<p>In this paper, we study the time-fractional Cahn-Hilliard (TFCH) equation, which models phase separation processes along with nonlocal memory effects. A fully-discrete numerical scheme is developed using the ultra-weak discontinuous Galerkin (UWDG) method in space with generalized numerical fluxes, coupled with a non-uniform L1 time-stepping scheme. The proposed method is demonstrated to preserve key physical properties of the continuous model, including mass conservation and energy dissipation. We establish the unconditional unique solvability of the fully-discrete scheme using the convex-concave splitting approach and derive stability bounds for the order parameter. An optimal a priori error estimate in the <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(L^2\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>L</mi> <mn>2</mn> </msup> </math></EquationSource> </InlineEquation>-norm is established, confirming the convergence of the proposed scheme. Finally, numerical experiments are performed to validate the theoretically predicted convergence order for different choices of numerical fluxes. Numerical simulations confirm the scheme’s accuracy, demonstrating energy dissipation, mass conservation, and phase separation phenomena.</p>

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

Analysis of an Ultra-Weak Discontinuous Galerkin Method with Generalized Numerical Fluxes for Time-Fractional Cahn-Hilliard Equation

  • Deeksha Singh,
  • Monirul Islam,
  • Rajen Kumar Sinha

摘要

In this paper, we study the time-fractional Cahn-Hilliard (TFCH) equation, which models phase separation processes along with nonlocal memory effects. A fully-discrete numerical scheme is developed using the ultra-weak discontinuous Galerkin (UWDG) method in space with generalized numerical fluxes, coupled with a non-uniform L1 time-stepping scheme. The proposed method is demonstrated to preserve key physical properties of the continuous model, including mass conservation and energy dissipation. We establish the unconditional unique solvability of the fully-discrete scheme using the convex-concave splitting approach and derive stability bounds for the order parameter. An optimal a priori error estimate in the \(L^2\) L 2 -norm is established, confirming the convergence of the proposed scheme. Finally, numerical experiments are performed to validate the theoretically predicted convergence order for different choices of numerical fluxes. Numerical simulations confirm the scheme’s accuracy, demonstrating energy dissipation, mass conservation, and phase separation phenomena.