叶片安装角槽式轮毂处理的三维数值模拟

采用三维定常不可压Ｎ－Ｓ方程对叶片安装角槽式的轮毂处理的叶栅通道进行了数值模拟,分析了叶片安装角槽式的轮毂处理的扩稳机理,在计算中采用了非交错网格的ＳＩＭＰＬＥ算法,并解决了因使用非交错网格离散方程而带来的压力场振荡的抑制问题,紊流模型采用ｋ－ε两方程模型。     

跨洞庭湖大桥主梁节段定常气动力研究

采用风洞实验与数值模拟的相结合的手段,对跨洞庭湖大桥主梁节段的定常气动力特性进行了研究;流场数值模拟是采用ｋ－ε湍流模型封闭的雷诺时均Ｎ－Ｓ方程,求解方法是建立在交错网格了的压力－速度场迭代求解的ＳＩＭＰＬＥ算法;实验是在国防科技大学ＫＤ－０３低速风洞进行,研究表明实验与计算结果较为一致。     

高速列车绕流的数值模拟

给出了用低速Ｎ－Ｓ方程求解高速列车绕流的数值计算方法，采用显隐式结合的方法来加快收敛速度，基于该方法编制了计算软件系统Ｔｒａｉｎ ３Ｄ，用该软件系统对高速列车外形的外流场进行了数值模拟，获得了该外形的空气动力特性和流场结构。     

翼型喷气增升效应的数值模拟

应用数值方法模拟了飞行器升力面下表面喷气的增升作用。数值模拟的出发方程是Ｎ－Ｓ方程,格式为Ｂｅａｍ－Ｗａｒｍｉｎｇ格式的改进型。模拟结果表明：喷气增升主要是通过改变上、下翼面的压力分布实现的。在一定范围内升力随气流出口速度增加而加大,气流出口速度超过一定值后升力开始下降。喷口位置后移增升效果较好。单一增大气流出口压力时增升效果无明显变化。     

旋翼近地飞行时诱导速度的N—S方程计算

采用有限体积法用三维Ｎ－Ｓ方程对旋翼近地飞行时的速度场进行了数值计算，对旋翼桨盘处的诱导速度进行了分析。旋翼对流场的作用由分布在特定区域内的动量源项模拟，动量源项的强度由旋翼的几何尺寸，叶片剖面的空气动力特性及当地流场的性质迭代定出。本文讨论了旋翼近地飞行时，在不同前进比及不同离地高度情况下的诱导速度随径向距离及方位角的变化，并将数值计算结果与实验结果进行了比较。     

高Ra数自然对流问题的张力样条数值模拟

本文将基于压力泊松方程技术的原参数Ｎ－Ｓ方程的张力样条数值解法推广应用于自然对流问题的数值模拟，提出了自然对流问题一个新的数值模拟方法，对一方内高Ｒａ数自然对流问题的数值模拟表明，结果是令人满意的。     

真空羽流场的 N-S 和 DSMC 耦合数值模拟

根据调姿发动机真空羽流场的特点,采用差分求解Ｎ－Ｓ方程数值模拟喷管内流场,用ＤＳＭＣ方法数值模拟喷管外羽流场,得到了符合流动规律的数值模拟结果。为姿控发动机喷管设计提供了参考,并为真空羽流污染研究打下了良好的基础。     

S弯进气道中自动定位叶片对气流流动影响的研究

对Ｓ弯扩压管道内安置自动定位叶片后的气流流动特性进行了实验研究。结果表明，位于模型进气口处的可转动叶片在气流作用下能有两个自动稳定的位置。如果忽略可转动叶片转动轴的摩擦力矩，当叶片处在这两个自动稳定位置时，气流对叶片的合作用力矩为零。研究结果还进一步表明，自动定位叶片对Ｓ弯进气道出口压力场与速度场起着重要作用。本研究为以后如何采取叶片技术抑制Ｓ弯进气道出口的旋流流动，以及模拟旋流流动提供了重要技术基础。     

超燃冲压发动机燃烧室三维化学反应流场的数值模拟

采用ＬＵ隐式方法强耦合求解有限体积法离散的完全Ｎ－Ｓ方程及组分方程，数值模拟超燃冲压发动机机燃烧室机的化学反应流流场。计算中和了ＡＵＳＭ通量分裂格式，两方程的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代化学反应流场。     

细长体简谐振动绕流的数值模拟

用数值模拟方法研究了椭圆锥细长体简谐振动非定常过程的流场特性。数值模拟的出发方程是可压缩流的全Ｎ－Ｓ方程，计算格式为ＴＶＤ格式。数值模拟中，研究了细长简谐振动绕流的非定常流场特性，不同振动缩减频率和压缩性效应的影响，揭示了非定常流场的特性和流动机理。     

Aerodynamics of smart flap under ground effect

