BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Zhenning Cai

Zhenning Cai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math-phmath.MPmath.NAphysics.flu-dynphysics.comp-phquant-phNumerical Analysiscond-mat.stat-mechhep-lat

Catalog footprint

What is connected

22works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Asymptotic analysis of diabatic surface hopping algorithm in the adiabatic and non-adiabatic limits

Surface hopping algorithms, as an important class of quantum dynamics simulation algorithms for non-adiabatic dynamics, are typically performed in the adiabatic representation, which can break down in the presence of ill-defined adiabatic potential energy surfaces (PESs) and adiabatic coupling term. Another issue of surface hopping algorithms is the difficulty in capturing the correct scaling of the transition rate in the Marcus (weak-coupling/non-adiabatic) regime. Though the first issue can be circumvented by exploiting the diabatic representation, diabatic surface hopping algorithms usually lack justification on the theoretical level. We consider the diabatic surface hopping algorithm proposed in [Fang, Lu. Multiscale Model. Simul. 16:4, 1603-1622, 2018] and provide the asymptotic analysis of the transition rate in the Marcus regime that justifies the correct scaling for the spin-boson model. We propose two conditions that guarantee the correctness for general potentials. In the opposite (strong-coupling/adiabatic) regime, we derive the asymptotic behavior of the algorithm that interestingly matches a type of mean-field description. The techniques used here may shed light on the analysis for other diabatic-based algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast Algorithms of Bath Calculations in Simulations of Quantum System-Bath Dynamics

We present fast algorithms for the summation of Dyson series and the inchworm Monte Carlo method for quantum systems that are coupled with harmonic baths. The algorithms are based on evolving the integro-differential equations where the most expensive part comes from the computation of bath influence functionals. To accelerate the computation, we design fast algorithms based on reusing the bath influence functionals computed in the previous time steps to reduce the number of calculations. It is proven that the proposed fast algorithms reduce the number of such calculations by a factor of $O(N)$, where $N$ is the total number of time steps. Numerical experiments are carried out to show the efficiency of the method and to verify the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Differential Equation Based Path Integral for System-Bath Dynamics

We propose the differential equation based path integral (DEBPI) method to simulate the real-time evolution of open quantum systems. In this method, a system of partial differential equations is derived based on the continuation of a classical numerical method called iterative quasi-adiabatic propagator path integral (i-QuAPI). While the resulting system has infinite equations, we introduce a reasonable closure to obtain a series of finite systems. New numerical schemes can be derived by discretizing these differential equations. It is numerically verified that in certain cases, by selecting appropriate systems and applying suitable numerical schemes, the memory cost required in the i-QuAPI method can be significantly reduced.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Regularization of Complex Langevin Method

The complex Langevin method, a numerical method used to compute the ensemble average with a complex partition function, often suffers from runaway instability. We study the regularization of the complex Langevin method via augmenting the action with a stabilization term. Since the regularization introduces biases to the numerical result, two approaches, named 2R and 3R methods, are introduced to recover the unbiased result. The 2R method supplements the regularization with regression to estimate the unregularized ensemble average, and the 3R method reduces the computational cost by coupling the regularization with a reweighting strategy before regression. Both methods can be generalized to the SU(n) theory and are assessed from several perspectives. Several numerical experiments in the lattice field theory are carried out to show the effectiveness of our approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Flows Between Parallel Plates: Analytical Solutions of Regularized 13-Moment Equations for Inverse-Power-Law Models

We study the structure of stationary channel flows predicted by the regularized 13-moment equations. Compared with the previous work [P. Taheri et al., Phys. Fluids, 21 (2009), 017102], we focus on gases whose molecules satisfy the general inverse power law. The analytical solutions are obtained for the semi-linear equations, and the structures of Couette, Fourier, and Poiseuille flows are solved by coupling the general solutions with newly derived boundary conditions. The results show excellent agreement with the reference solution in the slip-flow regime. Our results also show that the R13 equations derived from inverse power law models can have better accuracy than the R13 equations of Maxwell molecules with altered viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hyperbolic Model Reduction for Kinetic Equations

We make a brief historical review to the moment model reduction to the kinetic equations, particularly the Grad's moment method for Boltzmann equation. The focus is on the hyperbolicity of the reduced model, which is essential to the existence of its classical solution as a Cauchy problem. The theory of the framework we developed in last years is then introduced, which may preserve the hyperbolic nature of the kinetic equations with high universality. Some lastest progress on the comparison between models with/without hyperbolicity is presented to validate the hyperbolic moment models for rarefied gases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical Analysis of a High-Order Scheme for Nonlinear Fractional Differential Equations with Uniform Accuracy

We introduce a high-order numerical scheme for fractional ordinary differential equations with the Caputo derivative. The method is developed by dividing the domain into a number of subintervals, and applying the quadratic interpolation on each subinterval. The method is shown to be unconditionally stable, and for general nonlinear equations, the uniform sharp numerical order $3-ν$ can be rigorously proven for sufficiently smooth solutions at all time steps. The proof provides a general guide for proving the sharp order for higher-order schemes in the nonlinear case. Some numerical examples are given to validate our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

How does Gauge Cooling Stabilize Complex Langevin?

We study the mechanism of the gauge cooling technique to stabilize the complex Langevin method in the one-dimensional periodic setting. In this case, we find the exact solutions for the gauge transform which minimizes the Frobenius norm of link variables. Thereby, we derive the underlying stochastic differential equations by continuing the numerical method with gauge cooling, and thus provide a number of insights on the effects of gauge cooling. A specific case study is carried out for the Polyakov loop model in $SU(2)$ theory, in which we show that the gauge cooling may help form a localized distribution to guarantee there is no excursion too far away from the real axis.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the Holway-Weiss Debate: Convergence of the Grad-Moment-Expansion in Kinetic Gas Theory

Moment expansions are used as model reduction technique in kinetic gas theory to approximate the Boltzmann equation. Rarefied gas models based on so-called moment equations became increasingly popular recently. However, in a seminal paper by Holway [Phys. Fluids 7/6, (1965)] a fundamental restriction on the existence of the expansion was used to explain sub-shock behavior of shock profile solutions obtained by moment equations. Later, Weiss [Phys. Fluids 8/6, (1996)] argued that this restriction does not exist. We will revisit and discuss their findings and explain that both arguments have a correct and incorrect part. While a general convergence restriction for moment expansions does exist, it cannot be attributed to sub-shock solutions. We will also discuss the implications of the restriction and give some numerical evidence for our considerations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Regularized 13-Moment Equations for Inverse Power Law Models

We propose a systematic methodology to derive the regularized thirteen-moment equations in the rarefied gas dynamics for a general class of linearized collision models. Detailed expressions of the moment equations are written down for all inverse power law models as well as the hard-sphere model. By linear analysis, we show that the equations are stable near the equilibrium. The models are tested for shock structure problems to show its capability to capture the correct flow structure in strong nonequilibirum.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Preserving Hyperbolicity in Stochastic Galerkin Method for Uncertainty Quantification

We first investigate the structure of the systems derived from the gPC based stochastic Galerkin method for the nonlinear hyperbolic systems with random inputs. This method adopts a generalized Polynomial Chaos (gPC) approximations in the stochastic Galerkin framework, but such approximations to the nonlinear hyperbolic systems do not necessarily yield hyperbolic systems \cite{Lucor2013}. Thus based on the work in \cite{framework}, we propose a framework to carry out the model reduction for the general nonlinear hyperbolic system to derive a final global system. Within this framework, the nonlinear hyperbolic system in one space dimension and the symmetric hyperbolic system in multiple space dimensions are reduced into a symmetric hyperbolic system based on the stochastic Galerkin method. We note that the basis functions in the expansion are not restricted to the random-dependent polynomials as that in gPC method and there is no restriction on the dimensions of the random variables neither.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Framework on Moment Model Reduction for Kinetic Equation

By a further investigation on the structure of the coefficient matrix of the globally hyperbolic regularized moment equations for Boltzmann equation in [Z. Cai, Y. Fan and R. Li, Comm. Math. Sci., 11 (2013), pp. 547-571], we propose a uniform framework to carry out model reduction to general kinetic equations, to achieve certain moment system. With this framework, the underlying reason why the globally hyperbolic regularization in [Z. Cai, Y. Fan and R. Li, Comm. Math. Sci., 11 (2013), pp. 547-571] works is revealed. The even fascinating point is, with only routine calculation, existing models are represented and brand new models are discovered. Even if the study is restricted in the scope of the classical Grad's 13-moment system, new model with global hyperbolicity can be deduced.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Hyperbolicity of 13-Moment System

We point out that the thermodynamic equilibrium is not an interior point of the hyperbolicity region of Grad's 13-moment system. With a compact expansion of the phase density, which is compacter than Grad's expansion, we derived a modified 13-moment system. The new 13-moment system admits the thermodynamic equilibrium as an interior point of its hyperbolicity region. We deduce a concise criterion to ensure the hyperbolicity, thus the hyperbolicity region can be quantitatively depicted.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Globally Hyperbolic Regularization of Grad's Moment System

In this paper, we propose a globally hyperbolic regularization to the general Grad's moment system in multi-dimensional spaces. Systems with moments up to an arbitrary order are studied. The characteristic speeds of the regularized moment system can be analytically given and only depend on the macroscopic velocity and the temperature. The structure of the eigenvalues and eigenvectors of the coefficient matrix is fully clarified. The regularization together with the properties of the resulting moment systems is consistent with the simple one-dimensional case discussed in [1]. Besides, all characteristic waves are proven to be genuinely nonlinear or linearly degenerate, and the studies on the properties of rarefaction waves, contact discontinuities and shock waves are included.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Globally Hyperbolic Regularization of Grad's Moment System in One Dimensional Space

In this paper, we present a regularization to 1D Grad's moment system to achieve global hyperbolicity. The regularization is based on the observation that the characteristic polynomial of the Jacobian of the flux in Grad's moment system is independent of the intermediate moments. The method is not relied on the form of the collision at all, thus this regularization is applicable to the system without collision terms. Moreover, the proposed approach is proved to be the unique one if only the last moment equation is allowed to be alternated to match the condition of non-equilibrium independent characteristic speeds. The hyperbolic structure of the regularized system, including the signal speeds, Riemann invariants and the properties of the characteristic waves including the rarefaction wave, contact discontinuity and shock are provided in the perfect formations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantum Hydrodynamic Model by Moment Closure of Wigner Equation

In this paper, we derive the quantum hydrodynamics models based on the moment closure of the Wigner equation. The moment expansion adopted is of the Grad type firstly proposed in \cite{Grad}. The Grad's moment method was originally developed for the Boltzmann equation. In \cite{Fan_new}, a regularization method for the Grad's moment system of the Boltzmann equation was proposed to achieve the globally hyperbolicity so that the local well-posedness of the moment system is attained. With the moment expansion of the Wigner function, the drift term in the Wigner equation has exactly the same moment representation as in the Boltzmann equation, thus the regularization in \cite{Fan_new} applies. The moment expansion of the nonlocal Wigner potential term in the Wigner equation is turned to be a linear source term, which can only induce very mild growth of the solution. As the result, the local well-posedness of the regularized moment system for the Wigner equation remains as for the Boltzmann equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solving Vlasov Equations Using NRxx Method

In this paper, we propose a moment method to numerically solve the Vlasov equations using the framework of the NRxx method developed in [6, 8, 7] for the Boltzmann equation. Due to the same convection term of the Boltzmann equation and the Vlasov equation, it is very convenient to use the moment expansion in the NRxx method to approximate the distribution function in the Vlasov equations. The moment closure recently presented in [5] is applied to achieve the globally hyperbolicity so that the local well-posedness of the moment system is attained. This makes our simulations using high order moment expansion accessible in the case of the distribution far away from the equilibrium which appears very often in the solution of the Vlasov equations. With the moment expansion of the distribution function, the acceleration in the velocity space results in an ordinary differential system of the macroscopic velocity, thus is easy to be handled. The numerical method we developed can keep both the mass and the momentum conserved. We carry out the simulations of both the Vlasov-Poisson equations and the Vlasov-Poisson-BGK equations to study the linear Landau damping. The numerical convergence is exhibited in terms of the moment number and the spatial grid size, respectively. The variation of discretized energy as well as the dependence of the recurrence time on moment order is investigated. The linear Landau damping is well captured for different wave numbers and collision frequencies. We find that the Landau damping rate linearly and monotonically converges in the spatial grid size. The results are in perfect agreement with the theoretic data in the collisionless case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The NRxx Method for Polyatomic Gases

