Characteristic-based finite volume scheme for 1D Euler equations

In this paper, a high-order finite-volume scheme is presented for the one- dimensional scalar and inviscid Euler conservation laws. The Simpson＇s quadrature rule is used to achieve high-order accuracy in time. To get the point value of the Simpson＇s quadrature, the characteristic theory is used to obtain the positions of the grid points at each sub-time stage along the characteristic curves, and the third-order and fifth-order central weighted essentially non-oscillatory （CWENO） reconstruction is adopted to estimate the cell point values. Several standard one-dimensional examples are used to verify the high-order accuracy, convergence and capability of capturing shock.     

Non-oscillatory shock-capturing finite element methods for the one-dimensional compressible Euler equations

A class of shock-capturing Petrov–Galerkin finite element methods that use high-order non-oscillatory interpolations is presented for the one-dimensional compressible Euler equations. Modified eigenvalues which employ total variation diminishing (TVD), total variation bounded (TVB) and essentially non-oscillatory (ENO) mechanisms are introduced into the weighting functions. A one-pass Euler explicit transient algorithm with lumped mass matrix is used to integrate the equations. Numerical experiments with Burgers' equation, the Riemann problem and the two-blast-wave interaction problem are presented. Results indicate that accurate solutions in smooth regions and sharp and non-oscillatory solutions at discontinuities are obtainable even for strong shocks.     

Non-oscillatory central-upwind scheme for hyperbolic conservation laws

We introduce a new fourth order, semi-discrete, central-upwind scheme for solving systems of hyperbolic conservation laws. The scheme is a combination of a fourth order non-oscillatory reconstruction, a semi-discrete central-upwind numerical flux and the third order TVD Runge-Kutta method. Numerical results suggest that the new scheme achieves a uniformly high order accuracy for smooth solutions and produces non-oscillatory profiles for discontinuities. This is especially so for long time evolution problems. The scheme combines the simplicity of the central schemes and accuracy of the upwind schemes. The advantages of the new scheme will be fully realized when solving various examples.     

Some explicit triangular finite element schemes for the Euler equations

Several explicit schemes are presented for triangular P₀ and P₁ finite elements. A first-order accurate upwind P₀ scheme is compared to a FLIC type method. A second-order accurate Richtmyer scheme is constructed. Applications are given for the Euler system of conservation laws in the 2-dimensional case.     

Characteristic-based finite volume scheme for 1D Euler eouations

In this paper, a high-order finite-volume scheme is presented for the one-dimensional scalar and inviscid Euler conservation laws. The Simpson's quadrature rule is used to achieve high-order accuracy in time. To get the point value of the Simpson's quadrature, the characteristic theory is used to obtain the positions of the grid points at each sub-time stage along the characteristic curves, and the third-order and fifth-order central weighted essentially non-oscillatory (CWENO) reconstruction is adopted to estimate the cell point values. Several standard one-dimensional examples are used to verify the high-order accuracy, convergence and capability of capturing shock.     

Improved artificial dissipation schemes for the Euler equations

The influence of artificial dissipation schemes on the accuracy and stability of the numerical solution of compressible flow is extensively examined. Using an implicit central difference factored scheme, an improved form of artificial dissipation is introduced which highly reduces the errors due to numerical viscosity. A function of the local Mach number is used to scale the amount of numerical damping added into the solution according to the character of the flow in several flow regimes. The resulting scheme is validated through several inviscid flow test cases.     

A shock-capturing scheme for body-fitted meshes

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional Euler equations of gas dynamics in a generalized co-ordinate system. The scheme is based on numerical characteristic decomposition and solves locally linearized Riemann problems using upwind differencing. The decomposition is for a generalized co-ordinate system and a convex equation of state. This ensures good shock-capturing properties when incorporated with operator splitting and the advantage of using body-fitted co-ordinates. The resulting scheme is applied to supersonic flow of real air' past a circular cylinder.     

Symmetric Gauss-Seidel multigrid solution of the Euler equations on structured and unstructured grids

The efficient symmetric Gauss-Seidel (SGS) algorithm for solving the Euler equations of inviscid, compressible flow on structured grids, developed in collaboration with Jameson of Stanford University, is extended to unstructured grids. The algorithm uses a nonlinear formulation of an SGS solver, implemented within the framework of multigrid. The earlier form of the algorithm used the natural (lexicographic) ordering of the mesh cells available on structured grids for the SGS sweeps, but a number of features of the method that are believed to contribute to its success can also be implemented for computations on unstructured grids. The present paper reviews, the features of the SGS multigrid solver for structured gr0ids, including its nonlinear implementation, its use of “absolute” Jacobian matrix preconditioning, and its incorporation of multigrid, and then describes the incorporation of these features into an algorithm suitable for computations on unstructured grids. The implementation on unstructured grids is based on the agglomerated multigrid method developed by Sørensen, which uses an explicit Runge-Kutta smoothing algorithm. Results of computations for steady, transonic flows past two-dimensional airfoils are presented, and the efficiency of the method is evaluated for computations on both structured and unstructured meshes.     

