基于动态混合网格的多体分离数值模拟方法

在以往发展的动态混合网格技术的基础上,耦合六自由度运动方程计算和非定常流场计算,建立了一种多体分离问题的非定常数值模拟方法。首先采用刚体运动的六自由度运动方程,建立刚体运动轨迹计算模块,然后以该模块为纽带将以往建立的非定常计算方法与改进的基于Delaunay背景网格插值方法和局部网格重构方法相结合的动态混合网格生成方法耦合起来,形成了精度较高、效率较快的多体分离数值模拟方法。利用该方法对典型多体分离问题进行了数值模拟,计算结果与风洞实验结果吻合较好。     

网格拓扑对DLR-F6构型数值模拟的影响

基于TRIP3.0软件平台和多块结构网格技术,开展了控制方程和网格拓扑等因素对DLR-F6构型数值模拟的影响研究。数值计算采用与试验相同的参数,采用了RANS/TLNS方程,生成了O型和H型2种拓扑的粗、中、细网格进行模拟。分别从气动特性、压力系数分布曲线和表面流态3个方面对结果进行分析。通过与试验数据的对比表明,在模拟小分离流动时,采用H型拓扑网格和RANS方程模型,获得的计算结果会更为准确。     

低雷诺数翼型层流分离现象大涡模拟方法

为准确预测翼型低雷诺数条件下出现的层流分离流动现象,发展了基于结构化拼接网格的大涡模拟方法。控制方程为Favre滤波Naver-Stokes方程,并选用多种亚格子模型。空间离散采用AUSM格式以及高阶WENO格式,时间推进采用显式方法和隐式方法。以SD7003翼型在雷诺数60 000及4°迎角下的层流分离流动为研究对象,对比分析了数值格式、亚格子模型、网格尺度对流场预测结果的影响。数值结果表明:划分的计算网格能够有效解析小尺度流动结构,基于隐式亚格子模型、采用AUSM格式和双时间步推进的大涡模拟方法能够准确预测流动分离、转捩、再附等复杂流动现象,计算得到的平均压强系数与雷诺应力分布与文献数据吻合较好。与转捩模式计算结果对比进一步表明:发展的大涡模拟方法能够准确预测翼型低雷诺数层流分离流动非定常演化过程,为下一步的研究工作建立了有效的数值模拟手段。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

薄片分离过程流场特性的数值仿真研究

利用弹簧近似法和局部网格重构相结合的非结构动网格方法,耦合求解Navier-Stokes方程和六自由度刚体动力学方程,数值模拟了三维薄片非定常分离过程。数值仿真结果有效反映出薄片分离过程中运动速度、轨迹、受力变化规律及流场特性的发展与变化。根据数值仿真结果对薄片分离机理、多体气动干扰作用进行了分析。分离过程中单个薄片与主体段周围形成复杂流动,并伴随着尾迹涡的运动、发展等非定常过程。气动载荷的非线性变化是薄片分离的主要原因。     

多体全自由飞过程的数值模拟

以三维非结构粘性直角网格重构为构架,数值求解了固定空间网格瞬间运动边界的N-S方程,逼近了流场的非定常过程,结合六自由度的弹道方程,对单体或多体运动过程进行了全自由飞数值模拟,数值模拟结果和级间分离的风洞自由飞试验进行了对比,初步的计算结果说明,全自由飞数值模拟在一定程度上可以描述风洞试验的飞行轨迹,并且适应能力更强,整个全自由飞数值模拟系统可以为更复杂的单体机动飞行或多体运动过程提供快速实用的数值模拟方法.     

失速限制型风轮叶片流场的数值模拟

采用双时间步方法在旋转坐标系下求解非定常方程,构造了零偏航条件下风轮叶片流场的数值模拟方法,采用该方法对失速限制型风轮叶片NREL Phase-Ⅵ进行了模拟,其中网格采用重叠网格进行模拟,控制方程为可压缩N-S方程,同时加入了低速预处理技术,模拟的风速范围从13m/s到25m/s,这时候风轮叶片发生失速以限制扭矩输出.本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析,从风轮扭矩、载荷分布等方面对数值计算结果进行了对比验证.     

