Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Unconditionally stable monte carlo simulation for solving the multi-dimensional Allen–Cahn equation

  • Received: 13 April 2023 Revised: 27 June 2023 Accepted: 04 July 2023 Published: 17 July 2023
  • In this study, we present an efficient and novel unconditionally stable Monte Carlo simulation (MCS) for solving the multi-dimensional Allen–Cahn (AC) equation, which can model the motion by mean curvature flow of a hypersurface. We use an operator splitting method, where the diffusion and nonlinear terms are solved separately. The diffusion term is calculated using MCS for the stochastic differential equation, while the nonlinear term is locally computed for each particle in a virtual grid. Several numerical experiments are presented to demonstrate the performance of the proposed algorithm. The computational results confirm that the proposed algorithm can solve the AC equation more efficiently as the dimension of space increases.

    Citation: Youngjin Hwang, Ildoo Kim, Soobin Kwak, Seokjun Ham, Sangkwon Kim, Junseok Kim. Unconditionally stable monte carlo simulation for solving the multi-dimensional Allen–Cahn equation[J]. Electronic Research Archive, 2023, 31(8): 5104-5123. doi: 10.3934/era.2023261

    Related Papers:

    [1] Jingwen He, Shirley Simon, Feng-Kuang Chiang . A comparative study of pre-service teachers' perceptions on STEAM education in UK and China. STEM Education, 2022, 2(4): 318-344. doi: 10.3934/steme.2022020
    [2] Ana Barbosa, Isabel Vale, Dina Alvarenga . The use of Tinkercad and 3D printing in interdisciplinary STEAM education: A focus on engineering design. STEM Education, 2024, 4(3): 222-246. doi: 10.3934/steme.2024014
    [3] Vincent Theodore M. Balo, Joje Mar P. Sanchez . Evaluation of educational assessment module for flexible STEM education. STEM Education, 2025, 5(1): 130-151. doi: 10.3934/steme.2025007
    [4] Ana Barbosa, Isabel Vale . Rebuilding manipulatives through digital making in teacher education. STEM Education, 2025, 5(4): 515-545. doi: 10.3934/steme.2025025
    [5] Hyunkyung Kwon, Yujin Lee . A meta-analysis of STEM project-based learning on creativity. STEM Education, 2025, 5(2): 275-290. doi: 10.3934/steme.2025014
    [6] Usman Ghani, Xuesong Zhai, Riaz Ahmad . Mathematics skills and STEM multidisciplinary literacy: Role of learning capacity. STEM Education, 2021, 1(2): 104-113. doi: 10.3934/steme.2021008
    [7] Sasha Nikolic, Zach Quince, Anna Lidfors Lindqvist, Peter Neal, Sarah Grundy, May Lim, Faham Tahmasebinia, Shannon Rios, Josh Burridge, Kathy Petkoff, Ashfaque Ahmed Chowdhury, Wendy S.L. Lee, Rita Prestigiacomo, Hamish Fernando, Peter Lok, Mark Symes . Project-work Artificial Intelligence Integration Framework (PAIIF): Developing a CDIO-based framework for educational integration. STEM Education, 2025, 5(2): 310-332. doi: 10.3934/steme.2025016
    [8] Giuseppe Carbone, Elio Matteo Curcio, Stefano Rodinò, Francesco Lago . A Robot-Sumo student competition at UNICAL as a learning-by-doing strategy for STEM education. STEM Education, 2022, 2(3): 262-274. doi: 10.3934/steme.2022016
    [9] Ibrahim Khalil, Amirah AL Zahrani, Bakri Awaji, Mohammed Mohsen . Teachers' perceptions of teaching mathematics topics based on STEM educational philosophy: A sequential explanatory design. STEM Education, 2024, 4(4): 421-444. doi: 10.3934/steme.2024023
    [10] Tanya Evans, Heiko Dietrich . Inquiry-based mathematics education: a call for reform in tertiary education seems unjustified. STEM Education, 2022, 2(3): 221-244. doi: 10.3934/steme.2022014
  • In this study, we present an efficient and novel unconditionally stable Monte Carlo simulation (MCS) for solving the multi-dimensional Allen–Cahn (AC) equation, which can model the motion by mean curvature flow of a hypersurface. We use an operator splitting method, where the diffusion and nonlinear terms are solved separately. The diffusion term is calculated using MCS for the stochastic differential equation, while the nonlinear term is locally computed for each particle in a virtual grid. Several numerical experiments are presented to demonstrate the performance of the proposed algorithm. The computational results confirm that the proposed algorithm can solve the AC equation more efficiently as the dimension of space increases.



    We analyze the oscillation of linear differential equations in the form

    [r(t)x(t)]+s(t)x(t)=0, (1.1)

    where the coefficients r>0,s are continuous on an interval [τ,). We point out that Eq (1.1) is called oscillatory if all solutions are oscillatory (i.e., any solution has infinitely many zero points in any neighbourhood of ). In the other case (any non-trivial solution has the biggest zero point), we say that Eq (1.1) is non-oscillatory. Concerning the basics of the oscillation theory of Eq (1.1), we refer, e.g., to [3,51] with references cited therein. We focus on the study of conditionally oscillatory equations which are treated as "ideal testing equations" (it is described, e.g., in [21]).

    To recall the notion of the conditional oscillation, we consider Eq (1.1) in the following modified form

    [r(t)x(t)]+γs(t)x(t)=0, (1.2)

    where γR is a parameter. We say that Eq (1.2) is conditionally oscillatory if there exists the so-called critical oscillation constant Γ>0 such that Eq (1.2) is oscillatory for γ>Γ and non-oscillatory for γ<Γ. Note that the case γ=Γ is not covered by the definition above. Some simple equations are non-oscillatory in the borderline case (see, e.g., [27,35,37,48]). Nevertheless, many equations can be oscillatory in the borderline case (see, e.g., [11,13,14,33] and also [28,50] in the discrete case).

    In this paragraph, we mention examples of conditionally oscillatory equations together with the corresponding critical oscillation constants Γ. We add that log denotes the natural logarithm and p>0 stands for an arbitrarily given number in the whole paper. In [47] (see also [46,49]), there is shown that the equation

    x(t)+γsinttx(t)=0

    is conditionally oscillatory for Γ=1/2 and the equation

    x(t)+γsin(t2)x(t)=0

    is conditionally oscillatory for Γ=2. In [16,38], it is proved that the equation

    [r(t)x(t)]+γ1t2s(t)x(t)=0,

    where the coefficients r,s are α-periodic (for some α>0) and positive, is conditionally oscillatory for

    Γ=14(1αα0dτr(τ))1(1αα0s(τ)dτ)1. (1.3)

    Analogously (see [24]), for α-periodic and positive functions r,s, the equation

    [tr(t)x(t)]+γ1tlog2ts(t)x(t)=0

    is conditionally oscillatory for the critical oscillation constant in (1.3) and, from [8] (see also [27,Corollary 4.2]), it follows that the equation

    [tqr(t)x(t)]+γtqt2s(t)x(t)=0,qR{1},

    is conditionally oscillatory for

    Γ=(q1)24(1αα0dτr(τ))1(1αα0s(τ)dτ)1.

    For other relevant results about the considered conditional oscillation, we refer at least to [15,31,34,39] (and also [19,40]). However, the main motivation for our current research comes from [21], where the following theorem is proved.

    Theorem 1.1. Let us considerthe equation

    [logptr(t)x(t)]+logptt2s(t)x(t)=0, (1.4)

    where r,s:R(0,) are continuous and periodic functions with period α>0.

    (A) If

    4(1αα0r(τ)dτ)(1αα0s(τ)dτ)>1,

    then Eq (1.4) is oscillatory.

    (B) If

    4(1αα0r(τ)dτ)(1αα0s(τ)dτ)<1,

    then Eq (1.4) is non-oscillatory.

    The aim of this paper is to extend Theorem 1.1, i.e., to identify a more general class of conditionally oscillatory linear equations. In fact, we consider equations in the form of Eq (1.4) for non-periodic functions r,s, for periodic r,s having different periods and for periodic s changing its sign. All these cases are covered. To prove such a result, we use the modified Prüfer angle in this paper. Note that Theorem 1.1 is proved using the Riccati transformation, i.e., the methods are dissimilar.

    To complete the literature overview, we add that conditionally oscillatory equations are studied in the field of difference equations and dynamic equations on time scales (see, e.g., [4,12,20,45] and also [1,2,18,25,52,53] for possible research directions). Concerning more general half-linear conditionally oscillatory equations, we refer to [9,17,23,29] in the continuous case and to [22,26,36,44] in the discrete case; concerning non-linear equations, we point out at least [5,32,41,42] (see also [6,43] for possible research directions).

    This paper is organized as follows. In the next section, we present the main tool of our paper, i.e., the used version of the Prüfer angle. In Section 3, we collect all auxiliary results. The main result is proved in Section 4, where we also mention an example of conditionally oscillatory equations. This example demonstrates how our result substantially generalizes Theorem 1.1. In the last section, we present new corollaries together with simple illustrative examples of linear equations whose oscillation properties can be deduced from the corollaries and do not follow from any previously known results.

    Since we study the oscillation and non-oscillation of differential equations, it suffices to consider all equations only for t large enough. For simplicity, we consider te, where e is the base of the natural logarithm. We also put Re:=[e,).

    We consider linear second order differential equations of the form

    [R(t)x(t)]+S(t)x(t)=0, (2.1)

    where R,S:ReR are continuous functions and R is positive. Let x be a non-trivial solution of Eq (2.1). For x(t)0, the well-known Riccati transformation

    w(t)=R(t)x(t)x(t) (2.2)

    gives

    w(t)+S(t)+R1(t)w2(t)=0. (2.3)

    For details, we can refer, e.g., to [30].

    Applying the substitution

    v(t)=tlogptw(t), (2.4)

    from Eq (2.3), we have

    v(t)=(tlogpt)w(t)+tlogptw(t)=logtplogp+1tw(t)tlogpt[S(t)+R1(t)w2(t)]=logtptlogtv(t)tlogptS(t)R1(t)logpttv2(t). (2.5)

    Then, for a non-trivial solution x of Eq (2.1), we consider the modified Prüfer transformation in the form

    x(t)=ρ(t)sinφ(t),x(t)=ρ(t)R1(t)logpttcosφ(t). (2.6)

    We have (see (2.2), (2.4) and (2.6))

    v(t)=tlogptw(t)=tlogptR(t)x(t)x(t)=tlogptR(t)ρ(t)R1(t)logpttcosφ(t)ρ(t)sinφ(t)=cosφ(t)sinφ(t)=cotφ(t), (2.7)

    i.e.,

    v(t)=1sin2φ(t)φ(t). (2.8)

    Thus (see (2.5) and (2.8)), we obtain

    1sin2φ(t)φ(t)=logtptlogtv(t)tlogptS(t)R1(t)logpttv2(t) (2.9)

    and (see (2.7) and (2.9))

    φ(t)=logtptlogtcosφ(t)sinφ(t)+tlogptS(t)sin2φ(t)+R1(t)logpttcos2φ(t),

    i.e.,

    φ(t)=logpttR(t)cos2φ(t)logtptlogtcosφ(t)sinφ(t)+tS(t)logptsin2φ(t). (2.10)

    Let the functions R and S have the forms

    R(t)=logptr(t),S(t)=logptt2s(t),tRe, (2.11)

    i.e., let us consider the equation

    [logptr(t)x(t)]+logptt2s(t)x(t)=0, (2.12)

    where continuous functions r:Re(0,) and s:ReR satisfy

    limtt+1tr(τ)dτt=limtt+1t|s(τ)|dτt=0. (2.13)

    Note that (2.13) implies

    limtt+αtr(τ)dτt=limtt+αt|s(τ)|dτt=0 (2.14)

    for any α>0.

    Finally, for the coefficients in (2.11), Eq (2.10) has the form

    φ(t)=1t[r(t)cos2φ(t)logtplogtcosφ(t)sinφ(t)+s(t)sin2φ(t)], (2.15)

    which is the considered form of the equation for the modified Prüfer angle, i.e., we use Eq (2.15) to study Eq (2.12). In addition, for given α>0 and a solution φ of Eq (2.15) on Re, we use the auxiliary average function ψα:ReR defined by the formula

    ψα(t):=1αt+αtφ(τ)dτ,tRe. (2.16)

    In this section, we collect all used lemmas. At first, we mention a known result (in the form which plays a crucial role in the proof of our main result).

    Lemma 3.1. Let us considerEq (2.12) together with Eq (2.15). Let φ:ReR be asolution of Eq (2.15).

    (A) If limtφ(t)=, then Eq (2.12) is oscillatory.

    (B) If lim suptφ(t)<, then Eq (2.12) is non-oscillatory.

    Proof. The non-oscillation of Eq (2.12) is equivalent to the boundedness from above of the Prüfer angle φ given by Eq (2.15). See, e.g., [37] or directly consider the transformation in (2.6) and Eq (2.15) for sinφ(t)=0. In addition (see the form of Eq (2.15)), the set of all values of φ is unbounded if and only if limtφ(t)=.

    One can easily reformulate Lemma 3.1 as follows.

    Lemma 3.2. Let us considerEq (2.12) together with Eq (2.15). Let φ:ReR be asolution of Eq (2.15).

    (A) If Eq (2.12) is oscillatory, then limtφ(t)=.

    (B) If Eq (2.12) is non-oscillatory, then lim suptφ(t)<.

    Now we mention a simple corollary of Theorem 1.1.

    Lemma 3.3. Let C,D>0. Let us considerthe equation

    [logptCx(t)]+logptt2Dx(t)=0. (3.1)

    (A) If 4CD>1, then Eq (3.1) is oscillatory.

    (B) If 4CD<1, then Eq (3.1) is non-oscillatory.

    Proof. It suffices to consider Theorem 1.1 for constant functions r,s.

    The two previous lemmas give the following result.

    Lemma 3.4. Let C,D>0. Let η:ReR be a solution ofthe equation

    η(t)=1t[Ccos2η(t)logtplogtcosη(t)sinη(t)+Dsin2η(t)]. (3.2)

    (A) If 4CD>1, then limtη(t)=.

    (B) If 4CD<1, then lim suptη(t)<.

    Proof. The statement of the lemma follows immediately from Lemmas 3.2 and 3.3, because Eq (3.2) has the form of the equation for the modified Prüfer angle η which corresponds to Eq (3.1).

    Next, we deal with the auxiliary average function ψα.

    Lemma 3.5. Let φ:ReR be a solution of Eq (2.15) and ψα be defined in (2.16). The limit

    limtt|φ(τ)ψα(t)|=0 (3.3)

    exists uniformly with respect to τ[t,t+α].

    Proof. The continuity of φ gives

    |φ(τ)ψα(t)|maxτ1,τ2[0,α]|φ(t+τ1)φ(t+τ2)|,tRe,τ[t,t+α].

    Hence, we obtain (see (2.14) and Eq (2.15))

    lim suptt|φ(τ)ψα(t)|lim supttmaxτ1,τ2[0,α]|φ(t+τ1)φ(t+τ2)|lim supttt+αt|φ(τ)|dτ=lim supttt+αt|1τ[r(τ)cos2φ(τ)logτplogτcosφ(τ)sinφ(τ)+s(τ)sin2φ(τ)]|dτlim suptttt+αt[r(τ)+1+plogt+|s(τ)|]dτ=lim supt(1tt+αtr(τ)dτ+1tt+αt|s(τ)|dτ)=0

    uniformly with respect to τ[t,t+α].

    Lemma 3.6. Let φ:ReR be a solution of Eq (2.15) and ψα be defined in (2.16). Then, there exists a continuousfunction F:(e,)R satisfying

    limtF(t)=0 (3.4)

    and

    ψα(t)=1t[(1αt+αtr(τ)dτ)cos2ψα(t)logtplogtcosψα(t)sinψα(t)+(1αt+αts(τ)dτ)sin2ψα(t)+F(t)] (3.5)

    for all t>e.

    Proof. In fact, the considered continuous function F:(e,)R can be introduced as the function for which (3.5) holds, i.e., one can put

    F(t):=tψα(t)(1αt+αtr(τ)dτ)cos2ψα(t)+logtplogtcosψα(t)sinψα(t)(1αt+αts(τ)dτ)sin2ψα(t) (3.6)

    for all t>e. The aim is to prove (3.4) for this function F. For all t>e, we have (see Eq (2.15) and (2.16))

    ψα(t)=1α(φ(t+α)φ(t))=1αt+αtφ(τ)dτ=1αt+αt1τ[r(τ)cos2φ(τ)logτplogτcosφ(τ)sinφ(τ)+s(τ)sin2φ(τ)]dτ

    and

    |tψα(t)1αt+αt[r(τ)cos2φ(τ)logτplogτcosφ(τ)sinφ(τ)+s(τ)sin2φ(τ)]dτ|=tα|t+αt1τ[r(τ)cos2φ(τ)logτplogτcosφ(τ)sinφ(τ)+s(τ)sin2φ(τ)]dτt+αt1t[r(τ)cos2φ(τ)logτplogτcosφ(τ)sinφ(τ)+s(τ)sin2φ(τ)]dτ|tαt+αt(1t1τ)[r(τ)+1+plogt+|s(τ)|]dτ1tt+αt[r(τ)+1+p+|s(τ)|]dτ=1t[1tt+αtr(τ)dτ+α(1+p)t+1tt+αt|s(τ)|dτ],

    which implies (see (2.14))

    limt|tψα(t)1αt+αt[r(τ)cos2φ(τ)logτplogτcosφ(τ)sinφ(τ)+s(τ)sin2φ(τ)]dτ|=0.

    Therefore (see directly (3.6)), (3.4) can be proved by

    limt(1αt+αtr(τ)cos2ψα(t)dτ1αt+αtr(τ)cos2φ(τ)dτ)=0, (3.7)
    limt(logtplogtcosψα(t)sinψα(t)1αt+αtlogτplogτcosφ(τ)sinφ(τ)dτ)=0, (3.8)
    limt(1αt+αts(τ)sin2ψα(t)dτ1αt+αts(τ)sin2φ(τ)dτ)=0. (3.9)

    To obtain (3.7) and (3.9), we use Lemma 3.5 and the Lipschitz continuity of y=cos2x and y=sin2x which gives L>0 such that

    |cos2x1cos2x2|L|x1x2|,x1,x2R, (3.10)

    and that

    |sin2x1sin2x2|L|x1x2|,x1,x2R. (3.11)

    Thus, we have (see (2.14), (3.3) and (3.10))

    lim supt|1αt+αtr(τ)cos2ψα(t)dτ1αt+αtr(τ)cos2φ(τ)dτ|lim supt1αt+αtr(τ)|cos2ψα(t)cos2φ(τ)|dτlim supt1tt+αtr(τ)Lαt|ψα(t)φ(τ)|dτlim supt1tt+αtr(τ)dτ=0 (3.12)

    and (see (2.14), (3.3) and (3.11))

    lim supt|1αt+αts(τ)sin2ψα(t)dτ1αt+αts(τ)sin2φ(τ)dτ|lim supt1αt+αt|s(τ)||sin2ψα(t)sin2φ(τ)|dτlim supt1tt+αt|s(τ)|Lαt|ψα(t)φ(τ)|dτlim supt1tt+αt|s(τ)|dτ=0. (3.13)

    It is seen that (3.12) and (3.13) affirms (3.7) and (3.9), respectively.

    Similarly, we use the Lipschitz continuity of y=cosxsinx which implies that there exists M>0 for which

    |cosx1sinx1cosx2sinx2|M|x1x2|,x1,x2R. (3.14)

    Now (3.8) follows from

    lim supt|logtplogtcosψα(t)sinψα(t)1αt+αtlogτplogτcosφ(τ)sinφ(τ)dτ|lim supt|logtplogtcosψα(t)sinψα(t)1αt+αtlogtplogtcosφ(τ)sinφ(τ)dτ|+lim supt|1αt+αtlogtplogtcosφ(τ)sinφ(τ)dτ1αt+αtlogτplogτcosφ(τ)sinφ(τ)dτ|(p+1)lim supt1αt+αt|cosψα(t)sinψα(t)cosφ(τ)sinφ(τ)|dτ+lim supt1αt+αt|logtplogtlogτplogτ|dτM(p+1)lim supt1t(1αt+αtt|ψα(t)φ(τ)|dτ)+plim supt1αt+αt(1logt1log(t+α))dτmM(p+1)lim supt1t+plim suptαlogtlog(t+α)=0,

    where (3.3) (see Lemma 3.5) and (3.14) are applied.

    We prove the conditional oscillation of Eq (2.12), i.e., the announced generalization of Theorem 1.1. This main result follows.

    Theorem 4.1. Let us consider Eq (2.12) with continuous coefficients r:Re(0,) and s:ReR satisfying (2.13). Let A,B,α>0.

    (A) If 4AB>1 and if the inequalities

    1αt+αtr(τ)dτA,1αt+αts(τ)dτB (4.1)

    are valid for all large t, then Eq (2.12) is oscillatory.

    (B) If 4AB<1 and if the inequalities

    1αt+αtr(τ)dτA,1αt+αts(τ)dτB (4.2)

    are valid for all large t, then Eq (2.12) is non-oscillatory.

    Proof. In both parts of the proof, we use the equation for the modified Prüfer angle which corresponds to Eq (2.12), i.e., we consider Eq (2.15). Let φ:ReR be a solution of Eq (2.15) and let ψα:ReR be the corresponding average function introduced in (2.16). Based on Lemma 3.1, we analyze the limit behaviour of φ. In addition, due to Lemma 3.5, it suffices to show that ψα is unbounded in the first case and bounded from above in the second case.

    We begin with the first part (i.e., the oscillation part) of the theorem. Considering (3.5) in Lemma 3.6 and (4.1), we obtain

    ψα(t)=1t[(1αt+αtr(τ)dτ)cos2ψα(t)logtplogtcosψα(t)sinψα(t)+(1αt+αts(τ)dτ)sin2ψα(t)+F(t)]1t[Acos2ψα(t)logtplogtcosψα(t)sinψα(t)+Bsin2ψα(t)+F(t)]=1t[(A+F(t))cos2ψα(t)logtplogtcosψα(t)sinψα(t)+(B+F(t))sin2ψα(t)]

    for all large t, where limtF(t)=0 (see (3.4)). Let δ>0 be so small that

    Aδ,Bδ>0,4(Aδ)(Bδ)>1 (4.3)

    and t so large that |F(t)|<δ. Then,

    ψα(t)>1t[(Aδ)cos2ψα(t)logtplogtcosψα(t)sinψα(t)+(Bδ)sin2ψα(t)] (4.4)

    for all large t. Putting C=Aδ and D=Bδ, from (4.3), we have C,D>0 and 4CD>1. We apply Lemma 3.4, (A). Thus, any solution ζ:ReR of the equation

    ζ(t)=1t[(Aδ)cos2ζ(t)logtplogtcosζ(t)sinζ(t)+(Bδ)sin2ζ(t)] (4.5)

    satisfies limtζ(t)=. Let T>e be such that (4.4) is valid for all tT. Since y=sin2x, y=cos2x, y=cosxsinx are π-period functions, comparing Eq (4.5) with the right-hand side of (4.4), we have

    ψα(t)ζ(t)|ψα(T)ζ(T)|π

    for all tT. Therefore,

    lim inftψα(t)limtζ(t)=.

    Now we prove the second part. We proceed analogously as in the first part. Applying Lemma 3.6 and (4.2), we have

    ψα(t)=1t[(1αt+αtr(τ)dτ)cos2ψα(t)logtplogtcosψα(t)sinψα(t)+(1αt+αts(τ)dτ)sin2ψα(t)+F(t)]1t[Acos2ψα(t)logtplogtcosψα(t)sinψα(t)+Bsin2ψα(t)+F(t)]=1t[(A+F(t))cos2ψα(t)logtplogtcosψα(t)sinψα(t)+(B+F(t))sin2ψα(t)]

    for all large t, where limtF(t)=0 (see (3.4)). Let δ>0 satisfy 4(A+δ)(B+δ)<1 and t be such that |F(t)|<δ. Then,

    ψα(t)<1t[(A+δ)cos2ψα(t)logtplogtcosψα(t)sinψα(t)+(B+δ)sin2ψα(t)] (4.6)

    for all large t. We put C=A+δ, D=B+δ. Considering 4CD<1, Lemma 3.4, (B) says that any solution θ:ReR of the equation

    θ(t)=1t[(A+δ)cos2θ(t)logtplogtcosθ(t)sinθ(t)+(B+δ)sin2θ(t)] (4.7)

    satisfies lim suptθ(t)<. From (4.6) and (4.7), we have

    lim suptψα(t)lim suptθ(t)+|ψα(T)θ(T)|+π<,

    where T>e is a sufficiently large number.

    Remark 1. Concerning the statement of Theorem 4.1, we remark that the requirement (see (2.13))

    limtt+1tr(τ)dτt=0 (4.8)

    can be omitted (in contrast to the second requirement about s in (2.13)), i.e., it is evident that Theorem 4.1 is valid also without this limitation provided we have proved its statement with (4.8). Note that we use (4.8) in the proofs of Lemmas 3.5 and 3.6.

    Remark 2. Concerning the conditional oscillation of perturbed Euler type equations (see, e.g., [7,10]), we conjecture that it is not possible to decide the (non-)oscillation of Eq (2.12) for general functions r,s satisfying

    limt(1αt+αtr(τ)dτ)(1αt+αts(τ)dτ)=14,

    where α>0. We can also refer to [27].

    Theorem 4.1 covers equations with unbounded coefficients which oscillate non-trivially. To illustrate this fact, we provide the following example.

    Example 1. For μ>0 and νR, we define

    r(t):=μ[1+3n(tn)],t[n,n+12n),nN;r(t):=μ[1+3n(n+22nt)],t[n+12n,n+22n],nN;r(t):=μ,t(n+22n,n+1),nN,

    and

    s(t):=ν+4nnlog(n+1)(tn4i4n),t[n+4i4n,n+1+4i4n),i{0,1,,4n11},nN;s(t):=ν+4nnlog(n+1)(n+2+4i4nt),t[n+1+4i4n,n+3+4i4n),i{0,1,,4n11},nN;s(t):=ν+4nnlog(n+1)(tn4+4i4n),t[n+3+4i4n,n+4+4i4n),i{0,1,,4n11},nN.

    For these functions, let us consider Eq (2.12). Let α=1, i.e., let the average ψα be taken over intervals of the length α=1. One can easily verify that

    limtt+1tr(τ)dτ=μ,limtt+1ts(τ)dτ=ν,limtt+1t|s(τ)|dτt=0.

    Applying Theorem 4.1 (consider also Remark 1), we obtain the oscillation of the considered equation for 4μν>1 and its non-oscillation for 4μν<1.

    In this section, to explain the novelty of our main result, we mention its direct corollaries together with simple examples. We point out that none of the examples below is covered by previously known results and that all results below are new for any p>0. At first, we recall the concept of the mean value for continuous functions which is used in the first corollary.

    Definition 1. Let a continuous function f:ReR be such that the limit

    M(f):=lima1at+atf(τ)dτ

    is finite and exists uniformly with respect to tRe. The number M(f) is called the mean value of f.

    Corollary 5.1. Let us consider Eq (2.12), where continuous functions r:Re(0,) and s:ReR have mean values M(r),M(s) and

    limtt+1t|s(τ)|dτt=0.

    (A) If 4M(r)M(s)>1, then Eq (2.12) is oscillatory.

    (B) If 4M(r)M(s)<1, then Eq (2.12) is non-oscillatory.

    Proof. According to the well-known Sturm theory (see, e.g., [51]), we can assume that M(r)>0 and M(s)>0. The corollary follows from Theorem 4.1 (and from Remark 1). Indeed, for δ>0 such that M(r)δ,M(s)δ>0 and 4(M(r)δ)(M(s)δ)>1 or 4(M(r)+δ)(M(s)+δ)<1, from Definition 1, we obtain the existence of α>0 with the property that

    1αt+αtr(τ)dτM(r)δ,1αt+αts(τ)dτM(s)δ

    or

    1αt+αtr(τ)dτM(r)+δ,1αt+αts(τ)dτM(s)+δ

    for all t, respectively.

    There are known many types of continuous functions which have mean values. For example, we recall asymptotically almost periodic functions which are treated in the following example.

    Example 2. For μ2 and νR, let us consider the equation

    [logp(t+2)μ+sint+cos(2t)x(t)]+logptν+sin(3t)+cos(2t)+arctan t+1t2(t+t)2x(t)=0, (5.1)

    which has the form of Eq (2.12) for

    r(t)=(μ+sint+cos(2t))(logtlog(t+2))p,tRe,
    s(t)=(ν+sin(3t)+cos(2t)+arctan t+1t2)t2(t+t)2,tRe.

    Since

    M(r)=μ+M(sint)+M(cos(2t))=μ

    and

    M(s)=ν+M(sin(3t))+M(cos(2t))+M(arctan t+1t2)=ν,

    Corollary 5.1 guarantees the oscillation of Eq (5.1) for 4μν>1 and the non-oscillation of Eq (5.1) for 4μν<1.

    Now we concentrate on more concrete cases of the studied equations with periodic coefficients (nonpositive or with different periods). We highlight that these results are new as well.

    Corollary 5.2. Let us consider Eq (2.12), where continuous functions r:Re(0,) and s:ReR are α-periodic for α>0.

    (A) If

    4(1αe+αer(τ)dτ)(1αe+αes(τ)dτ)>1,

    then Eq (2.12) is oscillatory.

    (B) If

    4(1αe+αer(τ)dτ)(1αe+αes(τ)dτ)<1,

    then Eq (2.12) is non-oscillatory.

    Proof. It suffices to consider Corollary 5.1 for

    M(r)=1αe+αer(τ)dτ,M(s)=1αe+αes(τ)dτ

    and the boundedness of s.

    Corollary 5.3. Let us consider Eq (2.12), where r:Re(0,) is a continuousand α-periodic function and s:Re(0,) is a continuous and β-periodic function for α,β>0.

    (A) If

    4(1αe+αer(τ)dτ)(1βe+βes(τ)dτ)>1,

    then Eq (2.12) is oscillatory.

    (B) If

    4(1αe+αer(τ)dτ)(1βe+βes(τ)dτ)<1,

    then Eq (2.12) is non-oscillatory.

    Proof. It suffices to consider Corollary 5.1 for

    M(r)=1αe+αer(τ)dτ,M(s)=1βe+βes(τ)dτ

    and the boundedness of s.

    Using the above corollaries, one can identify the critical oscillation constants for several equations. We mention at least the examples below.

    Example 3. For μ>1, let us consider the equations

    [logptμ+sintx(t)]+logpt18+sintt2x(t)=0,[logptμ+sintx(t)]+logpt18+costt2x(t)=0,[logptμ+costx(t)]+logpt18+sintt2x(t)=0,[logptμ+costx(t)]+logpt18+costt2x(t)=0.

    For these equations, one can apply Corollary 5.2. Therefore, the equations are oscillatory if μ>2, and they are non-oscillatory if μ<2.

    Example 4. For ν>1, let us consider the equations

    [logptν+sintx(t)]+logpt2+sin(πt)16t2x(t)=0,[logptν+sintx(t)]+logpt2+cos(πt)16t2x(t)=0,[logptν+costx(t)]+logpt2+sin(πt)16t2x(t)=0,[logptν+costx(t)]+logpt2+cos(πt)16t2x(t)=0.

    Now one can apply Corollary 5.3. Thus, the equations are oscillatory for ν>2 and non-oscillatory for ν<2.

    The research presented in this paper was supported by the Czech Science Foundation (grant no. GA20-11846S).

    All authors declare no conflicts of interest in this paper.



    [1] 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
    [2] M. Olshanskii, X. Xu, V. Yushutin, A finite element method for Allen–Cahn equation on deforming surface, Comput. Math. Appl., 90 (2021), 148–158. https://doi.org/10.1016/j.camwa.2021.03.018 doi: 10.1016/j.camwa.2021.03.018
    [3] J. W. Choi, H. G. Lee, D. Jeong, J. Kim, An unconditionally gradient stable numerical method for solving the Allen–Cahn equation, Physica A, 388 (2009), 1791–1803. https://doi.org/10.1016/j.physa.2009.01.026 doi: 10.1016/j.physa.2009.01.026
    [4] Y. Li, H. G. Lee, D. Jeong, J. Kim, An unconditionally stable hybrid numerical method for solving the Allen–Cahn equation, Comput. Math. Appl., 60 (2010), 1591–1606. https://doi.org/10.1016/j.camwa.2010.06.041 doi: 10.1016/j.camwa.2010.06.041
    [5] H. D. Vuijk, J. M. Brader, A. Sharma, Effect of anisotropic diffusion on spinodal decomposition, Soft Matter, 15 (2019), 1319–1326. https://doi.org/10.1039/C8SM02017E doi: 10.1039/C8SM02017E
    [6] T. H. Fan, J. Q. Li, B. Minatovicz, E. Soha, L. Sun, S. Patel, et al., Phase-field modeling of freeze concentration of protein solutions, Polymers, 11 (2018), 10. https://doi.org/10.3390/polym11010010 doi: 10.3390/polym11010010
    [7] R. B. Marimont, M. B. Shapiro, Nearest neighbour searches and the curse of dimensionality, IMA J. Appl. Math., 24 (1979), 59–70. https://doi.org/10.1093/imamat/24.1.59 doi: 10.1093/imamat/24.1.59
    [8] X. Fang, L. Qiao, F. Zhang, F. Sun, Explore deep network for a class of fractional partial differential equations, Chaos Solitons Fractals, 172 (2023), 113528. https://doi.org/10.1016/j.chaos.2023.113528 doi: 10.1016/j.chaos.2023.113528
    [9] V. Charles, J. Aparicio, J. Zhu, The curse of dimensionality of decision-making units: A simple approach to increase the discriminatory power of data envelopment analysis, Eur. J. Oper. Res., 279 (2019), 929–940. https://doi.org/10.1016/j.ejor.2019.06.025 doi: 10.1016/j.ejor.2019.06.025
    [10] V. Berisha, C. Krantsevich, P. R. Hahn, S. Hahn, G. Dasarathy, P. Turaga, et al., Digital medicine and the curse of dimensionality, npj Digit. Med., 4 (2021), 153. https://doi.org/10.1038/s41746-021-00521-5 doi: 10.1038/s41746-021-00521-5
    [11] S. Koohy, G. Yao, K. Rubasinghe, Numerical solutions to low and high-dimensional Allen–Cahn equations using stochastic differential equations and neural networks, Partial Differ. Equations Appl. Math, 7 (2023), 100499. https://doi.org/10.1016/j.padiff.2023.100499 doi: 10.1016/j.padiff.2023.100499
    [12] S. Ham, Y. Hwang, S. Kwak, J. Kim, Unconditionally stable second-order accurate scheme for a parabolic sine–Gordon equation, AIP Adv., 12 (2022), 025203. https://doi.org/10.1063/5.0081229 doi: 10.1063/5.0081229
    [13] S. Ayub, H. Affan, A. Shah, Comparison of operator splitting schemes for the numerical solution of the Allen–Cahn equation, AIP Adv., 9 (2019), 125202. https://doi.org/10.1063/1.5126651 doi: 10.1063/1.5126651
    [14] H. G. Lee, A second-order operator splitting Fourier spectral method for fractional-in-space reaction-diffusion equations, J. Comput. Appl. Math., 333 (2018), 395–403. https://doi.org/10.1016/j.cam.2017.09.007 doi: 10.1016/j.cam.2017.09.007
    [15] D. Jeong, J. Kim, An explicit hybrid finite difference scheme for the Allen–Cahn equation, J. Comput. Appl. Math., 340 (2018), 247–255. https://doi.org/10.1016/j.cam.2018.02.026 doi: 10.1016/j.cam.2018.02.026
    [16] X. Yang, Z. Yang, C. Zhang, Stochastic heat equation: Numerical positivity and almost surely exponential stability, Comput. Math. Appl., 119 (2022), 312–318. https://doi.org/10.1016/j.camwa.2022.05.031 doi: 10.1016/j.camwa.2022.05.031
    [17] Y. Sun, M. Kumar, Numerical solution of high dimensional stationary Fokker–Planck equations via tensor decomposition and Chebyshev spectral differentiation, Comput. Math. Appl., 67 (2014), 1960–1977. https://doi.org/10.1016/j.camwa.2014.04.017 doi: 10.1016/j.camwa.2014.04.017
    [18] S. Shrestha, Monte carlo method to solve diffusion equation and error analysis, J. Nepal Math. Soc., 4 (2021), 54–60. https://doi.org/10.3126/jnms.v4i1.37113 doi: 10.3126/jnms.v4i1.37113
    [19] A. Medved, R. Davis, P. A. Vasquez, Understanding fluid dynamics from Langevin and Fokker–Planck equations, Fluids, 5 (2020), 40. https://doi.org/10.3390/fluids5010040 doi: 10.3390/fluids5010040
    [20] H. Naeimi, F. Kowsary, Finite Volume Monte Carlo (FVMC) method for the analysis of conduction heat transfer, J. Braz. Soc. Mech. Sci. Eng., 41 (2019), 1–10. https://doi.org/10.1007/s40430-019-1762-3 doi: 10.1007/s40430-019-1762-3
    [21] H. Naeimi, Monte carlo methods for heat transfer, Int. J. Math. Game Theory Algebra, 29 (2020), 113–170.
    [22] T. Nakagawa, A. Tanaka, On a Monte Carlo scheme for some linear stochastic partial differential equations, Monte Carlo Methods Appl., 27 (2021), 169–193. https://doi.org/10.1515/mcma-2021-2088
    [23] G. Venkiteswaran, M. Junk, Quasi-Monte Carlo algorithms for diffusion equations in high dimensions, Monte Carlo Methods Appl., 68 (2005), 23–41. https://doi.org/10.1016/j.matcom.2004.09.003 doi: 10.1016/j.matcom.2004.09.003
    [24] A. Novikov, D. Kuzmin, O. Ahmadi, Random walk methods for Monte Carlo simulations of Brownian diffusion on a sphere, Appl. Math. Comput., 364, (2020), 124670. https://doi.org/10.1016/j.amc.2019.124670
    [25] D. Lee, The numerical solutions for the energy-dissipative and mass-conservative Allen–Cahn equation, Comput. Math. Appl., 80 (2020), 263–284. https://doi.org/10.1016/j.camwa.2020.04.007 doi: 10.1016/j.camwa.2020.04.007
    [26] H. G. Lee, J. Y. Lee, A second order operator splitting method for Allen–Cahn type equations with nonlinear source terms, Physica A, 432 (2015), 24–34. https://doi.org/10.1016/j.physa.2015.03.012 doi: 10.1016/j.physa.2015.03.012
    [27] Y. Cheng, A. Kurganov, Z. Qu, T. Tang, Fast and stable explicit operator splitting methods for phase-field models, J. Comput. Phys., 303 (2015), 45–65. https://doi.org/10.1016/j.jcp.2015.09.005 doi: 10.1016/j.jcp.2015.09.005
    [28] A. Chertock, C. R. Doering, E. Kashdan, A. Kurganov, A fast explicit operator splitting method for passive scalar advection, J. Sci. Comput., 45 (2010), 200–214. https://doi.org/10.1007/s10915-010-9381-2 doi: 10.1007/s10915-010-9381-2
    [29] K. Poochinapan, B. Wongsaijai, Numerical analysis for solving Allen–Cahn equation in 1D and 2D based on higher-order compact structure-preserving difference scheme, Appl. Math. Comput., 434 (2022), 127374. https://doi.org/10.1016/j.amc.2022.127374 doi: 10.1016/j.amc.2022.127374
    [30] J. Kai, S. Wei, High-order energy stable schemes of incommensurate phase-field crystal model, Electron. Res. Arch., 28 (2020), 1077–1093. https://doi.org/10.3934/era.2020059 doi: 10.3934/era.2020059
    [31] W. Liupeng, H. Yunqing, Error estimates for second-order SAV finite element method to phase field crystal model, Electron. Res. Arch., 29 (2021), 1735–1752. https://doi.org/10.3934/era.2020089 doi: 10.3934/era.2020089
    [32] J. W. Choi, H. G. Lee, D. Jeong, J. Kim, An unconditionally gradient stable numerical method for solving the Allen–Cahn equation, Physica A, 388 (2009), 1791–1803. https://doi.org/10.1016/j.physa.2009.01.026 doi: 10.1016/j.physa.2009.01.026
    [33] M. Franken, M. Rumpf, B. Wirth, A phase field based PDE constraint optimization approach to time discrete Willmore flow, Int. J. Numer. Anal. Model., 2011 (2011).
    [34] X. Chen, C. M. Elliott, A. Gardiner, J. JING ZHAO, Convergence of numerical solutions to the Allen–Cahn equation, Appl. Anal., 69 (1998), 47–56.
    [35] G. B. Folland, Introduction to Partial Differential Equations, Princeton university press, 1976. https://doi.org/10.1515/9780691213033
    [36] R. A. Fisher, The wave of advance of advantageous genes, Ann. Eugen., 7 (1937), 355–369. https://doi.org/10.1111/j.1469-1809.1937.tb02153.x doi: 10.1111/j.1469-1809.1937.tb02153.x
    [37] A. Kolmogorov, Étude de l'équation de la diffusion avec croissance de la quantité de matière et son application à un problème biologique, Moscow Univ. Bull. Math., 1 (1937), 1–25.
    [38] P. Román-Román, F. Torres-Ruiz, Modelling logistic growth by a new diffusion process: application to biological systems, Biosystems, 110 (2012), 9–21. https://doi.org/10.1016/j.biosystems.2012.06.004 doi: 10.1016/j.biosystems.2012.06.004
    [39] X. Y. Wang, Exact and explicit solitary wave solutions for the generalised fisher equation, Phys. Lett. A, 131 (1988), 277–279. https://doi.org/10.1016/0375-9601(88)90027-8 doi: 10.1016/0375-9601(88)90027-8
  • This article has been cited by:

    1. Engin Can, Hakan Adiguzel, Oscillation analysis of conformable fractional generalized Lienard equations, 2022, 26, 0354-9836, 647, 10.2298/TSCI22S2647C
    2. Petr Hasil, Michal Pospíšil, Jiřina Šišoláková, Michal Veselý, Oscillation criterion for linear equations with coefficients containing powers of natural logarithm, 2024, 203, 0026-9255, 91, 10.1007/s00605-023-01910-6
    3. Jiřina Šišoláková, Non-oscillation of linear differential equations with coefficients containing powers of natural logarithm, 2024, 22, 2391-5455, 10.1515/math-2024-0012
    4. Peng E, Tingting Xu, Linhua Deng, Yulin Shan, Miao Wan, Weihong Zhou, Solutions of a class of higher order variable coefficient homogeneous differential equations, 2025, 20, 1556-1801, 213, 10.3934/nhm.2025011
  • Reader Comments
  • © 2023 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(1470) PDF downloads(64) Cited by(11)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog