Source author record

Yulong Xing

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

ResearcherUnclaimed source record

math.NA Numerical Analysis astro-ph.IM astro-ph.SR gr-qc physics.comp-ph

Catalog footprint

What is connected

10works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Discontinuous Galerkin methods for stochastic Maxwell equations with multiplicative noise

In this paper we propose and analyze finite element discontinuous Galerkin methods for the one- and two-dimensional stochastic Maxwell equations with multiplicative noise. The discrete energy law of the semi-discrete DG methods were studied. Optimal error estimate of the semi-discrete method is obtained for the one-dimensional case, and the two-dimensional case on both rectangular meshes and triangular meshes under certain mesh assumptions. Strong Taylor 2.0 scheme is used as the temporal discretization. Both one- and two-dimensional numerical results are presented to validate the theoretical analysis results.

preprint2022arXiv

Energy conserving and well-balanced discontinuous Galerkin methods for the Euler-Poisson equations in spherical symmetry

This paper presents high-order Runge-Kutta (RK) discontinuous Galerkin methods for the Euler-Poisson equations in spherical symmetry. The scheme can preserve a general polytropic equilibrium state and achieve total energy conservation up to machine precision with carefully designed spatial and temporal discretizations. To achieve the well-balanced property, the numerical solutions are decomposed into equilibrium and fluctuation components which are treated differently in the source term approximation. One non-trivial challenge encountered in the procedure is the complexity of the equilibrium state, which is governed by the Lane-Emden equation. For total energy conservation, we present second- and third-order RK time discretization, where different source term approximations are introduced in each stage of the RK method to ensure the conservation of total energy. A carefully designed slope limiter for spherical symmetry is also introduced to eliminate oscillations near discontinuities while maintaining the well-balanced and total-energy-conserving properties. Extensive numerical examples -- including a toy model of stellar core-collapse with a phenomenological equation of state that results in core-bounce and shock formation -- are provided to demonstrate the desired properties of the proposed methods, including the well-balanced property, high-order accuracy, shock capturing capability, and total energy conservation.

preprint2021arXiv

High order well-balanced asymptotic preserving finite difference WENO schemes for the shallow water equations in all Froude numbers

In this paper, high order semi-implicit well-balanced and asymptotic preserving finite difference WENO schemes are proposed for the shallow water equations with a non-flat bottom topography. We consider the Froude number ranging from O(1) to 0, which in the zero Froude limit becomes the "lake equations" for balanced flow without gravity waves. We apply a well-balanced finite difference WENO reconstruction, coupled with a stiffly accurate implicit-explicit (IMEX) Runge-Kutta time discretization. The resulting semi-implicit scheme can be shown to be well-balanced, asymptotic preserving (AP) and asymptotically accurate (AA) at the same time. Both one- and two-dimensional numerical results are provided to demonstrate the high order accuracy, AP property and good performance of the proposed methods in capturing small perturbations of steady state solutions.

preprint2021arXiv

Multi-symplectic discontinuous Galerkin methods for the stochastic Maxwell equations with additive noise

One- and multi-dimensional stochastic Maxwell equations with additive noise are considered in this paper. It is known that such system can be written in the multi-symplectic structure, and the stochastic energy increases linearly in time. High order discontinuous Galerkin methods are designed for the stochastic Maxwell equations with additive noise, and we show that the proposed methods satisfy the discrete form of the stochastic energy linear growth property and preserve the multi-symplectic structure on the discrete level. Optimal error estimate of the semi-discrete DG method is also analyzed. The fully discrete methods are obtained by coupling with symplectic temporal discretizations. One- and two-dimensional numerical results are provided to demonstrate the performance of the proposed methods, and optimal error estimates and linear growth of the discrete energy can be observed for all cases.

preprint2020arXiv

High order finite difference Hermite WENO fast sweeping methods for static Hamilton-Jacobi equations

In this paper, we propose a novel Hermite weighted essentially non-oscillatory (HWENO) fast sweeping method to solve the static Hamilton-Jacobi equations efficiently. During the HWENO reconstruction procedure, the proposed method is built upon a new finite difference fifth order HWENO scheme involving one big stencil and two small stencils. However, one major novelty and difference from the traditional HWENO framework lies in the fact that, we do not need to introduce and solve any additional equations to update the derivatives of the unknown function $ϕ$. Instead, we use the current $ϕ$ and the old spatial derivative of $ϕ$ to update them. The traditional HWENO fast sweeping method is also introduced in this paper for comparison, where additional equations governing the spatial derivatives of $ϕ$ are introduced. The novel HWENO fast sweeping methods are shown to yield great savings in both computational time and storage, which improves the computational efficiency of the traditional HWENO scheme. In addition, a hybrid strategy is also introduced to further reduce computational costs. Extensive numerical experiments are provided to validate the accuracy and efficiency of the proposed approaches.

preprint2019arXiv

On structure-preserving discontinuous Galerkin methods for Hamiltonian partial differential equations: Energy conservation and multi-symplecticity

