A LOW-REGULARITY FOURIER INTEGRATOR FOR THE DAVEY-STEWARTSON II SYSTEM WITH ALMOST MASS CONSERVATION

In this work,we propose a low-regularity Fourier integrator with almost mass conservation to solve the Davey-StewartsonⅡsystem(hyperbolic-elliptic case).Arbitrary order mass convergence could be achieved by the suitable addition of correction terms,while keeping the first order accuracy in H~γ×H^(γ+1)for initial data in H^(γ+1)×H^(γ+1)withγ>1.The main theorem is that,up to some fixed time T,there exist constantsτ_(0)and C depending only on T and‖u‖_(L^(∞)((0,T);H^(γ+1)))such that,for any 0<τ≤τ_(0),we have that‖u(t_(n),·)-u^(n)‖H_γ≤C_(τ),‖v(t_(n),·)-v^(n)‖_(Hγ+1)≤C_(τ),where u^(n)and v^(n)denote the numerical solutions at t_(n)=nτ.Moreover,the mass of the numerical solution M(u^(n))satisfies that|M(u^(n))-M(u_0)|≤Cτ~5.

DISCREPANCY OF THE GLOBAL AIR MASS AND WATER BUDGETS AMONG 20 CMIP5 CLIMATE MODELS

The consistency of global atmospheric mass and water budget performance in 20 state-of-the-art ocean-atmosphere Coupled Model Intercomparison Project Phase 5(CMIP5) coupled models has been assessed in a historical experiment. All the models realistically reproduce a climatological annual mean of global air mass(AM) close to the ERA-Interim AM during 1989-2005. Surprisingly, the global AM in half of the models shows nearly no seasonal variation,which does not agree with the seasonal processes of global precipitable water or water vapor, given the mass conservation constraint. To better understand the inconsistencies, we evaluated the seasonal cycles of global AM tendency and water vapor source(evaporation minus precipitation). The results suggest that the inconsistencies result from the poor balance between global AM tendency and water vapor source based on the global AM budget equation. Moreover, the cross-equatorial dry air mass flux, or hemispheric dry mass divergence, is not well represented in any of the 20 CMIP5 models, which show a poorly matched seasonal cycle and notably larger amplitude, compared with the hemispheric tendencies of dry AM in both the Northern Hemisphere and Southern Hemisphere. Pronounced erroneous estimations of tropical precipitation also occur in these models. We speculate that the large inaccuracy of precipitation and possibly evaporation in the tropics is one of the key factors for the inconsistent cross-equatorial mass flux. A reasonable cross-equatorial mass flux in well-balanced hemispheric air mass and moisture budgets remains a challenge for both reanalysis assimilation systems and climate modeling.

Effects of a Dry-Mass Conserving Dynamical Core on the Simulation of Tropical Cyclones

The accurate forecasting of tropical cyclones(TCs)is a challenging task.The purpose of this study was to investigate the effects of a dry-mass conserving(DMC)hydrostatic global spectral dynamical core on TC simulation.Experiments were conducted with DMC and total(moist)mass conserving(TMC)dynamical cores.The TC forecast performance was first evaluated considering 20 TCs in the West Pacific region observed during the 2020 typhoon season.The impacts of the DMC dynamical core on forecasts of individual TCs were then estimated.The DMC dynamical core improved both the track and intensity forecasts,and the TC intensity forecast improvement was much greater than the TC track forecast improvement.Sensitivity simulations indicated that the DMC dynamical core-simulated TC intensity was stronger regardless of the forecast lead time.In the DMC dynamical core experiments,three-dimensional winds and warm and moist cores were consistently enhanced with the TC intensity.Drier air in the boundary inflow layer was found in the DMC dynamical core experiments at the early simulation times.Water vapor mixing ratio budget analysis indicated that this mainly depended on the simulated vertical velocity.Higher updraft above the boundary layer yielded a drier boundary layer,resulting in surface latent heat flux(SLHF)enhancement,the major energy source of TC intensification.The higher DMC dynamical core-simulated updraft in the inner core caused a higher net surface rain rate,producing higher net internal atmospheric diabatic heating and increasing the TC intensity.These results indicate that the stronger DMC dynamical coresimulated TCs are mainly related to the higher DMC vertical velocity.

