三维化学非平衡无粘流场的数值模拟

本文用ＮＮＤ格式求解了化学非平衡流场，利用点隐式处理组分连续方程的刚性，采用全流场激波捕捉办法和特征理论，成功地进行了物面边界条件的处理，模拟了高马赫数来流条件下捻激波流场。化学反应采用了五组分五反应模型，建立起来的三维流场计算程序可用于复杂外形化学非平衡无粘流场的数值模拟。     

一次强对流天气过程伴随云场的数值模拟研究

本文利用中尺度数值预报模式MM5对2005年6月25—26日的一次由中尺度对流系统（MCS）引起的暴雨过程进行数值模拟研究，分析表明：MM5很好地模拟了此次MCS过程及伴随的云物理过程。地面强降水中心位置与云水含量高值中心相吻合，云场的分布及移动与地面降水分布一致：雨水含量高值区的分布与其上空的云水含量分布相似，高值中心的出现与云水向雨水转化的时间一致；云水向雨水的转化，加大了低层的雨水含量，对降水的组织和维持有重要作用。     

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

弹体穿越平面激波三维非定常流场的数值模拟

用数值方法求解三维非定常可压缩Navier-Stokes方程，计算了复杂弹体在相同马赫数、不同迎角下的流场和气动力特征，计算结果与实验结果吻合，证明本文采用的方法能够准确计算此类复杂流场；然后采用该方法模拟了弹体穿越正激波的非定常过程，得到了弹体周围流场及气动力特征随穿越过程的非定常变化，阻力特征穿越后急剧下降。     

高超声速喷管三维热化学非平衡流场的数值模拟

利用数值求解三维热化学非平衡Navier-Stokes方程的方法和多块网格技术,数值模拟了典型高焓激波风洞锥型喷管的三维非平衡流场,分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、组分N2和O2振动温度、马赫数和组分质量分数在轴对称线上的分布,也给出了喷管出口附近的速度和温度沿径向分布等结果,获得了喷管流场的详细信息,并与已经存在的采用轴对称Navier-Stokes方程计算的结果,进行了比较,二者吻合得很好.计算结果表明:目前的代码是可以模拟多块结构网格下的热化学非平衡流动的.     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

叶栅式反推力装置开启过程的三维非稳态数值模拟与分析

针对涵道比为8的涡扇发动机叶栅式反推力装置,计算分析了反推力装置运动部件在不同运动控制规律下的开启过程对外涵流场、风扇背压、阻流门受力等影响。结果表明:在反推力装置开启过程总时间一定时,随着阻流门开始旋转的时间点向后推移,风扇背压的脉动强度增大,而阻流门受到的气动载荷会减小,存在折中的阻流门开始旋转时间点,即移动外罩开启1/3后阻流门开始旋转;开启反推力装置总时间变化对风扇背压脉动强度和阻流门受力的影响较小;紧急停飞状态下开启反推力装置,风扇背压脉动强度最大值达到20%,超过允许值,而阻流门所受到的最大气动载荷达到4500N,相当于正常开启反推力状态下的4倍以上。      

可压缩后台阶流的数值模拟

基于可压缩Navier-Stokes方程,采用Beam-Warming矢通量分裂格式以及vanLeer的矢通量分裂方法,分别通过Johnson-King及Baldwin-Lomax湍流模型对亚音速和超音速的后台阶流进行了计算,计算结果与实验符合较好。     

自适应网格局部加密方法在可压缩流场模拟中的应用

本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除无需再加密的旧细网格,以最小的运算量达到给定的精度要求。网格可一层层加密下去,按层次覆盖,每个网格是有任意方向的矩形均匀网格。这套算法独立于求解所用的差分格式,很容易和各种格式结合。 我们从M.J.Berger和J.Oliger程序出发,参考M.J.Berger和P.Colella,对守恒型方程组,实现了粗细网格交接面上数值通量守恒,从而能计算间断解。对网格边界处理、误差估计和网格重新生成等做了一些改进,并且可以灵活地在任何感兴趣区域指定加密。 我们用不定常Euler方程组计算了带前台阶的二维管道中Mach-3流动,采用MacCormack显式附加人工粘性差分格式。共用了四层网格结构:Δ=1/10、1/20、1/80和1/320。其中基本网格是Δ=1/10;到台阶上第一次激波反射区域从Δ=1/20加密至Δ=1/80,计算到定常所需运算量仅是同等均匀网格△=1/80的28％;我们还用Δ=1/320网格对台阶拐角和切向间断处指定加密,结果表明可明显改善台阶上的激波Mach反射,切向间断清晰可辨,所需的计算量是上述均匀网格Δ=1/     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

注塑模熔体非牛顿流动的数值模拟

研究注塑模流动的数值模拟。用有限元-有限差分方法模拟非牛顿流体的压力场和温度场,用流体网格法（FAN）模拟熔体锋面的推进过程,这种方法具有通用性好和容易实现自动推进的优点。在486微机上对几个算例进行了模拟计算。比较了牛顿和非牛顿流体的流动特点和锋面推进自动处理的效果,得到了满意的结果。     

刷式封严流动和换热的数值模拟

对刷式封严的流动和换热进行了数值模拟。计算中对结构复杂的刷丝部分采用多孔介质处理,用Ｄａｒｃｙ定律取代动量方程,能量方程采用改进形式,而对其它部分采用简单的矩形网格进行完全Ｎ－Ｓ求解。Ｄａｒｃｙ定律的扩散系数由所提出的实验方法确定。典型的计算结果显示所提方法是正确、合理而可行的。     

铸件定向凝固过程的数值模拟

 本文探讨大温度梯度下铸件定向凝固过程数值模拟的一些规律性。 为了提高计算精度,讨论了在大温度梯度条件下热传导方程中项对计算结果的影响。对如何处理凝固潜热,铸件-铸型界面热性质和初始条件等问题也进行了分析。 对同一铸件,在初始条件、边界条件和热物理参数相同的情况下,采用有限差分和有限元两种方法进行计算,并比较其优缺点。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

导弹弹道仿真的模型建立

运用MATLAB的Simulink工具箱建立起了导弹三通道弹道仿真模型,并给出一定条件下的弹道仿真结果,以此证明模型建立的正确性、实用性及兼容性.     

可压流双方程紊流模型及其应用

 本文提出了可压流紊流模型。所用的方法类似于K-KL双方程模型,但考虑了气流压缩性的影响。方程中的常数已由实验确定。用本模型计算了管内超音速混合流动,其理论计算曲线与实验数据相当符合,证明该模型较不可压紊流模型更为适用。     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。