skip to main content


This content will become publicly available on November 1, 2024

Title: Uniform stability for local discontinuous Galerkin methods with implicit-explicit Runge-Kutta time discretizations for linear convection-diffusion equation

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. theL2L^2norm of the numerical solution does not increase in time, under the time step conditionτ<#comment/>≤<#comment/>F(h/c,d/c2)\tau \le \mathcal {F}(h/c, d/c^2), with the convection coefficientcc, the diffusion coefficientdd, and the mesh sizehh. The functionF\mathcal {F}depends on the specific IMEX temporal method, the polynomial degreekkof the discrete space, and the mesh regularity parameter. Moreover, the time step condition becomesτ<#comment/>≲<#comment/>h/c\tau \lesssim h/cin the convection-dominated regime and it becomesτ<#comment/>≲<#comment/>d/c2\tau \lesssim d/c^2in 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
Award ID(s):
2010107
NSF-PAR ID:
10498236
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
AMS
Date Published:
Journal Name:
Mathematics of Computation
Volume:
92
Issue:
344
ISSN:
0025-5718
Page Range / eLocation ID:
2475 to 2513
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. By discretizing an argument of Kislyakov, Naor and Schechtman proved that the 1-Wasserstein metric over the planar grid{0,1,…<#comment/>,n}2\{0,1,\dots , n\}^2hasL1L_1-distortion bounded below by a constant multiple oflog⁡<#comment/>n\sqrt {\log n}. We provide a new “dimensionality” interpretation of Kislyakov’s argument, showing that if{Gn}n=1∞<#comment/>\{G_n\}_{n=1}^\inftyis a sequence of graphs whose isoperimetric dimension and Lipschitz-spectral dimension equal a common numberδ<#comment/>∈<#comment/>[2,∞<#comment/>)\delta \in [2,\infty ), then the 1-Wasserstein metric overGnG_nhasL1L_1-distortion bounded below by a constant multiple of(log⁡<#comment/>|Gn|)1δ<#comment/>(\log |G_n|)^{\frac {1}{\delta }}. We proceed to compute these dimensions for⊘<#comment/>\oslash-powers of certain graphs. In particular, we get that the sequence of diamond graphs{Dn}n=1∞<#comment/>\{\mathsf {D}_n\}_{n=1}^\inftyhas isoperimetric dimension and Lipschitz-spectral dimension equal to 2, obtaining as a corollary that the 1-Wasserstein metric overDn\mathsf {D}_nhasL1L_1-distortion bounded below by a constant multiple oflog⁡<#comment/>|Dn|\sqrt {\log | \mathsf {D}_n|}. This answers a question of Dilworth, Kutzarova, and Ostrovskii and exhibits only the third sequence ofL1L_1-embeddable graphs whose sequence of 1-Wasserstein metrics is notL1L_1-embeddable.

     
    more » « less
  2. IfIIis an ideal in a Gorenstein ringSS, andS/IS/Iis Cohen-Macaulay, then the same is true for any linked idealII’; but such statements hold for residual intersections of higher codimension only under restrictive hypotheses, not satisfied even by ideals as simple as the idealLnL_{n}of minors of a generic2×<#comment/>n2 \times nmatrix whenn>3n>3.

    In this paper we initiate the study of a different sort of Cohen-Macaulay property that holds for certain general residual intersections of the maximal (interesting) codimension, one less than the analytic spread ofII. For example, suppose thatKKis the residual intersection ofLnL_{n}by2n−<#comment/>42n-4general quadratic forms inLnL_{n}. In this situation we analyzeS/KS/Kand show thatIn−<#comment/>3(S/K)I^{n-3}(S/K)is a self-dual maximal Cohen-MacaulayS/KS/K-module with linear free resolution overSS.

    The technical heart of the paper is a result about ideals of analytic spread 1 whose high powers are linearly presented.

     
    more » « less
  3. We show that for primesN,p≥<#comment/>5N, p \geq 5withN≡<#comment/>−<#comment/>1modpN \equiv -1 \bmod p, the class number ofQ(N1/p)\mathbb {Q}(N^{1/p})is divisible bypp. Our methods are via congruences between Eisenstein series and cusp forms. In particular, we show that whenN≡<#comment/>−<#comment/>1modpN \equiv -1 \bmod p, there is always a cusp form of weight22and levelΓ<#comment/>0(N2)\Gamma _0(N^2)whoseℓ<#comment/>\ellth Fourier coefficient is congruent toℓ<#comment/>+1\ell + 1modulo a prime abovepp, for all primesℓ<#comment/>\ell. We use the Galois representation of such a cusp form to explicitly construct an unramified degree-ppextension ofQ(N1/p)\mathbb {Q}(N^{1/p}).

     
    more » « less
  4. Let(R,m)(R,\mathfrak {m})be a Noetherian local ring of dimensiond≥<#comment/>2d\geq 2. We prove that ife(R^<#comment/>red)>1e(\widehat {R}_{red})>1, then the classical Lech’s inequality can be improved uniformly for allm\mathfrak {m}-primary ideals, that is, there existsε<#comment/>>0\varepsilon >0such thate(I)≤<#comment/>d!(e(R)−<#comment/>ε<#comment/>)ℓ<#comment/>(R/I)e(I)\leq d!(e(R)-\varepsilon )\ell (R/I)for allm\mathfrak {m}-primary idealsI⊆<#comment/>RI\subseteq R. This answers a question raised by Huneke, Ma, Quy, and Smirnov [Adv. Math. 372 (2020), pp. 107296, 33]. We also obtain partial results towards improvements of Lech’s inequality when we fix the number of generators ofII.

     
    more » « less
  5. This is the first of our papers on quasi-split affine quantum symmetric pairs(U~<#comment/>(g^<#comment/>),U~<#comment/>ı<#comment/>)\big (\widetilde {\mathbf U}(\widehat {\mathfrak g}), \widetilde {{\mathbf U}}^\imath \big ), focusing on the real rank one case, i.e.,g=sl3\mathfrak g = \mathfrak {sl}_3equipped with a diagram involution. We construct explicitly a relative braid group action of typeA2(2)A_2^{(2)}on the affineı<#comment/>\imathquantum groupU~<#comment/>ı<#comment/>\widetilde {{\mathbf U}}^\imath. Real and imaginary root vectors forU~<#comment/>ı<#comment/>\widetilde {{\mathbf U}}^\imathare constructed, and a Drinfeld type presentation ofU~<#comment/>ı<#comment/>\widetilde {{\mathbf U}}^\imathis then established. This provides a new basic ingredient for the Drinfeld type presentation of higher rank quasi-split affineı<#comment/>\imathquantum groups in the sequels.

     
    more » « less