建筑室内环境通风模拟研究

自然通风是一种简单、经济而有效的通风方式，既能节省有限的能源，又不会造成环境污染问题。通过介绍建筑室内环境模拟过程中的多区通风网络宏观模型及计算流体力学CFD微观模型，并对两种模型的原理及计算方法进行了论述，在此基础上对一典型建筑户型做了通风模拟及结果分析。     

刻槽球形轴承气膜运动方程的一个计算方法

一、引言研究气体动压轴承时,必须计算气膜的运动。在设计轴承时,为了提高其承载能力和稳定性,一般在轴承的转子(或定子)上刻一些槽,这更增加了计算的复杂性。对这类问题已有一些计算方法,这些方法中有的是假定有无限多条槽(窄槽理论),有的是给出了稳态解,或在稳态解上加一个小扰动,这些方法都有一定的局限性。本文根据对气膜运动周期性的分析,给出了一个计算气膜运动方程的渐近周期解的方法,它不必假定窄槽,可以计算大偏心等一般情况。试算表明,计算结果和实验结果符合较好。     

燃油箱燃油连通管的出流能力研究

本文研究了直升机加油过程中燃油箱通油连通管出流能力的两种计算方法:一是将连通管的流动型式全部假定为孔口出流进行计算;一是将连通管的流动型式细分为堰流及孔口出流,并建立弓形堰的出流计算公式进行计算,对这两种计算结果的差异进行了比较,供实际应用时参考.     

D φ—ψ流线坐标系中求解跨音速绕流问题

本文提出了在一种新的流线坐标系，位函数一流函数坐系系中求解跨音速流问题的数值方法。由于跨音速流动对边界条件的变化比较敏感，现在数值计算跨音速流动已大量采用边界拟合坐标系，流线会标系也是一种边界拟合坐标系，且由于流线随着计算的进而而自动调整。因此，本方法与近来发展较快的自适应网格技术有许多相似之处。数值计算表明，本方法在稳定性和收敛性方面都比较好，与实验结果及其它数值计算结果相比，相互吻合程度令人满     

超视距空战仿真中的策略识别

针对超视距空战仿真中敌机策略的识别问题,研究了一种基于案例的策略识别方法。该方法通过构建包含假定的对手任务目标、观测数据、策略等内容的案例库,采用相似度计算选择与新的观测相似的案例库子集,并计算策略概率分布来识别对手策略。实验证明,相比与常规的基于案例推理方法,在案例中增加敌机任务目标提高了策略识别准确率,并且在假定敌机目标不正确时,能修复错误假定并进行策略识别,改善了基于案例策略识别方法的性能。     

弹道再入飞行器的无粘绕流及表面加热率的一种工程数值解法

RVSPHR是一种确定弹道再入飞行器无粘气动力和表面热流率的近似数值计算程序。在头部亚跨音速区内,由工程方法给出激波角及压力分布,并假定激波层内的物面法向压力和法向速度剖面,而相应的激波形状则按质量守恒条件迭代调整。下游超音速区仍由严格运动方程差分数值求解。在无粘流计算的基础上,用无需迭代的熵吞技术,计算表面热流率。所有计算结果与严格数值解及实验值作了比较,其一致性较好。     

喷雾运动初始条件问题的研究

根据所完成的流动环境中颗粒在初始区域内速度分布的测量结果，对颗粒运动初始条件的问题进行了全面的研究，提出了统计性，随机性和确定性的三种初始条件及其数学处理方法，在此基础上，应用喷雾两相流的数值方法，计算了喷雾的运动，比较了不同初始条件对计算结果的影响，并将所得计算结果与实验结果做了比较，两者符合较好。     

湍流射流扩散火焰的层流火焰面模拟

采用双尺度κ-ε湍流模型，标量联合的概率密度函数（PDF）输运方程和层流火争面模型相结合，模拟湍流射流扩散火焰，比较了几种不同κ-ε湍流模型对轴对称湍流射流速度场模拟的结果，建立了混合物分数和湍流频率的联合PDF输运方程，并假定湍流频率和标量 线性关系，得到标量耗散率的概率分布，化学组分的浓度和温度由层流火焰面系综的平均得到，计算结果与直接用化学反应动力学计算的结果及实验数据做了比较。     

用神经网络法实现反向工程中的曲面插补

提出了一种用神经网络方法对反向工程中曲面的坐标．点进行插补，将插补点与NURBS计算结果进行了比较，结果证明这是一和非常有效的构造曲面方法．并对这种方法进行了探讨。     

民航飞机起飞决断速度的一种计算方法

起飞决断速度是民航飞机起飞过程中一个至关重要的参考速度．其准确与否直接关系到飞机的起飞安全。传统的查表方法确定起飞决断速度有时不够准确，所以提出了一种迭代求解方法．并通过B757—200大型运输机的基本参数计算了其起飞决断速度，将计算结果与查表所得的结果进行了分析比较，证明了该算法是切实有效的．既满足适航要求，又利于工程计算。     

基于数值模拟的流场附面层边缘识别方法

附面层边缘通常取在速度达到主流速度0.99倍的位置,而复杂流场中主流流动往往并不均匀,给附面层边缘的准确识别造成了困难。为解决此问题,提出了用"参考主流"代替实际主流识别附面层边缘的方法:通过零剪切力滑移壁面边界条件下数值模拟得到不受附面层干扰的参考主流,在根据附面层定义确定附面层边缘时以该参考主流中的速度代替实际的主流速度。通过斜楔压缩和弯曲压缩两个超声速压缩流场对该识别方法进行了验证,所得到的斜楔压缩出口截面上附面层厚度与采用实际主流速度判断得到的厚度相对误差仅4.1%。根据该方法的识别结果对弯曲压缩型面设计进行附面层修正后,弯曲激波高度与无黏设计值之间的误差从修正前的2.0%降低至0.3%,压缩面末端压力的相对误差从修正前的6.6%降低至2.3%。该方法避免了指定主流速度的主观性,识别结果较为准确。     

