An explicit numerical method for the conservative Allen–Cahn equation on a cubic surface

Youngjin Hwang; Jyoti; Soobin Kwak; Hyundong Kim; Junseok Kim; Youngjin Hwang; Jyoti; Soobin Kwak; Hyundong Kim; Junseok Kim

doi:10.3934/math.20241641

AIMS Mathematics

2024, Volume 9, Issue 12: 34447-34465. doi: 10.3934/math.20241641

Previous Article Next Article

Research article Special Issues

An explicit numerical method for the conservative Allen–Cahn equation on a cubic surface

1.
Department of Mathematics, Korea University, Seoul, 02841, Republic of Korea
2.
The Institute of Basic Science, Korea University, Seoul, 02841, Republic of Korea
3.
Department of Mathematics and Physics, Gangneung-Wonju National University, Gangneung 25457, Republic of Korea
4.
Institute for Smart Infrastructure, Gangneung-Wonju National University, Gangneung 25457, Republic of Korea

Received: 08 October 2024 Revised: 27 November 2024 Accepted: 04 December 2024 Published: 06 December 2024
MSC : 80M20, 39A14, 65M06

We introduced a fully explicit finite difference method (FDM) designed for numerically solving the conservative Allen–Cahn equation (CAC) on a cubic surface. In this context, the cubic surface refers to the combined areas of the six square faces that enclose the volume of a cube. The proposed numerical solution approach is structured into two sequential steps. First, the Allen–Cahn (AC) equation was solved by applying the fully explicit FDM, which is computationally efficient. Following this, the conservation term is resolved using the updated solution from the AC equation to ensure consistency with the underlying conservation principles. To evaluate the effectiveness of the proposed scheme, computational tests are performed to verify that the resulting numerical solution of the CAC equation successfully conserves the discrete mass. Additionally, the solution is examined for its ability to exhibit the property of constrained motion by mass conserving mean curvature, a critical characteristic of the CAC equation. These two properties are fundamental to the integrity and accuracy of the CAC equation.

Keywords:

Citation: Youngjin Hwang, Jyoti, Soobin Kwak, Hyundong Kim, Junseok Kim. An explicit numerical method for the conservative Allen–Cahn equation on a cubic surface[J]. AIMS Mathematics, 2024, 9(12): 34447-34465. doi: 10.3934/math.20241641

Related Papers:

[1]	Mansour Shrahili, Mohamed Kayid . Uncertainty quantification based on residual Tsallis entropy of order statistics. AIMS Mathematics, 2024, 9(7): 18712-18731. doi: 10.3934/math.2024910
[2]	Ramy Abdelhamid Aldallal, Haroon M. Barakat, Mohamed Said Mohamed . Exploring weighted Tsallis extropy: Insights and applications to human health. AIMS Mathematics, 2025, 10(2): 2191-2222. doi: 10.3934/math.2025102
[3]	H. M. Barakat, M. A. Alawady, I. A. Husseiny, M. Nagy, A. H. Mansi, M. O. Mohamed . Bivariate Epanechnikov-exponential distribution: statistical properties, reliability measures, and applications to computer science data. AIMS Mathematics, 2024, 9(11): 32299-32327. doi: 10.3934/math.20241550
[4]	Mohamed Said Mohamed, Najwan Alsadat, Oluwafemi Samson Balogun . Continuous Tsallis and Renyi extropy with pharmaceutical market application. AIMS Mathematics, 2023, 8(10): 24176-24195. doi: 10.3934/math.20231233
[5]	Mohamed Said Mohamed, Haroon M. Barakat, Aned Al Mutairi, Manahil SidAhmed Mustafa . Further properties of Tsallis extropy and some of its related measures. AIMS Mathematics, 2023, 8(12): 28219-28245. doi: 10.3934/math.20231445
[6]	Refah Alotaibi, Mazen Nassar, Zareen A. Khan, Wejdan Ali Alajlan, Ahmed Elshahhat . Entropy evaluation in inverse Weibull unified hybrid censored data with application to mechanical components and head-neck cancer patients. AIMS Mathematics, 2025, 10(1): 1085-1115. doi: 10.3934/math.2025052
[7]	Alaa M. Abd El-Latif, Hanan H. Sakr, Mohamed Said Mohamed . Fractional generalized cumulative residual entropy: properties, testing uniformity, and applications to Euro Area daily smoker data. AIMS Mathematics, 2024, 9(7): 18064-18082. doi: 10.3934/math.2024881
[8]	M. Nagy, H. M. Barakat, M. A. Alawady, I. A. Husseiny, A. F. Alrasheedi, T. S. Taher, A. H. Mansi, M. O. Mohamed . Inference and other aspects for $q-$ Weibull distribution via generalized order statistics with applications to medical datasets. AIMS Mathematics, 2024, 9(4): 8311-8338. doi: 10.3934/math.2024404
[9]	Areej M. AL-Zaydi . On concomitants of generalized order statistics arising from bivariate generalized Weibull distribution and its application in estimation. AIMS Mathematics, 2024, 9(8): 22002-22021. doi: 10.3934/math.20241069
[10]	G. M. Mansour, M. A. Abd Elgawad, A. S. Al-Moisheer, H. M. Barakat, M. A. Alawady, I. A. Husseiny, M. O. Mohamed . Bivariate Epanechnikov-Weibull distribution based on Sarmanov copula: properties, simulation, and uncertainty measures with applications. AIMS Mathematics, 2025, 10(5): 12689-12725. doi: 10.3934/math.2025572

Abstract

1. Introduction

Chaos theory is the branch of mathematics that deals with complex behavior emerging in dynamical systems ^[1]. Chaotic systems possess special features such as extreme sensitivity to initial conditions, random-like unpredictable behavior and ergodicity. Improving these characteristics enhance applications in almost every field of science and engineering, such as: Biology ^[2], astrophysics ^[3], mechanics ^[4], economics ^[5], secure communications ^[6], cryptography ^[7], robotics ^[8], control ^[9] and so on. With the study of fractional calculus, it has been found that Fractional-Order (FO) derivatives serve as an extra-degree freedom for enhancing chaos-based applications, this thanks to the FO derivatives that have inheritance or memory properties that preserve complex phenomena with more accuracy ^[10,11]. For instance, in cryptography and secure communications, FO chaotic systems can generate larger keyspace than their integer counterpart ^[12,13].

According to the literature, attractors of dynamical systems can be divided into two categories depending on the localization of their equilibrium points in the corresponding basin of attraction ^[14]. On the one hand, chaotic systems with self-excited attractors are those whose basin of attraction intersects with an unstable equilibrium. On the other hand, chaotic systems with hidden attractors have a basin of attraction that does not intersect with any open neighborhoods of equilibria. Since the first hidden attractor was found in the classical Chua's circuit ^[15], much more attention has been shifted towards this new type of attractors. Firstly, the dynamics of hidden attractors are highly complex to analyze since they can lead to unexpected behavior. For instance, complex phenomena such as "Multistability" ^[16], which refers to when two or more different attractors coexist under the same parameters are likely to appear in dynamical systems with hidden attractors. The multistability phenomena depend on the initial conditions to switch between one attractor to a totally different one, and this property can be used as an additional source of randomness for potential applications. For instance, but not limited to, they might enhance the performance on cryptography ^[17], and secure communications ^[18].

There exist few studies related to FO chaotic systems with hidden attractors, but FO derivatives can generate several families of hidden chaotic attractors in their commensurate and incommensurate models where multistability phenomena can also be observed ^[19,20]. Therefore, the research effort oriented to FO chaotic systems with hidden attractors is an area of opportunity to enhance chaos-based applications. One of the critical factors that guarantees the success of those applications, is the complexity degree of the chaotic attractor that can be quantified by its Kaplan-Yorke dimension ( $D_{KY}$ ). In this manner, the application of metaheuristics, such as: Differential Evolution (DE) and Accelerated Particle Swarm Optimization (APSO) algorithms, has gained interest in recent years for solving optimization problems of FO chaotic systems ^[21]. In ^[22] the optimization through metaheuristics of the $D_{KY}$ is achieved on four FO chaotic systems with self-excited attractors, a similar approach is presented in ^[23]. To the best of the authors' knowledge, there are no works that explore the optimization of FO chaotic systems with hidden attractors.

Motivated by the above discussion, this paper aims to introduce the optimization of FO chaotic systems with hidden attractors. One of the challenges in the optimization process, is that hidden attractors have small basin of attraction ^[24,25,26], and thus, tiny variations in the parameters, and fractional-order derivatives or initial conditions can easily shift the chaotic behavior to a periodic dynamics. Hence, careful considerations needs to be taken to ensure that the optimization process is performing correctly. The literature reports several methods for detecting Chaos in dynamical systems, such as phase portraits, bifurcation diagrams, equilibrium points analysis through Shilnikov theorem, Lyapunov exponents (LEs), Poincaré maps, among others. It is well-known that LEs provide quantitative information over the stretching and folding of the chaotic attractors. However, accurate computation of Lyapunov exponents requires the analysis of large time series by specialized library packages such as TISEAN ^[27], which can be high computational expensive, but LEs are necessary for calculating the $D_{KY}$ of a chaotic attractor. For such problem, one can introduce constraints that ensures the chaotic motion for its LEs computation through TISEAN. In this manner, Poincaré maps provide a qualitative understating of the dynamical characteristics of a $n$ -dimensional system into a $(n-1)$ representation ^[28]. This is done by the intersection analysis of a hyper-plane transverse to the flow of the dynamical system, then the Poincaré map corresponding to a chaotic system should display enough distribution over a 2-D plane. In this regard, the contribution of this work is to detect chaotic motion in FO chaotic systems with hidden attractors through the distribution quantification of Poincaré maps. Moreover, to save execution time, the optimization algorithms are implemented by using the Numba just-in-time (jit) compiler in the Python programming language.

The rest of the paper is organized as follows: Section 2 summarizes the background of this work. DE and APSO algorithms are described in Section 3. Section 4 shows the numerical analysis of the FO chaotic systems. Section 5 details the optimization process whereas the discussion of the optimization results is presented in Section 6. Finally, the authors' conclusions are given in Section 7.

2. Mathematical background

2.1. Grünwald-Letnikov method

Although one can cite various definitions of FO derivatives such as as the Grünwald-Letnikov ^[29], Riemman Liouville ^[30] and Caputo definitions ^[31], these mathematical interpretations do not allow a direct understanding of the derivative solution. Therefore, the solution of FO differential equations are often obtained through numerical approximations accordingly the used fractional derivative definition ^{[32,33,34,35]}. For instance, the Adams Bashforth Moulton (ABM) method is a predictor-corrector scheme highly recognized as a valuable fractional calculus tool due to the good approximation that it produces ^[36]. Drawbacks of this scheme are related to the complexity for representing the memory properties of FO derivatives, that requires from long time computation. The Grünwald Letnikov (GL) method is another numerical scheme very suitable for solving FO systems. It can be considered as a generalization of the Euler method for solving ordinary differential equations ^[37]. The GL method is of low computational cost and has approximately the same order of accuracy and good match of numerical solutions as ABM ^[38]. Therefore, this work uses the Grünwald-Letnikov method for solving the FO differential equations that comprises the FO chaotic systems that are the case study. The mathematical description is given below.

According to the explicit numerical approximation of the Grünwald-Letnikov derivative ^[39], the FO derivative $0 < q < 1$ for a discrete function $f(t_k)$ , can be described by (2.1), $h$ is the integration step-size, and $C^q_j$ are binomial coefficients which are recursively calculated by (2.2).

$\begin{equation} D^q_{t_k}f(t_k)\approx h^{-q}\sum\limits_{j = 0}^kC^q_jf(t_{k-j}), \end{equation}$

(2.1)

$\begin{equation} C^q_0 = 1,\;C^q_j = \bigg( 1 - \frac{1+q}{j}\bigg)C^q_{j-1},\;j = 1,2,...,k. \end{equation}$

(2.2)

The solution of a FO system in the form $D^q_{t_k}x(t_k) = f(x(t_k))$ , can be obtained by (2.3), which ideally requires infinite memory length for numerical simulation since the summation term depends on the discretized time $t_k$ .

$\begin{equation} x(t_k) = f(x(t_{k-1}))h^q - \sum\limits_{j = 1}^kC^q_jx(t_{k-j}). \end{equation}$

(2.3)

2.2. Poincaré maps

One of the features of chaos behavior is the existence of a dense set of periodic orbits, which implies that any periodic trajectory of the orbit visits an arbitrarily small neighborhood of a non-periodic one ^[40]. In this manner, and since Poincaré maps preserve many of the orbits of dynamical systems, it is a viable tool for detecting chaos in both integer and FO systems ^[41,42]. Below is described the elaboration of a Poincaré map.

Lets consider a FO system of the form given in (2.4), where $0 < q < 1$ , $x(t) = [x_1(t), x_2(t), ..., x_n(t)]^T\in\mathbb{R}^n$ and $f\big(x(t)\big) = [f_1\big(x(t)\big), f_2\big(x(t)\big), ..., f_n\big(x(t)\big)]^T$ .

$\begin{equation} D^q_{t}x(t) = f\big(x(t)\big), \end{equation}$

(2.4)

The Poincaré map of the system given in (2.4) is obtained as follows:

$(1)$ Without loss of generality, let $\Sigma\in\mathbb{R}^{n-1}$ be the Poincaré section of system (2.4) such that $\Sigma\cap x(t) = [x_1(t), x_2(t), ..., x_{n-1}(t)]^T$ .

$(2)$ Next, the Poincaré section is settled on a given $x_n = \sigma$ .

$(3)$ Finally, all states $[x_1(t), x_2(t), ..., x_{n-1}(t)]^T$ are captured each time the dynamical flow intersects the Poincaré section.

The captured points by the above procedure portrait the Poincaré map of system (2.4), then, the behavior of the system is obtained by analyzing the distribution of the distinct points of the Poincaré map. If the Poincaré map consists neither of a small number of points or filling a continuous line but instead a large and irregular distribution, then it is a strong indicator of chaos ^[43].

2.3. Kaplan Yorke dimension

The Kaplan Yorke dimension $D_{KY}$ describes the fractal dimension of an attractor by means of its Lyapunov exponents as follows ^[44]:

$\begin{equation} D_{KY} = j + \frac{1}{|\lambda_{j+1}|} \sum\limits_{k = 1}^j \lambda_k, \end{equation}$

(2.5)

herein, the Lyapunov exponents are sorted from highest to lowest $\lambda_1\geq\lambda_2\geq\cdots\geq\lambda_{j+1}$ and $j\in\mathbb{N}$ is the largest integer for which $\lambda_1 + \lambda_2 +\cdots + \lambda_j\geq0$ . The $D_{KY}$ gives a meaningful measure to quantify the complexity of a chaotic attractor ^[45,46]. Therefore, it is the interest of this optimization work to determine a set of parameters and fractional orders that leads to a higher $D_{KY}$ .

3. DE and APSO algorithms

3.1. Differential Evolution algorithm

The Differential Evolution (DE) algorithm belongs to the family of algorithms based on biological evolution, which uses genetic cross and mutation operators to generate new candidate solutions ^[47]. The structure of DE is illustrated in the chart flow of Figure 1, whose description is as follows:

