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

Ruo Li

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

ResearcherUnclaimed source record
math.NAmath-phmath.MPNumerical Analysisphysics.comp-phmath.APcond-mat.stat-mechmath.OCphysics.flu-dynphysics.plasm-phquant-ph

Catalog footprint

What is connected

41works
11topics
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

41 published item(s)

preprint2026arXiv

A reconstructed discontinuous approximation for distributed elliptic control problems

In this paper, we present and analyze an interior penalty discontinuous Galerkin method for the distributed elliptic optimal control problems. It is based on a reconstructed discontinuous approximation which admits arbitrarily high-order approximation space with only one unknown per element. Applying this method, we develop a proper discretization scheme that approximates the state and adjoint variables in the approximation space. Our main contributions are twofold: (1) the derivation of both a priori and a posteriori error estimates of the $L^2$-norm and the energy norms, and (2) the implementation of an efficiently solvable discrete system, which is solved via a linearly convergent projected gradient descent method. Numerical experiments are provided to verify the convergence order in a priori error estimate and the efficiency of a posteriori error estimate.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A general improvement in the WENO-Z-type schemes

A new type of finite volume WENO schemes for hyperbolic problems was devised in [36] by introducing the order-preserving (OP) criterion. In this continuing work, we extend the OP criterion to the WENO-Z-type schemes. We firstly rewrite the formulas of the Z-type weights in a uniform form from a mapping perspective inspired by extensive numerical observations. Accrodingly, we build the concept of the locally order-preserving (LOP) mapping which is an extension of the order-preserving (OP) mapping and the resultant improved WENO-Z-type schemes are denoted as LOP-GMWENO-X. There are four major advantages of the LOP-GMWENO-X schemes superior to the existing WENO-Z-type schemes. Firstly, the new schemes can amend the serious drawback of the existing WENO-Z-type schemes that most of them suffer from either producing severe spurious oscillations or failing to obtain high resolutions in long calculations of hyperbolic problems with discontinuities. Secondly, they can maintain considerably high resolutions on solving problems with high-order critical points at long output times. Thirdly, they can obtain evidently higher resolution in the region with high-frequency but smooth waves. Finally, they can significantly decrease the post-shock oscillations for simulations of some 2D problems with strong shock waves. Extensive benchmark examples are conducted to illustrate these advantages.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Coupling conditions for linear hyperbolic relaxation systems in two-scales problems

This work is concerned with coupling conditions for linear hyperbolic relaxation systems with multiple relaxation times. In the region with small relaxation time, an equilibrium system can be used for computational efficiency. Under the assumption that the relaxation system satisfies the structural stability condition and the interface is non-characteristic, we derive a coupling condition at the interface to couple the two systems in a domain decomposition setting. We prove the validity by the energy estimate and Laplace transform, which shows how the error of the domain decomposition method depends on the smaller relaxation time and the boundary layer effects. In addition, we propose a discontinuous Galerkin (DG) scheme for solving the interface problem with the derived coupling condition and prove the L2 stability. We validate our analysis on the linearized Carleman model and the linearized Grad's moment system and show the effectiveness of the DG scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On Well-posed Boundary Conditions for the Linear Non-homogeneous Moment Equations in Half-space

We propose a necessary and sufficient condition for the well-posedness of the linear non-homogeneous Grad moment equations in half-space. The Grad moment system is based on Hermite expansion and regarded as an efficient reduction model of the Boltzmann equation. At a solid wall, the moment equations are commonly equipped with a Maxwell-type boundary condition named the Grad boundary condition. We point out that the Grad boundary condition is unstable for the non-homogeneous half-space problem. Thanks to the proposed criteria, we verify the well-posedness of a class of modified boundary conditions. The technique to make sure the existence and uniqueness mainly includes a well-designed preliminary simultaneous transformation of the coefficient matrices and Kreiss' procedure about the linear boundary value problem with characteristic boundaries. The stability is established by a weighted estimate. At the same time, we obtain the analytical expressions of the solution, which may help solve the half-space problem efficiently.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A modified adaptive improved mapped WENO method

We propose a new family of mapped WENO schemes by using several adaptive control functions and a smoothing approximation of the signum function. The proposed schemes introduce the adaptivity and admit an extensive permitted range of the parameters in the mapping functions. Consequently, they have the capacity to achieve optimal convergence rates, even near critical points. Particularly, one of these new schemes with fine-tuned parameters illustrates a significant advantage when solving problems with discontinuities. It produces numerical solutions with high resolution without generating spurious oscillations, especially for long output times.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A new mapped WENO scheme using order-preserving mapping

