A dynamically optimized finite difference scheme for Large-Eddy Simulation

A low-dispersive dynamic finite difference scheme for Large-Eddy Simulation is developed. The dynamic scheme is constructed by combining Taylor series expansions on two different grid resolutions. The scheme is optimized dynamically through the real-time adaption of a dynamic coefficient according to the spectral content of the flow, such that the global dispersion error is minimal. In the case of DNS-resolution, the dynamic scheme reduces to the standard Taylor-based finite difference scheme with formal asymptotic order of accuracy. When going to LES-resolution, the dynamic scheme seamlessly adapts to a dispersion-relation preserving scheme. The scheme is tested for Large-Eddy Simulation of Burgers equation. Very good results are obtained.     

Mathematical proof of closed form expressions for finite difference approximations based on Taylor series

Taylor series based finite difference approximations of derivatives of a function have already been presented in closed forms, with explicit formulas for their coefficients. However, those formulas were not derived mathematically and were based on observation of numerical results. In this paper, we provide a mathematical proof of those formulas by deriving them mathematically from the Taylor series.     

High order Hybrid central—WENO finite difference scheme for conservation laws

In this article we present a high resolution hybrid central finite difference—WENO scheme for the solution of conservation laws, in particular, those related to shock–turbulence interaction problems. A sixth order central finite difference scheme is conjugated with a fifth order weighted essentially non-oscillatory WENO scheme in a grid-based adaptive way. High order multi-resolution analysis is used to detect the high gradients regions of the numerical solution in order to capture the shocks with the WENO scheme while the smooth regions are computed with the more efficient and accurate central finite difference scheme. The application of high order filtering to mitigate the dispersion error of central finite difference schemes is also discussed. Numerical experiments with the 1D compressible Euler equations are shown.     

An Engquist-Osher type finite difference scheme with a discontinuous flux function in space

An Engquist-Osher type finite difference scheme is derived for dealing with scalar conservation laws having a flux that is spatially dependent through a possibly discontinuous coefficient. The new monotone difference scheme is based on introducing a new interface numerical flux function, which is called a generalized Engquist-Osher flux. By means of this scheme, the existence and uniqueness of weak solutions to the scalar conservation laws are obtained and the convergence theorem is established. Some numerical examples are presented and the corresponding numerical results are displayed to illustrate the efficiency of the methods.     

Uniformly convergent novel finite difference methods for singularly perturbed reaction–diffusion equations

In this paper, we construct a kind of novel finite difference (NFD) method for solving singularly perturbed reaction–diffusion problems. Different from directly truncating the high‐order derivative terms of the Taylor's series in the traditional finite difference method, we rearrange the Taylor's expansion in a more elaborate way based on the original equation to develop the NFD scheme for 1D problems. It is proved that this approach not only can highly improve the calculation accuracy but also is uniformly convergent. Then, applying alternating direction implicit technique, the newly deduced schemes are extended to 2D equations, and the uniform error estimation based on Shishkin mesh is derived, too. Finally, numerical experiments are presented to verify the high computational accuracy and theoretical prediction.     

High‐order finite difference schemes for numerical solutions of the generalized Burgers–Huxley equation

In this article, up to tenth‐order finite difference schemes are proposed to solve the generalized Burgers–Huxley equation. The schemes based on high‐order differences are presented using Taylor series expansion. To establish the numerical solutions of the corresponding equation, the high‐order schemes in space and a fourth‐order Runge‐Kutta scheme in time have been combined. Numerical experiments have been conducted to demonstrate the high‐order accuracy of the current algorithms with relatively minimal computational effort. The results showed that use of the present approaches in the simulation is very applicable for the solution of the generalized Burgers–Huxley equation. The current results are also seen to be more accurate than some results given in the literature. The proposed algorithms are seen to be very good alternatives to existing approaches for such physical applications. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1313‐1326, 2011     

An optimal 25-point finite difference scheme for the Helmholtz equation with PML

