skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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. the L 2 L^2 norm of the numerical solution does not increase in time, under the time step condition τ<#comment/> ≤<#comment/> F ( h / c , d / c 2 ) \tau \le \mathcal {F}(h/c, d/c^2) , with the convection coefficient c c , the diffusion coefficient d d , and the mesh size h h . The function F \mathcal {F} depends on the specific IMEX temporal method, the polynomial degree k k of the discrete space, and the mesh regularity parameter. Moreover, the time step condition becomes τ<#comment/> ≲<#comment/> h / c \tau \lesssim h/c in the convection-dominated regime and it becomes τ<#comment/> ≲<#comment/> d / c 2 \tau \lesssim d/c^2 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
Award ID(s):
2010107
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. We show that if L 1 \mathcal {L}_1 and L 2 \mathcal {L}_2 are linear transformations from Z d \mathbb {Z}^d to Z d \mathbb {Z}^d satisfying certain mild conditions, then, for any finite subset A A of Z d \mathbb {Z}^d , | L 1 A + L 2 A | ≥<#comment/> ( | det ( L 1 ) | 1 / d + | det ( L 2 ) | 1 / d ) d | A | −<#comment/> o ( | A | ) . \begin{equation*} |\mathcal {L}_1 A+\mathcal {L}_2 A|\geq \left ( |\det (\mathcal {L}_1)|^{1/d}+|\det (\mathcal {L}_2)|^{1/d} \right )^d|A|- o(|A|). \end{equation*} This result corrects and confirms the two-summand case of a conjecture of Bukh and is best possible up to the lower-order term for certain choices of L 1 \mathcal {L}_1 and L 2 \mathcal {L}_2 . As an application, we prove a lower bound for | A + λ<#comment/> ⋅<#comment/> A | |A + \lambda \cdot A| when A A is a finite set of real numbers and λ<#comment/> \lambda is an algebraic number. In particular, when λ<#comment/> \lambda is of the form ( p / q ) 1 / d (p/q)^{1/d} for some p , q , d ∈<#comment/> N p, q, d \in \mathbb {N} , each taken as small as possible for such a representation, we show that | A + λ<#comment/> ⋅<#comment/> A | ≥<#comment/> ( p 1 / d + q 1 / d ) d | A | −<#comment/> o ( | A | ) . \begin{equation*} |A + \lambda \cdot A| \geq (p^{1/d} + q^{1/d})^d |A| - o(|A|). \end{equation*} This is again best possible up to the lower-order term and extends a recent result of Krachun and Petrov which treated the case λ<#comment/> = 2 \lambda = \sqrt {2}
    more » « less
  2. We consider minimizing harmonic maps u u from Ω<#comment/> ⊂<#comment/> R n \Omega \subset \mathbb {R}^n into a closed Riemannian manifold N \mathcal {N} and prove: 1. an extension to n ≥<#comment/> 4 n \geq 4 of Almgren and Lieb’s linear law. That is, if the fundamental group of the target manifold N \mathcal {N} is finite, we have\[ H n −<#comment/> 3 ( sing ⁡<#comment/> u ) ≤<#comment/> C ∫<#comment/> ∂<#comment/> Ω<#comment/> | ∇<#comment/> T u | n −<#comment/> 1 d H n −<#comment/> 1 ; \mathcal {H}^{n-3}(\operatorname {sing} u) \le C \int _{\partial \Omega } |\nabla _T u|^{n-1} \,\mathrm {d}\mathcal {H}^{n-1}; \]2. an extension of Hardt and Lin’s stability theorem. Namely, assuming that the target manifold is N = S 2 \mathcal {N}=\mathbb {S}^2 we obtain that the singular set of u u is stable under small W 1 , n −<#comment/> 1 W^{1,n-1} -perturbations of the boundary data. In dimension n = 3 n=3 both results are shown to hold with weaker hypotheses, i.e., only assuming that the trace of our map lies in the fractional space W s , p W^{s,p} with s ∈<#comment/> ( 1 2 , 1 ] s \in (\frac {1}{2},1] and p ∈<#comment/> [ 2 , ∞<#comment/> ) p \in [2,\infty ) satisfying s p ≥<#comment/> 2 sp \geq 2 . We also discuss sharpness. 
    more » « less
  3. Let f f be analytic on [ 0 , 1 ] [0,1] with | f ( k ) ( 1 / 2 ) | ⩽<#comment/> A α<#comment/> k k ! |f^{(k)}(1/2)|\leqslant A\alpha ^kk! for some constants A A and α<#comment/> > 2 \alpha >2 and all k ⩾<#comment/> 1 k\geqslant 1 . We show that the median estimate of μ<#comment/> = ∫<#comment/> 0 1 f ( x ) d x \mu =\int _0^1f(x)\,\mathrm {d} x under random linear scrambling with n = 2 m n=2^m points converges at the rate O ( n −<#comment/> c log ⁡<#comment/> ( n ) ) O(n^{-c\log (n)}) for any c > 3 log ⁡<#comment/> ( 2 ) / π<#comment/> 2 ≈<#comment/> 0.21 c> 3\log (2)/\pi ^2\approx 0.21 . We also get a super-polynomial convergence rate for the sample median of 2 k −<#comment/> 1 2k-1 random linearly scrambled estimates, when k / m k/m is bounded away from zero. When f f has a p p ’th derivative that satisfies a λ<#comment/> \lambda -Hölder condition then the median of means has error O ( n −<#comment/> ( p + λ<#comment/> ) + ϵ<#comment/> ) O( n^{-(p+\lambda )+\epsilon }) for any ϵ<#comment/> > 0 \epsilon >0 , if k →<#comment/> ∞<#comment/> k\to \infty as m →<#comment/> ∞<#comment/> m\to \infty . The proof techniques use methods from analytic combinatorics that have not previously been applied to quasi-Monte Carlo methods, most notably an asymptotic expression from Hardy and Ramanujan on the number of partitions of a natural number. 
    more » « less
  4. In this paper we derive the best constant for the following L ∞<#comment/> L^{\infty } -type Gagliardo-Nirenberg interpolation inequality ‖<#comment/> u ‖<#comment/> L ∞<#comment/> ≤<#comment/> C q , ∞<#comment/> , p ‖<#comment/> u ‖<#comment/> L q + 1 1 −<#comment/> θ<#comment/> ‖<#comment/> ∇<#comment/> u ‖<#comment/> L p θ<#comment/> , θ<#comment/> = p d d p + ( p −<#comment/> d ) ( q + 1 ) , \begin{equation*} \|u\|_{L^{\infty }}\leq C_{q,\infty ,p} \|u\|^{1-\theta }_{L^{q+1}}\|\nabla u\|^{\theta }_{L^p},\quad \theta =\frac {pd}{dp+(p-d)(q+1)}, \end{equation*} where parameters q q and p p satisfy the conditions p > d ≥<#comment/> 1 p>d\geq 1 , q ≥<#comment/> 0 q\geq 0 . The best constant C q , ∞<#comment/> , p C_{q,\infty ,p} is given by C q , ∞<#comment/> , p = θ<#comment/> −<#comment/> θ<#comment/> p ( 1 −<#comment/> θ<#comment/> ) θ<#comment/> p M c −<#comment/> θ<#comment/> d , M c ∫<#comment/> R d u c , ∞<#comment/> q + 1 d x , \begin{equation*} C_{q,\infty ,p}=\theta ^{-\frac {\theta }{p}}(1-\theta )^{\frac {\theta }{p}}M_c^{-\frac {\theta }{d}},\quad M_c≔\int _{\mathbb {R}^d}u_{c,\infty }^{q+1} dx, \end{equation*} where u c , ∞<#comment/> u_{c,\infty } is the unique radial non-increasing solution to a generalized Lane-Emden equation. The case of equality holds when u = A u c , ∞<#comment/> ( λ<#comment/> ( x −<#comment/> x 0 ) ) u=Au_{c,\infty }(\lambda (x-x_0)) for any real numbers A A , λ<#comment/> > 0 \lambda >0 and x 0 ∈<#comment/> R d x_{0}\in \mathbb {R}^d . In fact, the generalized Lane-Emden equation in R d \mathbb {R}^d contains a delta function as a source and it is a Thomas-Fermi type equation. For q = 0 q=0 or d = 1 d=1 , u c , ∞<#comment/> u_{c,\infty } have closed form solutions expressed in terms of the incomplete Beta functions. Moreover, we show that u c , m →<#comment/> u c , ∞<#comment/> u_{c,m}\to u_{c,\infty } and C q , m , p →<#comment/> C q , ∞<#comment/> , p C_{q,m,p}\to C_{q,\infty ,p} as m →<#comment/> + ∞<#comment/> m\to +\infty for d = 1 d=1 , where u c , m u_{c,m} and C q , m , p C_{q,m,p} are the function achieving equality and the best constant of L m L^m -type Gagliardo-Nirenberg interpolation inequality, respectively. 
    more » « less
  5. Let K / Q p K/\mathbb {Q}_p be a finite unramified extension, ρ<#comment/> ¯<#comment/> : G a l ( Q ¯<#comment/> p / K ) →<#comment/> G L n ( F ¯<#comment/> p ) \overline {\rho }:\mathrm {Gal}(\overline {\mathbb {Q}}_p/K)\rightarrow \mathrm {GL}_n(\overline {\mathbb {F}}_p) a continuous representation, and τ<#comment/> \tau a tame inertial type of dimension n n . We explicitly determine, under mild regularity conditions on τ<#comment/> \tau , the potentially crystalline deformation ring R ρ<#comment/> ¯<#comment/> η<#comment/> , τ<#comment/> R^{\eta ,\tau }_{\overline {\rho }} in parallel Hodge–Tate weights η<#comment/> = ( n −<#comment/> 1 , ⋯<#comment/> , 1 , 0 ) \eta =(n-1,\cdots ,1,0) and inertial type τ<#comment/> \tau when theshapeof ρ<#comment/> ¯<#comment/> \overline {\rho } with respect to τ<#comment/> \tau has colength at most one. This has application to the modularity of a class of shadow weights in the weight part of Serre’s conjecture. Along the way we make unconditional the local-global compatibility results of Park and Qian [Mém. Soc. Math. Fr. (N.S.) 173 (2022), pp. vi+150]. 
    more » « less