基于人工粘性的间断Galerkin有限元方法

发展了一种基于人工粘性的间断Galerkin有限元方法,作出改进以增强适应性,并向非均匀网格推广.选取了典型算例对方法进行验证.一维激波管算例表明,改进的方法在保持计算结果高分辨率的同时,能够更好地抑制非物理振荡.分析得知,当人工粘性方法用于二阶DG格式时,所得计算结果的数值耗散较大,而当格式精度大于二阶时,采用人工粘性方法所得的结果的分辨率较高.通过计算圆柱高超声速粘性绕流,将三阶DG格式与三阶MUSCL格式和五阶WENO格式的结果进行对比,结果表明,该人工粘性方法对于高超声速流动计算也具有一定的优势.     

预处理AUSMDV格式在微型扑翼绕流数值模拟中的应用

由于微型扑翼的低速低雷诺特性,对其粘性绕流的数值模拟存在较大困难。传统的中心空间离散格式不具有格式耗散而需要加入带有自由度的人工耗散,给计算带来不确定因素;此外,在低速低雷诺条件下,RANS方程组的条件数太大,导致计算收敛缓慢甚至不收敛,需要引入预处理方法加速收敛。本文使用具有低耗散特性的AUSMDV空间离散格式结合预处理方法对扑翼的非定常粘性绕流进行了数值模拟,非定常推进采用了含双时间迭代的LU-SGS隐式时间推进,湍流模型采用了BL模型。计算结果表明,加入预处理能够克服AUSMDV格式难以收敛的缺点,并且与中心格式相比,基于预处理的AUSMDV格式能够更准确地模拟扑翼的非定常粘性绕流,通过将本文计算结果与文献结果对比验证了本文方法的有效性。在此基础上研究了扑翼相关参数对气动特性的影响,对进一步了解扑翼气动力产生机理,指导扑翼的气动设计具有参考价值。     

一类中心型三阶格式及应用

通过分析格式的耗散关系和重线性化技术,发展相应格式构造技术,提出了一类中心型的三阶非线性格式,格式具有与NND格式相当的网格模板.由于格式的线性形式具有很好的耗散特性,同时重线性化技术的使用使格式在光滑区尽量采用线性形式,从而理论上提升了计算对旋涡等结构的分辨和保持能力;另一方面,格式的非线性机制使格式在非光滑区退化为类NND格式,实现耗散和稳定机制.验证计算表明,格式能够开展粘性和旋涡及脉动流场的计算应用,算法形式简单,便于工程应用,容易进行有限体积推广.     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

样条Lax-Wendroff格式和样条跳蛙格式

1.引言 用三次样条求解粘性流体力学问题的数值方法已有了相当的发展。为了探讨用该方法求解无粘流动问题的可能性,本文尝试把差分方法中一些典型的计算格式扩展到样条方法,提出了样条Lax-Wendroff格式和样条跳蛙格式。以一维模型方程为例,按照Warming等的方法导出了它们的修正方程,分析了格式的稳定条件及数值弥散和耗散效应。然后扩展到对一般非线性双曲方程求解的研究。最后,讨论了某些数值解例。     

基于N-S方程的剪刀式尾桨前飞状态气动力计算研究

建立了一个适用于计算剪刀式尾桨前飞状态气动力的数值方法。在该方法中,采用惯性坐标系下的非定常Navier-Stokes方程作为主控方程,使用低数值耗散的Roe格式进行空间离散,并应用高效的隐式LU-SGS格式进行时间推进;采用双时间法以模拟尾桨前飞时的非定常现象;使用运动嵌套网格方法以计入尾桨桨叶的旋转运动。应用建立的方法,首先对有试验数据可作对比的Lynx直升机尾桨悬停状态气动性能进行了计算,验证了所建立方法对尾桨气动性能计算的有效性;同时,针对Helishape 7A旋翼进行了前飞状态气动载荷计算,验证了该方法对旋翼(尾桨)前飞气动载荷计算的有效性。最后,对剪刀式尾桨前飞状态进行了气动力计算,并分析了构型参数对剪刀式尾桨气动力的影响。计算结果表明,前飞状态时,相同总距下剪刀式尾桨的平均拉力要大于常规尾桨,且剪刀式尾桨在90°-180°-270°方位角内的桨-涡干扰现象较常规尾桨严重。     

增升状态下的全机位流数值模拟

研究了势流理论下的全机增升装置复杂外形气动力数值计算方法,该方法以低阶面元法为基础,具有计算量小、速度快的优点.复杂流动的粘性效应主要通过附面层粘性修正技术加以体现.采用该方法对某增升装置翼身组合体外形压力分布特性进行计算,并和风洞实验数据相比较,验证了方法的正确性和实用性.在此基础上对某型飞机真实复杂外形的压力分布特性进行了计算,为复杂外形飞行器增升装置的初步设计、方案选型与评估提供了非常理想的工程技术手段.     

微型飞行器小展弦比机翼的低雷诺数气动力特性分析

以固定机翼式微型飞行器的气动设计为研究背景,对小展弦比薄机翼的低雷诺数气动力特性进行了数值分析.不可压粘性绕流的求解采用了人工压缩性方法,其中对流项的离散应用了三阶迎风格式.湍流模型为Baldwin-Barth一方程模型.一无弯度薄翼的数值计算结果表明本文的数值方法是成功的.通过数值分析,揭示了弯度和展弦比对气动力特性的影响,也给出了不同弯度机翼流动结构的差别.     

适用于气动学问题的DRP类有限差分方法