In this paper, we present an optimal 25-point finite difference scheme for solving the Helmholtz equation with perfectly matched layer (PML) in two dimensional domain. Based on minimizing the numerical dispersion, we propose the refined choice strategy for choosing optimal parameters of the 25-point finite difference scheme. Numerical experiments are given to illustrate the improvement of the accuracy and the reduction of the numerical dispersion.     

An optimal 25-point finite difference scheme for the Helmholtz equation with PML

In this paper, we present an optimal 25-point finite difference scheme for solving the Helmholtz equation with perfectly matched layer (PML) in two dimensional domain. Based on minimizing the numerical dispersion, we propose the refined choice strategy for choosing optimal parameters of the 25-point finite difference scheme. Numerical experiments are given to illustrate the improvement of the accuracy and the reduction of the numerical dispersion.     

A generalized finite difference method for solid mechanics

Procedures are developed that improve the applicability of the finite difference method to problems in solid mechanics. This is accomplished by formulating the coefficients of the Taylor series expansion used to approximate derivative quantities in terms of physically interpretable strain gradients. Improvements realized include modeling of boundary conditions that has intuitive appeal and the use of irregular grids in a natural manner. These developments are demonstrated for the analysis of plane stress problems with traction boundary conditions. The results compare well with finite element solutions. The approach suggests further generalization of the finite difference method.     

Optimized composite finite difference schemes for atmospheric flow modeling

In this paper, we use some finite difference methods in order to solve an atmospheric flow problem described by an advection–diffusion equation. This flow problem was solved by Clancy using forward‐time central space (FTCS) scheme and is challenging to simulate due to large errors in phase and amplitude which are generated especially over long propagation times. Clancy also derived stability limits for FTCS scheme. We use Von Neumann stability analysis and the approach of Hindmarsch et al. which is an improved technique over that of Clancy in order to obtain the region of stability of some methods such as FTCS, Lax–Wendroff (LW), Crank–Nicolson. We also construct a nonstandard finite difference (NSFD) scheme. Properties like stability and consistency are studied. To improve the results due to significant numerical dispersion or numerical dissipation, we derive a new composite scheme consisting of three applications of LW followed by one application of NSFD. The latter acts like a filter to remove the dispersive oscillations from LW. We further improve the composite scheme by computing the optimal temporal step size at a given spatial step size using two techniques namely; by minimizing the square of dispersion error and by minimizing the sum of squares of dispersion and dissipation errors.     

Discontinuous Galerkin and mimetic finite difference methods for coupled Stokes–Darcy flows on polygonal and polyhedral grids

We study locally mass conservative approximations of coupled Darcy and Stokes flows on polygonal and polyhedral meshes. The discontinuous Galerkin (DG) finite element method is used in the Stokes region and the mimetic finite difference method is used in the Darcy region. DG finite element spaces are defined on polygonal and polyhedral grids by introducing lifting operators mapping mimetic degrees of freedom to functional spaces. Optimal convergence estimates for the numerical scheme are derived. Results from computational experiments supporting the theory are presented.     

An accurate semi-analytic finite difference scheme for two-dimensional elliptic problems with singularities

A high-order semi-analytic finite difference scheme is presented to overcome degradation of numerical performance when applied to two-dimensional elliptic problems containing singular points. The scheme, called Least-Square Singular Finite Difference Scheme (L-S SFDS), applies an explicit functional representation of the exact solution in the vicinity of the singularities, and a conventional finite difference scheme on the remaining domain. It is shown that the L-S SFDS is “pollution” free, i.e., no degradation in the convergence rate occurs because of the singularities, and the coefficients of the asymptotic solution in the vicinity of the singularities are computed as a by-product with a very high accuracy. Numerical examples for the Laplace and Poisson equations over domains containing re-entrant corners or abrupt changes in the boundary conditions are presented. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 281–296, 1998     

A finite difference scheme for the nonlinear time-fractional partial integro-differential equation

