Perfect plane-wave source for a high-order symplectic finite-difference time-domain scheme

The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.

High-order targeted essentially non-oscillatory scheme for two-fluid plasma model

The weakly ionized plasma flows in aerospace are commonly simulated by the single-fluid model,which cannot describe certain nonequilibrium phenomena by finite collisions of particles,decreasing the fidelity of the solution.Based on an alternative formulation of the targeted essentially non-oscillatory(TENO)scheme,a novel high-order numerical scheme is proposed to simulate the two-fluid plasmas problems.The numerical flux is constructed by the TENO interpolation of the solution and its derivatives,instead of being reconstructed from the physical flux.The present scheme is used to solve the two sets of Euler equations coupled with Maxwell's equations.The numerical methods are verified by several classical plasma problems.The results show that compared with the original TENO scheme,the present scheme can suppress the non-physical oscillations and reduce the numerical dissipation.

High-Order Semi-Lagrangian WENO Schemes Based on Non-polynomial Space for the Vlasov Equation

In this paper,we present a semi-Lagrangian(SL)method based on a non-polynomial function space for solving the Vlasov equation.We fnd that a non-polynomial function based scheme is suitable to the specifcs of the target problems.To address issues that arise in phase space models of plasma problems,we develop a weighted essentially non-oscillatory(WENO)scheme using trigonometric polynomials.In particular,the non-polynomial WENO method is able to achieve improved accuracy near sharp gradients or discontinuities.Moreover,to obtain a high-order of accuracy in not only space but also time,it is proposed to apply a high-order splitting scheme in time.We aim to introduce the entire SL algorithm with high-order splitting in time and high-order WENO reconstruction in space to solve the Vlasov-Poisson system.Some numerical experiments are presented to demonstrate robustness of the proposed method in having a high-order of convergence and in capturing non-smooth solutions.A key observation is that the method can capture phase structure that require twice the resolution with a polynomial based method.In 6D,this would represent a signifcant savings.

A staggered-grid high-order finite-difference modeling for elastic wave field in arbitrary tilt anisotropic media

The paper presents a staggered-grid any even-order accurate finite-difference scheme for two-dimensional (2D), three-component (3C), first-order stress-velocity elastic wave equation and its stability condition in the arbitrary tilt anisotropic media; and derives a perfectly matched absorbing layer (PML) boundary condition and its stag- gered-grid any even-order accurate difference scheme in the 2D arbitrary tilt anisotropic media. The results of nu- merical modeling indicate that the modeling precision is high, the calculation efficiency is satisfactory and the absorbing boundary condition is better. The wave-front shapes of elastic waves are complex in the anisotropic media, and the velocity of qP wave is not always faster than that of qS wave. The wave-front triplication of qS wave and its events in both reflected domain and propagated domain, which are not commonly hyperbola, is a common phenomenon. When the symmetry axis is tilted in the TI media, the phenomenon of S-wave splitting is clearly observed in the snaps of three components and synthetic seismograms, and the events of all kinds of waves are asymmetric.

A high-order splitting scheme for the advection diffusion equation

A high-order splitting scheme for the advection-diffusion equation of pollutants is proposed in this paper. The multidimensional advection-diffusion equation is splitted into several one-dimensional equations that are solved by the scheme. Only three spatial grid points are needed in each direction and the scheme has fourth-order spatial accuracy. Several typically pure advection and advection-diffusion problems are simulated. Numerical results show that the accuracy of the scheme is much higher than that of the classical schemes and the scheme can he efficiently solved with little programming effort.

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 High-Order Spatiotemporal Precision-Matching Taylor–Li Scheme for Time-Dependent Problems

Based on the Taylor series method and Li’s spatial differential method, a high-order hybrid Taylor–Li scheme is proposed.The results of a linear advection equation indicate that, using the initial values of the square-wave type, a result with thirdorder accuracy occurs. However, using initial values associated with the Gaussian function type, a result with very high precision appears. The study demonstrates that, when the order of the time integral is more than three, the corresponding optimal spatial difference order could be higher than six. The results indicate that the reason for why there is no improvement related to an order of spatial difference above six is the use of a time integral scheme that is not high enough. The author also proposes a recursive differential method to improve the Taylor–Li scheme’s computation speed. A more rapid and highprecision program than direct computation of the high-order space differential item is employed, and the computation speed is dramatically boosted. Based on a multiple-precision library, the ultrahigh-order Taylor–Li scheme can be used to solve the advection equation and Burgers’ equation.

High-order schemes for predicting computational aeroacoustic propagation with adaptive mesh refinement

High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical issues associated with high-order schemes on an adaptively refined mesh, such as stability and accuracy are addressed. Several CAA benchmark problems are used to demonstrate the feasibility and efficiency of the approach.

Convergence to Steady-State Solutions of the New Type of High-Order Multi-resolution WENO Schemes: a Numerical Study

A new type of high-order multi-resolution weighted essentially non-oscillatory(WENO)schemes(Zhu and Shu in J Comput Phys,375:659-683,2018)is applied to solve for steady-state problems on structured meshes.Since the classical WENO schemes(Jiang and Shu in J Comput Phys,126:202-228,1996)might suffer from slight post-shock oscillations(which are responsible for the residue to hang at a truncation error level),this new type of high-order finite-difference and finite-volume multi-resolution WENO schemes is applied to control the slight post-shock oscillations and push the residue to settle down to machine zero in steady-state simulations.This new type of multi-resolution WENO schemes uses the same large stencils as that of the same order classical WENO schemes,could obtain fifth-order,seventh-order,and ninth-order in smooth regions,and could gradually degrade to first-order so as to suppress spurious oscillations near strong discontinuities.The linear weights of such new multi-resolution WENO schemes can be any positive numbers on the condition that their sum is one.This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finitedifference and finite-volume WENO schemes for solving steady-state problems.In comparison with the classical fifth-order finite-difference and finite-volume WENO schemes,the residue of these new high-order multi-resolution WENO schemes can converge to a tiny number close to machine zero for some benchmark steady-state problems.

A High-Order Scheme for Fractional Ordinary Differential Equations with the Caputo-Fabrizio Derivative

In this paper, we consider numerical solutions of fractional ordinary diferential equations with the Caputo-Fabrizio derivative, and construct and analyze a high-order time-stepping scheme for this equation. The proposed method makes use of quadratic interpolation function in sub-intervals, which allows to produce fourth-order convergence. A rigorous stability and convergence analysis of the proposed scheme is given. A series of numerical examples are presented to validate the theoretical claims. Traditionally a scheme having fourth-order convergence could only be obtained by using block-by-block technique. The advantage of our scheme is that the solution can be obtained step by step, which is cheaper than a block-by-block-based approach.

High-Order Local Discontinuous Galerkin Algorithm with Time Second-Order Schemes for the Two-Dimensional Nonlinear Fractional Diffusion Equation

In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.

The Best Finite-Difference Scheme for the Helmholtz Equation

The best finite-difference scheme for the Helmholtz equation is suggested. A method of solving obtained finite-difference scheme is developed. The efficiency and accuracy of method were tested on several examples.

Optimization of a global seventh-order dissipative compact finite-difference scheme by a genetic algorithm

A global seventh-order dissipative compact finite-difference scheme is optimized in terms of time stability. The dissipative parameters appearing in the boundary closures are assumed to be different, resulting in an optimization problem with several parameters determined by applying a generic algorithm. The optimized schemes are analyzed carefully from the aspects of the eigenvalue distribution, the ε-pseudospectra, the short time behavior, and the Fourier analysis. Numerical experiments for the Euler equations are used to show the effectiveness of the final recommended scheme.

Properties of High-Order Finite Difference Schemes and Idealized Numerical Testing

Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.

A FAMILY OF HIGH-ORDER ACCURACY EXPLICIT DIFFERENCE SCHEMES WITH BRANCHING STABILITY FOR SOLVING 3-D PARABOLIC PARTIAL DIFFERENTIAL EQUATION

A family of high-order accuracy explict difference schemes for solving 3-dimension parabolic P. D. E. is constructed. The stability condition is r = Deltat/Deltax(2) Deltat/Deltay(2) = Deltat/Deltaz(2) < 1/2 ,and the truncation error is 0(<Delta>t(2) + Deltax(4)).

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.

Stability of High-Order Staggered-Grid Schemes for 3D Elastic Wave Equation in Heterogeneous Media

In this paper, we firstly derive the stability conditions of high-order staggered-grid schemes for the three-dimensional (3D) elastic wave equation in heterogeneous media based on the energy method. Moreover, the plane wave analysis yields a sufficient and necessary stability condition by the von Neumann criterion in homogeneous case. Numerical computations for 3D wave simulation with point source excitation are given.

SEVERAL NEW TYPES OF FINITE-DIFFERENCE SCHEMES FOR SHALLOW-WATER EQUATION WITH INITIAL-BOUNDARY VALUE AND THEIR NUMERICAL EXPERIMENT

This paper proposed several new types of finite-difference methods for the shallow water equation in absolute coordinate system and put forward an effective two-step predictor-corrector method, a compact and iterative algorithm for five diagonal matrix. Then the iterative method was used for a multi-grid procedure for shallow water equation. A t last, an initial-boundary value problem was considered, and the numerical results show that the linear sinusoidal wave would successively evolve into conoidal wave.

A NEW HIGH-ORDER ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING THREE-DIMENSIONAL PARABOLIC EQUATIONS

In this paper, a new three-level explicit difference scheme with high-order accuracy is proposed for solving three-dimensional parabolic equations. The stability condition is r = Delta t/Delta x(2) = Delta t/Delta gamma(2) = Delta t/Delta z(2) less than or equal to 1/4, and the truncation error is O(Delta t(2) + Delta x(4)).

High-order discontinuous Galerkin method for applications to multicomponent and chemically reacting flows

This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been successfully employed, DG simulations of chemically reacting flows introduce challenges that arise from flow unsteadiness, combustion, heat release, compressibility effects, shocks, and variations in thermodynamic properties. To address these challenges, algorithms are developed, including an entropy-bounded DG method, an entropy-residual shock indicator, and a new formulation of artificial viscosity. The performance and capabilities of the resulting DG method are demonstrated in several relevant applications, including shock/bubble interaction, turbulent combustion, and detonation. It is concluded that the developed DG method shows promising performance in application to multicomponent reacting flows. The paper concludes with a discussion of further research needs to enable the application of DG methods to more complex reacting flows.