Two Energy-Preserving Compact Finite Difference Schemes for the Nonlinear Fourth-Order Wave Equation

In this paper,two fourth-order compact finite difference schemes are derived to solve the nonlinear fourth-order wave equation which can be viewed as a generalized model from the nonlinear beam equation.Differing from the existing compact finite difference schemes which preserve the total energy in a recursive sense,the new schemes are proved to per-fectly preserve the total energy in the discrete sense.By using the standard energy method and the cut-off function technique,the optimal error estimates of the numerical solutions are established,and the convergence rates are of O(h^(4)+τ^(2))with mesh-size h and time-step τ.In order to improve the computational efficiency,an iterative algorithm is proposed as the outer solver and the double sweep method for pentadiagonal linear algebraic equations is introduced as the inner solver to solve the nonlinear difference schemes at each time step.The convergence of the iterative algorithm is also rigorously analyzed.Several numerical results are carried out to test the error estimates and conservative properties.

Combinations of nonstandard finite difference schemes and composition methods with complex time steps for population models

We propose an efficient numerical method for two population models, based on the nonstandard finite difference （NSFD） schemes and composition methods with complex time steps. The NSFD scheme is able to give positive numerical solutions that satisfy the conservation law, which is a key property for biological population models. The accuracy is improved by using the composition methods with complex time steps. Numerical tests on the plankton nutrient model and whooping cough model are presented to show the efficiency and advantage of the proposed numerical method.

Two finite difference schemes for the phase field crystal equation

The phase field crystal(PFC) model is a nonlinear evolutionary equation that is of sixth order in space.In the first part of this work,we derive a three level linearized difference scheme,which is then proved to be energy stable,uniquely solvable and second order convergent in L_2 norm by the energy method combining with the inductive method.In the second part of the work,we analyze the unique solvability and convergence of a two level nonlinear difference scheme,which was developed by Zhang et al.in 2013.Some numerical results with comparisons are provided.

High-Order Accurate Entropy Stable Finite Difference Schemes for One- and Two-Dimensional Special Relativistic Hydrodynamics

This paper develops the high-order accurate entropy stable finite difference schemes for one-and two-dimensional special relativistic hydrodynamic equations.The schemes are built on the entropy conservative flux and the weighted essentially non-oscillatory(WENO)technique as well as explicit Runge-Kutta time discretization.The key is to technically construct the affordable entropy conservative flux of the semi-discrete second-order accurate entropy conservative schemes satisfying the semi-discrete entropy equality for the found convex entropy pair.As soon as the entropy conservative flux is derived,the dissipation term can be added to give the semidiscrete entropy stable schemes satisfying the semi-discrete entropy inequality with the given convex entropy function.The WENO reconstruction for the scaled entropy variables and the high-order explicit Runge-Kutta time discretization are implemented to obtain the fully-discrete high-order entropy stable schemes.Several numerical tests are conducted to validate the accuracy and the ability to capture discontinuities of our entropy stable schemes.

AWeighted Average Finite Difference Scheme for the Numerical Solution of Stochastic Parabolic Partial Differential Equations

In the present paper,the numerical solution of It?type stochastic parabolic equation with a timewhite noise process is imparted based on a stochastic finite difference scheme.At the beginning,an implicit stochastic finite difference scheme is presented for this equation.Some mathematical analyses of the scheme are then discussed.Lastly,to ascertain the efficacy and accuracy of the suggested technique,the numerical results are discussed and compared with the exact solution.

A Class of Stable and Conservative Finite Difference Schemes for the Cahn-Hilliard Equation

In this paper, we propose a class of stable finite difference schemes for the initial-boundary value problem of the Cahn-Hilliard equation. These schemes are proved to inherit the total mass conservation and energy dissipation in the discrete level. The dissipation of the total energy implies boundness of the numerical solutions in the discrete H1 norm. This in turn implies boundedness of the numerical solutions in the maximum norm and hence the stability of the difference schemes. Unique existence of the numerical solutions is proved by the fixed-point theorem. Convergence rate of the class of finite difference schemes is proved to be O（h2 ＋ r2） with time step T and mesh size h. An efficient iterative algorithm for solving these nonlinear schemes is proposed and discussed in detail.

