Numerical Approximation for Fractional Neutron Transport Equation

Zhao, Zhengang; Zheng, Yunying

doi:https://doi.org/10.1155/2021/6676640

Journal of Mathematics

On this page

Abstract Introduction Results Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 6676640 | https://doi.org/10.1155/2021/6676640

Numerical Approximation for Fractional Neutron Transport Equation

Zhengang Zhao¹and Yunying Zheng²

Academic Editor: Imtiaz Ahmad

Received21 Dec 2020

Revised09 Jan 2021

Accepted20 Jan 2021

Published15 Mar 2021

Abstract

Fractional neutron transport equation reflects the anomalous transport processes in nuclear reactor. In this paper, we will construct the fully discrete methods for this type of fractional equation with Riesz derivative, where the generalized WENO5 scheme is used in spatial direction and Runge–Kutta schemes are adopted in temporal direction. The linear stabilities of the generalized WENO5 schemes with different stages and different order ERK are discussed detailed. Numerical examples show the combinations of forward Euler/two-stage, second-order ERK and WENO5 are unstable and the three-stage, third-order ERK method with generalized WENO5 is stable and can maintain sharp transitions for discontinuous problem, and its convergence reaches fifth order for smooth boundary condition.

1. Introduction

The fractional differential operations, such as the natural generalizations of the classical operators of integer order, have become more and more popular. Many phenomena in physics, engineering, biology, chemistry, and other sciences have been described successfully by fractional calculus and fractional differential equations [1]. Recently, several numerical methods to solve fractional differential equations have been developing rapidly. For example, the finite difference method [2–5], the finite element method [6–13], the discontinuous finite element method [14–20], the spectral method [21–24], the finite volume method [25, 26], etc. However, most of them are focusing on the fractional parabolic equation. It seems that there are few works regarding the numerical methods for fractional hyperbolic equation, where the physical phenomena often refer to vibration and oscillation of the complex material or complex progress. Compared to the classical hyperbolic equation, the solutions of fractional hyperbolic equation may also develop discontinuously even when the initial condition is smooth. It is impossible to deduce the analytical solutions in most cases, and one thus has to rely on numerical methods. So, it is very meaningful to study the numerical methods for fractional hyperbolic equation.

As a typical representative of hyperbolic equation, the classical neutron transport equation studies the motions and interactions of neutrons with materials in a nuclear reactor. Besides, the transport equation also plays an important role within the transport theory. Compared to the classical neutron transport equation, the fractional neutron transport equation denotes the anomalous transport processes [27], which was studied by many researchers. Das and Biswas [28] developed a constitutive neutron diffusion equation for describing the neutron flux profile by using the concept of fractional divergence. Vyawahare et al. [29] studied a fractional-order model of a nuclear reactor, which predicts subdiffusive behaviour for long time. Vyawahare and Nataraj [30] dealt with a fractional-order model of neutron transport process and pointed out that it is a more faithful and realistic representation of neutron movements than the classical model.

Weighted essentially nonoscillatory (WENO) schemes [31–34] are the effective numerical methods to solve the hyperbolic equation, specially, the fifth-order WENO (WENO5), which is a robust scheme for solving the classical neutron transport equation. Naturally, applying the WENO schemes to solve the fractional differential equation numerically is an interesting work. However, it seems that there are few works about the WENO schemes used to approximate the fractional differential equation except [3]. In [3], Deng et al. designed a WENO6 scheme for the fractional diffusion equation in the sense of Caputo derivative.

In this work, we will design the generalized WENO5 schemes to solve the following spatial fractional neutron transport equation:where is a given constant coefficient and represents the Riesz derivative given byin which .

The basic idea of numerical methods in this paper is to split the order fractional derivative into a fractional integral and a classical first derivative. Next, we apply the Gauss–Jacobi quadrature to approximate the fractional integral and WENO5 to approximate the classical first derivative separately. Then, the three-stage, third-order SSP (3,3) explicit Runge–Kutta (ERK) method or the third-order TVD ERK method [35] is selected as the discrete time method to solve the ODE systems. The numerical experiments confirm fifth-order convergence order and show the maintenance of sharp discontinuous transitions. Furthermore, the linear stability of various ERK methods together with generalized WENO5 to approximate the fractional neutron transport equation is also discussed. From the linear stability analyses, we get that the combinations of the forward Euler scheme with WENO5 and any two-stage, second-order ERK schemes with WENO5 are linearly unstable in the fractional neutron transport equation. And, the combination of any three-stage, third-order ERK method with WENO5 for the fractional neutron transport equation is linearly stable. All of these theoretical results are demonstrated by numerical experiments.

