高精度格式WCNS-E-5的Fourier分析与应用

对高精度加权紧致非线性格式WCNS-E-5进行了理论分析和应用研究。首先采用Fourier方法分析了WCNS-E-5与WCNS-5及其它高阶格式(迎风偏置格式EUW-5与标准格式)的差分误差特性,研究了它们在分辨效率方面的特性和相位误差在多维方向上的各向异性特性。分析结果显示WCNS-E-5与WCNS-5在色散与耗散特性方面优于EUW-5,分辨效率也普遍高于EUW-5和格式,而且它们的相位速度在多维方向上传播所表现的各向异性特性与其它高阶格式一致。WCNS-E-5的高精度特性与WCNS-5的一致,但在计算过程中少了三对角矩阵求解因而它的效率更高,于是采用WCNS-E-5数值模拟了二维/三维复杂流场,算例包括单涡斜向运动、二维Rie-mann问题以及存在分离的大攻角钝锥流动问题。计算结果体现了WCNS-E-5具有很低的数值耗散误差,它对激波、膨胀波和接触面等间断具有很好的捕捉能力,得到的图像清晰光滑,准确再现了真实流动现象。     

尺度无关高阶WCNS格式

高速流场的数值模拟中, 既要保证对小尺度结构的高保真分辨, 又要实现对激波稳定、无振荡地捕捉.当前工程中广泛应用的高精度数值格式虽然都能一定程度地满足上述两种要求, 但仍与理想目标存在较大差距.例如, 模拟雷诺应力模型等小尺度问题时, 高精度格式在间断解附近易产生数值振荡.基于高精度格式所存在的上述问题, 本文引入去尺度函数, 探索了一种更加简单稳定的非线性权重构造方法, 并将其应用于7阶精度加权紧致非线性格式WCNS, 提出了一种尺度无关的7阶WCNS格式.该格式的性能与灵敏度参数和尺度因子的选择无关, 并且在小尺度下仍可以有效捕捉流场激波.同时, 该格式在间断处具有基本无振荡性质, 且在任意尺度函数下保持尺度无关, 并且在极值点处也能保持最优精度.本文还推导了7阶D权函数的形式.最后, 在一维线性对流方程中验证了新格式在流场光滑区能够达到设计精度, 并通过一系列数值实验证明了尺度无关的7阶WCNS格式在激波捕捉能力上具有良好表现, 为WCNS格式改进和解决可压缩湍流等非线性问题提供了一种新途径.      

浅水波方程的高精度隐式WCNS格式

设计了高精度隐式WCNS格式,用于数值求解浅水波方程。方程通量导数采用四阶隐式紧致差分格式计算;半节点通量采用五阶非线性加权插值方法计算。对于含源项浅水波方程,采用特殊的源项处理,使得隐式WCNS格式具有和谐性,即在静水条件下保持流速为零,水位为常数。通过线性算例及几个经典的无源项和含源项数值算例对WCNS格式进行测试,结果表明:该格式在光滑区具有四阶精度,在间断处具有良好的捕捉激波能力,同时对于含源项问题具有保持静水条件下定常解的和谐性。     

一种面向激波噪声计算的加权优化紧致格式及非线性效应分析

在超声速流动中, 激波与湍流结构的相互作用会产生高强度的激波噪声. 激波噪声的高保真计算要求激波捕捉格式具有高阶精度、低耗散和低色散特性, 同时还要尽可能地减弱格式的非线性效应. 现有的六阶精度迎风/对称混合型加权非线性紧致格式CCSSR-HW-6在基于对称模板构造网格中心处的数值通量时引入了两级加权, 且两级加权都需要构造非线性的权系数, 因而非线性效应较强. 本文以修正波数的误差积分函数为优化目标函数, 优化了CCSSR-HW-6格式的非线性特性, 建立了加权优化紧致格式WOCS. 精度验证表明WOCS格式的精度高于5阶. 谱特性分析表明, 与原方法相比, WOCS格式的耗散误差和非线性效应显著降低. 典型激波噪声问题数值实验表明: WOCS格式不仅提高了对高频波的分辨能力, 而且显著地消除了数值解中因格式的非线性效应所导致的非物理振荡.      

基于HLL-HLLC的高阶WENO格式及其应用研究

HLL-HLLC格式能够克服HLLC在强激波附近的激波不稳定现象,并且保持了HLLC的低耗散特性,是一种适合更大马赫数范围的近似黎曼求解器。本文从RANS方程出发,将HLL-HLLC近似黎曼求解器结合五阶WENO重构,实现了对无粘通量的高阶离散;同时,采用完全守恒形式的四阶中心差分格式处理粘性项,建立了RANS 方程的高阶数值求解格式。通过对四个经典算例,钝头体、 ONERA M6机翼、DLR F6-WB翼身组合体和DLR F6-WBNP复杂外形的数值模拟,考察了两种WENO改进格式在复杂流场中的表现,研究了高阶格式的收敛特性;给出了在复杂流动中 WENO自由参数的推荐值,以增强求解的收敛性。算例结果表明,本文构造的高阶格式鲁棒性好,能够显著改善激波位置和激波强度,捕获更丰富的流场细节,满足复杂工程应用需求。     

