An efficient TVD‐WENO method for conservation laws

In this article, we present a high‐resolution hybrid scheme for solving hyperbolic conservation laws in one and two dimensions. In this scheme, we use a cheap fourth order total variation diminishing (TVD) scheme for smooth region and expensive seventh order weighted nonoscillatory (WENO) scheme near discontinuities. To distinguish between the smooth parts and discontinuities, we use an efficient adaptive multiresolution technique. For time integration, we use the third order TVD Runge‐Kutta scheme. The accuracy of the resulting hybrid high order scheme is comparable with these of WENO, but with significant decrease of the CPU cost. Numerical demonstrates that the proposed scheme is comparable to the high order WENO scheme and superior to the fourth order TVD scheme. Our scheme has the added advantage of simplicity and computational efficiency. Numerical tests are presented which show the robustness and effectiveness of the proposed scheme.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Development and tending to steady state of subsonic gas condensation on a plane condensed phase

The kinetic equation for a monatomic gas with a model collision operator (S-model) is used to study the development and tending to steady state of one-dimensional unsteady half-space gas condensation on a plane condensed phase. Initially, the gas is at rest and in equilibrium with the body’s surface and, then, the body temperature suddenly drops to a constant value. The problem is solved using an implicit second-order accurate quasi-monotone scheme. The process of reaching a steady flow regime is of primary interest. The effect of the evaporation (condensation) coefficient on the flow pattern is analyzed.     

Comparison with experiment for TVD calculations of blast waves from a shock tube

Harten's second-order-accurate total-variation-diminishing (TVD) scheme is applied to calculation of flow from the open end of a shock tube. Comparison of numerical results with available experimental data for overpressure at selected points around the shock tube exit shows good agreement. Numerically indicated positions of the moving shock front and Mach stem also compare well with flow shadowgraph data. Where the problem geometry is sufficiently simple and rectangular gridding can be used, Harten's method affords a good choice for blast wave calculations.     

TVD格式在超音速喷管三维粘性流动求解中的应用

详细给出了任意三维曲线坐标系中Ｎｏｖｉｅｒ－Ｓｔｏｋｅｓ方程的对流项ＴＶＤ格式的构造过程，建立了数值求解三维粘性流动的计算方法，应用该方法对三维超音速喷管中有激波及无激波情况下的两种工况的层流流场进行了数值求解，并与实验做了对比。结果表明本文建立的计算方法具有较高的精度，同时也证明ＴＶＤ格式具有分辩率高，稳定收敛等优点，为进一步开展叶栅流场及紊流的研究打下了基础。     

Numerical Calculation of Flows and Long-Range Transport of Contaminants in Lowland River Reservoirs

Flows and contaminant transport in the Novosibirsk reservoir are calculated on the basis of a two-dimensional (plane) nonstationary model with Saint Venant's equations. The model allows for the presence of a large number of islands. Coefficients of horizontal exchange (dispersion) are calculated by the formula taking into account dynamic velocity at the bottom. Numerical implementation of the model employs a semi-implicit conservative finite-difference TVD scheme on a distributed grid and procedures allowing for the flow past these islands. Model examples of calculations and computation results for dynamics of long-range transport of contaminants along the Novosibirsk reservoir are given.     

A fully implicit time accurate method for hypersonic combustion: application to shock-induced combustion instability

A new fully implicit, time accurate algorithm suitable for chemically reacting, viscous flows in the transonic-to-hypersonic regime is described. The method is based on a class of Total Variation Diminishing (TVD) schemes and uses successive Gauss-Seidel relaxation sweeps. The inversion of large matrices is avoided by partitioning the system into reacting and nonreacting parts, but a fully coupled interaction is still maintained. As a result, the matrices that have to be inverted are of the same size as those obtained with the commonly used point implicit methods. In this paper we illustrate applicability of the new algorithm to hypervelocity unsteady combustion applications. We present a series of numerical simulations of the periodic combustion instabilities observed in ballistic-range experiments of blunt projectiles flying at subdetonative speeds through hydrogenair mixtures. The computed frequencies of oscillation are in excellent agreement with experimental data.     

NND schemes and numerical simulation of axial symmetric free jet flows

Through a study on one-dimensional Navier-Stokes equations, it was found that the spurious oscillations occuring near shock waves with finite difference equations are related to the dispersion term in the corresponding modified differential equations. If the sign of dispersion coefficient is properly adjusted so that the sign changes across shock waves, the undesirable oscillations can be totally suppressed. Based on this finding, the non-oscillatory, containing no free parameters and dissipative shheme (NND scheme) is developed. This scheme is one of “TVD”. The axisymmetric free jet flows are simulated numerically using this scheme. The results obtained by the present scheme are compared with the experimental picture. It is shown that the agreement is very good, and that this scheme has advantages of high resolution for capturing shocks and contact discontinuities. Project supported by National Science Foundation of China     

高超声速三维化学非平衡流动N－S方程数值解

详细地叙述了模拟三维高超声速电离空气化学非平衡粘性流全流场（前体和底部近尾迹）的数值方法，化学反应粘性流求解是基于带有化学源项的Ｎ－Ｓ方程的守恒形式，总的连续方程由组元守恒方程组所代替．对于高温电离空气流动，存在如下７个主要组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ和ｅ￣－．研究发展了求解完全耦合的，与时间相关的偏微分方程组的数值方法．利用一类新的迎风ＴＶＤ激波捕捉有限差分格式求解守恒形式的控制方程组，对高超声速层流有攻角钝锥体绕流流场，在非催化壁条件下、对包括底部近尾迹在内的全流场各参数，如组元浓度，电子密度和温度，以及物面压力分布和热流率分布得出了计算结果．对化学反应气体和完全气体模型的计算结果进行了比较，计算的流场电子密度与飞行实验结果符合很好．     

前向爆轰驱动变截面激波管特性的数值模拟

讨论的激波管由氢氧前向爆轰产生高压驱动气体、并配有收缩截面段,它能产生高焓超高速气流.对此进行了模拟.既对变截面角度的作用,两端截面积比的作用又对比正向反向爆轰有何不同效果做了探讨.正向爆轰时主膜处的收缩段产生的汇聚作用既加强主激波又产生反向激波,缓解了爆轰波阵面后紧跟着的稀疏波导致主激波衰减偏快的不利影响.当收缩角度在30°和45°之间时,主激波的强度较高,衰减较小;当β接近90°时在主激波传播过一小段距离之后,主激波后高温高速气流较均匀,且主激波的衰减最小,具有实用价值.驱动段与被驱动段的截面积比越大,主激波的强度越高,但是最初阶段衰减也越快.反向爆轰时产生的主激波衰减最缓,但是同样的主激波强度需要的驱动段与被驱动段初始压力比前向爆轰高1个量级.     

Efficient solving method for unsteady incompressible interfacial flow problems

Unsteady interfacial problems, considered in an Eulerian form, are studied. The phenomena are modeled using the incompressible viscous Navier–Stokes equations to get the velocity field and an advection equation to predict interface evolutions. The momentum equation is solved by means of an implicit hybrid augmented Lagrangian–Projection method, whereas an explicit characteristic method coupled with a TVD SUPERBEE scheme is applied to the advection equation. The velocity components and the pressure are discretized on staggered grids with finite volumes. Emphasis is on the accuracy and robustness of the techniques described before. A precise explanation on the validation phase will be given, which uses such tests as the advection of a step function or Zalesak's problem to improve the calculation of the interface. The global approach is used on a physically hard interfacial test with strong disparities between viscosities and densities. Copyright © 1999 John Wiley & Sons, Ltd.