Aerodynamic characteristics of two-dimensional smart flap under the ground effect have been assessed by a numerical simulation. In this process, a pressure-based implicit procedure to solve Navier–Stokes equations on a nonorthogonal mesh with collocated finite volume formulation is used. The boundedness criteria for this procedure are determined from the Normalized Variable Diagram (NVD) scheme. The procedure incorporates the k–ε eddy–viscosity turbulence model. Cantilever beam with uniformly varying load with roller support at the free end is considered for the configuration of the smart flap. The method is first validated against experimental data. Then, the algorithm is applied for turbulent aerodynamic flows around airfoil with smart and conventional flaps for different attack angle, flap angle and ground clearance where the results of two flaps are compared. The comparisons show that the quality of the solution is considerable.     

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

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

激光功率和工质对连续激光推力器性能影响的数值模拟

对连续激光维持的等离子体加热推力器内流场建立一种计算模型,考虑的物理机制包括化学反应、高温气体性质、激光吸收、黏性、扩散、热传导以及辐射效应,模型方程形式为可压缩轴对称Navier-Stokes方程,对方程采用SIM-PLEC算法求解.在吸收室压强、聚焦光束形状参数及喷管构型相同的条件下,模拟了氩为工质的推力器比冲随激...     

任意非正交曲线坐标系中湍流的数值计算

本文旨在运用张量分析理论,结合成熟的SIMPLE算法,发展一种适用于复杂几何形状内湍流流场数值求解的方法。假设流动为二维单连通域内的等温、不可压、定常、无化学反应、单组分湍流流动,采用双方程湍流模型。由雷诺平均的湍流控制方程的基本形式[1,2],结合张量分析理论。      

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

叶轮机械三维粘性动静叶干涉的数值模拟

本文发展了 1套模拟叶轮机械三维粘性动静叶干涉的数值方法。在成熟的空间离散方法基础上 ,提出了 1种高效的双时间步隐式时间推进方法 ,并构造了一种保证守恒的动静叶交界面处的数值通量计算方法。本文对某三维亚音速涡轮的动静叶干涉进行了数值模拟 ,结果与实验符合较好。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

超声速小攻角下的圆锥转捩攻角效应

采用高阶差分格式求解三维可压缩滤波Navier-Stokes方程,对超声速零攻角和小攻角圆锥边界层转捩进行了空间大涡模拟研究.在三维Tollmien-Schlichting(T-S)波扰动作用下,计算得到了流场非定常转捩过程.计算结果清晰地展现了非定常大尺度结构演化特征和攻角效应引起的非对称转捩现象,流场发卡涡、横流涡等转捩主控大尺度结构与实验一致.     

基于黏性涡模型的旋翼流场数值方法

 建立了一种适用于旋翼非定常流场特性分析的黏性涡数值方法。在该方法中:流场中的大尺度涡被离散为若干微小的涡元,通过求解涡量-速度形式的Navier-Stokes方程模拟涡元的输运等过程;黏性扩散效应采用高精度的粒子强度交换法进行计算,而桨叶附着涡以及新生涡环量采用了Weissinger-L升力面理论进行求解;为显著提高计算效率,在诱导速度及其梯度的计算中还引入了快速多极子算法(FMM)。应用上述方法,对悬停和前飞状态下的多个旋翼流场算例进行了计算,通过对比旋翼尾迹涡量特征和诱导速度分布等,验证了该方法的有效性。此外,还将本方法与旋翼计算流体力学(CFD)方法及传统的自由尾迹方法进行了比较,结果表明黏性涡方法在兼顾效率的同时,还能够更好地捕捉旋翼尾迹运动。     

高超声速飞行器表面温度分布与气动热耦合数值研究

针对高超声速飞行器热防护设计中的高温气体非平衡效应问题和气动热环境精确预测问题,基于流场的非平衡Navier-Stokes方程、表面的能量守恒方程和内部的热传导方程,考虑流场的非平衡效应、表面的热辐射效应、催化效应和烧蚀效应以及热防护层内部的热传导效应,建立了初步的表面温度分布与气动热的耦合计算方法,完善了高超声速飞行器气动物理流场计算软件(AEROPH_Flow)。在表面材料为碳-碳(C-C)条件下,对飞行高度为65km和飞行速度为8,10km/s的半球以及飞行高度为50km和飞行速度为8km/s的球锥模型,开展了表面温度分布与气动热的耦合计算,验证了计算方法和计算软件,分析了表面温度分布对气动热环境的影响。研究结果表明:表面温度分布对气动热的计算结果有较大影响,在气动热环境的预测中,不仅要考虑热化学非平衡效应和表面催化效应的影响,还要考虑表面温度分布的影响,最好是采用表面温度分布与气动热耦合计算的方法,以减小表面温度分布对气动热计算结果的影响。为此,需要发展完善非平衡流场/表面催化和烧蚀/热传导温度场(气/表/固)的计算模型、耦合求解技术和计算软件,实现对高超声速飞行器的真实飞行条件下高温气体非平衡效应和气动热环境的精确模拟。