In this paper, we present and study discontinuous Galerkin (DG) methods for one-dimensional multi-symplectic Hamiltonian partial differential equations. We particularly focus on semi-discrete schemes with spatial discretization only, and show that the proposed DG methods can simultaneously preserve the multi-symplectic structure and energy conservation with a general class of numerical fluxes, which includes the well-known central and alternating fluxes. Applications to the wave equation, the Benjamin-Bona-Mahony equation, the Camassa-Holm equation, the Korteweg-de Vries equation and the nonlinear Schrödinger equation are discussed. Some numerical results are provided to demonstrate the accuracy and long time behavior of the proposed methods. Numerically, we observe that certain choices of numerical fluxes in the discussed class may help achieve better accuracy compared with the commonly used ones including the central fluxes.

preprint2015arXiv

Analysis of mixed interior penalty discontinuous Galerkin methods for the Cahn-Hilliard equation and the Hele-Shaw flow

This paper proposes and analyzes two fully discrete mixed interior penalty discontinuous Galerkin (DG) methods for the fourth order nonlinear Cahn-Hilliard equation. Both methods use the backward Euler method for time discretization and interior penalty discontinuous Galerkin methods for spatial discretization. They differ from each other on how the nonlinear term is treated, one of them is based on fully implicit time-stepping and the other uses the energy-splitting time-stepping. The primary goal of the paper is to prove the convergence of the numerical interfaces of the DG methods to the interface of the Hele-Shaw flow. This is achieved by establishing error estimates that depend on $ε^{-1}$ only in some low polynomial orders, instead of exponential orders. Similar to [14], the crux is to prove a discrete spectrum estimate in the discontinuous Galerkin finite element space. However, the validity of such a result is not obvious because the DG space is not a subspace of the (energy) space $H^1$ and it is larger than the finite element space. This difficult is overcome by a delicate perturbation argument which relies on the discrete spectrum estimate in the finite element space proved in \cite{Feng_Prohl04}. Numerical experiment results are also presented to gauge the theoretical results and the performance of the proposed fully discrete mixed DG methods.

preprint2015arXiv

On the advantage of well-balanced schemes for moving-water equilibria of the shallow water equations

This note aims at demonstrating the advantage of moving-water well-balanced schemes over still-water well-balanced schemes for the shallow water equations. We concentrate on numerical examples with solutions near a moving-water equilibrium. For such examples, still-water well-balanced methods are not capable of capturing the small perturbations of the moving-water equilibrium and may generate significant spurious oscillations, unless an extremely refined mesh is used. On the other hand, moving- water well-balanced methods perform well in these tests. The numerical examples in this note clearly demonstrate the importance of utilizing moving-water well-balanced methods for solutions near a moving-water equilibrium.

preprint2014arXiv

Bound-Preserving Discontinuous Galerkin Methods for Conservative Phase Space Advection in Curvilinear Coordinates

We extend the positivity-preserving method of Zhang & Shu (2010, JCP, 229, 3091-3120) to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for non-equilibrium transport problems in general relativity and also in science and engineering applications with specific geometries. The method achieves high-order accuracy using Discontinuous Galerkin (DG) discretization of phase space and strong stability-preserving, Runge-Kutta (SSP-RK) time integration. Special care in taken to ensure that the method preserves strict bounds for the phase space distribution function $f$; i.e., $f\in[0,1]$. The combination of suitable CFL conditions and the use of the high-order limiter proposed in Zhang & Shu (2010) is sufficient to ensure positivity of the distribution function. However, to ensure that the distribution function satisfies the upper bound, the discretization must, in addition, preserve the divergence-free property of the phase space flow. Proofs that highlight the necessary conditions are presented for general curvilinear coordinates, and the details of these conditions are worked out for some commonly used coordinate systems (i.e., spherical polar spatial coordinates in spherical symmetry and cylindrical spatial coordinates in axial symmetry, both with spherical momentum coordinates). Results from numerical experiments --- including one example in spherical symmetry adopting the Schwarzschild metric --- demonstrate that the method achieves high-order accuracy and that the distribution function satisfies the maximum principle.

preprint2010arXiv

Absolutely stable local discontinuous Galerkin methods for the Helmholtz equation with large wave number

Two local discontinuous Galerkin (LDG) methods using some non-standard numerical fluxes are developed for the Helmholtz equation with the first order absorbing boundary condition in the high frequency regime. It is shown that the proposed LDG methods are absolutely stable (hence well-posed) with respect to both the wave number and the mesh size. Optimal order (with respect to the mesh size) error estimates are proved for all wave numbers in the preasymptotic regime. To analyze the proposed LDG methods, they are recasted and treated as (non-conforming) mixed finite element methods. The crux of the analysis is to establish a generalized {\em inf-sup} condition, which holds without any mesh constraint, for each LDG method. The generalized {\em inf-sup} conditions then easily infer the desired absolute stability of the proposed LDG methods. In return, the stability results not only guarantee the well-posedness of the LDG methods but also play a crucial role in the derivation of the error estimates. Numerical experiments, which confirm the theoretical results and compare the proposed two LDG methods, are also presented in the paper.