Figure 1. Chart flow of the Differential Evolution algorithm.

DownLoad: Full-Size Img PowerPoint

$(1)$ In the first step, the initial population is generated within the range of interest.

$\begin{equation} X_i\in\mathbb{R}^n,\;i = 1,2,...,p, \end{equation}$

(3.1)

$(2)$ Next, the objective function is evaluated for each individual in the population

$\begin{equation} \forall X_i,\; do\; f(X_i)\in\mathbb{R}. \end{equation}$

(3.2)

$(3)$ In the next step, the mutation operation is carried out:

● (3.1) First, three randomly individuals $X_1$ , $X_2$ and $X_3$ are selected.

● (3.2) Then, the following operation is performed:

$\begin{equation} w_i = X_1 + \mu (X_2 - X_3) \end{equation}$

(3.3)

where $\mu\in\mathbb{R}$ is a scaling parameter.

$(4)$ In the next step, individuals are crossed according to the condition

$\begin{equation} F_i(j) = \left\{ \begin{array}{ll} w_i(j), & \mathrm{\ } {\rm{rand}} \leq Cr, \\ X_i(j), & \; {\rm{other}}. \\ \end{array} \right. \end{equation}$

(3.4)

Specifically, if the random number (rand) is less than or equal to the cross parameter $Cr$ , then the term $F_i (j)$ takes the value $w_i (j)$ , otherwise it will be $X_i (j)$ . Next, the first generation of individuals generates new individuals through mutation and cross operations. Subsequently, the objective function is evaluated in the new individual.

$(5)$ Steps 2–5 are repeated for a number of $k$ generations.

$(6)$ In the last step, the lowest cost function is printed.

3.2. Accelerated Particle Swarm Optimization algorithm

In the Accelerated Particle Swarm Optimization (APSO) algorithm, the tentative solutions are called particles. This algorithm consists of disturbing the behavior of the particles in such a way that they can move within a search space and finding the solution to the optimization problem ^[48]. The algorithm is inspired by the flight of birds to find food in an unknown field with the lowest possible energy expenditure, which can be translated into finding the minimum value of a function. The structure of the APSO is illustrate in the chart flow of Figure 2, and its description is given as follows:

Figure 2. Chart flow of the Accelerated Particle Swarm Optimization algorithm.

DownLoad: Full-Size Img PowerPoint

$(1)$ In the first step, the initial population is generated, meaning the particles position within a range of interest.

$\begin{equation} X_i\in\mathbb{R}^n,\;i = 1,2,...,p. \end{equation}$

(3.5)

$(2)$ Next, the objective function is evaluated for each particle

$\begin{equation} \forall X_i,\; do\; f(X_i)\in \mathbb{R}. \end{equation}$

(3.6)

$(3)$ In the next step, the position of the particles is disturbed:

● (3.1) First, the best particles are selected

$\begin{equation} g^* = \min\big(f(X_i)\big). \end{equation}$

(3.7)

● (3.2) Then, the position is updated

$\begin{equation} X_i (k+1) = (1-\beta) + \beta g^* + \alpha\epsilon_n . \end{equation}$

(3.8)

Specifically, the new population is generated from the change in position of the particles. The value for $\epsilon_n$ is randomly initialized within the range $[0, 1]$ whereas $\alpha$ and $\beta$ are estimated experimentally.

$(4)$ Steps 2–5 are repeated for a number of $k$ generations.

$(5)$ In the last step, the lowest cost function is printed.

4. Fractional order chaotic systems

This section shows the numerical simulations and Poincaré maps of the FO chaotic systems that are case study. For instance, the authors in ^[19] introduced the 3D FO chaotic system given in (4.1), where $0 < q_i < 1$ , $(i = 1, 2, 3)$ and $a\in\mathbb{R}$ . According to expression (4.2), it can be seen that by setting the parameter $a > 1/4$ , the 3D FO chaotic system does not posses an equilibrium point ( $E_p = (\bar{x}, \bar{y}, \bar{z})$ ), and thus, any existing basin of attraction is considered as hidden.

$\begin{equation} \begin{split} D^{q_1}_{t}x & = yz + x(y-a), \\ D^{q_2}_{t}y & = 1 - |x|,\\ D^{q_3}_{t}z & = -xy-z,\\ \end{split} \end{equation}$

(4.1)

$\begin{equation} \begin{split} E_{p_1}& = \big( 1, (1+\sqrt{1-4a})/2, -(1+\sqrt{1-4a})/2 \big), \\ E_{p_2}& = \big( -1, (1-\sqrt{1-4a})/2, (1-\sqrt{1-4a})/2 \big), \\ E_{p_3}& = \big( 1, (1-\sqrt{1-4a})/2, -(1-\sqrt{1-4a})/2 \big), \\ E_{p_4}& = \big( -1, (1+\sqrt{1-4a})/2, (1-\sqrt{1-4a})/2 \big). \\ \\ \end{split} \end{equation}$

(4.2)

By taking into account (2.3), the solution of the 3D FO chaotic system by means of the Grünwald-Letnikov method is given in (4.3), whereas Figure 3 shows the simulation results. Figure 4 shows the LEs spectra by varying the fractional derivative.

$\begin{equation} \begin{split} x(t_k) & = \bigg(y(t_{k-1})z(t_{k-1}) + x(t_{k-1})\big(y(t_{k-1})-a\big)\bigg)h^{q_1} - \sum\limits_{j = 1}^kC^{q_1}_jx(t_{k-j}),\\ y(t_k)& = \big( 1 - |x(t_{k-1})|\big)h^{q_2} - \sum\limits_{j = 1}^kC^{q_2}_jy(t_{k-j}),\\ z(t_k)& = \big(-x(t_{k-1})y(t_{k-1})-z(t_{k-1})\big)h^{q_3} - \sum\limits_{j = 1}^kC^{q_3}_jz(t_{k-j}). \end{split} \end{equation}$

(4.3)

Figure 3. Numerical simulation of the 3D FO chaotic system with hidden attractors, with step size

$h = 0.015$ , simulation time

$T_{sim} = 2200s$ , initial conditions

$(x(t_0), y(t_0), z(t_0)) = (1, 1, 1)$ , parameter

$a = 0.35$ , and fractional order

$q_i = 0.97$ ,

$(i = 1, 2, 3)$ . Hidden chaotic attractor depicted in (a)

$xy$ phase portrait, (b)

$xz$ phase portrait, (c)

$yz$ phase portrait and (d) the corresponding time series.

DownLoad: Full-Size Img PowerPoint

Figure 4. LEs spectra of the 3D FO system with hidden attractors, with

$a = 0.35$ and varying the fractional-order derivative

$q$ .

$\lambda_1$ is in blue color,

$\lambda_2$ is in orange color and

$\lambda_3$ is in green color.

DownLoad: Full-Size Img PowerPoint

To demonstrate the usefulness of Poincaré maps in detecting chaotic behavior, a Poincaré section $\Sigma\in\mathbb{R}^2$ on $y = 2$ is defined as depicted in , then, by following the steps of Section 2.2, the Poincaré map of is obtained. It can be seen that the displayed points resembles the $xz$ phase portrait of , the markers illustrates if the flow intersects the Poincaré section on $y > 2$ (blue, in-direction) and if it intersects on $y < 2$ (red, out-direction).

Figure 5. (a) Chaotic flow of the 3D FO chaotic system intersecting the Poincaré section

$\Sigma\in\mathbb{R}^2$ on

$y = 2$ , and (b) its corresponding Poincaré map.

DownLoad: Full-Size Img PowerPoint

It is worthy to mention that the number of points portrayed in the Poincaré map is in relation to the amount of used data. However, the number of points itself says nothing about the dynamical behavior. For instance, lets us consider the dynamical behavior of the system (4.1) illustrated in , then, a Poincaré section $\Sigma\in\mathbb{R}^2$ on $y = 2$ is defined as depicted in Figure 7 (a). The obtained Poincaré map in Figure 7 (b) displays only two points, indicating the flow intersects twice in different in-out directions of the Poincaré section, and thus, resembling a periodic behavior. Nevertheless, this is only apparent, due to the fact of the non-existence of exactly periodic solutions in FO systems ^[49,50], thereby, the importance of the distribution quantization of the Poincaré map. The next step is the distribution quantization of the Poincaré map. Herein, a simple procedure is proposed as follows:

Figure 6. Numerical simulation of the 3D FO system (4.1) with hidden attractors, with step size

$h = 0.015$ , simulation time

$T_{sim} = 2200s$ , initial conditions

$(x(t_0), y(t_0), z(t_0)) = (1, 0, 1)$ , parameter

$a = 0.35$ , and fractional order

$q_i = 0.996$ ,

$(i = 1, 2, 3)$ . Hidden chaotic attractor depicted in (a)

$xy$ phase portrait, (b)

$xz$ phase portrait, (c)

$yz$ phase portrait and (d) the corresponding time series.

DownLoad: Full-Size Img PowerPoint

Figure 7. (a) Quasi-periodic flow of the 3D FO system intersecting the Poincaré section

$\Sigma\in\mathbb{R}^2$ on

$y = 2$ , and (b) its corresponding Poincaré map.

DownLoad: Full-Size Img PowerPoint

$(1)$ Verify that the solution of a system (4.1) remains bounded but does not converge to a fixed point.

$(2)$ Next, the first 10 % samples of the time series are deleted to avoid the initial transient.

$(3)$ Then, the Poincaré map is portrayed with the obtained points when the flow intersects on $y < 2$ (out-direction).

$(4)$ In the next step, only the points that satisfies the condition (4.4) are considered, where $(u_j, v_j)$ and $(u_k, v_k)$ denotes the coordinates of the $P_j$ and $P_k$ points respectively, $p$ is the total number of points intersecting the Poincaré section and $\epsilon = 0.1$ , this ensures to preserve only the points distributed over the plane and remove those located closer than $(\epsilon, \epsilon)$ ,

$\begin{equation} |P_j(u_j,v_j)-P_k(u_k,v_k)|\geq(\epsilon,\epsilon),\;\;k = j+1,\cdots,p,\;\forall j = 1,2,\cdots p-1 . \end{equation}$

(4.4)

$(5)$ Finally, the number of points remaining in the Poincaré map is co-related with the dynamical behavior.

For instance, Figures 8 and illustrates the hidden attractors and its corresponding Poincaré maps when system (4.1) behaves quasi-periodic. After applying 1–5 from above procedure, the obtained results in shows that the remaining points are 10 and 3 respectively. On the other hand, when system (4.1) behaves chaotic such as in , the remaining points in the Poincaré map must portrait $\approx 100$ points.

Figure 8. (a) Quasi-periodic flow of the 3D FO sytem (4.1) with configuration

$a = 0.436$ ,

$(q_1, q_2, q_3) = (0.998, 0.989, 0.986)$ and initial conditions

$(x(t_0), y(t_0), z(t_0)) = (1, 0, 1)$ , intersecting the Poincaré section

$\Sigma\in\mathbb{R}^2$ on

$y = 2$ , (b) shows its corresponding Poincaré map.

DownLoad: Full-Size Img PowerPoint

Figure 9. (a) Quasi-periodic flow of the 3D FO system (4.1) with configuration

$a = 0.277$ ,

$(q_1, q_2, q_3) = (0.988, 0.948, 0.881)$ and initial conditions

$(x(t_0), y(t_0), z(t_0)) = (1, 0, 1)$ , intersecting the Poincaré section

$\Sigma\in\mathbb{R}^2$ on

$y = 2$ , (b) shows its corresponding Poincaré map.

DownLoad: Full-Size Img PowerPoint

Table 1. Relation between the remainder number of points in the Poincaré map of the 3D FO system (4.1) and its dynamical behavior. Obtained results with step size

$h = 0.025$ , simulation time

$T_{sim} = 2200s$ , and initial conditions

$(x(t_0), y(t_0), z(t_0)) = (1, 0, 1)$ .

Parameters	Fractional order derivative	Number of points	LEs	$D_{KY}$	Dynamical behavior
$a = 0.35$	$q_1 = 0.97$ , $q_2 = 0.97$ , $q_3 = 0.97$	107	$\lambda_1 = 0.0011$ , $\lambda_2 = -0.0006$ , $\lambda_3 = -0.0036$	2.1504	chaotic
$a = 0.3$	$q_1 = 0.993$ , $q_2 = 1$ , $q_3 = 0.75$	96	$\lambda_1 = 0.0010$ , $\lambda_2 = -0.0009$ , $\lambda_3 = -0.0047$	2.0208	chaotic
$a = 0.325$	$q_1 = 1$ , $q_2 = 0.99$ , $q_3 = 0.85$	104	$\lambda_1 = 0.0018$ , $\lambda_2 = -0.0006$ , $\lambda_3 = -0.0057$	2.2071	chaotic
$a = 0.337$	$q_1 = 0.966$ , $q_2 = 0.968$ , $q_3 = 0.926$	105	$\lambda_1 = 0.0013$ , $\lambda_2 = -0.0008$ , $\lambda_3 = -0.0036$	2.1452	chaotic
$a = 0.6$	$q_1 = 0.992$ , $q_2 = 0.938$ , $q_3 = 0.968$	94	$\lambda_1 = 0.0022,$ $\lambda_2 = -0.0009,$ $\lambda_3 = -0.0051$	2.2452	chaotic
$a = 0.437$	$q_1 = 0.993$ , $q_2 = 0.956$ , $q_3 = 0.874$	105	$\lambda_1 = 0.0010$ , $\lambda_2 = -0.0009$ , $\lambda_3 = -0.0043$	2.003	chaotic
$a = 0.322$	$q_1 = 0.982$ , $q_2 = 0.955$ , $q_3 = 0.985$	96	$\lambda_1 = 0.0017$ , $\lambda_2 = -0.0009$ , $\lambda_3 = -0.0052$	2.1593	chaotic
$a = 0.29$	$q_1 = 0.986$ , $q_2 = 0.967$ , $q_3 = 0.981$	17	$\lambda_1 = 0.0008$ , $\lambda_2 = -0.0017$ , $\lambda_3 = -0.0018$	1.5219	quasi-periodic
$a = 0.436$	$q_1 = 0.998$ , $q_2 = 0.989$ , $q_3 = 0.986$	10	$\lambda_1 = 0.0007$ , $\lambda_2 = -0.0012$ , $\lambda_3 = -0.0028$	1.8214	quasi-periodic
$a = 0.277$	$q_1 = 0.988$ , $q_2 = 0.948$ , $q_3 = 0.881$	3	$\lambda_1 = 0.0009$ , $\lambda_2 = -0.0038$ , $\lambda_3 = -0.0183$	1.8415	quasi-periodic
$a = 0.382$	$q_1 = 0.994$ , $q_2 = 0.922$ , $q_3 = 97$	2	$\lambda_1 = -0.0003$ , $\lambda_2 = -0.0011$ , $\lambda_3 = -0.0015$	1.06	quasi-periodic
$a = 0.387$	$q_1 = 0.976$ , $q_2 = 0.968$ , $q_3 = 0.804$	2	$\lambda_1 = 0.0001$ , $\lambda_2 = -0.0011$ , $\lambda_3 = -0.0016$	1.375	quasi-periodic

| Show Table

DownLoad: CSV

The second case study is a 4D FO chaotic system introduced in ^[51] as a 4D integer chaotic system with coexisting hidden attractors. Later on, in ^[52] the FO version given in (4.5) was proposed, where $0 < q_i < 1$ , $(i = 1, 2, 3)$ and $a\in\mathbb{R}$ . According to expression (4.6), it can be seen that for $a\neq 0$ , the 4D FO chaotic system does not present an equilibrium point ( $E_p = (\bar{x}, \bar{y}, \bar{z}, \bar{w})$ ), and thus, any existing basin of attraction is considered as hidden.

$\begin{equation} \begin{split} D^{q_1}_{t}x & = y - x, \\ D^{q_2}_{t}y & = -xz + w,\\ D^{q_3}_{t}z & = xy - a,\\ D^{q_4}_{t}z & = -by,\\ \end{split} \end{equation}$

(4.5)

$\begin{equation} \begin{split} y - x & = 0, \\ -xz + w & = 0,\\ xy - a & = 0,\\ -by & = 0.\\ \end{split} \end{equation}$

(4.6)

By taking into account (2.3), the solution of the 4D FO chaotic system by means of the Grünwald-Letnikov method is given in (4.7), whereas Figure 10 shows the simulation results. Figure 11 shows the LEs spectra by varying the fractional derivative.

$\begin{equation} \begin{split} x(t_k) & = \big(y(t_{k-1})- x(t_{k-1})\big)h^{q_1} - \sum\limits_{j = 1}^kC^{q_1}_jx(t_{k-j}),\\ y(t_k)& = \big(-x(t_{k-j})z(t_{k-j}) + w(t_{k-j}) \big)h^{q_2} - \sum\limits_{j = 1}^kC^{q_2}_jy(t_{k-j}),\\ z(t_k)& = \big(x(t_{k-1})y(t_{k-1})-a)\big)h^{q_3} - \sum\limits_{j = 1}^kC^{q_3}_jz(t_{k-j}), \\ z(t_k)& = \big(-by(t_{k-1})\big)h^{q_3} - \sum\limits_{j = 1}^kC^{q_4}_jw(t_{k-j}). \\ \end{split} \end{equation}$

(4.7)

Figure 10. Numerical simulation of the 4D FO chaotic system with hidden attractors, with step size

$h = 0.003$ , simulation time

$T_{sim} = 400s$ , initial conditions

$(x(t_0), y(t_0), z(t_0), w(t_0)) = (0.1, 5, 0, 2)$ , parameters

$a = 3.5$ ,

$b = 0.1$ and fractional order

$q_i = 0.98$ ,

$(i = 1, \cdots, 4)$ . Hidden chaotic attractor depicted in (a)

$xy$ phase portrait, (b)

$xz$ phase portrait, (c)

$yz$ phase portrait, (d)

$xw$ phase portrait and (e) the corresponding time series.

DownLoad: Full-Size Img PowerPoint

Figure 11. LEs spectra of the 4D FO system with hidden attractors (4.5), with

$a = 0.35$ and varying the fractional-order derivative

$q$ .

$\lambda_1$ is in blue color,

$\lambda_2$ is in orange color,

$\lambda_3$ is in green color and

$\lambda_4$ is in red color.

DownLoad: Full-Size Img PowerPoint

In the same manner, a Poincaré map is elaborated for analyzing the behavior of the 4D FO chaotic system. Herein, the Poincaré section $\Sigma\in\mathbb{R}^2$ is settled on $z = 0$ as depicted in Figure 12 (a), then, by following the steps of Section 2.2, the Poincaré map of Figure 12 (b) is obtained.

Figure 12. (a) Chaotic flow of the 4D FO chaotic system intersecting the Poincaré section

$\Sigma\in\mathbb{R}^2$ on

$z = 2$ , and (b) its corresponding Poincaré map.

DownLoad: Full-Size Img PowerPoint

The distribution quantization of the Poincaré map of the 4D FO chaotic system is quite similar to the quantization of the Poincaré map of the 3D FO chaotic system. Herein, only the obtained points when the flow intersects on $z > 0$ (in-direction) are considered with a separation of $\epsilon = 0.1$ . This procedure is repeated for different parameters and fractional order for the system (4.5). The obtained results with a number of data of $1.20\times 10^{5}$ are shown in . It can be seen that for chaotic behavior to exist in the system (4.1), its corresponding Poincaré map must portrait $\approx 40-60$ points.

Table 2. Relation between the remainder number of points in the Poincaré map of the 4D FO system (4.1) and its dynamical behavior. Obtained results with step size

$h = 0.003$ , simulation time

$T_{sim} = 400s$ , and initial conditions

$(x(t_0), y(t_0), z(t_0), w(t_0)) = (0.1, 5, 0, 2)$ .

Parameters	Fractional order derivative	Number of points	LEs	$D_{KY}$	Dynamical behavior
$a = 3.5$ , $b = 0.1$	$q_1 = 0.98$ , $q_2 = 0.98$ , $q_3 = 0.98$ , $q_4 = 0.98$	43	$\lambda_1 = 0.0011$ , $\lambda_2 = 0$ , $\lambda_3 = -0.0001$ , $\lambda_4 = -0.0024$	3.4257	chaotic
$a = 3.5$ , $b = 0.1$	$q_1 = 1.0$ , $q_2 = 0.997$ , $q_3 = 1.0$ , $q_4 = 0.45$	50	$\lambda_1 = 0.0030$ , $\lambda_2 = 0$ , $\lambda_3 = -0.0017$ , $\lambda_4 = -0.0040$	3.3368	chaotic
$a = 2.9$ , $b = 0.1$	$q_1 = 0.45$ , $q_2 = 0.95$ , $q_3 = 1.0$ , $q_4 = 0.985$	40	$\lambda_1 = 0.0008$ , $\lambda_2 = 0$ , $\lambda_3 = -0.0003$ , $\lambda_4 = -0.0030$	3.1666	chaotic
$a = 3.162$ , $b = 0.182$	$q_1 = 0.980$ , $q_2 = 0.942$ , $q_3 = 0.978$ , $q_4 = 0.975$	42	$\lambda_1 = 0.0008$ , $\lambda_2 = 0.0001$ , $\lambda_3 = 0$ , $\lambda_4 = -0.0018$	3.5029	hyper-chaotic
$a = 3.128$ , $b = 0.285$	$q_1 = 0.946$ , $q_2 = 0.964$ , $q_3 = 0.957$ , $q_4 = 0.923$	41	$\lambda_1 = 0.0003$ , $\lambda_2 = 0.0001$ , $\lambda_3 = 0$ , $\lambda_4 = -0.0006$	3.5103	hyper-chaotic
$a = 3.131$ , $b = 0.285$	$q_1 = 0.945$ , $q_2 = 0.963$ , $q_3 = 0.956$ , $q_4 = 0.922$	41	$\lambda_1 = 0.0004$ , $\lambda_2 = 0.0002$ , $\lambda_3 = 0$ , $\lambda_4 = -0.0011$	3.5282	hyper-chaotic
$a = 3.088$ , $b = 0.12$	$q_1 = 0.925$ , $q_2 = 0.948$ , $q_3 = 0.933$ , $q_4 = 0.932$	38	$\lambda_1 = 0.0016$ , $\lambda_2 = 0$ , $\lambda_3 = -0.0004$ , $\lambda_4 = -0.0031$	3.387	chaotic
$a = 3.601$ , $b = 0.3$	$q_1 = 0.999$ , $q_2 = 0.999$ , $q_3 = 0.999$ , $q_4 = 0.962$	48	$\lambda_1 = 0.0005$ , $\lambda_2 = 0$ , $\lambda_3 = -0.0001$ , $\lambda_4 = -0.0022$	3.2014	chaotic
$a = 9.4$ , $b = 0.01$	$q_1 = 0.936$ , $q_2 = 0.935$ , $q_3 = 0.958$ , $q_4 = 0.946$	60	$\lambda_1 = 0.0009$ , $\lambda_2 = 0$ , $\lambda_3 = -0.0001$ , $\lambda_4 = -0.0013$	3.5972	chaotic
$a = 8.4$ , $b = 0.1$	$q_1 = 0.990$ , $q_2 = 0.980$ , $q_3 = 0.980$ , $q_4 = 0.981$	8	$\lambda_1 = 0$ , $\lambda_2 = -0.0002$ , $\lambda_3 = -0.0014$ , $\lambda_4 = -0.0039$	2.5897	quasi-periodic
$a = 7.8$ , $b = 0.1$	$q_1 = 0.990$ , $q_2 = 0.980$ , $q_3 = 0.970$ , $q_4 = 0.981$	10	$\lambda_1 = 0$ , $\lambda_2 = -0.0002$ , $\lambda_3 = -0.0011$ , $\lambda_4 = -0.0041$	2.6829	quasi-periodic
$a = 7.3$ , $b = 0.08$	$q_1 = 0.989$ , $q_2 = 0.998$ , $q_3 = 0.967$ , $q_4 = 0.981$	11	$\lambda_1 = -0.0001$ , $\lambda_2 = -0.0002$ , $\lambda_3 = -0.0007$ , $\lambda_4 = -0.0042$	2.7619	quasi-periodic

| Show Table

DownLoad: CSV

5. Description of the optimization process

This section details the maximization process of the $D_{KY}$ for the 3D and 4D FO chaotic systems. Given the nature of the DE and APSO, the optimization problem is given as follows:

$\begin{equation} {\rm{Maximize}}\big(f(X)\big) = {\rm{Minimize}} \big(-f(X)\big), \;f(X)\in\mathbb{R}, \end{equation}$

(5.1)

where $X\in\mathbb{R}^n$ are the parameters and fractional order derivative values of the FO chaotic system of interest and $f(X)$ returns the $D_{KY}$ .

The optimization process is performed applying DE and APSO who are programmed in Python, because it is a multi-paradigm programming language highly versatile for creating all kind of applications and, therefore, it has become one of the most used in recent years. However, Python based applications pays a penalty in the form of speed execution, thereby complex projects can benefit from the use of the Numba JIT compiler ^[53], which allows to translate Python functions to optimized machine code and thus obtaining better speed performance ^[54]. Moreover, the implementation of the Numba JIT compiler only requires from slight modifications in the original Python code and, therefore, maintaining the flexibility of the language. shows a comparison of the speed performance between a normal Python code and one Python code using the Numba JIT compiler. In both executions, the argument is to solve the proposed FO chaotic system and elaborate its corresponding Poincaré maps for the chaotic behavior validation. It can be seen, that the optimized Python code performs $\approx 375$ times faster than the non-optimized. Thereby the interest of using this compiler in the optimization process.

Table 3. Execution time in seconds, between a Python code and a Python code using the Numba Jit compiler.

FO chaotic system	Python code	Code using the Numba Jit compiler	Number of generated data
3D	$1.025\times 10^{3}s$	$2.61s$	$146,666$
4D	$2.085\times 10^{3}s$	$5.55s$	$133,332$

| Show Table

DownLoad: CSV

Algorithms 1 and 2 show the pseudocodes for the implementation of DE and APSO, respectively. It can be seen that the implementation is quite similar in both algorithms. The Numba JIT compiler is used in most of the Python code with the exception of the Lyapunov exponents (LEs) computation $D_{KY}$ computation, in this step, the computation of the Lyapunov exponents $\lambda_i$ , $(i = 1, 2, ..., n)$ is made with the TISEAN package by removing the first $10\%$ samples of the time series.

Algorithm 1 Maximizing

$D_{KY}$ applying DE

Require:

$\mu$ ,

$C_r$ ,

$p$ ,

$gen$ .
Ensure: best

$D_j$ ,

$j=0, 2, ..., gen$
Initialize population

$X_i\in\mathbb{R}^n, \; i=1, 2, ..., p$
for

$i=1:p$ do
Compute numerical solution for each

$X_i$ using (2.3)
Calculate

$D_{KY_i}$ using (2.5)
store

$D_0 \gets$ best

$D_{KY_i}$
for

$j=1:gen$ do
for

$k=1:p$ do
generate

$X_k$ using (3.3) and (3.4)
Compute numerical solution for each new

$X_i$ using (2.3)
    Detect chaotic behavior using Poincaré map
   Compute LEs using Tisean
   Calculate

$D_{KY_k}$ using (2.5)
store

$D_j \gets$ best

$D_{KY_k}$

| Show Table

DownLoad: CSV

Algorithm 2 Maximizing

$D_{KY}$ applying APSO

Require:

$\alpha$ ,

$\beta$ ,

$p$ ,

$gen$ .
Ensure: Best

$D_j$ , for

$j=0, 2, \dots, gen$ .
Initialize population

$X_i\in\mathbb{R}^n, \; i=1, 2, ..., p$
for

$i=1:p$ do
Compute numerical solution for each

$X_i$ using (2.3)
Calculate

$D_{KY_i}$ using (2.5)
Store

$D_0 \gets$ best

$D_{KY_i}$
for

$j=1:gen$ do
for

$k=1:p$ do
Generate

$X_k$ using (3.8)
Compute numerical solution for each new

$X_i$ using (2.3)
    Detect chaotic behavior using Poincaré map
    Compute LEs using Tisean
    Calculate

$D_{KY_k}$ using (2.5)
Store

$D_j \gets$ best

$D_{KY_k}$

| Show Table

DownLoad: CSV

6. Optimization results

The DE and APSO were executed for a 100 population and for 50–100 generations. shows the ranges of interest and the configuration of the metaheuristics for the optimization of the 3D and 4D FO chaotic systems, respectively. The values of the parameters and fractional-order derivatives in the search space, were encoded within three decimal places. For convenience, the initial population is set with parameters and fractional-order derivatives that lead to chaotic behavior and, therefore, the elaboration of Poincaré maps is not necessary for this step. In and are shown the optimization results of the 3D and 4D FO chaotic systems by DE and APSO, respectively. It can be seen that in all cases, the obtained results lead to higher $D_{KY}$ than those shown in Tables 1 and 2, thus validating the optimization process. The best results for the 3D FO chaotic system are those obtained APSO by using 50 generations whereas the best results for the 4D FO chaotic system are obtained with DE by using 100 generations.

Table 4. Range of parameters and fractional order derivative values for the optimization of the 3D and 4D FO chaotic systems.

FO chaotic system	Parameters	FO derivative	DE parameters	APSO parameters
3D	$a\in[0.251, 0.5]$	$q_i\in[0.3, 0.999]$	$\mu=0.5$ , $C_r = 0.6$	$\alpha=0.1$ , $\beta=0.1$
4D	$a\in[0.251, 0.45]$ , $b\in[-0.1, 0.3]$	$q_i\in[0.3, 0.999]$	$\mu=0.5$ , $C_r = 0.6$	$\alpha=0.3$ , $\beta=0.1$

| Show Table

DownLoad: CSV

Table 5. Optimization results of the 3D and 4D FO chaotic systems by DE.

