激波结构内流动问题的气体运动论描述

从介观Boltzmann速度分布函数理论出发,发展计及分子粘性碰撞截面与扩散碰撞截面,可描述各流域一维气体流动问题的Boltzmann简化速度分布函数方程及其气体运动论数值计算方法。通过对不同Knudsen数下非定常激波管流动及不同马赫数定常正激波结构问题数值模拟,研究分析不同流区的激波突跃变化过程以及近连续流、稀薄过渡流特有的分子输运现象,揭示不同马赫数、不同分子模型的激波内流动与传热变化规律,证实基于Bo-ltzmann模型方程的气体运动论数值计算方法用于激波结构内流动研究的准确可靠性。     

基于动态混合网格的鱼体巡游数值模拟

对鱼类等柔性变形体的数值模拟需要解决如下几个问题(1) 柔性体动态网格生成;(2) 基于动态网格的不可压流非定常算法.本文基于Delaunay背景网格插值方法,发展了一套快速生成柔性体变形的动态混合网格生成方法;对于不可压流问题,采用虚拟压缩方法进行求解,结合双时间步方法来处理非定常问题;最后数值模拟了鱼体巡游的非定常流动,分析了鱼体的摆动频率对力、功耗等流体动力学特性的影响.     

超声速二维湍流分离流动的数值解法

本文采用层流中引导简化N-S方程的方法,建立了简化的湍流运动方程组,分别用零方程、一方程模型计算了绕二维压缩拐角的湍流分离流动。数值求解用显隐式的反扩散格式。     

飞船返回舱俯仰振荡的动态稳定性研究

本文研究了以平衡攻角为中心作单自由度俯仰振荡的飞船返回舱 ,其动态稳定形态随来流M∞ 的变化。设θ(t)是由平衡攻角起算的俯仰振荡角 ,Cm 是作用在飞船返回舱上的气动俯仰力矩系数 ,Cμ(θ, θ)· θ是机械阻尼力矩 (自由飞行时为零 ,实验时要计入其影响 ) ,文中给出飞船返回舱在平衡攻角处的俯仰振荡动态稳定性判据 ,并证明λ =λ(M∞) = Cm θ0+Cμ(0 ,0 ) 1 - Cm ¨θ0 是决定动稳定形态的重要参数。如果随M∞ 的变化 ,λ(M∞)由λ<0经过λ =0变化到λ>0 ,则飞船返回舱将由稳定的点吸引子形态 (即稳定在平衡攻角状态 )演化为周期吸引子形态 (即作周期振荡 )。对应于λ(M∞) =0的马赫数就是飞船返回舱的俯仰运动出现Hopf分叉的临界马赫数Mcr 。本文首先分析了飞船返回舱所受动态气动俯仰力矩的依赖状态变量 ,然后应用非线性动力学理论对飞船返回舱的俯仰运动进行了定性理论分析 ;最后耦合求解俯仰振荡方程和非定常Navier Stokes方程 ,数值模拟了飞船返回舱俯仰振荡随来流马赫数变化的Hopf分叉过程 ,验证了分析结论     

求解N-S方程的一个简单隐式算法

1.引言 近十几年来,求解分离流动的N-S方程的数值解法,取得了很大进展。显式有限差分解已发展得较为成熟。但由于稳定条件的限制,获得稳定解的计算时间较长。为了克服这一困难,发展了隐式方法,如隐式因式分解法和交替方向的隐式方法。但是这些方法都要求解三对角矩阵方程,同显式方法比较,程序复杂,每步计算时间较长。为了避免三对角矩阵迭代的困难,MacCormack提出了新的隐式算法,其计算迭代过程为两对角矩阵追杆。本文在[1]的基础上,提出了一种更为简单的隐式算法,其迭代过程     

基于Boltzmann模型方程不同流区复杂三维绕流HPF并行计算

随着现代航空航天快速发展,是否能建立一套可有效模拟稀薄流到连续流各流域绕流问题统一算法,已成为工程应用部门和学术研究领域关心的问题。通过开展基于Boltzmann简化速度分布函数方程数值求解,发展起从稀薄流到连续流统一算法。借助区域分解并行化方法研究建立气体运动论统一算法并行方案,通过对统一算法进行HPF并行化程序设计及算法考验,拟定不同流域三维球体绕流及类"神舟"返回舱外形体绕流算例,进行HPF并行计算,并将计算结果与有关实验数据、DSMC模拟值进行定量比较、分析。研究表明,所发展的统一算法具有很好的并行独立性,基本达到了线性加速的并行效果,算法负载平衡和并行可扩展性较好,可望建立起新型的能可靠模拟不同流域三维绕流问题的HPF并行算法研究方向。     

三维控制翼周围粘性分离流动的数值模拟

本文首先分析了中心型隐式格式和迎风型隐式格式的优缺点,在此基础上发展了一种自由参数可确定的反扩散隐式格式。用此格式计算了绕三维控制翼的层流和湍流分离流动。出发方程是三维简化N-S方程,用Baldwin-Lomax湍流模型。结果表明,这种隐式格式效率高,能满意地刻划流场。     

二维超声速混合层流动稳定性的数值分析与并行计算

本文用DNS方法对两股空气的二维超声速混合流动(Mc=0.5)的稳定性问题进行了计算模拟.通过对采集的流动参数在时空两个方面进行的分析,发现流动失稳后出现了间歇结构.有限差分计算采用了并行计算技术,在微机并行环境下可高效率地完成大规模计算任务.     

差分计算中激波上、下游解出现波动的探讨

本文探讨了差分计算中激波上、下游解出现波动的原因。研究指出:在NS激波方程中,附加二阶耗散项,可平滑激波;但附加三阶弥散项或四阶耗散项,在一定条件下,可引起解在激波上、下游波动。因此在求解NS激波方程时,若采用二阶或三阶差分格式,由于格式弥散项和耗散项的存在,将引起解的波动。     

自由分子流气动实验的模拟理论

本文给出了自由分子流气面作用的相似准则;分析了分子束流动与无限空间自由分子来流的相似性;提出了自由分子流气动实验的模拟理论。