液氢/液氧火箭发动机喷雾燃烧过程三维数值模拟

用三维湍流Ｎ－Ｓ方程及拟流体模型描述液体火箭发动机内部喷雾两相燃烧流动过程,两相之间的质量交换采用蒸发模型计算,气相化学反应速率采用Ａｒｒｈｅｎｉｕｓ公式计算。采用ＳＩＭＰＬＥＲ与ＩＰＳＡ算法求解两相控制方程并对氢氧发动机燃烧流动过程进行了三维数值分析,得到了发动机推力室内燃气参数的详细分布情况。结果表明对氢氧发动机而言,蒸发过程非常快,推进剂混合过程是燃烧过程的决定因素。     

二维超音速湍流混合特性研究

着重讨论顺流喷射情况下的超音速混合的流场特征及影响因素。简要介绍了顺流喷射混合的研究现状。同时，探讨了不同的来流参数对流场特征的影响及混合机理。计算采用了一种二阶ＴＶＤ格式有限体积法及Ｂ－Ｌ湍流模型求解二维湍流Ｎ－Ｓ方程。结果表明：所采用的计算方法能成功地抑制非物理振荡，适用于二维超音速内流的湍流混合特性研究。     

k-ω SST两方程湍流模型中参数影响的初步分析

Menter的k-ωSST两方程湍流模型在流体力学计算中有良好表现。利用均匀试验设计方法对k-ωSST模型中八个参数在小攻角亚声速、跨声速和大攻角三种典型的翼型绕流流场计算中的影响规律进行了分析。分析结果表明在附着流中参数a1是对流场的影响最大的参数;而在分离流中,参数a1对流场的影响与湍流产生与耗散等项的作用相当,模型中的四个参数:σω2、β2、β*和a1的取值都会对计算结果有较大的影响。     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

共轭梯度方法在求解不可压N-S方程中的应用

采用非结构网格上的SIMPLE(Semi-implicit Method for Pressure-linked Equations)算法、k-ε湍流模型求解了三维不可压Navier-stokes方程,对低速不可压粘性流场进行数值模拟,在求解压强方程时采用共轭梯度方法,求解其它变量方程时,采用预处理的BICG(biconjugate gradients)算法.对NACA0012翼型绕流流场和飞艇绕流流场进行数值模拟并对结果进行分析,取得较好的结果.     

液体推进剂火箭发动机推力室再生冷却通道三维流动与传热数值计算

应用湍流模型对液体推进剂火箭发动机再生冷却推力室通道的流动与传热进行了三维数值模拟,冷却工质为氢气,其密度、导热系数、动力粘度随着温度和压力而变化,冷却剂比热容及金属固体物性随着温度而变化。计算采用标准k-ε双方程湍流模型及气-固耦合算法。结果表明:推力室燃气侧壁面的温度和热流密度的最高点均发生在喉部附近,喉部横截面固体区域最大温度梯度靠近燃气,喉部附近氢气在垂直主流方向的截面上产生了二次流。气固耦合面最大热流密度及最大对流换热系数同样位于推力室喉部附近。      

双燃式冲压发动机中超燃燃烧室冷态流场数值模拟

采用Ｍａｃ Ｃｏｒｍ ａｃｋ 矢通分裂格式, 求解雷诺平均的Ｎ－Ｓ方程组, 采用Ｂａｌｄｗ ｉｎ－Ｌｏ－ｍ ａｘ 代数湍流模型和混合长度模型, 模拟超燃燃烧室内两股平行超声速来流的冷态流场, 取得了与试验相一致的结果, 计算表明回流区和射流流动参数对超燃燃烧室流场内波系的变化有着较大的影响。     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

半波GLM方法封闭的湍流方程对转捩过程的数值模拟

ＴＲＡＮＳＩＴＩＯＮＳＩＭＵＬＡＴＩＯＮＵＳＩＮＧＧＬＭ－ＥＭＢＡＳＥＤＣＬＯＳＥＤＥＱＵＡＴＩＯＮＳＯＦＴＵＲＢＵＬＥＮＣＥＧａｏＧｅ（４ｔｈＤｅｐｔ．ＢｅｉｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏｎａｕｔｉｃｓ，Ｂｅ...     

Analysis of pressure characteristics under laminar and turbulent flow states inside the pilot stage of a deflection flapper servo-valve: Mathematical modeling with CFD study and experimental validation

Electro-hydraulic servo-valves are widely used components in the mechanical industry, aerospace and aerodynamic devices which precisely control the airplane or missile wings. Due to the small size and complex structure in the pilot stage of deflection flapper servo-valves, accurate mathematical models for the flow and pressure characteristics have always been very difficult to be built. In this paper, mathematical models for the pilot stage of deflection flapper servo-valve are investigated to overcome some gaps between the theoretical formulation and overall performance of the valve by considering different flow states. Here, a mathematical model of the velocity distribution at the flapper groove exit is established by using Schlichting velocity equations for in-compressible laminar fluid flow. Moreover, when the flow becomes turbulent, a mathematical model of pressure characteristics in the receiving ports is built on the basis of the assumption of the collision between the liquid and the jet as the impact of the jet on a moving block of fluid particles. To verify the analytical models for both laminar and turbulent flows, the pressure characteristics of the deflection flapper pilot stage are calculated and tested by using numerical simulation and experiment. Experimental verification of the theory is also presented. The computed numerical and analytical results show a good agreement with experimental data.     

双对置二冲程柴油机扫气过程仿真研究

基于双对置二冲程柴油机直流扫气过程流动特点,建立了该发动机扫气过程瞬态流动三维计算流体动力学(CFD)仿真模型,应用AVL FIRE软件对该发动机额定工况(功率160 kW,转速2 500 r/min)下的扫气过程进行了仿真研究,结果显示该工况下发动机扫气效率可达94%.通过扫气过程不同时刻流场分析发现,采用直流扫气方式可较好避免新鲜充量与废气的相互掺混,有效提高高速二冲程发动机扫气质量.      

计算两相流的一种改进方法

针对两相流动中密度及颗粒浓度的变化，建立了一种考虑颗粒相体积浓度值变化的颗粒相压力修正方程及考虑可压缩性影响的流体相压力修正方程，在不可压两相流的计算方法上，发展了一种计算可压两相流的方法，以速度协变分量为求解变量，对不可压两相密相流及可压稀相两相流进行了计算，并和已有的实验数据及有关文献进行了比较，吻合性较好，表明计算方法可行。     

径向开槽药柱发动机内流场数值模拟

采用时间相关的显示ＭａｃＣｏｒｍａｃｋ预估－校正两步格式和可压缩全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了轴对称径向开槽药柱发动机内流场。计算中采用了贴体坐标系统和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。边界参数计算采用参考平面上的特征线法和物理边界条件。通过大量算例的分析表明，方法有效，计算结果合理。当有集中质量加入时，燃烧室头部压强升高，喷射强度和开槽口的倾斜方向对头部压强也有较大的影响。     

旋转坐标系中可压缩湍流的基本方程和湍流模型

本文在可压缩湍流瞬时量基本方程的基础上,用质量加权平均的概念,推导了旋转坐标系中可压湍流流动的连续方程、动量方程和能量方程。进而导出了旋转坐标系中可压湍流的应力输运方程、脉动平均动能方程和湍流耗散率输运方程。这些方程中包括了由旋转引起的 Coriolis力和离心力的作用项。这些方程是在旋转坐标系中研究平均 N-S方程求解湍流流场和研究湍流模型的基础。      

轴对称弯道流动的N—S方程数值解

利用涡量—流函数方法求解离心压气机进气弯道中的轴对称粘性流场。本文推导出了q—τ湍流模型方程在正交曲线坐标系下的表达式,并用以预测湍流流场。文中在涡量方程及q—τ方程中引进非定常项,采用时间推进法求解。通过对二维直通道及复杂边界的轴对称弯道内的层流、湍流流动的计算,获得了令人满意的结果。      

近壁湍流相干结构的数值模拟

提出了近壁湍流对称单个相干结构的理论模型,采用紧致有限差分和Fourier谱展开相结合的数值方法研究它的演化问题,得到雷诺应力分布及高剪切层的形成过程与实验结果一致。     

Unsteady flow phenomena associated with leading-edge vortices

This paper presents selected results from extensive experimental investigations on turbulent flow fields and unsteady surface pressures caused by leading-edge vortices, in particular, for vortex breakdown flow. Such turbulent flows may cause severe dynamic aeroelastic problems like wing and/or fin buffeting on fighter-type aircraft. The wind tunnel models used include a generic delta wing as well as a detailed aircraft configuration of canard-delta wing type. The turbulent flow structures are analyzed by root-mean-square and spectral distributions of velocity and pressure fluctuations. Downstream of bursting local maxima of velocity fluctuations occur in a limited radial range around the vortex center. The corresponding spectra exhibit significant peaks indicating that turbulent kinetic energy is channeled into a narrow band. These quasi-periodic velocity oscillations arise from a helical mode instability of the breakdown flow. Due to vortex bursting there is a characteristic increase in surface pressure fluctuations with increasing angle of attack, especially when the burst location moves closer to the apex. The pressure fluctuations also show dominant frequencies corresponding to those of the velocity fluctuations. Using the measured flow field data, scaling parameters are derived for design purposes. It is shown that a frequency parameter based on the local semi-span and the sinus of angle of attack can be used to estimate the frequencies of dynamic loads evoked by vortex bursting.     

Gao-Yong理性湍流方程

对近十年来Gao-Yong理性湍流方程的研究进行综述,给出矢量形式的可压湍流封闭方程组的推导过程。侧偏统计平均保留了湍流脉动量的一阶统计信息,在引入加权漂移速度对称性及正交各向异性后,导出了漂移流的连续方程、动量方程及机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组。方程组不含任何经验系数,不使用壁面函数,保留了NS方程的均化的非线性特性。其级数形式的能量方程与非线性现象多尺度层次现象相对应,具备了描述湍流统计平均流动及拟序结构流动的双重功能。已进行了大量算例验证,验证结果证明了Gao-Yong湍流方程对广泛范围的复杂湍流问题的适定性。     

S弯管中可压三元紊流控制方程和数值计算

本文通过对瞬态N-S方程的质量加权平均处理,对不可压缩流的k-ε两方程紊流模型做了可压缩性的修正。用有限差分法进行了数值计算。计算结果与实验数据吻合很好,这预示着数值方法可能会发展成为研究S弯管中高亚音速三元流动的一种手段。     

EXPERIMENTAL STUDY OF TURBULENT WINGBODY JUNCTION AND WAKE FLOW

ＥＸＰＥＲＩＭＥＮＴＡＬ　ＳＴＵＤＹ　ＯＦ　ＴＵＲＢＵＬＥＮＴ　ＷＩＮＧＢＯＤＹ　ＪＵＮＣＴＩＯＮ　ＡＮＤ　ＷＡＫＥＦＬＯＷＫａｎｇＨｏｎｇｗｅｎ；ＦａｎＭｅｎｇ；ＷｅｉＺｈｏｎｇｌｅｉ（ＤｅｐａｒｔｍｅｎｔｏｆＭｅｃｈａｎｉｃｓ，ＰｅｋｉｎｇＵｎｉ...