Second‐order entropy diminishing scheme for the Euler equations

This paper presents a method of controlling the water levels in a conduit system by employing optimal control theory and the finite element method. A shallow‐water equation is employed for the analysis of flow behaviour. Optimal control theory is utilized to obtain a control value for the target state value. The Sakawa–Shindo method is employed as a minimization technique. For the computational storage requirements, the time domain decomposition method is applied. The Crank–Nicolson method is used for temporal discretization. In addition to a method for optimally controlling water level, a method is presented for determining transversality conditions, the terminal condition of the Lagrange multiplier. Copyright © 2006 John Wiley & Sons, Ltd.     

Approximate symmetries and conservation laws of the geodesic equations for the Schwarzschild metric

Approximate symmetries have been defined in the context of differential equations and systems of differential equations. They give approximately, conserved quantities for Lagrangian systems. In this paper, the exact and the approximate symmetries of the system of geodesic equations for the Schwarzschild metric, and in particular for the radial equation of motion, are studied. It is noted that there is an ambiguity in the formulation of approximate symmetries that needs to be clarified by consideration of the Lagrangian for the system of equations. The significance of approximate symmetries in this context is discussed.     

An adaptive multigrid finite-volume scheme for incompressible Navier–Stokes equations

An algorithm for the solutions of the two-dimensional incompressible Navier–Stokes equations is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. It is based on an artificial compressibility method and uses higher-order upwind finite-volume techniques for the convective terms and a second-order finite-volume technique for the viscous terms. Three upwind schemes for discretizing convective terms are proposed here. An interesting result is that the solutions computed by one of them is not sensitive to the value of the artificial compressibility parameter. A second-order, two-step Runge–Kutta integration coupling with an implicit residual smoothing and with a multigrid method is used for achieving fast convergence for both steady- and unsteady-state problems. The numerical results agree well with experimental and other numerical data. A comparison with an analytically exact solution is performed to verify the space and time accuracy of the algorithm.     

求解一维理想磁流体方程的移动网格熵稳定格式

磁流体方程的数值求解在等离子体物理学、天体物理研究以及流动控制等领域具有重要意义,本文构造了用于求解理想磁流体动力学方程的基于移动网格的熵稳定格式,此方法将Roe型熵稳定格式与自适应移动网格算法结合,空间方向采用熵稳定格式对磁流体动力学方程进行离散,利用变分法构造网格演化方程并通过Gauss-Seidel迭代法对其迭代求解实现网格的自适应分布,在此基础上采用守恒型插值公式实现新旧节点上的量值传递,利用三阶强稳定Runge-Kutta方法将数值解推进到下一时间层。数值实验表明,该算法能有效捕捉解的结构(特别是激波和稀疏波),分辨率高,通用性好,具有强鲁棒性。     

Numerical solutions of incompressible Euler and Navier-Stokes equations by efficient discrete singular convolution method

An efficient discrete singular convolution (DSC) method is introduced to the numerical solutions of incompressible Euler and Navier-Stokes equations with periodic boundary conditions. Two numerical tests of two-dimensional Navier-Stokes equations with periodic boundary conditions and Euler equations for doubly periodic shear layer flows are carried out by using the DSC method for spatial derivatives and fourth-order Runge-Kutta method for time advancement, respectively. The computational results show that the DSC method is efficient and robust for solving the problems of incompressible flows, and has the potential of being extended to numerically solve much broader problems in fluid dynamics. The project supported by the National Natural Science Foundation of China (No.19902010).     

一般形式的二维时—空守恒格式

本文将经作者改进后的一维时－空守恒格式推广到了二维情形，得到了一个一般形式的二维Ｅｕｌｅｒ方程时－空守恒格式，该格式对各种不规则几何区域内的流动问题具有很强的适应性，同时它还保留了一维格式的优点。几个典型算例的计算结果表明，本文格式不仅精度高，通用性好，而且对激波等间断具有很高的分辨率。     

Computation of three-dimensional flow using the Euler equations and a multiple-grid scheme