The structure of this paper is as follows. In Section 2, the generalized WENO5 schemes of the spatial fractional neutron transport equation are designed, including the Gauss quadrature formula relative to Jacobi weights, to compute the fractional integral terms. Then, the linear stability analyses of the designed schemes are discussed in Section 3. Finally, the numerical examples are presented to verify the convergence and stability of the designed schemes in Section 4.

2. The Generalized WENO5 Schemes

To achieve our method, we need to introduce the fractional integrals and fractional derivatives.

Firstly, we define the th order left- and right-fractional integrals of function as follows:where . Next, we give the th order left- and right-Caputo derivatives:in which .

According to [36], when , we can easily derive the following relations between two fractional derivatives:

Therefore, equation (1) can be rewritten as

By using the definition of the Caputo derivative, we havewhere .

In order to get a high-order accuracy of the numerical solution, we apply Gauss-type quadrature formula relative to Jacobi weights to compute the fractional integral terms in (7). Therefore, we transfer the integral intervals and to the standard interval , respectively:

Now, we denote as the spatial step size, , where is the numbers of subintervals divided in the spatial directions. Then, by using points Gauss–Jacobi quadrature and the Lagrange interpolation, the left- and right-fractional integrals in the right-hand side of (8) and (9) at the mesh points can be, respectively, approximated bywhere is the Gauss nodes, and are Jacobi weights, and is the Lagrange interpolation basis function associated with the points , which can be replaced by any other kind of interpolation basis functions.

In this paper, we denote and as the exact solution and the numerical solution at mesh point , respectively. For the intermediate function , we use WENO5 to do the discretization as follows:where the term is the numerical flux and can be split into positive and negative parts:

On stencil, , .

For and , we define them on each stencil as follows:where, for parameters , stencil , and the other details of the WENO5 method, one can refer to [31, 34, 37].

Next, we will use the forward Euler method, SSP (3,3) schemes, and some other Runge–Kutta schemes to solve the ODE system in the temporal direction after doing the spatial discretization.

3. Linear Stability Analysis

In this section, we will discuss the linear stability properties of the generalized WENO5 method when coupled with different kinds of ERK methods for solving the fractional neutron transport equation (1). Here, we define CFL number as and only consider the corresponding homogeneous problem:under homogeneous boundary conditions .

In the following, we will show that the combination of the generalized WENO5 and forward Euler scheme for the fractional neutron transport equation (15) is linearly unstable.

Lemma 1 (see [37]). When the forward Euler method is applied with the classical WENO5 scheme, the function , which is defined as , can be approximated bywhere the Taylor expansions of and are used, is a small positive number, only related to the spatial step , and is the imaginary unit, i.e., .

Theorem 1. The combination of the generalized WENO5 method and the forward Euler method for the fractional neutron transport equation (15) is linearly unstable.

Proof. Here, we only consider the Lax–Frichrichs flux splitting, which means that the negative part of the numerical flux is zero [37]. And, the same analysis and conclusions can be shifted to the other flux splitting.
According to the Fourier transform, we assumewhere , and is the spatial frequency. By using (13) and (17), we further conclude thatWhen the forward Euler method is adopted, the fully discrete generalized WENO5 method becomeswhere , which can be defined as .
Then, equation (15) becomesDefine the amplification factor:By Lemma 1, we get the amplification factor as follows:where and are denoted by and correspondingly.
According to a simple calculation, we can get thatThe proof is complete.
Next, we consider the following stage ERK method to replace the forward Euler method:where the coefficients satisfy , for , and .

Lemma 2 (see [38, 39]). The amplification factor of the stage, order ERK method iswhere are the coefficients of the stage, order ERK method, , and , where is caused by the character of the spatial operator.

Remark that, from Theorem 1, we can see that the operator should be adjusted to in the fractional neutron transport equation, which is caused by the left and right spatial fractional integral terms.

Theorem 2. The combination of the generalized WENO5 method with any two-stage, second-order ERK method for the fractional neutron transport equation (15) is linearly unstable.

Proof. According to Lemmas 1 and 2, the amplification factor can be written aswhere the constant is denoted by .
Then, we have thatThe proof is finished.

Theorem 3. The combination of the generalized WENO5 method with any three-stage, three-order ERK method for the fractional neutron transport equation (15) is linearly stable.

Proof. According to Lemmas 1 and 2 and the proof of Theorem 2, we can get the amplification factor as follows:Then, we have thatThe proof is finished.
From the above proof, we can see that the conclusion of the linear stabilities for the combination of the generalized WENO5 method with the ERK method in the fractional neutron transport equation is the same as in the classical neutron transport equation. The reason is that the added fractional integral terms are approximated by a linear combination, in which the weights are bounded in the Gauss–Jacobi quadrature formulas. Therefore, all of the conclusions about the linear stabilities for the classical hyperbolic conservation law [37] can be shifted to problem (1), where the proofs are omitted.

4. Numerical Examples and Results

In this section, we will design three numerical examples to examine the effectiveness of the numerical algorithms. The first example with periodic boundary conditions is designed to verify the linear stability of the proposed schemes. The second example with two points’ boundary conditions is computed to verify the convergence order. And, the third example is designed for numerical solving piecewise smooth solution and discontinuous initial condition.

Example 1. Consider the following fractional neutron transport equation:where is chosen as the exact solution.
Firstly, we use the forward Euler method with the WENO5 scheme. Next, we use the two-stage, second-order SSP ERK method with the WENO5 scheme. We compute the two cases for , which are shown in Figures 1 and 2 separately. According to the theoretical analyses, spurious oscillations due to linear instabilities are present.
Then, we use the SSP (3,3) ERK method with the WENO5 scheme. error and convergence order are shown in Table 1. From Table 1, we can see that the numerical scheme performs well and shows the fifth-order convergence accuracy.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

Example 2. Consider the following fractional neutron transport equation:with the initial condition, two points’ boundary conditions,and the source term iswhere . Then, the exact solution is .
Applying the SSP (3,3) method with the WENO5 scheme, we get error and convergence order in Table 2. And the fifth-order convergence of the combination of the SSP (3,3) method with the WENO5 scheme is still observed.

Example 3. Suppose that . Consider the fractional neutron transport equation as follows:with the boundary conditions , and the initial conditionIn the simulation of Figure 3, we apply the SSP (3, 3) ERK method with the WENO5 scheme and take . Figure 3 shows numerical solutions with different orders of fractional derivative, at . It can be seen that sharp transitions are retained in the numerical simulations for different .
Figure 4 shows the numerical solutions and their absolute errors for and 200, and at . It can be seen that the bigger the value, the sharper the transitions in the simulations.
Figure 5 shows the numerical solutions when the parameters are chosen as , and separately.
Figure 6 shows the numerical solutions when the parameters are chosen as , at separately.

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

(a)

(b)

(c)

(d)

5. Conclusion

This paper develops the generalized WENO5 methods for numerically solving the spatial fractional neutron transport equation with smooth solutions or piecewise smooth solutions, which is a new attempt in the work of numerical solution for fractional differential equation. By using the classical WENO5 schemes and the Gauss–Jacobi quadrature, we perform the spatial discretization. Then, we get the nonlinear ODE systems, which are computed by the SSP (3,3) Runge–Kutta method and some other Runge–Kutta method. The linear stabilities of the numerical schemes are proved. And, the numerical results are present to verify the effectiveness and reliability of the numerical schemes.

Generally speaking, the WENO scheme is widely used to solve the problem that the solution of the equation is not very smooth and can maintain a high-order convergence rate as well. Therefore, for the upcoming work, we will mainly discuss the discontinuous problems and the nonlinear problems with the WENO method in the future. Maybe, by using the nonuniform mesh method or meshless method together with the WENO scheme for the fractional problems is a good choice.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors have equally contributed to the article and read and approved the final version of the paper.

Acknowledgments

This work was supported by the Natural Science Foundation of Shanghai (no. 19ZR1422000), Natural Science Foundation of Anhui Province (no. 2008085MA11), and Natural Science Foundation of Education Department of Anhui Province (no. KJ2018A0385).

References

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, Amsterdam, Netherlands, 2006.
M. M. Meerschaert and C. Tadjeran, “Finite difference approximations for fractional advection-dispersion flow equations,” Journal of Computational and Applied Mathematics, vol. 172, no. 1, pp. 65–77, 2004.
View at: Publisher Site | Google Scholar
W. Deng, S. Du, and Y. Wu, “High order finite difference WENO schemes for fractional differential equations,” Applied Mathematics Letters, vol. 26, no. 3, pp. 362–366, 2013.
View at: Publisher Site | Google Scholar
H. Sun, Z.-Z. Sun, and G.-H. Gao, “Some temporal second order difference schemes for fractional wave equations,” Numerical Methods for Partial Differential Equations, vol. 32, no. 3, pp. 970–1001, 2016.
View at: Publisher Site | Google Scholar
Z. Hao and W. Cao, “An improved algorithm based on finite difference schemes for fractional boundary value problems with nonsmooth solution,” Journal of Scientific Computing, vol. 73, no. 1, pp. 395–415, 2017.
View at: Publisher Site | Google Scholar
W. H. Deng, “Finite element method for the space and time fractional Fokker-Planck equation,” SIAM Journal on Numerical Analysis, vol. 47, pp. 204–226, 2008.
View at: Google Scholar
V. J. Ervin, N. Heuer, and J. P. Roop, “Numerical approximation of a time dependent, nonlinear, space‐fractional diffusion equation,” SIAM Journal on Numerical Analysis, vol. 45, no. 2, pp. 572–591, 2007.
View at: Publisher Site | Google Scholar
C. Li, Z. Zhao, and Y. Chen, “Numerical approximation of nonlinear fractional differential equations with subdiffusion and superdiffusion,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 855–875, 2011.
View at: Publisher Site | Google Scholar
W. Bu, Y. Tang, and J. Yang, “Galerkin finite element method for two-dimensional Riesz space fractional diffusion equations,” Journal of Computational Physics, vol. 276, pp. 26–38, 2014.
View at: Publisher Site | Google Scholar
B. Jin, R. Lazarov, J. Pasciak, and W. Rundell, “Variational formulation of problems involving fractional order differential operators,” Mathematics of Computation, vol. 84, no. 296, pp. 2665–2700, 2015.
View at: Publisher Site | Google Scholar
Z. Zhang, W. Deng, and G. E. Karniadakis, “A Riesz basis Galerkin method for the tempered fractional Laplacian,” SIAM Journal on Numerical Analysis, vol. 56, no. 5, pp. 3010–3039, 2018.
View at: Publisher Site | Google Scholar
Z. Zhao, Y. Zheng, and P. Guo, “A Galerkin finite element scheme for time-space fractional diffusion equation,” International Journal of Computer Mathematics, vol. 93, no. 7, pp. 1212–1225, 2015.
View at: Publisher Site | Google Scholar
Z. Zhao, Y. Zheng, and P. Guo, “A galerkin finite element method for a class of time-space fractional differential equation with nonsmooth data,” Journal of Scientific Computing, vol. 70, no. 1, pp. 386–406, 2017.
View at: Publisher Site | Google Scholar
Q. Xu and J. S. Hesthaven, “Discontinuous Galerkin method for fractional convection-diffusion equations,” SIAM Journal on Numerical Analysis, vol. 52, no. 1, pp. 405–423, 2014.
View at: Publisher Site | Google Scholar
X. Tian and Q. Du, “Nonconforming discontinuous Galerkin methods for nonlocal variational problems,” SIAM Journal on Numerical Analysis, vol. 53, no. 2, pp. 762–781, 2015.
View at: Publisher Site | Google Scholar
Y. Du, Y. Liu, H. Li, Z. Fang, and S. He, “Local discontinuous Galerkin method for a nonlinear time-fractional fourth-order partial differential equation,” Journal of Computational Physics, vol. 344, pp. 108–126, 2017.
View at: Publisher Site | Google Scholar
Y. Zheng and Z. Zhao, “The discontinuous Galerkin finite element method for fractional cable equation,” Applied Numerical Mathematics, vol. 115, pp. 32–41, 2017.
View at: Publisher Site | Google Scholar
C. Li and Z. Wang, “The local discontinuous Galerkin finite element methods for Caputo-type partial differential equations: numerical analysis,” Applied Numerical Mathematics, vol. 140, pp. 1–22, 2019.
View at: Publisher Site | Google Scholar
Y. Zheng, Z. Zhao, and Y. Cui, “The discontinuous Galerkin finite element approximation of the multi-order fractional initial problems,” Applied Mathematics and Computation, vol. 348, pp. 257–269, 2019.
View at: Publisher Site | Google Scholar
Y. Zheng and Z. Zhao, “The time discontinuous space-time finite element method for fractional diffusion-wave equation,” Applied Numerical Mathematics, vol. 150, pp. 105–116, 2020.
View at: Publisher Site | Google Scholar
X. J. Li and C. J. Xu, “Existence and uniqueness of the weak solution of the space-time fractional diffusion equation and a spectral method approximation,” Communications in Computational Physics, vol. 8, pp. 1016–1051, 2010.
View at: Google Scholar
F. Zeng, F. Liu, C. Li, K. Burrage, I. Turner, and V. Anh, “A crank--nicolson ADI spectral method for a two-dimensional Riesz space fractional nonlinear reaction-diffusion equation,” SIAM Journal on Numerical Analysis, vol. 52, no. 6, pp. 2599–2622, 2014.
View at: Publisher Site | Google Scholar
Z. Mao and G. E. Karniadakis, “A spectral method (of exponential convergence) for singular solutions of the diffusion equation with general two-sided fractional derivative,” SIAM Journal on Numerical Analysis, vol. 56, no. 1, pp. 24–49, 2018.
View at: Publisher Site | Google Scholar
N. Wang, Z. Mao, C. Huang, and G. E. Karniadakis, “A spectral penalty method for two-sided fractional differential equations with general boundary conditions,” SIAM Journal on Scientific Computing, vol. 41, no. 3, pp. A1840–A1866, 2019.
View at: Publisher Site | Google Scholar
H. Hejazi, T. Moroney, and F. Liu, “Stability and convergence of a finite volume method for the space fractional advection-dispersion equation,” Journal of Computational and Applied Mathematics, vol. 255, pp. 684–697, 2014.
View at: Publisher Site | Google Scholar
A. Simmons, Q. Yang, and T. Moroney, “A finite volume method for two-sided fractional diffusion equations on non-uniform meshes,” Journal of Computational Physics, vol. 335, pp. 747–759, 2017.
View at: Publisher Site | Google Scholar
R. Klages, G. Radons, and I. M. Sokolov, Anomalous Transport, Wiley-VCH Verlag GmbH & Co., New York, NY, USA, 2008.
S. Das and B. B. Biswas, “Fractional divergence for neutron flux profile in nuclear reactors,” International Journal of Nuclear Energy Science and Technology, vol. 3, no. 2, pp. 139–159, 2007.
View at: Publisher Site | Google Scholar
V. A. Vyawahare, P. S. V. Nataraj, G. Espinosa-Paredes, and R.-I. Cázares-Ramírez, “Nuclear reactor with subdiffusive neutron transport: development of linear fractional-order models,” International Journal of Dynamics and Control, vol. 5, no. 4, pp. 1182–1200, 2017.
View at: Publisher Site | Google Scholar
V. A. Vyawahare and P. S. V. Nataraj, “Fractional-order modeling of neutron transport in a nuclear reactor,” Applied Mathematical Modelling, vol. 37, no. 23, pp. 9747–9767, 2013.
View at: Publisher Site | Google Scholar
C. W. Shu, “Efficient implementation of essentially non-oscillatory shock-captring schemes,” Journal of Computational Physics, vol. 77, pp. 439–471, 1998.
View at: Google Scholar
X.-D. Liu, S. Osher, and T. Chan, “Weighted essentially non-oscillatory schemes,” Journal of Computational Physics, vol. 115, no. 1, pp. 200–212, 1994.
View at: Publisher Site | Google Scholar
G.-S. Jiang and C.-W. Shu, “Efficient implementation of weighted ENO schemes,” Journal of Computational Physics, vol. 126, no. 1, pp. 202–228, 1996.
View at: Publisher Site | Google Scholar
C.-W. Shu, “High order weighted essentially nonoscillatory schemes for convection dominated problems,” SIAM Review, vol. 51, no. 1, pp. 82–126, 2009.
View at: Publisher Site | Google Scholar
R. J. Spiteri and S. J. Ruuth, “A new class of optimal high-order strong-stability-preserving time discretization methods,” SIAM Journal on Numerical Analysis, vol. 40, no. 2, pp. 469–491, 2002.
View at: Publisher Site | Google Scholar
I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, CA, USA, 1999.
R. Wang and R. J. Spiteri, “Linear instability of the fifth-order WENO method,” SIAM Journal on Numerical Analysis, vol. 45, no. 5, pp. 1871–1901, 2007.
View at: Publisher Site | Google Scholar
E. Tadmor, “Stability analysis of finite difference, pseudospectral and Fourier-Galerkin approximations for time-dependent problems,” SIAM Review, vol. 29, pp. 469–491, 1987.
View at: Publisher Site | Google Scholar
D. Levy and E. Tadmor, “From semidiscrete to fully discrete: stability of runge--kutta schemes by the energy method,” SIAM Review, vol. 40, no. 1, pp. 40–73, 1998.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Zhengang Zhao and Yunying Zheng. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies