一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

热平衡等离子体的二维复合流动数值计算

采用中性气体流动欧拉方程的近似因子分解法，耦合等离子体的热平衡模型方程，对用微波产生的等离子体在喷管内二维轴对称复合流动进行了数值计算。计算结果表明：微波等离子体电子数密度不超出临界值，1kW以下功率微波产生的是弱电离等离子体；气流的马赫数分布不受微波加热功率、喷管入口压强的影响，喷管内等离子体中电子数密度随微波功率的加大而增加，随入口压强的增加而减小。     

NEW TIME-EXPLICIT ASYMPTOTIC METHOD FOR THE SOLUTION OF AN EXACT CONTROLLABILITY PROBLEM OF SCATTERING WAVES

In this paper a new time-explicit asymptotic method is presented through introducing an auxiliary parameter for the solution of an exact controllability problem of scattering waves. Based on an exact control function, the asymptotic iteration is controlled by the auxiliary parameter and the optimal auxiliary parameter is updated during the iteration based on the existing or current iterated solutions of the wave equation. The numerical results show that the new method presented has a significant advantage over the purely asymptotic method in the history of convergence and has the ability to solve the scattering by the multi-bodies.     

基于N-S方程的地面效应模拟

模拟了低速情况下二维翼型的地面效应问题，其中翼面流动利用一种改进的、基于不可压N－S方程的离散涡方法来模拟，地面模拟则利用一系列连续分布的涡板块来实现。利用上述模型研究了地面对翼型表面流动和压力分布的影响，并对有地面存在条件下的尾涡发展规律进行了讨论。计算结果表明，该方法能够有效地模拟低速情况下翼型的近地面流动，具有一定的工程应用价值。     

直升机非线性运动方程及其数值分析

基于牛顿法建立了单旋翼带尾桨直升机的非线性运动方程，提出了一种求解直升机非线性运动方程的数值方法。建立方程时，用欧拉角描述直升机在空中的姿态，旋翼采用有挥舞铰且在铰上带有弹性约束的模型，桨盘处的入流采用线性模型，通过对一算例直升机的操纵响应进行数值仿真，将非线性与小扰动线性化的直升机运动方程进行了对比分析。结果表明：从操纵输入到随后的2s左右，非线性和小扰动线性化的直升机运动方程的计算结果非常接近，但在此之后，非线性运动方程的计算结果更能反映直升机的运动规律。     

绕大攻角翼型低速非定常流动的SUPG有限元数值模拟

采用既有良好稳定性又有较高精度的ＳＵＰＧ有限元方法,通过数值求解涡量一流函数形式的不可压缩流ＮＳ方程组,模拟大攻角翼型在高Ｒｅ数下几种典型情况的非定常绕流流场,同时验证了ＳＵＰＧ有限元方法在求解对流占优的二维非定常对流的－扩散方程时的有铲性。     

高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

元件疲劳可靠性估算的剩余强度模型

 提出了一个估算元件疲劳可靠性的剩余强度模型,讨论了剩余强度模型的一般形式,给出了元件在常幅和变幅载荷作用下的剩余强度的退化模型。在此基础上系统阐述了元件在常幅和多级载荷下的疲劳可靠性分析方法。算例的分析结果与试验结果符合很好,表明所提出的剩余强度模型是合理可行的。     

基于N-S方程的翼型气动优化设计

本文简要回顾了气动优化设计的最新发展，并采用连续伴随方法对粘性条件下的翼型气动外形进行了优化设计。设计过程中，采用B样条曲线的控制顶点作为设计变量，通过求解伴随方程得到目标函数关系设计变量的梯度，最后用最速下降法和BFGS变尺度法对二维翼型进行了外形优化，验证了本方法直接、快速的设计能力。     

基于分区的隐式求解二维不可压NS方程的并行实现

本文通过对目前常用并行算法的分析,基于区域分解的思想,给出了适用于隐式求解原始变量二维不可压NS方程压力修正方法的并行算法,并在网络异构编程PVM环境下进行了数值验证,讨论了计算结果,结果显示对任务级的并行,并行效率也同样随粒度的增加而提高.另外此方法还可推广到三维复杂流场的计算问题.