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

Yingda Cheng

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

ResearcherUnclaimed source record
math.NANumerical AnalysisMachine Learningmath-phmath.APmath.MPphysics.comp-ph

Catalog footprint

What is connected

18works
7topics
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

18 published item(s)

preprint2022arXiv

Machine learning moment closure models for the radiative transfer equation I: directly learning a gradient based closure

In this paper, we take a data-driven approach and apply machine learning to the moment closure problem for radiative transfer equation in slab geometry. Instead of learning the unclosed high order moment, we propose to directly learn the gradient of the high order moment using neural networks. This new approach is consistent with the exact closure we derive for the free streaming limit and also provides a natural output normalization. A variety of benchmark tests, including the variable scattering problem, the Gaussian source problem with both periodic and reflecting boundaries, and the two-material problem, show both good accuracy and generalizability of our machine learning closure model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An adaptive multiresolution discontinuous Galerkin method with artificial viscosity for scalar hyperbolic conservation laws in multidimensions

In this paper, we develop an adaptive multiresolution discontinuous Galerkin (DG) scheme for scalar hyperbolic conservation laws in multidimensions. Compared with previous work for linear hyperbolic equations \cite{guo2016transport, guo2017adaptive}, a class of interpolatory multiwavelets are applied to efficiently compute the nonlinear integrals over elements and edges in DG schemes. The resulting algorithm, therefore can achieve similar computational complexity as the sparse grid DG method for smooth solutions. Theoretical and numerical studies are performed taking into consideration of accuracy and stability with regard to the choice of the interpolatory multiwavelets. Artificial viscosity is added to capture the shock and only acts on the leaf elements taking advantages of the multiresolution representation. Adaptivity is realized by auto error thresholding based on hierarchical surplus. Accuracy and robustness are demonstrated by several numerical tests.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An adaptive multiresolution interior penalty discontinuous Galerkin method for wave equations in second order form

In this paper, we propose a class of adaptive multiresolution (also called adaptive sparse grid) discontinuous Galerkin (DG) methods for simulating scalar wave equations in second order form in space. The two key ingredients of the schemes include an interior penalty DG formulation in the adaptive function space and two classes of multiwavelets for achieving multiresolution. In particular, the orthonormal Alpert's multiwavelets are used to express the DG solution in terms of a hierarchical structure, and the interpolatory multiwavelets are further introduced to enhance computational efficiency in the presence of variable wave speed or nonlinear source. Some theoretical results on stability and accuracy of the proposed method are presented. Benchmark numerical tests in 2D and 3D are provided to validate the performance of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An adaptive multiresolution ultra-weak discontinuous Galerkin method for nonlinear Schrodinger equations

This paper develops a high order adaptive scheme for solving nonlinear Schrodinger equations. The solutions to such equations often exhibit solitary wave and local structures, which makes adaptivity essential in improving the simulation efficiency. Our scheme uses the ultra-weak discontinuous Galerkin (DG) formulation and belongs to the framework of adaptive multiresolution schemes. Various numerical experiments are presented to demonstrate the excellent capability of capturing the soliton waves and the blow-up phenomenon.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A sparse grid discontinuous Galerkin method for high-dimensional transport equations and its application to kinetic simulations

In this paper, we develop a sparse grid discontinuous Galerkin (DG) scheme for transport equations and applied it to kinetic simulations. The method uses the weak formulations of traditional Runge-Kutta DG (RKDG) schemes for hyperbolic problems and is proven to be $L^2$ stable and convergent. A major advantage of the scheme lies in its low computational and storage cost due to the employed sparse finite element approximation space. This attractive feature is explored in simulating Vlasov and Boltzmann transport equations. Good performance in accuracy and conservation is verified by numerical tests in up to four dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Adaptive Multiresoluton Discontinuous Galerkin Method for Time-Dependent Transport Equations in Multi-dimensions

In this paper, we develop an adaptive multiresolution discontinuous Galerkin (DG) scheme for time-dependent transport equations in multi-dimensions. The method is constructed using multiwavlelets on tensorized nested grids. Adaptivity is realized by error thresholding based on the hierarchical surplus, and the Runge-Kutta DG (RKDG) scheme is employed as the reference time evolution algorithm. We show that the scheme performs similarly to a sparse grid DG method when the solution is smooth, reducing computational cost in multi-dimensions. When the solution is no longer smooth, the adaptive algorithm can automatically capture fine local structures. The method is therefore very suitable for deterministic kinetic simulations. Numerical results including several benchmark tests, the Vlasov-Poisson (VP) and oscillatory VP systems are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Asymptotic Preserving Maxwell Solver Resulting in the Darwin Limit of Electrodynamics

