High-resolution central difference scheme for the shallow water equations

A two-dimensional nonoscillatory central difference scheme was extended to the shallow water equations. A high-resolution numerical method for solving the shallow water equations was presented. In order to prevent oscillation, the nonlinear limiter is employed to approximate the discrete slopes. The main advantage of the presented method is simplicity comparable with the upwind schemes. This method does not require Riemann solvers or some form of flux difference splitting methods. Furthermore, the discrete derivatives of flux can be approximated by the component-wise approach and thus the computation of Jacobian can be avoided. The method retains high resolution and high accuracy similar to the upwind results. It is applied to simulating several tests, including circular dam-break problem, shock focusing problem and partial dam-break problem. The results are in good agreement with the numerical results obtained by other methods. The simulated results also demonstrate that the presented method is stable and efficient.

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.

Multidomain Hybrid Direct DG and Central Difference Methods for Viscous Terms in Hyperbolic-Parabolic Equations

A class of multidomain hybrid methods of direct discontinuous Galerkin(DDG)methods and central difference(CD)schemes for the viscous terms is pro-posed in this paper.Both conservative and nonconservative coupling modes are dis-cussed.To treat the shock wave,the nonconservative coupling mode automatically switch to conservative coupling mode to preserve the conservative property when discontinuities pass through the artificial interface.To maintain the accuracy of the hybrid methods,the Lagrange interpolation polynomials and their derivatives are reconstructed to handle the coupling cells in the DDG subdomain,while the values of ghost points for the CD subdomain are calculated by the approximate polynomials from the DDG methods.The linear stabilities of these methods are demonstrated in detail through von-Neumann analysis.The multidomain hybrid DDG and CD meth-ods are then extended to one-and two-dimensional hyperbolic-parabolic equations.Numerical results validate that the multidomain hybrid methods are high-order ac-curate in the smooth regions,robust for viscous shock simulations and capable to save computational cost.

UNSTEADY/STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.

A numerical study of 1-D nonlinear P-wave propagation in solid

Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlinear solid media by the numerical method. Distinctly different from the case of numerical modeling of linear elastic wave, there may be several difficulties in the numerical treatment to the nonlinear partial differential equation, such as the steep gradients, shocks and unphysical oscillations. All of them are the great obstacles to the stability and conver-gence of numerical calculation. Fortunately, the comparative study on the modeling of nonlinear wave by the two FD schemes presented in the paper can provide us with an easy method to keep the stability and convergence in the calculation field when the product of the absolute value of nonlinear coefficient and the value of u/x are small enough, namely, the value of bu/x is much smaller than 1. Several results are founded in the numerical study of nonlinear P-wave propagation, such as the waveform aberration, the generation and growth of harmonic wave and the energy redistribution among different frequency components. All of them will be more violent when the initial amplitude A0 is larger or the nonlinearity of medium is stronger. Correspondingly, we have found that the nonlinear P-wave propagation velocity will change with different initial frequency f of source wave or the wave velocity c (equal to the P-wave velocity in the same medium without considering nonlinearity).