Adaptive hierarchic transformations for dynamically p-enriched slope-limiting over discontinuous Galerkin systems of generalized equations

We study a family of generalized slope limiters in two dimensions for Runge–Kutta discontinuous Galerkin (RKDG) solutions of advection-diffusion systems. We analyze the numerical behavior of these limiters applied to a pair of model problems, comparing the error of the approximate solutions, and discuss each limiter’s advantages and disadvantages. We then introduce a series of coupled p-enrichment schemes that may be used as standalone dynamic p-enrichment strategies, or may be augmented via any in the family of variable-in-p slope limiters presented.     

2阶精度时间后差隐式格式TVD充分条件

非定常流动问题计算中常用到含三个时间层的二阶精度时间后差隐式格式,并且希望构成TVD格式,然而理论上的问题多年一直没有解决。本文找到了解决办法,构造了这种类型的隐式TVD格式,证明了其为TVD的充分条件。理论结果为计算所验证,并表明通常未采取本文对时间差分处理方法的格式尚不具备TVD性质。     

高精度TVD有限差分格式的自由梯度条件

1自由流条件守恒性气动方程组的高精度TVD有限差分格式一般通过使用通量限制器加入非线性的人工粘性来构成。在曲线坐标系中，通量限制器比较的两个因子与各自所在单元的界面雅可比转换行列式相关。为了消除畸变网格单元界面雅可比转换行列式的变化对通量限制器比较结果的影响，本文提出了高精度TVD格式的自由梯度性质，并分别就Osher和HartenTVD格式讨论了同时满足自由流和自由梯度的度量张量计算方法，使高精度TVD格式在畸变网格上得以实现。任意曲线坐标系中守恒型二维欧拉方程表示为其中，E、F为笛卡尔坐标系中的通量；J为坐标变…     

Résolution des équations de dispersion–diffusion en présence de forts gradientsResolution of advection–diffusion equation in presence of sharp gradient

In order to resolve the convection terms in the simulation of coastal flows, we propose a numerical method based on the Total Variation Diminishing scheme completed with the Artificial Compression Method to reduce the artificial diffusion and to sharpen the tracer profile in strong gradient regions. This method was applied successfully at discontinuous regions, but out of these regions a limitation of its application is necessary. An application is made to estimate the width of the fresh water band in the Eastern English Channel. To cite this article: C. Sart et al., C. R. Mecanique 330 (2002) 159–165.     

Axisymmetric impingement of a hot jet on a flat plate: equilibrium flow analysis of high-temperature air

The flow structure of an underexpanded supersonic jet with high reservoir temperature impinging on a flat plate has been numerically investigated using a Total Variation Diminishing (TVD) scheme. When the temperature of the flow field is high enough to cause chemical reaction, the specific heat ratio,, is no longer equal to 1.4, nor constant. This explains the difference found in the literature between the flow properties of the calorically perfect gas and that of the chemically reacting flow. Under the equilibrium flow assumption the effect of high temperature gas on the impinging jet has been taken into account in the present paper by using specific heat ratio and speed of sound given by correlation polynomials of thermodynamic variables. The limiting case of cold jet calculation in the present numerical results agreed well with the existing experimental data. For the equilibrium jet with high reservoir temperature,T _o=1000K, qualitative support of the present result has been provided by means of the approximation theory.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

运动激波和气泡串相互作用的初步数值模拟

通过对激波和流体界面相互作用诱导的大变形界面演化的数值模拟,验证Level set方法精确模拟多个流体界面的有效性.采用2阶迎风TVD求解欧拉方程得到流场解,采用5阶WENO求解Level set方程追踪多流体界面,采用GFM方法处理流体内界面.利用文[1]的计算结果校核本文程序.在此基础上,对运动激波和气泡串相互作用过程进行了初步数值模拟,得到了不同时刻运动激波和圆管内的两个气泡作用后的演化图象,包括压力和密度等值线分布.计算结果表明:针对推广后的多界面Level set方程,该方法仍可高质量地捕捉多个流体界面.     

一类基于通量分裂的二阶精度TVD差分格式

基于通量分裂和逆风特性,选取单元交界面上的正、负数值通量,并结合高阶Runge-Kutta TVD时间离散方法,构造了一维非线性双曲型守恒律的一类二阶精度的差分格式,证明了差分格式的TVD特性。按分量形式推广到方程组。通过几个典型的数值算例验证了格式的有效性。     

Unsteady supersonic flow over an oscillatory aerofoil using a second-order TVD scheme

The unsteady flow over an oscillatory NACA0012 aerofoil has been simulated by the calculation with Euler equations. The equations are discretized by an implicit Euler in time, and a second-order space-accurate TVD scheme based on flux vector splitting with van Leer's limiter. Modified eigenvalues are proposed to overcome the slope discontinuities of split eigenvalues at Mach = 0·0 and ± 1·0, and to generate a bow shock in front of the aerofoil. A moving grid system around the aerofoil is generated by Sorenson's boundary fitted co-ordinates for each time step. The calculations have been done for two angles of attack θ = 5·0° sin (ωt) and θ = 3·0° + 3·0° sin (ωt) for the free-stream Mach numbers 2·0 and 3·0. The results show that pressure and Mach cells flow along characteristic lines. To examine unsteady effects, the responses of wall pressure and normal force coefficients are analysed by a Fourier series expansion.     

Comparisons of TVD schemes for turbulent transonic projectile aerodynamics computations with a two-equation model of turbulence

The development of a computer program to solve the axisymmetric full Navier--Stokes equations with k-ε two-equation model of turbulence using various total variation diminishing (TVD) schemes is the primary interest of this study. The computations are performed for the turbulent, transonic, viscous flow over a projectile with/without supporting sting at zero angle of attack. The predicted results, as well as the convergence characteristics, by various TVD schemes are compared with each other. The results show that the TVD schemes of higher-order accuracy do have influence on the regions of high gradients such as shock, base corner and base flow. However, the schemes of third-order accuracy do not necessarily improve the agreement with measured data (which is not available on the base) than that of second-order accuracy, but surely generate apparent different result of base flow. The supporting sting on the projectile base will complicate the base flow and the existence of the sting will slightly shift the shock location and slightly change the flow field after the shock. More iteration steps are needed to get the converged results in the computation for the projectile with sting.     

A comparison of spatial discretization schemes for differential solution methods of the radiative transfer equation

A comparison of discretization schemes required to evaluate the radiation intensity at the cell faces of a control volume in differential solution methods of the radiative transfer equation is presented. Several schemes developed using the normalized variable diagram and the total variation diminishing formalisms are compared along with essentially non-oscillatory schemes and genuinely multidimensional schemes. The calculations were carried out using the discrete ordinates method, but the analysis is equally valid for the finite-volume method. It is shown that the S schemes of the genuinely multidimensional family perform quite well, particularly in problems with discontinuous radiation intensity fields. However, they are time consuming, and so they do not always become more attractive regarding the trade-off between accuracy and computational requirements, in comparison with other high-order schemes that, although being less accurate, are also more economical.