In plasma simulations, where the speed of light divided by a characteristic length is at a much higher frequency than other relevant parameters in the underlying system, such as the plasma frequency, implicit methods begin to play an important role in generating efficient solutions in these multi-scale problems. Under conditions of scale separation, one can rescale Maxwell's equations in such a way as to give a magneto static limit known as the Darwin approximation of electromagnetics. In this work, we present a new approach to solve Maxwell's equations based on a Method of Lines Transpose (MOL$^T$) formulation, combined with a fast summation method with computational complexity $O(N\log{N})$, where $N$ is the number of grid points (particles). Under appropriate scaling, we show that the proposed schemes result in asymptotic preserving methods that can recover the Darwin limit of electrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence of discontinuous Galerkin schemes for front propagation with obstacles

We study semi-Lagrangian discontinuous Galerkin (SLDG) and Runge-Kutta discontinuous Galerkin (RKDG) schemes for some front propagation problems in the presence of an obstacle term, modeled by a nonlinear Hamilton-Jacobi equation of the form $\min(u_t + c u_x, u - g(x))=0$, in one space dimension. New convergence results and error bounds are obtained for Lipschitz regular data. These "low regularity" assumptions are the natural ones for the solutions of the studied equations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

One dimensional dissipative Boltzmann equation: measure solutions, cooling rate and self-similar profile

This manuscript investigates the following aspects of the one dimensional dissipative Boltzmann equation associated to variable hard-spheres kernel: (1) we show the optimal cooling rate of the model by a careful study of the system satisfied by the solution's moments, (2) give existence and uniqueness of measure solutions, and (3) prove the existence of a non-trivial self-similar profile, i.e. homogeneous cooling state, after appropriate scaling of the equation. The latter issue is based on compactness tools in the set of Borel measures. More specifically, we apply a dynamical fixed point theorem on a suitable stable set, for the model dynamics, of Borel measures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sparse Grid Discontinuous Galerkin Methods for High-Dimensional Elliptic Equations

This paper constitutes our initial effort in developing sparse grid discontinuous Galerkin (DG) methods for high-dimensional partial differential equations (PDEs). Over the past few decades, DG methods have gained popularity in many applications due to their distinctive features. However, they are often deemed too costly because of the large number of degrees of freedom of the approximation space, which are the main bottleneck for simulations in high dimensions. In this paper, we develop sparse grid DG methods for elliptic equations with the aim of breaking the \emph{curse of dimensionality}. Using a hierarchical basis representation, we construct a sparse finite element approximation space, reducing the degrees of freedom from the standard {$O(h^{-d})$ to $O(h^{-1}|\log_2 h|^{d-1})$} for $d$-dimensional problems, where $h$ is the uniform mesh size in each dimension. Our method, based on the interior penalty (IP) DG framework, can achieve accuracy of $O(h^{k}|\log_2 h|^{d-1})$ in the energy norm, where $k$ is the degree of polynomials used. Error estimates are provided and confirmed by numerical tests in multi-dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy-conserving Discontinuous Galerkin Methods for the Vlasov-Maxwell System

In this paper, we generalize the idea in our previous work for the Vlasov-Ampère (VA) system \cite{cheng_va} and develop energy-conserving discontinuous Galerkin (DG) methods for the Vlasov-Maxwell (VM) system. The VM system is a fundamental model in the simulation of collisionless magnetized plasmas. Compared to \cite{cheng_va}, additional care needs to be taken for both the temporal and spatial discretizations to achieve similar type of conservation when the magnetic field is no longer negligible. Our proposed schemes conserve the total particle number and the total energy at the same time, and therefore can obtain accurate, yet physically relevant solutions. The main components of our methods include second order and above, explicit or implicit energy-conserving temporal discretizations, and DG methods for Vlasov and Maxwell's equations with carefully chosen numerical fluxes. Benchmark numerical tests such as the streaming Weibel instability are provided to validate the accuracy and conservation of the schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Study of the Two-Species Vlasov-Ampère System: Energy-Conserving Schemes and the Current-Driven Ion-Acoustic Instability

