Processing math: 58%
Research article Special Issues

On an efficient numerical procedure for the Functionalized Cahn-Hilliard equation

  • The Functionalized Cahn Hilliard (FCH) equation was used to model micro-phase separation in mixtures of amphiphilic molecules in solvent. In this paper, we proposed a Tri-Harmonic Modified (THM) numerical approach for efficiently solving the FCH equation with symmetric double well potential by extending the ideas of the Bi-harmonic Modified (BHM) method. THM formulation allowed for the nonlinear terms in the FCH equation to be computed explicitly, leading to fast evaluations at every time step. We investigated the convergence properties of the new approach by using benchmark problems for phase-field models, and we directly compared the performance of the THM method with the recently developed scalar auxiliary variable (SAV) schemes for the FCH equation. The THM modified scheme was able to produce smaller errors than those obtained from the SAV formulation. In addition to this direct comparison with the SAV schemes, we tested the adaptability of our scheme by using an extrapolation technique which allows for errors to be reduced for longer simulation runs. We also investigated the adaptability of the THM method to other 6th order partial differential equations (PDEs) by considering a more complex form of the FCH equation with nonsymmetric double well potential. Finally, we also couple the THM scheme with a higher order time-stepping method, (implicit-explicit) IMEX schemes, to demonstrate the robustness and adaptability of the new scheme. Numerical experiments are presented to investigate the performance of the new approach.

    Citation: Saulo Orizaga, Ogochukwu Ifeacho, Sampson Owusu. On an efficient numerical procedure for the Functionalized Cahn-Hilliard equation[J]. AIMS Mathematics, 2024, 9(8): 20773-20792. doi: 10.3934/math.20241010

    Related Papers:

    [1] Monica-Felicia Bota, Liliana Guran . Existence of a solution of fractional differential equations using the fixed point technique in extended b-metric spaces. AIMS Mathematics, 2022, 7(1): 518-535. doi: 10.3934/math.2022033
    [2] Aftab Hussain . Fractional convex type contraction with solution of fractional differential equation. AIMS Mathematics, 2020, 5(5): 5364-5380. doi: 10.3934/math.2020344
    [3] Muhammad Sarwar, Aiman Mukheimer, Syed Khayyam Shah, Arshad Khan . Existence of solutions of fractal fractional partial differential equations through different contractions. AIMS Mathematics, 2024, 9(5): 12399-12411. doi: 10.3934/math.2024606
    [4] Abdullah Eqal Al-Mazrooei, Jamshaid Ahmad . Fixed point approach to solve nonlinear fractional differential equations in orthogonal F-metric spaces. AIMS Mathematics, 2023, 8(3): 5080-5098. doi: 10.3934/math.2023255
    [5] Rajagopalan Ramaswamy, Gunaseelan Mani . Application of fixed point result to solve integral equation in the setting of graphical Branciari -metric spaces. AIMS Mathematics, 2024, 9(11): 32945-32961. doi: 10.3934/math.20241576
    [6] Müzeyyen Sangurlu Sezen . Interpolative best proximity point results via γ-contraction with applications. AIMS Mathematics, 2025, 10(1): 1350-1366. doi: 10.3934/math.2025062
    [7] Mohammed Shehu Shagari, Akbar Azam . Integral type contractions of soft set-valued maps with application to neutral differential equations. AIMS Mathematics, 2020, 5(1): 342-358. doi: 10.3934/math.2020023
    [8] Gonca Durmaz Güngör, Ishak Altun . Fixed point results for almost (ζθρ)-contractions on quasi metric spaces and an application. AIMS Mathematics, 2024, 9(1): 763-774. doi: 10.3934/math.2024039
    [9] Ahmed Morsy, Kottakkaran Sooppy Nisar, Chokkalingam Ravichandran, Chandran Anusha . Sequential fractional order Neutral functional Integro differential equations on time scales with Caputo fractional operator over Banach spaces. AIMS Mathematics, 2023, 8(3): 5934-5949. doi: 10.3934/math.2023299
    [10] Mohammed Shehu Shagari, Qiu-Hong Shi, Saima Rashid, Usamot Idayat Foluke, Khadijah M. Abualnaja . Fixed points of nonlinear contractions with applications. AIMS Mathematics, 2021, 6(9): 9378-9396. doi: 10.3934/math.2021545
  • The Functionalized Cahn Hilliard (FCH) equation was used to model micro-phase separation in mixtures of amphiphilic molecules in solvent. In this paper, we proposed a Tri-Harmonic Modified (THM) numerical approach for efficiently solving the FCH equation with symmetric double well potential by extending the ideas of the Bi-harmonic Modified (BHM) method. THM formulation allowed for the nonlinear terms in the FCH equation to be computed explicitly, leading to fast evaluations at every time step. We investigated the convergence properties of the new approach by using benchmark problems for phase-field models, and we directly compared the performance of the THM method with the recently developed scalar auxiliary variable (SAV) schemes for the FCH equation. The THM modified scheme was able to produce smaller errors than those obtained from the SAV formulation. In addition to this direct comparison with the SAV schemes, we tested the adaptability of our scheme by using an extrapolation technique which allows for errors to be reduced for longer simulation runs. We also investigated the adaptability of the THM method to other 6th order partial differential equations (PDEs) by considering a more complex form of the FCH equation with nonsymmetric double well potential. Finally, we also couple the THM scheme with a higher order time-stepping method, (implicit-explicit) IMEX schemes, to demonstrate the robustness and adaptability of the new scheme. Numerical experiments are presented to investigate the performance of the new approach.


    One of the concepts which have significant impact in Clifford analysis is studying the approximation of a Cliffordian function h(m)(x), xRm+1 as a series of the form:

    n=0Q(m)n(x) a(m)n,  a(m)nCm, (1.1)

    where {Q(m)n(x), xRm+1} is a prescribed base of Cllifordian polynomials and Cm is the real Clifford algebra of dimension 2m. In 1990, this problem was addressed by the authors of [1]. The polynomials are assumed to form a Hamel basis of Sm(x) (the Clifford linear space of all special monogenic polynomials (SMPs) with Clifford algebra coefficients). The series (1.1) is called the basic series associated with the base of SMPs. Many results about the approximation of SMFs and ESMFs by basic series, which can be associated with bases of SMPs [2,3].

    The theory of basic series in the case of one complex variable was originally discovered by Whittaker and Cannon [4,5,6,7] about 90 years ago. As we have mentioned earlier, the attempt done by authors of [1] were the first to extend the notion of basic series in the case of Clifford analysis.

    In the case of a single complex variable, the approximation properties of the derivative and integral bases of a certain base of polynomials of in a disk of center origin have been studied by many authors, of whom we may mention Makar [8], Mikhail [9], and Newns [10]. In the case of of several complex variables the domains of representation are hyperspherical, hyperelliptical and polycylinderical regions (see [11,12]). Afterwards, the authors of [13,14] generalized this problem in Clifford analysis, which is called hypercomplex primitive and derivative bases of SMPs and their representations is in closed hyperballs.

    Approximation theory is a rich topic which has numerous applications in various scientific disciplines such that mathematical analysis, statistics, engineering and physics. Recently, order moment of the wind power time series has been studied in [15]. Although our study is narrowed to theoretical aspects, the basic sets (bases) of polynomials proved its efficiency in as solutions to important partial differential equations, such as the heat equation [16] and wave equation [17,18].

    The authors of [19] introduced an expansion of a SMF by basic series of generalized Bessel SMPs. They proved that the GBSMPs are solutions of second order homogeneous differential equations. Furthermore, in [3], the authors proved an extended version of Hadamard's three-hyperballs theorem to study the overconvergence properties. One of the recent fascinating research findings can be found in [20] where the authors of used the Hadamard's three-hyperballs theorem to generalize the Whittaker-Cannon theorem in open hyperballs in Rm+1. Precisely, they proved that the hypercomplex Cannon functions preserved the effectiveness properties of both Cannon and non-Cannon bases. In the very recent paper [21] the authors derived a new base of SMPs in F-modules, named the equivalent base. They have also studied the convergence properties (effectiveness, order and type, Tρ-property) of these base.

    In 2017, a study based primarily on combination of Clifford analysis and functional analysis [26] when the considered bases {Q(m)n(x)} are not necessarily consisting of polynomials. The convergence properties of these general bases had been studied in F-modules. Precisely, a general criterion for effectiveness of basic series in F-modules was constructed.

    Recently in [22], the authors have studied a new base called hypercomplex Ruscheweyh derivative bases (HRDBs). They investigated the effectiveness properties of HRDBs of a given base of SMPs in different regions of convergence in F-modules. The above treatment is considered to extend and improve the results in Clifford and complex given in [8,9,10,13,14].

    Motivated by the preceding discussion, the current work introduces a modified generalization of the Hasse derivative operator (HDO). Acting by hypercomplex HDO on bases, we derive a base of SMPs, which we may call the hypercomplex Hasse derivative bases of SMPs (HHDBSMPs). Consequently, we discuss the effectiveness properties, mode of increase, and the Tρ-property of such a base in several regions: closed and open hyperballs, open regions surrounding closed hyperball, at the origin, and for all entire SMFs. Some applications on the HHD of Bernoulli SMPs (BSMPs), Euler SMPs (ESMPs), proper Bessel SMPs (PBSMPs), general Bessel SMPs (GBSMPs) and Chebyshev SMPs (CSMPs) are also provided. The obtained results offer new generalizations of existing work concerning the convergence properties of polynomials bases in both complex and Clifford settings.

    This section collects some notations and basic results which are needed throughout the paper. More details can be found in the literature, see [1,23,24,25,26]. The real Clifford algebra over R is defined as

    Cm={b=B{1,...m} bBeB, bBR},

    where ei=e{i}, i=1,,m, e0=eϕ=1 and eB=eβ1...eβh, with 1β1<β2<<βhm. The product in Cm is determined by the relations eiej+ejei=2δij where δij denotes the Kronecker delta and e0=1 for 1ijm (for details on the main concepts about Cm, see [25]). The norm of a Clifford number is given by |b|=(BN|bB|2)12 where N stands for {1,,m}.

    Since Cm is isomorphic to R2m we may provide it with the R2m-norm |b| and one sees easily that for any b,cCm, |bc|2m2 |b||c|.

    The elements (x0,x)=(x0,x1,...,xm)Rm+1 will be identified with the Clifford numbers x0+x_=x0+mj=1ejxj. Note that if x=x0+x_Rm+1, ¯x=x0x_.

    Definition 2.1. Let xRm+1 and ΩRm+1 be an open set, then the function h(m)(x) is called left monogenic in Ω if D[h(m)(x)]=0 where

    D=mi=0eixi

    is the generalized Cauchy-Riemann operator. Similarly, h(m)(x) is a right monogenic function if [h(m)(x)]D=0.

    Definition 2.2. A polynomial Q(m)(x) is SMP iff DQ(m)(x)=0 and has the form:

    Q(m)(x)=finitei,j¯xixja(m)i,j,

    where a(m)i,jCm.

    Let Sm[x] be the space of all SMPs is the right Cm-module defined by

    Sm[x]=spanCm{Q(m)n(x): nN},

    where Q(m)n(x) was given in [1] as follows:

    Q(m)n(x)=n!(m)nr+s=n(m12)r(m+12)sr!s!¯xrxs, (2.1)

    where for βR, (β)l=β(β+1)(β+l1) is the Pochhamer symbol.

    Definition 2.3. Let ΩRm+1 be a connected open containing 0 and h(m) is monogenic in Ω. The function h(m) is said to be SMF in Ω if and only if its Taylor series near zero exists and cab be expressed as: h(m)(x)=n=0Q(m)n(x) a(m)n for some SMPs Q(m)n(x).

    If Q(m)n(x) is a homogeneous SMP has degree n in x, (see [1])

    Q(m)n(x)=Q(m)n(x) β(m)n,

    where β(m)nCm is a constant. Accordingly, it follows that:

    Q(m)nR=sup¯B(R)|Q(m)n(x)|=Rn.

    Next, we recall the definition of F-module.

    Definition 2.4. An F-module E over Cm is a complete Hausdorff topological vector space by countable family of a proper system of semi-norms Q={.s}s0 such that s<th(m)sh(m)t; (h(m)E), Hence WE is open iff h(m)W, ϵ>0, M0 such that {g(m)E:h(m)g(m)s)ϵ}W,sM.

    Definition 2.5. A sequence {h(m)n} in an F-module E converges to g(m) in E if

    limn h(m)ng(m)s=0

    for all .sQ.

    The domains of representation adopted here are the open hyperball B(R), the closed hyperball ¯B(R) and B+(R); R>0, where B+(R) any open hyperball enclosing closed hyperball, these are the sets defined by

    B(R)={xRm+1:|x|<R},
    ¯B(R)={xRm+1:|x|R},
    B+(R)={xRm+1:|x|<R+}.

    Table 1 summarizes certain classes of SMFs which represent F-modules where xRm+1 and each space is associated with the a proper countable system of semi-norms as follows.

    Table 1.  F-modules.
    Space The Associated Semi-Norms
    M[B(R)]: Class of SMFs in B(R) h(m)r=sup¯B(r)|h(m)(x)|,  r<R, h(m)M[B(R)],
    M[¯B(R)]: Classe of SMFs in ¯B(R) h(m)R=sup¯B(R)|h(m)(x)|,   h(m)M[¯B(R)],
    M[B+(R)]: Class of SMFs in B+(R) h(m)r=sup¯B(r)|h(m)(x)|,  R<r, h(m)M[B+(R)],
    M[0+]: Class of SMFs at the origin h(m)ϵ=sup¯B(ϵ)|h(m)(x)|, ϵ>0  h(m)M[0+] ,
    M[]: Class of ESMFs on Rm+1 h(m)n=sup¯B(n)|h(m)(x)|, n<  h(m)M[].

     | Show Table
    DownLoad: CSV

    Now, let {Q(m)n(x)} be a base of an F-module E such that

    Q(m)n(x)=k=0Q(m)k(x) Q(m)n,k,  Q(m)n,kCm, (2.2)
    Q(m)n(x)=k=0 Q(m)k(x) π(m)n,k,  π(m)n,kCm, (2.3)
    Q(m)nR=sup¯B(R)|Q(m)n(x)|, (2.4)
    ΨQ(m)n(R)=kQ(m)k π(m)n,kR, (2.5)

    this sum is called hypercomplex Cannon sum, where

    Q(m)k π(m)n,kR=sup¯B(R)|Q(m)k(x) π(m)n,k|,ΨQ(m)(R)=lim supn{ΨQ(m)n(R)}1n, (2.6)

    where ΨQ(m)(R) is called the hypercomplex Cannon function of the base {Q(m)n(x)} in closed hyperball ¯B(R).

    Let Dn is the degree of the polynomial of highest degree in the representation (2.3) the following restrictions are imposed.

    limn{Dn}1n=1, (2.7)
    Dn=O[na],  a1,  (2.8)
    Dn=o(nlogn). (2.9)

    If dk is the degree of the polynomials {Q(m)k(x)}, then dkDn for all kn (see [1]).

    If Q(m)=(Q(m)n,k) and Π(m)=(π(m)n,k) are the Clifford matrices of coefficients and operators respectively of the set {Q(m)n(x)}. Thus according to [1] the set {Q(m)n(x)} will be base iff

    Q(m)Π(m)=Π(m)Q(m)=I, (2.10)

    where I is the unit matrix.

    Let h(m)(x)=n=0Q(m)n(x) an(h(m)) be any function which is SMF at the origin, substituting for Q(m)n(x) from (2.3) we obtain the basic series

    h(m)(x)n=0Q(m)n(x) Πn(h(m)), (2.11)

    where

    Πn(h(m))=k=0π(m)k,n ak(h(m)) . (2.12)

    The authors in [22,26] introduced the idea of effectiveness for the class M[¯B(R)]. A base {Q(m)n(x)} is effective for the class M[¯B(R)] If the basic series (2.11) converges normally to every function h(m)(x)M[¯B(R)] which is SMF in ¯B(R). Similar definitions are used for the classes M[B(R)], M[B+(R)], M[] and M[0+].

    They also proved:

    Theorem 2.1. A base {Q(m)n(x)} is effective for the classes M[¯B(R)], M[B(R)], M[B+(R)], M[] or M[0+] if and only if ΨQ(m)(R)=R, ΨQ(m)(r)<R   r<R, ΨQ(m)(R+)=R, ΨQ(m)(R)<  R<, or ΨQ(m)(0+)=0.

    For the definition of bases of SMPs and theorems governing the effectiveness properties of bases of SMPs, the reader is referred to the authors [21,22,26].

    The complex Hasse derivative operator (CHDO) of order i is defined in [28,29,30]. Using the definition of the complex Hasse derivative, we can define a new operator in the case of Clifford setting called the hypercomplex Hasse derivative (HHD) as follows:

    Definition 3.1. For each integer i0 the HHD H(i) of order i is defined by

    H(i)(Q(m)n(x))=ζn,i Q(m)ni(x), (3.1)

    where

    ζn,i=nii!i1j=1(1jn)

    and H(i) is closely related to the higher hypercomplex derivative (12¯D)i: H(i)=1i!(12¯D)i.

    The set {Q(m)n(x)} is an Appell sequence with respect to x0 or 12¯D : 12¯DQ(m)n(x)=nQ(m)n1(x).

    Remark 3.1. If xC1 then (3.1) is reduced to the ordinary Hasse derivative of order i (see [28,29,30]),

    Definition 3.2. Let {Q(m)n(x)} be a base. By acting on both sides of Eq (2.2) with the operator H(i), we get

    H(i)Q(m)n(x)=k  ζk,i Q(m)ki(x)Q(m)n,k. (3.2)

    The set {H(i)Q(m)n(x)} = {H(i,m)(x)} is defined the Hypercomplex Hasse derivative bases (HHDBs).

    The present work deals principally with the convergence properties of certain classes of bases, namely HHDBs. In fact we shall study the convergence of the expansion of certain classes of functions as series of HHDBs. This study will be based on the already established theorems dealing with the convergence of basic series of HHDBs. The convergence properties of HHDBs are mainly classified as follows:

    (1) The region of effectiveness of HHDBs for the classes M[B(R)], M[¯B(R)], M[B+(R)], M[0+], and M[].

    (2) The mode of increase of HHDBs which determined by the order and type.

    (3) The Tρ-property of HHDBs.

    In the following sections, we will investigated all of these problems.

    In the current section, the property of effectiveness concerning the HHDBs in several regions such as M[B(R)], M[0+], M[] and M[B+(R)] are demonstrated.

    Theorem 4.1. If {Q(m)n(x)} is a base, then the HHD set {H(i,m)(x)} is also base.

    Proof. We form the coefficient matrix H(i,m) by defining the HHDBs in (2.2)

    H(i,m)n(x)=k Q(m)ki(x) ζk,iQ(m)n,k.

    Hence, the coefficients matrix H(i,m) is given by the following:

    H(i,m)=(H(i,m)n,k)=(ζk,i Q(m)n,k).

    Also, the operators matrix Π(i,m) follows from the effect H(i) on both sides of the representation (2.3) where

    Q(m)ki(x)=1ζn,ikπ(m)n,kH(i,m)k(x),

    and

    Π(i,m)=(π(i,m)n,k)=(1ζn,i π(m)n,k).

    Consequently,

    H(i,m)Π(i,m)=(kH(i,m)n,kπ(i,m)k,h)=(kQ(m)n,kπ(m)k,h)=(δn,h)=I.

    Moreover,

    Π(i,m)H(i,m)=(kπ(i,m)n,kH(i,m)k,h)=(k1ζn,iπ(m)n,k ζh,i Q(m)k,h)=(ζh,iζn,iδn,h)=I.

    We easily obtain from (2.10) that the set {H(i,m)n(x)} is a base.

    Theorem 4.2. The base {Q(m)n(x)} and its HHDBs {H(i,m)n(x)} have the same region of effectiveness for the class M[B(R)].

    Proof. If Q(m)n(x) is a base, Q(m)nr=sup¯B(r)|Q(m)n(x)| and H(i,m)nr=sup¯B(r)|H(i,m)n(x)|, then

    H(i,m)nr=sup¯B(r)|H(i,m)n(x)|=sup¯B(r)|jζj,i Q(m)ji(x)Qn,j|2m/2jζj,i rjiQ(m)nRRj=2m/2riΥ(r,R) Q(m)nR=K1 Q(m)nR    for all r<R, (4.1)

    where K1=2m/2riΥ(r,R) and Υ(r,R)=j=0ζj,i(rR)j<.

    Using (2.5) and (4.1), it follows that the hypercomplex Cannon sum of the HHDBs {H(i,m)n(x)} is given by

    ΨH(i,m)n(r)=kH(i,m)kπ(i,m)n,krK1kQ(m)k π(i,m)n,kR=K1ζn,ikQ(m)kπ(m)n,kR=K1ζn,iΨQ(m)n(R). (4.2)

    Using (2.6) and (4.2), we obtain that the hypercomplex Cannon function of the HHDBs is given by:

    ΨH(i,m)(r)ΨQ(m)(R),   r<R. (4.3)

    Now, suppose that the base {Q(m)n(x)} is effective for M[B(R)], we can apply Theorem 2.1, we have

    ΨQ(m)(r)<R,   r<R. (4.4)

    Hence there is a number r1 such that r<r1<R, then from (4.3) and (4.4), we deduce that

    ΨH(i,m)(r)ΨQ(m)(r1)<R,   r<R,

    that is to say the base {H(i,m)n(x)} is effective for M[B(R)].

    Theorem 4.3. The base {Q(m)n(x)} and its HHDBs {H(i,m)n(x)} have the same region of effectiveness for the class M[0+] or M[].

    Proof. Suppose that the base {Q(m)n(x)} is effective for M[0+], we can apply Theorem 2.1, it follows that ΨQ(m)(0+)=0. Making R,r0+ in (4.3), we have ΨH(i,m)(0+)ΨQ(m)(0+)=0 but we know that ΨH(i,m)(0+)0, thus, ΨH(i,m)(0+)=0. Therefore, the base {H(i,m)n(x)} is effective for M[0+].

    Now, suppose that the base {Q(m)n(x)} is effective for M[]. Applying Theorem 2.1 we conclude that

    ΨQ(m)(r)<, r<. (4.5)

    Thus if we choose the number r2 such that r<r2<R, making R in (4.3). Then, by using (4.5), we obtain that

    ΨH(i,m)(r)ΨQ(m)(r2)<, r<,

    and, the base {H(i,m)n(x)} will be effective for M[].

    Theorem 4.4. The base {Q(m)n(x)} and its HHDBs {H(i,m)n(x)} have the same region of effectiveness for the class M[B+(R)].

    Proof. If the base {Q(m)n(x)} is effective for M[B+(r3)] and r3 is any positive number such that r3<r, we can apply Theorem 2.1, we obtain

    ΨQ(m)(r+3)=r3,  r3<r<R. (4.6)

    Making Rr+3 in (4.3), we easily obtain, from (4.6) that ΨH(i,m)(r+3)ΨQ(m)(r+3)=r3, but ΨH(i,m)(r+3)r3 which implies that ΨH(i,m)(r+3)=r3. Hence, the base {H(i,m)n(x)} is indeed effective for M[B+(r3)] as required.

    When the representation (2.3) is finite then the base is called SMPs. In this section we will discuss the region of effectiveness of HHDBSMPs for the class of SMFs in ¯B(R). The following result states the purpose of this section.

    Theorem 5.1. The base {Q(m)n(x)} for which the condition (2.7) is satisfied and its HHDBSMPs {H(i,m)n(x)} have the same region of effectiveness for the class M[¯B(R)].

    Proof. If Q(m)n(x) is a base of SMPs, Q(m)nR=sup¯B(R)|Q(m)n(x)| and H(i,m)nR=sup¯B(R)|H(i,m)n(x)|, then

    H(i,m)nR=sup¯B(R)|H(i,m)n(x)|=sup¯B(R)|jQ(m)n,j ζj,i Q(m)ji(x)|jQ(m)nRRj ζj,i Rji=Q(m)nRRij ζj,iQ(m)nRRi ζdn,i (ζdn,i+1), (5.1)

    where dn is the degree of the polynomial Q(m)n(x), dnDn. Applying (2.5) and (5.1), it follows that

    ΨH(i,m)n(R)=k H(i,m)k π(i,m)n,kR1ζn,iRikQ(m)k π(m)n,kR ζdk,i (ζdk,i+1)1ζn,iRiζDn,i (ζDn,i+1) ΨQ(m)n(R). (5.2)

    A combination of (2.6), (2.7) and (5.2), gives ΨH(i,m)n(R)ΨP(m)(R)R.

    But ΨH(i,m)n(R)R. We finally deduce that

    ΨH(i,m)n(R)=R. (5.3)

    and the HHDBSMPs {H(i,m)n(x)} is effective for M[¯B(R)].

    The following example shows that the condition (2.7) imposed on the class of the base {Q(m)n(x)} cannot be relaxed.

    Example 5.1. Theorem 5.1 is not always correct if the condition (2.7) is not satisfied. Let

    Q(m)n(x)={Q(m)n(x),n is even,Q(m)n(x)+ Q(m)b(x), b=2n,n is odd.

    When n is even, we have Q(m)n(x)=Q(m)n(x) and hence ΨQ(m)n(R)=Rn. Thus, by taking R=1, then ΨQ(m)n(1)=1, and limn{ΨQ(m)2n(1)}12n=1.

    Furthermore, Q(m)n(x)=Q(m)n(x)Q(m)b(x), when n is odd, then

    ΨQ(m)n(R)=Rn+2 Rb.

    So that when R=1, ΨQ(m)n(1)=3, we get

    limn {ΨQ(m)2n+1(1)}12n+1=1.

    Consequently, ΨQ(m)(1)=lim supn {ΨQ(m)n(1)}1n=1, and the base {Q(m)n(x)} is effective for M[¯B(1)].

    Forming the HHDBSMPs {H(i,m)n(x)}, we easily get

    H(i,m)n(x)={ζn,i Q(m)ni(x),n is even, and 2,ζn,i Q(m)ni(x)+ζb,iQ(m)bi(x), n is odd.

    Since Q(m)ni(x)=(1ζn,i) H(i,m)n(x), when n is even, then ΨH(i,m)n(R)=Rni, taking R=1, ΨH(i,m)n(1)=1. Hence,

    limn {ΨH(i,m)2n(1)}12n=1.

    When n is odd, Q(m)ni(x)=(1ζn,i)[H(i,m)Pn(x)ζb,iH(i,m)b(x)]. Hence we have, ΨH(i,m)n(R)=(1ζn,i)[ζn,iRni+2 ζb,iRbi].

    Taking R=1, then we get

    ΨH(i,m)(1)=lim supn {ΨH(i,m)2n+1(1)}12n+1=2>1,

    and the HHDBSMPs {H(i,m)n(x)} is not effective for M[¯B(1)].

    For a simple base of SMPs (Dn=n) (see [1]), we obtain the following corollary.

    Corollary 5.1. When the simple base {Q(m)n(x)} of SMPs is effective for M[¯B(R)], so also will be the HHDBSMPs {H(i,m)n(x)}.

    In [1,23], the idea of the order and type of the base {Q(m)n(x} of SMPs was introduced as follows:

    ρQ(m)=limRlim supnlog ΨQ(m)n(R)n log n. (6.1)

    and

    τQ(m)=limReρQ(m)lim supn{ΨQ(m)n(R)}1n ρQ(m)n. (6.2)

    Importantly, if the base {Q(m)n(x)} has finite order ρQ(m) and finite type τQ(m), then it can represent every ESMF of order less than 1ρQ(m) and type less than 1τQ(m) in any finite hyperball. Rich investigation on the order of certain classes of bases can be found in [31,32].

    Now, we explore the relation between the order and type of SMPs {Q(m)n(x)} and our constructed base; {H(i,m)n(x)} as follows.

    Theorem 6.1. Let ρQ(m) and τQ(m) be the order and type of the base of SMPs {Q(m)n(x)} satisfying the condition (2.8). Then the HHDBSMPs {H(i,m)n(x)} will be of order ρH(i,m)ρQ(m) and type τH(i,m)τQ(m) whenever ρH(i,m)=ρQ(m). The values of ρQ(m) and τQ(m) are attainable.

    Proof. The proof of this theorem denoted on the inequality (5.2), since

    ΨH(i,m)n(R)1ζn,iRαζDn,i (ζDn,i+1) ΨQ(m)n(R).

    Then

    limRlim supnlog ΨH(i,m)n(R)n log nlimRlim supnlog ζDn,i (ζDn,i+1)+log ΨQ(m)n(R)n log n.

    It follows, in view of (6.1), that the HHDBSMPs is at most ρQ(m).

    If ρH(i,m)=ρQ(m), we have

    limReρH(i,m)lim supn{Ψ(i)Hn(R)}1n(ρH(i,m))nlimReρQ(m)lim supn{ΨQn(R)}1n(ρQ(m))n,

    and the type of the HHDBSMPs is at most τQ(m).

    Note that the upper bound given in this theorem is attainable. We will illustrate this fact by introducing the following example:

    Example 6.1. Let {Q(m)n(x)} be the base of SMPs given by Q(m)n(x)=nn+Q(m)n(x), Q(m)0(x)=1, for which

    ΨQ(m)n(R)=nn[2+(Rn)n].

    It is easily seen that the base {Q(m)n(x)} is of order ρQ(m)=1 and type τQ(m)=e. Construct now the base {H(i,m)n(x)} such that

    H(i,m)n(x)=nn+ζn,i Q(m)ni(x),  Q(m)0(x)=1.

    Hence,

    ΨH(i,m)n(R)=nnζn,i[2+ζn,iRi(Rn)n].

    Therefore, the base {H(i,m)n(x)} is of order ρH(i,m)=1 and type τH(i,m)=e.

    The following example illustrates the best possibility of condition (2.8).

    Example 6.2. Let the base {Q(m)n(x)} of SMPs be defined by

    Q(m)n(x)={Q(m)n(x),n is even,Q(m)n(x)+μb2μ Q(m)2μ(x),n is odd and μ=nn,b>1.

    Hence,

    Q(m)n(x)=Q(m)n(x)μb2μQ(m)2μ(x),

    and

    ΨQ(m)n(R)=Rn+2μ(Rb)2μ.

    It is easy to see that the base Q(m)n(x) is of order ρQ(m)=1.

    For the HHDBSMPs {H(i,m)n(x)} it can verified that

    H(i,m)n(x)={ζn,i Q(m)ni(x),n is even,ζn,i Q(m)ni(x)+μb2μ ζ2μ,i Q(m)2μi(x),n is odd.

    Thus,

    \mathcal{Q}_{n-i}^{(m)}(x) = \frac{1}{\zeta_{n, i}}\mathbb{H}_n^{(i, m)}(x)-\frac{\mu}{b^{2\mu} }\ \frac{\zeta_{2\mu, i}}{\zeta_{n, i}} \ \mathbb{H}_{2\mu}^{(i, m)}(x).

    Consequently,

    \Psi_{\mathbb{H}_{n}^{(i, m)}}(R) = R^{n-i}+\frac{2\mu}{b^{i}} \ \frac{\zeta_{2\mu, i}}{\zeta_{n, i}} \ \left(\frac{R}{b}\right)^{2\mu-i}.

    Therefore, \rho_{\mathbb{H}^{(i, m)}} = 2 and \rho_{\mathbb{H}^{(i, m)}} > \rho_{Q^{(m)}}. This completes the proof.

    If the base of SMPs \{Q_{n}^{(m)}(x)\} is simple base ( D_{n} = n ) (see [1]), then the following corollary is a special case of Theorem 6.1.

    Corollary 6.1. When the simple base \{Q_{n}^{(m)}(x)\} of SMPs is of order \rho_{Q^{(m)}} and type \tau_{Q^{(m)}} , then the HHDBSMPs \{\mathbb{H}_n^{(i, m)}(x)\} will be of order \rho_{\mathbb{H}^{(i, m)}}\leq \rho_{Q^{(m)}} and type \tau_{\mathbb{H}^{(i, m)}}\leq \tau_{Q^{(m)}} whenever \rho_{\mathbb{H}^{(i, m)}} = \rho_{Q^{(m)}} .

    In the following, we determine the \mathbb{T}_{\rho_{Q^{(m)}}} -property of the HHDBs. The authors of [2] deduced \mathbb{T}_{\rho_{Q^{(m)}}} -property of the base \{Q^{(m)}_{n}(x)\} in Clifford analysis in open hyperball B(R) , closed hyperball \overline{B}(R) and at the origin are defined as follows:

    Definition 6.1 If the base \{Q^{(m)}_{n}(x)\} represents all ESMFs of order less than \rho_{Q^{(m)}} in \overline{B}(R) , B(R) or at the origin, then it is said to have property T_{\rho_{Q^{(m)}}} in \overline{B}(R) , B(R) or at the origin.

    Let

    \Psi_{Q^{(m)}}(R) = \limsup\limits_{n\rightarrow \infty}\ \frac{\log \ \Psi_{Q^{(m)}_{n}}(R)}{n \ \log \ n}.

    The following theorem concerning the property \mathbb{T}_{\rho_{Q^{(m)}}} of the base \{Q^{(m)}_{n}(x)\} (see [2]).

    Theorem 6.2. A base \{Q^{(m)}_{n}(x)\} to have the property T_{\rho_{Q^{(m)}}} for all ESMF of order less than \rho_{Q^{(m)}} in closed hyperball \overline{B}(R) , open hyperball B(R) or at the origin iff, \Psi_{Q^{(m)}}(R)\leq \frac{1}{\rho_{Q^{(m)}}} , \Psi_{Q^{(m)}}(r)\leq \frac{1}{\rho_{Q^{(m)}}} for all r < R or \Psi_{Q^{(m)}}(0^{+})\leq \frac{1}{\rho_{Q^{(m)}}}.

    Next, we construct the \mathbb{T}_{\rho_{\mathbb{H}^{(i, m)}}} -property of the HHDBSMPs in the closed hyperball \overline{B}(R) , for R > 0 .

    Theorem 6.3. Let \{Q_{n}^{(m)}(x)\} be the base of SMPs have \mathbb{T}\rho_{Q^{(m)}} -property in \overline{B}(R) , where R > 0 and for which the condition (2.9) is satisfied. Then the HHDBSMPs \{\mathbb{H}_n^{(i, m)}(x)\} have the same property.

    Proof. Suppose that the function \Psi_{\mathbb{H}^{(i, m)}}(R) given by:

    \begin{equation} \Psi_{\mathbb{H}^{(i, m)}}(R) = \limsup\limits_{n\rightarrow \infty}\ \frac{log \ \Psi_{\mathbb{H}^{(i, m)}_{n}}(R)}{n \ log \ n}, \end{equation} (6.3)

    where \Psi_{\mathbb{H}_{n}^{(i, m)}}(R) is the Cannon sum of the HHDBSMPs \{\mathbb{H}_n^{(i, m)}(x)\}. Then by using (2.9), (5.2) and (6.3), we obtain that

    \begin{equation} \Psi_{\mathbb{H}^{(i, m)}}(R)\leq \limsup\limits_{n\rightarrow \infty } \frac{\text{log}\zeta_{D_{n}, i} \ (\zeta_{D_{n}, i}+1)+ \text{log} \ \Psi_{Q_{n}^{(m)}}(R)}{n \ \text{log} \ n}\leq \Psi_{Q^{(m)}}(R) . \end{equation} (6.4)

    Since the base \{Q_{n}^{(m)}(x)\} has the property \mathbb{T}_{\rho_{Q^{(m)}}} in \overline{B}(R) , R > 0 . Hence by inequality (6.4) and Theorem 6.2, we have

    \Psi_{\mathbb{H}^{(i, m)}}(R)\leq \Psi_{Q^{(m)}}(R)\leq \frac{1}{\rho_{Q^{(m)}}},

    and the base \{\mathbb{H}_n^{(i, m)}(x)\} has the property \mathbb{T}_{\rho_{Q^{(m)}}} in \overline{B}(R) , R > 0 .

    The fact that HHDBSMPs \{\mathbb{H}_n^{(i, m)}(x)\} does not have the property \mathbb{T}_{\rho_{Q^{(m)}}} in \overline{B}(R) if the condition (2.9) is not satisfied is illustrated by the following example.

    Example 6.3. Let \{Q_{n}^{(m)}(x)\} be the base of SMPs, is defined by:

    \begin{equation*} Q_{n}^{(m)}(x ) = \begin{cases} \mathcal{Q}_{n}^{(m)}(x), &{ n\ is\ even, }\\ \mathcal{Q}_{n}^{(m)}(x)+ \ \frac{\mathcal{Q}_{s(n)}^{(m)}(x)}{2^{(n^{n})}}, &{ n\ is\ odd, } \end{cases} \end{equation*}

    where s(n) is the nearest even integer to n\log n+n^{n}.

    When n is odd, we obtain:

    \mathcal{Q}_{n}^{(m)}(x) = Q_{n}^{(m)}(x)-\frac{Q_{t(n)}^{(m)}(x)}{2^{(n^{n})}}.

    Hence,

    \Psi_{Q_{n}^{(m)}}(R) = R^{n}+2\frac{R^{t(n)}}{2^{(n^{n})}}.

    Putting R = 2 , it follows that

    \Psi_{Q_{n}^{(m)}}(2) = 2^{n}+2^{n\log n+1},

    so that

    \Psi_{Q^{(m)}}(2) = \limsup\limits_{n\rightarrow \infty}\ \frac{log \ \Psi_{Q_{n}^{(m)}}(2)}{n \ log \ n}\leq \log 2.

    It follows that, the base {Q_{n}^{(m)}(x)} has the \mathbb{T}_{\frac{1}{\log2}} -property in \overline{B}(2). The HHDBSMPs \{\mathbb{H}_n^{(i, m)}(x)\} is

    \begin{equation*} \mathbb{H}_n^{(i, m)}(x) = \begin{cases} \zeta_{n, i} \mathcal{Q}_{n-i}^{(m)}(x), & n \ \mbox{is even, }\\ \zeta_{n, i} \mathcal{Q}_{n-i}^{(m)}(x)+\zeta_{t(n), i} \frac{\mathcal{Q}_{t(n)-i}^{(m)}(x)}{2^{(n^{n})}}, &n \ \mbox{is odd.} \end{cases} \end{equation*}

    Hence, when n is odd, we obtain

    \Psi_{\mathbb{H}_{n}^{(i, m)}}(R) = R^{n-i}+2\frac{\zeta_{t(n), i}}{\zeta_{n, i}} \ \frac{R^{t(n)-i}}{2^{(n^{n})}},

    so that when R = 2 ,

    \Psi_{\mathbb{H}_{n}^{(i, m)}}(2) = 2^{n-i}+2\frac{\zeta_{t(n), i}}{\zeta_{n, i}} \ \frac{2^{t(n)-i}}{2^{(n^{n})}}.

    Thus,

    \Psi_{\mathbb{H}^{(i, m)}}(2) = \limsup\limits_{n\rightarrow \infty}\ \frac{\log \ \Psi_{\mathbb{H}_{n}^{(i, m)}}(2)}{n \ log \ n}\leq 1+\log 2,

    and the HHDBSMPs \mathbb{H}_n^{(i, m)}(x), does not have the \mathbb{T}_{\frac{1}{\log2}} -property in \overline{B}(2) as required.

    If the base of SMPs \{Q_{n}^{(m)}(x)\} is simple base ( D_{n} = n ) (see [1]), then the following corollary is a special case of Theorem 6.3.

    Corollary 6.2. When the simple base \{Q_{n}^{(m)}(x)\} of SMPs have \mathbb{T}\rho_{Q^{(m)}} -property in \overline{B}(R) , R > 0 . Then the HHDBSMPs \{\mathbb{H}_n^{(i, m)}(x)\} is also have the \mathbb{T}\rho_{Q^{(m)}} -property.

    In the following, we deduce that the base \{Q_{n}^{(m)}(x)\} and the HHDBs \{\mathbb{H}_{n}^{(i, m)}(x)\} have the same \mathbb{T}\rho_{Q^{(m)}} in an open hyperball B(R) , where R > 0 or at the origin.

    Theorem 6.4. Let \{Q_{n}^{(m)}(x)\} be a base of SMPs have the \mathbb{T}\rho_{Q^{(m)}} -property in B(R) , R > 0 or at the origin. Then the HHDBs \{\mathbb{H}_{n}^{(i, m)}(x)\} have the same property.

    Proof. Let \{Q_{n}^{(m)}(x)\} be have the property \mathbb{T}\rho_{Q^{(m)}} in B(R) , R > 0 , then

    \begin{equation} \Psi_{Q^{(m)}}(r)\leq \frac{1}{\rho_{Q^{(m)}}} \ \ \forall r < R . \end{equation} (6.5)

    It follows from (4.2) that

    \begin{equation} \Psi_{\mathbb{H}^{(i, m)}}(r) = \limsup\limits_{n\rightarrow \infty}\ \frac{\log \ \Psi_{\mathbb{H}_{n}^{(i, m)}}(r)}{n \log \ n}\leq \Psi_{Q^{(m)}}(r_{1}) , \end{equation} (6.6)

    such that r < r_{1} < R. Using (6.5) and (6.6), we have \Psi_{\mathbb{H}_{n}^{(i, m)}}(r)\leq \frac{1}{\rho_{Q^{(m)}}} \forall r < R and the base \{\mathbb{H}_{n}^{(i, m)}(x)\} has the property \mathbb{T}\rho_{Q^{(m)}} in an open hyperball B(R) , R > 0.

    Suppose that the base \{Q_{n}^{(m)}(x)\} has the property \mathbb{T}\rho_{Q^{(m)}} at the origin, then we get

    \begin{equation} \Psi_{Q^{(m)}}(o^{+})\leq \frac{1}{\rho_{Q^{(m)}}}. \end{equation} (6.7)

    Let r_{1}\rightarrow 0^{+} in (6.6), then by (6.7), we have

    \Psi_{\mathbb{H}^{(i, m)}}(o^{+}) \leq \Psi_{Q^{(m)}}(o^{+})\leq \frac{1}{\rho_{Q^{(m)}}},

    and the base \{\mathbb{H}_{n}^{(i, m)}(x)\} has the property \mathbb{T}\rho_{Q^{(m)}} at the origin.

    The problem of classical special functions can be considered as an application of bases of SMPs. Recently, the authors in [19,33] proved that the proper Bessel SMPs (PBSMPs) \{P_{n}^{(m)}(x)\} and the general Bessel SMPs (GBSMPs) \{G_{n}^{(m)}(x)\} are effective for \mathcal{M}[\overline{B}(R)] . Furthermore, recently in [34], the authors proved that the Chebyshev SMPs (CSMPs) \{T_{n}(x)\} is effective for \mathcal{M}[\overline{B}(1)] .

    The following results follows directly by applying Theorem 5.1.

    Corollary 7.1. The base of PBSMPs \{P_{n}^{(m)}(x)\} and the HHD of PBSMPs \{\mathbb{P}_n^{(i, m)}(x)\} have the same region of effectiveness for the class \mathcal{M}[\overline{B}(R)].

    Corollary 7.2. The base of GBSMPs \{G_{n}^{(m)}(x)\} and the HHD of GBSMPs \{\mathbb{G}_n^{(i, m)}(x)\} have the same region of effectiveness for the class \mathcal{M}[\overline{B}(R)].

    Corollary 7.3. The base of CSMPs \{T_{n}(x)\} and the HHD of CSMPs \{\mathbb{T}_n^{(i, m)}(x)\} have the same region of effectiveness for the class \mathcal{M}[\overline{B}(1)].

    In [27] the authors proved that the Bernoulli SMPs (BSMPs) \{B_{n}^{(m)}(x)\} is of order 1 and type \frac{1}{2\pi} and the Euler SMPs (ESMPs)) \{E_{n}^{(m)}(x)\} is of order 1 and type \frac{1}{\pi} .

    According to Theorem 6.1, we obtain the following corollaries:

    Corollary 7.4. The base of BSMPs \{B_{n}^{(m)}(x)\} and the HHD of BSMPs \{\mathbb{B}_{n}^{(i, m)}(x)\} are of the same order 1 and type \frac{1}{2\pi} .

    Corollary 7.5. The base of ESMPs \{B_{n}^{(m)}(x)\} and the HHD of ESMPs \{\mathbb{B}_{n}^{(i, m)}(x)\} are of the same order 1 and type \frac{1}{\pi} .

    Moreover, in [27], the BSMPs \{B_{n}^{(m)}(x)\} and the ESMPs \{E_{n}^{(m)}(x)\} have the property \mathbb{T}_{1} . According to Theorem 6.3, we conclude directly the following corollary:

    Corollary 7.6. If the BSMPs \{B_{n}^{(m)}(x)\} and the ESMPs \{E_{n}^{(m)}(x)\} have the property \mathbb{T}_{1} , then the HHD of BSMPs \{\mathbb{B}_{n}^{(i, m)}(x)\} and ESMPs \{\mathbb{E}_{n}^{(i, m)}(x)\} have the same property, respectively.

    Now, suppose that J_{N}(\mathbb{H}^{(i)}) is a polynomial of the operator \mathbb{H}^{(i)} as given in (3.1) such that

    J_{N}(\mathbb{H}^{(i)}) = \sum\limits_{j = 1}^{N}\lambda_{j} \ (\mathbb{H}^{(i)})^{j}, \ \ \lambda_{i}\in \mathcal{C}_{m},

    where (\mathbb{H}^{(i)})^{j} = (\mathbb{H}^{(i)})^{j-1}\mathbb{H}^{(i)}. Obviously that Theorems 4.1–4.4, 5.1, 6.1, 6.3 and 6.4 will be valid when we replace the base \{\mathbb{H}^{(i)} Q_{n}^{(m)}(x)\} by the base \{J_{N}(\mathbb{H}^{(i)}) Q_{n}^{(m)}(x)\}

    Similar results for the generalized hypercomplex Ruscheweyh derivative base \{J_{N}(\mathfrak{R}^{(i)}) Q_{n}^{(m)}(x)\} , where \mathfrak{R}^{(i)} is the hypercomplex Ruscheweyh derivative. These results generalize the result in [22].

    This work is mainly devoted to derive a generalized form for the Hasse operator in the Clifford setting. Using the defined operator, we accordingly construct the hypercomplex Hasse derivative bases (HHDBs). The approximation properties (effectiveness, order and type, the Property of \mathbb{T}\rho_{Q^{(m)}} ) have been describe for the derived HHDBSMPs in multiple regions in F-modules. Our results are considered as a modified generalization to those given in [8,9,10]. It is clear that that when x\in \mathcal C_{1} in Theorems 4.1–4.4, 5.1, 6.1, 6.3 and 6.4 results obtained by [8,9,10] yield. Additionally considering x to be an element of \mathcal C_{2} in Theorems 4.1–4.4, 5.1, 6.1, 6.3 and 6.4, our results coincide with the quaternion analysis \mathbb{H} . Our results improve and extend the corresponding ones in the Clifford analysis with regards to the region of effectiveness and the mode of increase of HDB (see [13,14]).

    As a result of the growing interest in fractional calculus and its numerous real-world applications, recent contributions were placed on representing analytic functions in terms of complex conformable fractional derivatives and integral bases in different domains in Fréchet spaces [35]. In [36], the authors investigated uncertain barrier swaption pricing problems based on the fractional differential equation in Caputo sense. Relevantly, the fraction Dirac operator constructed using Caput derivative in the case of Clifford variables were studied in [37]. Furthermore, in [38], the authors introduced a new class of time-fractional Dirac type operators with time-variable coefficients. It will be of great interest in the future to explore the convergence properties of fractional derivative bases in the context of Clifford analysis.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through large group Research Project under grant number RGP2/429/44.

    The authors declare no conflict of interest to disclose.



    [1] J. W. Cahn, J. E. Hilliard, Free energy of a nonuniform system. I. interfacial free energy, J. Chem. Phys., 28 (1958), 258–267. https://doi.org/10.1063/1.1744102 doi: 10.1063/1.1744102
    [2] S. M. Allen, J. W. Cahn, A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening, Acta Metall., 27 (1979), 1085–1095. https://doi.org/10.1016/0001-6160(79)90196-2 doi: 10.1016/0001-6160(79)90196-2
    [3] A. Miranville, The Cahn-Hilliard equation and some of its variants, AIMS Math., 2 (2017), 479–544. https://doi.org/10.3934/Math.2017.2.479 doi: 10.3934/Math.2017.2.479
    [4] S. M. Wise, J. S. Lowengrub, H. B. Frieboes, V. Cristini, Three-dimensional multispecies nonlinear tumor growth-i: Model and numerical method, J. Theor. Biol., 253 (2008), 524–543. https://doi.org/10.1016/j.jtbi.2008.03.027 doi: 10.1016/j.jtbi.2008.03.027
    [5] H. Garcke, K. F. Lam, V. Styles, Cahn-hilliard inpainting with the double obstacle potential, SIAM J. Imaging Sci., 11 (2018), 2064–2089. https://doi.org/10.1137/18M1165633 doi: 10.1137/18M1165633
    [6] A. L. Bertozzi, S. Esedoglu, A. Gillette, Inpainting of binary images using the cahn-hilliard equation, IEEE T. Image Proc., 16 (2007), 285–291. https://doi.org/10.1109/TIP.2006.887728 doi: 10.1109/TIP.2006.887728
    [7] L. Q. Chen, Phase-field models for microstructure evolution, Annu. Rev. Mater. Res., 32 (2002), 113–140. https://doi.org/10.1146/annurev.matsci.32.112001.132041 doi: 10.1146/annurev.matsci.32.112001.132041
    [8] J. W. Barrett, J. F. Blowey, Finite element approximation of the Cahn-Hilliard equation with concentration dependent mobility, Math. Comp., 68 (1999), 487–517.
    [9] K. Glasner, S. Orizaga, Improving the accuracy of convexity splitting methods for gradient flow equations, J. Comput. Phys., 315 (2016), 52–64. https://doi.org/10.1016/j.jcp.2016.03.042 doi: 10.1016/j.jcp.2016.03.042
    [10] H. Song, Energy SSP-IMEX Runge-Kutta methods for the Cahn-Hilliard equation, J. Comput. Appl. Math., 292 (2016), 576–590. https://doi.org/10.1016/j.cam.2015.07.030 doi: 10.1016/j.cam.2015.07.030
    [11] J. M. Church, Z. Guo, P. K. Jimack, A. Madzvamuse, K. Promislow, B. Wetton, et al., High accuracy benchmark problems for allen-cahn and cahn-hilliard dynamics, Commun. Comput. Phys., 26 (2019), 947–972. https://doi.org/10.4208/cicp.OA-2019-0006 doi: 10.4208/cicp.OA-2019-0006
    [12] J. Shen, J. Xu, J. Yang, The scalar auxiliary variable (SAV) approach for gradient flows, J. Comput. Phys., 353 (2018), 407–416. https://doi.org/10.1016/j.jcp.2017.10.021 doi: 10.1016/j.jcp.2017.10.021
    [13] S. Orizaga, K. Glasner, Instability and reorientation of block copolymer microstructure by imposed electric fields, Phys. Rev. E, 93 (2016), 052504. https://doi.org/10.1103/PhysRevE.93.052504 doi: 10.1103/PhysRevE.93.052504
    [14] K. R. Elder, M. Grant, Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals, Phys. Rev. E, 70 (2004), 051605. https://doi.org/10.1103/PhysRevE.70.051605 doi: 10.1103/PhysRevE.70.051605
    [15] H. Gomez, X. Nogueira, An unconditionally energy-stable method for the phase field crystal equation, Comput. Meth. Appl. Mech. Eng., 249-252 (2012), 52–61. https://doi.org/10.1016/j.cma.2012.03.002 doi: 10.1016/j.cma.2012.03.002
    [16] H. Gomez, M. Bures, A. Moure, A review on computational modelling of phasetransition problems, Philos. T. Royal Soc. A: Math. Phys. Eng. Sci., 377 (2019), 20180203. https://doi.org/10.1098/rsta.2018.0203 doi: 10.1098/rsta.2018.0203
    [17] G. Gompper, M. Schick, Correlation between structural and interfacial properties of amphiphilic systems, Phys. Rev. Lett., 65 (1990), 1116–1119. https://doi.org/10.1103/PhysRevLett.65.1116 doi: 10.1103/PhysRevLett.65.1116
    [18] W. Feng, Z. Guan, J. Lowengrub, C. Wang, S. M. Wise, Y. Chen, A uniquely solvable, energy stable numerical scheme for the functionalized cahn–hilliard equation and its convergence analysis, J. Sci. Comput., 76 (2018), 1938–1967. https://doi.org/10.1007/s10915-018-0690-1 doi: 10.1007/s10915-018-0690-1
    [19] C. Zhang, J. Ouyang, C. Wang, S. M. Wise, Numerical comparison of modified-energy stable SAV-type schemes and classical BDF methods on benchmark problems for the functionalized cahn-hilliard equation, J. Comput. Phys., 423 (2020), 109772. https://doi.org/10.1016/j.jcp.2020.109772 doi: 10.1016/j.jcp.2020.109772
    [20] A. L. Bertozzi, N. Ju, H. W. Lu, A biharmonic-modified forward time stepping method for fourth order nonlinear diffusion equations, Discrete Contin. Dyn. Syst., 29 (2011), 1367–1391. https://doi.org/10.3934/dcds.2011.29.1367 doi: 10.3934/dcds.2011.29.1367
    [21] D. J. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation, MRS Online Proceedings Library (OPL), 529 (1998), 39. https://doi.org/10.1557/PROC-529-39 doi: 10.1557/PROC-529-39
    [22] S. Dai, Q. Du, Computational studies of coarsening rates for the Cahn-Hilliard equation with phase-dependent diffusion mobility, J. Comput. Phys., 310 (2016), 85–108. https://doi.org/10.1016/j.jcp.2016.01.018 doi: 10.1016/j.jcp.2016.01.018
    [23] H. D. Ceniceros, C. J. Garcia-Cervera, A new approach for the numerical solution of diffusion equations with variable and degenerate mobility, J. Comput. Phys., 246 (2013), 1–10. https://doi.org/10.1016/j.jcp.2013.03.036 doi: 10.1016/j.jcp.2013.03.036
    [24] L. N. Trefethen, Spectral Methods in MatLab, New York: SIAM, 2000.
    [25] G. Sheng, T. Wang, Q. Du, K. G. Wang, Z. K. Liu, L. Q. Chen, Coarsening kinetics of a two phase mixture with highly disparate diffusion mobility, Commun. Comput. Phys., 8 (2010), 249–264. https://doi.org/10.4208/cicp.160709.041109a doi: 10.4208/cicp.160709.041109a
    [26] J. Zhu, L. Q. Chen, J. Shen, V. Tikare, Coarsening kinetics from a variablemobility cahn-hilliard equation: application of a semi-implicit fourier spectral method, Phys. Rev. E, 60 (1999), 3564. https://doi.org/10.1103/PhysRevE.60.3564 doi: 10.1103/PhysRevE.60.3564
    [27] D. Li, Z. Qiao, T. Tang, Characterizing the stabilization size for semi-implicit Fourier-spectral method to phase field equations, SIAM J. Numer. Anal., 54 (2016), 1653–1681. https://doi.org/10.1137/140993193 doi: 10.1137/140993193
    [28] D. Li, Z. Qiao, On second order semi-implicit fourier spectral methods for 2d Cahn-Hilliard equations, J. Sci. Comput., 70 (2017), 301–341. https://doi.org/10.1007/s10915-016-0251-4 doi: 10.1007/s10915-016-0251-4
    [29] Z. Xu, H. Zhang, Stabilized semi-implicit numerical schemes for the Cahn-Hilliard-like equation with variable interfacial parameter, J. Comput. Appl. Math., 346 (2019), 307–322. https://doi.org/10.1016/j.cam.2018.06.031 doi: 10.1016/j.cam.2018.06.031
    [30] L. Duchemin, J. Eggers, The explicit-implicit-null method: Removing the numerical instability of PDEs, J. Comput. Phys., 263 (2014), 37–52. https://doi.org/10.1016/j.jcp.2014.01.013 doi: 10.1016/j.jcp.2014.01.013
    [31] Y. Yan, W. Chen, C. Wang, S. M. Wise, A second-order energy stable BDF numerical scheme for the Cahn-Hilliard equation, Commun. Comput. Phys., 23 (2018), 572–602. https://doi.org/10.4208/cicp.OA-2016-0197 doi: 10.4208/cicp.OA-2016-0197
    [32] C. Canuto, M. Y. Hussaini, A. Quarteroni, T. A. Zang, Spectral Methods: Fundamentals in Single Domains, Berlin: Springer Science Business Media, 2007.
    [33] J. Shen, J. Xu, J. Yang, A new class of efficient and robust energy stable schemes for gradient flows, SIAM Rev., 61 (2019), 474–506. https://doi.org/10.1137/17M1150153 doi: 10.1137/17M1150153
    [34] A. Christlieb, J. Jones, K. Promislow, B. Wetton, M. Willoughby, High accuracy solutions to energy gradient flows from material science models, J. Comput. Phys., 257 (2014), 193–215. https://doi.org/10.1016/j.jcp.2013.09.049 doi: 10.1016/j.jcp.2013.09.049
    [35] A. Doelman, G. Hayrapetyan, K. Promislow, B. Wetton, Meander and pearling of single-curvature bilayer interfaces in the functionalized cahn-hilliard equation, SIAM J. Math. Anal., 46 (2014), 3640–3677. https://doi.org/10.1137/13092705X doi: 10.1137/13092705X
    [36] C. Zhang, J. Ouyang, Unconditionally energy stable second-order numerical schemes for the functionalized cahn–hilliard gradient flow equation based on the SAV approach, Comput. Math. Appl., 84 (2021), 16–38. https://doi.org/10.1016/j.camwa.2020.12.003 doi: 10.1016/j.camwa.2020.12.003
    [37] K. Promislow, Q. Wu, Existence of pearled patterns in the planar functionalized cahn-hilliard equation, J. Diff. Equ., 259 (2015), 3298–3343. https://doi.org/10.1016/j.jde.2015.04.022 doi: 10.1016/j.jde.2015.04.022
    [38] Z. Tan, L. Chen, J. Yang, Generalized allen-cahn-type phase-field crystal model with fcc ordering structure and its conservative high-order accurate algorithm, Comput. Phys. Commun., 286 (2023), 108656. https://doi.org/10.1016/j.cpc.2023.108656 doi: 10.1016/j.cpc.2023.108656
    [39] N. Gavish, J. Jones, Z. Xu, A. Christlieb, K. Promislow, Variational models of network formation and ion transport: Applications to perfluorosulfonate ionomer membranes, Polymers, 4 (2012), 630–655. https://doi.org/10.3390/polym4010630 doi: 10.3390/polym4010630
    [40] J. Burkardt, M. Gunzburger, W. Zhao, High-precision computation of the weak galerkin methods for the fourth-order problem, Numer. Algor., 84 (2020), 181–205. https://doi.org/10.1007/s11075-019-00751-5 doi: 10.1007/s11075-019-00751-5
  • This article has been cited by:

    1. Mohra Zayed, Gamal Hassan, Kronecker product bases and their applications in approximation theory, 2025, 33, 2688-1594, 1070, 10.3934/era.2025048
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1045) PDF downloads(62) Cited by(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog