A high‐order compact boundary value method for solving one‐dimensional heat equations

We combine fourth‐order boundary value methods (BVMs) for discretizing the temporal variable with fourth‐order compact difference scheme for discretizing the spatial variable to solve one‐dimensional heat equations. This class of new compact difference schemes achieve fourth‐order accuracy in both temporal and spatial variables and are unconditionally stable due to the favorable stability property of BVMs. Numerical results are presented to demonstrate the accuracy and efficiency of the new compact difference scheme, compared to the standard second‐order Crank‐Nicolson scheme. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 846–857, 2003.     

A finite difference scheme for solving parabolic two‐step micro‐heat transport equations in a double‐layered micro‐sphere heated by ultrashort‐pulsed lasers

Ultrashort‐pulsed lasers with pulse durations of the order of sub‐picosecond to femtosecond domain possess exclusive capabilities in limiting the undesirable spread of the thermal process zone in the heated sample. Parabolic two‐step micro heat transport equations have been widely applied for thermal analysis of thin metal films exposed to picosecond thermal pulses. In this study, we develop a three‐level finite difference scheme for solving the micro heat transport equations in a double‐layered micro sphere. It is shown by the discrete energy method that the scheme is unconditionally stable. Numerical results for thermal analysis of a gold layer coated on a chromium padding layer are obtained. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A high‐order ADI finite difference scheme for a 3D reaction‐diffusion equation with neumann boundary condition

In this article, we extend the fourth‐order compact boundary scheme in Liao et al. (Numer Methods Partial Differential Equations 18 (2002), 340–354) to a 3D problem and then combine it with the fourth‐order compact alternating direction implicit (ADI) method in Gu et al. (J Comput Appl Math 155 (2003), 1–17) to solve the 3D reaction‐diffusion equation with Neumann boundary condition. First, the reaction‐diffusion equation is solved with a compact fourth‐order finite difference method based on the Padé approximation, which is then combined with the ADI method and a fourth‐order compact scheme to approximate the Neumann boundary condition, to obtain fourth order accuracy in space. The accuracy in the temporal dimension is improved to fourth order by applying the Richardson extrapolation technique, although the unconditional stability of the numerical method is proved, and several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed new algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Adequacy of finite difference schemes for convection‐diffusion equations

Finite difference schemes for the numerical solution of singularly perturbed convection problems on uniform grids are studied in the limit case where the viscosity and the meshsize approach zero at the same time. The present error estimates are given in terms of order of magnitude in the above limit process and are useful in a priori choosing adequate schemes and meshsizes for boundary‐layer problems and problems with closed characteristics. Published 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 280–295, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10007     

Technical note: A fourth‐order finite difference scheme for a system of a 2D nonlinear elliptic partial differential equations

In this article we present a fourth‐order finite difference scheme, for a system of two‐dimensional, second‐order, nonlinear elliptic partial differential equations with mixed spatial derivative terms, using 13‐point stencils with a uniform mesh size h on a square region R subject to Dirichlet boundary conditions. The scheme of order h⁴ is derived using the local solution of the system on a single stencil. The resulting system of algebraic equations can be solved by iterative methods. The difference scheme can be easily modified to obtain formulae for grid points near the boundary. Computational results are given to demonstrate the performance of the scheme on some problems including Navier‐Stokes equations. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 43–53, 2001     

A new fifth order finite difference WENO scheme for Hamilton‐Jacobi equations

In this continuing paper of (Zhu and Qiu, J Comput Phys 318 (2016), 110–121), a new fifth order finite difference weighted essentially non‐oscillatory (WENO) scheme is designed to approximate the viscosity numerical solution of the Hamilton‐Jacobi equations. This new WENO scheme uses the same numbers of spatial nodes as the classical fifth order WENO scheme which is proposed by Jiang and Peng (SIAM J Sci Comput 21 (2000), 2126–2143), and could get less absolute truncation errors and obtain the same order of accuracy in smooth region simultaneously avoiding spurious oscillations nearby discontinuities. Such new WENO scheme is a convex combination of a fourth degree accurate polynomial and two linear polynomials in a WENO type fashion in the spatial reconstruction procedures. The linear weights of three polynomials are artificially set to be any random positive constants with a minor restriction and the new nonlinear weights are proposed for the sake of keeping the accuracy of the scheme in smooth region, avoiding spurious oscillations and keeping sharp discontinuous transitions in nonsmooth region simultaneously. The main advantages of such new WENO scheme comparing with the classical WENO scheme proposed by Jiang and Peng (SIAM J Sci Comput 21 (2000), 2126–2143) are its efficiency, robustness and easy implementation to higher dimensions. Extensive numerical tests are performed to illustrate the capability of the new fifth WENO scheme. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1095–1113, 2017     

