Abstract
We present general analytical solutions to the nuclear dynamics-related Neutron Boltzmann Transport Equation inside nanoenergy reactors. Finding a particular solution to the neutron equation by making use of boundary conditions and initial conditions may be too much for the present study and reduce the generality of the solutions. Some simple assumptions have been introduced in the main system thanks to the Boubaker Polynomial Expansion Scheme, BPES, in order to make the general analytical procedure simple and adaptable for solving similar real-life problems.
1. Introduction
Several choices are possible for describing neutron behavior in a medium filled with nuclei. A neutron is a subatomic particle called a baryon having the characteristics strong force of the standard model. Thus, a quantum mechanical description seems appropriate, leading to an involved system of Schrodinger equations describing neutron motion between and within nuclei [1–7]. A neutron is also a relativistic particle with variation of its mass over time when travelling near the speed of light. Additionally, a neutron posses wave and classical particle properties simultaneously and therefore a collective description like that of Maxwell’s equations also seems appropriate. In reality, a neutron displays all of the above characteristics at one time or another. When a neutron collides with the nucleus, its strong force interacts with all of the individual nucleons [2, 4].
However, between nuclear collisions, neutrons move ballistically. Neutrons with energies above 20 Mev with speeds of more than 20% the speed of light exhibit relativistic motion, but most in a reactor are rarely above 0.17 c. The neutron wavelength is most important for ultra-low-energy neutrons mainly existing in the laboratory. Fortunately, classical neutral particle description with quantum mechanics describing collisions energies is the most appropriate for the investigation of neutron motion within a nuclear reactor. Neutron transport equation, which is also called neutron Boltzmann equation, is an equation which characterizes a relatively small number of neutrons colliding in a vast sea of nuclei [2, 4].
Mathematically, a neutron is a neutral point particle, experiencing deflection from or capture by a nucleus at the center of an atom. If the conditions are just right, the captured neutron causes a fissile nucleus to fission, producing more neutrons [2, 4]. The statistically large number of neutrons interacting in a reactor allows for a continuum-like description through averaging resulting in the linear Boltzmann equation. Also, a statistical mechanics formulation, first attempted by Boltzmann for interacting gases, provides appropriate descriptions.
2. Neutron Transport Equation
The neutron transport equation is a linearized version of the Boltzmann equation with wide applications in physics, geophysics, and astrophysics. The neutron transport equation models the transport of neutral particles in a scattering, fission, and absorption events with no self-interactions [4, 8]. It is used in radiation shielding and nanoenergy reactor calculations, as well as in radiative transfer of stellar and planetary atmospheres, and it also describes dispersion of light, the passage of γ-rays through dispersive media, and so forth [4, 8, 9].
The resolution of problems dealing with transport phenomena is the subject of several works, especially in the context of transfer multidimensional problems based on analytical and numerical approaches. One can refer, for example, to Fourier transform [10–12] and many others using the Laplace transform [13].
Chebyshev spectral methods for radiative transfer problems are also studied, for example, by Kim and Ishimaru in [14], by Kim and Moscoso in [15], and by Kadem in [16]. For more detailed study on Chebyshev spectral method and also approximations by the spectral methods, we refer the reader to monographs by Boyd [17] and Bernardi and Maday [18].
The neutron transport problem has been studied analytically [19] and recently, analytical solutions were sought using spectral methods [20]. However, most of these studies deal with the reduced neutron transport equation, in which one or two of the variables are left out (especially, the time-dependent part). Hence, we shall present an analytical procedure of solving the general neutron transport equation, while keeping all the fundamental variables.
3. Part I
3.1. Analytical Procedure
In this study, all particles including nuclei are in motion along with random collisions [3, 20]. Although several forms of the neutron transport equation exist, the integrodifferential formulation, arguably the most popular form in neutron transport and reactor physics applications, is presented and would be solved analytically for general applications. This equation is given as follows [2, 3, 20, 21]: where is the neutron speed, and are macroscopic cross-sections, is the scattering cross-section, is the distribution function, is the neutron angular flux, and are energies, and and are neutron directions, while is the source function [3].
Employing the method of variable separation (which is somehow similar to the multigroup theory), we shall make the following assumptions: Introducing these expressions in (3.1), we write represents the average number of neutrons per fission. Now, since is actually a vector that represents (or general) coordinates, we may write If we also write It follows that As far as we are concerned, the functions in would be taken as constant variables in (3.8). We may simplify the above problem further if we write Equation (3.8) becomes It is assumed that is separable such that Equation (3.10) can be written in the form If the function is separable such that and also, if the possibility of expressing the other functions, as follows, exists, the problems become very easy to handle: Hence, Since both sides of (3.16) are independent of one another, they must be equal to a constant , leading to the following equations: Using the method of integrating factor or by direct integration, the solutions to (3.17) are as follows: where and are constants.
Therefore, from (3.13), we could write From (3.2), the expression for the flux is given as This solution is very interesting because it leaves us with the choice of different expression for the functions , , , , , and . Hence, the nature of a nuclear reactor, the intended use, and different applications may be easily imposed on the analytical expression for the neutron flux.
However, there is a change that must be overcome, that is, separating the function into a sum of and . Once the possibility of obtaining such separation exists, then our solution is very easy to apply in any real life situations (nuclear reactor design, transport in porous media, fractured/random media analyses, etc.). However, we shall attempt a way around the challenge.
The macroscopic cross-section is given usually as follows: for a given nuclide and reaction type , where is the nuclear atomic density and is the microscopic cross-section.
Comparing (3.21) and (3.4), we see that the function is actually the microscopic cross-section. Hence, if we assume that nuclear atomic density is independent of time: Then This is true for the second macroscopic cross-section of the initial energy, that is, . Therefore, if is taken to be 1, it follows that Hence, (3.9) becomes It would be observed that (3.25) is independent of time and this means that Therefore, the expression for the angular flux becomes
3.2. Expression for the Distribution Function
From (3.26), we may write [2, 3, 20, 21]
4. Part II
4.1. Criticality and Analytical Solutions
The neutron transport equation without delayed neutrons is given as [2, 3, 20, 21] where is the external sources of neutrons, is the microscopic cross-section, is the neutron speed, and is the average number of neutrons per fission.
The equation assumes that all neutrons are emitted instantaneously at the time of fission. In fact, small fraction of neutrons is emitted later due to certain fission products [2, 3].
Now, if we seek as asymptotic solutions to (4.1) in the following form: where the solution satisfies the boundary conditions, and if the system is source free , (4.1) becomes We shall make the following assumptions: Equation (4.3) then becomes In the differential equation (4.6), all variables other than functions of are taken to be constant. Hence, we may write If we suppose that the integral is convergent, and taking into account the characteristics of a given nuclear reactor with spherical symmetry (), within radial range , where and are core reactor characteristic constants.
For solving (4.8), the Boubaker Polynomials Expansion Scheme, BPES, [22–35] is proposed. This scheme is applied through setting the following expression: where are the -order Boubaker polynomials, is the radius (), are minimal positive roots, is a prefixed integer, and are unknown pondering real coefficients.
The main advantage of this step lies in (4.10), which ensures verifying the four boundary conditions in (4.9), at the earliest stage of resolution protocol. In fact, due to the properties of the Boubaker polynomials [25–33], and since are roots of , (4.9) is reduced to The solution is then assigned to the set of pondering real coefficients , which minimizes the Minimum Square functional : which gives the following solution to (4.8): with constant.
From (4.8) and our earlier assumptions, we write Equation (4.16) is the neutron flux which defined the distribution neutron in the reactor core.
4.2. The -Eigenvalue
The -eigenvalue from of the critically problem is formulated by assuming that , the average number of neutrons per fission, can be adjusted to obtain a time-independent solution such as the expression of (4.16). Hence, we replace by [2, 3, 20, 21] so that such a solution is equivalent to (4.11). Therefore, it follows that We may write where
5. Conclusion
The analytical solution of the neutron transport equation applying the method of separation of variable has been presented. It is very interesting to note from (4.18) that the determination of the criticality of a nuclear nanoreactor depends on how term weighs compared to . Three interesting cases can be illustrated as follows [2, 3, 20, 21].(1)If , the system is supercritical.(2)If , the system is critical.(3)If , the system is subcritical.
Although working with these analytical solutions may not be an easy task, it is very crucial to note that the advantages are quite enormous.(1)For instance, we could now see how each parameters involved in nuclear reactor modeling interact with one another and we can now understand the in-depth physics better, giving us the opportunity to manipulate the parameters within natural and experimental restrictions.(2)The solutions have a lot of undefined functions and hence, a nuclear scientist is at liberty of defining such functions as the field of application of neutron transport requires. This may be very helpful in applications related to the Boltzmann equation.(3)We now have analytical expressions that may be used to access the reliability of the computational/numerical tools, which are used in reactor modeling presently. Algorithms designed from any of these expressions are definitely going to be very fast compared to Monte Carlo analysis.(4)Boundary conditions are now easier to use because of the separation of the main functions in the neutron transport equation.
These advantages will be explored one at a time in our next investigation.
Acknowledgment
The authors wish to thank Dr. Teresa Kulikowska for making some important materials available to us in the course of this study.