ترغب بنشر مسار تعليمي؟ اضغط هنا

Highly efficient and energy dissipative schemes for the time fractional Allen-Cahn equation

138   0   0.0 ( 0 )
 نشر من قبل Dianming Hou
 تاريخ النشر 2021
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English




اسأل ChatGPT حول البحث

In this paper, we propose and analyze a time-stepping method for the time fractional Allen-Cahn equation. The key property of the proposed method is its unconditional stability for general meshes, including the graded mesh commonly used for this type of equations. The unconditional stability is proved through establishing a discrete nonlocal free energy dispassion law, which is also true for the continuous problem. The main idea used in the analysis is to split the time fractional derivative into two parts: a local part and a history part, which are discretized by the well known L1, L1-CN, and $L1^{+}$-CN schemes. Then an extended auxiliary variable approach is used to deal with the nonlinear and history term. The main contributions of the paper are: First, it is found that the time fractional Allen-Chan equation is a dissipative system related to a nonlocal free energy. Second, we construct efficient time stepping schemes satisfying the same dissipation law at the discrete level. In particular, we prove that the proposed schemes are unconditionally stable for quite general meshes. Finally, the efficiency of the proposed method is verified by a series of numerical experiments.



قيم البحث

اقرأ أيضاً

A time-fractional Allen-Cahn equation with volume constraint is first proposed by introducing a nonlocal time-dependent Lagrange multiplier. Adaptive linear second-order energy stable schemes are developed for the proposed model by combining invarian t energy quadratization and scalar auxiliary variable approaches with the recent L1$^{+}$ formula. The new developed methods are proved to be volume-preserving and unconditionally energy stable on arbitrary nonuniform time meshes. The accelerated algorithm and adaptive time strategy are employed in numerical implement. Numerical results show that the proposed algorithms are computationally efficient in multi-scale simulations, and appropriate for accurately resolving the intrinsically initial singularity of solution and for efficiently capturing the fast dynamics away initial time.
In this paper, we study diagonal dominance of the stiffness matrix resulted from the piecewise linear finite element discretisation of the integral fractional Laplacian under global homogeneous Dirichlet boundary condition in one spatial dimension. W e first derive the exact form of this matrix in the frequency space which is extendable to multi-dimensional rectangular elements. Then we give the complete answer when the stiffness matrix can be strictly diagonally dominant. As one application, we apply this notion to the construction of maximum principle preserving schemes for the fractional-in-space Allen-Cahn equation, and provide ample numerical results to verify our findings.
We present several first-order and second-order numerical schemes for the Cahn-Hilliard equation with discrete unconditional energy stability. These schemes stem from the generalized Positive Auxiliary Variable (gPAV) idea, and require only the solut ion of linear algebraic systems with a constant coefficient matrix. More importantly, the computational complexity (operation count per time step) of these schemes is approximately a half of those of the gPAV and the scalar auxiliary variable (SAV) methods in previous works. We investigate the stability properties of the proposed schemes to establish stability bounds for the field function and the auxiliary variable, and also provide their error analyses. Numerical experiments are presented to verify the theoretical analyses and also demonstrate the stability of the schemes at large time step sizes.
We develop and analyze a class of maximum bound preserving schemes for approximately solving Allen--Cahn equations. We apply a $k$th-order single-step scheme in time (where the nonlinear term is linearized by multi-step extrapolation), and a lumped m ass finite element method in space with piecewise $r$th-order polynomials and Gauss--Lobatto quadrature. At each time level, a cut-off post-processing is proposed to eliminate extra values violating the maximum bound principle at the finite element nodal points. As a result, the numerical solution satisfies the maximum bound principle (at all nodal points), and the optimal error bound $O(tau^k+h^{r+1})$ is theoretically proved for a certain class of schemes. These time stepping schemes under consideration includes algebraically stable collocation-type methods, which could be arbitrarily high-order in both space and time. Moreover, combining the cut-off strategy with the scalar auxiliary value (SAV) technique, we develop a class of energy-stable and maximum bound preserving schemes, which is arbitrarily high-order in time. Numerical results are provided to illustrate the accuracy of the proposed method.
This article is mainly devoted to the asymptotic analysis of a fractional version of the (elliptic) Allen-Cahn equation in a bounded domain $Omegasubsetmathbb{R}^n$, with or without a source term in the right hand side of the equation (commonly calle d chemical potential). Compare to the usual Allen-Cahn equation, the Laplace operator is here replaced by the fractional Laplacian $(-Delta)^s$ with $sin(0,1/2)$, as defined in Fourier space. In the singular limit $varepsilonto 0$, we show that arbitrary solutions with uniformly bounded energy converge both in the energetic and geometric sense to surfaces of prescribed nonlocal mean curvature in $Omega$ whenever the chemical potential remains bounded in suitable Sobolev spaces. With no chemical potential, the notion of surface of prescribed nonlocal mean curvature reduces to the stationary version of the nonlocal minimal surfaces introduced by L.A. Caffarelli, J.M. Roquejoffre, and O. Savin. Under the same Sobolev regularity assumption on the chemical potential, we also prove that surfaces of prescribed nonlocal mean curvature have a Minkowski codimension equal to one, and that the associated sets have a locally finite fractional $2s^prime$-perimeter in $Omega$ for every $s^primein(0,1/2)$.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا