跨声速欧拉方程的一种隐式有限元解法

本文用有限体积Ｇａｌｅｒｋｉｎ法求解了跨声速二维欧拉方程。求通量运用了Ｏｓｈｅｒ格式，通过外插把原来只有一阶精度的逆风格式改进成二阶精度。并采用基于一性化流量的隐式格式和非结构化网格，使程序具有较高的收敛效率。     

定常可压缩无粘流非线性方程解的完全边界积分表示式

 本文以ρV为因变量将可压缩等熵无粘流非线性方程的解表示成完全边界积分形式,从而将问题的计算维数降低一维。     

熵变量欧拉方程组跨超声速流SUPG有限元数值解

通过变量代换将拟线性守恒变量欧拉方程组转换为具有对称性的熵变量欧拉方程组，对该方程组构造SUPG的弱解式，其中权函数摄动项中包含流向算子和间断捕获算子。对翼型的跨超声速流动算例表明，本法具有较好的稳定性和能有效地抑制激波两侧振荡。应用带预处理的GMRES迭代法求解非对称方程组能显著地提高收敛速度。     

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

二维可压无粘流的自适应流量修正有限元解

从非定常Ｅｕｌｅｒ方程组出发,采用非结构三角形网格,研究二维无粘流动问题的自适应流量修正有限元解。将ＦＥＭ－ＦＣＴ与网格自适应相结合计算了ＮＡＣＡ００１２翼型的三个算例。计算结果与实验及其他计算结果比较,亚音速和超音速绕流符合很好;跨音速绕流情况存在一些差别。     

具有轴向流的可压缩非定常旋涡流

本文研究埋定常轴对称自由旋涡流动结构的可压缩效应，探索旋涡的密度变化，粘性扩散和热传导的耦合作用机制。在回顾了前人某些典型结果之后，首次给邮了具有轴向拉伸率的可压缩非定常轴对称自由旋涡在高Ｒｅ数和小Ｍ数下的近似解析解。     

钝锥超声速粘性绕流的隐式推进求解方法研究

本文提出了一种新的推进求解方法。通过对模型方程的稳定性分析和钝锥超声速绕流的数值模拟,表明该方法既克服了文献中推进格式步长有下界和稳定性差的不足,也克服了亚声速区压力不好处理的困难。文中建立的一阶、二阶隐式推进格式,稳定性好,推进步长无下界,计算结果的精度较高。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

可压缩非平行流边界层稳定性研究

 采用线性稳定性理论和多重尺度方法,研究三维可压缩的非平行流边界层稳定性问题。分析了可解条件的特征,导出精确计算所需的伴随问题方程渐近外边界条件的矩阵表达式,给出有控制的重正化方法,以有效地克服刚性方程在积分求解中的困难。探讨与非平行性作用相关的方程和影响因素,特别是新的特征函数畸变对扰动增加率的作用。通过算例,清楚地显示了流动的非平行性对边界层稳定性的影响。     

二维可压粘流N—S方程的边界元解法

本文给出了一种求解二维可压Navier-Stokes(简称N-S)方程的边界元法。使用线化技术把控制方程变成变系数的线性偏微分方程,其基本解由组合代数的方法构造,从而给出了N-S方程的积分表达式。本文计算了几例低R_e数的跨声速翼型绕流,结果表明边界元法可以应用到可压粘流中。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

欧拉方程隐式迭代矢通量分裂跨音流场计算

 本文用矢通量分裂格式的有限体积法求解时间相关欧拉方程组。文中证明,矢通量分裂的有限体积法不能再采用物面反射边界条件。二维问题用隐式迭代法解代数方程组,第三维方向上则用强隐格式往返扫描。此法CFL数可达100的数量级,收敛快,占内存少。跨音流场计算结果与实验符合良好,激波只占一两个网格。     

数值求解可压缩N—S方程的增量迭代法

本文在Beam-Warming格式的基础上,引入通量分裂的概念,应用行高斯迭代技术提出了一个二阶隐式的求解BFC坐标变换下的N-S方程的增量迭代法。它消除了因式分解带来的寄生误差,合理地处理了控制方程中的粘性项,具有较好的收敛性和稳定性。二维Poissuill流动和Couette流动的计算结果和精确解吻合良好。S形进气道的数值解反映出的现象和实验结果是一致的。     

任意截面细长弯曲管道内的可压缩势流解

利用摄动法给出了任意截面细长弯曲管道内的可压缩定常势流的解。引入了一个表示管道细长程度的小参数ε,并将流场的解对小参数进行渐近展开。文中给出了对密度和轴向速度具有两项展开式、对横向速度具有一项展开式的结果。理论分析表明,在细长管道的情况下,问题的解可降低一维。同时,在上述近似条件下,求横向速度分量相当于求解一个平面Laplace方程的Neumann边值问题。     

三维机翼可压缩层湍流边界层计算

 从三维时间相依可压缩边界层动量积分方程和平均流动能积分方程出发,将二维可压缩层、湍流边界层积分方程算法推广到求解有限翼展后掠机翼的边界层流动。利用四步Runge-Kutta时间步进方案数值求解积分方程组,并利用当地时间步长加速迭代收敛。文中分析了数值方法的稳定性与收敛性,并考查了横向耗散项对计算结果的影响。算例表明,能获得令人满意的三维机翼定常可压缩层、湍流边界层的计算结果。