Next Article in Journal
Self-Similar Growth and Synergistic Link Prediction in Technology-Convergence Networks: The Case of Intelligent Transportation Systems
Next Article in Special Issue
Higher-Order Dispersive and Nonlinearity Modulations on the Propagating Optical Solitary Breather and Super Huge Waves
Previous Article in Journal
On the Estimation of the Persistence Exponent for a Fractionally Integrated Brownian Motion by Numerical Simulations
Previous Article in Special Issue
Fractional Study of the Non-Linear Burgers’ Equations via a Semi-Analytical Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional View Study of the Brusselator Reaction–Diffusion Model Occurring in Chemical Reactions

by
Saleh Alshammari
1,
M. Mossa Al-Sawalha
1,* and
Jamal R. Humaidi
2
1
Department of Mathematics, College of Science, University of Ha’il, P.O. Box 2440, Ha’il 81451, Saudi Arabia
2
Department of Chemistry, Faculty of Science, University of Ha’il, P.O. Box 2440, Ha’il 81451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(2), 108; https://doi.org/10.3390/fractalfract7020108
Submission received: 30 November 2022 / Revised: 12 January 2023 / Accepted: 13 January 2023 / Published: 20 January 2023

Abstract

:
In this paper, we study a fractional Brusselator reaction–diffusion model with the help of the residual power series transform method. Specific reaction–diffusion chemical processes are modeled by applying the fractional Brusselator reaction–diffusion model. It should be mentioned that many problems in nonlinear science are characterized by fractional differential equations, where an unknown term occurs when a fractional-order derivative is operating on it. The analytic method of this problem is rarely discussed in the literature, despite numerous scholars having researched its application and usefulness. To validate our proposed method’s accuracy, we compare the numerical results of the residual power series transform method and the exact result with different fractional orders. The solution shows that the introduced approach is a good tool for solving linear and nonlinear fractional system differential equations. Finally, we provide two and three-dimensional graphical plots to support the impact of the fractional derivative on the behavior of the achieved profile results to the proposed equations.

1. Introduction

Fractional calculus, which includes fractional order derivatives and fractional integrals, has attracted interest in the areas of natural sciences and technology due to its vast implementations in electromagnetic waves, viscoelasticity, diffusion of biological population, nanotechnology, image and signal processing, random walk, control theory of dynamical systems, anomalous transport, and anomalous diffusion [1,2,3]. Nonlinear ordinary and partial differential equations (PDEs) are used to represent the majority of real-world scientific and engineering phenomena [4]. In the mathematical modeling of numerous natural phenomena, fractional nonlinear differential equations are particularly helpful. The essential quality of fractional derivatives is their nonlocality, which emphasizes the fact that the future state derives from the current state and all prior states. The study of nonlinear ordinary differential equations is required to observe the behavior of a system. However, evaluating the analytical solution for these types of situations is difficult. As a result, they are dealt with by utilizing approximate numeric techniques [5,6,7].
Chemical reactions play an essential function, as shown by the Brusselsator reaction–diffusion system. The reaction terms in this model come from the modeling of the chemical structures, which includes plasma and laser physics, enzyme reactions, and other processes. This nonlinear partial differential equation develops since chemical reactions involve two unstable intermediates with controllable output and input compounds. Ozone generation by atomic oxygen is one of the physical issues the Brusselsator reaction–diffusion system solves. Specific chemical reaction–diffusion processes are modeled using the fractional Brusselator reaction–diffusion system. It should be mentioned that many models in nonlinear science are characterized by fractional differential equations, where an unknown function occurs when a fractional-order derivative is operating on it. The approximate analysis of this system is rarely discussed in the literature, despite numerous scholars having researched its application and usefulness [8,9,10,11]. Recently, the development of numerical schemes for the numerical solutions of reaction–diffusion Brusselator systems has received considerable attention, such as the decomposition method [12,13], the finite-difference method [14], the Runge–Kutta–Chebyshev technique [15], the dual-reciprocity boundary element technique [16], the differential quadrature technique [17], the collocation technique applying the radial basis function [18], the homotopy perturbation method [19], the residual power series method [20] and the variational iteration method [21].
The residual power series method is a semianalytical tool for investigating many types of fuzzy, partial, ordinary, and fractional-order integrodifferential problems. It is an effective optimization method, since it delivers solutions in the closed type of available features. Various fractional differential and integral equations and fuzzy fractional differential equations have been effectively solved using the residual power series technique, for example, the Newell–Whitehead–Segel fractional order equations [22] and the fractional Fokker–Planck equations [23]. Researchers have coupled the fractional resonant Schrodinger equation [24], massive fractional Thirring and Kundu–Eckhaus problems [25], and particular categories of fractional fuzzy differential equations [26,27]. The Yang transform (YT), which Xiao-Jun Yang introduced, can be used to resolve a variety of differential equations with constant coefficients. This article aims to analyze the analytic and approximated solution for a nonlinear system’s fractional partial differential equations by applying the residual power series transform method. The residual power series transform method is a combination of the Yang transformation and the residual power series method, which provides the analytical results as a easily fractional power series solution by converting the suggested problems to Yang space and creating the results for the equations. Finally, the result can be achieved by the inverse Yang transform of the achieved solutions. Using the limit concept, the unknown coefficients in a Yang expansion can be determined, as opposed to the residual power series transform method, which relies on the fractional derivative and may require time to determine the fractional derivatives to identify the results. The residual power series transform method requires fewer processing resources and less time and provides more precision.
This work is arranged as follows: Section 2 contains some useful basic definitions of the Yang transform. The suggested method is explained in terms of achieving the results for the fractional partial differential equations schemes in Section 3. In Section 4, numerical solutions are presented for some problems. Lastly, concluding remarks are drawn in Section 2.

2. Preliminaries

This section is concerned with the fundamental concept of fractional calculus along with the Yang transform theory.

2.1. Definition

The fractional Caputo derivative is given as [28]
D ϑ U ( y , ϑ ) = 1 Γ ( k ) 0 ϑ ( ϑ ) k 1 U ( k ) ( y , ) d , k 1 < k , k N .

2.2. Definition

The Yang transform is stated as [28]
Y { U ( ϑ ) } = M ( u ) = 0 e ϑ u U ( ϑ ) d ϑ , ϑ > 0 , u ( ϑ 1 , ϑ 2 ) ,
and the inverse Yang transform is
Y 1 { M ( u ) } = U ( ϑ ) .

2.3. Definition

The nth derivative Yang transform is stated as [28]
Y { U n ( ϑ ) } = M ( u ) u n k = 0 n 1 U k ( 0 ) u n k 1 , n = 1 , 2 , 3 ,

2.4. Definition

The Yang transform in connection with the fractional-order derivative is stated as [28]
Y { U ( ϑ ) } = M ( u ) u k = 0 n 1 U k ( 0 ) u ( k + 1 ) , n 1 < n .

3. Methodology

In this section, we present the suggested methodology of the fractional partial differential equations [29]
D τ ω 1 ( φ , τ ) = 1 ω 1 , ω 2 , , ω m + f 1 ( φ , τ ) D τ ω 2 ( φ , τ ) = 2 ω 1 , ω 2 , , ω m + f 2 ( φ , τ ) D τ ω m ( φ , τ ) = m ω 1 , ω 2 , , ω m + f m ( φ , τ )
subject to:
ω i ( φ , 0 ) = g i ( φ ) , i = 1 , 2 , , m ,
where D τ refers to the CD of order , ( 0 , 1 ] , i , i = 1 , 2 , , m are well-known nonlinear analytic functions, and f i , i = 1 , 2 , , m can be represented as
Now, the technique for the RPSM solution of the fractional partial differential equations in (6) and (7) can be constructed through the following steps.
Step 1. We write the fractional nonlinear differential equation as in Equations (6) and (7) in the following short form:
D τ ω j ( φ , τ ) j ω 1 , ω 2 , , ω m f j ( φ , τ ) = 0 , j = 1 , 2 , , m ,
subject to:
ω j ( φ , 0 ) = g j ( φ ) , j = 1 , 2 , , m .
Step 2. Applying the Yang transform to both sides of Equation (8), it becomes:
W j ( φ , s ) s ω j ( φ , 0 ) s N j ( φ , s ) s F j ( φ , s ) = 0 ,
where W j ( φ , s ) = Y τ ω j ( φ , τ ) , N j ( φ , s ) = Y τ N j ω 1 , ω 2 , , ω m , and F j ( φ , s ) = Y τ f j ( φ , τ ) .
Step 3. Suppose that the solutions of fractional partial differential equations in Equation (10) are in the form:
W j ( φ , s ) = n = 0 s n + 1 g j n ( φ ) , s > 0 , j = 1 , 2 , , m .
Applying Equation (9), we have g j ( φ ) = ω j ( φ , 0 ) = lim s s W j ( φ , s ) = g j 0 ( φ ) . Therefore, the series form solutions can be represented as:
W j ( φ , s ) = s g j ( φ ) + n = 1 s n + 1 g j n ( φ ) , s > 0 , j = 1 , 2 , , m .
Step 4. We define the k th-truncated series of W j ( φ , s ) as follows:
W j k ( φ , s ) = s g j ( φ ) + n = 1 k s n φ + 1 g j n ( φ ) , s > 0 j = 1 , 2 , , m .
Step 5. We define the Yang residual functions of Equation (10) and the k th-Yang residual functions, respectively, as follows:
Y Res j ( φ , s ) = W j ( φ , s ) s ω j ( φ , 0 ) s N j ( φ , s ) s F j ( φ , s ) .
Y Res j k ( φ , s ) = W j k ( φ , s ) s ω j ( φ , 0 ) s N j ( φ , s ) s F j ( φ , s ) .
Step 6. We substitute the form of W j k ( φ , s ) as in Equation (13) into Equation (15).
Step 7. We multiply Equation (15) by s k a + 1 , k = 1 , 2 , , j = 1 , 2 , , m .
Step 8. We solve the following algebraic equations recursively:
lim s s k + 1 L Res j k ( φ , s ) = 0 , k = 1 , 2 , 3 , j = 1 , 2 , , m ;
then, we obtain the forms of the unknown coefficients, g j 1 ( φ ) , g j 2 ( φ ) , , g j ( k 1 ) ( φ ) , respectively.
Step 9. We substitute the obtained forms of g j i ( φ ) , i = 1 , 2 , , m to obtain the k t h series form solutions of Equation (10).
Step 10. We operate the inverse YT on the final form of W j k ( φ , s ) , j = 1 , 2 , , m to obtain the k t h series form solution, ω j k ( φ , τ ) , j = 1 , 2 , , m , of Equations (6) and (7).

4. Results

