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: A linear second-order maximum bound principle-preserving BDF scheme for the Allen-Cahn equation with a general mobility
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 H 1 H^{1} error estimate and energy stability for the classic constant mobility case and the L ∞<#comment/> L^{\infty } 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
Award ID(s):
2109633
PAR ID:
10544595
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:
2515-2542
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We formulate and prove a Conner–Floyd isomorphism for the algebraic K-theory of arbitrary qcqs derived schemes. To that end, we study a stable ∞<#comment/> \infty -category of non- A 1 \mathbb {A}^1 -invariant motivic spectra, which turns out to be equivalent to the ∞<#comment/> \infty -category of fundamental motivic spectra satisfying elementary blowup excision, previously introduced by the first and third authors. We prove that this ∞<#comment/> \infty -category satisfies P 1 \mathbb {P}^1 -homotopy invariance and weighted A 1 \mathbb {A}^1 -homotopy invariance, which we use in place of A 1 \mathbb {A}^1 -homotopy invariance to obtain analogues of several key results from A 1 \mathbb {A}^1 -homotopy theory. These allow us in particular to define a universal oriented motivic E ∞<#comment/> \mathbb {E}_\infty -ring spectrum M G L \mathrm {MGL} . We then prove that the algebraic K-theory of a qcqs derived scheme X X can be recovered from its M G L \mathrm {MGL} -cohomology via a Conner–Floyd isomorphism\[ M G L ∗<#comment/> ∗<#comment/> ( X ) ⊗<#comment/> L Z [ β<#comment/> ±<#comment/> 1 ] ≃<#comment/> K ∗<#comment/> ∗<#comment/> ( X ) , \mathrm {MGL}^{**}(X)\otimes _{\mathrm {L}{}}\mathbb {Z}[\beta ^{\pm 1}]\simeq \mathrm {K}{}^{**}(X), \]where L \mathrm {L}{} is the Lazard ring and K p , q ( X ) = K 2 q −<#comment/> p ( X ) \mathrm {K}{}^{p,q}(X)=\mathrm {K}{}_{2q-p}(X) . Finally, we prove a Snaith theorem for the periodized version of M G L \mathrm {MGL}
    more » « less
  2. 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
  3. 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
  4. 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 = s l 3 \mathfrak g = \mathfrak {sl}_3 equipped with a diagram involution. We construct explicitly a relative braid group action of type A 2 ( 2 ) A_2^{(2)} on the affine ı<#comment/> \imath quantum group U ~<#comment/> ı<#comment/> \widetilde {{\mathbf U}}^\imath . Real and imaginary root vectors for U ~<#comment/> ı<#comment/> \widetilde {{\mathbf U}}^\imath are constructed, and a Drinfeld type presentation of U ~<#comment/> ı<#comment/> \widetilde {{\mathbf U}}^\imath is then established. This provides a new basic ingredient for the Drinfeld type presentation of higher rank quasi-split affine ı<#comment/> \imath quantum groups in the sequels. 
    more » « less
  5. 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