THE COMPACT SECOND-ORDER UPWIND FINITE DIFFERENCE SCHEMES FOR THE CONVECTION-DIFFUSION EQUATIONS

The compact second-order upwind finite difference schemes free of ceil Reynolds number limitation are developed in this paper for the one-to three-dimensional steady convection- diffusion equations,using a perturbational technique applied to the classical first-order upwind schemes.The present second-order schemes take essentially the same form as those of the first- order schemes,but involve a simple modification to the diffusive coefficients.Numerical exam- ples including one-to three-dimensional model equations of fluid flow and a problem of natural convection with boundary-layer effect are given to illustrate the excellent behavior of the present schemes.

THE REMAINDER－EFFECT ANALYSIS OF FINITE DIFFERENCE SCHEMES AND THE APPLICATIONS

In the present paper two contents are enclosed .First ,the Fourier analysis approach of the dispersion relation and group velocity effect of finite difference schemes is discussed.the defects of the approach is pointed out and the correction is made;Second,a new systematic analysis method -remaider -effect analysis (abbr.REAM)is proposed by means of the modified partial differential equations (abbr MPDE)of finite difference schemes.The analysis is based on the synthetical study of the rational dispersion-and dissipation relations of finite difference schemes.And the method clearly possesses constructivity

HIGH RESOLUTION POSITIVITY-PRESERVING DIFFERENCE SCHEMES FOR TWO DIMENSIONAL EULER EQUATIONS

A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.

THE UPWIND OPERATOR SPLITTING FINITE DIFFERENCE METHOD FOR COMPRESSIBLE TWO-PHASE DISPLACEMENT PROBLEM AND ANALYSIS

For compressible two-phase displacement problem, a kind of upwind operator splitting finite difference schemes is put forward and make use of operator splitting, of calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L 2 norm are derived to determine the error, in the approximate solution.

The Stability Research for the Finite Difference Scheme of a Nonlinear Reaction-diffusion Equation

In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.

A Finite Difference Method and Effective Modification of Gradient Descent Optimization Algorithm for MHD Fluid Flow Over a Linearly Stretching Surface

Present contribution is concerned with the construction and application of a numerical method for the fluid flow problem over a linearly stretching surface with the modification of standard Gradient descent Algorithm to solve the resulted difference equation.The flow problem is constructed using continuity,and Navier Stoke equations and these PDEs are further converted into boundary value problem by applying suitable similarity transformations.A central finite difference method is proposed that gives third-order accuracy using three grid points.The stability conditions of the present proposed method using a Gauss-Seidel iterative procedure is found using Von-Neumann stability criteria and order of the finite difference method is proved by applying the Taylor series on the discretised equation.The comparison of the presently modified optimisation algorithm with the Gauss-Seidel iterative method and standard Newton’s method in optimisation is also made.It can be concluded that the presently modified optimisation Algorithm takes a few iterations to converge with a small value of the parameter contained in it compared with the standard descent algorithm that may take millions of iterations to converge.The present modification of the steepest descent method converges faster than Gauss-Seidel method and standard steepest descent method,and it may also overcome the deficiency of singular hessian arise in Newton’s method for some of the cases that may arise in optimisation problem(s).

Optimized finite difference iterative scheme based on POD technique for 2D viscoelastic wave equation

This study develops an optimized finite difference iterative （OFDI） scheme for the two-dimensional （2D） viscoelastic wave equation. The OFDI scheme is obtained using a proper orthogonal decomposition （POD） method. It has sufficiently high accuracy with very few unknowns for the 2D viscoelastic wave equation. Existence, stability, and convergence of the OFDI solutions are analyzed. Numerical simulations verify efficiency and feasibility of the proposed scheme.

DIFFERENCE SCHEME AND NUMERICAL SIMULATION BASED ON MIXED FINITE ELEMENT METHOD FOR NATURAL CONVECTION PROBLEM

The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.

Reduced finite difference scheme and error estimates based on POD method for non-stationary Stokes equation

The proper orthogonal decomposition （POD） is a model reduction technique for the simulation Of physical processes governed by partial differential equations （e.g., fluid flows）. It has been successfully used in the reduced-order modeling of complex systems. In this paper, the applications of the POD method are extended, i.e., the POD method is applied to a classical finite difference （FD） scheme for the non-stationary Stokes equation with a real practical applied background. A reduced FD scheme is established with lower dimensions and sufficiently high accuracy, and the error estimates are provided between the reduced and the classical FD solutions. Some numerical examples illustrate that the numerical results are consistent with theoretical conclusions. Moreover, it is shown that the reduced FD scheme based on the POD method is feasible and efficient in solving the FD scheme for the non-stationary Stokes equation.

Convergence on Finite Difference Solution for Semilinear Wave Equation in One Space Variable

In this paper,the convergence and stability of the ’Leap-frog’ finite difference scheme for the semilinear wave equation are proved by using of the bounded extensive method under more generalized condition for the nonlinear term. The more complex standard a priori estimates are avoided so that the theoretical results are complemented for the scheme which was presented by Perring and Skyrne (1962).

LARGE TIME STEP GENERALIZATION OF RANDOM CHOICE FINITE DIFFERENCE SCHEME FOR HYPERBOLIC CONSERVATION LAWS

A natural generalization of random choice finite difference scheme of Harten and Lax for Courant number larger than 1 is obtained. We handle interactions between neighboring Riemann solvers by linear superposition of their conserved quantities. We show consistency of the scheme for arbitrarily large Courant numbers. For scalar problems the scheme is total variation diminishing.A brief discussion is given for entropy condition.

Finite Difference Approach for Estimating the Thermal Conductivity by 6-point Crank-Nicolson Scheme

Based on inverse heat conduction theory, a theoretical model using 6-point Crank-Nicolson finite difference scheme was used to calculate the thermal conductivity from temperature distribution, which can be measured experimentally. The method is a direct approach of second-order and the key advantage of the present method is that it is not required a priori knowledge of the functional form of the unknown thermal conductivity in the calculation and the thermal parameters are estimated only according to the known temperature distribution. Two cases were numerically calculated and the influence of experimental deviation on the precision of this method was discussed. The comparison of numerical and analytical results showed good agreement.

Numerical study of convective heat transfer in static arrangements of particles with arbitrary shapes:A monolithic hybrid lattice Boltzmann-finite difference-phase field solver

A compressible lattice Boltzmann-finite difference method is extended by the phase-field approach into a monolithic scheme to study fluid flow and heat transfer through regular arrangements of solid bodies of circular,elliptical and irregular shapes.The advantage of using the phase-field method is demon-strated both in its simplicity of accounting for flow and thermal boundary conditions at solid surfaces with irregular shapes and in the capability of generating such complex-shaped objects.For an array of discs,numerical results for the overall solid-to-gas heat transfer rate are validated via experiments on flow through arrays of hot cylinders.The thus validated compressible LB-FD-PF hybrid scheme is used to study the dependence of heat transfer on flow and thermal boundary conditions(Reynolds number,temperature difference between the hot solid bodies and the inlet gas),porosity as well as on the shape of solid objects.Results are rationalized in terms of the residence time of the gas close to the solid body and downstream variations of gas velocity and temperature.Perspective for further applications of the proposed methodology are also discussed.

High-order maximum-principle-preserving and positivity-preserving weighted compact nonlinear schemes for hyperbolic conservation laws

In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.