<p>In this paper, we propose and analyze a linear, structure-preserving method for solving the Allen–Cahn equation based on the second-order backward differentiation formula (BDF2) with variable time steps and the scalar auxiliary variable (SAV) approach. To this end, we first design a novel and essential auxiliary functional that serves twofold functions: (i) ensuring that a first-order approximation to the auxiliary variable, which is essentially important for deriving the unconditional energy dissipation law, does not affect the second-order temporal accuracy of the phase function <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_3077_Article_IEq1.gif" Format="GIF" Height="17" Rendition="HTML" Resolution="72" Type="Linedraw" Width="16" /> </InlineMediaObject> <EquationSource Format="TEX">\(\phi \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ϕ</mi> </math></EquationSource> </InlineEquation>; and (ii) enabling the development of effective stabilization terms that facilitate the construction of maximum bound principle (MBP)-preserving linear methods. Combining this novel functional with a standard central difference stencil, we then propose a linear, second-order variable-step BDF2 type stabilized exponential SAV scheme, referred to as BDF2-sESAV-I. The scheme is proven to preserve both the discrete modified energy dissipation law under the temporal stepsize ratio <InlineEquation ID="IEq2"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_3077_Article_IEq2.gif" Format="GIF" Height="19" Rendition="HTML" Resolution="72" Type="Linedraw" Width="216" /> </InlineMediaObject> <EquationSource Format="TEX">\( 0&lt; r_{k} := \tau _{k}/\tau _{k-1} &lt; 4.864 - \delta \)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>0</mn> <mo>&lt;</mo> <msub> <mi>r</mi> <mi>k</mi> </msub> <mo>:</mo> <mo>=</mo> <msub> <mi>τ</mi> <mi>k</mi> </msub> <mo stretchy="false">/</mo> <msub> <mi>τ</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>&lt;</mo> <mn>4.864</mn> <mo>-</mo> <mi>δ</mi> </mrow> </math></EquationSource> </InlineEquation> with a positive constant <InlineEquation ID="IEq3"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_3077_Article_IEq3.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="12" /> </InlineMediaObject> <EquationSource Format="TEX">\(\delta \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>δ</mi> </math></EquationSource> </InlineEquation> and the MBP under <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_3077_Article_IEq4.gif" Format="GIF" Height="19" Rendition="HTML" Resolution="72" Type="Linedraw" Width="120" /> </InlineMediaObject> <EquationSource Format="TEX">\( 0&lt; r_{k} &lt; 1 + \sqrt{2} \)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>0</mn> <mo>&lt;</mo> <msub> <mi>r</mi> <mi>k</mi> </msub> <mo>&lt;</mo> <mn>1</mn> <mo>+</mo> <msqrt> <mn>2</mn> </msqrt> </mrow> </math></EquationSource> </InlineEquation>. Moreover, the approximation of the original energy by the modified energy is analyzed. By leveraging the kernel recombination technique, optimal <InlineEquation ID="IEq5"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_3077_Article_IEq5.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="23" /> </InlineMediaObject> <EquationSource Format="TEX">\( H^{1}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>H</mi> <mn>1</mn> </msup> </math></EquationSource> </InlineEquation>- and <InlineEquation ID="IEq6"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_3077_Article_IEq6.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="26" /> </InlineMediaObject> <EquationSource Format="TEX">\( L^{\infty }\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>L</mi> <mi>∞</mi> </msup> </math></EquationSource> </InlineEquation>-norm error estimates for the variable-step BDF2-sESAV-I scheme are rigorously established. Numerical examples are carried out to verify the theoretical results and demonstrate the effectiveness and efficiency of the proposed scheme.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Energy Dissipation Law and Maximum Bound Principle-Preserving Linear BDF2 Schemes with Variable Steps for the Allen–Cahn Equation

  • Bingyin Zhang,
  • Hongfei Fu,
  • Rihui Lan,
  • Shusen Xie

摘要

In this paper, we propose and analyze a linear, structure-preserving method for solving the Allen–Cahn equation based on the second-order backward differentiation formula (BDF2) with variable time steps and the scalar auxiliary variable (SAV) approach. To this end, we first design a novel and essential auxiliary functional that serves twofold functions: (i) ensuring that a first-order approximation to the auxiliary variable, which is essentially important for deriving the unconditional energy dissipation law, does not affect the second-order temporal accuracy of the phase function \(\phi \) ϕ ; and (ii) enabling the development of effective stabilization terms that facilitate the construction of maximum bound principle (MBP)-preserving linear methods. Combining this novel functional with a standard central difference stencil, we then propose a linear, second-order variable-step BDF2 type stabilized exponential SAV scheme, referred to as BDF2-sESAV-I. The scheme is proven to preserve both the discrete modified energy dissipation law under the temporal stepsize ratio \( 0< r_{k} := \tau _{k}/\tau _{k-1} < 4.864 - \delta \) 0 < r k : = τ k / τ k - 1 < 4.864 - δ with a positive constant \(\delta \) δ and the MBP under \( 0< r_{k} < 1 + \sqrt{2} \) 0 < r k < 1 + 2 . Moreover, the approximation of the original energy by the modified energy is analyzed. By leveraging the kernel recombination technique, optimal \( H^{1}\) H 1 - and \( L^{\infty }\) L -norm error estimates for the variable-step BDF2-sESAV-I scheme are rigorously established. Numerical examples are carried out to verify the theoretical results and demonstrate the effectiveness and efficiency of the proposed scheme.