CONSERVATION OF MASS FOR A PARTICLEMOVING WITH HIGH VELOCITY

By using the revision of the momentum for a particle moving with high velocity and by investigating the famous Bucherer's experiment of an electron deflecting with high velocity in the electromagnetic fields in 1908, the paper determines that mass of the electron with high velocity is still to observe the law of conservation of mass.

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CONTINUUM THEORIES (Ⅵ)—CONSERVATION LAWS OF MASS AND INERTIA

The purpose is to reestablish the coupled conservation laws, the local conservation equations and the jump conditions of mass and inertia for polar continuum theories. In this connection the new material derivatives of the deformation gradient, the line element, the surface element and the volume element were derived and the generalized Reynolds transport theorem was presented. Combining these conservation laws of mass and inertia with the balance laws of momentum, angular momentum and energy derived in our previous papers of this series, a rather complete system of coupled basic laws and principles for polar continuum theories is constituted on the whole. From this system the coupled nonlocal balance equations of mass, inertia, momentum, angular momentum and energy may be obtained by the usual localization.

Theoretical analysis of minimum metal foil thickness achievable by asymmetric rolling with fixed identical roll diameters

A novel approach is proposed for computing the minimum thickness of a metal foil that can be achieved by asymmetric rolling using rolls with identical diameter. This approach is based on simultaneously solving Tselikov equation for the rolling pressure and the modified Hitchcock equation for the roller flattening. To minimize the effect of the elastic deformation on the equal flow per second during the ultrathin foil rolling process, the law of conservation of mass was employed to compute the proportions of the forward slip, backward slip, and the cross shear zones in the contact arc, and then a formula was derived for computing the minimum thickness for asymmetric rolling. Experiment was conducted to find the foil minimum thickness for 304 steel by asymmetric rolling under the asymmetry ratios of 1.05, 1.15 and 1.30. The experimental results are in good agreement with the calculated ones. It was validated that the proposed formula can be used to calculate the foil minimum thickness under the asymmetric rolling condition.

A local energy-preserving scheme for Zakharov system

In this paper, we propose a local conservation law for the Zakharov system. The property is held in any local time- space region which is independent of the boundary condition and more essential than the global energy conservation law. Based on the rule that the numerical methods should preserve the intrinsic properties as much as possible, we propose a local energy-preserving （LEP） scheme for the system. The merit of the proposed scheme is that the local energy conservation law can be conserved exactly in any time-space region. With homogeneous Dirchlet boundary conditions, the proposed LEP scheme also possesses the discrete global mass and energy conservation laws. The theoretical properties are verified by numerical results.

Stabilized Crouzeix-Raviart element for the coupled Stokes and Darcy problem

This paper introduces a new stabilized finite element method for the coupled Stokes and Darcy problem based on the nonconforming Crouzeix-Raviart element. Optimal error estimates for the fluid velocity and pressure are derived. A numerical example is presented to verify the theoretical predictions.

Finite Element Methods for Coupled Stokes and Darcy Problems

We derived and analyzed a new numerical scheme for the coupled Stokes and Darcy problems by using H（div） conforming elements in the entire domain. The approach employs the mixed finite element method for the Darcy equations and a stabilized H（div） finite element method for the Stokes equations. Optimal error estimates for the fluid velocity and pressure are derived. The finite element solutions from the new scheme not only feature a full satisfaction of the continuity equation, which is highly demanded in scientific computing, but also satisfy the mass conservation.

A BLOCK-CENTERED UPWIND APPROXIMATION OF THE SEMICONDUCTOR DEVICE PROBLEM ON A DYNAMICALLY CHANGING MESH

The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.

A MIXED FINITE ELEMENT AND CHARACTERISTIC MIXED FINITE ELEMENT FOR INCOMPRESSIBLE MISCIBLE DARCY-FORCHHEIMER DISPLACEMENT AND NUMERICAL ANALYSIS

In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order.The concentration equation is solved by the method of mixed finite element,where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method.The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation.In actual computations the characteristics adopts a large time step without any loss of accuracy.The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux.In order to derive the optimal 3/2-order error estimate in L^(2) norm,a post-processing technique is included in the approximation to the scalar unknowns.Numerical experiments are illustrated finally to validate theoretical analysis and efficiency.This method can be used to solve such an important problem.

CONVERGENCE ANALYSIS OF MIXED VOLUME ELEMENT-CHARACTERISTIC MIXED VOLUME ELEMENT FOR THREE-DIMENSIONAL CHEMICAL OIL-RECOVERY SEEPAGE COUPLED PROBLEM

The physical model is described by a seepage coupled system for simulating numerically three-dimensional chemical oil recovery, whose mathematical description includes three equations to interpret main concepts. The pressure equation is a nonlinear parabolic equation, the concentration is defined by a convection-diffusion equation and the saturations of different components are stated by nonlinear convection-diffusion equations. The transport pressure appears in the concentration equation and saturation equations in the form of Darcy velocity, and controls their processes. The flow equation is solved by the conservative mixed volume element and the accuracy is improved one order for approximating Darcy velocity. The method of characteristic mixed volume element is applied to solve the concentration, where the diffusion is discretized by a mixed volume element method and the convection is treated by the method of characteristics. The characteristics can confirm strong computational stability at sharp fronts and it can avoid numerical dispersion and nonphysical oscillation. The scheme can adopt a large step while its numerical results have small time-truncation error and high order of accuracy. The mixed volume element method has the law of conservation on every element for the diffusion and it can obtain numerical solutions of the concentration and adjoint vectors. It is most important in numerical simulation to ensure the physical conservative nature. The saturation different components are obtained by the method of characteristic fractional step difference. The computational work is shortened greatly by decomposing a three-dimensional problem into three successive one-dimensional problems and it is completed easily by using the algorithm of speedup. Using the theory and technique of a priori estimates of differential equations, we derive an optimal second order estimates in 12 norm. Numerical examples are given to show the effectiveness and practicability and the method is testified as a powerful tool to solve the important problems.

Spline Solution for the Nonlinear Schrödinger Equation

We develop an exponential spline interpolation method to solve the nonlinear Schr&oumldinger equation. The truncation error and stability analysis of the method are investigated and the method is shown to be unconditionally stable. The conservation quantities are computed to determine the conservation properties of the problem. We will describe the method and present numerical tests by two problems. The numerical simulations results demonstrate the well performance of the proposed method.

Utilizing neural networks to supplant chemical kinetics tabulation through mass conservation and weighting of species depletion

Artificial Neural Networks(ANNs)have emerged as a powerful tool in combustion simulations to replace memory-intensive tabulation of integrated chemical kinetics.Complex reaction mechanisms,however,present a challenge for standard ANN approaches as modeling multiple species typically suffers from inaccurate predictions on minor species.This paper presents a novel ANN approach which can be applied on complex reaction mechanisms in tabular data form,and only involves training a single ANN for a complete reaction mechanism.The approach incorporates a network architecture that automatically conserves mass and employs a particular loss weighting based on species depletion.Both modifications are used to improve the overall ANN performance and individual prediction accuracies,especially for minor species mass fractions.To validate its effectiveness,the approach is compared to standard ANNs in terms of performance and ANN complexity.Four distinct reaction mechanisms(H_(2),C_(7)H_(16),C_(12)H_(26),OME_(34))are used as a test cases,and results demonstrate that considerable improvements can be achieved by applying both modifications.

Solving coupled nonlinear Schrödinger equations via a direct discontinuous Galerkin method

In this work,we present the direct discontinuous Galerkin(DDG) method for the one-dimensional coupled nonlinear Schrdinger(CNLS) equation.We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system.The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge-Kutta method.Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.

A sharp-interface immersed smoothed point interpolation method with improved mass conservation for fluid-structure interaction problems

To solve the problem of inaccurate boundary identification and to eliminate the spurious pressure oscillation in the previously developed immersed smoothed point interpolation method(IS-PIM),a new sharp-interface IS-PIM combining mass conservation algorithm,called Sharp-ISPIM-Mass,is proposed in this work.Based on the so called sharp-interface method,the technique of quadratic local velocity reconstruction has been developed by combining with the mass conservation algorithm,which enables the present method improve the accuracy of the velocity field and satisfy the mass conservation condition near the boundary field.So the proposed method would not encounter the problem of spurious mass flux.In addition,a new form of FSI force evaluation considering pressure and viscous force to perform a whole function from the fluid domain to fictitious fluid domain is introduced,which makes the present method obtain more accurate results of FSI force than the original one.Through the numerical studies of a number of benchmark examples,the performance of the Sharp-ISPIM-Mass has been examined and illustrated.

Improved conservative level set method for free surface flow simulation

By coupling the standard and the conservative level set methods, an improved conservative level set method is proposed to capture the free surface smoothly with excellent mass conservation properties. The improvement lies in the fact that the surface normal is computed from a signed distance function instead of the Heaviside function. Comparing with the conservative level set method, the inevitable numerical discretization errors to point the surface normal in arbitrary directions could be eliminated, and the instability of the numerical solution could be improved efficiently. The advantage is clear in the straightforward combination of the standard level set and the conservative level set and a little effort is taken in coding compared with other coupled methods. The present method is validated with several well-known benchmark problems, including the 2-D Zalesak＇s disk rotating, the 3-D sphere stretching in deformation vortex and the dam break flow simulation. The results are shown to be in good agreement with the published experimental data and numerical results.

A Conservative Modification to the Ghost Fluid Method for Compressible Multiphase Flows

A conservative modification to the ghost fluid method(GFM)is developed for compressible multiphase flows.The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance.We track the conservative variables near the material interface and use this information to modify the numerical solution for an interfacing cell when the interface has passed the cell.The modification procedure can be used on the GFM with any base schemes.In this paper we use the fifth order finite difference WENO scheme for the spatial discretization and the third order TVD Runge-Kutta method for the time discretization.The level set method is used to capture the interface.Numerical experiments show that the method is at least mass and momentum conservative and is in general comparable in numerical resolution with the original GFM.

A Conservative Gradient Discretization Method for Parabolic Equations

In this paper,we propose a new conservative gradient discretization method(GDM)for one-dimensional parabolic partial differential equations(PDEs).We use the implicit Euler method for the temporal discretization and conservative gradient discretization method for spatial discretization.The method is based on a new cellcentered meshes,and it is locally conservative.It has smaller truncation error than the classical finite volume method on uniform meshes.We use the framework of the gradient discretization method to analyze the stability and convergence.The numerical experiments show that the new method has second-order convergence.Moreover,it is more accurate than the classical finite volume method in flux error,L2 error and L¥error.

Third-order unconditional positivity-preserving schemes for reactive flows keeping both mass and mole balance

In this paper,the previously proposed second-order process-based modified Patankar Runge-Kutta schemes are extended to the third order of accuracy.Owing to the process-based implicit handling of reactive source terms,the mass conservation,mole balance and energy conservation are kept simultaneously while the positivity for the density and pressure is preserved unconditionally even with stiff reaction networks.It is proved that the first-order truncation terms for the Patankar coefficients must be zero to achieve a prior third order of accuracy for most cases.A twostage Patankar procedure for each Runge-Kutta step is designed to eliminate the first-order truncation terms,accomplish the prior third order of accuracy and maximize the Courant number which the total variational diminishing property requires.With the same approach as the second-order schemes,the third-order ones are applied to Euler equations with chemical reactive source terms.Numerical studies including both 1D and 2D ordinary and partial differential equations are conducted to affirm both the prior order of accuracy and the positivity-preserving property for the density and pressure.