The hh-transforms of positivity preserving semigroups and their associated Markov processes are investigated in this paper. In particular, it is shown that any quasi-regular positivity preserving coercive form is hh-a...The hh-transforms of positivity preserving semigroups and their associated Markov processes are investigated in this paper. In particular, it is shown that any quasi-regular positivity preserving coercive form is hh-associated with a pair of special standard processes which are in weak duality.展开更多
This paper presents an improved well-balanced Godunov-type 2-D finite volume model with structured grids to simulate shallow flows with wetting and drying fronts over an irregular topography. The intercell flux is com...This paper presents an improved well-balanced Godunov-type 2-D finite volume model with structured grids to simulate shallow flows with wetting and drying fronts over an irregular topography. The intercell flux is computed using a central upwind scheme, which is a Riemann-problem-solver-free method for hyperbolic conservation laws. The nonnegative reconstruction method for the water depth is implemented to resolve the stationary or wet/dry fronts. The bed slope source term is discretized using a central difference method to capture the static flow state over the irregular topography. Second-order accuracy in space is achieved by using the slope limited linear reconstruction method. With the proposed method, the model can avoid the partially wetting/drying cell problem and maintain the mass conservation. The proposed model is tested and verified against three theoretical benchmark tests and two experimental dam break flows. Further, the model is applied to predict the maximum water level and the flood arrival time at different gauge points for the Malpasset dam break event. The predictions agree well with the numerical results and the measurement data published in literature, which demonstrates that with the present model, a well-balanced state can be achieved and the water depth can be nonnegative when the Courant number is kept less than 0.25.展开更多
A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By usi...A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.展开更多
In order to meet the needs of practical design, an interpolation technique is employed to constrain the shape of surfaces. The method of preserving positivity on the interpolation surface and constraint on interpolati...In order to meet the needs of practical design, an interpolation technique is employed to constrain the shape of surfaces. The method of preserving positivity on the interpolation surface and constraint on interpolating data is also developed. The advantage of this new method is that it can be used to constrain the shape of an interpolating surface only by selecting suitable parameters, and numerical examples are presented to show the performance of the method.展开更多
In this work,spatial second order positivity preserving characteristic blockcentered finite difference methods are proposed for solving convection dominated diffusion problems.By using a conservative piecewise parabol...In this work,spatial second order positivity preserving characteristic blockcentered finite difference methods are proposed for solving convection dominated diffusion problems.By using a conservative piecewise parabolic interpolation with positive constraint,the temporal first order scheme is shown to conserve mass exactly and preserve the positivity property of solution.Taking advantage of characteristics,there is no strict restriction on time steps.The scheme is extended to temporal second order by using a particular extrapolation along the characteristics.To restore solution positivity,a mass conservative local limiter is introduced and verified to keep second order accuracy.Numerical examples are carried out to demonstrate the performance of proposed methods.展开更多
In this paper we propose and analyze a second order accurate numericalscheme for the Cahn-Hilliard equation with logarithmic Flory Huggins energy potential. A modified Crank-Nicolson approximation is applied to the l...In this paper we propose and analyze a second order accurate numericalscheme for the Cahn-Hilliard equation with logarithmic Flory Huggins energy potential. A modified Crank-Nicolson approximation is applied to the logarithmic nonlinear term, while the expansive term is updated by an explicit second order AdamsBashforth extrapolation, and an alternate temporal stencil is used for the surface diffusion term. A nonlinear artificial regularization term is added in the numerical scheme,which ensures the positivity-preserving property, i.e., the numerical value of the phasevariable is always between -1 and 1 at a point-wise level. Furthermore, an unconditional energy stability of the numerical scheme is derived, leveraging the special formof the logarithmic approximation term. In addition, an optimal rate convergence estimate is provided for the proposed numerical scheme, with the help of linearizedstability analysis. A few numerical results, including both the constant-mobility andsolution-dependent mobility flows, are presented to validate the robustness of the proposed numerical scheme.展开更多
This paper generalizes the exponential Runge-Kutta asymptotic preserving(AP)method developed in[G.Dimarco and L.Pareschi,SIAM Numer.Anal.,49(2011),pp.2057–2077]to compute the multi-species Boltzmann equation.Compared...This paper generalizes the exponential Runge-Kutta asymptotic preserving(AP)method developed in[G.Dimarco and L.Pareschi,SIAM Numer.Anal.,49(2011),pp.2057–2077]to compute the multi-species Boltzmann equation.Compared to the single species Boltzmann equation that the method was originally applied on,this set of equation presents a new difficulty that comes from the lack of local conservation laws due to the interaction between different species.Hence extra stiff nonlinear source terms need to be treated properly to maintain the accuracy and the AP property.The method we propose does not contain any nonlinear nonlocal implicit solver,and can capture the hydrodynamic limit with time step and mesh size independent of the Knudsen number.We prove the positivity and strong AP properties of the scheme,which are verified by two numerical examples.展开更多
We present a new Dirichlet boundary condition for the rate-type non-Newtonian diffusive constitutive models. The newly proposed boundary condition is compared with two such well-known and popularly used boundary condi...We present a new Dirichlet boundary condition for the rate-type non-Newtonian diffusive constitutive models. The newly proposed boundary condition is compared with two such well-known and popularly used boundary conditions as the pure Neumann condition and the Dirichlet condition by Sureshkumar and Beris. Our condition is demonstrated to be more stable and robust in a number of numerical test cases. A new Dirichlet boundary condition is implemented in the framework of the finite difference Marker and Cell （MAC） method. In this paper, we also present an energy-stable finite difference MAC scheme that preserves the positivity for the conformation tensor and show how the addition of the diffusion helps the energy-stability in a finite difference MAC scheme-setting.展开更多
In this paper, we propose and analyze a second order accurate in time, masslumped mixed finite element scheme for the Cahn-Hilliard equation with a logarithmic Flory-Huggins energy potential. The standard backward dif...In this paper, we propose and analyze a second order accurate in time, masslumped mixed finite element scheme for the Cahn-Hilliard equation with a logarithmic Flory-Huggins energy potential. The standard backward differentiation formula (BDF)stencil is applied in the temporal discretization. In the chemical potential approximation,both the logarithmic singular terms and the surface diffusion term are treatedimplicitly, while the expansive term is explicitly updated via a second-order Adams-Bashforth extrapolation formula, following the idea of the convex-concave decompositionof the energy functional. In addition, an artificial Douglas-Dupont regularizationterm is added to ensure the energy dissipativity. In the spatial discretization, the masslumped finite element method is adopted. We provide a theoretical justification of theunique solvability of the mass lumped finite element scheme, using a piecewise linearelement. In particular, the positivity is always preserved for the logarithmic argumentsin the sense that the phase variable is always located between -1 and 1. In fact, thesingular nature of the implicit terms and the mass lumped approach play an essentialrole in the positivity preservation in the discrete setting. Subsequently, an unconditionalenergy stability is proven for the proposed numerical scheme. In addition, theconvergence analysis and error estimate of the numerical scheme are also presented.Two numerical experiments are carried out to verify the theoretical properties.展开更多
We discuss the development,verification,and performance of a GPU accelerated discontinuous Galerkin method for the solutions of two dimensional nonlinear shallow water equations.The shallow water equations are hyperbo...We discuss the development,verification,and performance of a GPU accelerated discontinuous Galerkin method for the solutions of two dimensional nonlinear shallow water equations.The shallow water equations are hyperbolic partial differential equations and are widely used in the simulation of tsunami wave propagations.Our algorithms are tailored to take advantage of the single instruction multiple data(SIMD)architecture of graphic processing units.The time integration is accelerated by local time stepping based on a multi-rate Adams-Bashforth scheme.A total variational bounded limiter is adopted for nonlinear stability of the numerical scheme.This limiter is coupled with a mass and momentum conserving positivity preserving limiter for the special treatment of a dry or partially wet element in the triangulation.Accuracy,robustness and performance are demonstrated with the aid of test cases.Furthermore,we developed a unified multi-threading model OCCA.The kernels expressed in OCCA model can be cross-compiled with multi-threading models OpenCL,CUDA,and OpenMP.We compare the performance of the OCCA kernels when cross-compiled with these models.展开更多
The stochastic comparison and preservation of positive correlations for Levy-type processes on R^d are studied under the condition that Levy measure v satisfies f{0〈｜z｜≤1）｜z｜｜v（x, dz） - v（x, d（-z））｜ 〈 ...The stochastic comparison and preservation of positive correlations for Levy-type processes on R^d are studied under the condition that Levy measure v satisfies f{0〈｜z｜≤1）｜z｜｜v（x, dz） - v（x, d（-z））｜ 〈 ∞, x∈ R^d, while the sufficient conditions and necessary ones for them are obtained. In some cases the conditions for stochastic comparison are not only sufficient but also necessary.展开更多
文摘The hh-transforms of positivity preserving semigroups and their associated Markov processes are investigated in this paper. In particular, it is shown that any quasi-regular positivity preserving coercive form is hh-associated with a pair of special standard processes which are in weak duality.
基金Project supported by Natural Science Foundation of Zhejiang Province(Grant No.LR16E090001)the Research Funding of Shenzhen City(Grant No.JCYJ20160425164642646)the Zhejiang Province Science and Technology Research Funding(Grant No.2015C03015)
文摘This paper presents an improved well-balanced Godunov-type 2-D finite volume model with structured grids to simulate shallow flows with wetting and drying fronts over an irregular topography. The intercell flux is computed using a central upwind scheme, which is a Riemann-problem-solver-free method for hyperbolic conservation laws. The nonnegative reconstruction method for the water depth is implemented to resolve the stationary or wet/dry fronts. The bed slope source term is discretized using a central difference method to capture the static flow state over the irregular topography. Second-order accuracy in space is achieved by using the slope limited linear reconstruction method. With the proposed method, the model can avoid the partially wetting/drying cell problem and maintain the mass conservation. The proposed model is tested and verified against three theoretical benchmark tests and two experimental dam break flows. Further, the model is applied to predict the maximum water level and the flood arrival time at different gauge points for the Malpasset dam break event. The predictions agree well with the numerical results and the measurement data published in literature, which demonstrates that with the present model, a well-balanced state can be achieved and the water depth can be nonnegative when the Courant number is kept less than 0.25.
文摘A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.
基金Supported by National nature Science Foundation of China(10771125)Nature Science Foundation of the Shandong Province(Y2007A20)
文摘In order to meet the needs of practical design, an interpolation technique is employed to constrain the shape of surfaces. The method of preserving positivity on the interpolation surface and constraint on interpolating data is also developed. The advantage of this new method is that it can be used to constrain the shape of an interpolating surface only by selecting suitable parameters, and numerical examples are presented to show the performance of the method.
基金Partial support for this research was provided by the National Natural Science Foundation of China(Grant No.11601497)the Natural Science Foundation of Shandong Province(Grant No.ZR2016AB16)Fundamental Research Funds for the Central Universities(No.201964008).
文摘In this work,spatial second order positivity preserving characteristic blockcentered finite difference methods are proposed for solving convection dominated diffusion problems.By using a conservative piecewise parabolic interpolation with positive constraint,the temporal first order scheme is shown to conserve mass exactly and preserve the positivity property of solution.Taking advantage of characteristics,there is no strict restriction on time steps.The scheme is extended to temporal second order by using a particular extrapolation along the characteristics.To restore solution positivity,a mass conservative local limiter is introduced and verified to keep second order accuracy.Numerical examples are carried out to demonstrate the performance of proposed methods.
基金This work is supported in part by the grants NSFC 12071090(W.Chen)NSF DMS-2012669(C.Wang)+2 种基金NSFC 11871159Guangdong Provincial Key Laboratory for Computational Science and Material Design 2019B030301001(X.Wang)NSF DMS-1719854,DMS-2012634(S.Wise).C.Wang also thanks the Key Laboratory of Mathematics for Nonlinear Sciences,Fudan University,for the support.
文摘In this paper we propose and analyze a second order accurate numericalscheme for the Cahn-Hilliard equation with logarithmic Flory Huggins energy potential. A modified Crank-Nicolson approximation is applied to the logarithmic nonlinear term, while the expansive term is updated by an explicit second order AdamsBashforth extrapolation, and an alternate temporal stencil is used for the surface diffusion term. A nonlinear artificial regularization term is added in the numerical scheme,which ensures the positivity-preserving property, i.e., the numerical value of the phasevariable is always between -1 and 1 at a point-wise level. Furthermore, an unconditional energy stability of the numerical scheme is derived, leveraging the special formof the logarithmic approximation term. In addition, an optimal rate convergence estimate is provided for the proposed numerical scheme, with the help of linearizedstability analysis. A few numerical results, including both the constant-mobility andsolution-dependent mobility flows, are presented to validate the robustness of the proposed numerical scheme.
基金supported by the NSF grant DMS-1114546 and NSF Research Network in Mathematical Sciences“KI-Net:Kinetic description of emerging challenges in multiscale problems of natural sciences”X.Y.was partially supported by the startup funding of University of California,Santa Barbara。
文摘This paper generalizes the exponential Runge-Kutta asymptotic preserving(AP)method developed in[G.Dimarco and L.Pareschi,SIAM Numer.Anal.,49(2011),pp.2057–2077]to compute the multi-species Boltzmann equation.Compared to the single species Boltzmann equation that the method was originally applied on,this set of equation presents a new difficulty that comes from the lack of local conservation laws due to the interaction between different species.Hence extra stiff nonlinear source terms need to be treated properly to maintain the accuracy and the AP property.The method we propose does not contain any nonlinear nonlocal implicit solver,and can capture the hydrodynamic limit with time step and mesh size independent of the Knudsen number.We prove the positivity and strong AP properties of the scheme,which are verified by two numerical examples.
文摘We present a new Dirichlet boundary condition for the rate-type non-Newtonian diffusive constitutive models. The newly proposed boundary condition is compared with two such well-known and popularly used boundary conditions as the pure Neumann condition and the Dirichlet condition by Sureshkumar and Beris. Our condition is demonstrated to be more stable and robust in a number of numerical test cases. A new Dirichlet boundary condition is implemented in the framework of the finite difference Marker and Cell （MAC） method. In this paper, we also present an energy-stable finite difference MAC scheme that preserves the positivity for the conformation tensor and show how the addition of the diffusion helps the energy-stability in a finite difference MAC scheme-setting.
基金NSFC(No.12071090)the National Key R&D Program of China(No.2019YFA0709502)+2 种基金Z.R.Zhang is partially supported by NSFC No.11871105 and Science Challenge Project No.TZ2018002C.Wang is partially supported by the NSF DMS-2012269S.M.Wise is partially supported by the NSF DMS-1719854,DMS-2012634.
文摘In this paper, we propose and analyze a second order accurate in time, masslumped mixed finite element scheme for the Cahn-Hilliard equation with a logarithmic Flory-Huggins energy potential. The standard backward differentiation formula (BDF)stencil is applied in the temporal discretization. In the chemical potential approximation,both the logarithmic singular terms and the surface diffusion term are treatedimplicitly, while the expansive term is explicitly updated via a second-order Adams-Bashforth extrapolation formula, following the idea of the convex-concave decompositionof the energy functional. In addition, an artificial Douglas-Dupont regularizationterm is added to ensure the energy dissipativity. In the spatial discretization, the masslumped finite element method is adopted. We provide a theoretical justification of theunique solvability of the mass lumped finite element scheme, using a piecewise linearelement. In particular, the positivity is always preserved for the logarithmic argumentsin the sense that the phase variable is always located between -1 and 1. In fact, thesingular nature of the implicit terms and the mass lumped approach play an essentialrole in the positivity preservation in the discrete setting. Subsequently, an unconditionalenergy stability is proven for the proposed numerical scheme. In addition, theconvergence analysis and error estimate of the numerical scheme are also presented.Two numerical experiments are carried out to verify the theoretical properties.
基金The authors gratefully acknowledge travel grants from Pan-American Advanced Studies Institute,grant from DOE and ANL(ANL Subcontract No.1F-32301 on DOE grant No.DE-AC02-06CH11357)grant from ONR(Award No.N00014-13-1-0873)fellowships from Ken Kennedy Institute of technology at Rice University and support from Shell(Shell Agreement No.PT22584),NVIDIA,and AMD.The authors also acknowledge Dr.Jesse Chan for fruitful discussions during the preparation of this manuscript.
文摘We discuss the development,verification,and performance of a GPU accelerated discontinuous Galerkin method for the solutions of two dimensional nonlinear shallow water equations.The shallow water equations are hyperbolic partial differential equations and are widely used in the simulation of tsunami wave propagations.Our algorithms are tailored to take advantage of the single instruction multiple data(SIMD)architecture of graphic processing units.The time integration is accelerated by local time stepping based on a multi-rate Adams-Bashforth scheme.A total variational bounded limiter is adopted for nonlinear stability of the numerical scheme.This limiter is coupled with a mass and momentum conserving positivity preserving limiter for the special treatment of a dry or partially wet element in the triangulation.Accuracy,robustness and performance are demonstrated with the aid of test cases.Furthermore,we developed a unified multi-threading model OCCA.The kernels expressed in OCCA model can be cross-compiled with multi-threading models OpenCL,CUDA,and OpenMP.We compare the performance of the OCCA kernels when cross-compiled with these models.
基金Supported in part by Creative Research Group Fund of the National Natural Science Foundation of China (No. 10121101)the "985" Project from the Ministry of Education in China
文摘The stochastic comparison and preservation of positive correlations for Levy-type processes on R^d are studied under the condition that Levy measure v satisfies f{0〈｜z｜≤1）｜z｜｜v（x, dz） - v（x, d（-z））｜ 〈 ∞, x∈ R^d, while the sufficient conditions and necessary ones for them are obtained. In some cases the conditions for stochastic comparison are not only sufficient but also necessary.
