Improvement of convergence to steady state solutions of Euler equations with weighted compact nonlinear schemes

The convergence to steady state solutions of the Euler equations for weighted compact nonlinear schemes (WCNS) [Deng X. and Zhang H. (2000), J. Comput. Phys. 165, 22-44 and Zhang S., Jiang S. and Shu C.-W. (2008), J. Comput. Phys. 227, 7294-7321] is studied through numerical tests. Like most other shock capturing schemes, WCNS also suffers from the problem that the residue can not settle down to machine zero for the computation of the steady state solution which contains shock waves but hangs at the truncation error level. In this paper, the techniques studied in [Zhang S. and Shu. C.-W. (2007), J. Sci. Comput. 31, 273-305 and Zhang S., Jiang S and Shu. C.-W. (2011), J. Sci. Comput. 47, 216-238], to improve the convergence to steady state solutions for WENO schemes, are generalized to the WCNS. Detailed numerical studies in one and two dimensional cases are performed. Numerical tests demonstrate the effectiveness of these techniques when applied to WCNS. The residue of various order WCNS can settle down to machine zero for typical cases while the small post-shock oscillations can be removed.     

High-order accurate difference schemes for solving gasdynamic equations by the Godunov method with antidiffusion

A technique is proposed for improving the accuracy of the Godunov method as applied to gasdynamic simulations in one dimension. The underlying idea is the reconstruction of fluxes arsoss cell boundaries (“large” values) by using antidiffusion corrections, which are obtained by analyzing the differential approximation of the schemes. In contrast to other approaches, the reconstructed values are not the initial data but rather large values calculated by solving the Riemann problem. The approach is efficient and yields higher accuracy difference schemes with a high resolution.     

Unstructured-grid finite-volume discretization of the Navier-Stokes equations based on high-resolution difference schemes

An unstructured-grid discretization of the Navier-Stokes equations based on the finite volume method and high-resolution difference schemes in time and space is described as applied to fluid dynamics problems in two and three dimensions. The control volume is defined as the cell-vertex median dual control volume. The fluxes through the faces of internal and boundary control volumes are written identically, which simplifies their software implementation. The gradient and the pseudo-Laplacian are calculated at the midpoint of a control volume face by using relations adapted to the computations on a strongly stretched grid in the boundary layer.     

High-order accurate monotone compact running scheme for multidimensional hyperbolic equations

Monotone absolutely stable conservative difference schemes intended for solving quasilinear multidimensional hyperbolic equations are described. For sufficiently smooth solutions, the schemes are fourth-order accurate in each spatial direction and can be used in a wide range of local Courant numbers. The order of accuracy in time varies from the third for the smooth parts of the solution to the first near discontinuities. This is achieved by choosing special weighting coefficients that depend locally on the solution. The presented schemes are numerically efficient thanks to the simple two-diagonal (or block two-diagonal) structure of the matrix to be inverted. First the schemes are applied to system of nonlinear multidimensional conservation laws. The choice of optimal weighting coefficients for the schemes of variable order of accuracy in time and flux splitting is discussed in detail. The capabilities of the schemes are demonstrated by computing well-known two-dimensional Riemann problems for gasdynamic equations with a complex shock wave structure.     

Optimal first- to sixth-order accurate Runge-Kutta schemes

An optimal choice of free parameters in explicit Runge-Kutta schemes up to the sixth order is discussed. A sixth-order seven-stage scheme that is immediately ahead of Butcher’s second barrier is constructed. The study is performed in the most general form, and its results are applicable to both autonomous and nonautonomous problems.     

Conservative compact difference schemes for the coupled nonlinear schrödinger system

In this article, some conservative compact difference schemes are explored for the strongly coupled nonlinear schrödinger system. After transforming the scheme into matrix form, we prove the existence and uniqueness, convergence and stability of the difference solutions for one nonlinear scheme in the norm by using some techniques of matrix theory. Numerical results show that one nonlinear scheme is the most efficient of all the compact schemes constructed here. It allows much larger time steps than the others. The second most efficient compact scheme is a linear one. We then give numerical simulations to two soliton interactions for the two most efficient compact schemes. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 749–772, 2014     

New higher-order compact finite difference schemes for 1D heat conduction equations

In this paper, we present two higher-order compact finite difference schemes for solving one-dimensional (1D) heat conduction equations with Dirichlet and Neumann boundary conditions, respectively. In particular, we delicately adjust the location of the interior grid point that is next to the boundary so that the Dirichlet or Neumann boundary condition can be applied directly without discretization, and at the same time, the fifth or sixth-order compact finite difference approximations at the grid point can be obtained. On the other hand, an eighth-order compact finite difference approximation is employed for the spatial derivative at other interior grid points. Combined with the Crank–Nicholson finite difference method and Richardson extrapolation, the overall scheme can be unconditionally stable and provides much more accurate numerical solutions. Numerical errors and convergence rates of these two schemes are tested by two examples.     

Application of compact and multioperator schemes to the numerical simulation of acoustic fields generated by instability waves in supersonic jets

Acoustic fields generated by instability waves in supersonic jets were numerically simulated. A seventh-order multioperator scheme was used to solve the Euler equations linearized about the mean flow field in an axisymmetric turbulent jet. The mean field was computed using fifth-order compact approximations of the convective terms under conditions similar to experimental data. The numerical results were found to agree well with the experiment.     

一维非线性Schrödinger 方程的两个无条件收敛的守恒紧致差分格式

本文对一维非线性 Schrödinger 方程给出两个紧致差分格式, 运用能量方法和两个新的分析技 巧证明格式关于离散质量和离散能量守恒, 而且在最大模意义下无条件收敛. 对非线性紧格式构造了 一个新的迭代算法, 证明了算法的收敛性, 并在此基础上给出一个新的线性化紧格式. 数值算例验证 了理论分析的正确性, 并通过外推进一步提高了数值解的精度.     

On convergence and performance of iterative methods with fourth-order compact schemes

We study the convergence and performance of iterative methods with the fourth-order compact discretization schemes for the one- and two-dimensional convection–diffusion equations. For the one-dimensional problem, we investigate the symmetrizability of the coefficient matrix and derive an analytical formula for the spectral radius of the point Jacobi iteration matrix. For the two-dimensional problem, we conduct Fourier analysis to determine the error reduction factors of several basic iterative methods and comment on their potential use as the smoothers for the multilevel methods. Finally, we perform numerical experiments to verify our Fourier analysis results. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14:263–280, 1998     

Remarks on the nonlinear instability of incompressible Euler equations

In this paper, we consider the nonlinear instability of incompressible Euler equations. If a steady density is non-monotonic, then the smooth steady state is a nonlinear instability. First, we use variational method to find a dominant eigenvalue which is important in the construction of approximate solutions, then by energy technique and analytic method, we obtain the dynamical instability result.     

Fourth‐order compact schemes of a heat conduction problem with Neumann boundary conditions

In this article, a set of fourth‐order compact finite difference schemes is developed to solve a heat conduction problem with Neumann boundary conditions. It is derived through the compact difference schemes at all interior points, and the combined compact difference schemes at the boundary points. This set of schemes is proved to be globally solvable and unconditionally stable. Numerical examples are provided to verify the accuracy.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Parameterized families of finite difference schemes for the wave equation

This article is devoted to an analysis of simple families of finite difference schemes for the wave equation. These families are dependent on several free parameters, and methods for obtaining stability bounds as a function of these parameters are discussed in detail. Access to explicit stability bounds such as those derived here may, it is hoped, lead to optimization techniques for so‐called spectral‐like methods, which are difference schemes dependent on many free parameters (and for which maximizing the order of accuracy may not be the defining criterion). Though the focus is on schemes for the wave equation in one dimension, the analysis techniques are extended to two dimensions; implicit schemes such as ADI methods are examined in detail. Numerical results are presented. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 463–480, 2004.     

Compact almost automorphic solutions for some nonlinear dissipative differential equations in Banach spaces

The aim of this work is to prove the existence and uniqueness of compact almost automorphic solutions for some dissipative differential equations in Banach spaces when the input function is only almost automorphic in the sense of Stepanov. Examples and a numerical simulation are provided to illustrate the theoretical findings.     

Note on a proposed weighted residual method for solving nonlinear differential boundary problems

A proposed iterative method for solving nonlinear differential two-point boundary value problems is generalized and shown to be equivalent to using the standard method of adjoints on a linearized form of the nonlinear boundary-value problem.     

Lattice-Boltzmann type relaxation systems and high order relaxation schemes for the incompressible Navier-Stokes equations

A relaxation system based on a Lattice-Boltzmann type discrete velocity model is considered in the low Mach number limit. A third order relaxation scheme is developed working uniformly for all ranges of the mean free path and Mach number. In the incompressible Navier-Stokes limit the scheme reduces to an explicit high order finite difference scheme for the incompressible Navier-Stokes equations based on nonoscillatory upwind discretization. Numerical results and comparisons with other approaches are presented for several test cases in one and two space dimensions.

     

Method of weighted residuals as applied to nonlinear differential equations

A new method for solving a class of nonlinear boundary-value problems is presented. In this method, the nonlinear equation is linearized by guessing an initial solution and using it to evaluate the nonlinear terms. Next, a method of weighted residuals is applied to transform the linearized form of the boundary value problem to an initial value problem. The second (improved) solution is obtained by integrating the initial value problem by a fourth order Runge-Kutta scheme. The entire process is repeated until a desired convergence criterion is achieved.