面向降落伞稳态CFD计算的网格生成方法研究

降落伞稳态流体动力学(ComputationalFluidDynamics,CFD)计算的网格生成是一个带边界层的三维复杂域网格生成问题,目前用于边界层计算的混合网格生成方法,往往存在计算繁琐、应用范围窄、自动化差、以及难以适应复杂外形的缺点。文章提出了一种结合约束德洛内(Delaunay)网格生成和网格前沿推进技术的降落伞稳态CFD计算的流场网格生成方法,实现了降落伞网格和包括边界层区域在内的流场网格一体化、全自动、高品质的生成。该方法采用网格前沿推进法来生成边界层区域的网格节点;算法整体上采用三维约束Delaunay网格生成技术,边界层层节点集合在网格生成过程中作为约束Delaunay三角化的约束条件,避免了复杂的网格求交计算和拓扑处理。网格生成实例表明,该方法能全自动生成降落伞稳态CFD网格,生成的网格品质、网格规模等满足降落伞稳态CFD计算的需求,同时该方法具有一定的通用性。     

应用TTM网格研究冲压增程弹丸进气道内外流场

基于TTM贴体网格生成法,研究了某冲压增程弹丸超音速进气道内外流场数值模拟中单区域网格生成问题,流场求解采用了二阶隐式TVD格式。NACA0012翼型跨音速流场算例结果表明,所采用的算法和编制的程序可靠,能够适用于具有亚、跨、超音速三种流态共存的复杂流场数值计算。对超音速进气道数值模拟结果显示,得到的流场复杂波系结构是合理的,来流马赫数和攻角变化显著影响着进气道性能。     

基于Van Leer+AUSM混合格式超声速流场的并行数值算法研究

基于CLUSTER集群并行体系结构与消息传递库MPI并行环境,开展了超声速流场的并行数值算 法研究。计算网格采用多块对接粘性结构网格,控制方程为积分形式的Navier\|Stoke s(N-S)方程,空间离散采用基于Van Leer+AUSM混合格式的有限体积法,粘性项采用中心 格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,应用一方程S -A模型计算湍流的粘性。并行计算采用 分区并行的思想,设计了满足负载平衡要求的区域分解和通讯方法,编程上采用单控制流多 数据流(SPMD)模型以及Master／Slave并行模式。最后,通过对某弹丸和某制导火箭弹超声 速流动的数值模拟, 验证了并行算法的准确性、高效性和处理复杂几何外形的能力。
     

球形贮箱中三维液体大幅晃动数值模拟

提出一种新颖有效的网格移动算法,使用任意拉格朗日-欧拉（ALE）有限元方法对球形贮箱中三维液体大幅晃动问题进行数值模拟。对接触线网格及自由液面网格指定其运动轨道,对湿壁面网格采用代数更新算法,对域内网格采用Laplace方程对网格速度进行平滑。同时采用基于特征线的时间分裂（CBS）格式来处理对流效应引起的数值振荡。数值计算结果与已有试验数据进行了对比,两者吻合良好,证明了网格移动算法与计算格式的可行性与准确性。     

Rankine涡作用下激波变形和分叉的数值研究

数值模拟了激波和Rankine涡相互作用，重点研究了波涡相互作用引起的激波结构的变化过程。应用非结构化适应网格下的二阶精度Godunov型的PLM格式来求解Elder方程，对一平面激波与不同涡强Rankine涡相互作用下瞬态激波的运动、发展进行了研究。漩涡与激波间相对强度的大小对激波结构的影响明显，弱相互作用诱导激波变形，强相互作用诱导激波变形并产生分叉。计算结果表明，应用非结构化适应网格的PLM格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息。     

Rankine涡作用下激波变形和分叉的数值研究

数值模拟了激波和Rankine涡相互作用,重点研究了波涡相互作用引起的激波结构的变化过程.应用非结构化适应网格下的二阶精度Godunov型的PLM格式来求解Euler方程,对一平面激波与不同涡强Rankine涡相互作用下瞬态激波的运动、发展进行了研究.漩涡与激波间相对强度的大小对激波结构的影响明显,弱相互作用诱导激波变形,强相互作用诱导激波变形并产生分叉.计算结果表明,应用非结构化适应网格的PLM格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息.     

非结构网格高超声速MHD流场逆风格式数值模拟

应用逆风通量分裂格式在非结构混合网格上对二维高超声速理想磁流体绕钝头体的流场进行了数值模拟。控制方程为Euler方程耦合Maxwell方程的理想MHD方程,空间离散采用AUSM格式,时间推进采用显式5步Runge-Kutta格式。引用双曲型散度清除技术加强·Ｂ＝0的条件。计算模型为二维钝头体,在高超声速来流条件下,对不同磁感应强度的均匀分布磁场干扰下的流场进行了计算,得到了较满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法可用于高超声速MHD方程的求解。     

基于预处理的非稳态流场计算方法

采用预处理方法可提高Navier-Stokes方程或Eu ler方程计算低马赫数流场的准确度和收敛性能。采用双时间步和预处理相结合的方法计算非稳态流场,保证了采用预处理方法后计算结果的时间精确性。在每一物理时间步,控制方程虚拟时间导数项采用LU隐式时间格式,对流通量采用AUSM 和MUSCL格式进行求解。通过对前台阶流和圆柱绕流的计算,验证了该方法的广泛适用性和计算结果的可靠性。     

高超声速飞行器外形热流密度分布计算的高精度方法研究

为提高高超声速条件下飞行器表面热环境的计算精度,发展了高精度算法WCNS-E-5,并开展了热流密度分布的数值计算研究.该算法采用的是五阶精度加权显式非线性格式,结合四阶精度的二阶偏导数差分近似、四阶精度的边界格式,同时还对网格偏导数及温度梯度也离散为四阶精度.通过求解三维Navier-Stokes方程,数值研究了高超声速飞行器外形如钝锥和双椭球体的物面热流密度分布.模拟结果表明,这种保证全流场高精度的WCNS-E-5得到的流场图像清晰、真实.与通常的二阶算法MUSCL对比,WCNS-E-5具有更高的流线分辨率,获得的热流密度更接近于实验测得分布.     

基于全隐式紧耦合算法的气动弹性数值仿真

为进行气动弹性问题的计算,提出了一种全隐式紧耦合算法,在子迭代过程中分别采用LU\|SGS隐格式和隐式线性多步法交替求解气动和结构方程,以获得物理时间域的高精度解。一种径向基函数和超限插值结合的方法被用来进行气动网格的快速变形。运用该算法,进行了Isogai wing和AGARD 445.6 wing的颤振分析,颤振边界的计算结果与文献值和实验值较符合,表明该全隐式紧耦合算法能够有效地计算气动弹性问题。