The unsteady Euler equations are numerically solved using the finite volume one-step scheme recently developed by Ron-Ho Ni. The multiple-grid procedure of Ni is also implemented. The flows are assumed to be homo-enthalpic; the energy equation is eliminated and the static pressure is determined by the steady Bernoulli equation; a local time-step technique is used. Inflow and outflow boundaries are treated with the compatibility relations method of ONERA. The efficiency of the multiple-grid scheme is demonstrated by a two-dimensional calculation (transonic flow past the NACA 12 aerofoil) and also by a three-dimensional one (transonic lifting flow past the M6 wing). The third application presented shows the ability of the method to compute the vortical flow around a delta wing with leading-edge separation. No condition is applied at the leading-edge; the vortex sheets are captured in the same sense as shock waves. Results indicate that the Euler equations method is well suited for the prediction of flows with shock waves and contact discontinuities, the multiple-grid procedure allowing a substantial reduction of the computational time.     

A multigrid solver for the vorticity-velocity Navier-Stokes equations

This paper provides a multigrid incremental line-Gauss-Seidel method for solving the steady Navier-Stokes equations in two and three dimensions expressed in terms of the vorticity and velocity variables. The system of parabolic and Poisson equations governing the scalar components of the vector unknowns is solved using centred finite differences on a non-staggered grid. Numerical results for the two-dimensional driven cavity problem indicate that the spatial discretization of the equation defining the value of the vorticity on the boundary is extremely critical to obtaining accurate solutions. In fact, a standard one-sided three-point second-order-accurate approximation produces very inaccurate results for moderate-to-high values of the Reynolds number unless an exceedingly fine mesh is employed. On the other hand, a compact two-point second-order-accurate discretization is found to be always satisfactory and provides accurate solutions for Reynolds number up to 3200, a target impossible heretofore using this formulation and a non-staggered grid.     

A new flux-limiting approach–based kinetic scheme for the Euler equations of gas dynamics

This paper proposes a new kinetic-theory-based high-resolution scheme for the Euler equations of gas dynamics. The scheme uses the well-known connection that the Euler equations are suitable moments of the collisionless Boltzmann equation of kinetic theory. The collisionless Boltzmann equation is discretized using Sweby's flux-limited method and the moment of this Boltzmann level formulation gives a Euler level scheme. It is demonstrated how conventional limiters and an extremum-preserving limiter can be adapted for use in the scheme to achieve a desired effect. A simple total variation diminishing criteria relaxing parameter results in improving the resolution of the discontinuities in a significant way. A 1D scheme is formulated first and an extension to 2D on Cartesian meshes is carried out next. Accuracy analysis suggests that the scheme achieves between first- and second-order accuracy as is expected for any second-order flux-limited method. The simplicity and the explicit form of the conservative numerical fluxes add to the efficiency of the scheme. Several standard 1D and 2D test problems are solved to demonstrate the robustness and accuracy.     

Use of the splitting scheme and multigrid method to compute flow separation

The splitting difference scheme is used to study flow separation. Flows behind a circular cylinder are computed as a model problem. In view of the nature of the flow, the variables are transformed. The boundary condition for the pressure is given from an intermediate velocity. The free-slip velocity boundary conditions on the rigid wall are given by interpolation. The multigrid algorithm is applied to the pressure iteration. We also choose better initial values for the model problem by means of the multigrid algorithm idea.     

An assessment of the finite termination property of the defect correction method

This paper investigates the use of defect correction procedures for the solution of finite volume approximations to systems of conservation laws. Particular emphasis is laid on the order of accuracy obtained after a fixed finite number of iterations. It is shown that a high order of accuracy may be achieved after only one defect correction iteration, involving two inversions of a stable lower-order-accurate operator. However, this result is found to be critically dependent on the consistency of the lower-order operator, a property which does not always hold for conservative finite volume discretizations. Through numerical experiments, the lack of consistency of these schemes is found to inhibit severely the finite termination property of the defect correction process. Results are presented for linear advection, Poisson's equation, and the Euler equations.     

A second-order positivity-preserving well-balanced finite volume scheme for Euler equations with gravity for arbitrary hydrostatic equilibria

We present a well-balanced finite volume scheme for the compressible Euler equations with gravity, where the approximate Riemann solver is derived using a Suliciu relaxation approach. Besides the well-balanced property, the scheme is robust with respect to the physical admissible states. General hydrostatic solutions are captured up to machine precision by deriving, for a given initial value problem, suitable time-independent functions and using them in the discretization of the source term. The first-order scheme is extended to a second-order scheme by reconstructing in equilibrium variables while preserving the well-balanced and robustness properties. Numerical examples are performed to demonstrate the accuracy, well-balanced, and robustness properties of the presented scheme for up to three space dimensions.