A note on the importance of mass conservation in long-time stability of Navier–Stokes simulations using finite elements

We prove a long-time stability result for the finite element in space, linear extrapolated Crank–Nicolson in time discretization of the Navier–Stokes equations (NSE). From this result and a numerical experiment, we show the importance of discrete mass conservation in long-time simulations of the NSE. That is, we show that using elements that strongly enforce mass conservation can provide significantly more accurate solutions over long times, compared to those that enforce it weakly.     

Conservation laws of turbulence models

The conservation of mass, momentum, energy, helicity, and enstrophy in fluid flow are important because these quantities organize a flow, and characterize change in the flow's structure over time. In turbulent flow, conservation laws remain important in the inertial range of wave numbers, where viscous effects are negligible. It is in the inertial range where energy, helicity (3d), and enstrophy (2d) must be accurately cascaded for a turbulence model to be qualitatively correct. A first and necessary step for an accurate cascade is conservation; however, many turbulent flow simulations are based on turbulence models whose conservation properties are little explored and might be very different from those of the Navier-Stokes equations.We explore conservation laws and approximate conservation laws satisfied by LES turbulence models. For the Leray, Leray deconvolution, Bardina, and Nth order deconvolution models, we give exact or approximate laws for a model mass, momentum, energy, enstrophy and helicity. The possibility of cascades for model quantities is also discussed.     

Efficient unconditionally stable numerical schemes for a modified phase field crystal model with a strong nonlinear vacancy potential

We consider numerical approximations for a modified phase field crystal model with a strong nonlinear vacancy potential. Based on the invariant energy quadratization approach and stabilized strategies, we develop linear, unconditionally energy stable numerical schemes using the first-order Euler method, the second-order backward differentiation formulas and the second-order Crank–Nicolson method, respectively. We rigorously prove the unconditional energy stability, the mass conservation of these three numerical schemes and carry out error estimates in time for the first-order numerical scheme. Various numerical experiments in 2D and 3D are carried out to validate the accuracy, energy stability, mass conservation, and efficiency of the proposed schemes.     

The heat and moisture transfer balance theory of garment simulation

To establish the human body model to analyze the heat and moisture transfer on body surface, a new explicit definition of rational L-recursion surface is given and the L-recursion surfaces, in Grassmann spaces, are constructed by using blossom method of the homogeneous normal pyramid form. Based on our human body model, the balance theory of garment simulation, the heat and moisture transfer balance equations, called ICAD-balance equations are obtained. The balance theory of garment simulation integrally studies the complex system of human body-fabric-environment. At the same time, the method of obtaining the heat and moisture transfer balance equations is also based on the mass conservation law, the energy conservation law and the Fish law of capillarity. A finite volume method is employed to solve the ICAD-balance equations.     

The heat and moisture transfer balance theory of garment simulation

To establish the human body model to analyze the heat and moisture transfer on body surface, a new explicit definition of rational L-recursion surface is given and the L-recursion surfaces, in Grassmann spaces, are constructed by using blossom method of the homogeneous normal pyramid form. Based on our human body model, the balance theory of garment simulation, the heat and moisture transfer balance equations, called ICAD-balance equations are obtained. The balance theory of garment simulation integrally studies the complex system of human body–fabric–environment. At the same time, the method of obtaining the heat and moisture transfer balance equations is also based on the mass conservation law, the energy conservation law and the Fish law of capillarity. A finite volume method is employed to solve the ICAD-balance equations.     

On the blowing up of solutions to quantum hydrodynamic models on bounded domains

The blow-up in finite time for the solutions to the initial-boundary value problem associated to the multi-dimensional quantum hydrodynamic model in a bounded domain is proved. The model consists on conservation of mass equation and a momentum balance equation equivalent to a compressible Euler equations corrected by a dispersion term of the third order in the momentum balance. The proof is based on a priori estimates for the energy functional for a new observable constructed with an auxiliary function, and it is shown that, under suitable boundary conditions and assumptions on the initial data, the solution blows up after a finite time. I.M. Gamba is supported by NSF-DMS0507038. M.P. Gualdani acknowledges partial support from the Deutsche Forschungsgemeinschaft, grants JU359/5 and was partially supported under the Feodor Lynen Research fellowship. P. Zhang is partially supported by the NSF of China under Grant 10525101 and 10421101, and the innovation grant from the Chinese Academy of Sciences. Part of the work was done when P. Zhang visited the Department of Mathematics of Texas University at Austin, the author would like to thank the hospitality of the department. Support from the Institute for Computational Engineering and Sciences at the University of Texas at Austin is also gratefully acknowledged.     

A conservative local discontinuous Galerkin method for the solution of nonlinear Schrdinger equation in two dimensions

