More Like this

1. Semi-implicit hybrid discrete (H$^T_N$) approximation of thermal radiative transfer

   https://doi.org/https://link.springer.com/article/10.1007/s10915-021-01686-7  

   R.G. McClarren, J.A. Rossmanith ( November 2021 , Journal of scientific computing)

   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 expensivemore »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.« less
2. A local discontinuous Galerkin method for nonlinear parabolic SPDEs

   https://doi.org/10.1051/m2an/2020026  

   Li, Yunzhang ; Shu, Chi-Wang ; Tang, Shanjian ( January 2021 , ESAIM: Mathematical Modelling and Numerical Analysis)

   In this paper, we propose a local discontinuous Galerkin (LDG) method for nonlinear and possibly degenerate parabolic stochastic partial differential equations, which is a high-order numerical scheme. It extends the discontinuous Galerkin (DG) method for purely hyperbolic equations to parabolic equations and shares with the DG method its advantage and flexibility. We prove the L 2 -stability of the numerical scheme for fully nonlinear equations. Optimal error estimates ( O ( h (k+1) )) for smooth solutions of semi-linear stochastic equations is shown if polynomials of degree k are used. We use an explicit derivative-free order 1.5 time discretization scheme to solve the matrix-valued stochastic ordinary differential equations derived from the spatial discretization. Numerical examples are given to display the performance of the LDG method.
3. Quantitative stability and error estimates for optimal transport plans

   https://doi.org/10.1093/imanum/draa045  

   Li, Wenbo ; Nochetto, Ricardo H ( July 2020 , IMA Journal of Numerical Analysis)

   Abstract Optimal transport maps and plans between two absolutely continuous measures $\mu$ and $\nu$ can be approximated by solving semidiscrete or fully discrete optimal transport problems. These two problems ensue from approximating $\mu$ or both $\mu$ and $\nu$ by Dirac measures. Extending an idea from Gigli (2011, On Hölder continuity-in-time of the optimal transport map towards measures along a curve. Proc. Edinb. Math. Soc. (2), 54, 401–409), we characterize how transport plans change under the perturbation of both $\mu$ and $\nu$. We apply this insight to prove error estimates for semidiscrete and fully discrete algorithms in terms of errors solely arising from approximating measures. We obtain weighted $L^2$ error estimates for both types of algorithms with a convergence rate $O(h^{1/2})$. This coincides with the rate in Theorem 5.4 in Berman (2018, Convergence rates for discretized Monge–Ampère equations and quantitative stability of optimal transport. Preprint available at arXiv:1803.00785) for semidiscrete methods, but the error notion is different.
4. High Spatial Order Energy Stable FDTD Methods for Maxwell’s Equations in Nonlinear Optical Media in One Dimension

   Bokil, Vrushali A ; Cheng, Yingda ; Jiang, Yan ; Li, Fengyan ( April 2018 , Journal of scientific computing)

   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 CFLmore »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.« less
5. Stability analysis and error estimates of arbitrary Lagrangian–Eulerian discontinuous Galerkin method coupled with Runge–Kutta time-marching for linear conservation laws

   https://doi.org/10.1051/m2an/2018069  

   Zhou, Lingling ; Xia, Yinhua ; Shu, Chi-Wang ( January 2019 , ESAIM: Mathematical Modelling and Numerical Analysis)

   In this paper, we discuss the stability and error estimates of the fully discrete schemes for linear conservation laws, which consists of an arbitrary Lagrangian–Eulerian discontinuous Galerkin method in space and explicit total variation diminishing Runge–Kutta (TVD-RK) methods up to third order accuracy in time. The scaling arguments and the standard energy analysis are the key techniques used in our work. We present a rigorous proof to obtain stability for the three fully discrete schemes under suitable CFL conditions. With the help of the reference cell, the error equations are easy to establish and we derive the quasi-optimal error estimates in space and optimal convergence rates in time. For the Euler-forward scheme with piecewise constant elements, the second order TVD-RK method with piecewise linear elements and the third order TVD-RK scheme with polynomials of any order, the usual CFL condition is required, while for other cases, stronger time step restrictions are needed for the results to hold true. More precisely, the Euler-forward scheme needs τ ≤ ρh 2 and the second order TVD-RK scheme needs $ \tau \le \rho {h}^{\frac{4}{3}}$ for higher order polynomials in space, where τ and h are the time and maximum space step, respectively, and ρmore »is a positive constant independent of τ and h .« less