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.展开更多
Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms.In our earlier work[31–33],we designed high order well-balanced schemes to a cl...Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms.In our earlier work[31–33],we designed high order well-balanced schemes to a class of hyperbolic systems with separable source terms.In this paper,we present a different approach to the same purpose:designing high order well-balanced finite volume weighted essentially non-oscillatory(WENO)schemes and RungeKutta discontinuous Galerkin(RKDG)finite element methods.We make the observation that the traditional RKDG methods are capable of maintaining certain steady states exactly,if a small modification on either the initial condition or the flux is provided.The computational cost to obtain such a well balanced RKDG method is basically the same as the traditional RKDG method.The same idea can be applied to the finite volume WENO schemes.We will first describe the algorithms and prove the well balanced property for the shallow water equations,and then show that the result can be generalized to a class of other balance laws.We perform extensive one and two dimensional simulations to verify the properties of these schemes such as the exact preservation of the balance laws for certain steady state solutions,the non-oscillatory property for general solutions with discontinuities,and the genuine high order accuracy in smooth regions.展开更多
The lack of descriptions regarding the order of precedence between the local laws of cities with subordinate districts and the regulations of provincial governments in Legislation Law of the People's Republic of C...The lack of descriptions regarding the order of precedence between the local laws of cities with subordinate districts and the regulations of provincial governments in Legislation Law of the People's Republic of China(Legislation Law) has led to two divergent views. One holds that "the local laws of cities with subordinate districts should take precedence over the regulations of provincial governments," while the other supports the exact opposite. This is a value judgment issue in legislation. To reach a solution, we need to clarify the premises based on the characteristics of the laws in question so that a basic common ground can be established for discussion. The first premise for traditional legislation is that a law should be based on experience as well as logic; the second is that the experience of authority subjects, plus the three aspects of logic should outweigh the experience of social subjects, plus the three aspects of logic. With respect to postmodern legislation, the first premise is that experience should override logic, and the second is that the experience of the authority subject should take precedence over that of social subject, with no requirements for logical consistency. Since Legislation Law fal s into the category of postmodern legislation, according to the premises, the argument that the local laws of cities with subordinate districts should take precedence enjoys wider acceptance, but the view is logically challenged in terms of conceptual consistency, system consistency and principle consistency. More studies must be conducted to facilitate the discussion.展开更多
In this paper we present a 2D/3D high order accurate finite volume scheme in the context of direct Arbitrary-Lagrangian-Eulerian algorithms for general hyperbolic systems of partial differential equations with non-con...In this paper we present a 2D/3D high order accurate finite volume scheme in the context of direct Arbitrary-Lagrangian-Eulerian algorithms for general hyperbolic systems of partial differential equations with non-conservative products and stiff source terms.This scheme is constructed with a single stencil polynomial reconstruction operator,a one-step space-time ADER integration which is suitably designed for dealing even with stiff sources,a nodal solver with relaxation to determine the mesh motion,a path-conservative integration technique for the treatment of non-conservative products and an a posteriori stabilization procedure derived from the so-called Multidimensional Optimal Order Detection(MOOD)paradigm.In this work we consider the seven equation Baer-Nunziato model of compressible multi-phase flows as a representative model involving non-conservative products as well as relaxation source terms which are allowed to become stiff.The new scheme is validated against a set of test cases on 2D/3D unstructured moving meshes on parallel machines and the high order of accuracy achieved by the method is demonstrated by performing a numerical convergence study.Classical Riemann problems and explosion problems with exact solutions are simulated in 2D and 3D.The overall numerical code is also profiled to provide an estimate of the computational cost required by each component of the whole algorithm.展开更多
In this paper, an infinite sequence of conservation laws for a generalized variable-coefficient fifth-order Korteweg-de Vries equation in fluids are constructed based on the Backlund transformation. Hirota bilinear fo...In this paper, an infinite sequence of conservation laws for a generalized variable-coefficient fifth-order Korteweg-de Vries equation in fluids are constructed based on the Backlund transformation. Hirota bilinear form and symbolic computation are applied to obtain three kinds of solutions. Variable coefficients can affect the conserved density, associated flux, and appearance of the characteristic lines. Effects of the wave number on the soliton structures are also discussed and types of soliton structures, e.g., the double-periodic soliton, parallel soliton and soliton complexes, are presented.展开更多
Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas.However,exact local energy–momentum conservation laws for the electromagnetic gyrokinetic system...Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas.However,exact local energy–momentum conservation laws for the electromagnetic gyrokinetic system have not been found despite continuous effort.Without such local conservation laws,energy and momentum can be instantaneously transported across spacetime,which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory.The standard Noether procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic field reside on different manifolds.To overcome this difficulty,we develop a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact,local conservation laws,in particular the energy–momentum conservation laws,for the electromagnetic gyrokinetic system.A weak Euler–Lagrange(EL)equation is established to replace the standard EL equation for the particles.It is discovered that an induced weak EL current enters the local conservation laws,and it is the new physics captured by the high-order field theory on heterogeneous manifolds.A recently developed gauge-symmetrization method for high-order electromagnetic field theories using the electromagnetic displacement-potential tensor is applied to render the derived energy–momentum conservation laws electromagnetic gauge-invariant.展开更多
In this paper,we consider high order multi-domain penalty spectral Galerkin methods for the approximation of hyperbolic conservation laws.This formulation has a penalty parameter which can vary in space and time,allow...In this paper,we consider high order multi-domain penalty spectral Galerkin methods for the approximation of hyperbolic conservation laws.This formulation has a penalty parameter which can vary in space and time,allowing for flexibility in the penalty formulation.This flexibility is particularly advantageous for problems with an inhomogeneous mesh.We show that the discontinuous Galerkin method is equivalent to the multi-domain spectral penalty Galerkin method with a particular value of the penalty parameter.The penalty parameter has an effect on both the accuracy and stability of the method.We examine the numerical issues which arise in the implementation of high order multi-domain penalty spectral Galerkin methods.The coefficient truncation method is proposed to prevent the rapid error growth due to round-off errors when high order polynomials are used.Finally,we show that an inconsistent evaluation of the integrals in the penalty method may lead to growth of errors.Numerical examples for linear and nonlinear problems are presented.展开更多
In this paper,the modulation instability(MI),rogue waves(RWs)and conseryation laws of the coupled higher-order nonlinear Schrodinger equation are investigated.According to MI and the 2×2 Lax pair,Darboux-dressing...In this paper,the modulation instability(MI),rogue waves(RWs)and conseryation laws of the coupled higher-order nonlinear Schrodinger equation are investigated.According to MI and the 2×2 Lax pair,Darboux-dressing transformation with an asymptotic expansion method,the existence and properties of the one-,second-,and third-order RWs for the higher-order nonlinear Schrodinger equation are constructed.In addition,the main characteristics of these solutions are discussed through some graphics,which are draw widespread attention in a variety of complex systems such as optics,Bose-Einstein condensates,capillary fow,superfluidity,fluid dynamics,and finance.In addition,infinitely-many conservation laws are established.展开更多
This paper presents a novel high-order space-time method for hyperbolic conservation laws.Two important concepts,the staggered space-time mesh of the space-time conservation element/solution element(CE/SE)method and t...This paper presents a novel high-order space-time method for hyperbolic conservation laws.Two important concepts,the staggered space-time mesh of the space-time conservation element/solution element(CE/SE)method and the local discontinuous basis functions of the space-time discontinuous Galerkin(DG)finite element method,are the two key ingredients of the new scheme.The staggered spacetime mesh is constructed using the cell-vertex structure of the underlying spatial mesh.The universal definitions of CEs and SEs are independent of the underlying spatial mesh and thus suitable for arbitrarily unstructured meshes.The solution within each physical time step is updated alternately at the cell level and the vertex level.For this solution updating strategy and the DG ingredient,the new scheme here is termed as the discontinuous Galerkin cell-vertex scheme(DG-CVS).The high order of accuracy is achieved by employing high-order Taylor polynomials as the basis functions inside each SE.The present DG-CVS exhibits many advantageous features such as Riemann-solver-free,high-order accuracy,point-implicitness,compactness,and ease of handling boundary conditions.Several numerical tests including the scalar advection equations and compressible Euler equations will demonstrate the performance of the new method.展开更多
基金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.
基金supported by ARO grant W911NF-04-1-0291,NSF grant DMS-0510345 and AFOSR grant FA9550-05-1-0123.
文摘Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms.In our earlier work[31–33],we designed high order well-balanced schemes to a class of hyperbolic systems with separable source terms.In this paper,we present a different approach to the same purpose:designing high order well-balanced finite volume weighted essentially non-oscillatory(WENO)schemes and RungeKutta discontinuous Galerkin(RKDG)finite element methods.We make the observation that the traditional RKDG methods are capable of maintaining certain steady states exactly,if a small modification on either the initial condition or the flux is provided.The computational cost to obtain such a well balanced RKDG method is basically the same as the traditional RKDG method.The same idea can be applied to the finite volume WENO schemes.We will first describe the algorithms and prove the well balanced property for the shallow water equations,and then show that the result can be generalized to a class of other balance laws.We perform extensive one and two dimensional simulations to verify the properties of these schemes such as the exact preservation of the balance laws for certain steady state solutions,the non-oscillatory property for general solutions with discontinuities,and the genuine high order accuracy in smooth regions.
基金part of the results(presented in stages)of"Research on the Legislative System of Cities with Subordinate Districts"(16XFX004)-a program of National Social Sciences Fund in Western China"Empirical Research on Local Legislation"(16XW16)-a research focus of Sichuan Academy of Social Sciences under a key program launched by the Publicity Department of the CPC Sichuan Provincial Committee
文摘The lack of descriptions regarding the order of precedence between the local laws of cities with subordinate districts and the regulations of provincial governments in Legislation Law of the People's Republic of China(Legislation Law) has led to two divergent views. One holds that "the local laws of cities with subordinate districts should take precedence over the regulations of provincial governments," while the other supports the exact opposite. This is a value judgment issue in legislation. To reach a solution, we need to clarify the premises based on the characteristics of the laws in question so that a basic common ground can be established for discussion. The first premise for traditional legislation is that a law should be based on experience as well as logic; the second is that the experience of authority subjects, plus the three aspects of logic should outweigh the experience of social subjects, plus the three aspects of logic. With respect to postmodern legislation, the first premise is that experience should override logic, and the second is that the experience of the authority subject should take precedence over that of social subject, with no requirements for logical consistency. Since Legislation Law fal s into the category of postmodern legislation, according to the premises, the argument that the local laws of cities with subordinate districts should take precedence enjoys wider acceptance, but the view is logically challenged in terms of conceptual consistency, system consistency and principle consistency. More studies must be conducted to facilitate the discussion.
基金W.B.has been financed by the European Research Council(ERC)under the European Union’s Seventh Framework Programme(FP7/2007-2013)with the research project STiMulUs,ERC Grant agreement no.278267R.L.has been partially funded by the ANR under the JCJC project“ALE INC(ubator)3D”JS01-012-01the“International Centre for Mathematics and Computer Science in Toulouse”(CIMI)partially supported by ANR-11-LABX-0040-CIMI within the program ANR-11-IDEX-0002-02.The authors would like to acknowledge PRACE for awarding access to the SuperMUC supercomputer based in Munich,Germany at the Leibniz Rechenzentrum(LRZ).Parts of thematerial contained in this work have been elaborated,gathered and tested while W.B.visited the Mathematical Institute of Toulouse for three months and R.L.visited the Dipartimento di Ingegneria Civile Ambientale e Meccanica in Trento for three months.
文摘In this paper we present a 2D/3D high order accurate finite volume scheme in the context of direct Arbitrary-Lagrangian-Eulerian algorithms for general hyperbolic systems of partial differential equations with non-conservative products and stiff source terms.This scheme is constructed with a single stencil polynomial reconstruction operator,a one-step space-time ADER integration which is suitably designed for dealing even with stiff sources,a nodal solver with relaxation to determine the mesh motion,a path-conservative integration technique for the treatment of non-conservative products and an a posteriori stabilization procedure derived from the so-called Multidimensional Optimal Order Detection(MOOD)paradigm.In this work we consider the seven equation Baer-Nunziato model of compressible multi-phase flows as a representative model involving non-conservative products as well as relaxation source terms which are allowed to become stiff.The new scheme is validated against a set of test cases on 2D/3D unstructured moving meshes on parallel machines and the high order of accuracy achieved by the method is demonstrated by performing a numerical convergence study.Classical Riemann problems and explosion problems with exact solutions are simulated in 2D and 3D.The overall numerical code is also profiled to provide an estimate of the computational cost required by each component of the whole algorithm.
基金Supported by the National Natural Science Foundation of China under Grant No.60772023by the Slpported Project under Grant No.SKLSDE-2010ZX-07 of the State Key Laboratory of Software Development Environment,Beijing University of Aeronautics and As tronautics+2 种基金by the Specialized Research Fund for the Doctoral Program of Higher Educatioi under Grant No.200800130006Chinese Ministry of Education,and by the Innovation Foundation for Ph.D.Graduates under Grant Nos.30-0350 and 30-0366Beijing University of Aeronautics and Astronautics
文摘In this paper, an infinite sequence of conservation laws for a generalized variable-coefficient fifth-order Korteweg-de Vries equation in fluids are constructed based on the Backlund transformation. Hirota bilinear form and symbolic computation are applied to obtain three kinds of solutions. Variable coefficients can affect the conserved density, associated flux, and appearance of the characteristic lines. Effects of the wave number on the soliton structures are also discussed and types of soliton structures, e.g., the double-periodic soliton, parallel soliton and soliton complexes, are presented.
基金supported by the Chinese Scholarship Council(CSC)(No.201806340074)Shenzhen Clean Energy Research Institute and National Natural Science Foundation of China(No.12005141)+3 种基金supported by the US Department of Energy(No.DE-AC02-09CH11466)supported by the National MC Energy R&D Program(No.2018YFE0304100)National Key Research and Development Program(Nos.2016YFA0400600,2016YFA0400601 and 2016YFA0400602)the National Natural Science Foundation of China(Nos.11905220 and 11805273)。
文摘Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas.However,exact local energy–momentum conservation laws for the electromagnetic gyrokinetic system have not been found despite continuous effort.Without such local conservation laws,energy and momentum can be instantaneously transported across spacetime,which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory.The standard Noether procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic field reside on different manifolds.To overcome this difficulty,we develop a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact,local conservation laws,in particular the energy–momentum conservation laws,for the electromagnetic gyrokinetic system.A weak Euler–Lagrange(EL)equation is established to replace the standard EL equation for the particles.It is discovered that an induced weak EL current enters the local conservation laws,and it is the new physics captured by the high-order field theory on heterogeneous manifolds.A recently developed gauge-symmetrization method for high-order electromagnetic field theories using the electromagnetic displacement-potential tensor is applied to render the derived energy–momentum conservation laws electromagnetic gauge-invariant.
基金The work of both authors has been supported by the NSF under Grant No.DMS-0608844.
文摘In this paper,we consider high order multi-domain penalty spectral Galerkin methods for the approximation of hyperbolic conservation laws.This formulation has a penalty parameter which can vary in space and time,allowing for flexibility in the penalty formulation.This flexibility is particularly advantageous for problems with an inhomogeneous mesh.We show that the discontinuous Galerkin method is equivalent to the multi-domain spectral penalty Galerkin method with a particular value of the penalty parameter.The penalty parameter has an effect on both the accuracy and stability of the method.We examine the numerical issues which arise in the implementation of high order multi-domain penalty spectral Galerkin methods.The coefficient truncation method is proposed to prevent the rapid error growth due to round-off errors when high order polynomials are used.Finally,we show that an inconsistent evaluation of the integrals in the penalty method may lead to growth of errors.Numerical examples for linear and nonlinear problems are presented.
基金This work is supported by the National Natural Science Foundation of China(Grant No.71690242,No.11731014,No.12001241)the Basic Research Program of Jiangsu Province(Grant No.BK20200885).
文摘In this paper,the modulation instability(MI),rogue waves(RWs)and conseryation laws of the coupled higher-order nonlinear Schrodinger equation are investigated.According to MI and the 2×2 Lax pair,Darboux-dressing transformation with an asymptotic expansion method,the existence and properties of the one-,second-,and third-order RWs for the higher-order nonlinear Schrodinger equation are constructed.In addition,the main characteristics of these solutions are discussed through some graphics,which are draw widespread attention in a variety of complex systems such as optics,Bose-Einstein condensates,capillary fow,superfluidity,fluid dynamics,and finance.In addition,infinitely-many conservation laws are established.
基金This work is supported by the U.S.Air Force Office of Scientific Research(AFOSR)Computational Mathematics Program under the Award No.FA9550-08-1-0122.
文摘This paper presents a novel high-order space-time method for hyperbolic conservation laws.Two important concepts,the staggered space-time mesh of the space-time conservation element/solution element(CE/SE)method and the local discontinuous basis functions of the space-time discontinuous Galerkin(DG)finite element method,are the two key ingredients of the new scheme.The staggered spacetime mesh is constructed using the cell-vertex structure of the underlying spatial mesh.The universal definitions of CEs and SEs are independent of the underlying spatial mesh and thus suitable for arbitrarily unstructured meshes.The solution within each physical time step is updated alternately at the cell level and the vertex level.For this solution updating strategy and the DG ingredient,the new scheme here is termed as the discontinuous Galerkin cell-vertex scheme(DG-CVS).The high order of accuracy is achieved by employing high-order Taylor polynomials as the basis functions inside each SE.The present DG-CVS exhibits many advantageous features such as Riemann-solver-free,high-order accuracy,point-implicitness,compactness,and ease of handling boundary conditions.Several numerical tests including the scalar advection equations and compressible Euler equations will demonstrate the performance of the new method.
