全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

类航天飞机前身结构与高超声速流场的耦合传热模拟分析

开展了高超声速流场与结构温度场的耦合数值计算.流场部分求解了三维非定常全Navier-Stokes(N-S)方程,空间差分采用Harten-Yee的TVD(总变差衰减)格式,时间离散采用双时间步推进.固体结构传热部分求解了三维非稳态的热传导方程.通过流固交界面,流体从固体部分得到温度边界条件,固体从流体部分得到热流边界条件,从而实现流场和固体温度场的紧耦合计算.用绕无限长圆柱的气动加热计算验证了该算法的有效性,并对类航天飞机前身结构在气动加热过程中的温度变化做了比较详细的分析.计算结果表明,固体结构在遭遇到气动加热后的一段时间内,壁面温升对壁面热流的影响是很大的,由于一体化计算能很好的综合考虑热壁的影响,因此,开展一体化计算是很有必要的.      

高超声速流场与结构温度场一体化计算方法

在对国内外流场与结构温度场一体化计算方法的不足进行细致讨论的基础上,提出了一种高超声速流场与结构温度场一体化计算方法.采用统一的积分方程组作为气动加热和结构传热物理过程的控制方程,对整个物理场进行统一的迎风格式有限体积方法离散,给出了流场与结构交界面上温度、温度梯度及导热系数等参数的计算方法.在时间推进方面,定常状态采用多步龙格库塔迭代格式,非定常状态则采用双时间步长方法.采用发展的一体化算法对二维圆管模型的气动加热和结构传热问题进行了数值模拟.结果表明:2s时驻点物面温度为390K,与其他文献的误差在3.1K范围内;稳定时驻点物面温度为647K.      

时空守恒元和解元方法的爆震波一维数值模拟

运用基元化学反应模型和流体力学的基本原理，建立了脉冲爆震发动机一维非定常反应流数学模型。应用一种改进的时空守恒元和解元(CESE)方法建立控制方程的差分格式，针对化学反应源项的刚性问题，在对非线性源项进行线性化处理的基础上，应用隐式梯形方法实现源项的积分。对氢氧爆震波的形成和传播过程进行数值计算，结果表明，应用CESE方法和基元化学反应模型能够成功地模拟氢氧爆震波的精细结构。     

双微射流作动器合成流场数值模拟

采用精度高和稳定性强的隐式高阶紧致差分格式结合Beam-Warming近似因式分解法求解全N—S方程，对二维、粘性、非定常、可压双微射流作动器合成流场进行数值模拟。内置牛顿似子迭代用来消除因近似因式分解、线性化、显式边界条件及隐式一边采用低阶空间离散所带来的误差。与其它人工粘性方法相比，隐式高阶数值过滤方法对许多情况，特别是对于马赫数非常低的流场计算有明显优越性。计算结果揭示了双微射流合成流场的特点及其产生、发展与耗散的过程。此外，双微射流合成流场的形成受双微射流作动器相位差的影响，合成射流偏向相位超前的一方。相位差愈大，偏转愈剧烈。     

固体火箭发动机复合结构喷管温度场与应力场理论计算

对某长时间工作的固体火箭发动机复合结构喷管的二维轴对称瞬态温度场及二维轴对称准静态应力场进行了理论预估。温度场的计算采用时间域上的有限差分及空间域上的有限元相结合的方法。应力场的计算采用以位移作为未知量的有限元法。计算中考虑了材料的方向性及物性参数随温度的变化，采用变带宽一维压缩存贮刚度矩阵，节约了大量的内存。计算结果与试验结果吻合较好。     

三维方柱不可压缩绕流的大涡模拟计算

用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究.求解程序采用虚拟压缩方法和双时间步的时间推进方案,对流项的差分采用高阶精度的数值差分格式.通过与实验和文献结果比较分析,验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的.     

带操纵面机翼的非定常跨声速流数值解法

本文介绍一种带操纵面机翼的非定常跨声速流的有限差分计算方法。采用的方程是三元非定常跨声速修正小扰动位势方程,使用时间积分法,求解格式是近似因式分解交替方向隐式(ADI)格式。使用这种方法计算了F-5机翼在来流马赫数为0.9和0.925时的定常气动力和操纵面振荡的非定常气动力。计算结果与国外的NLR试验结果和XTRAN3S方法的计算结果进行了比较,表明计算是成功的。     

侧向脉冲喷流瞬态干扰流场探讨

当姿态控制发动机工作时间很短时, 喷流干扰产生的非定常气动特性就显得非常重要。为了研究脉冲发动机起动到关闭过程中引起的非定常气动特性对导弹控制效应的影响, 采用了二阶时间精度的LU隐式时间推进格式和双时间计算技术, 数值求解了飞行马赫数M∞=5, 迎角α=0°情况下, 拦截导弹姿态控制脉冲发动机喷流的瞬态过程和所产生的非定常气动特性。以0 001ms间隔观察了导弹表面在喷流出口附近非定常干扰流动分离区的发展过程。给出了详细的喷流瞬态干扰流场结构, 以及喷流瞬态干扰区随时间改变的流场细节特性。研究表明: 喷流前的高压区和喷流后的低压区对非定常效应非常敏感。并且, 当喷流已完全关闭时刻还存在喷流羽流的残余干扰量。     

固体火箭发动机喷管喉部沉积过程的传热分析

推进剂含金属填加剂的固体火箭发动机在工作过程中有时会发生喷管喉部沉积现象,可能对发动机性能带来严重的影响。本文根据有关实验现象,建立了喉部沉积过程的传热模型。假设喷管喉部结构由沉积层、耐熔喉衬、绝热衬套及外壳等导热性能不同的材料组成的,其中沉积层的厚度是随时间变化的,是多层壁变边界的不稳定传热问题。根据喉部沉积的传热模型建立了偏微分方程组,采用有限差分完全隐式格式用电算机进行数值分析计算。计算分析结果给出了喉部截面温度场及其变化规律,并从理论上预示有关因素对喉部沉积的影响,与实验规律一致。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法（Ⅱ）

建立了二维振荡机翼的非定常反问题的有限差分解法,并通过若干算例对求解方法进行了验证.首先根据基于欧拉方程的非定常气动反问题数学模型确定了具体求解方案,接着重点讨论了壁面非定常反问题边界条件和远场非定常无反射边界条件的处理方法,编写了相应的计算程序,并以绕1/4弦长点做小幅俯仰周期振动的NACA-65-0012翼型为初始机翼,在不同目标压强的条件下进行了气动反问题计算得到新翼型,它们的压强分布与预设目标相吻合,这表明本文的反问题求解方案与计算方法都是可行的.     

基于时间精确自由尾迹方法的旋翼非定常 气动特性分析

 给出了一种显式修正的二阶向后差分格式(PC2B)的时间精确自由尾迹算法,用于分析旋翼非定常气动特性。由于旋翼自由尾迹控制方程具有非线性特征,求解过程中易产生数值不稳定性问题,因此本文推导了修正PC2B算法的修正方程,分析了该算法的求解精度和数值稳定性,并通过数值试验对其稳定性和计算效率加以验证及分析。基于该算法,首先对悬停和总距突增时旋翼非定常气动特性进行了计算,通过与试验结果的对比验证了模型的有效性;其次研究了总距操纵时旋翼非定常气动特性的变化规律。研究结果表明,旋翼自由尾迹几何形状的瞬态变化是造成总距操纵条件下旋翼非定常气动特性波动的主要原因。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

直升机旋翼非定常气动力及结构动力响应计算

以二维非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程为模型方程,采用Ｊａｍｅｓｏｎ有限体积并引入变系数隐式残值光顺及固接网格计算了直升机旋翼不同展向位置处翼剖面的非定常气动力,进而采用插值的方法模拟直升机在悬停及前飞状态下,旋翼的桨叶所受的非常气动力载荷,并将之与桨叶的结构振动方程相耦合,求解其时域内的结构响应,非定常下气动力计算结果与实验结果做了对比,两者比较吻合,旋翼结构响应计算结果合理。     

烧蚀端头的瞬变外形及内部温度分布

本文通过对传热、烧蚀和形状变化的耦合计算取得了各向同性和各向异性石墨再入端头帽的温度场。直接从热传导偏微分方程出发,通过座标变换,建立三层无条件稳定的显式差分格式,计算简捷。在碳表面,包括十七组元的平衡烧蚀计算表明,壁面的主要化学成分有N_2、CO、C_2N、C_3、CN等几种。对非线性形状变化方程,采用可显式计算的隐式差分格式,稳定性良好。二种石墨端头的温度场差别很大,但算例条件下的总烧蚀量相差无几。     

欧拉方程计算非定常流动的一种数值方法

提出了利用欧拉方程计算非定常流动的一种数值方法。在利用中心差分的有限体积法和双时间推进法求解非定常欧拉方程过程中，通过在固定物面边界上满足完全的非定常边界条件，计算出非定常流场。算例结果表明，本文方法具有高效、简便的特点，而且能描述激波的非定常运动和变化情况。     

高超音速流场与结构温度场耦合计算

采用松耦合方法对高超音速流动与结构温度场进行耦合计算。流场空间离散采用AUSM+格式,时间离散采用隐式LU-SGS方法。结构温度场求解采用控制容积法。分析了高超音速圆柱表面热流和结构温度场的变化。通过与实验结果对比表明了本文流固热耦合计算方法的准确性。     

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

跨声速操纵面嗡鸣的数值研究

我国的一些高速歼击机，在飞行中曾多次出现操纵面嗡鸣问题，它严重影响了飞机的性能和飞行安全，本文采用数值计算方法，对此问题做较深入的分析和研究。气动力计算采用Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，并采用ＬＵ－ＴＶＤ混合差分格式求解由此气动力解操纵面的弹性振动方程，求出现嗡鸣的Ｍ数及振动的频率、振幅等参数。计算中采用Ｂ－Ｌ代数湍流模型，网格生成采用代数方法及保角变换方法，文中讨论了隐式算法在非定常气动力计算     

固体火箭发动机喷管喉衬动边界条件下的温度场计算

本文探讨了固体火箭发动机喷管喉村在烧蚀退移、燃气中固体微粒的沉积以及沉积物吹除所引起的动边界条件下的温度场数值计算模型。模型的建立以有限体积能量平衡的基本原理为依据,采用有限差分的数值计算方法及稳定性较好的交替方向隐式格式求解。为确保数值计算结果在边界退移及网格脱落时的稳定性,采用了插值与迭代相结合的方法,以全隐式格式求解表层网格节点和受热壁面的温度。模型还考虑了多层背村材料结构、变物性材料特性等诸种因素,使之具有一定的普遍性,并可应用于工程设计和研究中。数次发动机实验证明,本文提出的温度场计算模型能够较真实地反映喉衬温度响应,具有较为满意的计算精度。