BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ruo Li

Ruo Li contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phcond-mat.stat-mechmath-phmath.APmath.MPmath.OCphysics.flu-dynphysics.plasm-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

20 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