Existing mapped WENO schemes can hardly prevent spurious oscillations while preserving high resolutions at long output times. We reveal in this paper the essential reason for such phenomena. It is actually caused by that the mapping function in these schemes can not preserve the order of the nonlinear weights of the stencils. The nonlinear weights may be increased for non-smooth stencils and be decreased for smooth stencils. It is then indicated to require the set of mapping functions to be Order-Preserving in mapped WENO schemes. Therefore, we propose a new mapped WENO scheme with a set of mapping functions to be order-preserving which exhibits a remarkable advantage over the mapped WENO schemes in references. For long output time simulations, the new scheme has the capacity to attain high resolutions and avoid spurious oscillations near discontinuities meanwhile.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An efficient mapped WENO scheme using approximate constant mapping

We present a novel mapping approach for WENO schemes through the use of an approximate constant mapping function which is constructed by employing an approximation of the classic signum function. The new approximate constant mapping function is designed to meet the overall criteria for a proper mapping function required in the design of the WENO-PM6 scheme. The WENO-PM6 scheme was proposed to overcome the potential loss of accuracy of the WENO-M scheme which was developed to recover the optimal convergence order of the WENO-JS scheme at critical points. Our new mapped WENO scheme, denoted as WENO-ACM, maintains almost all advantages of the WENO-PM6 scheme, including low dissipation and high resolution, while decreases the number of mathematical operations remarkably in every mapping process leading to a significant improvement of efficiency. The convergence rates of the WENO-ACM scheme have been shown through one-dimensional linear advection equation with various initial conditions. Numerical results of one-dimensional Euler equations for the Riemann problems, the Mach 3 shock-density wave interaction and the Woodward-Colella interacting blastwaves are improved in comparison with the results obtained by the WENO-JS, WENO-M and WENO-PM6 schemes. Numerical experiments with two-dimensional problems as the 2D Riemann problem, the shock-vortex interaction, the 2D explosion problem, the double Mach reflection and the forward-facing step problem modeled via the two dimensional Euler equations have been conducted to demonstrate the high resolution and the effectiveness of the WENO-ACM scheme. The WENO-ACM scheme provides significantly better resolution than the WENO-M scheme and slightly better resolution than the WENO-PM6 scheme, and compared to the WENO-M and WENO-PM6 schemes, the extra computational cost is reduced by more than 83% and 93%, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An extension of the order-preserving mapping to the WENO-Z-type schemes

In our latest studies, by introducing the novel order-preserving (OP) criterion, we have successfully addressed the widely concerned issue of the previously published mapped weighted essentially non-oscillatory (WENO) schemes that it is rather difficult to achieve high resolutions on the premise of removing spurious oscillations for long-run simulations of the hyperbolic systems. In the present study, we extend the OP criterion to the WENO-Z-type schemes as the forementioned issue has also been extensively observed numerically for these schemes. Firstly, we innovatively present the concept of the generalized mapped WENO schemes by rewriting the Z-type weights in a uniform formula from the perspective of the mapping relation. Then, we naturally introduce the OP criterion to improve the WENO-Z-type schemes, and the resultant schemes are denoted as MOP-GMWENO-X. Finally, extensive numerical experiments have been conducted to demonstrate the benefits of these new schemes. We draw the conclusion that, the convergence propoties of the proposed schemes are equivalent to the corresponding WENO-X schemes. The major benefit of the new schemes is that they have the capacity to achieve high resolutions and simultaneously remove spurious oscillations for long simulations. The new schemes have the additional benefit that they can greatly decrease the post-shock oscillations on solving 2D Euler problems with strong shock waves.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Direct Flux Gradient Approximation to Close Moment Model for Kinetic Equations

