Next Article in Journal
Characterization of Aeroelastic Behavior in a High Aspect Ratio Wing Using Computational and Wind Tunnel Experiments
Next Article in Special Issue
Fractional p-Laplacian Coupled Systems with Multi-Point Boundary Conditions
Previous Article in Journal
“Computational Mathematics and Mathematical Physics”—Editorial I (2021–2023)
Previous Article in Special Issue
Analysis of a Fractional-Order Quadratic Functional Integro-Differential Equation with Nonlocal Fractional-Order Integro-Differential Condition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Auxiliary Function Method for Analyzing Nonlinear System of Belousov–Zhabotinsky Equation with Caputo Operator

1
Department of Mathematical Sciences, Faculty of Sciences, Princess Nourah Bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
2
Department of Basic Sciences, Preparatory Year Deanship, King Faisal University, Al-Ahsa 31982, Saudi Arabia
3
Department of Mathematics, Abdul Wali Khan University Mardan, Mardan 23200, Pakistan
*
Author to whom correspondence should be addressed.
Submission received: 29 July 2023 / Revised: 25 August 2023 / Accepted: 25 August 2023 / Published: 28 August 2023
(This article belongs to the Special Issue Recent Advances in Fractional Differential Equations and Inequalities)

Abstract

:
This paper introduces the optimal auxiliary function method (OAFM) for solving a nonlinear system of Belousov–Zhabotinsky equations. The system is characterized by its complex dynamics and is treated using the Caputo operator and concepts from fractional calculus. The OAFM provides a systematic approach to obtain approximate analytical solutions by constructing an auxiliary function. By optimizing the parameters of the auxiliary function, an approximate solution is derived that closely matches the behavior of the original system. The effectiveness and accuracy of the OAFM are demonstrated through numerical simulations and comparisons with existing methods. Fractional calculus enhances the understanding and modeling of the nonlinear dynamics in the Belousov–Zhabotinsky system. This study contributes to fractional calculus and nonlinear dynamics, offering a powerful tool for analyzing and solving complex systems such as the Belousov–Zhabotinsky equation.

1. Introduction

Fractional nonlinear systems of partial differential equations (PDEs) have garnered substantial research interest owing to their capacity to describe diverse phenomena across areas such as engineering, physics, and biology. Fractional derivatives, non-local operator that account for historical impacts, and complicated non-linearities are all part of these schemes. The usefulness of fractional calculus, which deals with fractional derivatives and integrals, has been demonstrated in real-world circumstances by simulating memory effects, long-range interactions, and anomalous diffusion. Simultaneously, nonlinear dynamics investigates the complicated behavior of complex systems that exhibit nonlinear interactions among their basic components [1,2,3]. The combination of these two domains has resulted in fractional nonlinear partial differential equation (PDE) systems capable of capturing the complicated dynamics found in several physical and biological systems. The analysis of fractional nonlinear PDE systems has provided a multitude of intriguing discoveries and practical applications within this rapidly developing field of study. Current research efforts are mostly focused on diving into the stability and bifurcation characteristics of solutions, as well as their uniqueness and existence, which is frequently supported by numerical simulations. These systems have been applied to a variety of disciplines, including fluid dynamics, electrochemistry, and population dynamics, each of which provides distinct insights and opportunities for investigation [4,5,6,7,8,9].
One class of nonlinear systems that has attracted significant attention is the nonlinear system of Belousov–Zhabotinsky equation (NSBZE). The NSBZE is a mathematical model that describes the complex behavior of reaction–diffusion systems exhibiting chemical oscillations and pattern formation. It finds applications in diverse areas such as chemistry, biology, and material science [10,11]. The nonlinear system of Belousov–Zhabotinsky equation, initially rooted in chemical reactions, has broadened its impact across multiple scientific domains. It is a cornerstone for understanding chemical kinetics, reaction–diffusion dynamics, and pattern formation. This equation finds applications in diverse fields, including materials science, biological systems, neural networks, artificial life, optics, and unconventional computing. Its insights inform the study of self-organizing materials, neural pattern formation, brain dynamics, and innovative information processing methods. Moreover, the equation’s visually captivating patterns and oscillations make it a valuable tool for educational outreach and fostering a deeper understanding of complex behaviors and nonlinear dynamics in various scientific contexts. To tackle the NSBZE, fractional calculus provides a powerful mathematical framework. Fractional calculus extends the concept of differentiation and integration to non-integer orders, allowing for analyzing systems with memory and long-range dependencies. The Caputo operator is commonly employed in fractional calculus, a fractional derivative operator widely used to describe the dynamics of physical systems governed by fractional differential equations [12,13,14,15,16,17].
The optimal auxiliary function method (OAFM) has proven to be a powerful tool for solving various types of nonlinear differential equations. This method introduces an auxiliary function that satisfies an ordinary differential equation optimized to yield approximate solutions for the original problem. The OAFM has been used effectively in various areas, including fluid dynamics, chemical reactions, and population dynamics [18,19,20,21]. In this paper, we aim to combine the optimal auxiliary function method with the Caputo operator in the context of the nonlinear system of Belousov–Zhabotinsky equation. We aim to develop an effective numerical technique for obtaining approximate solutions to this challenging nonlinear system. By utilizing the capabilities of the OAFM and the mathematical foundation provided by fractional calculus, we seek to enhance the accuracy and efficiency of the solution process.

