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. 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. This work is a continuation of our previous
research on energy-stable 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 leap-frog or
trapezoidal time integrators over different frequency ranges using a variety of measures
more »
« less
High Spatial Order Energy Stable FDTD Methods for Maxwell’s Equations in Nonlinear Optical Media in One Dimension
In this paper, we consider electromagnetic (EM) wave propagation in nonlinear
optical media in one spatial dimension. We model the EM wave propagation by the time-
dependent Maxwell’s equations coupled with a system of nonlinear ordinary differential
equations (ODEs) for the response of the medium to the EM waves. The nonlinearity in the
ODEs describes the instantaneous electronic Kerr response and the residual Raman molecular
vibrational response. The ODEs also include the single resonance linear Lorentz dispersion.
For such model, we will design and analyze fully discrete finite difference time domain
(FDTD) methods that have arbitrary (even) order in space and second order in time. It is
challenging to achieve provable stability for fully discrete methods, and this depends on
the choices of temporal discretizations of the nonlinear terms. In Bokil et al. (J Comput
Phys 350:420–452, 2017), we proposed novel modifications of second-order leap-frog and
trapezoidal temporal schemes in the context of discontinuous Galerkin methods to discretize
the nonlinear terms in this Maxwell model. Here, we continue this work by developing
similar time discretizations within the framework of FDTD methods. More specifically, we
design fully discrete modified leap-frog FDTD methods which are proved to be stable under
appropriate CFL conditions. These method can be viewed as an extension of the Yee-FDTD
scheme to this nonlinear Maxwell model. We also design fully discrete trapezoidal FDTD
methods which are proved to be unconditionally stable. The performance of the fully discrete
FDTD methods are demonstrated through numerical experiments involving kink, antikink
waves and third harmonic generation in soliton propagation.
more »
« less
- Award ID(s):
- 1720116
- NSF-PAR ID:
- 10065567
- Date Published:
- Journal Name:
- Journal of scientific computing
- ISSN:
- 0885-7474
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We provide a short introduction to the devising of a special type of methods for numerically approximating the solution of Hamiltonian partial differential equations. These methods use Galerkin space-discretizations which result in a system of ODEs displaying a discrete version of the Hamiltonian structure of the original system. The resulting system of ODEs is then discretized by a symplectic time-marching method. This combination results in high-order accurate, fully discrete methods which can preserve the invariants of the Hamiltonian defining the ODE system. We restrict our attention to linear Hamiltonian systems, as the main results can be obtained easily and directly, and are applicable to many Hamiltonian systems of practical interest including acoustics, elastodynamics, and electromagnetism. After a brief description of the Hamiltonian systems of our interest, we provide a brief introduction to symplectic time-marching methods for linear systems of ODEs which does not require any background on the subject. We describe then the case in which finite-difference space-discretizations are used and focus on the popular Yee scheme (1966) for electromagnetism. Finally, we consider the case of finite-element space discretizations. The emphasis is placed on the conservation properties of the fully discrete schemes. We end by describing ongoing work.more » « less
-
more » « less
Abstract This work focuses on modeling the interaction between an incompressible, viscous fluid and a poroviscoelastic material. The fluid flow is described using the time‐dependent Stokes equations, and the poroelastic material using the Biot model. The viscoelasticity is incorporated in the equations using a linear Kelvin–Voigt model. We introduce two novel, noniterative, partitioned numerical schemes for the coupled problem. The first method uses the second‐order backward differentiation formula (BDF2) for implicit integration, while treating the interface terms explicitly using a second‐order extrapolation formula. The second method is the Crank–Nicolson and Leap‐Frog (CNLF) method, where the Crank–Nicolson method is used to implicitly advance the solution in time, while the coupling terms are explicitly approximated by the Leap‐Frog integration. We show that the BDF2 method is unconditionally stable and uniformly stable in time, while the CNLF method is stable under a CFL condition. Both schemes are validated using numerical simulations. Second‐order convergence in time is observed for both methods. Simulations over a longer period of time show that the errors in the solution remain bounded. Cases when the structure is poroviscoelastic and poroelastic are included in numerical examples.
-
The thermal radiative transfer (TRT) equations form an integro-differential system that describes the propagation and collisional interactions of photons. Computing accurate and efficient numerical solutions TRT are challenging for several reasons, the first of which is that TRT is defined on a high-dimensional phase space that includes the independent variables of time, space, and velocity. In order to reduce the dimensionality of the phase space, classical approaches such as the P$_N$ (spherical harmonics) or the S$_N$ (discrete ordinates) ansatz are often used in the literature. In this work, we introduce a novel approach: the hybrid discrete (H$^T_N$) approximation to the radiative thermal transfer equations. This approach acquires desirable properties of both P$_N$ and S$_N$, and indeed reduces to each of these approximations in various limits: H$^1_N$ $\equiv$ P$_N$ and H$^T_0$ $\equiv$ S$_T$. We prove that H$^T_N$ results in a system of hyperbolic partial differential equations for all $T\ge 1$ and $N\ge 0$. Another challenge in solving the TRT system is the inherent stiffness due to the large timescale separation between propagation and collisions, especially in the diffusive (i.e., highly collisional) regime. This stiffness challenge can be partially overcome via implicit time integration, although fully implicit methods may become computationally expensive due to the strong nonlinearity and system size. On the other hand, explicit time-stepping schemes that are not also asymptotic-preserving in the highly collisional limit require resolving the mean-free path between collisions, making such schemes prohibitively expensive. In this work we develop a numerical method that is based on a nodal discontinuous Galerkin discretization in space, coupled with a semi-implicit discretization in time. In particular, we make use of a second order explicit Runge-Kutta scheme for the streaming term and an implicit Euler scheme for the material coupling term. Furthermore, in order to solve the material energy equation implicitly after each predictor and corrector step, we linearize the temperature term using a Taylor expansion; this avoids the need for an iterative procedure, and therefore improves efficiency. In order to reduce unphysical oscillation, we apply a slope limiter after each time step. Finally, we conduct several numerical experiments to verify the accuracy, efficiency, and robustness of the H$^T_N$ ansatz and the numerical discretizations.more » « less
-
null (Ed.)In this paper, we apply two fully-discrete local discontinuous Galerkin (LDG) methods to the compressible wormhole propagation. We will prove the stability and error estimates of the schemes. Traditional LDG methods use the diffusion term to control of convection term to obtain the stability for some linear equations. However, the variables in wormhole propagation are coupled together and the whole system is highly nonlinear. Therefore, it is extremely difficult to obtain the stability for fully-discrete LDG methods. To fix this gap, we introduce a new auxiliary variable including both the convection and diffusion terms. Moreover, we also construct a special time integration for the porosity, leading to physically relevant numerical approximations and controllable growth rate of the porosity. With a reasonable growth rate, it is possible to handle the time level mismatch in the first-order fully discrete scheme and obtain the stability of the scheme. For the whole system, we will prove that under weak temporal-spatial conditions, the optimal error estimates for the pressure, velocity, porosity and concentration under different norms can be obtained. Numerical experiments are also given to verify the theoretical results.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
