 Award ID(s):
 2012699
 NSFPAR ID:
 10355684
 Date Published:
 Journal Name:
 Journal of scientific computing
 Volume:
 90
 ISSN:
 15737691
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

In this paper, we consider Maxwell’s equations in linear dispersive media described by a singlepole 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 semiimplicit leapfrog and implicit trapezoidal temporal schemes. By performing detailed dispersion analysis for the semidiscrete 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. This work is a continuation of our previous research on energystable numerical schemes for nonlinear dispersive optical media [6,7]. The results for the numerical dispersion analysis of the reduced linear model, considered in the present paper, 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 leapfrog or trapezoidal time integrators over different frequency ranges using a variety of measuresmore » « less

In this paper, we consider the linear convectiondiffusion equation in one dimension with periodic boundary conditions, and analyze the stability of fully discrete methods that are defined with local discontinuous Galerkin (LDG) methods in space and several implicitexplicit (IMEX) RungeKutta methods in time. By using the forward temporal differences and backward temporal differences, respectively, we establish two general frameworks of the energymethod based stability analysis. From here, the fully discrete schemes being considered are shown to have monotonicity stability, i.e. the
$L^2$ norm of the numerical solution does not increase in time, under the time step condition$\tau \le \mathcal {F}(h/c, d/c^2)$ , with the convection coefficient$c$ , the diffusion coefficient$d$ , and the mesh size$h$ . The function$\mathcal {F}$ depends on the specific IMEX temporal method, the polynomial degree$k$ of the discrete space, and the mesh regularity parameter. Moreover, the time step condition becomes$\tau \lesssim h/c$ in the convectiondominated regime and it becomes$\tau \lesssim d/c^2$ in the diffusiondominated regime. The result is improved for a first order IMEXLDG method. To complement the theoretical analysis, numerical experiments are further carried out, leading to slightly stricter time step conditions that can be used by practitioners. Uniform stability with respect to the strength of the convection and diffusion effects can especially be relevant to guide the choice of time step sizes in practice, e.g. when the convectiondiffusion equations are convectiondominated in some subregions. 
Abstract We present the implementation of a twomomentbased generalrelativistic multigroup radiation transport module in the
G eneralrelativisticmu ltigridnu merical (Gmunu ) code. On top of solving the generalrelativistic magnetohydrodynamics and the Einstein equations with conformally flat approximations, the code solves the evolution equations of the zeroth and firstorder moments of the radiations in the Eulerianframe. An analytic closure relation is used to obtain the higher order moments and close the system. The finitevolume discretization has been adopted for the radiation moments. The advection in spatial space and frequencyspace are handled explicitly. In addition, the radiation–matter interaction terms, which are very stiff in the optically thick region, are solved implicitly. The implicit–explicit Runge–Kutta schemes are adopted for time integration. We test the implementation with a number of numerical benchmarks from frequencyintegrated to frequencydependent cases. Furthermore, we also illustrate the astrophysical applications in hot neutron star and corecollapse supernovae modelings, and compare with other neutrino transport codes. 
In this article, we develop several first order fully discrete Galerkin finite element schemes for the Oldroyd model and establish the corresponding stability results for these numerical schemes with smooth and nonsmooth initial data. The stable mixed finite element method is used to the spatial discretization, and the temporal treatments of the spatial discrete Oldroyd model include the first order implicit, semi‐implicit, implicit/explicit, and explicit schemes. The
‐stability results of the different numerical schemes are provided, where the first‐order implicit and semi‐implicit schemes are the ‐unconditional stable, the implicit/explicit scheme is the ‐almost unconditional stable, and the first order explicit scheme is the ‐conditional stable. Finally, some numerical investigations of the ‐stability results of the considered numerical schemes for the Oldroyd model are provided to verify the established theoretical findings. 
Abstract We present a new class of discontinuous Galerkin methods for the space discretization of the timedependent Maxwell equations whose main feature is the use of time derivatives and/or time integrals in the stabilization part of their numerical traces.These numerical traces are chosen in such a way that the resulting semidiscrete schemes exactly conserve a discrete version of the energy.We introduce four model ways of achieving this and show that, when using the midpoint rule to march in time, the fully discrete schemes also conserve the discrete energy.Moreover, we propose a new threestep technique to devise fully discrete schemes of arbitrary order of accuracy which conserve the energy in time.The first step consists in transforming the semidiscrete scheme into a Hamiltonian dynamical system.The second step consists in applying a symplectic timemarching method to this dynamical system in order to guarantee that the resulting fully discrete method conserves the discrete energy in time.The third and last step consists in reversing the abovementioned transformation to rewrite the fully discrete scheme in terms of the original variables.more » « less