基于结构网格点的Euler方程无网格自适应求解算法

提出了一种自适应无网格Euler方程求解算法。构造了基于无网格点云结构的布点加密技术,实现了借助压力梯度变化准则的流场局部加密;借助流场结构化剖分技术,给出初始无网格节点分布,并采用四步显式Runge—Kutta时间格式推进计算,求解了绕翼型的典型流动问题。数值算例表明,本文提出的方法能有效地提高流场的分辨率,如捕捉激波等流动特征清晰明了。     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

处理动边界问题的无网格/直角网格混合算法

把无网格/直角网格混合算法发展用于涉及动边界的非定常流动计算问题.采用无网格/直角网格混合算法,计算区域整体采用直角网格,只在物面附近嵌入局部无网格.对于动边界问题,本文让主体直角网格保持固定,用原混合算法中的局部无网格区跟随运动边界,使得动边界物面附近有不变的点云结构,避免了点云重构;无网格区与直角网格交界面随边界运动跟踪确定,由于主体基于直角网格,方便了交界面的跟踪.用发展的混合算法,结合非定常Euler方程双时间推进方法,先对翼型作俯仰振荡的非定常跨声速流动进行了数值模拟,计算结果与现有文献和实验数据进行了比较与验证;结合刚体运动方程,成功模拟了二维外挂物投放过程,算例展示出发展的混合算法适合处理大位移运动边界问题.     

基于混合动网格的二维非定常粘性流动数值模拟

利用基于Delaunay三角化的动网格方法,在混合网格上求解了带运动边界的二维非定常Navier-Stokes方程。在有限体积法格心格式的基础上,采用带人工耗散项的Jameson中心格式以及双时间步长推进方法对Navier-Stokes方程中的空间项和时间项分别进行了离散。湍流模型为Spalart-Allmaras方程模型。对在跨声速粘性绕流中进行俯仰振动以及后缘摆动的翼型进行了数值模拟,得到了令人满意的结果。     

一种适合全流速的流场数值方法

发展一种适合于所有流速的压力修正方法，在考虑了密度的可压缩性后，压力修正方程夸为对流扩散模型的，从而与般标量一样使用二阶中心差分离散。整个算法具有二阶精度，进一步使用多步压力修正来增强网格非正交时的收敛性态，在对流项离散中，设计了独特的人工耗散，也同样适用于所有流速，利用本算法求解全三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，计算采用了非交错网格方案，使用Ｒｈｉｅ－Ｃｈｏｗ方法消除了速度－压力的不关联，     

基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。     

非结构网格高超声速绕流数值模拟及涵道构型减阻特性分析

在非结构网格上对二维高超声速化学非平衡粘性绕流进行了数值模拟，并应用此方法对涵道构型高超声速减阻特性进行了数值分析。空间离散采用VanLeer逆风通量分裂格式，时间推进采用显式的Runge—Kutta格式。化学非平衡动力学模型为五组元五反应模型，对化学反应源项进行了点隐式格式处理，温度场的计算采用牛顿迭代法。分别对二维类弹头体涵道构型、二维半圆柱模型的高超声速绕流流场进行了数值模拟，得到了数值结果，并与完全气体绕流计算结果进行了对比分析。结果表明本文方法可以应用于高超声速绕流计算，并且涵道构型具有高超声速减阻特性。     

基于SNGR方法的二维空腔噪声数值模拟

基于SNGR方法,采用有限差分法求解带源项的线化欧拉方程。数值方法采用色散关系保持(DRP)格式。低耗散色散龙格-库塔格式显式时间推进,采用无反射远场边界条件。对二维空腔,长深比为9∶1,在亚声速、跨声速和超声速多个马赫数下进行噪声预测,并与实验结果对比,两者吻合较好。     

绘制等值线的一种新算法

本文以非矩形离散网格和任意离散数据点集的情况为背景,首次将数据结构引进等值线的求解中,给出在三角形网格上绘制等值线的一种快速、简洁、通用的算法。     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

动边界问题无网格算法及其动态点云

提出基于Delaunay图映射的点云运动处理技术.该技术使点云随物体运动再生成问题变成简单的映射关系,无需迭代,且能保持点云结构不变,故有效地避免了因点云结构变动而产生的附加计算量.运用该技术,结合非定常Euler方程双时间步无网格算法,成功地模拟了翼型的跨声速非定常流动.其中物理时间迭代采用二阶隐式格式,伪时间遮代采用四步龙格一库塔显式格式.文中给出了不同翼型的计算结果,并与实验结果进行了比较,取得了令人满意的结果.     