组合抽吸对大转角扩压叶栅性能的影响

实验研究了低速条件下在端壁和吸力面同时进行附面层吸除对某大转角扩压叶栅性能的影响.在叶片表面及端壁进行了墨迹流动显示,并对叶栅出口截面进行了测量.基于实验的抽吸槽布置方式:端壁吸气主要影响区域是吸力面角区;吸力面抽吸可以减小角区范围,延迟叶片吸力面附面层转捩,改善中径处流场;组合抽吸则优化了叶栅整体流场,使流动更加均一高效,由于削弱了端壁和吸力面附面层间的相互作用,组合抽吸在大吸气量下优于前两种吸气方式.      

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动。该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递。分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法。流场求解采用Spalart-Allmaras带湍流模型的Navier-Stokes方程组,其中对流项采用流通矢量分裂和5阶WENO(weighted essentially non-oscillatory)-Z格式离散,黏性项采用6阶中心差分格式离散,时间项采用龙格库塔显式格式离散。数值验证算例表明:该方法具备高于4阶的空间求解精度,并适用于刚体动网格非定常流场模拟,且无需更新网格形状。等距面贴体网格生成过程简单,避免了繁琐的人工调节过程,与笛卡儿网格结合,可提供足够的壁面附近网格密度,同时有效减少了网格总量需求。     

嵌边法出流条件在可压缩流直接数值模拟中的应用

以嵌边法为出口边界条件试算了亚、超声速平板边界层扰动波的演化、超声速流中涡的传播、亚声速平板边界层尾缘扰动的演化等,以验证嵌边法作为出流条件在可压缩流中的有效性。然后,探讨了嵌边法有关参数的选取问题。     

Adjoint boundary sensitivity method to assess the effect of nonuniform boundary conditions

The performance of turbomachinery is largely affected by the nonuniform boundary conditions caused by the coupling between neighboring parts, such as the inlet distortion and hot streak. Existing works study this problem by comparing the flow fields with uniform and nonuniform boundary conditions, which is cost extensive. In this work a new and efficient method is developed by computing the sensitivities of arbitrary performance metric relative to the boundary conditions and this method quantita...     

凹型粗糙元对边界层稳定性的影响

为了能够利用eN方法对带有凹型粗糙元的平板边界层进行转捩预测.通过数值模拟和流动稳定性的方法研究了马赫数为4.5的超声速边界层中凹型粗糙元对扰动演化的影响.对两种尺度的凹型粗糙元下扰动沿流向的幅值和增长率分布进行了计算.结果表明:凹型粗糙元对基本流的影响只是局限在凹型粗糙元及其附近很小的范围内.凹型粗糙元对扰动幅值的演化有抑制作用,尺度越大抑制作用越明显.在凹型粗糙元的后面,凹型粗糙元对扰动幅值的增长率的影响很小,相对变化量在2%以内.验证了对壁面有凹型粗糙元的情况,可以通过流动稳定性分析加N值修正的方法进行转捩预测.      

压气机转子通道内尖区三维平均流场

在低速大尺寸压气机试验台上,借助旋转四坐标全电动探针位移机构,用锥形五孔压力探针分别测量压气机设计状态和近失速状态转子通道内尖区的三维平均流场,揭示压气机转子通道内尖部的流动结构及其变化     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

高超声速边界层的转捩问题

取40km高空处的气体参数,对来流马赫数为8、10和12等三种不同情况的平板边界层和楔角为20°的楔型体边界层,做了转捩发生位置的预测。所用方法为改进的eN方法。结果发现,壁面条件对是第一还是第二模态波决定转捩位置有很大关系。如对于传统的eN方法,如果取N为10作为转捩判断标准,对马赫数为12的楔体边界层,等温壁条件下由第二模态波决定的转捩位置距前缘13m左右,而平板边界层中两模态波所决定的转捩位置距前缘均超过50m。而采用改进的eN方法,平板边界层中两模态波所决定的转捩位置距前缘均超过10m,而对马赫数为12的楔体边界层,两模态波所决定的转捩位置距前缘均出现小于5m的情况。     

Experimental investigation of a counter-rotating compressor with boundary layer suction

This paper presents the aspirated modification of a dual-stage counter-rotating compressor which contains inlet guide vanes(IGVs), two counter-rotating rotors(R1, R2), and outlet guide vanes(OGVs).Uniform circular holes are circumferentially distributed over the rotors' tips on the shroud casing which are designed to remove the low-energy fluid near the shroud casing.OGVs are hollow blades with two slots designed on the suction side which can better control the flow on the suction surface through boundary layer suction.Related works about the experiments have been carried out since June 2012 and the effect of flow suction on the performance of the compressor is investigated in detail.Characteristic lines at a 70% corrected rotating speed are tested and those with higher rotating speeds will be studied in the near future.Experimental results indicate that boundary layer suction can improve the compressor characteristics and the best suction methodology varies along the operating line.At the near stall condition, suction from the R2 tip region can obviously increase the efficiency and the total pressure ratio, as well as improve the flow capacity.Isentropic efficiency can be maximally increased by 4.24% with an increase of 1.94% in massflow under a suction flow of 160 m3/h.Suction at the R1 position with a suction rate below 0.35% in a high flow situation can make the performance of the compressor better than others.Around the peak efficiency point, boundary layer suction from the slots of OGVs is the best choice in improving the efficiency, but it causes a 0.1% loss in the total pressure ratio.