一种求解欧拉方程改进的通量分裂方法

本文以Steger-Warming提出求解Euler方程的通量分裂方法为基础,利用改进的特征值分裂方法,建立了一种连续的通量分裂方法。该方法克服了Steger-Warming通量分裂法存在的分裂后的通量在个别点导数不连续的问题,同时方法简单,便于推广。对于通量分裂后的Euler方程,采用迎风格式,建立了相应的二阶隐式有限体积方程。为了减少求解所得块矩阵方程组所需要的工作量,本文利用特征变换,在不改变原方程组离散化精度的前提下,将所得块矩阵方程组化为分离的代数方程组进行求解,从而大大减少了计算量。文中分别对一维管流,二维绕NACA-0012和NACA64A-10翼型在亚声速和超声速来流的跨声速流场进行了计算,并将所得结果与有关实验和计算结果进行了比较,均符合良好。 本文所建立的方法具有占机内存小,计算时间短的优点。在普通中小型计算机上亦可使用。     

多级风扇三维粘性流场的数值模拟

 进行了多级跨音速风扇粘性流动的数值模拟研究。求解全三维N-S方程,通过高分辨率格式提高跨音速流场计算的准确性,利用LU-SGS隐式解法加快计算收敛速度。给出了一双级高负荷跨音速风扇设计和非设计工况下的详细计算结果,各工况下的计算与实验数据均吻合较好。     

多体系统的计算机自动建模方法

文章基于D’Alembert虚功原理，系统地论述了计算机自动建立链状多体系统动力学方程的方法，并给出了在系统中加入约束条件的矩阵方法，因而使所推导的计算机自动建模方法也适用于约束多体链系统。     

非惯性非正交曲线坐标系中的边界层方程

应用张量分析方法将 N-S方程写于作非惯性运动的非正交曲线坐标系中 ,然后用匹配渐近展开法导出此坐标系中的边界层方程。消去压强项后 ,与惯性坐标系中的边界层方程相比 ,非惯性系中的边界层方程在 x1和 x2 方向的动量方程中分别只多了一哥氏力项 ,它们是由物体转动角速度的物面法向分量产生的 ,其它惯性力项均不进入边界层方程。方程中所含曲线坐标系的拉梅系数及其导数只需在物面上取值 ,物面曲率影响不进入边界层方程     

高应变率条件下38CrMoAl钢的动态力学行为及失效模型

38CrMoAl钢是飞机某频繁经受冲击载荷结构的主要材料,为了获得该钢在冲击载荷下的力学行为和失效参数,本文对38CrMoAl钢开展了准静态试验、应变率从650 s^-1至5 500 s^-1的动态力学试验以及应力三轴度试验,对材料的金相组织和断口形貌进行了微观观察分析,最后基于试验结果确定了材料的修正强度模型和失效模型参数。试验结果表明,38CrMoAl钢具有较强的应变率敏感性,材料屈服强度随着应变率的增加具有正应变率强化效应,应变率在10^-3~5 500 s^-1的范围变化时,屈服应力从450 MPa上升到了1 085 MPa,应变率为5 500 s^-1时材料的屈服应力达到了准静态条件下的2.41倍。对材料组织的微观观察分析发现材料主要由呈现颗粒状的回火索氏体组成,材料在高应变速率加载下其塑性发生了下降,断口由韧窝型塑性断口向局部解理断面发生转变,韧窝尺寸也随着应变率的升高而变小。对Johnson-Cook（J-C）强度模型中的应变率强化项进行了修正,使修正后的本构模型可以更有效地描述38CrMoAl钢的动态力学行为,最后结合数值模拟方法确定了材料的J-C失效准则参数。     

固体火箭发动机燃气射流流场和声场数值计算

燃气射流噪声是固体火箭发动机工作过程中的主要噪声源之一,射流流场的参数对其产生的射流噪声有重要影响。通过大涡模拟(LES)对不同尺寸喷管形成的超声速高温射流进行了三维非稳态数值模拟,随后在合适的声源面中,采用FW-H(Ffowcs Williams-Hawkings)声学模型和傅里叶变换得到了燃气射流噪声声压级的空间分布。计算结果表明,随着喷管尺寸增大,超声速射流核心区变大,喷口流场波节数增加,对喷管尾流场的影响域扩大,其产生的射流噪声也增强;燃气射流噪声辐射有较强的指向性,在射流轴向30°角方向噪声声压级最大,与相关文献中的试验结果比较吻合。研究结果可为后续固体火箭发动机降噪设计提供参考。     

抛物化 Navier-Stokes 方程(PNS)和高氏 PNS 理论及其应用的评述

抛物化 NS 方程得到广泛应用,已经成为工业标准气动计算的基础。现有的八种抛物化 NS 方程有不同的名称,方程中粘性项的形式略有不同,其中的 PNS 和薄层（TL）NS 方程应用最多。但是这些方程都具有类似的数学性质,例如,当流向方向上马赫数大于1时,他们都是抛物型方程,可以采用空间推进算法（SMA）进行求解。与采用时间推进算法求解的 NS 方程或雷诺平均（RA）NS 方程相比,PNS-SMA 计算降低了空间的维数,节省了大量的存储空间和 CPU 计算时间。PNS-SMA 算法也获得了巨大的进展。但是,早期 PNS 研究在理论上是相当模糊的,高智在1990年提出的粘性／无粘干扰剪切流理论（ISF）弥补了这一不足。ISF 理论概括了 PNS 方程所能描述的基本流动,提出了其流动的运动规律及数学定义式,所导出的 ISF 方程组也属于 PNS 方程的一种。为了不增加新的名称,我们将 ISF 方程组也称为高氏 PNS 理论和方程组。这一理论在 NS 方程和RANS 方程的计算中均有重要的应用。例如,计算最优坐标系的选择以减少伪扩散,网格尺度选择及局部网格加密设计以捕捉高超声速流动中物体表面热流等的急剧变化,壁面压力边界条件的选择以及由高 PNS 导出的壁面判据来进行 NS 和 RANS 近壁数值解可信度评估。本文评述了一些初步的应用,进一步的应用和综合 PNS-SMA,RANS-SMA 以及 PSE-SMA 计算值得深入研究,这里 PSE指抛物化稳定性方程。     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

数值模拟二维跨音速叶栅可压粘性流动

用有限体积法数值模拟了任意回转流面跨音速压气机叶栅中的粘性流动，四步龙格－库塔法用于N－S方程的时间推进。给出了一种抑制叶栅计算中“数值失速”的方法，并采用隐式残差光顺技术加速收敛。用所发展的方法计算了三种叶栅在宽广来流Ma数条件下的流动，计算结果与实验结果吻合良好，流场激波具有较高的分辨率。     

飞机结构疲劳试验载荷的优化设计

提出用对载荷误差求极值的方法实现最小误差控制的试验载荷计算,并通过精度加权系数综合控制误差的合理分配,完成试验载荷的优化设计。该方法已经应用于某型飞机疲劳试验的试验载荷计算,为该机全尺寸疲劳试验的成功奠定了基础。