用于非定常粘性流动计算的无网格算法(英文)

把无网格算法发展用于求解涉及动边界的非定常粘性流动问题。在处理粘性流动的布点问题时,通过在远离物面的区域采用无粘流动计算时所用的各向同性点云,而在物面附近引入各向异性点云,从而在保证物面法线方向布点较密以准确模拟边界层的同时,有效控制了布点总量,减少了计算时间。对于动边界问题,本文在上述方法获得的初始布点的基础上,采用基于扰动衰减规律的点云移动技术,快速得到物面运动到其他位置时流场求解所需点云。在所获得的点云上,采用一种带权系数的二次极小曲面逼近方法来离散Navier-Stokes方程的空间导数,权系数的引入使得流场求解时对各向同性点云和各向异性点云无需分开考虑,而是采用统一的求解方法,从而简化了编程求解。用发展的无网格算法,结合Navier-Stokes方程求解双时间推进方法,并耦合Spalart-Allmaras湍流模型,先成功地模拟出NLR7301翼型后缘摆动诱发的非定常流,并通过与实验比较,验证了本方法.接着进行了NACA0012失速模拟,给出了动态失速过程,展示出用本文方法处理复杂动边界问题的效果。     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

圆柱形空腔内气云爆炸波对壁面的作用

从非定常Ｅｕｌｅｒ方程出发，用高精度高分辨率的ＴＶＤ格式，对圆柱形空腔内的气云爆炸波对壁面作用问题进行了数值模拟。计算获得圆柱形空腔内气云爆炸场的复杂的波系结构及其演化过程、壁面任意指定点的压力时间曲线。本文还给出了同体积、不同形状的气云爆炸对壁面的作用效应。     

一种划分航空知识点类型的方法

基于航空知识提出一种划分知识点类型的方法,对各种知识点类型进行统计,给出知识点的存储结构。     

弹尾超声速轴对称喷流的正格式数值模拟(英文)

采用二阶正格式方法对超声速轴对称喷流流动进行数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对导弹尾部超声速伴随射流进行了数值计算。数值结果与实验照片所反映的流动特征吻合较好,与三阶精度间断有限元方法计算结果吻合,相比于二阶高分辨率TVD格式的计算结果,正格式方法的计算结果在间断点处具有较大的梯度变化。这表明该方法对激波具有较强的捕捉能力,在间断处不会出现数值振荡,对弹尾声音速伴随射流的数值模拟真实、有效。     

超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

尖前缘任意翼型的超音速和高超音速俯仰安定性导数

本文推广了文[1]的欧拉方程摄动法,讨论了尖前缘任意翼型俯仰安定性导数的计算方法,并应用该方法计算了双圆弧翼型的俯仰安定性导数,其数值结果与文[2]的超音速线化位流方程有限差分法的数值结果比较,更接近于实验值。在激波附体的条件下,本文的方法可用于计算大迎角尖前缘任意翼型的超音速和高超音速俯仰安定性导数。     

Roe格式在多组元燃烧流场数值模拟中的应用

对Roe通量分裂格式进行了详细讨论，并将其推广应用到基于非结构网格的高超声速多组元燃烧流场的数值模拟中。控制方程采用三维多组元Euler方程，重点针对Roe格式开展研究；考虑到燃烧流场求解过程中出现的刚性问题，对化学反应源项采用点隐式方法处理。对超声速条件下钝头体激波诱导燃烧流场进行了数值模拟，比较分析了3种经典化学动力学模型对流场结构的影响，获得与数值模拟、试验数据相符的结果。结果表明发展的方法可以用于多组元燃烧流场的数值分析，且化学动力学模型的选取直接影响激波诱导燃烧流场中诱导区厚度。     

Stiffened状态方程下γ-model方法跟踪二维多介质流动

在Stiffened状态方程下，运用Level—set方程跟踪界面运动变化，把界面捕捉的等效方程、Level-set函数和欧拉方程组耦合，求解包含比热比和材料参数的耦合形式的流体力学方程组。计算方法采用二阶精度Wave-propagation算法，通过对气体-液体两种流体的一维Riemann问题、激波和气泡相互作用以及气液两种流体Richtmyer—Meshkov界面不稳定性问题进行数值模拟，抑制了不同介质界面两侧的非物理振荡。