<p>We focus on establishing a novel convergence analysis of the numerical approximation methods applied to axisymmetric mean curvature flows (MCFs) with genus-1 surfaces, including both isotropic and anisotropic cases. This advancement builds upon and enhances the findings presented in a previous work in Barrett et al. (IMA J Numer Anal 41:1641–1667, 2021), which exists a strict restriction on the time-space step ratio, i.e. <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_2821_Article_IEq1.gif" Format="GIF" Height="21" Rendition="HTML" Resolution="72" Type="Linedraw" Width="71" /> </InlineMediaObject> <EquationSource Format="TEX">\(\Delta t \lesssim \sqrt{h}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi mathvariant="normal">Δ</mi> <mi>t</mi> <mo>≲</mo> <msqrt> <mi>h</mi> </msqrt> </mrow> </math></EquationSource> </InlineEquation>. In this work, we first adopt a novel time-space error splitting technique to remove above restriction under the isotropic case. Furthermore, we discover that this splitting also aids in addressing another unresolved issue present in anisotropic MCFs. Indeed, almost any restrictions on the time and space steps cannot derive the <InlineEquation ID="IEq2"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_2821_Article_IEq2.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>-boundedness of <InlineEquation ID="IEq3"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_2821_Article_IEq3.gif" Format="GIF" Height="26" Rendition="HTML" Resolution="72" Type="Linedraw" Width="52" /> </InlineMediaObject> <EquationSource Format="TEX">\(D_t\vec {X}_{h, \rho }^m\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi>D</mi> <mi>t</mi> </msub> <msubsup> <mover accent="true"> <mi>X</mi> <mo stretchy="false">→</mo> </mover> <mrow> <mi>h</mi> <mo>,</mo> <mi>ρ</mi> </mrow> <mi>m</mi> </msubsup> </mrow> </math></EquationSource> </InlineEquation>, which is required for convergence analysis in the anisotropic case. Given that <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_2821_Article_IEq4.gif" Format="GIF" Height="26" Rendition="HTML" Resolution="72" Type="Linedraw" Width="272" /> </InlineMediaObject> <EquationSource Format="TEX">\(\Vert D_t\vec {X}_{h, \rho }^m\Vert _{L^\infty } \lesssim \Delta t^{-1}h^{-1/2}(\Delta t + h) + 1\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mrow> <mo stretchy="false">‖</mo> </mrow> <msub> <mi>D</mi> <mi>t</mi> </msub> <msubsup> <mover accent="true"> <mi>X</mi> <mo stretchy="false">→</mo> </mover> <mrow> <mi>h</mi> <mo>,</mo> <mi>ρ</mi> </mrow> <mi>m</mi> </msubsup> <msub> <mrow> <mo stretchy="false">‖</mo> </mrow> <msup> <mi>L</mi> <mi>∞</mi> </msup> </msub> <mo>≲</mo> <mi mathvariant="normal">Δ</mi> <msup> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <mi>h</mi> <mrow> <mo>-</mo> <mn>1</mn> <mo stretchy="false">/</mo> <mn>2</mn> </mrow> </msup> <mrow> <mo stretchy="false">(</mo> <mi mathvariant="normal">Δ</mi> <mi>t</mi> <mo>+</mo> <mi>h</mi> <mo stretchy="false">)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> </math></EquationSource> </InlineEquation>, we need to ensure that <InlineEquation ID="IEq5"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_2821_Article_IEq5.gif" Format="GIF" Height="21" Rendition="HTML" Resolution="72" Type="Linedraw" Width="155" /> </InlineMediaObject> <EquationSource Format="TEX">\(h^{-1/2}(1 + \Delta t^{-1}h) \lesssim 1\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msup> <mi>h</mi> <mrow> <mo>-</mo> <mn>1</mn> <mo stretchy="false">/</mo> <mn>2</mn> </mrow> </msup> <mrow> <mo stretchy="false">(</mo> <mn>1</mn> <mo>+</mo> <mi mathvariant="normal">Δ</mi> <msup> <mi>t</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>h</mi> <mo stretchy="false">)</mo> </mrow> <mo>≲</mo> <mn>1</mn> </mrow> </math></EquationSource> </InlineEquation>. This condition is not feasible unless <InlineEquation ID="IEq6"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10915_2025_2821_Article_IEq6.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="42" /> </InlineMediaObject> <EquationSource Format="TEX">\(h &gt; 1\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>h</mi> <mo>&gt;</mo> <mn>1</mn> </mrow> </math></EquationSource> </InlineEquation> is required. Within the scope of this research, we demonstrate that the error analytical technique established for the isotropic case, is also helpful for the anisotropic cases. Finally, we showcase a variety of numerical experiments, including the convergence tests to validate the theoretical analysis, as well as numerical simulations to to further support our findings.</p>

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

Error Analysis of Finite Element Approximation for Mean Curvature Flows in Axisymmetric Geometry

  • Meng Li

摘要

We focus on establishing a novel convergence analysis of the numerical approximation methods applied to axisymmetric mean curvature flows (MCFs) with genus-1 surfaces, including both isotropic and anisotropic cases. This advancement builds upon and enhances the findings presented in a previous work in Barrett et al. (IMA J Numer Anal 41:1641–1667, 2021), which exists a strict restriction on the time-space step ratio, i.e. \(\Delta t \lesssim \sqrt{h}\) Δ t h . In this work, we first adopt a novel time-space error splitting technique to remove above restriction under the isotropic case. Furthermore, we discover that this splitting also aids in addressing another unresolved issue present in anisotropic MCFs. Indeed, almost any restrictions on the time and space steps cannot derive the \(L^\infty \) L -boundedness of \(D_t\vec {X}_{h, \rho }^m\) D t X h , ρ m , which is required for convergence analysis in the anisotropic case. Given that \(\Vert D_t\vec {X}_{h, \rho }^m\Vert _{L^\infty } \lesssim \Delta t^{-1}h^{-1/2}(\Delta t + h) + 1\) D t X h , ρ m L Δ t - 1 h - 1 / 2 ( Δ t + h ) + 1 , we need to ensure that \(h^{-1/2}(1 + \Delta t^{-1}h) \lesssim 1\) h - 1 / 2 ( 1 + Δ t - 1 h ) 1 . This condition is not feasible unless \(h > 1\) h > 1 is required. Within the scope of this research, we demonstrate that the error analytical technique established for the isotropic case, is also helpful for the anisotropic cases. Finally, we showcase a variety of numerical experiments, including the convergence tests to validate the theoretical analysis, as well as numerical simulations to to further support our findings.