Fengyan Li

Yingda Cheng, Irene M. Gamba, Fengyan Li et al.

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.

Vrushali A. Bokil, Yingda Cheng, Yan Jiang et al.

The propagation of electromagnetic waves in general media is modeled by the time-dependent Maxwell's partial differential equations (PDEs), coupled with constitutive laws that describe the response of the media. In this work, we focus on nonlinear optical media whose response is modeled by a system of first order nonlinear ordinary differential equations (ODEs), which include a single resonance linear Lorentz dispersion, and the nonlinearity comes from the instantaneous electronic Kerr response and the residual Raman molecular vibrational response. To design efficient, accurate, and stable computational methods, we apply high order discontinuous Galerkin discretizations in space to the hybrid PDE-ODE Maxwell system with several choices of numerical fluxes, and the resulting semi-discrete methods are shown to be energy stable. Under some restrictions on the strength of the nonlinearity, error estimates are also established. When we turn to fully discrete methods, the challenge to achieve provable stability lies in the temporal discretizations of the nonlinear terms. To overcome this, novel strategies are proposed to treat the nonlinearity in our model within the framework of the second-order leap-frog and implicit trapezoidal time integrators. The performance of the overall algorithms are demonstrated through numerical simulations of kink and antikink waves, and third-harmonic generation in soliton propagation.

Anqi Chen, Fengyan Li, Yingda Cheng

In this paper, we develop an ultra-weak discontinuous Galerkin (DG) method to solve the one-dimensional nonlinear Schrödinger equation. Stability conditions and error estimates are derived for the scheme with a general class of numerical fluxes. The error estimates are based on detailed analysis of the projection operator associated with each individual flux choice. Depending on the parameters, we find out that in some cases, the projection can be defined element-wise, facilitating analysis. In most cases, the projection is global, and its analysis depends on the resulting $2\times2$ block-circulant matrix structures. For a large class of parameter choices, optimal $\textit{a priori}$ $L^2$ error estimates can be obtained. Numerical examples are provided verifying theoretical results.

Kimberly Matsuda, Fengyan Li

Kinetic transport models are mesoscopic mathematical descriptions of the transport of particles as well as their interactions with the background media or among themselves, and they have wide applications in many areas of mathematical physics such as nuclear and biomedical engineering, rarefied gas dynamics, and plasma physics. They are often multi-scale, with different characteristics (e.g. hyperbolic, diffusive) depending on the material properties. As our continuing effort to design and analyze numerical methods for accurate and robust simulation of the multi-scale kinetic transport models, in this work, we consider a linear kinetic transport model, a simplified radiative transfer equation, in a diffusive scaling, and propose and analyze three families of asymptotic preserving (AP) methods. Numerical methods with the AP property, that is to preserve the asymptotic behavior of the models at the discrete level on under-resolved meshes, can work uniformly well to simulate multi-scale models across a wide range of scales. The proposed methods start from the micro-macro decomposition of the model, and involve discontinuous Galerkin (DG) methods in space, the discrete ordinates method (i.e. $S_N$ method) in velocity, and implicit-explicit (IMEX) BDF methods in time, with three different IMEX partitionings. A systematic study, both analytically and computationally, is presented regarding their difference in stability, accuracy, computational complexity and AP property. These methods, with multi-step time integrators, are also compared in terms of their accuracy and efficiency with the ones that only differ in using certain IMEX Runge-Kutta methods in time. Together with our previous developments, the present work further contributes to high order DG AP methods for multi-scale kinetic simulation, especially by utilizing the structure of the micro-macro decomposition of the models.

Yan Jiang, Puttha Sakkaplangkul, Vrushali A. Bokil et al.

In this paper, we consider Maxwell's equations in linear dispersive media described by a single-pole Lorentz model for electronic polarization. We study two classes of commonly used spatial discretizations: finite difference methods (FD) with arbitrary even order accuracy in space and high spatial order discontinuous Galerkin (DG) finite element methods. Both types of spatial discretizations are coupled with second order semi-implicit leap-frog and implicit trapezoidal temporal schemes studied in our previous research [5,6]. By performing detailed dispersion analysis for the semi-discrete and fully discrete schemes, we obtain rigorous quantification of the dispersion error for Lorentz dispersive dielectrics. In particular, comparisons of dispersion error can be made taking into account the model parameters, and mesh sizes in the design of the two types of schemes. The results for the numerical dispersion analysis can guide us in the optimal choice of discretization parameters for the more complicated and nonlinear models. The numerical dispersion analysis of the fully discrete FD and DG schemes, for the dispersive Maxwell model considered in this paper, clearly indicate the dependence of the numerical dispersion errors on spatial and temporal discretizations, their order of accuracy, mesh discretization parameters and model parameters. The results obtained here cannot be arrived at by considering discretizations of Maxwell's equations in free space. In particular, our results contrast the advantages and disadvantages of using high order FD or DG schemes and leap-frog or trapezoidal time integrators over different frequency ranges using a variety of measures of numerical dispersion errors. Finally, we highlight the limitations of the second order accurate temporal discretizations considered.