2. Preliminaries

Definition 1.
The Caputo fractional derivative of a function R ( ζ , τ ) with a fractional order α is expressed as [18]
C D t α R ( ζ , t ) = J t m α R m ( ζ , t ) , m 1 < α m , t > 0 .
Definition 2.
The expression for the Riemann fractional integral can be presented as follows [18]:
J t α R ( ζ , t ) = 1 Γ ( α ) 0 t ( t r ) α 1 R ( ζ , r ) d r
Lemma 1.
For n 1 < α n , p > 1 , t 0 and λ R , we have [18]:
1. 
D t α t p = Γ ( α + 1 ) Γ ( p α + 1 ) t p α
2. 
D t α λ = 0
3. 
D t α I t α R ( ζ , t ) = R ( ζ , t )
4. 
I t α D t α R ( ζ , t ) = R ( α , t ) i = 0 n 1 i R ( ζ , 0 ) t i i !

3. General Procedure of OAFM

In this section, we consider the most general type of a nonlinear fractional order differential equation [18]:
α R ( ξ , τ ) τ α = g ( ξ , τ ) + N ( R ( ξ , τ ) )
with the initial condition
k R ( ξ , 0 ) τ k = h k ( ξ ) , k = 0 ,   1 ,   2 ,   .
where α τ α represents the fractional derivative of order α , g ( ξ , τ ) is the unknown function, and N ( R ( ξ , τ ) ) is the nonlinear function involving R, spatial variable ξ , and time variable τ .
  • Step 1: To solve Equation (3), we will use an approximate solution that has two components, such as:
    R ( ξ , τ ) = R 0 ( ξ , τ ) + R 1 ( ξ , τ , C i ) , i = 1 ,   2 ,   3 ,   ,   p
  • Step 2: To determine the zero and first-order solutions, we substitute Equation (5) into Equation (3), which results in:
    α R 0 ( ξ , τ ) τ α + α R 1 ( ξ , τ ) τ α + g ( ξ , τ ) + N α R 0 ( ξ , τ ) τ α + α R 1 ( ξ , τ , C i ) τ α = 0
  • Step 3: For the purpose of determining the first approximation R 0 ( ξ , τ ) based on the linear equation.
    α R 0 ( ξ , τ ) τ α + g ( ξ , τ ) = 0
Using the inverse operator, we arrive at R 0 ( ξ , τ ) , which is expressed as follows:
R 0 ( ξ , τ ) = g ( ξ , τ )
  • Step 4: The nonlinear term seen in expanding form (6) is,
    N ρ R 0 ( ξ , τ ) τ α + R 1 α ( ξ , τ , C i ) t α = N R 0 ( ξ , τ ) + k = 1 R 1 k k ! N ( k ) R 0 ( ξ , τ )
  • Step 5: To quickly solve Equation (9) and accelerate the convergence of the first-order approximations R ( ξ , τ ) , we propose the following alternative equation:
    α R 1 ( ξ , τ , C i ) τ α = A 1 R 0 ( ξ , τ ) N R 0 ( ξ , τ ) + A 2 R 0 ( ξ , τ ) , C j .
  • Step 6: We find a first-order solutions R 1 ( ξ , τ ) by using the inverse operator after putting the Auxiliary function into Equation (10).
  • Step 7: Several methods are used to find the numerical values of the convergence control parameters C i , including least squares, Galerkin’s, Ritz, and collocation. To eliminate mistakes, we employ the least squares approach.
    J ( C i , C j ) = 0 t Ω 2 ( y , t ; C i , C j ) d ξ d τ .
    where denotes the residual,
    ( ξ , τ , C i , C j ) = R ( ξ , τ , C i , C j ) τ + g ( ξ , τ ) + N [ R ( ξ , τ , C i , C j ) ] , i = 1 ,   2 ,   ,   s , j = S + 1 ,   S + 2 ,   ,   p

4. Applications

4.1. Problem 1

Let us consider the Belousov–Zhabotinsky System of fractional order, which is represented by
D τ α R ( ξ , τ ) R ( ξ , τ ) 2 R ( ξ , τ ) ξ 2 + R 2 ( ξ , τ ) + λ T ( ξ , τ ) = 0 , D τ α T ( ξ , τ ) 2 T ( ξ , τ ) ξ 2 + γ T ( ξ , τ ) T ( ξ , τ ) = 0 , where , 0 < α 1
Subject to the following IC’s:
R ( ξ , 0 ) = 1 2 1 t a n h 2 ( ξ 2 ) , T ( ξ , 0 ) = 1 2 + t a n h ( ξ 2 ) + 1 2 t a n h 2 ( ξ 2 )
Consider linear terms in Equation (13)
L ( R 0 ( ξ , τ ) ) = α R 0 ( ξ , τ ) τ α , L ( T 0 ( ξ , τ ) ) = α T 0 ( ξ , τ ) τ α
Consider nonlinear terms in Equation (13)
N ( R 0 ( ξ , τ ) ) = R 0 ( ξ , τ ) 2 R 0 ( ξ , τ ) ξ 2 + R 0 2 ( ξ , τ ) + λ T 0 ( ξ , τ ) , N ( T 0 ( ξ , τ ) ) = 2 T 0 ( ξ , τ ) ξ 2 + γ T 0 ( ξ , τ ) R 0 ( ξ , τ )
Zeroth-order approximation
α R 0 ( ξ , τ ) τ α = 0 α T 0 ( ξ , τ ) τ α = 0
By apply of the inverse operator, we obtain result as:
R 0 ( ξ , τ ) = 1 2 1 t a n h 2 ( ξ 2 ) , T 0 ( ξ , τ ) = 1 2 + t a n h ( ξ 2 ) + 1 2 t a n h 2 ( ξ 2 )
By making use of Equation (18) in Equation (16), the system of nonlinear terms becomes
N [ R 0 ( ξ , τ ) ] = 1 14 sech 4 ξ 2 ( λ ( 2 sinh ( ξ ) + sinh ( 2 x ) 2 ) 2 ( λ 2 ) cosh ( ξ ) ) N [ T 0 ( ξ , τ ) ] = 1 4 sech 4 ξ 2 ( γ ( sinh ( ξ ) ) + γ + sinh ( ξ ) + cosh ( ξ ) 2 )
We choose the auxiliary function as:
A 1 = c 1 1 2 1 t a n h 2 ( ξ 2 ) , A 2 = c 2 1 2 1 t a n h 2 ( ξ 2 ) 3 , A 3 = c 3 1 2 + t a n h ( ξ 2 ) + 1 2 t a n h 2 ( ξ 2 ) 5 A 4 = c 4 1 2 + t a n h ( ξ 2 ) + 1 2 t a n h 2 ( ξ 2 ) 7
First-order approximation using the OAFM approach, as described in Section 3.
α R 1 ( ξ , τ ) τ α = A 1 [ R 0 ( ξ , τ ) N [ R 0 ( ξ , τ ) ] + A 2 [ R 0 ( ξ , τ ) ; c j ] α T 1 ( ξ , τ ) τ α = A 3 [ T 0 ( ξ , τ ) N [ T 0 ( ξ , τ ) ] + A 4 [ T 0 ( ξ , τ ) ; c j ]
using Equations (19) and (20) in Equation (21).
α R 1 ( ξ , τ ) τ ρ = 1 8 sech 6 ξ 2 ( c 1 λ sinh ( ξ ) + c 1 cosh ( ξ ) ( λ sinh ( ξ ) λ + 2 ) c 1 λ + c 2 ) α T 1 ( ξ , τ ) τ ρ = 1 256 ( sinh ( ξ ) 1 ) 5 sech 14 ξ 2 ( 2 sinh ( ξ ) ( c 3 γ + c 3 2 c 4 ) + 2 c 3 ( γ 2 ) +   2 c 3 cosh ( ξ ) c 4 cosh ( 2 ξ ) + c 4 )
Using the inverse operator to Equation (22) we get
R 1 ( ξ , τ ) = τ α sech 6 ξ 2 ( c 1 λ sinh ( ξ ) + c 1 cosh ( ξ ) ( λ sinh ( ξ ) λ + 2 ) c 1 λ + c 2 ) 8 Γ ( 1 + α ) T 1 ( ξ , τ ) = τ α sech 14 ξ 2 ( sinh ( ξ ) 1 ) 5 ( 2 c 3 ( γ 2 ) + 2 c 3 cosh ( x ) + c 4 cosh ( 2 ξ ) + c 4 ) 256 Γ ( α + 1 )   τ α ( sinh ( ξ ) 1 ) 5 sech 14 ξ 2 ( 2 sinh ( ξ ) ( c 3 γ + c 3 2 c 4 ) ) 256 Γ ( α + 1 )
According to the OAFM procedure
R ( ξ , τ ) = R 0 ( ξ , τ ) + R 1 ( ξ , τ ) T ( ξ , τ ) = T 0 ( ξ , τ ) + T 1 ( ξ , τ )
In this discussion, we analyze the behavior of the Belousov–Zhabotinsky (BZ) system, as shown in Figure 1 and Figure 2. The BZ system is an example of a nonlinear chemical oscillator and is widely studied due to its oscillatory behavior, which can be observed in various fields such as biology, chemistry, and physics. We focus on the response functions R ( ξ , τ ) and T ( ξ , τ ) , which describe the behavior of the system under varying conditions. The parameter α plays a crucial role in these response functions and is investigated in this discussion. Figure 1a presents a 2D plot of the response function R ( ξ , τ ) at different values of the parameter α . As the parameter α increases, the plot illustrates the change in the oscillatory behavior of the BZ system. The response function exhibits more complex oscillatory patterns with increasing values of α . This suggests that the parameter α has a significant impact on the system’s response and can be used to control its behavior. Similarly, Figure 1b shows a 2D plot of the response function T ( ξ , τ ) at different values of the parameter α . The plot reveals that the oscillatory behavior of the BZ system also changes with varying values of α in the case of T ( ξ , τ ) . This further supports the idea that the parameter α is a key factor in determining the overall behavior of the BZ system. To provide a more comprehensive understanding of the BZ system, Figure 2a,b display 3D plots of the response functions R ( ξ , τ ) and T ( ξ , τ ) , respectively, at different values of the parameter α . These 3D plots offer a more detailed perspective on the complex oscillatory behavior of the BZ system. Figure 2a demonstrates how the 3D plot of the response function R ( ξ , τ ) changes as the parameter α is varied. The plot exhibits intricate patterns, indicating the complex relationship between ξ , τ , and α . Similarly, Figure 2b showcases a 3D plot of the response function T ( ξ , τ ) at different values of α . This plot also displays intricate patterns, further highlighting the significance of the parameter α in the BZ system’s behavior. In Table 1 and Table 2, a comparison of the R ( ξ , τ ) and T ( ξ , τ ) OAFM solution and exact solution and their corresponding absolute error at fractional order α = 1 where λ = 2 and γ = 3 for Problem 1. In conclusion, the graphical analysis of Figure 1 and Figure 2 clearly demonstrates that the parameter α plays a crucial role in determining the behavior of the BZ system. The response functions R ( ξ , τ ) and T ( ξ , τ ) exhibit complex oscillatory patterns that change with varying values of α . This analysis is beneficial for understanding and controlling the oscillatory behavior of the BZ system and other similar nonlinear chemical oscillators.
The values of auxiliary constant obtained by collocation method.
c 1 = 2.358192481417809 , c 2 = 0.7753100423748968 , c 3 = 83.49971233901873 and c 4 = 33.24684800337086

4.2. Problem 2

Let us consider the system of a nonlinear fractional-order equation, which is represented by
D τ α R ( ξ , τ ) T ( ξ , τ ) 2 R ( ξ , τ ) ξ 2 R ( ξ , τ ) 2 T ( ξ , τ ) ξ 2 2 + 2 ξ 2 + 2 τ 2 = 0 , D τ α T ( ξ , τ ) T ( ξ , τ ) 2 T ( ξ , τ ) ξ 2 + R ( ξ , τ ) 2 R ( ξ , τ ) ξ 2 1 3 2 ξ 2 3 2 τ 2 = 0 , where , 1 < α 2
Subject to the following IC’s:
R ( ξ , 0 ) = ξ 2 , τ R ( ξ , 0 ) = 0 T ( ξ , 0 ) = 1 2 ξ 2 , τ T ( ξ , 0 ) = 0
Consider linear terms in Equation (25)
L ( R 0 ( ξ , τ ) ) = α R 0 ( ξ , τ ) τ α , L ( T 0 ( ξ , τ ) ) = α T 0 ( ξ , τ ) τ α
Consider nonlinear terms in Equation (25)
N ( R 0 ( ξ , τ ) ) = T 0 ( ξ , τ ) 2 R 0 ( ξ , τ ) ξ 2 R 0 ( ξ , τ ) 2 T 0 ( ξ , τ ) ξ 2 , N ( T 0 ( ξ , τ ) ) = T 0 ( ξ , τ ) 2 T 0 ( ξ , τ ) ξ 2 + R 0 ( ξ , τ ) 2 R 0 ( ξ , τ ) ξ 2
Zeroth-order approximation
α R 0 ( ξ , τ ) τ α = 0 α T 0 ( ξ , τ ) τ α = 0
By applying of the inverse operator, we obtain result as:
R 0 ( ξ , τ ) = ξ 2 , T 0 ( ξ , τ ) = 1 2 ξ 2
By making use of Equation (30) in Equation (28), the system of nonlinear terms becomes
N [ R 0 ( ξ , τ ) ] = 2 ξ 2 N [ T 0 ( ξ , τ ) ] = 3 2 ξ 2
We choose the auxiliary function as:
A 1 = ( c 1 + ξ c 2 ) ξ 2 , A 2 = c 3 ( ξ 2 ) 3 , A 3 = ( c 4 + ξ c 5 ) ( 3 2 ξ 2 ) 4 , A 4 = c 6 ( 3 2 ξ 2 ) 7
First-order approximation using the OAFM approach, as described in Section 3.
α R 1 ( ξ , τ ) τ α = A 1 [ R 0 ( ξ , τ ) N [ R 0 ( ξ , τ ) ] + A 2 [ R 0 ( ξ , τ ) ; c j ] α T 1 ( ξ , τ ) τ α = A 3 [ T 0 ( ξ , τ ) N [ T 0 ( ξ , τ ) ] + A 4 [ T 0 ( ξ , τ ) ; c j ]
using Equations (30) and (32) in Equation (33)
α R 1 ( ξ , τ ) τ ρ = 2 ξ 4 ( c 1 + c 2 ξ ) c 3 ξ 6 α T 1 ( ξ , τ ) τ ρ = 3 32 ξ 10 ( c 4 + c 5 ξ ) c 6 ξ 14 128
applying inverse operator to Equation (34) we get
R 1 ( ξ , τ ) = ξ 4 τ α ( 2 c 1 + ξ ( 2 c 2 c 3 ξ ) ) Γ ( α + 1 ) T 1 ( ξ , τ ) = ξ 10 τ α 12 c 4 + 12 c 5 ξ + c 6 ξ 4 128 Γ ( α + 1 )
According to the OAFM procedure
R ( ξ , τ ) = R 0 ( ξ , τ ) + R 1 ( ξ , τ ) T ( ξ , τ ) = T 0 ( ξ , τ ) + T 1 ( ξ , τ )
The values of auxiliary constant obtained by collocation method.
c 1 = 31,430.6586668333 , c 2 = 110,362.31997412356 , c 3 = 182,153.25836719057 , c 4 = 125,830.5231198654 , c 5 = 41,444.51349867818 and c 6 = 2.0600130170669805 × 10 6
Figure 3 and Figure 4 present the graphical representations of the reflection, R ( ξ , τ ) , and transmission, T ( ξ , τ ) , coefficients obtained for a nonlinear fractional system of partial differential equations using the optimal amplitude frequency modulation (OAFM) method. Figure 3 illustrates the 2D plots of the coefficients: (a) represents the 2D plot of the reflection coefficient, R ( ξ , τ ) , and (b) displays the 2D plot of the transmission coefficient, T ( ξ , τ ) . The horizontal axis in both plots represents the spatial variable, ξ , while the vertical axis represents the temporal variable, τ . The variation in color intensity indicates the magnitude of the coefficients. The 2D plots enable a clear visualization of the coefficients’ behavior with respect to the spatial and temporal variables, as well as how they change in the presence of nonlinearities. Figure 4, on the other hand, provides a more comprehensive representation of the coefficients through 3D plots. In this figure, (a) is the 3D plot of the reflection coefficient, R ( ξ , τ ) , while (b) is the 3D plot of the transmission coefficient, T ( ξ , τ ) . The horizontal plane represents the spatial variable, ξ , and the temporal variable, τ . The vertical axis illustrates the magnitude of the coefficients. By using a 3D representation, it becomes easier to visualize the distribution of the coefficients and how they evolve in the nonlinear fractional system. Table 3 and Table 4 are a comparison of the OAFM solution and exact solution and their corresponding absolute error at fractional order α = 2 for R ( ξ , τ ) and T ( ξ , τ ) for Problem 2. In conclusion, Figure 3 and Figure 4 offer a graphical representation of the reflection and transmission coefficients, R ( ξ , τ ) and T ( ξ , τ ) , for a nonlinear fractional system of partial differential equations using the OAFM method. The 2D and 3D plots allow for an in-depth understanding of the coefficients’ behavior and how they are affected by the nonlinearities present in the system.

5. Conclusions

In conclusion, the optimal Auxiliary function method (OAFM) has proven to be a powerful and effective technique for analyzing nonlinear systems, as demonstrated in this paper by applying it to the nonlinear system of the Belousov–Zhabotinsky equation. By utilizing the Caputo operator and concepts from fractional calculus, the OAFM provides a valuable approach to solving and understanding complex dynamic systems. The integration of Caputo operators further enhances the accuracy and authenticity of the model, enabling a more faithful representation of the underlying physical phenomena. The OAFM not only allows for accurate and efficient computations but also offers insights into the behavior and stability of the system. The results obtained using the OAFM highlight its effectiveness in addressing challenging problems in nonlinear dynamics. The combination of the OAFM, nonlinear systems, and fractional calculus opens up new avenues for future research, paving the way for further advancements and applications in various scientific and engineering domains.

Author Contributions

Conceptualization, A.S.A.; methodology, A.S.A.; software, H.Y.; validation, H.Y.; formal analysis, M.W.A.; investigation, M.W.A.; resources, A.K. and R.S.; data curation, A.K. and R.S.; Funding H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R183), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. This work was supported by the Deanship of Scientific Research, the Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant No. 4000).

Data Availability Statement

The numerical data used to support the findings of this study are included within the article. Mathematica codes for drawing the figures are available, which can be requested from R. Shah.

Acknowledgments

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R183), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. This work was supported by the Deanship of Scientific Research, the Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant No. 4000).

Conflicts of Interest

The authors declare that there are no conflict of interest regarding the publication of this article.

References

  1. Alderremy, A.A.; Aly, S.; Fayyaz, R.; Khan, A.; Shah, R.; Wyal, N. The Analysis of Fractional-Order Nonlinear Systems of Third Order KdV and Burgers Equations via a Novel Transform. Complexity 2022, 2022, 4935809. [Google Scholar] [CrossRef]
  2. Iqbal, N. The fractional-order system of singular and non-singular thermo-elasticity system in the sense of homotopy perturbation transform method. Fractals 2023. [Google Scholar] [CrossRef]
  3. Mukhtar, S.; Noor, S. The numerical investigation of a fractional-order multi-dimensional Model of Navier-Stokes equation via novel techniques. Symmetry 2022, 14, 1102. [Google Scholar] [CrossRef]
  4. Alderremy, A.A.; Shah, R.; Iqbal, N.; Aly, S.; Nonlaopon, K. Fractional Series Solution Construction for Nonlinear Fractional Reaction-Diffusion Brusselator Model Utilizing Laplace Residual Power Series. Symmetry 2022, 14, 1944. [Google Scholar] [CrossRef]
  5. Alqhtani, M.; Saad, K.M.; Shah, R.; Weera, W.; Hamanah, W.M. Analysis of the Fractional-Order Local Poisson Equation in Fractal Porous Media. Symmetry 2022, 14, 1323. [Google Scholar] [CrossRef]
  6. Iqbal, N.; Chughtai, M.T.; Ullah, R. Fractional Study of the Non-Linear Burgers’ Equations via a Semi-Analytical Technique. Fractal Fract. 2023, 7, 103. [Google Scholar] [CrossRef]
  7. Jin, H.-Y.; Wang, Z.-A. Global stabilization of the full attraction-repulsion Keller-Segel system. Discret. Contin. Dyn. Syst. A 2020, 40, 3509–3527. [Google Scholar] [CrossRef]
  8. Li, H.; Peng, R.; Wang, Z.-A. On a Diffusive Susceptible-Infected-Susceptible Epidemic Model with Mass Action Mechanism and Birth-Death Effect: Analysis, Simulations, and Comparison with Other Mechanisms. SIAM J. Appl. Math. 2018, 78, 2129–2153. [Google Scholar] [CrossRef]
  9. Jin, H.-Y.; Wang, Z.-A. Boundedness, blowup and critical mass phenomenon in competing chemotaxis. J. Differ. Equ. 2016, 260, 162–196. [Google Scholar] [CrossRef]
  10. Epstein, I.R.; Pojman, J.A. An Introduction to Nonlinear Chemical Dynamics: Oscillations, Waves, Patterns, and Chaos; Oxford University Press: Oxford, UK, 1998. [Google Scholar]
  11. Karaagac, B.; Esen, A.; Ucar, Y.; Yagmurlu, N.M. A new outlook for analysis of Noyes-Field model for the nonlinear Belousov-Zhabotinsky reaction using operator splitting method. Comput. Math. Appl. 2023, 136, 127–135. [Google Scholar] [CrossRef]
  12. Prakash, A.; Kaur, H.; Rahul; Goyal, M. A Reliable Numerical Approach for Liouville-Caputo Time-fractional Belousov-Zhabotinsky Equation Arising in Chemical Oscillatory Reaction. In Control and Measurement Applications for Smart Grid: Select Proceedings of SGESC 2021; Springer Nature: Singapore, 2022; pp. 383–395. [Google Scholar]
  13. Akinyemi, L. A fractional analysis of Noyes–Field model for the nonlinear Belousov—Zhabotinsky reaction. Comput. Appl. Math. 2020, 39, 175. [Google Scholar] [CrossRef]
  14. Li, Q.; Lin, H.; Tan, X.; Du, S. H Consensus for Multiagent-Based Supply Chain Systems Under Switching Topology and Uncertain Demands. IEEE Trans. Syst. Man, Cybern. Syst. 2020, 50, 4905–4918. [Google Scholar] [CrossRef]
  15. Wang, B.; Shen, Y.; Li, N.; Zhang, Y.; Gao, Z. An adaptive sliding mode fault-tolerant control of a quadrotor unmanned aerial vehicle with actuator faults and model uncertainties. Int. J. Robust Nonlinear Control. 2023. [Google Scholar] [CrossRef]
  16. Song, J.; Mingotti, A.; Zhang, J.; Peretto, L.; Wen, H. Accurate Damping Factor and Frequency Estimation for Damped Real-Valued Sinusoidal Signals. IEEE Trans. Instrum. Meas. 2022, 71, 6503504. [Google Scholar] [CrossRef]
  17. Ali, T.A.A.; Xiao, Z.; Jiang, H.; Li, B. A Class of Digital Integrators Based on Trigonometric Quadrature Rules. IEEE Trans. Ind. Electron. 2023, 1–11. [Google Scholar] [CrossRef]
  18. Liu, H.; Xu, Y.; Liu, J.; Chen, Y. Optimal Auxiliary Function Method for Fractional Differential Equations. J. Appl. Math. 2018, 2018, 1–13. [Google Scholar]
  19. He, J.-H. Optimal Auxiliary Function Method for Self-Adjoint Perturbed Nonlinear Problems. Nonlinear Sci. Lett. 2006, 1, 1–8. [Google Scholar]
  20. Chen, J.; Liu, Z.; Lomovtsev, F.E.; Obukhovskii, V. Optimal feedback control for a class of second-order evolution differential inclusions with Clarkes subdifferential. J. Nonlinear Var. Anal. 2022, 6, 551–565. [Google Scholar]
  21. Tang, Y. A characteristic mixed finite element method for bilinear convection-diffusion optimal control problems. J. Nonlinear Funct. Anal. 2022, 2022, 39. [Google Scholar]
Figure 1. The 2D-plots of R ( ξ , τ ) and T ( ξ , τ ) using the OAFM solution where λ = 1.5 and γ = 2.5 .
Figure 1. The 2D-plots of R ( ξ , τ ) and T ( ξ , τ ) using the OAFM solution where λ = 1.5 and γ = 2.5 .
Axioms 12 00825 g001
Figure 2. The 3D-plots of R ( ξ , τ ) and T ( ξ , τ ) using OAFM solution where the values λ = 1.5 and γ = 2.5 .
Figure 2. The 3D-plots of R ( ξ , τ ) and T ( ξ , τ ) using OAFM solution where the values λ = 1.5 and γ = 2.5 .
Axioms 12 00825 g002
Figure 3. The 2D-plots of R ( ξ , τ ) , T ( ξ , τ ) using OAFM solution at different values of fractional order.
Figure 3. The 2D-plots of R ( ξ , τ ) , T ( ξ , τ ) using OAFM solution at different values of fractional order.
Axioms 12 00825 g003
Figure 4. The 3D-plots of R ( ξ , τ ) , T ( ξ , τ ) using OAFM solution at different values of fractional order.
Figure 4. The 3D-plots of R ( ξ , τ ) , T ( ξ , τ ) using OAFM solution at different values of fractional order.
Axioms 12 00825 g004
Table 1. A comparison of the R ( ξ , τ ) OAFM solution and exact solution and their corresponding absolute error at fractional order α = 1 where λ = 2 and γ = 3 for Problem 1.
Table 1. A comparison of the R ( ξ , τ ) OAFM solution and exact solution and their corresponding absolute error at fractional order α = 1 where λ = 2 and γ = 3 for Problem 1.
ξ R ( ξ , τ ) OAFM R ( ξ , τ ) ExactAbs.error
0.1−0.4987−0.49875.1 × 10 5
0.2−0.495−0.49503.3 × 10 5
0.3−0.4889−0.48891.5 × 10 5
0.4−0.4805−0.48052.7 × 10 6
0.5−0.470−0.47002.0 × 10 5
0.6−0.4576−0.45753.7 × 10 5
0.7−0.4434−0.44345.2 × 10 5
0.8−0.4279−0.42786.6 × 10 5
0.9−0.411−0.4117.8 × 10 5
1−0.3933−0.39328.9 × 10 5
Table 2. Comparison of the T ( ξ , τ ) OAFM solution and exact solution and their corresponding absolute error at fractional order α = 1 for Problem 1.
Table 2. Comparison of the T ( ξ , τ ) OAFM solution and exact solution and their corresponding absolute error at fractional order α = 1 for Problem 1.
ξ T ( ξ , τ ) OAFM T ( ξ , τ ) ExactAbs.error
0.1−0.4487−0.44874.9 × 10 5
0.2−0.3953−0.39538.2 × 10 5
0.3−0.34−0.33991.0 × 10 4
0.4−0.2831−0.2831.1 × 10 4
0.5−0.2251−0.2251.2 × 10 4
0.6−0.1663−0.16611.2 × 10 4
0.7−0.1071−0.10691.2 × 10 4
0.8−0.0479−0.04781.2 × 10 4
0.90.01090.011021.2 × 10 4
10.068890.069011.1 × 10 4
Table 3. A comparison of the OAFM solution and exact solution and their corresponding absolute error at fractional order α = 2 for Problem 2.
Table 3. A comparison of the OAFM solution and exact solution and their corresponding absolute error at fractional order α = 2 for Problem 2.
ξ R ( ξ , τ ) OAFM R ( ξ , τ ) ExactAbs.error
0.10.010.014.06 × 10 7
0.150.02250.02252.72 × 10 6
0.20.040010.047.12 × 10 6
0.250.062510.06251.30 × 10 5
0.30.090020.091.90 × 10 5
0.350.122520.12252.10 × 10 5
0.40.160020.161.70 × 10 5
0.450.20250.20252.74 × 10 6
0.50.249980.252.20 × 10 5
Table 4. A comparison of the OAFM solution and exact solution and their corresponding absolute error at fractional order α = 2 for Problem 2.
Table 4. A comparison of the OAFM solution and exact solution and their corresponding absolute error at fractional order α = 2 for Problem 2.
ξ T ( ξ , τ ) OAFM T ( ξ , τ ) ExactAbs.error
0.100.0050.0052.000 × 10 8
0.150.011250.011252.0001 × 10 8
0.200.02000.02002.0025 × 10 8
0.250.031250.031252.0242 × 10 8
0.300.04500.045002.1515 × 10 8
0.350.061250.061252.7125 × 10 8
0.400.08000.08004.7134 × 10 8
0.450.101250.101251.0770 × 10 7
0.500.12500.125002.6870 × 10 7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alshehry, A.S.; Yasmin, H.; Ahmad, M.W.; Khan, A.; Shah, R. Optimal Auxiliary Function Method for Analyzing Nonlinear System of Belousov–Zhabotinsky Equation with Caputo Operator. Axioms 2023, 12, 825. https://doi.org/10.3390/axioms12090825

AMA Style

Alshehry AS, Yasmin H, Ahmad MW, Khan A, Shah R. Optimal Auxiliary Function Method for Analyzing Nonlinear System of Belousov–Zhabotinsky Equation with Caputo Operator. Axioms. 2023; 12(9):825. https://doi.org/10.3390/axioms12090825

Chicago/Turabian Style

Alshehry, Azzh Saad, Humaira Yasmin, Muhammad Wakeel Ahmad, Asfandyar Khan, and Rasool Shah. 2023. "Optimal Auxiliary Function Method for Analyzing Nonlinear System of Belousov–Zhabotinsky Equation with Caputo Operator" Axioms 12, no. 9: 825. https://doi.org/10.3390/axioms12090825

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop