Efficient Finite Difference WENO Scheme for Hyperbolic Systems withNon-conservativeProducts

Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of the comparable accuracy.This makes them quite attractive for several science and engineering applications.But,to the best of our knowledge,such schemes have not been extended to non-linear hyperbolic systems with non-conservative products.In this paper,we perform such an extension which improves the domain of the applicability of such schemes.The extension is carried out by writing the scheme in fluctuation form.We use the HLLI Riemann solver of Dumbser and Balsara(J.Comput.Phys.304:275-319,2016)as a building block for carrying out this extension.Because of the use of an HLL building block,the resulting scheme has a proper supersonic limit.The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme,thus expanding its domain of the applicability.Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions,making it very easy for users to transition over to the present formulation.For conservation laws,the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO,with two major advantages:(i)It can capture jumps in stationary linearly degenerate wave families exactly.(i)It only requires the reconstruction to be applied once.Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of the accuracy for smooth multidimensional flows.Stringent Riemann problems and several novel multidimensional problems that are drawn from compressible Baer-Nunziato multiphase flow,multiphase debris flow and twolayer shallow water equations are also shown to document the robustness of the method.For some test problems that require well-balancing we have even been able to apply the scheme without any modification and obtain good results.Many useful PDEs may have stiff relaxation source terms for which the finite difference formulation of WENO is shown to provide some genuine advantages.

THE FINITE DIFFERENCE STREAMLINE DIFFUSION METHODS FOR TIME-DEPENDENT CONVECTION-DIFFUSION EQUATIONS

In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.

A-high-order Accuraqcy Implicit Difference Scheme for Solving the Equation of Parabolic Type

In this paper,a implicit difference scheme is proposed for solving the equation of one_dimension parabolic type by undetermined paameters.The stability condition is r=αΔt/Δx 2 1/2 and the truncation error is o(Δt 4+Δx 4) It can be easily solved by double sweeping method.

CONSERVATIVE DIFFERENCE SCHEME BASED ON NUMERICAL ANALYSIS FOR NONLINEAR SCHRDINGER EQUATION WITH WAVE OPERATOR

A new conservative finite difference scheme is presented based on the numerical analysis for an initialboundary value problem of a class of Schroedinger equation with the wave operator. The scheme can be linear and implicit or explicit based on the parameter choice. The initial value after discretization has second-order accuracy that is consistent with the scheme accuracy. The existence and the uniqueness of the difference solution are proved. Based on the priori estimates and an inequality about norms, the stability and the convergence of difference solutions with the second-order are proved in the energy norm. Experimental results demonstrate the efficiency of the new scheme.

An Explicit Difference Scheme with High Accuracy and Branching Stability for Solving Parabolic Partial Differential Equation

This paper presents an explicit difference scheme with accuracy and branching stability for solving onedimensional parabolic type equation by the method of undetermined parameters and its truncation error is O(△t4+△x4). The stability condition is r=a△t/△x2<1/2.

A Family of High_order Accuracy Explicit Difference Schemes for Solving 2-D Parabolic Partial Differential Equation

A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx 4).

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.

Stability of Difference Schemes with Intrinsic Parallelism for Quasilinear Parabolic Systems

In this paper, the first boundary problem of quasilinear parabolic system of second order is studied by the finite difference method with intrinsic parallelism. for the problem, the stability of the difference schemes with intrinsic parallelism are justified in the sense of the continuous dependence of the discrete vector solution of the difference schemes on the discrete data of the original problem, without assuming the existence of the smooth solutions for the origillal problem.

A CONSERVATIVE DIFFERENCE SCHEME FOR CONSERVATIVE DIFFERENTIAL EQUATION WITH PERIODIC BOUNDARY

The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, it is proved that the scheme converges uniformly to the solution of differential equation with order one.

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 CLASS OF COMPACT UPWIND TVD DIFFERENCE SCHEMES