FO chaotic system	Generations	Parameters	Fractional order derivative	Number of points	LEs	$D_{KY}$	Dynamical behavior
3D	50	$a=0.351$	$q_1=0.976$ , $q_2=0.956$ , $q_3=0.906$	101	$\lambda_1=0.0023$ , $\lambda_2=-0.0004$ , $\lambda_3=-0.0034$	2.5282	chaotic
	100	$a=0.322$	$q_1=0.986$ , $q_2=0.927$ , $q_3=0.971$	92	$\lambda_1=0.0029$ , $\lambda_2=-0.0005$ , $\lambda_3=-0.0032$	2.6120	chaotic
4D	50	$a=3.291$ , $b = 0.099$	$q_1=0.936$ , $q_2=0.965$ , $q_3=0.936$ , $q_4=0.887$	42	$\lambda_1=0.0020$ , $\lambda_2=0$ , $\lambda_3=-0.0003$ , $\lambda_4=-0.0023$	3.6803	chaotic
	100	$a=2.688$ , $b = 0.152$	$q_1=0.951$ , $q_2=0.970$ , $q_3=0.910$ , $q_4=0.970$	41	$\lambda_1=0.0144$ , $\lambda_2=0.0027$ , $\lambda_3=0$ , $\lambda_4=-0.0201$	3.8543	hyperchaotic

| Show Table

DownLoad: CSV

Table 6. Optimization results of the 3D and 4D FO chaotic systems by APSO.

APSO optimization results
FO chaotic system	Generations	Parameters	Fractional order derivative	Number of points	LEs	$D_{KY}$	Dynamical behavior
3D	50	$a=0.378$	$q_1=0.999$ , $q_2=0.969$ , $q_3=0.902$	101	$\lambda_1=0.0026$ , $\lambda_2=-0.0005$ , $\lambda_3=-0.0032$	2.66	chaotic
	109	$a=0.38$	$q_1=0.998$ , $q_2=0.970$ , $q_3=0.902$	113	$\lambda_1=0.0034$ , $\lambda_2=-0.0006$ , $\lambda_3=-0.0045$	2.5978	chaotic
4D	50	$a=3.193$ , $b = 0.299$	$q_1=0.965$ , $q_2=0.980$ , $q_3=0.985$ , $q_4=0.967$	40	$\lambda_1=0.0007$ , $\lambda_2=0.0003$ , $\lambda_3=-0.0001$ , $\lambda_4=-0.0011$	3.7633	hyperchaotic
	100	$a=3.106$ , $b = 0.19$	$q_1=0.878$ , $q_2=0.957$ , $q_3=0.945$ , $q_4=0.970$	41	$\lambda_1=0.0008$ , $\lambda_2=0.0002$ , $\lambda_3=0$ , $\lambda_4=-0.0015$	3.6358	hyperchaotic

| Show Table

DownLoad: CSV

7. Conclusions

In this paper, the optimization of the $D_{KY}$ of a 3D and a 4D FO chaotic systems was performed by applying two metaheuristics, DE and APSO. The main contribution of this work is the use of Poincaré maps for validating the chaotic behavior in FO chaotic systems with hidden attractors. Herein a simple numerical procedure was proposed for the analysis of points distribution in Poincaré maps. It has been shown that with this simple procedure, it is possible to differentiate a chaotic attractor from a periodic or quasi-periodic attractor. Moreover, the use of the Numba JIT compiler in the Python programming language has been an essential tool to save execution time, while preserving the versatility in the implementation of the metaheuristics. This last part is essential and it has been demonstrated that it facilitates the optimization process of FO chaotic systems with hidden attractors. Finally, it is worthy to mention that Poincaré maps can be used in both FO and integer chaotic systems regardless of the attractor nature (hidden or self-excited). The proposed approach can be applied in FO chaotic systems with self-excited attractor as done in integer chaotic systems with hidden and self-excited attractors.

Acknowledgments

Daniel Clemente-López acknowledges CONACYT for PhD scholarship with CVU no. 859036.

Conflict of interest

The authors declare no conflicts of interests.

References

