An official website of the United States government
Abstract A local discontinuous Galerkin (LDG) method for approximating large deformations of prestrained plates is introduced and tested on several insightful numerical examples in Bonito et al. (2022, LDG approximation of large deformations of prestrained plates. J. Comput. Phys., 448, 110719). This paper presents a numerical analysis of this LDG method, focusing on the free boundary case. The problem consists of minimizing a fourth-order bending energy subject to a nonlinear and nonconvex metric constraint. The energy is discretized using LDG and a discrete gradient flow is used for computing discrete minimizers. We first show $$\varGamma $$-convergence of the discrete energy to the continuous one. Then we prove that the discrete gradient flow decreases the energy at each step and computes discrete minimizers with control of the metric constraint defect. We also present a numerical scheme for initialization of the gradient flow and discuss the conditional stability of it.
more »
« less
Analysis and application of a local discontinuous Galerkin method for the electromagnetic concentrator model
Abstract In this paper, we develop a local discontinuous Galerkin (LDG) method to simulate the wave propagation in an electromagnetic concentrator. The concentrator model consists of a coupled system of four partial differential equations and one ordinary differential equation. Discrete stability and error estimate are proved for both semi‐discrete and full‐discrete LDG schemes. Numerical results are presented to justify the theoretical analysis and demonstrate the interesting wave concentration property by the electromagnetic concentrator.
more »
« less
- Award ID(s):
- 2011943
- PAR ID:
- 10487128
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Numerical Methods for Partial Differential Equations
- Volume:
- 40
- Issue:
- 2
- ISSN:
- 0749-159X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract In this paper, we focus on a discrete physical model describing granular crystals, whose equations of motion can be described by a system of differential difference equations. After revisiting earlier continuum approximations, we propose a regularized continuum model variant to approximate the discrete granular crystal model through a suitable partial differential equation. We then compute, both analytically and numerically, its travelling wave and periodic travelling wave solutions, in addition to its conservation laws. Next, using the periodic solutions, we describe quantitatively various features of the dispersive shock wave (DSW) by applying Whitham modulation theory and the DSW fitting method. Finally, we perform several sets of systematic numerical simulations to compare the corresponding DSW results with the theoretical predictions and illustrate that the continuum model provides a good approximation of the underlying discrete one.more » « less
-
null (Ed.)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.more » « less
-
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