A new method was proposed for constructing total variation diminishing （TVD） upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can ensure the nonlinear compact schemes TVD property. Two compact TVD （CTVD） schemes were tested, one is thirdorder accuracy, and the other is fifth-order. The performance of the numerical algorithms was assessed by one-dimensional complex waves and Riemann problems, as well as a twodimensional shock-vortex interaction and a shock-boundary flow interaction. Numerical results show their high-order accuracy and high resolution, and low oscillations across discontinuities.

Second-order difference scheme for a nonlinear model of wood drying process

A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.

Analysis of Extended Fisher-Kolmogorov Equation in 2D Utilizing the Generalized Finite Difference Method with Supplementary Nodes

In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.

A High-order Accuracy Explicit Difference Scheme with Branching Stability for Solving Higher-dimensional Heat-conduction Equation

A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncation error is O（△t^2 ＋ △x^4）.

A SIMPLE PROOF OF STABILITY AND CONVERGENCE IN L＿2 FOR SOME DIFFERENCE SCHEMES AND THEIR EXTRAPOLATION METHOD FOR PARABOLIC EQUATIONS

ＡＳＩＭＰＬＥＰＲＯＯＦＯＦＳＴＡＢＩＬＩＴＹＡＮＤＣＯＮＶＥＲＧＥＮＣＥＩＮＬ＿２ＦＯＲＳＯＭＥＤＩＦＦＥＲＥＮＣＥＳＣＨＥＭＥＳＡＮＤＴＨＥＩＲＥＸＴＲＡＰＯＬＡＴＩＯＮＭＥＴＨＯＤＦＯＲＰＡＲＡＢＯＬＩＣＥＱＵＡＴＩＯＮＳＳｕｎＺｈｉｚｈｏｎｇ...

Improving the Stability Problem of the Finite Difference Scheme for Reaction-diffusion Equation

This paper deals with the special nonlinear reaction-diffusion equation. The finite difference scheme with incremental unknowns approximating to the differential equation （2.1） is set up by means of introducing incremental unknowns methods. Through the stability analyzing for the scheme, it was shown that the stability conditions of the finite difference schemes with the incremental unknowns are greatly improved when compared with the stability conditions of the corresponding classic difference scheme.

New optimized flux difference schemes for improving high-order weighted compact nonlinear scheme with applications

To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal orders,or extending stencil widths,are avoided in the new optimized schemes by utilizing fluxes from both cell-edges and cell-nodes.Optimizations are implemented with Fourier analysis for linear schemes and the approximate dispersion relation(ADR)for nonlinear schemes.Classical difference schemes are restored near discontinuities to suppress numerical oscillations with use of a shock sensor based on smoothness indicators.The results of several benchmark numerical tests indicate that the new optimized difference schemes outperform the classical schemes,in terms of accuracy and resolution for smooth wave and vortex,especially for long-time simulations.Using optimized schemes increases the total CPU time by less than 4%.

Computational Stability of the Explicit Difference Schemes of the Forced Dissipative Nonlinear Evolution Equations

The computational stability of the explicit difference schemes of the forced dissipative nonlinear evolution equations is analyzed and the computational quasi-stability criterion of explicit difference schemes of the forced dissipative nonlinear atmospheric equations is obtained on account of the concept of computational quasi-stability, Therefore, it provides the new train of thought and theoretical basis for designing computational stable difference scheme of the forced dissipative nonlinear atmospheric equations. Key words Computational quasi-stability - Computational stability - Forced dissipative nonlinear evolution equation - Explicit difference scheme This work was supported by the National Outstanding Youth Scientist Foundation of China (Grant No. 49825109), the Key Innovation Project of Chinese Academy of Sciences (KZCX1-10-07), the National Natural Science Foundation of China (Grant Nos, 49905007 and 49975020) and the Outstanding State Key Laboratory Project (Grant No. 40023001).

A CLASS OF TWO-LEVEL EXPLICIT DIFFERENCE SCHEMES FOR SOLVING THREE DIMENSIONAL HEAT CONDUCTION EQUATION

A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.

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.