Example 1.
We considered the system of the fractional reaction–diffusion equation
D τ u = u u 2 u v + u φ φ , D τ v = v φ φ u v , 0 < 1 , τ > 0 ,
with the initial conditions
u ( φ , 0 ) = e p φ 1 + e 0.5 p φ 2 , v ( φ , 0 ) = 1 1 + e 0.5 p φ ,
where p is constant.
Using the Yang transformation on Equation (17) and applying Equation (18), we obtained
Y [ D τ u ] = Y [ u u 2 u v + u φ φ ] , Y [ D τ v ] = Y [ v φ φ u v ] ,
U ( φ , ϱ ) = ϱ u ( φ , 0 ) + ϱ Y [ u u 2 u v + u φ φ ] , V ( φ , ϱ ) = ϱ v ( φ , 0 ) + ϱ Y [ v φ φ u v ] .
The K t h -truncated term series were
U ( φ , ϱ ) = ϱ e p φ 1 + e 0.5 p φ 2 + n = 1 ß ϱ n + 1 f n ( φ ) , V ( φ , ϱ ) = ϱ 1 + e 0.5 p φ + n = 1 ß ϱ n + 1 g n ( φ ) ,
and the k t h -Yang residual functions were
Y τ R e s u , k ( φ , ϱ ) = U ß ( φ , ϱ ) ϱ e p φ 1 + e 0.5 p φ 2 ϱ Y τ [ Y τ 1 U ß ( φ , ϱ ) Y τ 1 U ß 2 ( φ , ϱ ) Y τ 1 U ß ( φ , ϱ ) Y τ 1 V ß ( φ , ϱ ) + Y τ 1 U ß , φ φ ( φ , ϱ ) ] , Y τ R e s u , k ( φ , ϱ ) = V ß ( φ , ϱ ) ϱ 1 + e 0.5 p φ ϱ Y τ Y τ 1 V ß , φ φ ( φ , ϱ ) Y τ 1 U ß ( φ , ϱ ) Y τ 1 V ß ( φ , ϱ ) .
Now, to find f ß ( φ ) , and g ß ( φ ) , k = 1 , 2 , 3 , , we substituted the k t h -truncated series Equation (21) into the k t h -Yang residual function Equation (22), multiplied the achieved solutions by ϱ k + 1 , and then recursively investigated the relations lim ϱ ( ϱ k + 1 Y τ R e s u , k ( φ , ϱ ) ) = 0 and lim ϱ ( ϱ k + 1 Y τ R e s v , k ( φ , ϱ ) ) = 0 , k = 1 , 2 , 3 , . The first few functions were
f 0 ( φ ) = e p φ 1 + e 0.5 p φ 2 , g 0 ( φ ) = 1 1 + e 0.5 p φ , f 1 ( φ ) = ( p 2 e 0.5 p φ 2 p 2 2 e 0.5 p φ ) e p φ ( 1 + e 0.5 p φ ) 4 , g 1 ( φ ) = 0.25 ( 4 e p φ + p 2 e 0.5 p φ p 2 e p x ) ( 1 + e 0.5 p φ ) 3 , f 2 ( φ ) = 1 8 ( 1 + e 0.5 p φ ) 6 ( 32 p 2 e 2 p φ + 16 e 2 p φ + 28 e 1.5 p φ p 2 33 p 4 e 1.5 p φ + 4 p 2 e 2.5 p φ + 18 p 4 e 2 p φ p 4 e 2.5 p φ + 8 p 4 e p φ ) , g 2 ( φ ) = 1 16 ( 1 + e 0.5 p φ ) 5 ( p 4 e 2 p φ + 16 e 2 p φ 8 p 2 e 2 p φ + 11 p 2 e p φ 11 p 4 e 1.5 p φ + 40 p 2 e 1.5 p φ p 4 e 0.5 p φ 32 p 2 e p φ 16 e 1.5 p φ ) .
Putting the values of f ß ( φ ) and g ß ( φ ) , k = 1 , 2 , 3 , into Equation (21), we obtained
U ( φ , ϱ ) = ϱ e p φ 1 + e 0.5 p φ 2 + ϱ + 1 ( p 2 e 0.5 p φ 2 p 2 2 e 0.5 p φ ) e p φ ( 1 + e 0.5 p φ ) 4 + ( ϱ 2 + 1 8 ( 1 + e 0.5 p φ ) 6 ( 32 p 2 e 2 p φ + 16 e 2 p φ + 28 e 1.5 p φ p 2 33 p 4 e 1.5 p φ + 4 p 2 e 2.5 p φ + 18 p 4 e 2 p φ p 4 e 2.5 p φ + 8 p 4 e p φ ) ) , V ( φ , ϱ ) = ϱ 1 + e 0.5 p φ ϱ + 1 0.25 ( 4 e p φ + p 2 e 0.5 p φ p 2 e p x ) ( 1 + e 0.5 p φ ) 3 + ( ϱ 2 + 1 16 ( 1 + e 0.5 p φ ) 5 ( p 4 e 2 p φ + 16 e 2 p φ 8 p 2 e 2 p φ + 11 p 2 e p φ 11 p 4 e 1.5 p φ + 40 p 2 e 1.5 p φ p 4 e 0.5 p φ 32 p 2 e p φ 16 e 1.5 p φ ) ) .
Using the inverse Yang transform, we obtained
u ( φ , τ ) = e p φ 1 + e 0.5 p φ 2 + ( p 2 e 0.5 p φ 2 p 2 2 e 0.5 p φ ) e p φ ( 1 + e 0.5 p φ ) 4 τ Γ ( + 1 ) + ( 1 8 ( 1 + e 0.5 p φ ) 6 ( 32 p 2 e 2 p φ + 16 e 2 p φ + 28 e 1.5 p φ p 2 33 p 4 e 1.5 p φ + 4 p 2 e 2.5 p φ + 18 p 4 e 2 p φ p 4 e 2.5 p φ + 8 p 4 e p φ ) ) τ 2 Γ ( 2 + 1 ) , v ( φ , τ ) = 1 1 + e 0.5 p φ 0.25 ( 4 e p φ + p 2 e 0.5 p φ p 2 e p x ) ( 1 + e 0.5 p φ ) 3 τ Γ ( + 1 ) + ( 1 16 ( 1 + e 0.5 p φ ) 5 ( p 4 e 2 p φ + 16 e 2 p φ 8 p 2 e 2 p φ + 11 p 2 e p φ 11 p 4 e 1.5 p φ + 40 p 2 e 1.5 p φ p 4 e 0.5 p φ 32 p 2 e p φ 16 e 1.5 p φ ) ) τ 2 Γ ( 2 + 1 ) .
The residual power series transformation method was implemented to achieve the analytical and exact results for the strongly nonlinear fractional reaction–diffusion equation with the initial conditions. In Figure 1, the analytical solutions of (a) u ( φ , ( τ ) ) and (b) v ( φ , ( τ ) ) of Example 1 are shown. In Figure 2, (a) the fractional order 0.5 , (b) 0.7 , (c) the actual result and (d) analytical solution of Eaxample 1 of u ( φ , τ ) . Similarly, in Figure 3, (a) the fractional order 0.5 , (b) 0.7 , (c) the actual result and (d) analytical solution of Eaxample 1 of v ( φ , τ ) . In the above, the figures showed that in the residual power series transformation method, the analytical results were close to each other. Table 1 and Table 2 detail the numerical analysis of u ( φ , τ ) and v ( φ , τ ) of the given example from Example 1, with different values of φ, τ, and p = 2 3 .
Example 2.
We considered the fractional order reaction–diffusion Brusselator equation
D τ u ( φ , ρ , τ ) = u 2 ( φ , ρ , τ ) v ( φ , ρ , τ ) u ( φ , ρ , τ ) ( A + 1 ) + 1 500 u φ φ ( φ , ρ , τ ) + u ρ ρ ( φ , ρ , τ ) + B , D τ v ( φ , ρ , τ ) = u 2 ( φ , ρ , τ ) v ( φ , ρ , τ ) + A u ( φ , ρ , τ ) + 1 500 u φ φ ( φ , ρ , τ ) + u ρ ρ ( φ , ρ , τ ) ,
with the initial conditions
u ( φ , ρ , 0 ) = 2 + 1 4 ρ , v ( φ , ρ , 0 ) = 1 + 4 5 φ ,
where u ( φ , ρ , τ ) and v ( φ , ρ , τ ) represent the chemical concentrations of the intermediate reaction products, and A and B are the constant concentrations of the input reagents, where A = 17 5 and B = 1 .
Using the Yang transformation on Equation (26) and applying Equation (27), we obtained
Y [ D τ u ( φ , ρ , τ ) ] = Y [ u 2 ( φ , ρ , τ ) v ( φ , ρ , τ ) u ( φ , ρ , τ ) ( A + 1 ) + 1 500 u φ φ ( φ , ρ , τ ) + u ρ ρ ( φ , ρ , τ ) + B ] , Y [ D τ v ( φ , ρ , τ ) ] = Y [ u 2 ( φ , ρ , τ ) v ( φ , ρ , τ ) + A u ( φ , ρ , τ ) + 1 500 u φ φ ( φ , ρ , τ ) + u ρ ρ ( φ , ρ , τ ) ] ,
U ( φ , ρ , ϱ ) = ϱ u ( φ , ρ , 0 ) + ϱ Y [ u 2 ( φ , ρ , τ ) v ( φ , ρ , τ ) u ( φ , ρ , τ ) ( A + 1 ) + 1 500 u φ φ ( φ , ρ , τ ) + u ρ ρ ( φ , ρ , τ ) + B ] , V ( φ , ρ , ϱ ) = ϱ v ( φ , ρ , 0 ) + ϱ Y [ u 2 ( φ , ρ , τ ) v ( φ , ρ , τ ) + A u ( φ , ρ , τ ) + 1 500 u φ φ ( φ , ρ , τ ) + u ρ ρ ( φ , ρ , τ ) ] .
The k t h truncate series terms were
U ( φ , ϱ ) = ϱ ( 2 + 1 4 ρ ) + n = 1 ß ϱ n + 1 f n ( φ , ρ ) , V ( φ , ϱ ) = ϱ ( 1 + 4 5 φ ) + n = 1 ß ϱ n + 1 g n ( φ , ρ ) ,
and the k t h -Yang residual functions were
Y τ R e s u , k ( φ , ρ , ϱ ) = U ß ( φ , ρ , ϱ ) ϱ 2 + 1 4 ρ 1 ϱ Y τ [ Y τ 1 U ß 2 ( φ , ρ , ϱ ) Y τ 1 V ß ( φ , ρ , ϱ ) Y τ 1 U ß ( φ , ρ , ϱ ) ( A + 1 ) + 1 500 Y τ 1 U ß , φ φ ( φ , ρ , ϱ ) + Y τ 1 U ß , ρ ρ ( φ , ρ , ϱ ) + B ] , Y τ R e s v , k ( φ , ρ , ϱ ) = V ß ( φ , ρ , ϱ ) ϱ 1 + 4 5 φ 1 ϱ Y τ [ Y τ 1 U ß 2 ( φ , ρ , ϱ ) Y τ 1 V ß ( φ , ρ , ϱ ) + A Y τ 1 U ß ( φ , ρ , ϱ ) + 1 500 Y τ 1 U ß , φ φ ( φ , ρ , ϱ ) + Y τ 1 U ß , ρ ρ ( φ , ρ , ϱ ) ] .
Now, to calculate f ß ( φ , ρ ) and g ß ( φ , ρ ) , k = 1 , 2 , 3 , , we put the k t h -truncate series (Equation (30)) into the k t h -Yang residual function (Equation (31)), multiplied the solution equation by ϱ k + 1 , and then recursively analyzed the relations lim ϱ ( ϱ k + 1 Y τ R e s u , k ( φ , ψ , ϱ ) ) = 0 , and lim ϱ ( ϱ k + 1 Y τ R e s v , k ( φ , ψ , ϱ ) ) = 0 , k = 1 , 2 , 3 , . The first few terms were defined as:
f 0 ( φ , ρ ) = 2 + ρ 4 , g 0 ( φ , ρ ) = 1 + 4 φ 5 , f 1 ( φ , ρ ) = 1899 500 + 16 φ 5 ρ 10 + 4 φ ρ 5 + 1 ρ 2 16 + ρ 2 φ 20 , g 1 ( φ , ρ ) = 1401 500 16 φ 5 3 ρ 20 4 ρ φ 5 ρ 2 16 ρ 2 φ 20 , f 2 ( φ , ρ ) = 66261 5000 36667 φ 1250 + 256 ρ 2 251 13 ρ 3 320 ρ 4 256 ρ 3 φ 20 + ρ 2 + 96 φ 2 ρ 25 + 16 φ 5 + 12 φ 2 ρ 2 25 + ρ 3 φ 2 50 ρ 4 φ 320 33 ρ 2 φ 50 9199 ρ φ 1250 2 5 2399 ρ 2 8000 157 ρ 1000 + 2 ρ φ 5 , g 2 ( φ , ρ ) = 44021 5000 + ρ 3 φ 20 96 ρ φ 2 25 12 ρ 2 φ 2 25 + ρ 4 φ 320 ρ 3 φ 2 50 + 28917 φ 1250 256 φ 2 25 + 13 ρ 3 320 + ρ 4 256 + 1899 ρ 2 8000 + 7699 ρ φ 1250 + 61 ρ 2 φ 100 243 ρ 1000 , .
Putting the value of f ß ( φ , ρ ) and g ß ( φ , ρ ) , k = 1 , 2 , 3 , , into Equation (30), we obtained
U ( φ , ρ , ϱ ) = ϱ ( 2 + 1 4 ρ ) + 1899 500 + 16 φ 5 ρ 10 + 4 φ ρ 5 + 1 ρ 2 16 + ρ 2 φ 20 ϱ + 1 + [ 66261 5000 36667 φ 1250 + 256 ρ 2 251 13 ρ 3 320 ρ 4 256 ρ 3 φ 20 + ρ 2 + 96 φ 2 ρ 25 + 16 φ 5 + 12 φ 2 ρ 2 25 + ρ 3 φ 2 50 ρ 4 φ 320 33 ρ 2 φ 50 9199 ρ φ 1250 2 5 2399 ρ 2 8000 157 ρ 1000 + 2 ρ φ 5 ] ϱ 2 + 1 + , V ( φ , ρ , ϱ ) = ϱ ( 1 + 4 5 φ ) + 1401 500 16 φ 5 3 ρ 20 4 ρ φ 5 ρ 2 16 ρ 2 φ 20 ϱ + 1 + [ 44021 5000 + ρ 3 φ 20 96 ρ φ 2 25 12 ρ 2 φ 2 25 + ρ 4 φ 320 ρ 3 φ 2 50 + 28917 φ 1250 256 φ 2 25 + 13 ρ 3 320 + ρ 4 256 + 1899 ρ 2 8000 + 7699 ρ φ 1250 + 61 ρ 2 φ 100 243 ρ 1000 ] ϱ 2 + 1 + .
Now, applying the inverse Yang transform, we obtained
u ( φ , ρ , τ ) = 2 + 1 4 ρ + 1899 500 + 16 φ 5 ρ 10 + 4 φ ρ 5 + ρ 2 16 + ρ 2 φ 20 τ Γ ( + 1 ) + [ 66261 5000 36667 φ 1250 + 256 ρ 2 251 13 ρ 3 320 ρ 4 256 ρ 3 φ 20 + ρ 2 + 96 φ 2 ρ 25 + 16 φ 5 + 12 φ 2 ρ 2 25 + ρ 3 φ 2 50 ρ 4 φ 320 33 ρ 2 φ 50 9199 ρ φ 1250 2 5 2399 ρ 2 8000 157 ρ 1000 + 2 ρ φ 5 ] τ 2 Γ ( 2 + 1 ) + , v ( φ , ρ , τ ) = 1 + 4 5 φ + 1401 500 16 φ 5 3 ρ 20 4 ρ φ 5 ρ 2 16 ρ 2 φ 20 τ Γ ( + 1 ) + [ 44021 5000 + ρ 3 φ 20 96 ρ φ 2 25 12 ρ 2 φ 2 25 + ρ 4 φ 320 ρ 3 φ 2 50 + 28917 φ 1250 256 φ 2 25 + 13 ρ 3 320 + ρ 4 256 + 1899 ρ 2 8000 + 7699 ρ φ 1250 + 61 ρ 2 φ 100 243 ρ 1000 ] τ 2 Γ ( 2 + 1 ) + .
The residual power series transformation method was implemented to achieve the analytical and exact results of the strongly nonlinear fractional reaction–diffusion equation with the initial conditions. In Figure 4, the analytical solutions of (a) u ( φ , ( τ ) ) and (b) v ( φ , ( τ ) ) of Example 2 are shown. In Figure 5, (a) the fractional order 0.5 , (b) 0.7 , (c) the actual result and (d) analytical solution of Eaxample 2 of u ( φ , τ ) . Similarly, in Figure 6, (a) the fractional order 0.5 , (b) 0.7 , (c) the actual result and (d) analytical solution of Example 2 of v ( φ , τ ) . These figures showed that in using the residual power series transformation method, the analytical results were close to each other. Table 3 and Table 4 detail the numerical analysis of u ( φ , τ ) and v ( φ , τ ) from Example 2 with different values of φ, τ, and p = 2 3 .

5. Conclusions

This paper examined the fractional reaction–diffusion and reaction–diffusion Brusselator equations in terms of the Caputo fractional derivative. With the aid of the RPSM, the series solution of the investigated model was successfully attained. Plots were made of the numerical simulations of the proposed solution using various fractional values of . The technique was applied to the suggested fractional-order systems to investigate its analytical solutions. The solutions showed that the method was inexpensive and suitable for achieving analytical solutions. We, therefore, conclude that the suggested fractional models connected to the Caputo fractional derivative are very helpful in efficiently investigating the problems arising in science and engineering.

Author Contributions

Conceptualization, S.A. and M.M.A.-S.; methodology, J.R.H.; software, S.A.; validation, M.M.A.-S.; formal analysis, J.R.H.; investigation, S.A; resources, M.M.A.-S.; data curation, J.R.H.; writing—original draft preparation, M.M.A.-S.; writing—review and editing, J.R.H.; visualization, S.A.; supervision, M.M.A.-S.; project administration, J.R.H.; funding acquisition, S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was funded by the Deputy for Research & Innovation, the Ministry of Education, through the Initiative of Institutional Funding at the University of Ha’il-Saudi Arabia through project number IFP-22 127.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Magin, R. Fractional calculus in bioengineering, part 1. Crit. Rev. Biomed. Eng. 2004, 32, 104. [Google Scholar]
  2. Machado, J.T.; Kiryakova, V.; Mainardi, F. Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1140–1153. [Google Scholar] [CrossRef] [Green Version]
  3. Sabatier, J.A.T.M.J.; Agrawal, O.P.; Machado, J.T. Advances in Fractional Calculus; Springer: Dordrecht, The Netherlands, 2007; Volume 4, p. 9. [Google Scholar]
  4. Giusti, A.; Colombaro, I.; Garra, R.; Garrappa, R.; Polito, F.; Popolizio, M.; Mainardi, F. A practical guide to Prabhakar fractional calculus. Fract. Calc. Appl. Anal. 2020, 23, 9–54. [Google Scholar] [CrossRef] [Green Version]
  5. Baleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus: Models and Numerical Methods; World Scientific: Singapore, 2012; Volume 3. [Google Scholar]
  6. Debnath, L. Recent applications of fractional calculus to science and engineering. Int. J. Math. Math. Sci. 2003, 2003, 3413–3442. [Google Scholar] [CrossRef] [Green Version]
  7. Axtell, M.; Bise, M.E. Fractional calculus application in control systems. In Proceedings of the IEEE Conference on Aerospace and Electronics, Dayton, OH, USA, 21–25 May 1990; pp. 563–566. [Google Scholar]
  8. El-Saka, H.A.A.; Arafa, A.A.M.; Gouda, M.I. Dynamical analysis of a fractional SIRS model on homogenous networks. Adv. Differ. Equ. 2019, 2019, 1–15. [Google Scholar] [CrossRef] [Green Version]
  9. Alyousef, H.A.; Yasmin, H.; Shah, N.A.; El-Sherif, L.S.; El-Tantawy, S.A. Mathematical Modeling and Analysis of the Steady Electro-Osmotic Flow of Two Immiscible Fluids: A Biomedical Application. Coatings 2023, 13, 115. [Google Scholar] [CrossRef]
  10. Alyousef, H.A.; Nonlaopon, K.; El-Sherif, L.S.; El-Tantawy, S.A. An Efficient Analytical Method for Analyzing the Nonlinear Fractional Klein-Fock-Gordon Equations. Symmetry 2022, 14, 2640. [Google Scholar] [CrossRef]
  11. Singh, J.; Rashidi, M.M.; Kumar, D.; Swroop, R. A fractional model of a dynamical Brusselator reaction-diffusion system arising in triple collision and enzymatic reactions. Nonlinear Eng. 2016, 5, 277–285. [Google Scholar] [CrossRef]
  12. Adomian, G. The diffusion-Brusselator equation. Comput. Math. Appl. 1995, 29, 1–3. [Google Scholar] [CrossRef] [Green Version]
  13. Wazwaz, A.M. The decomposition method applied to systems of partial differential equations and to the reaction-diffusion Brusselator model. Appl. Math. Comput. 2000, 110, 251–264. [Google Scholar] [CrossRef]
  14. Twizell, E.H.; Gumel, A.B.; Cao, Q. A second-order scheme for the Brusselator reaction-diffusion system. J. Math. Chem. 1999, 26, 297–316. [Google Scholar] [CrossRef]
  15. Verwer, J.G.; Hundsdorfer, W.H.; Sommeijer, B.P. Convergence properties of the Runge-Kutta-Chebyshev method. Numer. Math. 1990, 57, 157–178. [Google Scholar] [CrossRef] [Green Version]
  16. Ang, W.T. The two-dimensional reaction-diffusion Brusselator system: A dual-reciprocity boundary element solution. Eng. Anal. Bound. Elem. 2003, 27, 897–903. [Google Scholar] [CrossRef]
  17. Mittal, R.C.; Jiwari, R. Numerical study of two-dimensional reaction-diffusion Brusselator system by differential quadrature method. Int. J. Comput. Methods Eng. Sci. Mech. 2011, 12, 14–25. [Google Scholar] [CrossRef]
  18. Ali, A.; Haq, S. A computational modeling of the behavior of the two-dimensional reaction-diffusion Brusselator system. Appl. Math. Model. 2010, 34, 3896–3909. [Google Scholar]
  19. Kumar, S.; Khan, Y.; Yildirim, A. A mathematical modeling arising in the chemical systems and its approximate numerical solution. Asia-Pac. J. Chem. Eng. 2012, 7, 835–840. [Google Scholar] [CrossRef]
  20. 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]
  21. Jafari, H.; Kadem, A.; Baleanu, D. Variational iteration method for a fractional-order Brusselator system. Abstr. Appl. Anal. 2014, 2014, 496323. [Google Scholar] [CrossRef] [Green Version]
  22. Saadeh, R.; Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Salma Din, U.K. Application of fractional residual power series algorithm to solve Newell-Whitehead-Segel equation of fractional order. Symmetry 2019, 11, 1431. [Google Scholar] [CrossRef] [Green Version]
  23. Freihet, A.; Hasan, S.; Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Salma Din, U.K. Toward computational algorithm for time-fractional Fokker-Planck models. Adv. Mech. Eng. 2019, 11, 1687814019881039. [Google Scholar] [CrossRef]
  24. Al-Smadi, M.; Arqub, O.A.; Momani, S. Numerical computations of coupled fractional resonant Schrodinger equations arising in quantum mechanics under conformable fractional derivative sense. Phys. Scr. 2020, 95, 075218. [Google Scholar] [CrossRef]
  25. Al-Smadi, M.; Arqub, O.A.; Hadid, S. Approximate solutions of nonlinear fractional Kundu-Eckhaus and coupled fractional massive Thirring equations emerging in quantum field theory using conformable residual power series method. Phys. Scr. 2020, 95, 105205. [Google Scholar] [CrossRef]
  26. Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Salma Din, U.K. Computational optimization of residual power series algorithm for certain classes of fuzzy fractional differential equations. Int. J. Differ. Equ. 2018, 2018, 8686502. [Google Scholar] [CrossRef]
  27. Aljahdaly, N.H.; Agarwal, R.P.; Botmart, T. The analysis of the fractional-order system of third-order KdV equation within different operators. Alex. Eng. J. 2022, 61, 11825–11834. [Google Scholar] [CrossRef]
  28. Dattu, M.K.U. New integral transform: Fundamental properties, investigations and applications. Iaetsd J. Adv. Res. Appl. Sci. 2018, 5, 534–539. [Google Scholar]
  29. Oqielat, M.A.N.; Eriqat, T.; Al-Zhour, Z.; Ogilat, O.; El-Ajou, A.; Hashim, I. Construction of fractional series solutions to nonlinear fractional reaction-diffusion for bacteria growth model via Laplace residual power series method. Int. J. Dyn. Control. 2022, 2022, 1–8. [Google Scholar] [CrossRef]
Figure 1. The analytical solution for u ( φ , ( τ ) ) and v ( φ , ( τ ) ) from Example 1.
Figure 1. The analytical solution for u ( φ , ( τ ) ) and v ( φ , ( τ ) ) from Example 1.
Fractalfract 07 00108 g001
Figure 2. The actual and analytical results for u ( φ , τ ) at k = 2 3 , c = 1 , and different values of from Example 1.
Figure 2. The actual and analytical results for u ( φ , τ ) at k = 2 3 , c = 1 , and different values of from Example 1.
Fractalfract 07 00108 g002
Figure 3. The actual and analytical solutions for v ( φ , τ ) at p = 2 3 , c = 1 , and different values of from Example 1.
Figure 3. The actual and analytical solutions for v ( φ , τ ) at p = 2 3 , c = 1 , and different values of from Example 1.
Fractalfract 07 00108 g003
Figure 4. The analytical solution (a) u ( φ , ( τ ) ) and subplot (b) v ( φ , ( τ ) ) at A = 17 5 , B = 1 , x = 0.1 , and y = 0.1 , respectively, from Example 2.
Figure 4. The analytical solution (a) u ( φ , ( τ ) ) and subplot (b) v ( φ , ( τ ) ) at A = 17 5 , B = 1 , x = 0.1 , and y = 0.1 , respectively, from Example 2.
Fractalfract 07 00108 g004
Figure 5. The analytical solution for u ( φ , τ ) at A = 17 5 , B = 1 , and y = 1 for distinct values of from Example 2.
Figure 5. The analytical solution for u ( φ , τ ) at A = 17 5 , B = 1 , and y = 1 for distinct values of from Example 2.
Fractalfract 07 00108 g005
Figure 6. The analytical solution for v ( φ , τ ) at A = 17 5 , B = 1 , and y = 1 for different values of from Example 2.
Figure 6. The analytical solution for v ( φ , τ ) at A = 17 5 , B = 1 , and y = 1 for different values of from Example 2.
Fractalfract 07 00108 g006
Table 1. Numerical analysis of u ( φ , τ ) of the example given in Example 1 with different values of φ , τ , and p = 2 3 .
Table 1. Numerical analysis of u ( φ , τ ) of the example given in Example 1 with different values of φ , τ , and p = 2 3 .
τ φ AE at = 0.5 AE at = 0.75 AE at = 1 AE at = 1 [20]
0.20.08056294520.05022040620.03250208750.0325020875
0.40.08254547180.05153310390.03342417340.0334241734
0.10.60.08391118440.05247875870.03411705700.0341170570
0.80.08464373830.05304486060.03457064780.0345706478
10.08474527710.05322952740.03478137640.0347813764
0.20.13861503910.10640554040.08162947030.0816294703
0.40.14210278460.10918052840.08385714300.0838571430
0.250.60.14450635720.11116238910.08550122050.0855012205
0.80.14579629270.11232514540.08653873090.0865387309
10.14597646950.11266625590.08696338340.0869633834
Table 2. Numerical analysis of v ( φ , τ ) of the example given in Example 1 with different values of φ , τ , and p = 2 3 .
Table 2. Numerical analysis of v ( φ , τ ) of the example given in Example 1 with different values of φ , τ , and p = 2 3 .
τ φ AEat = 0.5 AEat = 0.75 AEat = 1 AEat = 1
0.20.03654093700.01626475170.00444859400.0044485940
0.40.03723492640.01675765850.00474332640.0047433264
0.10.60.03781619280.01720175440.00502253920.0050225392
0.80.03828028950.01759336080.00528386290.0052838629
10.03862422760.01792944690.00552517280.0055251728
0.20.04938018870.02771027020.01098974760.0109897476
0.40.05021671210.02856741550.01168086540.0116808654
0.250.60.05087906420.02933019370.01233031040.0123303104
0.80.05136288570.02999256190.01293267540.0129326754
10.05166622230.03054981270.01348322050.0134832205
Table 3. Numerical analysis of u ( φ , τ ) of the example given in 2 with different values of φ , τ , and p = 2 3 .
Table 3. Numerical analysis of u ( φ , τ ) of the example given in 2 with different values of φ , τ , and p = 2 3 .
τ φ Numerical Simulation at = 0.5 Numerical Simulation at = 0.75 Numerical Simulation at = 1
0.21.8377073731.9170525471.980814044
0.41.6255769791.9757016572.036756119
0.10.61.5300865862.0453774872.098530194
0.81.5512361932.1260800362.166136269
11.6890258002.2178093052.239574344
0.21.7135778151.6748590931.702150273
0.41.5539286911.6701899501.748038242
0.250.61.5039580221.7228978761.830376211
0.81.5636658091.8329828691.949164180
11.7330520512.0004449312.104402148
Table 4. Numerical analysis of v ( φ , τ ) of the example given in 2 with different values of φ , τ , and p = 2 3 .
Table 4. Numerical analysis of v ( φ , τ ) of the example given in 2 with different values of φ , τ , and p = 2 3 .
τ φ Numerical Simulation at = 0.5 Numerical Simulation at = 0.75 Numerical Simulation at = 1
0.21.4592388981.4244478741.321163456
0.41.7543692911.5286224551.421371381
0.10.61.9328596841.6050503371.515747306
0.81.9947100771.6537315231.604291231
11.9399204711.6746660091.687003156
0.21.3246481761.5288612541.499959102
0.42.0880508641.7761060481.590008633
0.250.62.5598535521.9136723871.643608164
0.82.7400562381.9415602701.660757695
12.6286589251.8597696981.641457227
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

Alshammari, S.; Al-Sawalha, M.M.; Humaidi, J.R. Fractional View Study of the Brusselator Reaction–Diffusion Model Occurring in Chemical Reactions. Fractal Fract. 2023, 7, 108. https://doi.org/10.3390/fractalfract7020108

AMA Style

Alshammari S, Al-Sawalha MM, Humaidi JR. Fractional View Study of the Brusselator Reaction–Diffusion Model Occurring in Chemical Reactions. Fractal and Fractional. 2023; 7(2):108. https://doi.org/10.3390/fractalfract7020108

Chicago/Turabian Style

Alshammari, Saleh, M. Mossa Al-Sawalha, and Jamal R. Humaidi. 2023. "Fractional View Study of the Brusselator Reaction–Diffusion Model Occurring in Chemical Reactions" Fractal and Fractional 7, no. 2: 108. https://doi.org/10.3390/fractalfract7020108

APA Style

Alshammari, S., Al-Sawalha, M. M., & Humaidi, J. R. (2023). Fractional View Study of the Brusselator Reaction–Diffusion Model Occurring in Chemical Reactions. Fractal and Fractional, 7(2), 108. https://doi.org/10.3390/fractalfract7020108

Article Metrics

Back to TopTop