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.

COMPACT FINITE DIFFERENCE-FOURIER SPECTRAL METHOD FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

Comparison of finite difference and pseudo-spectral methods in forward modelling based on metal ore model of random media

With more applications of seismic exploration in metal ore exploration,forward modelling of seismic wave has become more important in metal ore. Finite difference method and pseudo-spectral method are two important methods of wave-field simulation. Results of previous studies show that both methods have distinct advantages and disadvantages: Finite difference method has high precision but its dispersion is serious; pseudospectral method considers both computational efficiency and precision but has less precision than finite-difference. The authors consider the complex structural characteristics of the metal ore,furthermore add random media in order to simulate the complex effects produced by metal ore for wave field. First,the study introduced the theories of random media and two forward modelling methods. Second,it compared the simulation results of two methods on fault model. Then the authors established a complex metal ore model,added random media and compared computational efficiency and precision. As a result,it is found that finite difference method is better than pseudo-spectral method in precision and boundary treatment,but the computational efficiency of pseudospectral method is slightly higher than the finite difference method.

Finite Difference Preconditioners for Legendre Based Spectral Element Methods on Elliptic Boundary Value Problems

Finite difference type preconditioners for spectral element discretizations based on Legendre-Gauss-Lobatto points are analyzed. The latter is employed for the approximation of uniformly elliptic partial differential problems. In this work, it is shown that the condition number of the resulting preconditioned system is bounded independently of both of the polynomial degrees used in the spectral element method and the element sizes. Several numerical tests verify the h-p independence of the proposed preconditioning.

Efficient Finite Difference/Spectral Method for the Time Fractional Ito Equation Using Fast Fourier Transform Technic

A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.

Implicit finite difference method for fractional percolation equation with Dirichlet and fractional boundary conditions

An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.

Numerical storm surge model with higher order finite difference method of lines for the coast of Bangladesh

In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.

Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.

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

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

High-Order Bound-Preserving Finite Difference Methods for Multispecies and Multireaction Detonations

In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical reaction,say zi,is between 0 and 1,withΣz_(i)=1.Due to the lack of maximum-principle,most of the previous bound-preserving technique cannot be applied directly.To preserve those bounds,we will use the positivity-preserving technique to all the zi'is and enforceΣz_(i)=1 by constructing conservative schemes,thanks to conservative time integrations and consistent numerical fluxes in the system.Moreover,detonation is an extreme singular mode of flame propagation in premixed gas,and the model contains a significant stiff source.It is well known that for hyperbolic equations with stiff source,the transition points in the numerical approximations near the shocks may trigger spurious shock speed,leading to wrong shock position.Intuitively,the high-order weighted essentially non-oscillatory(WENO)scheme,which can suppress oscillations near the discontinuities,would be a good choice for spatial discretization.However,with the nonlinear weights,the numerical fluxes are no longer“consistent”,leading to nonconservative numerical schemes and the bound-preserving technique does not work.Numerical experiments demonstrate that,without further numerical techniques such as subcell resolutions,the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.

Sloshing simulation of standing wave with time-independent finite difference method for Euler equations

The numerical solutions of standing waves for Euler equations with the nonlinear free surface boundary condition in a two-dimensional （2D） tank are studied. The irregular tank is mapped onto a fixed square domain through proper mapping functions. A staggered mesh system is employed in a 2D tank to calculate the elevation of the transient fluid. A time-independent finite difference method, which is developed by Bang- fuh Chen, is used to solve the Euler equations for incompressible and inviscid fluids. The numerical results agree well with the analytic solutions and previously published results. The sloshing profiles of surge and heave motion with initial standing waves are presented. The results show very clear nonlinear and beating phenomena.

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.

Finite Difference Methods for the Time Fractional Advection-diffusion Equation

In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.

Second-Order Finite Difference/Spectral Element Formulation for Solving the Fractional Advection-Diffusion Equation

The main aim of this paper is to analyze the numerical method based upon the spectral element technique for the numerical solution of the fractional advection-diffusion equa-tion.The time variable has been discretized by a second-order finite difference procedure.The stability and the convergence of the semi-discrete formula have been proven.Then,the spatial variable of the main PDEs is approximated by the spectral element method.The convergence order of the fully discrete scheme is studied.The basis functions of the spectral element method are based upon a class of Legendre polynomials.The numerical experiments confirm the theoretical results.

An efficient locally one-dimensional finite-difference time-domain method based on the conformal scheme

An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.

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.

THE DYNAMICAL BEHAVIOR OF FULLY DISCRETE SPECTRAL METHOD FOR NONLINEAR SCHRODINGER EQUATION WITH WEAKLY DAMPED

Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference to

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.

Analysis of an Implicit Finite Difference Scheme for Time Fractional Diffusion Equation

Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.

Numerical Analysis of Upwind Difference Schemes for Two-Dimensional First-Order Hyperbolic Equations with Variable Coefficients

In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.