In this paper, a finite difference scheme is proposed for solving the nonlinear time-fractional integro-differential equation. This model involves two nonlocal terms in time, ie, a Caputo time-fractional derivative and an integral term with memory. The existence of numerical solutions is shown by the Leray-Schauder theorem. And we obtain the discrete L₂ stability and convergence with second order in time and space by the discrete energy method. Then the uniqueness of numerical solutions is derived. Moreover, an iterative algorithm is designed for solving the derived nonlinear system. Numerical examples are presented to validate the theoretical findings and the efficiency of the proposed algorithm.     

A finite difference method for fractional partial differential equation

An implicit unconditional stable difference scheme is presented for a kind of linear space–time fractional convection–diffusion equation. The equation is obtained from the classical integer order convection–diffusion equations with fractional order derivatives for both space and time. First-order consistency, unconditional stability, and first-order convergence of the method are proven using a novel shifted version of the classical Grünwald finite difference approximation for the fractional derivatives. A numerical example with known exact solution is also presented, and the behavior of the error is examined to verify the order of convergence.     

The characteristic finite difference streamline diffusion method for convection-dominated diffusion problems

In this paper, we consider the characteristic finite difference streamline diffusion method for two-dimensional convection-dominated diffusion problems. The scheme is combined the method of characteristics with the finite difference streamline diffusion (FDSD) method to create the characteristic FDSD (C-FDSD) procedures. Stability analysis and error estimate of the C-FDSD method are deduced. The scheme not only realizes the purpose of lowering the time-truncation error, using larger time step for solving the convection-dominated diffusion problems, but also keeps the favorable stability and high precision of the FDSD method. Finally, numerical experiments are presented to illustrate the availability of the scheme.     

A monotone scheme for Hamilton–Jacobi equations via the nonstandard finite difference method

A usual way of approximating Hamilton–Jacobi equations is to couple space finite element discretization with time finite difference discretization. This classical approach leads to a severe restriction on the time step size for the scheme to be monotone. In this paper, we couple the finite element method with the nonstandard finite difference method, which is based on Mickens' rule of nonlocal approximation. The scheme obtained in this way is unconditionally monotone. The convergence of the new method is discussed and numerical results that support the theory are provided. Copyright © 2009 John Wiley & Sons, Ltd.     

A reduced finite difference scheme based on singular value decomposition and proper orthogonal decomposition for Burgers equation

In this article, a reduced optimizing finite difference scheme (FDS) based on singular value decomposition (SVD) and proper orthogonal decomposition (POD) for Burgers equation is presented. Also the error estimates between the usual finite difference solution and the POD solution of reduced optimizing FDS are analyzed. It is shown by considering the results obtained for numerical simulations of cavity flows that the error between the POD solution of reduced optimizing FDS and the solution of the usual FDS is consistent with theoretical results. Moreover, it is also shown that the reduced optimizing FDS is feasible and efficient.     

Monotone scheme for finite difference equations concerning steady-state prey-predator interactions

In this article a system of semilinear elliptic partial differential equations is studied. This system determines the equilibria of the Volterra-Lotka equations describing prey-predator interactions with diffusion. To analyze the system, a new monotone scheme is presented. A rigorous foundation is given for numerical calculations by adapting a suitable finite difference method to the new monotone scheme. Earlier theories in finite differences are not successful in solving the system without this scheme.     

A High-Order Finite Difference Scheme for 3D Unsteady Convection Diffusion Reaction Equations北大核心CSCD

针对三维非稳态对流扩散反应方程,构造了一种高精度紧致有限差分格式,对空间的离散采用四阶紧致差分方法,对时间的离散采用Taylor级数展开和余项修正技术,所提格式在时间上的精度为二阶、在空间上的精度为四阶。利用Fourier稳定性分析法证明了该格式是无条件稳定的。最后给出数值算例验证了理论结果。     

Three-layer finite difference method for solving linear differential algebraic systems of partial differential equations

A boundary value problem is examined for a linear differential algebraic system of partial differential equations with a special structure of the associate matrix pencil. The use of an appropriate transformation makes it possible to split such a system into a system of ordinary differential equations, a hyperbolic system, and a linear algebraic system. A three-layer finite difference method is applied to solve the resulting problem numerically. A theorem on the stability and the convergence of this method is proved, and some numerical results are presented.