In this study, we present a conservative local discontinuous Galerkin(LDG) method for numerically solving the two-dimensional nonlinear Schrdinger(NLS) equation. The NLS equation is rewritten as a firstorder system and then we construct the LDG formulation with appropriate numerical flux. The mass and energy conserving laws for the semi-discrete formulation can be proved based on different choices of numerical fluxes such as the central, alternative and upwind-based flux. We will propose two kinds of time discretization methods for the semi-discrete formulation. One is based on Crank-Nicolson method and can be proved to preserve the discrete mass and energy conservation. The other one is Krylov implicit integration factor(IIF) method which demands much less computational effort. Various numerical experiments are presented to demonstrate the conservation law of mass and energy, the optimal rates of convergence, and the blow-up phenomenon.     

Design of a new dynamical core for global atmospheric models based on some efficient numerical methods

A careful study on the integral properties of the primitive hydrostatic balance equations for baroclinic atmosphere is carried out, and a new scheme to design the global adiabatic model of atmospheric dynamics is presented. This scheme includes a method of weighted equal-area mesh and a fully discrete finite difference method with quadratic and linear conservations for solving the primitive equation system. Using this scheme, we established a new dynamical core with adjustable high resolution acceptable to the available computer capability, which can be very stable without any filtering and smoothing. Especially, some important integral properties are kept unchanged, such as the anti-symmetries of the horizontal advection operators and the vertical convection operator, the mass conservation, the effective energy conservation under the standard stratification approximation, and so on. Some numerical tests on the new dynamical core, respectively regarding its global conservations and its integrated performances in climatic modeling, incorporated with the physical packages from the Community Atmospheric Model Version 2 (CAM2) of National Center for Atmospheric Research (NCAR), are included.     

Design of a new dynamical core for global atmospheric models based on some efficient numerical methods

A careful study on the integral properties of the primitive hydrostatic balance equations for baroclinic atmosphere is carried out, and a new scheme to design the global adiabatic model of atmospheric dynamics is presented. This scheme includes a method of weighted equal-area mesh and a fully discrete finite difference method with quadratic and linear conservations for solving the primitive equation system. Using this scheme, we established a new dynamical core with adjustable high resolution acceptable to the available     

On the hydrodynamic limit of a Ginzburg-Landau lattice model

Summary A lattice system of interacting diffusion processes is investigated. The evolution is attractive and time reversible, the spin satisfies a conservation law. It is shown that the rescaled spin field converges in probability to the corresponding solution to a nonlinear diffusion equation.Supported in part by the Hungarian National Foundation for Scientific Research, grant No. 819/1, and by the Mathematical Department of Rutgers University, N.S.F. grant DMR 8612369     

Mass and momentum conservation of the least-squares spectral collocation method for the Navier-Stokes equations

From the literature, it is known that the Least-Squares Spectral Element Method (LSSEM) for the stationary Stokes equations performs poorly with respect to mass conservation but compensates this lack by a superior conservation of momentum. Furthermore, it is known that the Least-Squares Spectral Collocation Method (LSSCM) leads to superior conservation of mass and momentum for the stationary Stokes equations. In the present paper, we consider mass and momentum conservation of the LSSCM for time-dependent Stokes and Navier-Stokes equations. We observe that the LSSCM leads to improved conservation of mass (and momentum) for these problems. Furthermore, the LSSCM leads to the well-known time-dependent profiles for the velocity and the pressure profiles. To obtain these results, we use only a few elements, each with high polynomial degree, avoid normal equations for solving the overdetermined linear systems of equations and introduce the Clenshaw-Curtis quadrature rule for imposing the average pressure to be zero. Furthermore, we combined the transformation of Gordon and Hall (transfinite mapping) with the least-squares spectral collocation scheme to discretize the internal flow problems.     

Mass and momentum conservation of the least-squares spectral collocation method for the Navier–Stokes equations

From the literature, it is known that the Least-Squares Spectral Element Method (LSSEM) for the stationary Stokes equations performs poorly with respect to mass conservation but compensates this lack by a superior conservation of momentum. Furthermore, it is known that the Least-Squares Spectral Collocation Method (LSSCM) leads to superior conservation of mass and momentum for the stationary Stokes equations. In the present paper, we consider mass and momentum conservation of the LSSCM for time-dependent Stokes and Navier–Stokes equations. We observe that the LSSCM leads to improved conservation of mass (and momentum) for these problems. Furthermore, the LSSCM leads to the well-known time-dependent profiles for the velocity and the pressure profiles. To obtain these results, we use only a few elements, each with high polynomial degree, avoid normal equations for solving the overdetermined linear systems of equations and introduce the Clenshaw–Curtis quadrature rule for imposing the average pressure to be zero. Furthermore, we combined the transformation of Gordon and Hall (transfinite mapping) with the least-squares spectral collocation scheme to discretize the internal flow problems.     

