Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implic...Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implicit computations and full precision for the explicit computations.In this work,we analyze the stability properties of these methods and their sensitivity to the low-precision rounding errors,and demonstrate their performance in terms of accuracy and efficiency.We develop codes in FORTRAN and Julia to solve nonlinear systems of ODEs and PDEs using the mixed-precision additive Runge-Kutta(MP-ARK)methods.The convergence,accuracy,and runtime of these methods are explored.We show that for a given level of accuracy,suitably chosen MP-ARK methods may provide significant reductions in runtime.展开更多
This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either...This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.展开更多
In this paper we develop explicit fast exponential Runge-Kutta methods for the numerical solutions of a class of parabolic equations.By incorporating the linear splitting technique into the explicit exponential Runge-...In this paper we develop explicit fast exponential Runge-Kutta methods for the numerical solutions of a class of parabolic equations.By incorporating the linear splitting technique into the explicit exponential Runge-Kutta schemes,we are able to greatly improve the numerical stability.The proposed numerical methods could be fast implemented through use of decompositions of compact spatial difference operators on a regular mesh together with discrete fast Fourier transform techniques.The exponential Runge-Kutta schemes are easy to be adopted in adaptive temporal approximations with variable time step sizes,as well as applied to stiff nonlinearity and boundary conditions of different types.Linear stabilities of the proposed schemes and their comparison with other schemes are presented.We also numerically demonstrate accuracy,stability and robustness of the proposed method through some typical model problems.展开更多
This paper generalizes the exponential Runge-Kutta asymptotic preserving(AP)method developed in[G.Dimarco and L.Pareschi,SIAM Numer.Anal.,49(2011),pp.2057–2077]to compute the multi-species Boltzmann equation.Compared...This paper generalizes the exponential Runge-Kutta asymptotic preserving(AP)method developed in[G.Dimarco and L.Pareschi,SIAM Numer.Anal.,49(2011),pp.2057–2077]to compute the multi-species Boltzmann equation.Compared to the single species Boltzmann equation that the method was originally applied on,this set of equation presents a new difficulty that comes from the lack of local conservation laws due to the interaction between different species.Hence extra stiff nonlinear source terms need to be treated properly to maintain the accuracy and the AP property.The method we propose does not contain any nonlinear nonlocal implicit solver,and can capture the hydrodynamic limit with time step and mesh size independent of the Knudsen number.We prove the positivity and strong AP properties of the scheme,which are verified by two numerical examples.展开更多
In this work,the exponential approximation is used for the numerical simulation of a nonlinear SITR model as a system of differential equations that shows the dynamics of the new coronavirus(COVID-19).The SITR mathema...In this work,the exponential approximation is used for the numerical simulation of a nonlinear SITR model as a system of differential equations that shows the dynamics of the new coronavirus(COVID-19).The SITR mathematical model is divided into four classes using fractal parameters for COVID-19 dynamics,namely,susceptible(S),infected(I),treatment(T),and recovered(R).The main idea of the presented method is based on the matrix representations of the exponential functions and their derivatives using collocation points.To indicate the usefulness of this method,we employ it in some cases.For error analysis of the method,the residual of the solutions is reviewed.The reported examples show that the method is reasonably efficient and accurate.展开更多
The reliability assessment of unit-system near two levels is the mostimportant content in the reliability multi-level synthesis of complex systems. Introducing theinformation theory into system reliability assessment,...The reliability assessment of unit-system near two levels is the mostimportant content in the reliability multi-level synthesis of complex systems. Introducing theinformation theory into system reliability assessment, using the addible characteristic ofinformation quantity and the principle of equivalence of information quantity, an entropy method ofdata information conversion is presented for the system consisted of identical exponential units.The basic conversion formulae of entropy method of unit test data are derived based on the principleof information quantity equivalence. The general models of entropy method synthesis assessment forsystem reliability approximate lower limits are established according to the fundamental principleof the unit reliability assessment. The applications of the entropy method are discussed by way ofpractical examples. Compared with the traditional methods, the entropy method is found to be validand practicable and the assessment results are very satisfactory.展开更多
Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time (HB-AFT) do- main technique has some advantages in dealing with nonlinear ...Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time (HB-AFT) do- main technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.展开更多
In this article,we consider a new family of exponential type estimators for estimating the unknown population mean of the study variable.We propose estimators taking advantage of the auxiliary variable information und...In this article,we consider a new family of exponential type estimators for estimating the unknown population mean of the study variable.We propose estimators taking advantage of the auxiliary variable information under the first and second non-response cases separately.The required theoretical comparisons are obtained and the numerical studies are conducted.In conclusion,the results show that the proposed family of estimators is the most efficient estimator with respect to the estimators in literature under the obtained conditions for both cases.展开更多
We present the application of differential quadrature(DQ) method for the buckling analysis of nanobeams with exponentially varying stiffness based on four different beam theories of Euler-Bernoulli, Timoshenko, Redd...We present the application of differential quadrature(DQ) method for the buckling analysis of nanobeams with exponentially varying stiffness based on four different beam theories of Euler-Bernoulli, Timoshenko, Reddy, and Levison.The formulation is based on the nonlocal elasticity theory of Eringen. New results are presented for the guided and simply supported guided boundary conditions. Numerical results are obtained to investigate the effects of the nonlocal parameter,length-to-height ratio, boundary condition, and nonuniform parameter on the critical buckling load parameter. It is observed that the critical buckling load decreases with increase in the nonlocal parameter while the critical buckling load parameter increases with increase in the length-to-height ratio.展开更多
A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. ...A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. It spans simultaneously wide range of exponential decay rates with multi scaling and does not suffer from zero crossing. These two conditions are necessary for many physical problems. For comparison, the method is used to solve different problems and compared with analytical and published results. The comparison exhibits the strengths and accuracy of the presented basis set.展开更多
We present the new predictor-corrector methods for systems of nonlinear differential equations, based on themethod of esponential time differencing. We compare the present schemes with the explicit multistep exponenti...We present the new predictor-corrector methods for systems of nonlinear differential equations, based on themethod of esponential time differencing. We compare the present schemes with the explicit multistep exponential time differecing and Adams Bashforth-Moulton method. The numerical results show that the schemes are more accurate and more efficient than Adams predictor-corrector method. The exponential time differencing method has been developed and perfected by the present studies.展开更多
In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation (RLW) is integrated fully by using an exponentia...In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation (RLW) is integrated fully by using an exponential B-spline Galerkin method in space together with Crank-Nicolson method in time. Three numerical examples related to propagation of sin- gle solitary wave, interaction of two solitary waves and wave generation are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.展开更多
By the discussion of the formula and properties of (4,4) parametric form rational approximation to function exp(q), the fourth order derivative one_step exponentially fitted method and the third order derivative hybri...By the discussion of the formula and properties of (4,4) parametric form rational approximation to function exp(q), the fourth order derivative one_step exponentially fitted method and the third order derivative hybrid one_step exponentially fitted method are presented, their order p satisfying 6≤p≤8. The necessary and sufficient conditions for the two methods to be A_ stable are given. Finally, for the fourth order derivative method, the error bound and the necessary and sufficient conditions for it to be median are discussed.展开更多
The Filon-type quadrature is efficient for highly oscillatory functions - Fourier transforms. Based on Cox and Matthews' ETD schemes, the higher order single step exponential time differencing schemes are presente...The Filon-type quadrature is efficient for highly oscillatory functions - Fourier transforms. Based on Cox and Matthews' ETD schemes, the higher order single step exponential time differencing schemes are presented based on the Filon-type integration and the A-stability of the two-order Adams-Bashforth exponential time differencing scheme is considered. The effectiveness and accuracy of the schemes is tested.展开更多
Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Th...Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Though, this paper shows that recent advance methods can be more favored. In this work, we have incorporated, throughout numerical comparison experiments, spectral methods, for the space discretization, in conjunction with second and fourth-order time integrating methods for approximating the solution of the reaction-diffusion differential equations. The results have revealed that these methods have advantages over the conventional methods, some of which to mention are: the ease of implementation, accuracy and CPU time.展开更多
In this paper, the analytical solution of a viscous and incompressible fluid towards an exponentially stretching porous sheet with surface heat flux in porous medium, for the boundary layer and heat transfer flow, is ...In this paper, the analytical solution of a viscous and incompressible fluid towards an exponentially stretching porous sheet with surface heat flux in porous medium, for the boundary layer and heat transfer flow, is presented. The equations of continuity, momentum and the energy are transformed into non-linear ordinary differential by using similarity transformation. The solutions of these highly non-linear ordinary differential equations are found analytically by means of Homotopy Analysis Method (HAM). The result obtained by HAM is compared with numerical results presented in the literature. The accuracy of the HAM is indicated by close agreement of the two sets of results. By this method, an expression is obtained which is admissible for all values of effective parameters. This method has the ability to control the convergence of the solution.展开更多
Following a six-step flow chart, exponentially-fitted variant of the 2-step Simpson’s method suitable for solving ordinary differential equations with periodic/oscillatory behaviour is constructed. The qualitative pr...Following a six-step flow chart, exponentially-fitted variant of the 2-step Simpson’s method suitable for solving ordinary differential equations with periodic/oscillatory behaviour is constructed. The qualitative properties of the constructed methods are also investigated. Numerical experiments on standard problems confirming the theoretical expectations regarding the constructed methods compared with other existing standard methods are also presented. Our results unify and improve the existing classical 2-step Simpson’s method.展开更多
We propose a symplectic partitioned Runge-Kutta (SPRK) method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this te...We propose a symplectic partitioned Runge-Kutta (SPRK) method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete (NAD) operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction (LWC) method and the eighth-order staggered-grid (SG)method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.展开更多
基金supported by ONR UMass Dartmouth Marine and UnderSea Technology(MUST)grant N00014-20-1-2849 under the project S31320000049160by DOE grant DE-SC0023164 sub-award RC114586-UMD+2 种基金by AFOSR grants FA9550-18-1-0383 and FA9550-23-1-0037supported by Michigan State University,by AFOSR grants FA9550-19-1-0281 and FA9550-18-1-0383by DOE grant DE-SC0023164.
文摘Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implicit computations and full precision for the explicit computations.In this work,we analyze the stability properties of these methods and their sensitivity to the low-precision rounding errors,and demonstrate their performance in terms of accuracy and efficiency.We develop codes in FORTRAN and Julia to solve nonlinear systems of ODEs and PDEs using the mixed-precision additive Runge-Kutta(MP-ARK)methods.The convergence,accuracy,and runtime of these methods are explored.We show that for a given level of accuracy,suitably chosen MP-ARK methods may provide significant reductions in runtime.
基金supported by the NSF under Grant DMS-2208391sponsored by the NSF under Grant DMS-1753581.
文摘This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.
基金The work is supported in part by China Fundamental Research of Civil Aircraft under grant number MJ-F-2012-04the Fundamental Research Funds for the Central Universities(YWF-15-SXXY-017).
文摘In this paper we develop explicit fast exponential Runge-Kutta methods for the numerical solutions of a class of parabolic equations.By incorporating the linear splitting technique into the explicit exponential Runge-Kutta schemes,we are able to greatly improve the numerical stability.The proposed numerical methods could be fast implemented through use of decompositions of compact spatial difference operators on a regular mesh together with discrete fast Fourier transform techniques.The exponential Runge-Kutta schemes are easy to be adopted in adaptive temporal approximations with variable time step sizes,as well as applied to stiff nonlinearity and boundary conditions of different types.Linear stabilities of the proposed schemes and their comparison with other schemes are presented.We also numerically demonstrate accuracy,stability and robustness of the proposed method through some typical model problems.
基金supported by the NSF grant DMS-1114546 and NSF Research Network in Mathematical Sciences“KI-Net:Kinetic description of emerging challenges in multiscale problems of natural sciences”X.Y.was partially supported by the startup funding of University of California,Santa Barbara。
文摘This paper generalizes the exponential Runge-Kutta asymptotic preserving(AP)method developed in[G.Dimarco and L.Pareschi,SIAM Numer.Anal.,49(2011),pp.2057–2077]to compute the multi-species Boltzmann equation.Compared to the single species Boltzmann equation that the method was originally applied on,this set of equation presents a new difficulty that comes from the lack of local conservation laws due to the interaction between different species.Hence extra stiff nonlinear source terms need to be treated properly to maintain the accuracy and the AP property.The method we propose does not contain any nonlinear nonlocal implicit solver,and can capture the hydrodynamic limit with time step and mesh size independent of the Knudsen number.We prove the positivity and strong AP properties of the scheme,which are verified by two numerical examples.
文摘In this work,the exponential approximation is used for the numerical simulation of a nonlinear SITR model as a system of differential equations that shows the dynamics of the new coronavirus(COVID-19).The SITR mathematical model is divided into four classes using fractal parameters for COVID-19 dynamics,namely,susceptible(S),infected(I),treatment(T),and recovered(R).The main idea of the presented method is based on the matrix representations of the exponential functions and their derivatives using collocation points.To indicate the usefulness of this method,we employ it in some cases.For error analysis of the method,the residual of the solutions is reviewed.The reported examples show that the method is reasonably efficient and accurate.
文摘The reliability assessment of unit-system near two levels is the mostimportant content in the reliability multi-level synthesis of complex systems. Introducing theinformation theory into system reliability assessment, using the addible characteristic ofinformation quantity and the principle of equivalence of information quantity, an entropy method ofdata information conversion is presented for the system consisted of identical exponential units.The basic conversion formulae of entropy method of unit test data are derived based on the principleof information quantity equivalence. The general models of entropy method synthesis assessment forsystem reliability approximate lower limits are established according to the fundamental principleof the unit reliability assessment. The applications of the entropy method are discussed by way ofpractical examples. Compared with the traditional methods, the entropy method is found to be validand practicable and the assessment results are very satisfactory.
基金supported by the National Natural Science Foundation of China(No.10632040)
文摘Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time (HB-AFT) do- main technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.
文摘In this article,we consider a new family of exponential type estimators for estimating the unknown population mean of the study variable.We propose estimators taking advantage of the auxiliary variable information under the first and second non-response cases separately.The required theoretical comparisons are obtained and the numerical studies are conducted.In conclusion,the results show that the proposed family of estimators is the most efficient estimator with respect to the estimators in literature under the obtained conditions for both cases.
文摘We present the application of differential quadrature(DQ) method for the buckling analysis of nanobeams with exponentially varying stiffness based on four different beam theories of Euler-Bernoulli, Timoshenko, Reddy, and Levison.The formulation is based on the nonlocal elasticity theory of Eringen. New results are presented for the guided and simply supported guided boundary conditions. Numerical results are obtained to investigate the effects of the nonlocal parameter,length-to-height ratio, boundary condition, and nonuniform parameter on the critical buckling load parameter. It is observed that the critical buckling load decreases with increase in the nonlocal parameter while the critical buckling load parameter increases with increase in the length-to-height ratio.
文摘A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. It spans simultaneously wide range of exponential decay rates with multi scaling and does not suffer from zero crossing. These two conditions are necessary for many physical problems. For comparison, the method is used to solve different problems and compared with analytical and published results. The comparison exhibits the strengths and accuracy of the presented basis set.
基金The project supported by National Natural Science Foundation of China under Grant No.19902002
文摘We present the new predictor-corrector methods for systems of nonlinear differential equations, based on themethod of esponential time differencing. We compare the present schemes with the explicit multistep exponential time differecing and Adams Bashforth-Moulton method. The numerical results show that the schemes are more accurate and more efficient than Adams predictor-corrector method. The exponential time differencing method has been developed and perfected by the present studies.
基金supported by the Scientific and Technological Research Council of Turkey(Grant No.113F394)
文摘In this paper, an approximate function for the Galerkin method is composed using the combination of the exponential B-spline functions. Regularized long wave equation (RLW) is integrated fully by using an exponential B-spline Galerkin method in space together with Crank-Nicolson method in time. Three numerical examples related to propagation of sin- gle solitary wave, interaction of two solitary waves and wave generation are employed to illustrate the accuracy and the efficiency of the method. Obtained results are compared with some early studies.
基金the Science Technology Foundation of Ministry of Machine_ Buildin
文摘By the discussion of the formula and properties of (4,4) parametric form rational approximation to function exp(q), the fourth order derivative one_step exponentially fitted method and the third order derivative hybrid one_step exponentially fitted method are presented, their order p satisfying 6≤p≤8. The necessary and sufficient conditions for the two methods to be A_ stable are given. Finally, for the fourth order derivative method, the error bound and the necessary and sufficient conditions for it to be median are discussed.
基金Projects(02JJY2006, 03JJY2001) supported by Natural Science Foundation of Hunan Province project supported by JSPS Fellowship Research Program
文摘The Filon-type quadrature is efficient for highly oscillatory functions - Fourier transforms. Based on Cox and Matthews' ETD schemes, the higher order single step exponential time differencing schemes are presented based on the Filon-type integration and the A-stability of the two-order Adams-Bashforth exponential time differencing scheme is considered. The effectiveness and accuracy of the schemes is tested.
基金Supported by National Natural Science Foundation of China(10571036)the Key Discipline Development Program of Beijing Municipal Commission (XK100080537)
文摘Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Though, this paper shows that recent advance methods can be more favored. In this work, we have incorporated, throughout numerical comparison experiments, spectral methods, for the space discretization, in conjunction with second and fourth-order time integrating methods for approximating the solution of the reaction-diffusion differential equations. The results have revealed that these methods have advantages over the conventional methods, some of which to mention are: the ease of implementation, accuracy and CPU time.
文摘In this paper, the analytical solution of a viscous and incompressible fluid towards an exponentially stretching porous sheet with surface heat flux in porous medium, for the boundary layer and heat transfer flow, is presented. The equations of continuity, momentum and the energy are transformed into non-linear ordinary differential by using similarity transformation. The solutions of these highly non-linear ordinary differential equations are found analytically by means of Homotopy Analysis Method (HAM). The result obtained by HAM is compared with numerical results presented in the literature. The accuracy of the HAM is indicated by close agreement of the two sets of results. By this method, an expression is obtained which is admissible for all values of effective parameters. This method has the ability to control the convergence of the solution.
文摘Following a six-step flow chart, exponentially-fitted variant of the 2-step Simpson’s method suitable for solving ordinary differential equations with periodic/oscillatory behaviour is constructed. The qualitative properties of the constructed methods are also investigated. Numerical experiments on standard problems confirming the theoretical expectations regarding the constructed methods compared with other existing standard methods are also presented. Our results unify and improve the existing classical 2-step Simpson’s method.
基金This research was supported by the National Natural Science Foundation of China (Nos. 41230210 and 41204074), the Science Foundation of the Education Department of Yunnan Province (No. 2013Z152), and Statoil Company (Contract No. 4502502663).
文摘We propose a symplectic partitioned Runge-Kutta (SPRK) method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete (NAD) operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction (LWC) method and the eighth-order staggered-grid (SG)method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.