螺旋桨流场Euler/NS数值解对比分析

采用多块点对点网格生成技术,生成螺旋桨流场空间计算网格。针对不同前进比,分别求解Euler和NS方程对螺旋桨流场进行数值模拟,并与实验结果进行了对比分析。研究发现无粘和有粘的计算结果均接近于实验结果,因此基于Euler方程的流场数值模拟更适用于工程设计。     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

一种新型颤振激励系统的数值计算研究

针对一种新型颤振激励系统的特点,采用非结构旋转运动网格技术,通过数值求解非定常欧拉方程,计算了该激励系统工作时的非定常气动力特性,并与准定常计算结果进行了比较。研究结果表明:采用非定常方法所得结果与准定常计算结果有明显的差别;气动力的变化与激励系统中圆柱的旋转方向有关,旋转角速度越大,气动力变化的幅值越小;气动力的幅值与圆柱缝道宽度也有一定关系。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

非定常环境中叶栅总压恢复系数的数值研究

利用以Ni格式和隐式双时间步法求解二维非定常N-S方程为基础的数值模拟方法研究了流动的非定常性对涡轮叶栅总压恢复系数的影响,并与文中提出的定常模型和准定常模型进行了对比。结果发现,非定常流动能够提高叶栅的总压恢复系数,但这是建立在与准定常流相比较的层面上,而和定常流比较时,非定常流只会使总压恢复系数下降。      

跨音速机翼结构动力响应的计算研究

通过交替求解流体运动方程和结构动力方程来计算时域内二维机翼的结构动力响应,并研究了结构参数对系统稳定性的影响,使用控制中的相平面图对计算结果进行分析,来判断系统的结构稳定性。计算中流体运动方程采用二维非定常Ｅｕｌｅｒ方程,并用Ｊａｍｅｓｏｎ提出的有限体积方案,四步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进求解;结构模型采用二自由度ＰＡＰ模型。经分析,计算结果合理。     

基于准定常的飞行器动导数的高效计算方法

结合工程设计经验和工程计算方法,针对飞行器动导数计算方法展开研究。基于求解准定常欧拉方程发展了一种适合于亚跨超声速状态下,任意气动外形气动导数的高效计算方法,并将计算结果与一般非定常方法的求解结果以及国际标准实验数据进行了对比,研究结果表明该方法结果可靠、效率较高。     

风扇进气畸变三维非定常数值模拟技术研究

采用三维欧拉方程加体积力模型的方法, 发展了一套三维非定常进气畸变数值模拟程序, 以实现风扇/压气机在进气畸变条件下的全环面三维非定常数值模拟.利用该程序模拟计算了某航空发动机风扇进口稳态总压畸变条件下特性和流场结构, 并与试验结果进行了对比, 结果表明该程序用较小的计算机资源即可以实现对风扇全环面三维非定场畸变流场的模拟, 能够以较高的精度预估气流畸变对压气机/风扇特性的影响.      

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

基于闭环控制的带翼导弹虚拟飞行数值模拟

通过耦合求解非定常Navier-Stokes方程、刚体六自由度（6DOF）运动方程及舵偏控制律,发展了一套虚拟飞行数值模拟平台,对基本带翼导弹外形进行了虚拟飞行数值模拟研究。流场求解采用基于结构网格的有限体积法,时间推进采用双时间步法,使用刚性重叠动网格技术模拟舵与弹身的相对运动,采用四阶Runger-Kutta方法求解运动方程组,舵偏控制律采用基于比例-积分-微分（PID）的反馈控制。模拟结果表明：发展的虚拟飞行数值模拟平台能够处理包含复杂运动边界的非定常运动问题,可以为非定常气动力建模提供数据、预示飞行器稳定性及检验飞行控制律,该模拟平台具备一定的工程应用价值。