基于Riemann问题的解在自适应无结构网格上求解Euler方程

给出了一种根据 Riemann问题的解 ,计算网格单元边界处的守恒量通量的方法。该方法不受网格单元形状的限制 ,具有较好的通用性。实际计算表明 ,采用这一方法在自适应无结构网格上求解 Euler方程 ,对于复杂的流动细节具有较高的分辨率 ,并且使用多重网格法能够加速收敛     

圆轴承三维热流体动力润滑的研究

联立求解了广义雷诺方程、三维能量方程、三维固体热传导方程及载荷平衡方程 ,并考虑粘度随温度及压力、密度随粘度及压力的变化 ,在油膜与轴瓦界面使用热流连续性边界条件 ,得到了圆轴承油膜及轴瓦的三维温度场。结果表明最高温度并不出现在油膜与轴瓦交界面 ,不同圆周位置沿膜厚方向的最高温度出现位置不同。证实了考虑三维温度得到的轴承的润滑性能与考虑二维温度得到的轴承润滑性能有很大的差异。采用Newton-Raphson法求解广义雷诺方程和载荷平衡方程 ,采用快速扫描法求解由有限差分离散的能量方程和热传导方程 ,在计算中发现所使用的算法收敛速度快 ,大大节约了计算时间。     

一种高效的多重网格三维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有限元数值模拟

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

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

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

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

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