研究了一类低耗散低色散的高阶精度有限差分方法,目的是直接计算非定常欧拉方程用于气动声学问题.采用的数值方法是空间四阶、时间三阶精度的色散关系保持(DRP:dispersion-relation-preserving)类有限差分格式,通过波动方程算例验证了格式模拟波动问题的能力;采用加变幅值的高波数人工耗散项来抑制高频数值振荡,使得该格式可推广用于含激波、初始间断或非线性波动的问题.通过含激波、初始间断条件的算例验证了所用变幅值人工耗散的有效性,最后作为初步应用计算了二维亚、超音速均匀流中简单点声源辐射问题,得到了很好的结果,表明此类方法在气动声学问题计算上深有潜力.     

求解双曲型守恒律的五阶松弛格式

给出了一种求解一维双曲型守恒律的五阶松弛格式.该格式以五阶WENO重构和显隐式Runge-Kutta方法为基础.本文格式保持了松弛格式简单的优点,即不用Riemann解算器和计算非线性通量函数的雅可比矩阵.用该格式对一维Euler方程进行了数值试验,并与三阶和四阶松弛格式的计算结果进行了比较,结果表明本文的格式具有更低的数值耗散和更高的分辨率.     

Relativistic MHD Simulations of Jets with Toroidal Magnetic Fields

This paper presents an application of the recent relativistic HLLC approximate Riemann solver by Mignone & Bodo to magnetized flows with vanishing normal component of the magnetic field. The numerical scheme is validated in two dimensions by investigating the propagation of axisymmetric jets with toroidal magnetic fields. The selected jet models show that the HLLC solver yields sharper resolution of contact and shear waves and better convergence properties over the traditional HLL approach.     

基于伴随算子的气动布局优化技术及其在大飞机机翼减阻中的应用

为了提高大型飞机机翼的气动性能,发展了基于Bezier-Bernstein曲线的机翼气动布局参数化技术和网格变形技术,以及基于Navier-Stokes方程的流场解算器、伴随算子解算器和二次规划寻优算法,形成了气动布局优化软件.通过改变约束条件、设计参数范围和初始点,对机翼进行了一系列的极多参数和多约束的气动布局精细优化设计,得到了良好的优化结果,优化结果表明该优化方法是高效和可靠的.      

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

NND格式的推广及在粘流计算中的应用

本文把ＮＮＤ格式推广到了特征变量，守恒变量和原始变量形式。通过引入通量差分分裂方法，由守恒变量形式的ＮＮＤ格式得到了Ｒｏｅ的近拟Ｒｉｅｍａｎｎ解法，并讨论了限制器的作用。在隐式ＮＮＤ格式中引入了一种分解算法，提高了计算效率。本文的数值计算表明，ＮＮＤ格式对复杂波系干扰，边界层分离与再附等流动现象均有较高的分辨率。     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

A hybrid multidimensional Riemann solver to couple self-similar method with MULTV method for complex flows

Since proposed, the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows. However, it needs numerous logical operations in supersonic cases, which limit the method’s applicability in engineering problems greatly. In order to overcome this defect, a hybrid multidimensional Riemann solver, called HMTHS (Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures), is proposed. It simulates the strongly interacting zone by adopting the MHLLES (Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures) scheme at subsonic speeds, which is with a high resolution by considering the second moment in the similarity variables. Also, it adopts the MULTV (Multidimensional Toro and Vasquez) scheme, which is with a high resolution in capturing discontinuities, to simulate the flux at supersonic speeds. Systematic numerical experiments, including both one-dimensional cases and two-dimensional cases, are conducted. One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves, shock waves, and linear contact discontinuities. Two-dimensional cases, such as the double Mach reflection case, the supersonic shock / boundary layer interaction case, the hypersonic flow over the cylinder case, and the hypersonic viscous flow over the double-ellipsoid case, indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows. Therefore, it is promising to be widely applied in both scholar and engineering areas.     

高亚声速二维位势流动的气动力计算

本文提供了一个精确表达式,通过它可将亚声速位势流动仿射变换成不可压缩流动。同时给出了新的气动力计算公式。在流过Joukovski翼型的气动力计算中,它与用计算机提供的精确数值计算结果相比较,二者很符合。     

MLP高阶重构格式应用

在有限体积框架下,利用MLP(multi-dimensional limiting process)系列重构格式结合HLLHLLC(Harten-Lax-van Leer with contact)近似黎曼求解器,同时引入激波探测函数进一步降低MLP在光滑流动区域的数值耗散,数值模拟了超声速前台阶流动、结冰翼型的非定常流动、高超声速双楔流动和DLR F6-WB跨声速流场,研究了MLP系列格式在可压缩复杂流场中的表现.结果表明:在多维空间中,MLP格式能够在如强斜激波与网格线不重合等复杂流场数值模拟中保持严格的流场单调性;具有和传统MUSCL(monotone upstream-centered schemes for conservation laws)格式类似的计算效率,可以实现5阶,甚至更高阶重构;数值耗散更低,捕获更准确的激波位置,对航空工程数值模拟具有重要意义.     

LENS风洞试验返回器模型气动热特性模拟

国家数值风洞高超声速流动模拟软件HyFLOW的研制对打破国外同类软件的技术壁垒具有重要意义。与国外DPLR软件进行了对比研究，同时系统介绍了HyFLOW软件求解器的数值方法、物理化学模型以及壁面催化特性计算模型等主要方法，采用典型算例对有限催化模型进行了数值验证，最后基于LENS风洞试验146 mm返回器标模外形开展了高超声速气动热特性数值模拟。研究结果表明，HyFLOW软件在高超声速热化学非平衡流动模拟与评估方面的气动力计算精度高，与国外同类软件DPLR相当，同时其壁面催化条件下的气动热计算精度可靠，可信度高。