ORIGINAL RESEARCH article

Front. Phys., 11 September 2023

Sec. Mathematical Physics

Volume 11 - 2023 | https://doi.org/10.3389/fphy.2023.1229142

A solution for the neutron diffusion equation in the spherical and hemispherical reactors using the residual power series

  • 1. Faculty of Science, Al Balqa Applied University, Al-Salt, Jordan

  • 2. Faculty of Science, Zarqa University, Zarqa, Jordan

Abstract

A novel analytical solution to the neutron diffusion equation is proposed in this study using the residual power series approach for both spherical and hemispherical fissile material reactors. Various boundary conditions are investigated, including zero flux on the boundary, zero flux on the extrapolated boundary distance, and the radiation boundary condition (RBC). The study also shows how two hemispheres with opposing flat faces interact. We give numerical results for the same energy neutrons diffused in pure . By qualitative comparison with the homotopy perturbation method and Bessel function-based solutions, the residual power series method (RPSM) presents accurate series solutions that converge to the exact solutions, as shown in this study. Moreover, numerical results were shown to be improved by the computer implementation of the analytic solutions.

1 Introduction

Modern physics is taken into consideration after classical physics failed to explain various phenomena, and many branches of modern physics developed, including nuclear physics. To perform nuclear fission in a nuclear reactor, it is essential to have a neutron in the thermal state to break through the nucleus, as each fission produces more than one neutron, leading to a new reaction where a chain reaction is created [1].

The continuity equation, known as the time-dependent transfer equation, expresses the distribution of neutrons in the reactor, which is given aswhere is the neutron velocity, is the neutron flux, is the neutron source term, is an absorption macroscopic cross section, and is the position vector at the time .

The continuity equation involves both neutron current and neutron flux , where the relationship between current and the flux can be shown using Fick’s law, which is given in the following mathematical equation:where is the diffusion coefficient.

By applying Fick’s law, the time-dependent transport equation can be simplified to the following time-dependent neutron diffusion equation:whose solution is considered very important in studying the behavior of neutrons in reactors [2, 3].

The stability of the nuclear reactor has been studied in detail in the past [4]. The condition of criticality in a nuclear reactor is essential to ensure its stability, as the number of neutrons in the reactor must be constant at any time. The flow in the critical state must be time-independent, which can be mathematically expressed using a time-independent diffusion equation as follows:

Various techniques, analytical and numerical, have been employed to solve the neutron diffusion equation throughout the past few decades, such as the homotopy perturbation method, the Laplace transform process, and the Galerkin method with B-splines [521]. To the best of our knowledge, the residual power series method (RPSM) [2129], an analytical technique, has never been used to handle issues relating to the neutron diffusion equation.

RPSM is a relatively new technique that provides exact and approximate analytical solutions in the series form to many differential and integrodifferential equations of integer and fractional orders [2129]. The idea of the RPSM is based on using the residual function and the derivative operator concepts in determining the coefficients of the solution, which is assumed to be in a power series form. In RPSM, it is usual to use a classical derivative operator to determine the coefficients of a power series solution if the differential equation is of integer order. However, a fractional derivative operator is used to determine the coefficients if the differential equation is of fractional order.

Later, we will rewrite Eq. 1 using spherical coordinates in spherical and hemispherical reactors, which gives rise to single equations. This work aims to adapt and modify RPSM to solve Eq. 1 after being formulated in spherical coordinates. This modification is a regeneration of RPSM and is presented for the first time in this work. This modification relies on using a fractional derivative operator, such as the Caputo fractional derivative (CFD) [30, 31] or the conformable fractional derivative [21, 32], in differential equations of the integer order. Since the equations to be solved are singular, we impose the solution in a Frobenius series form. Therefore, the new modified method can be considered an alternative method to the Frobenius method in determining the power series solution coefficients.