基于HLL-HLLC的高阶WENO格式及其应用研究

HLL-HLLC格式能够克服HLLC在强激波附近的激波不稳定现象,并且保持了HLLC的低耗散特性,是一种适合更大马赫数范围的近似黎曼求解器。本文从RANS方程出发,将HLL-HLLC近似黎曼求解器结合五阶WENO重构,实现了对无粘通量的高阶离散;同时,采用完全守恒形式的四阶中心差分格式处理粘性项,建立了RANS方程的高阶数值求解格式。通过对四个经典算例,钝头体、ONERA M6机翼、DLR F6-WB翼身组合体和DLR F6-WBNP复杂外形的数值模拟,考察了两种WENO改进格式在复杂流场中的表现,研究了高阶格式的收敛特性;给出了在复杂流动中WENO自由参数的推荐值,以增强求解的收敛性。算例结果表明,本文构造的高阶格式鲁棒性好,能够显著改善激波位置和激波强度,捕获更丰富的流场细节,满足复杂工程应用需求。     

高精度加权紧致非线性格式的研究进展

综述了高精度加权紧致非线性格式WCNS在理论分析以及复杂流动应用方面的研究进展. 首先回顾了国内外高精度格式研究的概况, 然后介绍了WCNS的研究与发展历程. 在对WCNS进行了Fourier分析和渐近稳定性分析后, 给出了WCNS求解多维复杂流动的算例.      

无粘可压缩流动的改进高精度方法

为更准确捕捉复杂流场的流动细节,通过对WENO格式的光滑因子进行改进,发展了一种新的五阶WENO格式。对三阶ENO格式进行加权可以得到五阶WENO格式,但是不同的加权处理,WENO格式在极值处保持加权基本无振荡的效果不同,本文构造了二阶精度的局部光滑因子,及不含一阶二阶导数的高阶全局光滑因子,从而实现WENO格式在极值处有五阶精度。基于改进五阶WENO格式,对一维对流方程、一维和二维可压缩无粘问题进行算例验证,并与传统WENO-JS格式和WENO-Z格式进行比较。计算结果表明,改进五阶WENO格式有较高的精度和收敛速度,有较低的数值耗散,能有效捕捉间断、激波和涡等复杂流动。     

求解无粘可压Euler方程组的虚拟流方法

首先将三阶Godunov型半离散中心迎风格式推广到四阶,之后再将该新的四阶半离散中心迎风格式与Level Set方法以及虚拟流方法结合起来,成功地处理了非反应激波问题和多介质流中的爆轰间断问题。由于Level Set函数能隐式地追踪到界面的位置,而虚拟流的构造能隐式地捕捉到界面的边界条件,故而本文的方法可以很自然地推广到多维情况。     

用隐式方法和WENO格式计算悬停旋翼跨声速无粘流场

寻找一种能够准确计算以涡为主要特征的复杂流场和克服尾迹耗散问题的数值方法,一直是旋翼空气动力学研究的热点和难点。本文发展了一种基于高阶迎风格式计算悬停旋翼无粘流场的隐式数值方法。无粘通量采用Roe通量差分分裂格式,为提高精度,使用五阶WENO格式进行左右状态插值,并与MUSCL插值进行比较。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。用该方法对一跨声速悬停旋翼无粘流场进行了数值计算,数值结果表明WENO-Roe的激波分辨率高于MUSCL-Roe,体现出了格式精度的提高对计算结果的改善,LU-SGS隐式方法的计算效率比5步Runge-Kutta显式方法的高。     

Some basic problems in numerically simulating effective properties and local fields of composite materials

It is pointed out that to numerically estimate the effective properties and local fields of matrix-inclusion composites, a commonly adopted method is accompanied with some serious draw-backs. We call this method the nominal loading scheme (NLS), which considers the actual inclusion distribution inside a finite domain, Ω say, treats the external domain of Ω to be of the pure matrix material, and imposes the actural traction, σ^∞ say on the remote boundary. It thus gives rise to the following basic problems: (i) Can NLS be improved remarkably just by adjusting σ^∞? (ii) What is the relationship between the size of Ω and the scale of inclusions? (iii) Which choice is better in calculating the effective properties, the whole domain Ω or an appropriately selected sub-domain of Ω? Targeting these problems, the equivalent loading scheme (ELS) and equivalent matrix scheme (EMS) are proposed. It is theoretically analyzed that both ELS and EMS can be used to precisely simulating the effective properties and local fields of matrix-inclusion composities, and both ELS and EMS are self-approved. As an application, ELS combined with a so-called pseudo-dislocations method is used to evaluate the effective properties and local fields of two-dimensional two-phase composites with close-packed circular inclusions, or randomly distributed circular inclusions, or randomly distributed microcracks. The results show that substituting the remore traction σ^∞ with the effective stress field σ^ɛ suggested by IDD scheme is a simple and effective method, and the estimation of the effective properties and local fields is very close to the accurate solution. The present work is supported by the National Natural Science Foundation of China (No. 19525207).     

