SPACE-TIME FINITE ELEMENT METHOD FOR SCHRDINGER EQUATION AND ITS CONSERVATION

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.

Numerical simulation of oxide nanoparticle growth characteristics under the gas detonation chemical reaction by space-time conservation element-solution element method

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.

A Compact High Order Space-Time Method for Conservation Laws

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.

A Space-Time Conservative Method for Hyperbolic Systems with Stiff and Non Stiff Source Terms

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 CE/SE Method for Solving Reduced Two-Fluid Flow Model

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.

A New Family of High Order Unstructured MOOD and ADER Finite Volume Schemes for Multidimensional Systems of Hyperbolic Conservation Laws

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.