In this paper, we consider the linear convection-diffusion 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 implicit-explicit (IMEX) Runge-Kutta methods in time. By using the forward temporal differences and backward temporal differences, respectively, we establish two general frameworks of the energy-method based stability analysis. From here, the fully discrete schemes being considered are shown to have monotonicity stability, i.e. the norm of the numerical solution does not increase in time, under the time step condition , with the convection coefficient , the diffusion coefficient , and the mesh size . The function depends on the specific IMEX temporal method, the polynomial degree of the discrete space, and the mesh regularity parameter. Moreover, the time step condition becomes in the convection-dominated regime and it becomes in the diffusion-dominated regime. The result is improved for a first order IMEX-LDG 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 convection-diffusion equations are convection-dominated in some sub-regions.
more »
« less
This content will become publicly available on July 1, 2026
A local discontinuous Galerkin method for the Novikov equation
In this paper, we propose a local discontinuous Galerkin (LDG) method for the Novikov equation that contains cubic nonlinear high-order derivatives. Flux correction techniques are used to ensure the stability of the numerical scheme. The -norm stability of the general solution and the error estimate for smooth solutions without using any priori assumptions are presented. Numerical examples demonstrate the accuracy and capability of the LDG method for solving the Novikov equation.
more »
« less
- Award ID(s):
- 2309249
- PAR ID:
- 10581292
- Publisher / Repository:
- AMS
- Date Published:
- Journal Name:
- Mathematics of Computation
- Volume:
- 94
- Issue:
- 354
- ISSN:
- 0025-5718
- Page Range / eLocation ID:
- 1603 to 1631
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this paper, we propose and analyze a linear second-order numerical method for solving the Allen-Cahn equation with a general mobility. The proposed fully-discrete scheme is carefully constructed based on the combination of first and second-order backward differentiation formulas with nonuniform time steps for temporal approximation and the central finite difference for spatial discretization. The discrete maximum bound principle is proved of the proposed scheme by using the kernel recombination technique under certain mild constraints on the time steps and the ratios of adjacent time step sizes. Furthermore, we rigorously derive the discrete error estimate and energy stability for the classic constant mobility case and the error estimate for the general mobility case. Various numerical experiments are also presented to validate the theoretical results and demonstrate the performance of the proposed method with a time adaptive strategy.more » « less
-
We consider minimizing harmonic maps from into a closed Riemannian manifold and prove: 1. an extension to of Almgren and Lieb’s linear law. That is, if the fundamental group of the target manifold is finite, we have\[ \]2. an extension of Hardt and Lin’s stability theorem. Namely, assuming that the target manifold is we obtain that the singular set of is stable under small -perturbations of the boundary data. In dimension both results are shown to hold with weaker hypotheses, i.e., only assuming that the trace of our map lies in the fractional space with and satisfying . We also discuss sharpness.more » « less
-
We construct a nonlinear least-squares finite element method for computing the smooth convex solutions of the Dirichlet boundary value problem of the Monge-Ampère equation on strictly convex smooth domains in . It is based on an isoparametric finite element space with exotic degrees of freedom that can enforce the convexity of the approximate solutions.A priorianda posteriorierror estimates together with corroborating numerical results are presented.more » « less
-
The space of vectors of norm in has a natural pseudo-Riemannian metric and a compatible almost complex structure. The group of automorphisms of both of these structures is the split real form . In this paper we consider a class of holomorphic curves in which we call alternating. We show that such curves admit a so called Frenet framing. Using this framing, we show that the space of alternating holomorphic curves which are equivariant with respect to a surface group is naturally parameterized by certain -Higgs bundles. This leads to a holomorphic description of the moduli space as a fibration over Teichmüller space with a holomorphic action of the mapping class group. Using a generalization of Labourie’s cyclic surfaces, we then show that equivariant alternating holomorphic curves are infinitesimally rigid.more » « less
An official website of the United States government
