N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

二维横向喷流的N-S方程数值解

让一股射流横向注入超音速流产生所需要的干扰效应是一种飞行器飞行控制的常用手段,如航天飞机和机动导弹上的反作用控制系统(RCS).所以,对射流与主流相互干扰的某些流动特性应该有个基本的了解.本文通过应用MacCormack显格式和Baldwin-Lomax修正的代数湍流模型求解二维RANS方程对带横向射流的绕后台阶的M_∞为2.19和6.00的超音速外干扰流场进行了数值计算.与无喷射流的情形相比,回流涡旋区由2个增加到4个,且底部压力有大的回升,喷口宽度越大,回升效应越显著.另外,数值计算表明,在压力梯度变化比较平缓的区域附加一四阶人工粘性能有效地抑制数值振荡,加快收敛.     

湍流的通用物理方程

利用变形张量及摄动变形位移的点积分解湍流脉动速度,并使用密度加权的侧偏系综平均方法取代 Reynolds平均方法,在不引入任何物理假设的情况下,对 Navier-Stokes方程进行平均。未知相关量在三阶精度上才出现,截断之后,获得具有二阶精度的封闭的湍流方程组。该方程组无相关量,不引入任何经验系数,无需构造任何湍流模型,即可描述 Ma<3.6流速的任意复杂湍流场并解释湍流的物理本质。本文仅对推导思路及方程组做概要介绍。      

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

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

二元喷管跨音速流场的数值计算

本文采用有限体积法进行空间离散,Runge-Kutta法进行时间推进的办法,通过求解Euler方程来模拟二元收敛—扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡,在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解。计算结果与试验数据和其他数值方法吻合良好。      

基于ADE-ELM的涡轴发动机建模方法

提出了基于自适应微分进化-极端学习机(ADE-ELM)求解平衡方程的高精度涡轴发动机实时部件级模型建立方法.基于牛顿-拉夫逊(N-R)迭代模型,以迭代计算前模型平衡方程残差为输入,迭代收敛后平衡方程猜值修正量为输出,训练极端学习机,并采用自适应微分进化(ADE)算法优化极端学习机(ELM)参数,提高猜值修正量映射精度.ADE算法中采用sigmoid型自适应缩放因子,提高了微分进化算法的寻优能力.在涡轴发动机不同飞行状态下的测试结果表明,以N-R迭代算法模型为基准,基于ADE-ELM的发动机模型,最大建模误差约为一次通过算法的1/3,运算耗时约为一次通过算法的1/3,验证了算法的有效性.      

RP-1煤油凝胶的制备与性能表征

为了研究凝胶燃料的流变机理和配方设计的基础理论,根据凝胶特殊的固-液形态,寻找凝胶化航空煤油RP-1需要的无机和有机凝胶剂,对其成胶机理、流变表征和稳定性进行研究,通过Brookfield流变仪测量流变参数,并构建出无机和有机煤油凝胶的HBE本构方程。结果表明:采用高速剪切分散比三辊研磨分散制备的无机煤油凝胶质量好;加入无水乙醇有利于凝胶剂A对煤油的凝胶化;无机煤油凝胶的物理稳定性受凝胶剂含量影响显著,实际应用选择6%和7%无机凝胶;所制备的煤油凝胶比典型的剪切变稀流体(0n1)的稀化能力更强;无机和有机煤油凝胶的复凝性很弱且不受凝胶剂含量的影响。     

固体火箭发动机喷管燃气舵铰链力矩参变分析

以固体火箭发动机喷管燃气舵绕流流场为研究对象,采用数值模拟技术,通过离散求解流动方程及相关湍流模型方程,对包含舵基、舵片和发动机喷管在内的复杂流场进行了变参数分析,获得并揭示了舵基、舵片附近的流动现象和流动特点.验证了所引入的数值模拟方法的有效性,分析结果对舵基、舵片的结构安装以及控制机构的设计具有借鉴意义.     

考虑低温增压真实气体效应的运输机气动特性数值模拟研究

为了给低温风洞试验数据修正提供参考,本文利用数值模拟手段研究了低温真实气体效应相比较于完全气体对飞行器气动特性的影响,以及该影响与雷诺数影响相比所占比例的大小等问题。文章应用 Aungier-Redlich-Kwongz 方程,发展了适用于模拟氮气低温高压真实气体效应的 RANS 求解软件。与 NIST 数据的对比表明,该状态方程在5倍大气压下,比热比等参数误差在0.3%以下。同时,标模测试结果表明本文软件计算精度与国外软件相当。应用本文方法研究低温增压风洞中氮气真实气体效应：以典型运输机构型 DLRF6为模型,分别计算了高速、低速状态下各种不同温度和压力工况下的流动。计算结果表明,在低温增压情况下,真实气体效应引起的气动力差异很小,升力、阻力、力矩最大相对误差均在0.3%以下,与雷诺数效应引起的偏差相比可以忽略不计。因此,以氮气为介质的低温风洞试验研究可以采用完全气体假设。     

不同襟翼偏角梯形翼构型气动特性数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用二阶离散精度的单调迎风格式(MUSCL),结合剪切应力输运(SST)两方程湍流模型和γ-Reθ转捩模型,研究了梯形翼高升力构型襟翼偏角变化对气动特性的影响。主要目的是进一步确认Trisonic Platform(TRIP)软件模拟高升力梯形翼不同襟翼偏角引起微小气动特性变化的能力。首先,简要介绍了采用的计算方法;其次,介绍了两种襟翼偏角的梯形翼模型及风洞试验;最后,在网格收敛性研究的基础上,采用全湍流和转捩两种方式模拟了梯形翼构型不同襟翼偏角对气动特性的影响。与试验数据的对比结果表明,采用全湍流和转捩两种方式均可以较好地模拟不同襟翼偏角对气动特性的影响量,采用γ-Reθ转捩模拟方式可以提高梯形翼构型气动特性的模拟精度,失速迎角附近的气动特性模拟需要进一步研究。