Developing shock-capturing difference methods

A new shock-capturing method is proposed which is based on upwind schemes and flux-vector splittings. Firstly, original upwind schemes are projected along characteristic directions. Secondly, the amplitudes of the characteristic decompositions are carefully controlled by limiters to prevent non-physical oscillations. Lastly, the schemes are converted into conservative forms, and the oscillation-free shock-capturing schemes are acquired. Two explicit upwind schemes (2nd-order and 3rd-order) and three compact upwind schemes (3rd-order, 5th-order and 7th-order) are modified by the method for hyperbolic systems and the modified schemes are checked on several one-dimensional and two-dimensional test cases. Some numerical solutions of the schemes are compared with those of a WENO scheme and a MP scheme as well as a compact-WENO scheme. The results show that the method with high order accuracy and high resolutions can capture shock waves smoothly.     

Comments on Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD

The explicit compact difference scheme,proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al.,published in Applied Mathematics and Mechanics (English Edition),2007,28(7),943-953,has the same performance as the conventional finite difference schemes.It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless,we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations,especially for higher accurate schemes.     

Comments on <Emphasis Type="Italic">Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD</Emphasis>

The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics （English Edition）, 2007, 28（7）, 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.     

Uniform algorithm for all-speed shock-capturing schemes

Many ideas exist for the development of shock-capturing schemes, such as Roe, Harten–Lax–van Leer (HLL) and advection upstream splitting method (AUSM) families, and their extension for all-speed flow. A uniform algorithm that expresses the three families in the same framework is proposed in this study. The algorithm has an explicit physical meaning, provides new understanding and comparison of the mechanism of schemes, and may play a significant role in further research. As an example of applying the uniform algorithm, the low Mach number behaviour of the schemes is analysed. A clear and simple explanation is provided based on the wall boundary, and a concise rule is proposed to determine whether a scheme has satisfied low Mach number behaviour.     

一个基于非均匀网格上的二阶基本振荡差分格式

在基本无振荡格式的构造中，将通常的对流通量f的逼近方式推广到对通量导数的逼近，这一构造方法可以有效地应用到非均匀或非结构网格。直接基于非均匀网格上，构造了一个二阶的基本无振荡（ENO）差分格式，该格式具有形式简单，对网格的划分灵活，对传统格式相比不增加计算量等优点，几个数值算例证明了格式的有效性。     

Improved laminar predictions using a stabilised time-dependent simple scheme

A new scheme which can solve unsteady incompressible flows is described in this paper. The scheme is a variant of the SIMPLE methodology. Typically, a scheme of this type tends to suffer from stability problems, which this new scheme overcomes by taking small intermediate steps within a time step. The calculations made in the intermediate steps are damped to enhance the stability of the scheme. The new stabilised scheme is evaluated for laminar flow around a square cylinder, impulsively started laminar flow over a backward-facing step and fluctuating laminar flow over a backward-facing step. Comparisons are made with other numerical predictions and experimental data. In general, good agreement is found, except for the fluctuating laminar flow over a backward-facing step problem. The new scheme is found to have the same level of accuracy, stability and efficiency in comparison with the PISO scheme, but it is easier to code. © 1998 John Wiley & Sons, Ltd.     

Migration-collision scheme of Lattice Boltzmann method for heat conduction problems involving solidification

The phase-change problem is solved by the migration-collision scheme of lattice Boltzmann method. After formula derivation, we find that this method can give a rigorous numerical value for the phase-change temperature, which is of crucial importance. One-dimensional solidification in half-space and two-dimensional solidification in a corner are simulated. The phase change temperature and the liquid-solid interface are both obtained, and the results conform to the analytical solution.     

利用多小波自适应格式求解流体力学方程

高阶计算格式的高精度、高分辨率对提高复杂流场的计算水平有重要的意义,为了提高AUSMPW格式对流场计算中激波等间断的分辨率,减小数值振荡,在原有AUSMPW格式的基础之上,利用多小波对函数进行多尺度分解,并采取阈值的方法生成自适应网格,提出了一种新的基于多小波自适应算法的AUSMPW格式,理论上可以达到任意阶精度. 将所得的压强、密度与原格式、TVD格式及WENO格式的计算结果进行了比较分析. 结果表明改进后的AUSMPW格式较原格式具有更高的分辨率、更强的捕捉间断的能力及更低的数值耗散.