In fact, several types of fractional derivative operators can be used in our hypothetical modification of the RPSM, such as the CFD operator, the Riemann–Liouville fractional derivative operator [33], and the conformable fractional derivative operator. The CFD is employed in the modified method without bias or preference to determine the Frobenius series coefficients. The CFD operator is defined using the following formula:where , is the CFD operator and is the gamma function.

Moreover, the numerical results in this study are contrasted using pure for the first time in studying spherical and hemispherical neutron diffusion equations, whereas the previous works in history studied pure .

The following section presents the methodology of the RPSM and its application to spheres and hemispheres. The radiation boundary condition is also covered in depth to set the stage for the computations. Section 3 presents numerical results for the flux distribution in the studied nuclear reactors and critical radius calculations. The associated technical problems encountered throughout the computations are presented, along with the corresponding solutions.

2 Theory

2.1 The basic idea of the RPSM

RPSM is a technique to find the power series solution coefficients without having a recurrence relation. Indeed, the coefficients of the power series are calculated by the concept of residual error through a sequence of algebraic equations, and in the end, a truncated series solution (approximate solution) is obtained. The major merit of RPSM is that it can be implemented to the problem directly without linearization, perturbation, or discretization and without any transformation by selecting appropriate initial conditions.

Most of the problems resolved using RPSM were to provide series solutions about regular points and did not address situations where there are singular points. In this work, we introduce a simple modification to the construction of the RPSM to be able to solve issues containing these types of points.

The basic definition in addition to the basic theories of RPSM and its applicability to different types of differential equations is given in [

21

29

]. For the convenience of the reader, we will introduce a review of the RPSM by clarifying the following algorithm on an ordinary differential equation:

  • (1) Express the th-order ordinary differential equation in the following general form:

subject to the initial conditions

where

is an unknown function of x on an open interval

containing

and

is a nonlinear differential operator of order

, and its coefficients are polynomials.

  • (2) If is a regular point for Eq. 3, assume the solution has a power series expansion:

whereas if

is a singular point for Eq.

3

, assume the solution has the following expansion (this is the new case that the RPSM did not address):

  • (3) Define the residual functions as follows:

Indeed,

, for every

and so

at every interior point of

.

  • (4) Substitute the expansions (5) and (6) of into Eq. 7 to obtain, respectively, the following expansions:

where

is the gamma function.

  • (5) In the singular point case, we need to find the value of . Therefore, we find the CFD of order for the residual function in Eq. 9 as follows:

where

is the CFD of order

, which has the following two properties: 1)

is a constant 2)

.

  • (6) Since the CFD of a constant is 0, the value of is the solution of the following indicial equation:

Since

and

for all

in the interior of

, the following equations are true:

  • (7) If is a regular point for Eq. 3, then , and for can be determined by solving the following algebraic equations, sequentially:

whereas if

is a singular point for Eq.

3

, then

is arbitrary and can be determined later by the given initial conditions, and

for

can be specified by solving the following algebraic equations, sequentially:

  • (8) Replace the values obtained for , with their place in the expanded series of to obtain the th approximate solution to Eq. 3.

  • (9) If there is a pattern in the coefficients of the series in terms of some well-known elementary functions, then we have the exact solution .

2.2 Bare sphere

First, the application of RPSM in spherical reactors is studied. We consider first the rather simple example of a bare spherical reactor of radius . The spherical symmetry of the system implies that the flux is a function of only. The time-independent diffusion Eq. 1 can thus be written aswhereis the buckling of the reactor. Substituting the Laplacian in spherical coordinates and using , we obtain

To solve Eq. 18 using RPSM, we assume that the solution has a fractional power series representation at as follows:

Define the residual function of Eq. 18 as

Substitute Eq. 19 into the residual function (20), we obtain the following series:

According to Eq. 10, the CFD of the residual function in Eq. 21 will be as

The solution of the indicial equations gives . Thus, Eq. 22 becomes

Thus, the th derivative of the residual function has the following general form:

According to step (9) in Subsection 2.1, the solution of the algebraic equationscan be summarized in Table 1.

TABLE 1

1
2
3
4

Some values of the coefficients of the expansion (19).

If we continue in the same manner, then we conclude that the general form of the coefficients in series (19) will be in the following form:

Therefore, the solution for Eq. 18 has the following exact solution:

Consequently, the final solution for Eq. 16, in a closed form, is given by

This flux must satisfy the relevant boundary conditions.

2.3 Hemispherical reactor

In this subsection, we will consider the hemispherical reactor. This case is more complicated because the flux is a function in two variables and . Therefore, the time-independent diffusion Eq. 1 can be expressed, by writing the Laplacian in the spherical coordinates and using , as follows:where the buckling is again given by Eq. 17.

2.3.1 Solution by the separation of variables and residual power series method

Applying the separation of variables, we obtain

We obtain the eigenvalue equations as follows:

To solve the radial part using RPSM, we first, for simplicity, consider and rewrite Eq. 31 asand we assume the solution has a power series representation at ,

The residual function of Eq. 33 has the following form:

If we substitute Eq. 34 into the residual function (35), then we obtain the following series function:and so, the CFD of the residual function of order in Eq. 36 has the following expansion:

Substituting in Eq. 37 gives the indicial equation of series (34), whose solution is . Hence, the residual function in Eq. 37 can be rewritten as

Following the same manner as in Section 2.2, if we apply the sequential operator , on both sides of Eq. 38, we obtain

Consequently, the first few coefficients of the series in Eq. 34 are given in Table 2.

TABLE 2

1
2
3
4

Some values of the coefficients of the expansion (34).

It is clear that there is a pattern in the coefficients of our series. So, the th coefficient, , has a general form as follows:which can be reformulated as follows:

Therefore, the eigensolutions to the Sturm–Liouville equation in Eq. 31 appear in the following series:

Now, to solve the angular part, which is a Legendre equation, we construct the residual function of Eq. 32 asand since is a regular point of Eq. 32, we assume the solution of the following expansion:

Substitute Eq. 44 into Eq. 43, and we obtain

Applying the differential operator , on both sides of Eq. 45 gives

Thus, the first few coefficients of the series in Eq. 44 are given in Table 3.

TABLE 3

0
1
2
3

Some values of the coefficients of expansion (44).

Note that there are two linearly independent solutions one of them is an infinite series and the other is a polynomial of degree m. If m is an even number, then the solution that contains terms of even exponents will be the polynomial whereas the solution that contains terms of odd exponents will be infinite series and vice versa if m is an odd number. Anyway, we can keep track of the pattern in the coefficients and write the solutions as follows:

Legendre was able to combine both solutions in a single format, regardless of whether the number was odd or even as

By using Rodrigues’s formula for the Legendre polynomials, the solution can be expressed as follows:

Therefore, the general solution for Eq. 29 can be expressed as

2.4 The radiation boundary condition

Although the choice of suitable boundary conditions in solving the neutron diffusion equation is at the core of the nuclear reactor physicist mission, here, we will apply three boundary conditions, namely, zero flux boundary condition (ZFBC), extrapolated boundary condition (EBC), and radiation boundary condition (RBC). These boundary conditions will be applied on both spherical and hemispherical reactors.

ZFBC cures about the point at the surface of the reactor that the flux is 0, and for more reality, the flux will vanish at a point less than that point which is called extrapolated distance . To be more accurate, especially for the hemisphere, the mixed boundary condition gives the closest results when comparing to the benchmark; this boundary condition is called RBC

The use of the application of RBC on the solution of the hemispherical reactor using RPSM based on the separation of variable techniques can be performed as follows:where varies over the surface, and this RBC is reported in previous works [34, 35].

The RBC on the curved surface where will be