Discrete Schrödinger-Poisson systems preserving energy and mass

A discrete predictor-corrector Schrödinger-Poisson system is proposed which has the property of mass and energy conservation exactly on the discrete level. The discretization is based on the Cranck-Nicholson scheme, which preserves these invariants of the Schrödinger-Poisson system, but involves the solution of nonlinear equations at each time step. A modified linearized scheme is proposed where conservation is achieved by introducing a phase modulation in the corrector step.     

Finite Difference Scheme for Barotropic Gas Equations

An implicit finite difference scheme approximating the equations of barotropic gas flow is proposed. This scheme ensures the positivity of density and the validity of an energy inequality and the mass conservation law. The continuity equation is approximated implicitly. It is proved that the resulting system of nonlinear equations has a solution for any time and space stepsizes. An iterative method for solving the system of nonlinear equations at each time step is proposed.     

Shock waves for compressible navier-stokes equations are stable

It is shown that shock waves for the compressible Navier-Stokes equations are nonlinearly stable. A perturbation of a shock wave tends to the shock wave, properly translated in phase, as time tends to infinity. Through the consideration of conservation of mass, momentum and energy we obtain an a priori estimate of the amount of translation of the shock wave and the strength of the linear and nonlinear diffusion waves that arise due to the perturbation. Our techniques include the energy method for parabolic-hyperbolic systems, the decomposition of waves, and the energy-characteristic method for viscous conservation laws introduced earlier by the author.     

An Eulerian‐Lagrangian control volume scheme for two‐dimensional unsteady advection‐diffusion problems

We develop a mass conservative Eulerian‐Lagrangian control volume scheme (ELCVS) for the solution of the transient advection‐diffusion equations in two space dimensions. This method uses finite volume test functions over the space‐time domain defined by the characteristics within the framework of the class of Eulerian‐Lagrangian localized adjoint characteristic methods (ELLAM). It, therefore, maintains the advantages of characteristic methods in general, and of this class in particular, which include global mass conservation as well as a natural treatment of all types of boundary conditions. However, it differs from other methods in that class in the treatment of the mass storage integrals at the previous time step defined on deformed Lagrangian regions. This treatment is especially attractive for orthogonal rectangular Eulerian grids composed of block elements. In the algorithm, each deformed region is approximated by an eight‐node region with sides drawn parallel to the Eulerian grid, which significantly simplifies the integration compared with the approach used in finite volume ELLAM methods, based on backward tracking, while retaining local mass conservation at no additional expenses in terms of accuracy or CPU consumption. This is verified by numerical tests which show that ELCVS performs as well as standard finite volume ELLAM methods, which have previously been shown to outperform many other well‐received classes of numerical methods for the equations considered. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

The modified method of characteristics with adjusted advection

Summary. The MMOC procedure for approximating the solutions of transport-dominated diffusion problems does not automatically preserve integral conservation laws, leading to (mass) balance errors in many kinds of flow problems. The variant, called the MMOCAA, discussed herein preserves the conservation law at a minor additional computational cost. It is shown that its solution, in either Galerkin or finite difference form, converges at the same rates as were proved earlier by Dougl as and Russell for the standard MMOC procedure. Received June 25, 1997 / Revised version received October 6, 1998 / Published online: July 7, 1999     

重建极性连续统理论的基本定律和原理(Ⅶ)——增率型

目的是建立微极连续统增率型的较为完整的运动方程,边界条件和能率方程.为此,先给出较为完整的变形梯度及其逆的定义.接着推导出各种应力率和偶应力率间的新关系式.最后,作为一种特殊情形得到连续统力学的耦合的增率型运动方程、边界条件和能率方程.     

The Cauchy problem for the nonlinear Schrödinger equation in H1

We consider the initial value problem for the nonlinear Schrödinger equation in H¹(Rⁿ). We establish local existence and uniqueness for a wide class of subcritical nonlinearities. The proofs make use of a truncation argument, space-time integrability properties of the linear equation, anda priori estimates derived from the conservation of energy. In particular, we do not need any differentiability property of the nonlinearity with respect to x.Research supported by NSF grants DMS 8201639 and DMS 8703096.     

Solving chemical kinetics equations by explicit methods

The paper presents two versions of numerical methods for solving systems of chemical kinetics which are based on quasi-steady state and partial equilibrium approximations. These methods are explicit and guarantee that the solutions are positive and bounded for an arbitrary period of time. In contrast to the existing methods, they ensure the conservation of mass and the fulfillment of stoichiometric relations, which enables us to increase the timestep considerably while keeping the same accuracy.