In this paper, we propose a numerical regularized moment method to solve the Boltzmann equation with ES-BGK collision term to simulate polyatomic gas flows. This method is an extension to the polyatomic case of the method proposed in [9], which is abbreviated as the NRxx method in [8]. Based on the form of the Maxwellian, the Laguerre polynomials of the internal energy parameter are used in the series expansion of the distribution function. We develop for polyatomic gases all the essential techniques needed in the NRxx method, including the efficient projection algorithm used in the numerical flux calculation, the regularization based on the Maxwellian iteration and the order of magnitude method, and the linearization of the regularization term for convenient numerical implementation. Meanwhile, the particular integrator in time for the ES-BGK collision term is put forward. The shock tube simulations with Knudsen numbers from 0.05 up to 5 are presented to demonstrate the validity of our method. Moreover, the nitrogen shock structure problem is included in our numerical experiments for Mach numbers from 1.53 to 6.1.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

NRxx Simulation of Microflows with Shakhov Model

In this paper, we propose a method to simulate the microflows with Shakhov model using the NRxx method developed in [4, 5, 6]. The equation under consideration is the Boltzmann equation with force terms and the Shakhov model is adopted to achieve the correct Prandtl number. As the focus of this paper, we derive a uniform framework for different order moment systems on the wall boundary conditions, which is a major difficulty in the moment methods. Numerical examples for both steady and unsteady problems are presented to show the convergence in the number of moments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical Regularized Moment Method for High Mach Number Flow

This paper is a continuation of our earlier work \cite{NRxx} in which a numerical moment method with arbitrary order of moments was presented. However, the computation may break down during the calculation of the structure of a shock wave with Mach number $M_0 \geqslant 3$. In this paper, we concentrate on the regularization of the moment systems. First, we apply the Maxwell iteration to the infinite moment system and determine the magnitude of each moment with respect to the Knudsen number. After that, we obtain the approximation of high order moments and close the moment systems by dropping some high-order terms. Linearization is then performed to obtain a very simple regularization term, thus it is very convenient for numerical implementation. To validate the new regularization, the shock structure of low order systems is computed with different number of moments and different shock Mach numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An Efficient NRxx Method for Boltzmann-BGK Equation

In \cite{NRxx}, we proposed a numerical regularized moment method of arbitrary order (abbreviated as NRxx method) for Boltzmann-BGK equation, which makes numerical simulation using very large number of moments possible. In this paper, we are further exploring the efficiency of NRxx method with techniques including the 2nd order HLL flux with linear reconstruction to improve spatial accuracy, the RKC schemes to relieve the time step length constraint by the regularization terms, and the revised Strang splitting to calculate convective and diffusive terms only once without loss of accuracy. It is validated by the numerical results that the overall efficiency is significantly improved and the convergence order is kept well.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Regularized Moment Method of Arbitrary Order for Boltzmann-BGK Equation

We introduce a numerical method for solving Grad's moment equations or regularized moment equations for arbitrary order of moments. In our algorithm, we do not need explicitly the moment equations. As an instead, we directly start from the Boltzmann equation and perform Grad's moment method \cite{Grad} and the regularization technique \cite{Struchtrup2003} numerically. We define a conservative projection operator and propose a fast implementation which makes it convenient to add up two distributions and provides more efficient flux calculations compared with the classic method using explicit expressions of flux functions. For the collision term, the BGK model is adopted so that the production step can be done trivially based on the Hermite expansion. Extensive numerical examples for one- and two-dimensional problems are presented. Convergence in moments can be validated by the numerical results for different number of moments.

0 citations0 reviewsNo code linkedOpen access