In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
In both quantum and classical field systems,conservation laws such as the conservation of energy and momentum are widely regarded as fundamental properties.A broadly accepted approach to deriving conservation laws is ...In both quantum and classical field systems,conservation laws such as the conservation of energy and momentum are widely regarded as fundamental properties.A broadly accepted approach to deriving conservation laws is built using Noether's method.However,this procedure is still unclear for relativistic particle-field systems where particles are regarded as classical world lines.In the present study,we establish a general manifestly covariant or geometric field theory for classical relativistic particle-field systems.In contrast to quantum systems,where particles are viewed as quantum fields,classical relativistic particle-field systems present specific challenges.These challenges arise from two sides.The first comes from the mass-shell constraint.To deal with the mass-shell constraint,the Euler–Lagrange–Barut(ELB)equation is used to determine the particle's world lines in the four-dimensional(4D)Minkowski space.Besides,the infinitesimal criterion,which is a differential equation in formal field theory,is reconstructed by an integro-differential form.The other difficulty is that fields and particles depend on heterogeneous manifolds.To overcome this challenge,we propose using a weak version of the ELB equation that allows us to connect local conservation laws and continuous symmetries in classical relativistic particle-field systems.By applying a weak ELB equation to classical relativistic particle-field systems,we can systematically derive local conservation laws by examining the underlying symmetries of the system.Our proposed approach provides a new perspective on understanding conservation laws in classical relativistic particle-field systems.展开更多
We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniq...We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniques in the literature rely on the creation of ghost cells outside the boundary and on some form of extrapolation from the inside that,taking into account the boundary conditions,fills the ghost cells with appropriate values,so that a standard reconstruction can be applied also in the boundary cells.In Naumann et al.(Appl.Math.Comput.325:252–270.https://doi.org/10.1016/j.amc.2017.12.041,2018),motivated by the difficulty of choosing appropriate boundary conditions at the internal nodes of a network,a different technique was explored that avoids the use of ghost cells,but instead employs for the boundary cells a different stencil,biased towards the interior of the domain.In this paper,extending that approach,which does not make use of ghost cells,we propose a more accurate reconstruction for the one-dimensional case and a two-dimensional one for Cartesian grids.In several numerical tests,we compare the novel reconstruction with the standard approach using ghost cells.展开更多
With the development of numerical methods the numerical computations require higher and higher accuracy. This paper is devoted to the high-order local absorbing boundary conditions (ABCs) for heat equation. We prove...With the development of numerical methods the numerical computations require higher and higher accuracy. This paper is devoted to the high-order local absorbing boundary conditions (ABCs) for heat equation. We proved that the coupled system yields a stable problem between the obtained high-order local ABCs and the partial differential equation in the computational domain. This method has been used widely in wave propagation models only recently. We extend the spirit of the methodology to parabolic ones, which will become a basis to design the local ABCs for a class of nonlinear PDEs. Some numerical tests show that the new treatment is very efficient and tractable.展开更多
This paper is to investigate the J-selfadjointness of a class of high-order complex coefficients differential operators with transmission conditions.Using the Lagrange bilinear form of J-symmetric differential equatio...This paper is to investigate the J-selfadjointness of a class of high-order complex coefficients differential operators with transmission conditions.Using the Lagrange bilinear form of J-symmetric differential equations,the definition of J-selfadjoint differential operators and the method of matrix representation,we prove that the operator is J-selfadjoint operator,and the eigenvectors and eigen-subspaces corresponding to different eigenvalues are C-orthogonal.展开更多
The paper is concerned with the numerical solution of Schr¨odinger equations on an unbounded spatial domain.High-order absorbing boundary conditions for one-dimensional domain are derived,and the stability of the...The paper is concerned with the numerical solution of Schr¨odinger equations on an unbounded spatial domain.High-order absorbing boundary conditions for one-dimensional domain are derived,and the stability of the reduced initial boundary value problem in the computational interval is proved by energy estimate.Then a second order finite difference scheme is proposed,and the convergence of the scheme is established as well.Finally,numerical examples are reported to confirm our error estimates of the numerical methods.展开更多
In this paper,we study the high-order nonlinear Schrodinger equation with periodic initial conditions via the unified transform method extended by Fokas and Lenells.For the high-order nonlinear Schrodinger equation,th...In this paper,we study the high-order nonlinear Schrodinger equation with periodic initial conditions via the unified transform method extended by Fokas and Lenells.For the high-order nonlinear Schrodinger equation,the initial value problem on the circle can be expressed in terms of the solution of a Riemann–Hilbert problem.The related jump matrix can be explicitly expressed based on the initial data alone.Furthermore,we present the explicit solution,which corresponds to a one-gap solution.展开更多
Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis (VTI media), a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis (TTI me...Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis (VTI media), a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis (TTI media). This is achieved using projection transformation, which rotates the direction vector in the coordinate system of observation toward the direction vector for the coordinate system in which the z-component is parallel to the symmetry axis of the TTI media. The equation has a simple form, is easily calculated, is not influenced by the pseudo-shear wave, and can be calculated reliably when δ is greater than ε. The finite difference method is used to solve the equation. In addition, a perfectly matched layer (PML) absorbing boundary condition is obtained for the equation. Theoretical analysis and numerical simulation results with forward modeling prove that the equation can accurately simulate a quasi-P wave in TTI medium.展开更多
The mixed finite element method is used to solve the exterior Poisson equations with higher-order local artificial boundary conditions in 3-D space. New unknowns are introduced to reduce the order of the derivatives ...The mixed finite element method is used to solve the exterior Poisson equations with higher-order local artificial boundary conditions in 3-D space. New unknowns are introduced to reduce the order of the derivatives of the unknown to two. The result is an equivalent mixed variational problem which was solved using bilinear finite elements. The primary advantage is that special finite elements are not needed on the adjacent layer of the artificial boundary for the higher-order derivatives. Error estimates are obtained for some local artificial boundary conditions with prescibed orders. A numerical example demonstrates the effectiveness of this method.展开更多
According to Godunov theorem for numerical calculations of advection equations, there exist no high-er-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy ex...According to Godunov theorem for numerical calculations of advection equations, there exist no high-er-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy exceeding the first-order. In case of advection-diffusion equations, so far there have been not found stable schemes with positive difference coefficients in a family of numerical schemes exceeding the second-order accuracy. We propose a third-order computational scheme for numerical fluxes to guarantee the non-negative difference coefficients of resulting finite difference equations for advection-diffusion equations. The present scheme is optimized so as to minimize truncation errors for the numerical fluxes while fulfilling the positivity condition of the difference coefficients which are variable depending on the local Courant number and diffusion number. The feature of the present optimized scheme consists in keeping the third-order accuracy anywhere without any numerical flux limiter by using the same stencil number as convemtional third-order shemes such as KAWAMURA and UTOPIA schemes. We extend the present method into multi-dimensional equations. Numerical experiments for linear and nonlinear advection-diffusion equations were performed and the present scheme’s applicability to nonlinear Burger’s equation was confirmed.展开更多
We employ the Riemann-Hilbert(RH)method to study the Hirota equation with arbitrary order zero poles under zero boundary conditions.Through the spectral analysis,the asymptoticity,symmetry,and analysis of the Jost fun...We employ the Riemann-Hilbert(RH)method to study the Hirota equation with arbitrary order zero poles under zero boundary conditions.Through the spectral analysis,the asymptoticity,symmetry,and analysis of the Jost functions are obtained,which play a key role in constructing the RH problem.Then we successfully established the exact solution of the equation without reflection potential by solving the RH problem.Choosing some appropriate parameters of the resulting solutions,we further derive the soliton solutions with different order poles,including four cases of a fourthorder pole,two second-order poles,a third-order pole and a first-order pole,and four first-order points.Finally,the dynamical behavior of these solutions are analyzed via graphic analysis.展开更多
A new weak boundary procedure for hyperbolic problems is presented.We consider high order finite difference operators of summation-by-parts form with weak boundary conditions and generalize that technique.The new boun...A new weak boundary procedure for hyperbolic problems is presented.We consider high order finite difference operators of summation-by-parts form with weak boundary conditions and generalize that technique.The new boundary procedure is applied near boundaries in an extended domain where data is known.We show how to raise the order of accuracy of the scheme,how to modify the spectrum of the resulting operator and how to construct non-reflecting properties at the boundaries.The new boundary procedure is cheap,easy to implement and suitable for all numerical methods,not only finite difference methods,that employ weak boundary conditions.Numerical results that corroborate the analysis are presented.展开更多
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
基金supported by National Natural Science Foundation of China(No.12005141)supported by National Natural Science Foundation of China(No.11805273)+2 种基金supported by the Collaborative Innovation Program of Hefei Science Center,CAS(No.2021HSCCIP019)National MC Energy R&D Program(No.2018YFE0304100)National Natural Science Foundation of China(No.11905220)。
文摘In both quantum and classical field systems,conservation laws such as the conservation of energy and momentum are widely regarded as fundamental properties.A broadly accepted approach to deriving conservation laws is built using Noether's method.However,this procedure is still unclear for relativistic particle-field systems where particles are regarded as classical world lines.In the present study,we establish a general manifestly covariant or geometric field theory for classical relativistic particle-field systems.In contrast to quantum systems,where particles are viewed as quantum fields,classical relativistic particle-field systems present specific challenges.These challenges arise from two sides.The first comes from the mass-shell constraint.To deal with the mass-shell constraint,the Euler–Lagrange–Barut(ELB)equation is used to determine the particle's world lines in the four-dimensional(4D)Minkowski space.Besides,the infinitesimal criterion,which is a differential equation in formal field theory,is reconstructed by an integro-differential form.The other difficulty is that fields and particles depend on heterogeneous manifolds.To overcome this challenge,we propose using a weak version of the ELB equation that allows us to connect local conservation laws and continuous symmetries in classical relativistic particle-field systems.By applying a weak ELB equation to classical relativistic particle-field systems,we can systematically derive local conservation laws by examining the underlying symmetries of the system.Our proposed approach provides a new perspective on understanding conservation laws in classical relativistic particle-field systems.
基金MIUR-PRIN project 2017KKJP4X“Innovative numerical methods for evolutionary partial differential equations and applications”.Gabriella Puppo acknowledges also the support of 2019 Ateneo Sapienza research project no.RM11916B51CD40E1.
文摘We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniques in the literature rely on the creation of ghost cells outside the boundary and on some form of extrapolation from the inside that,taking into account the boundary conditions,fills the ghost cells with appropriate values,so that a standard reconstruction can be applied also in the boundary cells.In Naumann et al.(Appl.Math.Comput.325:252–270.https://doi.org/10.1016/j.amc.2017.12.041,2018),motivated by the difficulty of choosing appropriate boundary conditions at the internal nodes of a network,a different technique was explored that avoids the use of ghost cells,but instead employs for the boundary cells a different stencil,biased towards the interior of the domain.In this paper,extending that approach,which does not make use of ghost cells,we propose a more accurate reconstruction for the one-dimensional case and a two-dimensional one for Cartesian grids.In several numerical tests,we compare the novel reconstruction with the standard approach using ghost cells.
文摘With the development of numerical methods the numerical computations require higher and higher accuracy. This paper is devoted to the high-order local absorbing boundary conditions (ABCs) for heat equation. We proved that the coupled system yields a stable problem between the obtained high-order local ABCs and the partial differential equation in the computational domain. This method has been used widely in wave propagation models only recently. We extend the spirit of the methodology to parabolic ones, which will become a basis to design the local ABCs for a class of nonlinear PDEs. Some numerical tests show that the new treatment is very efficient and tractable.
基金supported by the National Natural Science Foundation of China(Grant No.12261066)。
文摘This paper is to investigate the J-selfadjointness of a class of high-order complex coefficients differential operators with transmission conditions.Using the Lagrange bilinear form of J-symmetric differential equations,the definition of J-selfadjoint differential operators and the method of matrix representation,we prove that the operator is J-selfadjoint operator,and the eigenvectors and eigen-subspaces corresponding to different eigenvalues are C-orthogonal.
基金supported by FRG of Hong Kong Baptist University,RGC of Hong Kong,Natural Science Foundation of China(Grant Number 10871044)Singapore AcRF RG59/08(M52110092)NRF 2007IDM-IDM002-010.
文摘The paper is concerned with the numerical solution of Schr¨odinger equations on an unbounded spatial domain.High-order absorbing boundary conditions for one-dimensional domain are derived,and the stability of the reduced initial boundary value problem in the computational interval is proved by energy estimate.Then a second order finite difference scheme is proposed,and the convergence of the scheme is established as well.Finally,numerical examples are reported to confirm our error estimates of the numerical methods.
基金funded by National Natural Science Foundation of China(No.11471215)。
文摘In this paper,we study the high-order nonlinear Schrodinger equation with periodic initial conditions via the unified transform method extended by Fokas and Lenells.For the high-order nonlinear Schrodinger equation,the initial value problem on the circle can be expressed in terms of the solution of a Riemann–Hilbert problem.The related jump matrix can be explicitly expressed based on the initial data alone.Furthermore,we present the explicit solution,which corresponds to a one-gap solution.
基金supported by the National Natural Science Foundation of China(No.41674118)the national science and technology major project(No.2016ZX05027-002)
文摘Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis (VTI media), a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis (TTI media). This is achieved using projection transformation, which rotates the direction vector in the coordinate system of observation toward the direction vector for the coordinate system in which the z-component is parallel to the symmetry axis of the TTI media. The equation has a simple form, is easily calculated, is not influenced by the pseudo-shear wave, and can be calculated reliably when δ is greater than ε. The finite difference method is used to solve the equation. In addition, a perfectly matched layer (PML) absorbing boundary condition is obtained for the equation. Theoretical analysis and numerical simulation results with forward modeling prove that the equation can accurately simulate a quasi-P wave in TTI medium.
基金Supported by the National Natural Science Foundationof China(No.19772 0 2 2 )
文摘The mixed finite element method is used to solve the exterior Poisson equations with higher-order local artificial boundary conditions in 3-D space. New unknowns are introduced to reduce the order of the derivatives of the unknown to two. The result is an equivalent mixed variational problem which was solved using bilinear finite elements. The primary advantage is that special finite elements are not needed on the adjacent layer of the artificial boundary for the higher-order derivatives. Error estimates are obtained for some local artificial boundary conditions with prescibed orders. A numerical example demonstrates the effectiveness of this method.
文摘According to Godunov theorem for numerical calculations of advection equations, there exist no high-er-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy exceeding the first-order. In case of advection-diffusion equations, so far there have been not found stable schemes with positive difference coefficients in a family of numerical schemes exceeding the second-order accuracy. We propose a third-order computational scheme for numerical fluxes to guarantee the non-negative difference coefficients of resulting finite difference equations for advection-diffusion equations. The present scheme is optimized so as to minimize truncation errors for the numerical fluxes while fulfilling the positivity condition of the difference coefficients which are variable depending on the local Courant number and diffusion number. The feature of the present optimized scheme consists in keeping the third-order accuracy anywhere without any numerical flux limiter by using the same stencil number as convemtional third-order shemes such as KAWAMURA and UTOPIA schemes. We extend the present method into multi-dimensional equations. Numerical experiments for linear and nonlinear advection-diffusion equations were performed and the present scheme’s applicability to nonlinear Burger’s equation was confirmed.
基金supported by the National Natural Science Foundation of China under Grant No.11975306the Natural Science Foundation of Jiangsu Province under Grant No.BK20181351+2 种基金the Six Talent Peaks Project in Jiangsu Province under Grant No.JY-059the Fundamental Research Fund for the Central Universities under the Grant Nos.2019ZDPY07 and 2019QNA35the Postgraduate Research&Practice Innovation Program of Jiangsu Province under Grant No.KYCX212152.
文摘We employ the Riemann-Hilbert(RH)method to study the Hirota equation with arbitrary order zero poles under zero boundary conditions.Through the spectral analysis,the asymptoticity,symmetry,and analysis of the Jost functions are obtained,which play a key role in constructing the RH problem.Then we successfully established the exact solution of the equation without reflection potential by solving the RH problem.Choosing some appropriate parameters of the resulting solutions,we further derive the soliton solutions with different order poles,including four cases of a fourthorder pole,two second-order poles,a third-order pole and a first-order pole,and four first-order points.Finally,the dynamical behavior of these solutions are analyzed via graphic analysis.
基金supported by the National Science Foundation under Award No.0948304 and by the Southern California Earthquake Center.SCEC is funded by NSF Cooperative Agreement EAR-0529922 and USGS Cooperative Agreement 07HQAG0008(SCEC contribution number 1806).The work by the last author was carried out within the Swedish e-science Research Centre(SeRC)and supported by the Swedish Research Council(VR).
文摘A new weak boundary procedure for hyperbolic problems is presented.We consider high order finite difference operators of summation-by-parts form with weak boundary conditions and generalize that technique.The new boundary procedure is applied near boundaries in an extended domain where data is known.We show how to raise the order of accuracy of the scheme,how to modify the spectrum of the resulting operator and how to construct non-reflecting properties at the boundaries.The new boundary procedure is cheap,easy to implement and suitable for all numerical methods,not only finite difference methods,that employ weak boundary conditions.Numerical results that corroborate the analysis are presented.