In this paper, we propose energy-conserving Eulerian solvers for the two-species Vlasov-Ampère (VA) system and apply the methods to simulate current-driven ion-acoustic instability. The algorithm is generalized from our previous work for the single-species VA system and Vlasov-Maxwell (VM) system. The main feature of the schemes is their ability to preserve the total particle number and total energy on the fully discrete level regardless of mesh size. Those are desired properties of numerical schemes especially for long time simulations with under-resolved mesh. The conservation is realized by explicit and implicit energy-conserving temporal discretizations, and the discontinuous Galerkin (DG) spatial discretizations. We benchmarked our algorithms on a test example to check the one-species limit, and the current-driven ion-acoustic instability. To simulate the current-driven ion-acoustic instability, a slight modification for the implicit method is necessary to fully decouple the split equations. This is achieved by a Gauss-Seidel type iteration technique. Numerical results verified the conservation and performance of our methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A New Discontinuous Galerkin Finite Element Method for Directly Solving the Hamilton-Jacobi Equations

In this paper, we improve upon the discontinuous Galerkin (DG) method for Hamilton-Jacobi (HJ) equation with convex Hamiltonians in (Y. Cheng and C.-W. Shu, J. Comput. Phys. 223:398-415,2007) and develop a new DG method for directly solving the general HJ equations. The new method avoids the reconstruction of the solution across elements by utilizing the Roe speed at the cell interface. Besides, we propose an entropy fix by adding penalty terms proportional to the jump of the normal derivative of the numerical solution. The particular form of the entropy fix was inspired by the Harten and Hyman's entropy fix (A. Harten and J. M. Hyman. J. Comput. Phys. 50(2):235-269, 1983) for Roe scheme for the conservation laws. The resulting scheme is compact, simple to implement even on unstructured meshes, and is demonstrated to work for nonconvex Hamiltonians. Benchmark numerical experiments in one dimension and two dimensions are provided to validate the performance of the method.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy-conserving discontinuous Galerkin methods for the Vlasov-Ampère system

In this paper, we propose energy-conserving numerical schemes for the Vlasov-Ampère (VA) systems. The VA system is a model used to describe the evolution of probability density function of charged particles under self consistent electric field in plasmas. It conserves many physical quantities, including the total energy which is comprised of the kinetic and electric energy. Unlike the total particle number conservation, the total energy conservation is challenging to achieve. For simulations in longer time ranges, negligence of this fact could cause unphysical results, such as plasma self heating or cooling. In this paper, we develop the first Eulerian solvers that can preserve fully discrete total energy conservation. The main components of our solvers include explicit or implicit energy-conserving temporal discretizations, an energy-conserving operator splitting for the VA equation and discontinuous Galerkin finite element methods for the spatial discretizations. We validate our schemes by rigorous derivations and benchmark numerical examples such as Landau damping, two-stream instability and bump-on-tail instability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Study of conservation and recurrence of Runge-Kutta discontinuous Galerkin schemes for Vlasov-Poisson systems

In this paper we consider Runge-Kutta discontinuous Galerkin (RKDG) schemes for Vlasov-Poisson systems that model collisionless plasmas. One-dimensional systems are emphasized. The RKDG method, originally devised to solve conservation laws, is seen to have excellent conservation properties, be readily designed for arbitrary order of accuracy, and capable of being used with a positivity-preserving limiter that guarantees positivity of the distribution functions. The RKDG solver for the Vlasov equation is the main focus, while the electric field is obtained through the classical representation by Green's function for the Poisson equation. A rigorous study of recurrence of the DG methods is presented by Fourier analysis, and the impact of different polynomial spaces and the positivity-preserving limiters on the quality of the solutions is ascertained. Several benchmark test problems, such as Landau damping, two-stream instability and the KEEN (Kinetic Electrostatic Electron Nonlinear) wave, are given.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A discontinuous Galerkin solver for Boltzmann Poisson systems in nano devices

In this paper, we present results of a discontinuous Galerkin (DG) scheme applied to deterministic computations of the transients for the Boltzmann-Poisson system describing electron transport in semiconductor devices. The collisional term models optical-phonon interactions which become dominant under strong energetic conditions corresponding to nano-scale active regions under applied bias. The proposed numerical technique is a finite element method using discontinuous piecewise polynomials as basis functions on unstructured meshes. It is applied to simulate hot electron transport in bulk silicon, in a silicon $n^+$-$n$-$n^+$ diode and in a double gated 12nm MOSFET. Additionally, the obtained results are compared to those of a high order WENO scheme simulation and DSMC (Discrete Simulation Monte Carlo) solvers.

0 citations0 reviewsNo code linkedOpen access