To close the moment model deduced from kinetic equations, the canonical approach is to provide an approximation to the flux function not able to be depicted by the moments in the reduced model. In this paper, we propose a brand new closure approach with remarkable advantages than the canonical approach. Instead of approximating the flux function, the new approach close the moment model by approximating the flux gradient. Precisely, we approximate the space derivative of the distribution function by an ansatz which is a weighted polynomial, and the derivative of the closing flux is computed by taking the moments of the ansatz. Consequently, the method provides us an improved framework to derive globally hyperbolic moment models, which preserve all those conservative variables in the low order moments. It is shown that the linearized system at the weight function, which is often the local equilibrium, of the moment model deduced by our new approach is automatically coincided with the system deduced from the classical perturbation theory, which can not be satisfied by previous hyperbolic regularization framework. Taking the Boltzmann equation as example, the linearlization of the moment model gives the correct Navier-Stokes-Fourier law same as that the Chapman-Enskog expansion gives. Most existing globally hyperbolic moment models are re-produced by our new approach, and several new models are proposed based on this framework.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A discontinuous least squares finite element method for time-harmonic Maxwell equations

We propose and analyze a discontinuous least squares finite element method for solving the indefinite time-harmonic Maxwell equations. The scheme is based on the $L^2$ norm least squares functional with the weak imposition of the continuity across the interior faces. We minimize the functional over the piecewise polynomial spaces to seek numerical solutions. The method is shown to be stable without any constraint on the mesh size. We prove the convergence orders under both the energy norm and the $L^2$ norm. Numerical results in two and three dimensions are presented to verify the error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Nonlinear Hyperbolic Model for Radiative Transfer Equation in Slab Geometry

Linear models for the radiative transfer equation have been well developed, while nonlinear models are seldom investigated even for slab geometry due to some essential difficulties. We have proposed a moment model in MPN for slab geometry which combines the ideas of the classical PN and MN model. Though the model is far from perfect, it was demonstrated to be quite efficient in numerically approximating the solution of the radiative transfer equation, that we are motivated to further improve this model. Consequently we propose in this paper a new model following the chartmap in MPN with some significant theoretic progresses. The new model is derived with global hyperbolicity, and meanwhile some necessary physical properties are preserved. We give a complete analysis to the characteristic structure and propose a numerical scheme for the new model. Numerical examples are presented to demonstrate the numerical performance of the new model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Nonlinear Moment Model for Radiative Transfer Equation

We derive a nonlinear moment model for radiative transfer equation in 3D space, using the method to derive the nonlinear moment model for the radiative transfer equation in slab geometry. The resulted 3D HMPN model enjoys a list of mathematical advantages, including global hyperbolicity, rotational invariance, physical wave speeds, spectral accuracy, and correct higher-order Eddington approximation. Simulation examples are presented to validate the new model numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A sequential least squares method for elliptic equations in non-divergence form

We develop a new least squares method for solving the second-order elliptic equations in non-divergence form. Two least-squares-type functionals are proposed for solving the equations in two steps. We first obtain a numerical approximation to the gradient in a piecewisely irrotational polynomial space. Then together with the numerical gradient, we seek a numerical solution of the primitive variable in continuous finite element space. The variational setting naturally provides a posteriori error which could be used in an adaptive refinement algorithm. The error estimates in $L^2$ norm and energy norms for both two unknowns are derived. By a series of numerical experiments, we verify the convergence rates and show the efficiency of the adaptive algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Two-fluid Model for Plasma with Prandtl Number Correction

A two-fluid model is derived from the plasma kinetic equations using the moment model reduction method. The moment method we adopt was recently developed with a globally hyperbolic regularization where the moment model attained is locally well-posed in time. Based on the hyperbolic moment model with well-posedness, the Maxwellian iteration method is utilized to get the closure relations for the resulted two-fluid model. By taking the Shakhov collision operator in the Maxwellian iteration, the two-fluid model inherits the correct Prandtl number from the plasma kinetic equations. The new model is formally the same as the five-moment two-fluid model except for the closure relations, where the pressure tensor is anisotropic and the heat flux is presented. This provides the model the capacity to depict problems with anisotropic pressure tensor and large heat flux.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Integrated Quadratic Reconstruction for Finite Volume Schemes to Scalar Conservation Laws in Multiple Dimensions

We proposed a piecewise quadratic reconstruction method in multiple dimensions, which is in an integrated style, for finite volume schemes to scalar conservation laws. This integrated quadratic reconstruction is parameter-free and applicable on flexible grids. We show that the finite volume schemes with the new reconstruction satisfy a local maximum principle with properly setup on time steplength. Numerical examples are presented to show that the proposed scheme attains a third-order accuracy for smooth solutions in both 2D and 3D cases. It is indicated by numerical results that the local maximum principle is helpful to prevent overshoots in numerical solutions.

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

Reconstructed Discontinuous Approximation to Stokes Equation in A Sequential Least Squares Formulation

We propose a new least squares finite element method to solve the Stokes problem with two sequential steps. The approximation spaces are constructed by patch reconstruction with one unknown per element. For the first step, we reconstruct an approximation space consisting of piecewise curl-free polynomials with zero trace. By this space, we minimize a least squaresfunctional toobtain thenumericalapproximationstothe gradientof thevelocityand the pressure. In the second step, we minimize another least squares functional to give the solution to the velocity in the reconstructed piecewise divergence-free space. We derive error estimates for all unknowns under L2 norms and energy norms. Numerical results in two dimensions and three dimensions verify the convergence rates and demonstrate the great flexibility of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Least Squares Method for Linear Elasticity using A Patch Reconstructed Space

We propose a discontinuous least squares finite element method for solving the linear elasticity. The approximation space is obtained by patch reconstruction with only one unknown per element. We apply the L 2 norm least squares principle to the stress-displacement formulation based on discontinuous approximation with normal continuity across the interior faces. The optimal convergence order under the energy norm is attained. Numerical results of linear elasticity are presented to verify the error estimates. In addition to enjoying the advantages of discontinuous Galerkin method, we illustrate the great simplicity in implementation, the robustness and the improved efficiency of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Nonlinear Three-Moment Model for Radiative Transfer in Spherical Symmetry

We study the approximation of the radiative transfer equation with a relatively few moments in the spherically symmetric case. We propose a three-moment model based on choosing the beta distribution as the ansatz for the specific intensity. This ansatz enables our model to capture the anisotropy in the distribution function. The characteristic structure of the Riemann problem of the model is studied in detail. Numerical simulations demonstrate its validity in approximating the radiative transfer equation in the spherically symmetric case and its advantage in approximating highly anisotropic distribution functions in comparison to the $P_n$ method.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Approximation to Singular Quadratic Collision Model in Fokker-Planck-Landau Equation

We propose a Hermite-Galerkin spectral method to numerically solve the spatially homogeneous Fokker-Planck-Landau equation with singular quadratic collision model. To compute the collision model, we adopt a novel approximation formulated by a combination of a simple linear term and a quadratic term very expensive to evaluate. Using the Hermite expansion, the quadratic term is evaluated exactly by calculating the spectral coefficients. To deal with singularities, we make use of Burnett polynomials so that even very singular collision model can be handled smoothly. Numerical examples demonstrate that our method can capture low-order moments with satisfactory accuracy and performance.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Second Order Time Homogenized Model for Sediment Transport

A multi-scale method for the hyperbolic systems governing sediment transport in subcritical case is developed. The scale separation of this problem is due to the fact that the sediment transport is much slower than flow velocity. We first derive a zeroth order homogenized model, and then propose a first order correction. It is revealed that the first order correction for hyperbolic systems has to be applied on the characteristic speed of slow variables in one dimensional case. In two dimensional case, besides the characteristic speed, the source term is also corrected. We develop a second order numerical scheme following the framework of heterogeneous multi-scale method. The numerical results in both one and two dimensional cases demonstrate the effectiveness and efficiency of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Acceleration for Microflow Simulations of High-Order Moment Models by Using Lower-Order Model Correction

We study the acceleration of steady-state computation for microflow, which is modeled by the high-order moment models derived recently from the steady-state Boltzmann equation with BGK-type collision term. By using the lower-order model correction, a novel nonlinear multi-level moment solver is developed. Numerical examples verify that the resulting solver improves the convergence significantly thus is able to accelerate the steady-state computation greatly. The behavior of the solver is also numerically investigated. It is shown that the convergence rate increases, indicating the solver would be more efficient, as the total levels increases. Three order reduction strategies of the solver are considered. Numerical results show that the most efficient order reduction strategy would be $m_{l-1} = \lceil m_{l} / 2 \rceil$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Parity-decomposition and Moment Analysis for Stationary Wigner Equation with Inflow Boundary Conditions

We study the stationary Wigner equation on a bounded, one-dimensional spatial domain with inflow boundary conditions by using the parity decomposition in (Barletti and weifel, Trans. Theory Stat. Phys., 507--520, 2001). The decomposition reduces the half-range, two-point boundary value problem into two decoupled initial value problems of the even part and the odd part. Without using a cutoff approximation around zero velocity, we prove that the initial value problem for the even part is well-posed. For the odd part, we prove the uniqueness of the solution in the odd $L^2$-space by analyzing the moment system. An example is provided to show that how to use the analysis to obtain the solution of the stationary Wigner equation with inflow boundary conditions.

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
preprint2015arXiv

An Approximate $M_2$ Model for Radiative Transfer in Slab Geometry

We propose an approximate second order maximum entropy ($M_2$) model for radiative transfer in slab geometry. The model is based on the ansatz of the specific intensity in the form of a $\Beta$-distribution. This gives us an explicit form in its closure. The closure is very close to that of the maximum entropy, thus an approximation of the $M_2$ model. We prove that the new model is globally hyperbolic, sharing most of the advantages of the maximum entropy closure. Numerical examples illustrate that it provides solutions with satisfactory agreement with the $M_2$ model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Model Reduction of Kinetic Equations by Operator Projection

By a further study of the mechanism of the hyperbolic regularization of the moment system for Boltzmann equation proposed in [Z. Cai, Y. Fan, R. Li, Comm. Math. Sci. 11(2): 547-571, 2013], we point out that the key point is treating the time and space derivative in the same way. Based on this understanding, a uniform framework to derive globally hyperbolic moment systems from kinetic equations using an operator projection method is proposed. The framework is so concise and clear that it can be treated as an algorithm with four inputs to derive hyperbolic moment system by routine calculations. Almost all existing globally hyperbolic moment system can be included in the framework, as well as some new moment system including globally hyperbolic regularized versions of Grad ordered moment system and a multidimensional extension of the quadrature-based moment system.

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

A Nonlinear Multigrid Steady-State Solver for Microflow

We develop a nonlinear multigrid method to solve the steady state of microflow, which is modeled by the high order moment system derived recently for the steady-state Boltzmann equation with ES-BGK collision term. The solver adopts a symmetric Gauss-Seidel iterative scheme nested by a local Newton iteration on grid cell level as its smoother. Numerical examples show that the solver is insensitive to the parameters in the implementation thus is quite robust. It is demonstrated that expected efficiency improvement is achieved by the proposed method in comparison with the direct time-stepping scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence of Semi-discrete Stationary Wigner Equation with Inflow Boundary Conditions

Making use of the Whittaker-Shannon interpolation formula with shifted sampling points, we propose in this paper a well-posed semi-discretization of the stationary Wigner equation with inflow BCs. The convergence of the solutions of the discrete problem to the continuous problem is then analysed, providing certain regularity of the solution of the continuous problem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Globally Hyperbolic Moment System by Generalized Hermite Expansion

In a recent paper [Z.-N. Cai, Y.-W. Fan, and R. Li. Tech Report, Institude of Math, Peking Univeristy(2013)], it was revealed that a modified 13-moment system taking intrinsic heat fluxes as variables, instead of the heat fluxes along the coordinate vectors which is adopted in the classical Grad 13-moment system, attains some additional advantages than the classical Grad 13-moment system, particularly including that the equilibrium is turned to be the interior point of its hyperbolicity region. The modified 13-moment system was actually derived from the generalized Hermite expansion of the distribution function, where the anisotropy of Hermite expansion is specified by the full temperature tensor. We extend the method therein in this paper to high order of generalized Hermite expansion to derive arbitrary order moment systems, and proposed a globally hyperbolic regularization to achieve locally well-posedness similar to the method in [Z. Cai, Y. Fan, and R. Li, Comm. Pure Appl. Math.(online)(2013)]. Furthermore, the structure of the eigen-system of the coefficient matrix and all characteristic waves are fully clarified. The obtained systems provide a systematic class of hydrodynamic models as the refined version of Euler equations, which is gradually approaching the Boltzmann equation with increasing order of the expansion.

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
preprint2013arXiv

Stationary Wigner Equation with Inflow Boundary Conditions: Will a Symmetric Potential Yield a Symmetric Solution?

Based on the well-posedness of the stationary Wigner equation with inflow boundary conditions given in (A. Arnold, H et al. J. Math. Phys., 41, 2000), we prove without any additional prerequisite conditions that the solution of the Wigner equation with symmetric potential and inflow boundary conditions will be symmetric. This improve the result in (D. Taj et al. Europhys. Lett., 74, 2006) which depends on the convergence of solution formulated in the Neumann series. By numerical studies, we present the convergence of the numerical solution to the symmetric profile for three different numerical schemes. This implies that the upwind schemes can also yield a symmetric numerical solution, on the contrary to the argument given in (D. Taj et al. Europhys. Lett., 74, 2006).

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