[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]	S. Zhai, Z. Weng, Y. Mo, X. Feng, Energy dissipation and maximum bound principle preserving scheme for solving a nonlocal ternary Allen–Cahn model, Comput. Math. Appl., 155 (2024), 150–164. https://doi.org/10.1016/j.camwa.2023.06.018 doi: 10.1016/j.camwa.2023.06.018
[3]	M. Emamjomeh, M. Nabati, A. Dinmohammadi, Numerical study of two operator splitting localized radial basis function method for Allen–Cahn problem, Eng. Anal. Bound. Elem., 163 (2024), 126–137. https://doi.org/10.1016/j.enganabound.2023.07.014 doi: 10.1016/j.enganabound.2023.07.014
[4]	L. Q. Chen, Phase-field models for microstructure evolution, Ann. 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
[5]	M. Fatima, R. P. Agarwal, M. Abbas, P. O. Mohammed, M. Shafiq, N. Chorfi, Extension of Cubic B-Spline for Solving the Time-Fractional Allen–Cahn Equation in the Context of Mathematical Physics, Computation, 12 (2024), 51. https://doi.org/10.3390/computation12030051 doi: 10.3390/computation12030051
[6]	Z. Lu, J. Wang, A novel and efficient multi-scale feature extraction method for EEG classification, AIMS Math., 9 (2024), 16605–16622. https://doi.org/10.3934/math.2024848 doi: 10.3934/math.2024848
[7]	S. Kim, J. Kim, Automatic Binary Data Classification Using a Modified Allen–Cahn Equation, Int. J. Pattern Recognit. Artif. Intell., 35 (2021), 2150013. https://doi.org/10.1142/S021800142150013X doi: 10.1142/S021800142150013X
[8]	D. Lee, S. Lee, Image segmentation based on modified fractional Allen–Cahn equation, Math. Probl. Eng., 2019 (2019), 3980181. https://doi.org/10.1155/2019/8059716 doi: 10.1155/2019/8059716
[9]	M. Beneš, V. Chalupeck'y, K. Mikula, Geometrical image segmentation by the Allen–Cahn equation, Appl. Numer. Math., 51 (2004), 187–205. https://doi.org/10.1016/j.apnum.2004.04.006 doi: 10.1016/j.apnum.2004.04.006
[10]	D. Lee, Gradient-descent-like scheme for the Allen–Cahn equation, AIP Adv., 13 (2023), 085010. https://doi.org/10.1063/5.0154657 doi: 10.1063/5.0154657
[11]	D. Lee, Computing the area-minimizing surface by the Allen–Cahn equation with the fixed boundary, AIMS Math., 8 (2023), 23352–23371. https://doi.org/10.3934/math.20231184 doi: 10.3934/math.20231184
[12]	J. Yang, Y. Li, C. Lee, Y. Choi, J. Kim, Fast evolution numerical method for the Allen–Cahn equation, J. King Saud Univ. Sci., 35 (2023), 102430. https://doi.org/10.1016/j.jksus.2022.102430 doi: 10.1016/j.jksus.2022.102430
[13]	B. Xia, R. Yu, X. Song, X. Zhang, J. Kim, An efficient data assimilation algorithm using the Allen–Cahn equation, Eng. Anal. Bound. Elem., 155 (2023), 511–517. https://doi.org/10.1016/j.enganabound.2023.07.005 doi: 10.1016/j.enganabound.2023.07.005
[14]	H. Zhang, X. Qian, S. Song, Third-order accurate, large time-stepping and maximum-principle-preserving schemes for the Allen–Cahn equation, Numer. Algorithms, 95 (2024), 1213–1250. https://doi.org/10.1007/s11075-023-01482-w doi: 10.1007/s11075-023-01482-w
[15]	J. Sun, H. Zhang, X. Qian, S. Song, Up to eighth-order maximum-principle-preserving methods for the Allen–Cahn equation, Numer. Algorithms, 92 (2023), 1041–1062. https://doi.org/10.1007/s11075-022-01324-w doi: 10.1007/s11075-022-01324-w
[16]	J. Choi, S. Ham, S. Kwak, Y. Hwang, J. Kim, Stability analysis of an explicit numerical scheme for the Allen–Cahn equation with high-order polynomial potentials, AIMS Math., 9 (2024), 19332–19344. https://doi.org/10.3934/math.2024803 doi: 10.3934/math.2024803
[17]	H. Kim, G. Lee, S. Kang, S. Ham, Y. Hwang, J. Kim, Hybrid numerical method for the Allen–Cahn equation on nonuniform grids, Comput. Math. Appl., 158 (2024), 167–178. https://doi.org/10.1016/j.camwa.2023.07.010 doi: 10.1016/j.camwa.2023.07.010
[18]	X. Chen, X. Qian, S. Song, Fourth-order structure-preserving method for the conservative Allen–Cahn equation, Adv. Appl. Math. Mech., 15 (2023), 159–181. https://doi.org/10.4208/aamm.OA-2021-0303 doi: 10.4208/aamm.OA-2021-0303
[19]	J. Rubinstein, P. Sternberg, Nonlocal reaction-diffusion equations and nucleation, IMA J. Appl. Math., 48 (1992), 249–264. https://doi.org/10.1093/imamat/48.3.249 doi: 10.1093/imamat/48.3.249
[20]	M. Brassel, E. Bretin, A modified phase field approximation for mean curvature flow with conservation of the volume, Math. Meth. Appl. Sci., 10 (2011), 1157–1180. https://doi.org/10.1002/mma.1348 doi: 10.1002/mma.1348
[21]	J. Yang, J. Kim, Efficient and structure-preserving time-dependent auxiliary variable method for a conservative Allen–Cahn type surfactant system, Eng. Comput., 38 (2022), 5231–5250. https://doi.org/10.1007/s00366-021-01583-5 doi: 10.1007/s00366-021-01583-5
[22]	L. Bronsard, B. Stoth, Volume-preserving mean curvature flow as a limit of a nonlocal Ginzburg–Landau equation, SIAM J. Math. Anal., 28 (1997), 769–807. https://doi.org/10.1137/S0036141096298974 doi: 10.1137/S0036141096298974
[23]	X. Yang, J. J. Feng, C. Liu, J. Shen, Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method, J. Comput. Phys., 218 (2006), 417–428. https://doi.org/10.1016/j.jcp.2006.02.010 doi: 10.1016/j.jcp.2006.02.010
[24]	Z. Zhang, H. Tang, An adaptive phase field method for the mixture of two incompressible fluids, Comput. Fluids, 36 (2007), 1307–1318. https://doi.org/10.1016/j.compfluid.2006.10.001 doi: 10.1016/j.compfluid.2006.10.001
[25]	B. Xia, Y. Li, Z. Li, Second-order unconditionally stable direct methods for Allen–Cahn and conservative Allen–Cahn equations on surfaces, Mathematics, 8 (2020), 1486. https://doi.org/10.3390/math8091486 doi: 10.3390/math8091486
[26]	Z. Sun, S. Zhang, A radial basis function approximation method for conservative Allen–Cahn equations on surfaces, Appl. Math. Lett., 143 (2023), 108634. https://doi.org/10.1016/j.aml.2023.108634 doi: 10.1016/j.aml.2023.108634
[27]	X. Liu, Q. Hong, H. L. Liao, Y. Gong, A multi-physical structure-preserving method and its analysis for the conservative Allen–Cahn equation with nonlocal constraint, Numer. Algorithms, 97 (2024), 1–21. https://doi.org/10.1007/s11075-023-01502-9 doi: 10.1007/s11075-023-01502-9
[28]	J. Yang, J. Kim, Numerical study of incompressible binary fluids on 3D curved surfaces based on the conservative Allen–Cahn–Navier–Stokes model, Comput. Fluids, 228 (2021), 105094. https://doi.org/10.1016/j.compfluid.2021.105094 doi: 10.1016/j.compfluid.2021.105094
[29]	J. Yang, J. Kim, A phase-field model and its efficient numerical method for two-phase flows on arbitrarily curved surfaces in 3D space, Comput. Meth. Appl. Mech. Eng., 372 (2020), 113382. https://doi.org/10.1016/j.cma.2020.113382 doi: 10.1016/j.cma.2020.113382
[30]	C. Lee, S. Kim, S. Kwak, Y. Hwang, S. Ham, S. Kang, J. Kim, Semi-automatic fingerprint image restoration algorithm using a partial differential equation, AIMS Math., 8 (2023), 27528–27541. https://doi:10.3934/math.20231408 doi: 10.3934/math.20231408
[31]	Y. Choi, J. Kim, Maximum principle preserving and unconditionally stable scheme for a conservative Allen–Cahn equation, Eng. Anal. Bound. Elem., 150 (2023), 111–119. https://doi.org/10.1016/j.enganabound.2023.05.005 doi: 10.1016/j.enganabound.2023.05.005
[32]	J. Kim, S. Lee, Y. Choi, A conservative Allen–Cahn equation with a space–time dependent Lagrange multiplier, Int. J. Eng. Sci., 84 (2014), 11–17. https://doi.org/10.1016/j.ijengsci.2014.06.004 doi: 10.1016/j.ijengsci.2014.06.004
[33]	Y. Hwang, J. Yang, G. Lee, S. Ham, S. Kang, S. Kwak, et al., Fast and efficient numerical method for solving the Allen–Cahn equation on the cubic surface, Math. Comput. Simul., 215 (2024), 338–356. https://doi.org/10.1016/j.matcom.2023.07.024 doi: 10.1016/j.matcom.2023.07.024
[34]	Y. Wang, X. Xiao, X. Feng, Numerical simulation for the conserved Allen–Cahn phase field model of two-phase incompressible flows by an efficient dimension splitting method, Commun. Nonlinear Sci. Numer. Simul., 131 (2024), 107874. https://doi.org/10.1016/j.cnsns.2024.107874 doi: 10.1016/j.cnsns.2024.107874
[35]	X. Yang, Efficient, second-order in time, and energy stable scheme for a new hydrodynamically coupled three components volume-conserved Allen–Cahn phase-field model, Math. Models Meth. Appl. Sci., 31 (2021), 753–787. https://doi.org/10.1142/S0218202521500184 doi: 10.1142/S0218202521500184
[36]	W. Cai, J. Ren, X. Gu, Y. Wang, Parallel and energy conservative/dissipative schemes for sine–Gordon and Allen–Cahn equations, Comput. Meth. Appl. Mech. Eng., 425 (2024), 116938. https://doi.org/10.1016/j.cma.2024.116938 doi: 10.1016/j.cma.2024.116938
[37]	J. Kim, D. Jeong, S. D. Yang, Y. Choi, A finite difference method for a conservative Allen–Cahn equation on non-flat surfaces, J. Comput. Phys., 334 (2017), 170–181. https://doi.org/10.1016/j.jcp.2016.12.060 doi: 10.1016/j.jcp.2016.12.060

This article has been cited by:

Ghada Mohammed Mansour, Haroon Mohamed Barakat, Islam Abdullah Husseiny, Magdy Nagy, Ahmed Hamdi Mansi, Metwally Alsayed Alawady, Measures of cumulative residual Tsallis entropy for concomitants of generalized order statistics based on the Morgenstern family with application to medical data, 2025, 22, 1551-0018, 1572, 10.3934/mbe.2025058

Reader Comments

Your name:*

Email:*
© 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)