However, RBC on the flat surface will bewhere where the values of g change from 0 to 1, and the reactor geometry changes from spherical (0) to hemispherical ( 1) geometry as mentioned in [34, 35].

Applying RBC on the flat surface in Eq. 50 yields

Using the following recurrence relationships of :and some Legendre polynomials properties, we obtain the following relationship:

We can use the previous result to write the odd term amplitude in terms of the even term amplitude aswhere

as reported in [ 36].

Similarly, applying the RBC on the curved surface, we obtainor it can be written aswhere

Following [35], we obtainwhereand

The alternating method of iswhere , which is used in these calculations.

3 Numerical results

The numerical data of neutrons diffusing in pure are taken from [36], which are shown in Table 4.

TABLE 4

ʋfs
b b b

Pure cross sections.

The important diffusion constant and buckling are

Depending on these data, to calculate the flux, applying different boundary conditions, to find the critical radii, and to study the flux distribution for both spherical and hemispherical reactors,, this work has been compared with benchmark calculations [34].

3.1 Critical radius calculations

The critical radius is the dimension of the reactor determined as the flux vanished. The determination of the critical radius depends on the boundary conditions, and the simplest way of calculating it is by using ZFBC, where the flux vanishes. In real reactors, an extrapolated distance must be considered ( +2D) at which the value of the flux is 0; this boundary condition is the EBC. The use of ZFBC and EBC is suitable for sphere, cylinder, and slab reactors where there is high symmetry, but for more complicated reactors like hemispherical and cone reactors, we need to use the RBC.

The first case we will study is the spherical reactor. Here, the critical radius using ZFBC and EBC is considered, and the results are compared with RBC. A computer code based on Eq. 28 has been written to find the first root of Eq. 28. The value of at which the flux is 0 is cm. Since the flux is assumed to converge at using EBC, then . These results are tabulated in Table 5.

TABLE 5

BCSphere%ErrorHemisphere%Error
ZFBC12.5674
EBC10.5249
RBC6.5021-9.0956-

Critical radius (in cm) for pure for a sphere and hemisphere using different boundary conditions: ZFBC, EBC, and RBC.

After taking the bare sphere, the hemispherical reactor is considered, and the neutron flux using ZFBC will vanish on the flat surface (µ = 0). Using the Legendre polynomial properties, the even amplitudes must be equal to 0, and on the curved surface, we need to find the first zero of the solution of the radial part. The unique available case is for . So Eq. 50 will reduce to the following form:

The critical radii in pure that are calculated using the code depending on Eq. 59 are tabulated in Table 5.

Unlike the spherical reactor, applying RBC is essential for the hemispherical reactor.

The calculation of the critical dimension using RBC is illustrated in the following equation, and the infinite summation in Eq. 62 will be cut off to the value N, where the dimension of the matrix is finite, and rewritten aswhere

This simplification is an important numerical calculation to make cut-off limit value N. Then, the flux in Eq. 50 will be simplified to

A third computer code has been used to calculate the critical radius using RBC. The code calculates the determinant of . The value of a affects the determinant through the parameters (Eq. 61). As increases, the value of the determinant will change its sign. On the other hand, the calculated depends on the value of the interaction parameter g through (Eq. 58). Hence, the results for the sphere and the infinitely separated hemispheres are obtained by using and , respectively.

The determination of the critical dimension depends on the cut-off value N, where the number of terms will be , and the value of the critical dimension with the number of terms is given in Table 6.

TABLE 6

Number of termsCritical radius
4
69.0955147
89.0955608
109.0955676
129.0955673
189.0955655
229.0955651
329.0955652
389.0955657
429.0955664

Convergence of the critical radius (in cm) for the hemisphere using RBC.

The effect of the number of terms is clear; the value of the critical radius is converged by increasing the number of terms. Now, the effect of changing the interaction parameter on determining the critical radius is clarified in Table 7, where the cut-off value and the number of terms (22) are considered.

TABLE 7

Critical radius
06.50206
0.16.96796
0.27.3793148
0.37.7345081
0.48.0369619
0.58.2927704
0.68.5088470
0.78.6918102
0.88.8474778
0.98.9807312
19.0955651

Variation in the critical radius with the interaction parameter . (22 terms).

The reactor is a sphere when the interaction parameter , whereas it is a hemisphere when the interaction parameter , where the values of the critical radii are shown in Table 4 and Table 5, respectively.

3.2 Flux distribution

Now, the flux destruction inside the hemispherical reactor will be considered. It is straightforward that the symmetry of the spherical reactor makes the flux distribution a simple case study in [35], so the flux in Eq. 69 will be written as

For computational calculation purpose, this equation can be expressed as

The amplitudes in Eq. 62 are needed to be written in the formwhere has been normalized to unity. So A0, and the flux values are normalized to the averaged flux.

The value of the average flux is

Based on these data, the averaged flux calculated using the code was determined to be 0.756965. The variation in the flux at with changes from zero, and the critical radius is given in Table 8 for different values of N (different number of terms).

TABLE 8

4 terms12 terms22 terms46 terms62 terms86 terms
00.964770.964770.964770.964770.964770.96477
0.11.24571.245591.245591.24561.24561.2456
0.21.469421.469051.469071.469081.469081.46908
0.31.621841.621281.621311.621331.621331.62133
0.41.693951.693461.693511.693531.693531.69353
0.51.682381.682491.682551.682581.682581.68258
0.61.589681.591131.59121.591231.591241.59124
0.71.4241.427731.427811.427841.427841.42785
0.81.198411.205481.205561.205591.205591.20559
0.90.9298330.9412990.9413650.9413940.9413970.941399
10.637630.6544690.6545140.6545390.6545420.654544

Flux variation across the critical hemisphere at

The convergence of the flux values is clearly given in Table 8.

As shown in Figure 1, the fast convergence of the flux on the curved surface is obtained with the increasing number of terms (for and 30).

FIGURE 1

It is important to clarify, from Figure 1, that the convergence is stable and acceptable, and when the number of terms is 30, the angular flux, which expresses the behavior of neutrons when they are diffusing in pure , is expected which is the same as in using the homotopy perturbation method in solving the neutron diffusion equation where neutrons diffuse in the hemispherical reactor [35].

After studying the effect of the number of terms on the flux, it is taken into consideration at a wide range of angle values in Figure 2.

FIGURE 2

As shown in Figure 2, it is noted that the radial flux does not decrease when increases and the maximum value of the flux shifted with the increase in the angle. One more time, the neutron distribution is in the same manner as in previous work [35].

To clarify the flux behavior in more detail, the flux values for different angles and radii are tabulated in Table 9 and depicted in Figure 3.

TABLE 9

0.00.964770.964770.964770.964770.964770.964770.964770.964770.964770.96477
0.11.362131.356151.338341.309141.269251.219671.161651.096731.026660.953433
0.21.680611.669841.637741.58491.512331.421461.314181.192841.060310.919961
0.31.90311.888861.846351.776161.679311.557271.412011.246081.062710.865935
0.42.019262.002901.953991.873041.760911.618861.448541.252171.032660.793886
0.52.026172.008991.95761.872371.753961.603261.421441.209990.9709130.707143
0.61.928371.911581.861291.777791.661511.512981.332811.121630.8801470.60964
0.71.737341.721971.675921.599341.492521.35571.189020.9922910.7648280.505694
0.81.470491.457391.418041.352541.261051.143651.000170.8298780.6309080.399781
0.91.149581.139471.108581.057080.9850750.892580.7793180.6443780.4853810.296389
1.00.7954680.7928330.7709630.7351940.6850560.6204510.5414520.4471780.3355930.200533

Flux variation across the critical hemisphere at different values of . (22 terms).

FIGURE 3

Furthermore, a three-dimensional plot of the flux distribution across the core of a bare hemisphere as a function of and is shown in Figure 3.

The flux in Table 6 is represented graphically in Figure 3. The data in Figure 1 and Figure 2 are merged in the same figure, where the full neutron behavior is represented; however, this three-dimensional flux representation is not figured out previously for the hemispherical reactor.

Although the neutrons diffusing in pure in this paper are studied using RPSM for the first time and the neutrons diffusing in were previously studied using the homotopy perturbation method [35, 36], it is clear that their behavior is the same.

4 Conclusion

In this article, we introduced the application of RPSM in solving the neutron diffusion equation in hemispherical symmetry. The method is efficient and can present the solution in the form of a rapid series that converge to the exact solution. Two important cases are discussed: the spherical reactor and hemispherical reactor. The results are presented in tables and illustrative figures using codes. To show the validity of the presented method, we made qualitative comparisons between the obtained results using RPSM and those computed using the homotopy perturbation method. The used method and its computational implementation provide superior results when applied to more complex and reflected reactor geometries for one-group problems. It is recommended, however, for future work to further investigate the methodology for two- and multi-group problems.

In future work, we attend to solve more partial differential equations with physical applications, using the proposed method that showed its accuracy in handling similar problems, and further investigation of the convergence-produced series should be the goal of future specialized research [3741].

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1.

    KraneKS. Introductory nuclear physics. John Wiley & Sons (1991).

  • 2.

    DuderstadtJJHamiltonLJ. Nuclear reactor analysis. Boston, MA, USA: John Wiley & Sons (1976).

  • 3.

    LamarshJR. Introduction to nuclear engineering. 2nd ed.Boston, MA, USA: Addison-Wesley (1983).

  • 4.

    KochenovAS. The stability of a nuclear power plant. J Nucl Energ A. Reactor Sci (1960) 12(4):194200. 10.1016/0368-3265(60)90103-1

  • 5.

    DababnehSKhasawnehKOdibatZ. An alternative solution of the neutron diffusion equation in cylindrical symmetry. Ann Nucl Energ (2011) 38(5):11403. 10.1016/j.anucene.2010.12.011

  • 6.

    NairatMShqairMAlhalholyT. Cylindrically symmetric fractional Helmholtz equation. Appl Math E-notes (2019) 19:70817.

  • 7.

    ShqairM. Developing a new approaching technique of homotopy perturbation method to solve two-group reflected cylindrical reactor. Results Phys (2019) 12:18807. 10.1016/j.rinp.2019.01.063

  • 8.

    ShqairMFarragEAAl-SmadiM. Solving multi-group reflected spherical reactor system of equations using the homotopy perturbation method. Mathematics (2022) 10(10):1784. 10.3390/math10101784

  • 9.

    Ala’yedOSaadehRQazzaA. Numerical solution for the system of Lane-Emden type equations using cubic B-spline method arising in engineering. AIMS Maths (2023) 8(6):1474766. 10.3934/math.2023754

  • 10.

    EdwanRSaadehRHadidSAl-SmadiMMomaniS. Solving time-space-fractional Cauchy problem with constant coefficients by finite-difference method. Comput Math Appl (2020) 2546.

  • 11.

    ShqairMGhabarIBurqanA. Using Laplace residual power series method in solving coupled fractional neutron diffusion equations with delayed neutrons system. Fractal Fract (2023) 7:219. 10.3390/fractalfract7030219

  • 12.

    AboanberAENahlaAAAljawaznehSM. Fractional two energy groups matrix representation for nuclear reactor dynamics with an external source. Ann Nucl Energ (2021) 153:108062. 10.1016/j.anucene.2020.108062

  • 13.

    SardarTSaha RaySBeraRKBiswasBBDasS. The solution of coupled fractional neutron diffusion equations with delayed neutrons. Int J Nucl Energ Sci Technol (2010) 5(2):10513. 10.1504/ijnest.2010.030552

  • 14.

    El-NabulsiRA. Fractal neutrons diffusion equation: Uniformization of heat and fuel burn-up in nuclear reactor. Nucl Eng Des (2021) 380:111312. 10.1016/j.nucengdes.2021.111312

  • 15.

    El-NabulsiRA. Nonlocal effects to neutron diffusion equation in a nuclear reactor. J Comput Theor Transport (2020) 49(6):26781. 10.1080/23324309.2020.1816551

  • 16.

    El-NabulsiRA. Neutrons diffusion variable coefficient advection in nuclear reactors. Int J Adv Nucl Reactor Des Technol (2021) 3:1027. 10.1016/j.jandt.2021.06.005

  • 17.

    PatraAShoneTTMishraBB. Natural decomposition approximation solution for first order nonlinear differential equations. Int J Eng Technol (2018) 7(4):4425. 10.14419/ijet.v7i4.5.20202

  • 18.

    NahlaAAAl-MalkiFARokayaM. Numerical techniques for the neutron diffusion equations in the nuclear reactors. Adv Stud Theor Phys (2012) 6(14):64964.

  • 19.

    TumeleroFLapaCMBodmannBEVilhenaMT. Analytical representation of the solution of the space kinetic diffusion equation in a one-dimensional and homogeneous domain. Braz J Radiat Sci (2019) 7(2B):7. 10.15392/bjrs.v7i2b.389

  • 20.

    KhaledSM. Exact solution of the one-dimensional neutron diffusion kinetic equation with one delayed precursor concentration in Cartesian geometry. AIMS Math (2022) 7:1236473. 10.3934/math.2022686

  • 21.

    El-AjouAAl-SmadiMOqielatMNMomaniSHadidS. Smooth expansion to solve high-order linear conformable fractional PDEs via residual power series method: Applications to physical and engineering equations. Ain Shams Eng J (2020) 11(4):124354. 10.1016/j.asej.2020.03.016

  • 22.

    El-AjouA. Taylor’s expansion for fractional matrix functions: Theory and applications. J Math Comput Sci (2020) 21(1):117. 10.22436/jmcs.021.01.01

  • 23.

    OqielatMNEriqatTAl-ZhourZOgilatMNEl-AjouAHashimI. Construction of fractional series solutions to nonlinear fractional reaction–diffusion for bacteria growth model via Laplace residual power series method. Int J Dyn Control (2023) 11(2):5207. 10.1007/s40435-022-01001-8

  • 24.

    OqielatMNEriqatTOgilatOEl-AjouAAlhazmiSEAl-OmariS. Laplace-residual power series method for solving time-fractional reaction–diffusion model. Fractal and Fractional (2023) 7(4):309. 10.3390/fractalfract7040309

  • 25.

    SaadehRBurqanAEl-AjouA. Reliable solutions to fractional Lane-Emden equations via Laplace transform and residual error function. Alexandria Eng J (2022) 61(12):1055162. 10.1016/j.aej.2022.04.004

  • 26.

    SaadehR. A reliable algorithm for solving system of multi-pantograph equations. WSEAS Trans Math (2022) 21:792800. 10.37394/23206.2022.21.91

  • 27.

    KhresatHEl-AjouAAl-OmariSAlhazmiSEOqielatMA. Exact and approximate solutions for linear and nonlinear partial differential equations via Laplace residual power series method. Axioms (2023) 12(7):694. 10.3390/axioms12070694

  • 28.

    KorpinarZBaleanuDİnçMAlmohsenB. Some applications of the least squares-residual power series method for fractional generalized long wave equations. J Ocean Eng Sci (2021). 10.1016/j.joes.2021.09.001

  • 29.

    SalahEQazzaASaadehREl-AjouA. A hybrid analytical technique for solving multi-dimensional time-fractional Navier-Stokes system. AIMS Maths (2023) 8(1):171336. 10.3934/math.2023088

  • 30.

    El-AjouAAl-ZhourZ. A vector series solution for a class of hyperbolic system of Caputo-time-fractional partial differential equations with variable coefficients. Front Phys (2021) 9:276. 10.3389/fphy.2021.525250

  • 31.

    BurqanAEl-AjouASaadehRAl-SmadiM. A new efficient technique using Laplace transforms and smooth expansions to construct a series solution to the time-fractional Navier-Stokes equations. Alexandria Eng J (2022) 61(2):106977. 10.1016/j.aej.2021.07.020

  • 32.

    El-AjouA. A modification to the conformable fractional calculus with some applications. Alexandria Eng J (2020) 59(4):223949. 10.1016/j.aej.2020.02.003

  • 33.

    El-NabulsiRATorresDF. Necessary optimality conditions for fractional action-like integrals of variational calculus with Riemann–Liouville derivatives of order (α, β). Math Methods Appl Sci (2007) 30(15):19319. 10.1002/mma.879

  • 34.

    CassellJSWilliamsMM. A solution of the neutron diffusion equation for a hemisphere with mixed boundary conditions. Ann Nucl Energ (2004) 31(17):19872004. 10.1016/j.anucene.2004.04.008

  • 35.

    KhasawnehKDababnehSOdibatZ. A solution of the neutron diffusion equation in hemispherical symmetry using the homotopy perturbation method. Ann Nucl Energ (2009) 36(11-12):17117. 10.1016/j.anucene.2009.09.001

  • 36.

    SoodAForsterRAParsonsDK. Analytical benchmark test set for criticality code verification. Prog Nucl Energ (2003) 42(1):55106. 10.1016/s0149-1970(02)00098-7

  • 37.

    YangLSuHZhongCMengZLuoHLiXet alHyperspectral image classification using wavelet transform-based smooth ordering. Int J Wavelets, Multiresolution Inf Process (2019) 17(06):1950050. 10.1142/s0219691319500504

  • 38.

    GuarigliaE. Harmonic Sierpinski gasket and applications. Entropy (2018) 20(9):714. 10.3390/e20090714

  • 39.

    ZhengXTangYYZhouJ. A framework of adaptive multiscale wavelet decomposition for signals on undirected graphs. IEEE Trans Signal Process (2019) 67(7):1696711. 10.1109/tsp.2019.2896246

  • 40.

    GuarigliaESilvestrovS. Fractional-Wavelet Analysis of Positive definite Distributions and Wavelets on D'(C) D′(C). In Engineering mathematics II: Algebraic, stochastic and analysis structures for networks, data classification and optimization 2016. Springer International Publishing (2017). p. 33753.

  • 41.

    GuidoRCPedrosoFContrerasRCRodriguesLCGuarigliaENetoJS. Introducing the Discrete Path Transform (DPT) and its applications in signal analysis, artefact removal, and spoken word recognition. Digital Signal Process. (2021) 117:103158. 10.1016/j.dsp.2021.103158

Summary

Keywords

neutron diffusion, power series, radial flux, radiation boundary condition, nuclear reactor physics

Citation

El-Ajou A, Shqair M, Ghabar I, Burqan A and Saadeh R (2023) A solution for the neutron diffusion equation in the spherical and hemispherical reactors using the residual power series. Front. Phys. 11:1229142. doi: 10.3389/fphy.2023.1229142

Received

25 May 2023

Accepted

23 August 2023

Published

11 September 2023

Volume

11 - 2023

Edited by

Rami Ahmad El-Nabulsi, Chiang Mai University, Thailand

Reviewed by

Marwan Alquran, Jordan University of Science and Technology, Jordan

Yusif Gasimov, Azerbaijan University, Azerbaijan

Updates

Copyright

*Correspondence: Rania Saadeh,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics