Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved thro...Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.展开更多
Under harsh conditions （such as high temperature, high pressure, and millisecond lifetime chemical reaction）, a long-standing challenge remains to accurately predict the growth characteristics of nanosize spherical ...Under harsh conditions （such as high temperature, high pressure, and millisecond lifetime chemical reaction）, a long-standing challenge remains to accurately predict the growth characteristics of nanosize spherical particles and to determine the rapid chemical reaction flow field characteristics, The growth characteristics of similar spherical oxide nanoparticles are further studied by successfully introducing the space-time conservation element-solution element （CE/SE） algorithm with the monodisperse Kruis model. This approach overcomes the nanosize particle rapid growth limit set and successfully captures the characteristics of the rapid gaseous chemical reaction process. The results show that this approach quantitatively captures the characteristics of the rapid chemical reaction, nanosize particle growth and size distribution. To reveal the growth mechanism for numerous types of oxide nanoparticles, it is very important to choose a rational numerical method and particle physics model.展开更多
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.展开更多
In this article we propose a higher-order space-time conservative method for hyperbolic systems with stiff and non stiff source terms as well as relaxation systems.We call the scheme a slope propagation(SP)method.It i...In this article we propose a higher-order space-time conservative method for hyperbolic systems with stiff and non stiff source terms as well as relaxation systems.We call the scheme a slope propagation(SP)method.It is an extension of our scheme derived for homogeneous hyperbolic systems[1].In the present inhomogeneous systems the relaxation time may vary from order of one to a very small value.These small values make the relaxation term stronger and highly stiff.In such situations underresolved numerical schemes may produce spurious numerical results.However,our present scheme has the capability to correctly capture the behavior of the physical phenomena with high order accuracy even if the initial layer and the small relaxation time are not numerically resolved.The scheme treats the space and time in a unified manner.The flow variables and their slopes are the basic unknowns in the scheme.The source term is treated by its volumetric integration over the space-time control volume and is a direct part of the overall space-time flux balance.We use two approaches for the slope calculations of the flow variables,the first one results directly from the flux balance over the control volumes,while in the second one we use a finite difference approach.The main features of the scheme are its simplicity,its Jacobian-free and Riemann solver-free recipe,as well as its efficiency and high of order accuracy.In particular we show that the scheme has a discrete analog of the continuous asymptotic limit.We have implemented our scheme for various test models available in the literature such as the Broadwell model,the extended thermodynamics equations,the shallow water equations,traffic flow and the Euler equations with heat transfer.The numerical results validate the accuracy,versatility and robustness of the present scheme.展开更多
The space-time conservation element and solution element(CE/SE)method is proposed for solving a conservative interface-capturing reducedmodel of compressible two-fluid flows.The flow equations are the bulk equations,c...The space-time conservation element and solution element(CE/SE)method is proposed for solving a conservative interface-capturing reducedmodel of compressible two-fluid flows.The flow equations are the bulk equations,combined with mass and energy equations for one of the two fluids.The latter equation contains a source term for accounting the energy exchange.The one and two-dimensional flow models are numerically investigated in this manuscript.The CE/SE method is capable to accurately capture the sharp propagating wavefronts of the fluids without excessive numerical diffusion or spurious oscillations.In contrast to the existing upwind finite volume schemes,the Riemann solver and reconstruction procedure are not the building block of the suggested method.The method differs from the previous techniques because of global and local flux conservation in a space-time domain without resorting to interpolation or extrapolation.In order to reveal the efficiency and performance of the approach,several numerical test cases are presented.For validation,the results of the current method are compared with other finite volume schemes.展开更多
In this paper,we investigate the coupling of the Multi-dimensional Optimal Order Detection(MOOD)method and the Arbitrary high order DERivatives(ADER)approach in order to design a new high order accurate,robust and com...In this paper,we investigate the coupling of the Multi-dimensional Optimal Order Detection(MOOD)method and the Arbitrary high order DERivatives(ADER)approach in order to design a new high order accurate,robust and computationally efficient Finite Volume(FV)scheme dedicated to solve nonlinear systems of hyperbolic conservation laws on unstructured triangular and tetrahedral meshes in two and three space dimensions,respectively.The Multi-dimensional Optimal Order Detection(MOOD)method for 2D and 3D geometries has been introduced in a recent series of papers for mixed unstructured meshes.It is an arbitrary high-order accurate Finite Volume scheme in space,using polynomial reconstructions with a posteriori detection and polynomial degree decrementing processes to deal with shock waves and other discontinuities.In the following work,the time discretization is performed with an elegant and efficient one-step ADER procedure.Doing so,we retain the good properties of the MOOD scheme,that is to say the optimal high-order of accuracy is reached on smooth solutions,while spurious oscillations near singularities are prevented.The ADER technique permits not only to reduce the cost of the overall scheme as shown on a set of numerical tests in 2D and 3D,but it also increases the stability of the overall scheme.A systematic comparison between classical unstructured ADER-WENO schemes and the new ADER-MOOD approach has been carried out for high-order schemes in space and time in terms of cost,robustness,accuracy and efficiency.The main finding of this paper is that the combination of ADER with MOOD generally outperforms the one of ADER and WENO either because at given accuracy MOOD is less expensive(memory and/or CPU time),or because it is more accurate for a given grid resolution.A large suite of classical numerical test problems has been solved on unstructured meshes for three challenging multi-dimensional systems of conservation laws:the Euler equations of compressible gas dynamics,the classical equations of ideal magneto-Hydrodynamics(MHD)and finally the relativistic MHD equations(RMHD),which constitutes a particularly challenging nonlinear system of hyperbolic partial differential equation.All tests are run on genuinely unstructured grids composed of simplex elements.展开更多
This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorpo...This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorporating the contravariant second Piola-Kirchhoff stress tensor, the covariant Green’s strain tensor, and its rates up to order n. This mathematical model permits the study of finite deformation and finite strain compressible deformation physics with an ordered rate dissipation mechanism. Constitutive theories are derived using conjugate pairs in entropy inequality and the representation theorem. The resulting mathematical model is both thermodynamically and mathematically consistent and has closure. The solution of the initial value problems (IVPs) describing evolutions is obtained using a variationally consistent space-time coupled finite element method, derived using space-time residual functional in which the local approximations are in hpk higher-order scalar product spaces. This permits accurate description problem physics over the discretization and also permits precise a posteriori computation of the space-time residual functional, an accurate measure of the accuracy of the computed solution. Model problem studies are presented to demonstrate tensile shock formation, propagation, reflection, and interaction. A unique feature of this research is that tensile shocks can only exist in solid matter, as their existence requires a medium to be elastic (presence of strain), which is only possible in a solid medium. In tensile shock physics, a decrease in the density of the medium caused by tensile waves leads to shock formation ahead of the wave. In contrast, in compressive shocks, an increase in density and the corresponding compressive waves result in the formation of compression shocks behind of the wave. Although these are two similar phenomena, they are inherently different in nature. To our knowledge, this work has not been reported in the published literature.展开更多
The paper addresses tensile shock physics in compressible thermoviscoelastic solids with rheology i.e., in compressible polymeric solids. Polymeric solids have elasticity, dissipation mechanisms and relaxation phenome...The paper addresses tensile shock physics in compressible thermoviscoelastic solids with rheology i.e., in compressible polymeric solids. Polymeric solids have elasticity, dissipation mechanisms and relaxation phenomena due to the presence of long chain molecules in the viscous medium. Thus, the main focus of this investigation is to study how the presence of rheology influences tensile shock physics compared to the tensile shock physics in thermoelastic solids with same elasticity and dissipation but without rheology. A traveling stress wave in polymeric solids leaves nonzero stress signature behind it that naturally influences density. Complete relaxation of the nonzero signatures of stress and associated density change depends upon viscosity of the medium, Deborah number, strength of the stress or velocity wave etc. These aspects of the tensile shock physics in TVES with rheology are investigated in the paper. The mathematical model for finite deformation, finite strain is derived using CBL and CCM, and constitutive theories are derived using conjugate pairs in the entropy inequality and the representation theorem. This mathematical model is thermodynamically and mathematically consistent and has closure. The solution of the IVPs described by this mathematical model related to tensile shock physics in TVES with memory are obtained using space-time coupled finite element method based on space-time residual functional for a space-time strip with time marching. p-version hierarchical space-time local approximations with higher order global differentiability in hpk-scalar product spaces and use of minimally conforming spaces ensure that all space-time integrals over space-time discretization are Riemann. This facilitates more accurate description of the physics in the computational process. Model problem studies are presented to illustrate various aspects of tensile shock physics in compressible TVES with rheology.展开更多
We present a newly developed global magnetohydrodynamic(MHD) model to study the responses of the Earth's magnetosphere to the solar wind. The model is established by using the space-time conservation element and s...We present a newly developed global magnetohydrodynamic(MHD) model to study the responses of the Earth's magnetosphere to the solar wind. The model is established by using the space-time conservation element and solution element(CESE) method in general curvilinear coordinates on a six-component grid system. As a preliminary study, this paper is to present the model's numerical results of the quasi-steady state and the dynamics of the Earth's magnetosphere under steady solar wind flow with due northward interplanetary magnetic field(IMF). The model results are found to be in good agreement with those published by other numerical magnetospheric models.展开更多
基金Project supported by the National Basic Research Program of China (973 program) (No.G1999032804)
文摘Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.
基金This research was financially supported by the National Natural Science Foundation of China （No. 11502282）, the China Scholarship Council Fund （No. 201506425040）, the Natural Science Foundation of Jiangsu Province （No. BK20140178）.
文摘Under harsh conditions （such as high temperature, high pressure, and millisecond lifetime chemical reaction）, a long-standing challenge remains to accurately predict the growth characteristics of nanosize spherical particles and to determine the rapid chemical reaction flow field characteristics, The growth characteristics of similar spherical oxide nanoparticles are further studied by successfully introducing the space-time conservation element-solution element （CE/SE） algorithm with the monodisperse Kruis model. This approach overcomes the nanosize particle rapid growth limit set and successfully captures the characteristics of the rapid gaseous chemical reaction process. The results show that this approach quantitatively captures the characteristics of the rapid chemical reaction, nanosize particle growth and size distribution. To reveal the growth mechanism for numerous types of oxide nanoparticles, it is very important to choose a rational numerical method and particle physics model.
基金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.
文摘In this article we propose a higher-order space-time conservative method for hyperbolic systems with stiff and non stiff source terms as well as relaxation systems.We call the scheme a slope propagation(SP)method.It is an extension of our scheme derived for homogeneous hyperbolic systems[1].In the present inhomogeneous systems the relaxation time may vary from order of one to a very small value.These small values make the relaxation term stronger and highly stiff.In such situations underresolved numerical schemes may produce spurious numerical results.However,our present scheme has the capability to correctly capture the behavior of the physical phenomena with high order accuracy even if the initial layer and the small relaxation time are not numerically resolved.The scheme treats the space and time in a unified manner.The flow variables and their slopes are the basic unknowns in the scheme.The source term is treated by its volumetric integration over the space-time control volume and is a direct part of the overall space-time flux balance.We use two approaches for the slope calculations of the flow variables,the first one results directly from the flux balance over the control volumes,while in the second one we use a finite difference approach.The main features of the scheme are its simplicity,its Jacobian-free and Riemann solver-free recipe,as well as its efficiency and high of order accuracy.In particular we show that the scheme has a discrete analog of the continuous asymptotic limit.We have implemented our scheme for various test models available in the literature such as the Broadwell model,the extended thermodynamics equations,the shallow water equations,traffic flow and the Euler equations with heat transfer.The numerical results validate the accuracy,versatility and robustness of the present scheme.
基金supported by Higher Education Commission(HEC)of Pakistan through grant No.1375.
文摘The space-time conservation element and solution element(CE/SE)method is proposed for solving a conservative interface-capturing reducedmodel of compressible two-fluid flows.The flow equations are the bulk equations,combined with mass and energy equations for one of the two fluids.The latter equation contains a source term for accounting the energy exchange.The one and two-dimensional flow models are numerically investigated in this manuscript.The CE/SE method is capable to accurately capture the sharp propagating wavefronts of the fluids without excessive numerical diffusion or spurious oscillations.In contrast to the existing upwind finite volume schemes,the Riemann solver and reconstruction procedure are not the building block of the suggested method.The method differs from the previous techniques because of global and local flux conservation in a space-time domain without resorting to interpolation or extrapolation.In order to reveal the efficiency and performance of the approach,several numerical test cases are presented.For validation,the results of the current method are compared with other finite volume schemes.
基金the European Research Council(ERC)under the European Union’s Seventh Framework Programme(FP7/2007-2013)the research project STiMulUs,ERC Grant agreement no.278267+1 种基金.R.L.has been partially funded by the ANR under the JCJC project“ALE INC(ubator)3D”the reference LA-UR-13-28795.The authors would like to acknowledge PRACE for awarding access to the SuperMUC supercomputer based in Munich,Germany at the Leibniz Rechenzentrum(LRZ)。
文摘In this paper,we investigate the coupling of the Multi-dimensional Optimal Order Detection(MOOD)method and the Arbitrary high order DERivatives(ADER)approach in order to design a new high order accurate,robust and computationally efficient Finite Volume(FV)scheme dedicated to solve nonlinear systems of hyperbolic conservation laws on unstructured triangular and tetrahedral meshes in two and three space dimensions,respectively.The Multi-dimensional Optimal Order Detection(MOOD)method for 2D and 3D geometries has been introduced in a recent series of papers for mixed unstructured meshes.It is an arbitrary high-order accurate Finite Volume scheme in space,using polynomial reconstructions with a posteriori detection and polynomial degree decrementing processes to deal with shock waves and other discontinuities.In the following work,the time discretization is performed with an elegant and efficient one-step ADER procedure.Doing so,we retain the good properties of the MOOD scheme,that is to say the optimal high-order of accuracy is reached on smooth solutions,while spurious oscillations near singularities are prevented.The ADER technique permits not only to reduce the cost of the overall scheme as shown on a set of numerical tests in 2D and 3D,but it also increases the stability of the overall scheme.A systematic comparison between classical unstructured ADER-WENO schemes and the new ADER-MOOD approach has been carried out for high-order schemes in space and time in terms of cost,robustness,accuracy and efficiency.The main finding of this paper is that the combination of ADER with MOOD generally outperforms the one of ADER and WENO either because at given accuracy MOOD is less expensive(memory and/or CPU time),or because it is more accurate for a given grid resolution.A large suite of classical numerical test problems has been solved on unstructured meshes for three challenging multi-dimensional systems of conservation laws:the Euler equations of compressible gas dynamics,the classical equations of ideal magneto-Hydrodynamics(MHD)and finally the relativistic MHD equations(RMHD),which constitutes a particularly challenging nonlinear system of hyperbolic partial differential equation.All tests are run on genuinely unstructured grids composed of simplex elements.
文摘This paper addresses tensile shock physics in thermoviscoelastic (TVE) solids without memory. The mathematical model is derived using conservation and balance laws (CBL) of classical continuum mechanics (CCM), incorporating the contravariant second Piola-Kirchhoff stress tensor, the covariant Green’s strain tensor, and its rates up to order n. This mathematical model permits the study of finite deformation and finite strain compressible deformation physics with an ordered rate dissipation mechanism. Constitutive theories are derived using conjugate pairs in entropy inequality and the representation theorem. The resulting mathematical model is both thermodynamically and mathematically consistent and has closure. The solution of the initial value problems (IVPs) describing evolutions is obtained using a variationally consistent space-time coupled finite element method, derived using space-time residual functional in which the local approximations are in hpk higher-order scalar product spaces. This permits accurate description problem physics over the discretization and also permits precise a posteriori computation of the space-time residual functional, an accurate measure of the accuracy of the computed solution. Model problem studies are presented to demonstrate tensile shock formation, propagation, reflection, and interaction. A unique feature of this research is that tensile shocks can only exist in solid matter, as their existence requires a medium to be elastic (presence of strain), which is only possible in a solid medium. In tensile shock physics, a decrease in the density of the medium caused by tensile waves leads to shock formation ahead of the wave. In contrast, in compressive shocks, an increase in density and the corresponding compressive waves result in the formation of compression shocks behind of the wave. Although these are two similar phenomena, they are inherently different in nature. To our knowledge, this work has not been reported in the published literature.
文摘The paper addresses tensile shock physics in compressible thermoviscoelastic solids with rheology i.e., in compressible polymeric solids. Polymeric solids have elasticity, dissipation mechanisms and relaxation phenomena due to the presence of long chain molecules in the viscous medium. Thus, the main focus of this investigation is to study how the presence of rheology influences tensile shock physics compared to the tensile shock physics in thermoelastic solids with same elasticity and dissipation but without rheology. A traveling stress wave in polymeric solids leaves nonzero stress signature behind it that naturally influences density. Complete relaxation of the nonzero signatures of stress and associated density change depends upon viscosity of the medium, Deborah number, strength of the stress or velocity wave etc. These aspects of the tensile shock physics in TVES with rheology are investigated in the paper. The mathematical model for finite deformation, finite strain is derived using CBL and CCM, and constitutive theories are derived using conjugate pairs in the entropy inequality and the representation theorem. This mathematical model is thermodynamically and mathematically consistent and has closure. The solution of the IVPs described by this mathematical model related to tensile shock physics in TVES with memory are obtained using space-time coupled finite element method based on space-time residual functional for a space-time strip with time marching. p-version hierarchical space-time local approximations with higher order global differentiability in hpk-scalar product spaces and use of minimally conforming spaces ensure that all space-time integrals over space-time discretization are Riemann. This facilitates more accurate description of the physics in the computational process. Model problem studies are presented to illustrate various aspects of tensile shock physics in compressible TVES with rheology.
基金supported by the National Basic Research Program of China(Grant Nos.2012CB825601,2014CB845903,2012CB825604)the National Natural Science Foundation of China(Grant Nos.41031066,41231068,41274192,41074121,41204127,41174122)+1 种基金the Knowledge Innovation Program of the Chinese Academy of Sciences(Grant No.KZZD-EW-01-4)the Specialized Research Fund for State Key Laboratories
文摘We present a newly developed global magnetohydrodynamic(MHD) model to study the responses of the Earth's magnetosphere to the solar wind. The model is established by using the space-time conservation element and solution element(CESE) method in general curvilinear coordinates on a six-component grid system. As a preliminary study, this paper is to present the model's numerical results of the quasi-steady state and the dynamics of the Earth's magnetosphere under steady solar wind flow with due northward interplanetary magnetic field(IMF). The model results are found to be in good agreement with those published by other numerical magnetospheric models.