A Newton linearized compact finite difference scheme for one class of Sobolev equations

In this article, a Newton linearized compact finite difference scheme is proposed to numerically solve a class of Sobolev equations. The unique solvability, convergence, and stability of the proposed scheme are proved. It is shown that the proposed method is of order 2 in temporal direction and order 4 in spatial direction. Moreover, compare to the classical extrapolated Crank‐Nicolson method or the second‐order multistep implicit–explicit methods, the proposed scheme is easier to be implemented as it only requires one starting value. Finally, numerical experiments on one and two‐dimensional problems are presented to illustrate our theoretical results.     

A study on a second order finite difference scheme for fractional advection–diffusion equations

Second order finite difference schemes for fractional advection–diffusion equations are considered in this paper. We note that, when studying these schemes, advection terms with coefficients having the same sign as those of diffusion terms need additional estimates. In this paper, by comparing generating functions of the corresponding discretization matrices, we find that sufficiently strong diffusion can dominate the effects of advection. As a result, convergence and stability of schemes are obtained in this situation.     

High‐order compact finite difference and laplace transform method for the solution of time‐fractional heat equations with dirchlet and neumann boundary conditions

The work presents a novel coupling of the Laplace Transform and the compact fourth‐order finite‐difference discretization scheme for the efficient and accurate solution of linear time‐fractional nonhomogeneous diffusion equations subject to both Dirichlet and Neumann boundary conditions. A translational transformation of the dependent variable ensures the Caputo derivative is aligned with the Riemann‐Louiville fractional derivative. The resulting scheme is computationally efficient and shown to be uniquely solvable in all cases, accurate and convergent to in the spatial domain. The convergence rates in the temporal domain are contour dependent but exhibit geometric convergence. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1184–1199, 2016     

A note on second‐order finite‐difference schemes on uniform meshes for advection‐‐diffusion equations

An artificial‐viscosity finite‐difference scheme is introduced for stabilizing the solutions of advection‐diffusion equations. Although only the linear one‐dimensional case is discussed, the method is easily susceptible to generalization. Some theory and comparisons with other well‐known schemes are carried out. The aim is, however, to explain the construction of the method, rather than considering sophisticated applications. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 581–588, 1999     

Nonstandard finite difference schemes for reaction‐diffusion equations

We construct finite difference schemes for a particular class of one‐space dimension, nonlinear reaction‐diffusion PDEs. The use of nonstandard finite difference methods and the imposition of a positivity condition constrain the schemes to be explicit and allow the determination of functional relations between the space and time step‐sizes. The general procedure is illustrated by applying it to several important model systems of PDEs © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 201–214, 1999     

A fourth‐order compact finite difference scheme for solving an N‐carrier system with Neumann boundary conditions

In this study, we develop a fourth‐order compact finite difference scheme for solving a model of energy exchanges in a generalized N‐carrier system with heat sources and Neumann boundary conditions, which extends the concept of the well‐known parabolic two‐step model for microheat transfer. By using the matrix analysis, the compact finite difference numerical scheme is shown to be unconditionally stable. The accuracy of the solution obtained by the scheme is tested by a numerical example. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A second‐order finite difference approximation for a mathematical model of erythropoiesis

We present a second‐order finite difference scheme for approximating solutions of a mathematical model of erythropoiesis, which consists of two nonlinear partial differential equations and one nonlinear ordinary differential equation. We show that the scheme achieves second‐order accuracy for smooth solutions. We compare this scheme to a previously developed first‐order method and show that the first order method requires significantly more computational time to provide solutions with similar accuracy. We also compare this numerical scheme with other well‐known second‐order methods and show that it has better capability in approximating discontinuous solutions. Finally, we present an application to recovery after blood loss. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A high‐order finite difference discretization strategy based on extrapolation for convection diffusion equations

We propose a new high‐order finite difference discretization strategy, which is based on the Richardson extrapolation technique and an operator interpolation scheme, to solve convection diffusion equations. For a particular implementation, we solve a fine grid equation and a coarse grid equation by using a fourth‐order compact difference scheme. Then we combine the two approximate solutions and use the Richardson extrapolation to compute a sixth‐order accuracy coarse grid solution. A sixth‐order accuracy fine grid solution is obtained by interpolating the sixth‐order coarse grid solution using an operator interpolation scheme. Numerical results are presented to demonstrate the accuracy and efficacy of the proposed finite difference discretization strategy, compared to the sixth‐order combined compact difference (CCD) scheme, and the standard fourth‐order compact difference (FOC) scheme. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 18–32, 2004.     

Mixed two‐grid finite difference methods for solving one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations

The aim of this paper is to propose mixed two‐grid finite difference methods to obtain the numerical solution of the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. The finite difference equations at all interior grid points form a large‐sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a family of finite difference methods for discretizing the spatial and time derivatives. The obtained system has been solved by two‐grid method, where the two‐grid method is used for solving the large‐sparse linear systems. Also, in the proposed method, the spectral radius with local Fourier analysis is calculated for different values of h and Δt. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional and two‐dimensional Fitzhugh–Nagumo equations. Copyright © 2016 John Wiley & Sons, Ltd.     

Gidas-Ni-Nirenberg results for finite difference equations: Estimates of approximate symmetry

Are positive solutions of finite difference boundary value problems Δhu=f(u) in Ωh, u=0 on ∂Ωh as symmetric as the domain? To answer this question we first show by examples that almost arbitrary non-symmetric solutions can be constructed. This is in striking difference to the continuous case, where by the famous Gidas-Ni-Nirenberg theorem [B. Gidas, Wei-Ming Ni, L. Nirenberg, Symmetry and related problems via the maximum principle, Comm. Math. Phys. 68 (1979) 209-243] positive solutions inherit the symmetry of the underlying domain. Then we prove approximate symmetry theorems for solutions on equidistantly meshed n-dimensional cubes: explicit estimates depending on the data are given which show that the solutions become more symmetric as the discretization gets finer. The quality of the estimates depends on whether or not f(0)<0. The one-dimensional case stands out in two ways: the proofs are elementary and the estimates for the defect of symmetry are O(h) compared to O(1/|log(h)|) in the higher-dimensional case.     

Nonstandard finite difference schemes for Michaelis–Menten type reaction‐diffusion equations

We compare and investigate the performance of the exact scheme of the Michaelis–Menten (M–M) ordinary differential equation with several new nonstandard finite difference (NSFD) schemes that we construct using Mickens' rules. Furthermore, the exact scheme of the M–M equation is used to design several dynamically consistent NSFD schemes for related reaction‐diffusion equations, advection‐reaction equations, and advection‐reaction‐diffusion equations. Numerical simulations that support the theory and demonstrate computationally the power of NSFD schemes are presented. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A multigrid compact finite difference method for solving the one‐dimensional nonlinear sine‐Gordon equation

The aim of this paper is to propose a multigrid method to obtain the numerical solution of the one‐dimensional nonlinear sine‐Gordon equation. The finite difference equations at all interior grid points form a large sparse linear system, which needs to be solved efficiently. The solution cost of this sparse linear system usually dominates the total cost of solving the discretized partial differential equation. The proposed method is based on applying a compact finite difference scheme of fourth‐order for discretizing the spatial derivative and the standard second‐order central finite difference method for the time derivative. The proposed method uses the Richardson extrapolation method in time variable. The obtained system has been solved by V‐cycle multigrid (VMG) method, where the VMG method is used for solving the large sparse linear systems. The numerical examples show the efficiency of this algorithm for solving the one‐dimensional sine‐Gordon equation. Copyright © 2014 John Wiley & Sons, Ltd.     

A finite difference scheme for solving the heat transport equation at the microscale

Heat transport at the microscale is of vital importance in microtechnology applications. In this study, we develop a finite difference scheme of the Crank‐Nicholson type by introducing an intermediate function for the heat transport equation at the microscale. It is shown by the discrete energy method that the scheme is unconditionally stable. Numerical results show that the solution is accurate. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 697–708, 1999     

High‐order difference scheme for the solution of linear time fractional klein–gordon equations

In this article, we apply a high‐order difference scheme for the solution of some time fractional partial differential equations (PDEs). The time fractional Cattaneo equation and the linear time fractional Klein–Gordon and dissipative Klein–Gordon equations will be investigated. The time fractional derivative which has been described in the Caputo's sense is approximated by a scheme of order , and the space derivative is discretized with a fourth‐order compact procedure. We will prove the solvability of the proposed method by coefficient matrix property and the unconditional stability and ‐convergence with the energy method. Numerical examples demonstrate the theoretical results and the high accuracy of the proposed scheme. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1234–1253, 2014