改进的五阶WENO-Z+格式

WENO-ZWENO-Z$+\!$格式的性能提升依赖于新增项的作用,该项的作用是在WENO-Z格式的基础上进一步增大欠光滑子模板上的权重. 系数$\lambda$被设置用来调控该项的作用, 以避免负耗散. 本文指出了WENO-Z$+\!$格式的缺陷,其所采用$\lambda $的取值方式既不能充分发挥格式的潜力, 也未完全消除负耗散;提出$\lambda $的值应随当地流场数据变化,方能充分发挥新增项在降低数值耗散、提高分辨率上的潜力. 基于此,本文重新设计了$\lambda $的计算公式,该公式能自适应地调控新增项的作用: 只在欠光滑子模板上的权重容易过度增大的地方削弱该项的作用,以避免负耗散; 在其他地方则充分发挥新增项的作用,最大限度增大欠光滑子模板上的权重, 提高格式的分辨率.将使用该系数公式的新格式命名为WENO-Z++, 并对其数值性能进行了系统的研究.理论分析表明, 新格式在间断处具有基本无振荡(essentially non-oscillatory,ENO)特性和更低的数值耗散. 对近似色散关系(approximate dispersion relation,ADR)的研究表明,新格式有效地避免了因过度增大欠光滑子模板上的权重而带来的负耗散,其频谱特性也得到了显著提升.本文还推导了使新格式在极值点处也能保持最优阶的精度的参数设置.一系列求解Euler方程的数值试验表明,新格式的激波捕捉能力和对复杂流场结构的分辨率都显著优于原WENO-Z$+\!$格式.}     

两种矩形单元的Babuka-Brezzi条件

§1.引言 用原始变量有限元方法求解不可压缩流体流动方程组,要求速度有限元子空间V_h和压力有限元子空间Q_h 满足Babuska-Brezzi条件:     

二维超声速混压式进气道的数值模拟

本文采用TVD格式求解二维可压非定常Euler方程组,对二维超声速混压式进气道进行了数值计算,首先模拟了后部压力对进气道的影响过程,当进气道的后部压力与入口压力的比值逐渐增加时,进气道内的正激波逐渐往前移动,当比值达到一个临界值,进气道将不起动.其次.通过对R52.1和R54.5两种进气道模型进行计算,对它们的总压恢复进行了对比.通过计算得到了喉道处曲线弧的曲率越大,进气道的总压损失越大,进气道也越难起动.     

高超声速非定常流动的数值模拟与气动热计算

高超声速飞行器研究中的一个重点问题是飞行器表面的气动加热,它对飞行器的气动、热特性及安全性有重要的影响.受到当前实验技术的限制,地面实验无法准确模拟真实飞行条件,所以采用数值模拟研究气动加热问题成为目前重要的研究手段.本文采用数值方法求解三维N-S方程,得到钝头体再入模型绕流的瞬态流场,驻点温度及表面热流沿轨道变化规律.计算中采用变边界条件模拟沿轨道飞行的非定常性.     

Low-diffusion preconditioning scheme for numerical simulation of low-speed flows past airfoil

The preconditioning technique can address the stiffness of a low Mach number flow, while its stability is poor. Based on the conventional preconditioning method of Roe's scheme, a low-diffusion scheme is proposed. An adjustable parameter is introduced to control numerical dissipation, especially over the dissipation in the boundary layer and extremely in a low speed region. Numerical simulations of the low Mach number and low Reynolds number flows past a cylinder and the low Mach number and high Reynolds number flows past NACA0012 and S809 airfoils are performed to validate the new scheme. Results of the three tests well agree with experimental data, showing the applicability of the proposed scheme to low Mach number flow simulations.     

定常Navier-Stokes方程组罚函数方法中的惩罚因子

§1 引言 近年来人们越来越多地采用有限元方法求解Navier-Stokes方程组。在这个方程组中,不可压缩条件diy u=0是一个约束条件。本文所讨论的罚函数有限元方法是用div u=-εp来代替上述不可压缩条件,这里ε>0是一个参数,称为惩罚因子。     

类Apollo返回舱外形肩部热流的数值影响因素分析

预测类Apollo返回舱外形在高焓来流下的气动热特性,研究网格Reynolds数、壁面温度、多种化学反应模型以及限制器对预测热流的影响.采用ESI-CFD-FASTRAN软件作为数值模拟平台,使用基于温度梯度及分子扩散效应的热流模型;空间离散采用Roe-FDS格式,时间推进采用点隐式;采用等温壁面条件.数值计算表明:（1）热流在返回舱头部驻点处达到一个极值,沿着壁面热流不断下降,经过返回舱肩部热流有突越上升;（2）满足网格Reynolds数小于10的网格获得的热流较为准确;（3）使用Gupta模型计算得到的热流与Park85模型得到的类似,但是获得的热流分布类似;（4）采用湍流模型获得的头部肩部热流结果与层流结果相同;（5）二阶minmod限制器实现了高阶格式,其计算得到的热流结果在肩部略高,但是整体分布略低于不带限制器的格式.因此,在计算中采用满足网格Reynolds数壁面网格,采用带限制器的高阶格式计算获得的热流分布更加准确;由于头部热流主要贡献并非来源于湍流,因此对于肩部热流采用层流模型足够准确.      

低Mach数翼型绕流数值模拟的低耗散预处理格式

低Mach数流动中,基于可压缩流动的数值模拟算法存在严重的刚性问题,预处理方法可以有效地解决这一问题,但其计算结果不稳定．基于原有的预处理Roe格式,引入可调节参数,得到一种新的低耗散格式．该格式可以减弱边界层以及极低速区域的过度耗散,使得整个流场计算稳定．低Mach数、低Reynolds数定常圆柱绕流和低Mach数、高Reynolds数翼型(NACA0012和S809)绕流3个验证算例表明,带可调节参数的低耗散预处理方法正确可靠,是低速流动数值模拟的有效方法．