Source author record

Fengyan Li

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

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

9 published item(s)

preprint2022arXiv

A learning-based projection method for model order reduction of transport problems

The Kolmogorov $n$-width of the solution manifolds of transport-dominated problems can decay slowly. As a result, it can be challenging to design efficient and accurate reduced order models (ROMs) for such problems. To address this issue, we propose a new learning-based projection method to construct nonlinear adaptive ROMs for transport problems. The construction follows the offline-online decomposition. In the offline stage, we train a neural network to construct adaptive reduced basis dependent on time and model parameters. In the online stage, we project the solution to the learned reduced manifold. Inheriting the merits from both deep learning and the projection method, the proposed method is more efficient than the conventional linear projection-based methods, and may reduce the generalization error of a solely learning-based ROM. Unlike some learning-based projection methods, the proposed method does not need to take derivatives of the neural network in the online stage.

preprint2022arXiv

An Energy-Based Discontinuous Galerkin Method with Tame CFL Numbers for the Wave Equation

We extend and analyze the energy-based discontinuous Galerkin method for second order wave equations on staggered and structured meshes. By combining spatial staggering with local time-stepping near boundaries, the method overcomes the typical numerical stiffness associated with high order piecewise polynomial approximations. In one space dimension with periodic boundary conditions and suitably chosen numerical fluxes, we prove bounds on the spatial operators that establish stability for CFL numbers $c \frac {Δt}{h} < C$ independent of order when stability-enhanced explicit time-stepping schemes of matching order are used. For problems on bounded domains and in higher dimensions we demonstrate numerically that one can march explicitly with large time steps at high order temporal and spatial accuracy.

preprint2020arXiv

Asymptotic preserving IMEX-DG-S schemes for linear kinetic transport equations based on Schur complement

We consider a linear kinetic transport equation under a diffusive scaling, that converges to a diffusion equation as the Knudsen number $\varepsilon\rightarrow0$. In [3, 21], to achieve the asymptotic preserving (AP) property and unconditional stability in the diffusive regime with $\varepsilon\ll 1$, numerical schemes are developed based on an additional reformulation of the even-odd or micro-macro decomposed version of the equation. The key of the reformulation is to add a weighted diffusive term on both sides of one equation in the decomposed system. The choice of the weight function, however, is problem-dependent and ad-hoc, and it can affect the performance of numerical simulations. To avoid issues related to the choice of the weight function and still obtain the AP property and unconditional stability in the diffusive regime, we propose in this paper a new family of AP schemes, termed as IMEX-DG-S schemes, directly solving the micro-macro decomposed system without any further reformulation. The main ingredients of the IMEX-DG-S schemes include globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta (RK) temporal discretizations with a new IMEX strategy, discontinuous Galerkin (DG) spatial discretizations, discrete ordinate methods for the velocity space, and the application of the Schur complement to the algebraic form of the schemes to control the overall computational cost. The AP property of the schemes is shown formally. With an energy type stability analysis applied to the first order scheme, and Fourier type stability analysis applied to the first to third order schemes, we confirm the uniform stability of the methods with respect to $\varepsilon$ and the unconditional stability in the diffusive regime. A series of numerical examples are presented to demonstrate the performance of the new schemes.

preprint2020arXiv

Stability-enhanced AP IMEX1-LDG method: energy-based stability and rigorous AP property

In our recent work [22], a family of high order asymptotic preserving (AP) methods, termed as IMEX-LDG methods, are designed to solve some linear kinetic transport equations, including the one-group transport equation in slab geometry and the telegraph equation, in a diffusive scaling. As the Knudsen number $\varepsilon$ goes to zero, the limiting schemes are implicit discretizations to the limiting diffusive equation. Both Fourier analysis and numerical experiments imply the methods are unconditionally stable in the diffusive regime when $\varepsilon\ll1$. In this paper, we develop an energy approach to establish the numerical stability of the IMEX1-LDG method, the sub-family of the methods that is first order accurate in time and arbitrary order in space, for the model with general material properties. Our analysis is the first to simultaneously confirm unconditional stability when $\varepsilon\ll1$ and the uniform stability property with respect to $\varepsilon$. To capture the unconditional stability, a novel discrete energy is introduced by better exploring the contribution of the scattering term in different regimes. A general form of the weight function, introduced to obtain the unconditional stability for $\varepsilon\ll1$, is also for the first time considered in such stability analysis. Based on the uniform stability, a rigorous asymptotic analysis is then carried out to show the AP property.

preprint2014arXiv

Analysis of Asymptotic Preserving DG-IMEX Schemes for Linear Kinetic Transport Equations in a Diffusive Scaling

In this paper, some theoretical aspects will be addressed for the asymptotic preserving DG-IMEX schemes recently proposed in [J. Jang, F. Li, J.-M. Qiu and T. Xiong, submitted, arxiv:1306.0227] for kinetic transport equations under a diffusive scaling. We will focus on the methods that are based on discontinuous Galerkin (DG) spatial discretizations with the $P^k$ polynomial space and a first order IMEX temporal discretization, and apply them to two linear models: the telegraph equation and the one-group transport equation in slab geometry. In particular, we will establish uniform numerical stability with respect to Knudsen number $\varepsilon$ using energy methods, as well as error estimates for any given $\varepsilon$. When $\varepsilon\rightarrow 0$, a rigorous asymptotic analysis of the schemes is also obtained. Though the methods and the analysis are presented for one dimension in space, they can be generalized to higher dimensions directly.

preprint2014arXiv

High Order Asymptotic Preserving Nodal Discontinuous Galerkin IMEX Schemes for the BGK Equation

In this paper, we develop high-order asymptotic preserving (AP) schemes for the BGK equation in a hyperbolic scaling, which leads to the macroscopic models such as the Euler and compressible Navier-Stokes equations in the asymptotic limit. Our approaches are based on the so-called micro-macro formulation of the kinetic equation which involves a natural decomposition of the problem to the equilibrium and the non-equilibrium parts. The proposed methods are formulated for the BGK equation with constant or spatially variant Knudsen number. The new ingredients for the proposed methods to achieve high order accuracy are the following: we introduce discontinuous Galerkin (DG) discretization of arbitrary order of accuracy with nodal Lagrangian basis functions in space; we employ a high order globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta (RK) scheme as time discretization. Two versions of the schemes are proposed: Scheme I is a direct formulation based on the micro-macro decomposition of the BGK equation, while Scheme II, motivated by the asymptotic analysis for the continuous problem, utilizes certain properties of the projection operator. Compared with Scheme I, Scheme II not only has better computational efficiency (the computational cost is reduced by half roughly), but also allows the establishment of a formal asymptotic analysis. Specifically, it is demonstrated that when $0<ε\ll 1$, Scheme II, up to $\mathcal{O}(ε^2)$, becomes a local DG discretization with an explicit RK method for the macroscopic compressible Navier-Stokes equations, a method in a similar spirit to the ones in [Bassi \& Rabey 1997, Cockburn \& Shu 1998]. Numerical results are presented for a wide range of Knudsen number to illustrate the effectiveness and high order accuracy of the methods.

preprint2013arXiv

Discontinuous Galerkin Methods for the Vlasov-Maxwell Equations

Discontinuous Galerkin methods are developed for solving the Vlasov-Maxwell system, methods that are designed to be systematically as accurate as one wants with provable conservation of mass and possibly total energy. Such properties in general are hard to achieve within other numerical method frameworks for simulating the Vlasov-Maxwell system. The proposed scheme employs discontinuous Galerkin discretizations for both the Vlasov and the Maxwell equations, resulting in a consistent description of the distribution function and electromagnetic fields. It is proven, up to some boundary effects, that charge is conserved and the total energy can be preserved with suitable choices of the numerical flux for the Maxwell equations and the underlying approximation spaces. Error estimates are established for several flux choices. The scheme is tested on the streaming Weibel instability: the order of accuracy and conservation properties of the proposed method are verified.

preprint2013arXiv

Error Estimates of Runge-Kutta Discontinuous Galerkin Methods for the Vlasov-Maxwell System

In this paper, error analysis is established for Runge-Kutta discontinuous Galerkin (RKDG) methods to solve the Vlasov-Maxwell system. This nonlinear hyperbolic system describes the time evolution of collisionless plasma particles of a single species under the self-consistent electromagnetic field, and it models many phenomena in both laboratory and astrophysical plasmas. The methods involve a third order TVD Runge-Kutta discretization in time and upwind discontinuous Galerkin discretizations of arbitrary order in phase domain. With the assumption that the exact solution has sufficient regularity, the $L^2$ errors of the particle number density function as well as electric and magnetic fields at any given time $T$ are bounded by $C h^{k+\frac{1}{2}}+Cτ^3$ under a CFL condition $τ/h \leq γ$. Here $k$ is the polynomial degree used in phase space discretization, satisfying $k \geq \left \lceil \frac{d_x + 1}{2} \right \rceil$ (the smallest integer greater than or equal to $\frac{d_x+1}{2}$, with $d_x$ being the dimension of spatial domain), $τ$ is the time step, and $h$ is the maximum mesh size in phase space. Both $C$ and $γ$ are positive constants independent of $h$ and $τ$, and they may depend on the polynomial degree $k$, time $T$, the size of the phase domain, certain mesh parameters, and some Sobolev norms of the exact solution. The analysis can be extended to RKDG methods with other numerical fluxes and to RKDG methods solving relativistic Vlasov-Maxwell equations.

preprint2013arXiv

High Order Asymptotic Preserving DG-IMEX Schemes for Discrete-Velocity Kinetic Equations in a Diffusive Scaling

In this paper, we develop a family of high order asymptotic preserving schemes for some discrete-velocity kinetic equations under a diffusive scaling, that in the asymptotic limit lead to macroscopic models such as the heat equation, the porous media equation, the advection-diffusion equation, and the viscous Burgers equation. Our approach is based on the micro-macro reformulation of the kinetic equation which involves a natural decomposition of the equation to the equilibrium and non-equilibrium parts. To achieve high order accuracy and uniform stability as well as to capture the correct asymptotic limit, two new ingredients are employed in the proposed methods: discontinuous Galerkin spatial discretization of arbitrary order of accuracy with suitable numerical fluxes; high order globally stiffly accurate implicit-explicit Runge-Kutta scheme in time equipped with a properly chosen implicit-explicit strategy. Formal asymptotic analysis shows that the proposed scheme in the limit of epsilon -> 0 is an explicit, consistent and high order discretization for the limiting equation. Numerical results are presented to demonstrate the stability and high order accuracy of the proposed schemes together with their performance in the limit.