Next Article in Journal
Analysis of the Fractional HIV Model under Proportional Hadamard-Caputo Operators
Next Article in Special Issue
Fractional Model of Multiple Trapping with Charge Leakage: Transient Photoconductivity and Transit–Time Dispersion
Previous Article in Journal
An Intelligent Model for Supercapacitors with a Graphene-Based Fractal Electrode to Investigate the Cyclic Voltammetry
Previous Article in Special Issue
Noise Spectral of GML Noise and GSR Behaviors for FGLE with Random Mass and Random Frequency
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Laplace Residual Power Series Method in Solving Coupled Fractional Neutron Diffusion Equations with Delayed Neutrons System

College of Science, Zarqa University, Zarqa 13110, Jordan
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(3), 219; https://doi.org/10.3390/fractalfract7030219
Submission received: 25 December 2022 / Revised: 23 January 2023 / Accepted: 31 January 2023 / Published: 27 February 2023
(This article belongs to the Special Issue Application of Fractional-Calculus in Physical Systems)

Abstract

:
In this paper, a system of coupled fractional neutron diffusion equations with delayed neutrons was solved efficiently by using a combination of residual power series and Laplace transform techniques, and the anomalous diffusion was considered by taking the non-Gaussian case with different values of fractional parameter α. The Laplace residual power series method (LRPSM) does not require differentiation, conversion, or discretization for the assumed conditions, so the approach is simple and suitable for solving higher-order fractional differential equations. To assure the theoretical results, two different neutron flux initial conditions were presented numerically, where the needed Mathematica codes were performed using essential nuclear reactor cross-section data, and the results for different values of times were tabulated and graphically figured out. Finally, it must be noted that the results align with the Adomian decomposition method.

1. Introduction

The revolution of using energy generated from nuclear reactors can be considered as an essential part of a clean renewable source of energy; today, nuclear reactors cover 10% of the use of the global needed energy.
The inspiring idea is that the creation of nuclear energy is released when heavy nuclei undergo fission while they are bombarded by (thermal) slow neutrons. The neutron diffusion equation, which is simplified from the neutron transport equation when using Fick’s law, describes the behavior of neutrons in the nuclear fission reactor [1,2,3]. Researchers have investigated their scientific ambitions by studying different issues in this field, where different cases of neutron diffusion equations have been solved by various techniques. The homotopy perturbation method (HPM) has been used in solving bare hemispherical symmetry and cylindrical symmetry one group neutron diffusion equation reactors [4,5]; the residual power series method was used to solve multi-energy groups of neutron diffusion equations [6,7]; then the HPM is used to solve reflected reactor one group different geometries neutron diffusion equations; finally, solving a multi-group reflected spherical reactor system of equations using the HPM is performed [8].
The distribution of neutrons in nuclear reactors is constructed in the same way as dissolving the solute in the chemical solutions, which can be represented mathematically by the continuity equation:
1 v c ψ r , t t = s r , t Σ a ψ r , t · J r , t ,
where ψ r , t is the neutron flux; J r , t is neutron current density; a is the macroscopic absorption cross-section; v c   is neutron velocity; and s r , t is a source of neutrons term. Using Fick’s law, which gives the relation between neutron flux ψ r ,   t and the current density J r , t :
J r ,   t = D ψ r ,   t
where D is the neutron diffusion coefficient, then the neutron diffusion equation can be expressed as
1 v c ψ r , t t = s r ,   t Σ a   ψ r , t + D 2 ψ r ,   t .
This is the well-known nuclear reactor time-dependent neutron diffusion equation.
One step forward, the fractional study of the time-dependent neutron diffusion equation is to make the situation more detailed by studying the physical effect of delayed neutrons. These delayed neutrons have an essential role in the reactor. Investigating a system of coupled fractional neutron diffusion equations with delayed neutrons was the target of this work.
The importance of delayed neutrons in reactors is the effective lifetime of such neutrons (which are given by their prompt lifetime, and additional delayed time), which is characterized by the β-decay of their precursor, which is considerably longer than the prompt neutron lifetime. Hence, delayed neutrons substantially increase the time constant of a reactor so that effective control is possible.
The contribution of the effect of the delayed neutron is included in Equation (3) in the source term s r , t , which can be written as
s r , t = γ f   ψ r , t + λ c r , t
In this study, the general three-dimension vector r will be simplified to the one-dimension vector x to be more suitable.
The system of time-dependent neutron diffusion equations is known as the neutron point kinetic model; this is in one dimension:
1 v c t ψ x , t = D 2 x 2 ψ x , t + a + γ f ψ x , t + λ c x , t ,                      
t c x , t = β ν f   ψ x , t λ c x , t ,
with initial conditions:
ψ x , 0 = ψ 0 x ,     c x , 0 = β γ f λ ψ 0 x .
where the delayed neutron density is c x , t , and the constants γ, λ, and β are the average number of neutrons produced per fission, the radioactive decay constant, and the fraction of the delayed fission neutrons, respectively, where a + γ f can be considered as the total macroscopic cross-section .
The kinetic equation describes the dynamics of space-time neutron distribution in the reactors. These equations are divided into point kinetic equations and space kinetic equations; the time change in neutron density is described by the kinetic point equations, which estimate the reactor power response; the fractional calculus governs an important part of the time-dependent diffusion equation in the presence of delayed neutrons, which is known as anomalous diffusion; the anomalous diffusion is called by the variation from the Gaussian distribution of the diffusion.
Anomalous diffusion is characterized by the asymptotical nonlinear (especially power-law type) dependence of the mean-squared displacement on the time, which is x 2 t ~ t α .
The parameter α feeds back in the determination of the type of fractional neutron diffusion equation, where the anomalous diffusion ( α 1 ) is classified into sub diffusion ( α < 1 ) and super diffusion ( α > 1 ).
In the case of α = 1: the diffusion is Gaussian, mathematically representing the non-fractional neutron diffusion equation.
The trapping of neutrons in the medium causes slower neutron movement considered as sub diffusion, while the fast movement of neutrons when they move in one direction for a very long time without any collision is characterized by super diffusion [9,10].
The fractional neutron diffusion equation with a delayed neutron system has been studied in many works where the fractional two energy group matrix representation for nuclear reactor dynamics with an external source is considered [11], while the coupled fractional neutron diffusion equations with delayed neutrons have been solved [12], and an exact solution in the case of the one-dimensional neutron diffusion kinetic equation with one delayed precursor concentration in Cartesian geometry was studied [13].
Here, the fractional neutron diffusion equations with one delayed neutron group can be written as follows:
1 v c D t α ψ x , t = D 2 ψ x , t + γ f a   ψ x , t + λ   c x , t ,
D t α c x , t = β γ f   ψ x , t λ   c x , t ,
ψ x , 0 = ψ 0 x , c x , 0 = β γ f λ ψ 0 x .
Numerous numerical and analytical techniques including Laplace transforms, the homotopy method, the Pade-variational iteration method, the finite difference method, the Adomian decomposition method, and others [14,15,16,17,18,19,20,21,22,23,24] have been developed to solve fractional differential equations. However, these methods have limitations such as the need for long and complex computations. Therefore, they need to be connected to a transform operator. The authors of [25] combined the residual power series method (RPSM) with the Laplace transform to introduce actual and series solutions of the linear and nonlinear FDEs. The new method, known as the Laplace-residual power series method (LRPSM), was adopted to construct series solutions of a variety of fractional deferential equations (FDEs). LRPSM does not rely on fractional derivation to determine the coefficients of the series as in RPSM, but depends on the concept of the limit, so few calculations generate the coefficients compared to the residual power series method. The current technique is quick, requires little computer memory, and is not influenced by computational round-off errors. Furthermore, this technique computes the coefficients of the power series using a chain of equations with more than one variable, indicating that the present method has a rapid convergence.
This paper solved a system of coupled fractional neutron diffusion equations with delayed neutrons efficiently by applying LRPSM. The accuracy and effectiveness of the method were clarified by displaying numerical examples and comparing the solutions with the results of some approved approaches.
The rest of paper is structured as follows. The concepts and theorems in fractional calculus and Laplace transforms that are required to reach our results are introduced in Section 2. The LRPSM algorithm is described in Section 3. Section 4 introduces the numerical examples and compares the solutions with the results of some approved approaches. In Section 5, we finally present our conclusions.

2. Fractional Expansion in Laplace Space

The fractional derivatives are non-local operators because they are defined using integrals. Therefore, the time-fractional derivative contains information about the function at earlier points, thus it possesses a memory effect. Such derivatives consider the history and non-local distributed effects, which are essential for better and more accurate descriptions and an understanding of complex and dynamic system behavior.
The fractional calculus in Caputo’s sense and the Laplace transform, which are necessary to create the LRPS solution for the fractional neutron diffusion equations with one delayed neutron group, are introduced in this part with some basic concepts and ideas.
Definition 1.
The time-fractional Caputo derivative of order  α  of  ψ x , t  sense is defined as
D t α ψ x , t = J t n α t n ψ x , t ,   n 1 < α n ,   x I .   t > 0 ,
where  D t α  denotes the time-fractional Caputo derivative operator of order  α ; n ; I  is an interval; and J t β  is the time-fractional Riemann–Liouville integral operator of order  β  that is defined as:
J t β ψ x , t = 1 Γ β 0 t t τ β 1 ψ x , τ d τ ,       β > 0   ,   t > τ 0 ψ x , t ,     β = 0 ,
In the following lemma, we introduce certain properties of the operator D t α that we will require for the rest of the work in this paper. The reader should see [26] for more features.
Lemma 1.
For   n 1 < α n , ε > 1 ,   t 0 , and   μ , we have:
( i )   D t α t ε = Γ ε + 1 Γ ε + 1 α t ε α , ( ii )   D t α μ = 0 , ( iii )   D t α J t α ψ x , t = ψ x , t , ( iv )   J t α D t α ψ x , t = ψ x , t j = 0 n 1 t j ψ x , 0 + t j j ! .
Definition 2.
Let ψ x , t be a continuous function on I × 0 , and of exponential order δ . Then, the Laplace transform of the function ψ x , t is denoted and defined as follows:
Ψ x , s = L ψ x , t 0 e s t ψ x , t d t ,   s > δ ,
whereas the inverse Laplace transform of the function Ψ x , s is defined as follows:
ψ x , t = L 1 Ψ x , s c i c + i e s t Ψ x , s d s , c = R e s > c 0 ,
where   c 0   lies in the right half plane of the absolute convergence of the Laplace integral.
In the next lemma, we introduce several crucial characteristics related to the Laplace transform and the fractional derivative in the Caputo sense.
Lemma 2.
Let  ψ x , t  be a continuous function on I × 0 , and of exponential orders δ , and  Ψ x , s = L ψ x , t . Then,
( i )   lim s s Ψ x , s = ψ x , 0 ,   x I ; ( ii )   L   J t α ψ x , t = s α 1 Ψ x , s ,   α > 0 ; ( iii )   L D t α ψ x , t = s α Ψ x , s k = 0 n 1 s α k 1 t k ψ x , 0 ,   n 1 < α < n ; , ( iv )   L D t m α ψ x , t = s m α Ψ x , s k = 0 m 1 s m k α 1 D t k α ψ x , 0 , 0 < α < 1 .
where  D t m α = D t α · D t α D t α   ( m -times) .
Now, as shown in the following result, we add a new fractional expansion that serves as the foundation for creating a LRPS solution for the partial differential equations (PDEs).
Theorem 1.
Let  ψ x , t  be continuous on  I × 0 ,  and of exponential order  δ . Suppose that the function Ψ x , s = L ψ x , t has the following fractional expansion:
Ψ x , s = n = 0 f n x s n α + 1 ,   0 < α 1 ,   x I ,   s > δ ,
then  f n x = D t n α ψ x , 0 .
The following theorem shows the conditions that guarantee the convergence of the series in the expansion (17).
Theorem 2.
Let  ψ x , t  be a continuous function on  I × 0 ,  and of exponential order  δ  and  Ψ x , s = L ψ x , t  can be represented as the fractional expansion in Theorem 1. If  s   L D t n + 1 α ψ x , t M x , on  I × δ , γ  where  0 < α 1 , then the remainder  R n x , s  of the fractional expansion as in Equation (17) satisfies the following inequality:
R n x , s M x ,   s n + 1 α + 1   ,   x I ,   δ < s γ
.

3. Exact Solution of the Fractional Neutron Diffusion Equations with One Delayed Neutron Group

The main objective of this section is to implement the LRPSM to construct a series solution to the fractional neutron diffusion equations with one delayed neutron group
1 v c D t α ψ x , t = D 2 ψ x , t + γ f a   ψ x , t + λ   c x , t , D t α c x , t = β γ f   ψ x , t λ   c x , t ,
with the initial conditions
ψ x , 0 = ψ 0 x ,     c x , 0 = β ν f λ ψ 0 x ,
Rewriting the coupled equations in (18) by replacing B = β γ f and = γ f a , we obtain
D t α ψ x , t = v c   D   ψ x x x , t + v c   ψ x , t + v c λ   c x , t , D t α c x , t = B ψ x , t λ   c x , t ,
with the initial conditions
ψ x , 0 = ψ 0 x ,     c x , 0 = B λ ψ 0 x ,  
Operating the Laplace transform on the coupled equations in (20), we obtain
L D t α ψ x , t = v c   D   L   ψ x x x , t + v c     L ψ x , t + v c λ   L c x , t , L D t α c x , t = B ψ x , t λ   L c x , t .
According to Lemma 2 part (iii) and using the initial conditions in Equation (21), we can write equations in (22) in the following form:
s α Ψ x , s s α 1 ψ 0 x = v c   D   Ψ x x x , s + v c     Ψ x , s + v c λ   C x , s , s α C x , s s α 1 B λ ψ 0 x = B Ψ x , s λ   C x , s ,
where Ψ x , s = L ψ x , t ,   C x , s = L c x , t .
Equation (23) can be written as
Ψ x , s 1 s ψ 0 x v c   D s α   Ψ x x x , s v c   s α   Ψ x , s v c λ s α   C x , s = 0 , C x , s B λ   s ψ 0 x B s α   Ψ x , s + λ s α   C x , s = 0 .
The coupled equations in Equation (24) represent a linear system of partial differential equations that contain derivatives relative to the independent variable x . According to the Laplace residual power series method, we assume that the series solution of the system in Equation (24) has the form
Ψ x , s = n = 0 f n x s n α + 1 , x I ,   s > δ 0 , C x , s = n = 0 g n x s n α + 1 , x I ,   s > δ 0 .
Using the fact in Lemma 2 part (i), of the kth truncated series of Ψ x , s , C x , s can be expressed as
Ψ k x , s = ψ 0 x s + n = 1 k f n x s n α + 1 , x I ,   s > δ 0 , C k x , s = B λ   s ψ 0 x + n = 1 k g n x s n α + 1 , x I ,   s > δ 0 .
In the next step, we define the Laplace residual functions of the coupled equations in (24) in order to find the unknown coefficients of the series in Equation (26):
L R e s Ψ x , s = Ψ x , s 1 s ψ 0 x v c   D s α   Ψ x x x , s v c   s α   Ψ x , s v c λ s α   C x , s , L R e s C x , s = C x , s B λ   s ψ 0 x B s α   Ψ x , s + λ s α   C x , s .  
and the kth Laplace-residual functions as
L R e s k Ψ x , s = Ψ k x , s 1 s ψ 0 x v c   D s α   Ψ k x x x , s v c   s α   Ψ k x , s v c λ s α   C k x , s , L R e s k C x , s = C k x , s B λ   s ψ 0 x B s α   Ψ k x , s + λ s α   C k x , s .  
Since L R e s Ψ x , s = 0 and L R e s C x , s = 0 , we have s k α + 1 L R e s Ψ x , s = 0 ,   s k α + 1 L R e s C x , s = 0 . Therefore,
lim s s k α + 1 L R e s k Ψ x , s = 0 ,   lim s s k α + 1 L R e s k C x , s = 0   for   k = 0 , 1 , 2 ,   .  
To find f 1 x and g 1 x in Equation (26), we substituted Ψ 1 x , s = ψ 0 x s + f 1 x s α + 1 and C 1 x , s = B λ   s ψ 0 x + g 1 x s α + 1   in the first Laplace-residual functions to obtain
L R e s 1 Ψ x , s = f 1 x s α + 1 v c   D s α   ψ 0 x s + f 1 x s α + 1 v c   s α   ψ 0 x s + f 1 x s α + 1 v c λ s α   B λ   s ψ 0 x + g 1 x s α + 1 , L R e s 1 C x , s = g 1 x s α + 1 B s α   ψ 0 x s + f 1 x s α + 1 + λ s α   B λ   s ψ 0 x + g 1 x s α + 1 .  
Next, by solving l i m s s α + 1 L R e s 1 Ψ x , s = 0   ,   l i m s s α + 1 L R e s 1 C x , s = 0 , one can obtain
f 1 x = v c   D   ψ 0 x + v c   ψ 0 x + v c B ψ 0 x ,  
g 1 x = 0 .  
To find f 2 x and g 2 x in Equation (26), we substituted Ψ 2 x , s = ψ 0 x s + f 1 x s α + 1 + f 2 x s 2 α + 1 and C 2 x , s = B λ   s ψ 0 x + g 2 x s 2 α + 1 in the second Laplace-residual functions to obtain
L R e s 2 Ψ x , s = f 1 x s α + 1 + f 2 x s 2 α + 1 v c   D s α   ψ 0 x s + f 1 x s α + 1 + f 2 x s 2 α + 1 v c   s α   ψ 0 x s + f 1 x s α + 1 + f 2 x s 2 α + 1 v c λ s α   B λ   s ψ 0 x + g 2 x s 2 α + 1 , L R e s 2 C x , s = g 2 x s 2 α + 1 B s α   ψ 0 x s + f 1 x s α + 1 + f 2 x s 2 α + 1 + λ s α   B λ   s ψ 0 x + g 2 x s 2 α + 1 .  
Next, by solving l i m s s 2 α + 1 L R e s 2 Ψ x , s = 0 ,   l i m s s 2 α + 1 L R e s 2 C x , s = 0 , one can obtain
f 2 x = v c   D   f 1 x + v c     f 1 x ,  
g 2 x = B f 1 x .
In the third Laplace, the residual functions to obtain
L R e s 3 Ψ x , s = f 1 x s α + 1 + f 2 x s 2 α + 1 + f 3 x s 3 α + 1 v c   D s α   ψ 0 x s + f 1 x s α + 1 + f 2 x s 2 α + 1 + f 3 x s 3 α + 1 v c   s α   ψ 0 x s + f 1 x s α + 1 + f 2 x s 2 α + 1 + f 3 x s 3 α + 1 v c λ s α   B λ   s ψ 0 x + g 2 x s 2 α + 1 + g 3 x s 3 α + 1 , L R e s 3 C x , s = g 2 x s 2 α + 1 + g 3 x s 3 α + 1 B s α   ψ 0 x s + f 1 x s α + 1 + f 2 x s 2 α + 1 + f 3 x s 3 α + 1 + λ s α   B λ   s ψ 0 x + g 2 x s 2 α + 1 + g 3 x s 3 α + 1 .  
By solving l i m s s 3 α + 1 L R e s 3 Ψ x , s = 0 ,   l i m s s 3 α + 1 L R e s 3 C x , s = 0 , one obtains
f 3 x = v c   D   f 2 x + v c   f 2 x + v c λ g 2 x ,  
g 3 x = B f 2 x λ g 2 x .  
If we continue in the same manner, we can substitute the kth truncated series Ψ k x , s ,   C k x , s into the kth Laplace-residual function L R e s k Ψ x , s ,   L R e s k C x , s . Then, we can multiply the resulting equations by s k α + 1 and taking the limit as s , f k + 1 x , g k + 1 x for k 2 can be obtained by the following recurrence relation:
f n x = v c   D   f n 1 x + v c   f n 1 x + v c λ g n 1 x ,  
g n x = B f n 1 x λ g n 1 x ,   n = 1 , 2 , 3 ,
where g 0 x = B λ   ψ 0 x .
According to what was introduced, the series solution of system (24) is
Ψ x , s = ψ 0 x s + n = 1 f n x   s n α + 1 ,   x I ,   s > δ 0 . C x , s = B λ   s   ψ 0 x + n = 1 g n x   s n α + 1 ,   x I ,   s > δ 0 .
Therefore, the series solution of system (18) can be obtained by transforming the solution in (41) into the original space by using the inverse Laplace transform. Therefore, the Laplace residual power series solution of system (18) is given by:
ψ x , t = ψ 0 x + n = 1 f n x   t n α Γ n α + 1 , t 0 , x I . c x , t = B λ     ψ 0 x + k = 3 g n x   t n α Γ n α + 1 , t 0 , x I ,
where
f n x = v c   D   f n 1 x + v c   f n 1 x + v c λ g n 1 x ,   g n x = B f n 1 x λ g n 1 x , g 0 x = B λ   ψ 0 x , n = 1 , 2 , 3 ,   .

4. Numerical Results and Discussion

To verify the driven theory, we solved the time-dependent fractional neutron diffusion equations with delayed neutrons using the Laplace residual power series method, and the results were tested with the numerical values of the following nuclear reactor cross-section data [10]
v c = 220 , 000   cm / s , B = 0.000735   cm 1 , D = 0.356 ,   λ = 0.08   s 1 ,   Σ = 0.005   cm 2 .
It must be noted that these cross-sections are defined in Section 1.
The values of f n x   and   g n x were obtained using MATHEMATICA in 11th Gen Intel(R) Core (TM) i7-1165G7 2.80 GHz with a 16 GB RAM computer. Once f n x   and   g n x are calculated, the data to construct the tables and figures are achieved.
The numerical results of neutron flux at different values of α are shown in Table 1 for ψ 0 x = 1 .
The flux was calculated at different times, and the anomalous diffusion, α < 1 (sub diffusion) was considered. Additionally, it is clear from Table 1 that the neutron flux increased with time, as shown in Figure 1.
For more clarification, the neutron fluxes ψ x , t are plotted in Figure 1 with the initial condition ψ 0 x = x 2 .
Now, the different initial value problem ψ x , 0 = x 2 was used in the study of the following case, where the results are clarified in Figure 2.
The plotted profile of the neutron fluxes for different values of α ( α = 0.1 ,   0.2 , 0.3, 0.4 and α = 0.5 ) in both cases { ψ x , 0 = 1   and ψ x , 0 = x 2   } , as presented in Table 1 and Figure 1 and Figure 2, shows that the results were in agreement with the Adomian decomposition method [10].
As observed from the tables and figures, the Laplace residual power series method used to solve the time-dependent fractional neutron diffusion equations with delayed neutrons provided us with more details about the neutron flux in non-Gaussian diffusion for different selected times.
The importance of this work was manifested by uncovering the behavior of the flux in the presence of the delayed neutrons that affected them (this phenomenon is known as anomalous diffusion). This flux can be expressed mathematically by the solution of the time-dependent fractional neutron diffusion equations, as shown in Equation (42).
It must be noted that the time-dependent flux increased with time in both the studied numerical cases [9,10,12].

5. Conclusions

The Laplace residual power series method was presented efficiently by studying the coupled time-dependent fractional neutron diffusion equations with delayed neutrons to create an approximate fractional series solution, and the validity of the theoretical presentation of the method was confirmed by two different nuclear physics numerical case studies. The obtained fluxes for both cases were calculated and the anomalous diffusion, α < 1 (sub diffusion) was considered where the neutron fluxes increased with time; moreover, the results were successfully tested when they were compared with the Adomian decomposition method. The newly introduced concept of anomalous diffusion, which describes non-Gaussian diffusion, can be explained mathematically by fractional calculus using the Laplace residual power series method. It must be noted that the success of the combination of the residual power series and Laplace transformation in the proposed method has the potential for use in studying different nuclear reactor theories and other scientific phenomena.

Author Contributions

Conceptualization, M.S.; Methodology, A.B.; Software, I.G.; Investigation, M.S., I.G. and A.B.; Writing—original draft, A.B, M.S. and I.G.; Writing—review & editing, M.S. and I.G.; Visualization, I.G.; Supervision, A.B.; Project administration, M.S.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stacey, W.M. Nuclear Reactor Physics; John Wiley & Sons: Boston, MA, USA, 2001. [Google Scholar]
  2. Duderstadt, J.J.; Hamilton, L.J. Nuclear Reactor Analysis; John Wiley & Sons: Boston, MA, USA, 1976. [Google Scholar]
  3. Lamarsh, J.R. Introduction to Nuclear Engineering, 2nd ed.; Addison-Wesley: Boston, MA, USA, 1983. [Google Scholar]
  4. Khasawneh, K.; Dababneh, S.; Odibat, Z. A solution of the neutron diffusion equation in hemispherical symmetry using the homotopy perturbation method. Ann. Nucl. Energy 2009, 36, 1711–1717. [Google Scholar] [CrossRef]
  5. Dababneh, S.; Khasawneh, K.; Odibat, Z. An alternative solution of the neutron diffusion equation in cylindrical symmetry. Ann. Nucl. Energy 2010, 38, 1140–1143. [Google Scholar] [CrossRef]
  6. Shqair, M. Solution of different geometries reflected reactors neutron diffusion equation using the homotopy perturbation method. Results Phys. 2019, 12, 61–66. [Google Scholar] [CrossRef]
  7. Shqair, M. Developing a new approaching technique of homotopy perturbation method to solve two-group reflected cylindrical reactor. Results Phys. 2019, 12, 1880–1887. [Google Scholar] [CrossRef]
  8. Shqair, M.; Farrag, E.; Al-Smadi, M. Solving Multi-Group Reflected Spherical Reactor System of Equations Using the Homotopy Perturbation Method. Mathematics 2022, 10, 1784. [Google Scholar] [CrossRef]
  9. Vyawahare, V.; Nataraj, P. Fractional-Order Modeling of Nuclear Reactor: From Subdiffusive Neutron Transport to Control-Oriented Models; Springer Nature: Singapore, 2018. [Google Scholar]
  10. Ray, S. Fractional Calculus with Applications for Nuclear Reactor Dynamics; CRC Press: Broken, NW, USA, 2016. [Google Scholar]
  11. Aboanber, A.; Nahla, A.; Aljawazneh, S. Fractional two energy groups matrix representation for nuclear reactor dynamics with an external source. Ann. Nucl. Energy 2021, 153, 108062. [Google Scholar] [CrossRef]
  12. Sardar, T.; Ray, S.; Bera, R. The solution of coupled fractional neutron diffusion equations with delayed neutrons. Int. J. Nuclear Energy Sci. Technol. 2010, 5, 105–133. [Google Scholar] [CrossRef]
  13. Khaled, S. Exact solution of the one-dimensional neutron diffusion kinetic equation with one delayed precursor concentration in Cartesian geometry. AIMS Math. 2022, 7, 12364–12373. [Google Scholar] [CrossRef]
  14. Burqan, A.; El-Ajou, A.; Saadeh, R.; Al-Smadi, M. A new efficient technique using Laplace transforms and smooth expansions to construct a series solution to the time-fractional Navier-Stokes equations. Alex. Eng. J. 2022, 61, 1069–1077. [Google Scholar] [CrossRef]
  15. Sarhan, A.; Burqan, A.; Saadeh, R.; Al-Zhour, Z. Analytical Solutions of the Nonlinear Time-Fractional Coupled Boussinesq-Burger Equations Using Laplace Residual Power Series Technique. Fractal Fract. 2022, 6, 631. [Google Scholar] [CrossRef]
  16. Kazem, S. Exact solution of some linear fractional differential equations by Laplace transform. Int. J. Nonlinear Sci. 2013, 16, 3–11. [Google Scholar]
  17. Yavuz, M.; Özdemir, N. Numerical inverse Laplace homotopy technique for fractional heat equations. Therm. Sci. 2018, 22, 185–194. [Google Scholar] [CrossRef]
  18. Wang, Q. Numerical solutions for fractional KdV–Burgers equation by Adomian decomposition method. Appl. Math. Comput. 2006, 182, 1048–1055. [Google Scholar] [CrossRef]
  19. Cao, Y.; Nikan, O.; Avazzadeh, Z. A localized meshless technique for solving 2D nonlinear integro-differential equation with multi-term kernels. Appl. Numer. Math. 2023, 183, 140–156. [Google Scholar] [CrossRef]
  20. Akram, T.; Abbas, M.; Ali, A. A numerical study on time fractional Fisher equation using an extended cubic B-spline approximation. J. Math. Comput. Sci. 2021, 22, 85–96. [Google Scholar] [CrossRef]
  21. Oderinu, R.A.; Owolabi, J.A.; Taiwo, M. Approximate solutions of linear time-fractional differential equations. J. Math. Comput. Sci. 2023, 29, 60–72. [Google Scholar] [CrossRef]
  22. Ali, A.; Abbas, M.; Akram, T. New group iterative schemes for solving the two-dimensional anomalous fractional sub-diffusion equation. J. Math. Comp. Sci. 2021, 22, 119–127. [Google Scholar] [CrossRef]
  23. Liao, H.; Zhang, D.; Zhao, Y. Sharp error estimate of a nonuniform L1 formula for time-fractional reaction subdiffusion equations. SIAM J. Numer. Anal. 2018, 56, 1122–1133. [Google Scholar] [CrossRef]
  24. Lyu, P.; Vong, S. A symmetric fractional-order reduction method for direct nonuniform approximations of semilinear diffusion-wave equations. J. Sci. Comput. 2022, 93, 34. [Google Scholar] [CrossRef]
  25. Eriqat, T.; El-Ajou, A.; Moa’ath, N.O.; Al-Zhour, Z.; Momani, S. A new attractive analytic approach for solutions of linear and nonlinear neutral fractional pantograph equations. Chaos Solitons Fractals 2020, 138, 109957. [Google Scholar] [CrossRef]
  26. Oldham, K.; Spanier, J. The Fractional Calculus Theory and Applications of Differentiation and Integration to Arbitrary Order, 1st ed.; Elsevier: Amsterdam, The Netherlands, 1974. [Google Scholar]
Figure 1. Two-dimensional graphs of neutron flux ψ x ,   t for α = 0.1, 0.2, 0.3, 0.4, and 0.5 (colored online).
Figure 1. Two-dimensional graphs of neutron flux ψ x ,   t for α = 0.1, 0.2, 0.3, 0.4, and 0.5 (colored online).
Fractalfract 07 00219 g001
Figure 2. Three-dimensional graphs of neutron flux ψ x ,   t with ψ 0 x = x 2 for α = 0.1, 0.2, 0.3, 0.4 and 0.5 (colored online). Two-dimensional graphs of neutron flux ψ x ,   t for α = 0.1, 0.2, 0.3, 0.4, and 0.5 (colored online).
Figure 2. Three-dimensional graphs of neutron flux ψ x ,   t with ψ 0 x = x 2 for α = 0.1, 0.2, 0.3, 0.4 and 0.5 (colored online). Two-dimensional graphs of neutron flux ψ x ,   t for α = 0.1, 0.2, 0.3, 0.4, and 0.5 (colored online).
Fractalfract 07 00219 g002
Table 1. Numerical results of neutron flux at different values of α ,   ψ 0 x = 1 .
Table 1. Numerical results of neutron flux at different values of α ,   ψ 0 x = 1 .
t α = 0.1 α = 0.2 α = 0.3 α = 0.4 α = 0.5
0.000101.08E+086,841,627405,001.422,937.431302.468
0.000391.62E+0815,459,5761,371,299115,295.39415.568
0.000681.91E+0821,570,6042,258,076223,540.721,346.35
0.000972.13E+0826,687,1593,105,913341,437.736,086.62
0.001262.3E+0831,216,2113,927,982466,526.553,165.07
0.001552.45E+0835,342,3444,730,899597,429.772,290.41
0.001842.58E+0839,168,5945,518,637733,25393,259.25
0.002132.69E+0842,759,6876,293,845873,363.2115,919.1
0.002422.8E+0846,159,5517,058,4061,017,286140,150
0.002712.89E+0849,399,6837,813,7281,164,650165,854.5
0.003002.98E+0852,503,6118,560,9021,315,159192,951.2
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

Shqair, M.; Ghabar, I.; Burqan, A. Using Laplace Residual Power Series Method in Solving Coupled Fractional Neutron Diffusion Equations with Delayed Neutrons System. Fractal Fract. 2023, 7, 219. https://doi.org/10.3390/fractalfract7030219

AMA Style

Shqair M, Ghabar I, Burqan A. Using Laplace Residual Power Series Method in Solving Coupled Fractional Neutron Diffusion Equations with Delayed Neutrons System. Fractal and Fractional. 2023; 7(3):219. https://doi.org/10.3390/fractalfract7030219

Chicago/Turabian Style

Shqair, Mohammed, Ibrahim Ghabar, and Aliaa Burqan. 2023. "Using Laplace Residual Power Series Method in Solving Coupled Fractional Neutron Diffusion Equations with Delayed Neutrons System" Fractal and Fractional 7, no. 3: 219. https://doi.org/10.3390/fractalfract7030219

APA Style

Shqair, M., Ghabar, I., & Burqan, A. (2023). Using Laplace Residual Power Series Method in Solving Coupled Fractional Neutron Diffusion Equations with Delayed Neutrons System. Fractal and Fractional, 7(3), 219. https://doi.org/10.3390/fractalfract7030219

Article Metrics

Back to TopTop