In order to describe the distribution of neutrons in a highly heterogeneous configuration, it was necessary to extend the classical neutron diffusion equation. Numerical solution of the neutron diffusion equation. Pdf the solution of twodimensional neutron diffusion equation. Figure 6 shows an illustrative 5 group approximation. The neutron diffusion equation can be solved analytically in academic cases or using standard numerical analysis techniques such as the. The neutron diffusion equation describes the neutron population in a nuclear re actor core. They can be used to solve for the diffusion coefficient, d. Transport crosssection the effect of the scattering angular distribution on. A thermoelectric generator containing plutonium dioxide, used as a source of thermal and electric power in spacecraft, was studied. Davison provides an authoritative, and in many places an elegant description of nearly all the known methods of solving the transport equation. Apr 06, 2008 because diffusion is composed by a sum of different terms, and if your temrs assorbtion and generation that describe the nuclear reaction are 0, there still a term that describe the neutron collision and this term involves the flick law so there still neutron diffusion also without nuclear reacion.
This will yield the required neutron balance equation for the reactor. Each term represents a gain or a loss of a neutron, and the balance, in essence, claims that neutrons gained equals neutrons lost. Multilevel method to compute the lambda modes of the. Modifying the neutron diffusion equation using spatial. Diffusion approximation to neutron transport equation with. It must be added the constants a and c cannot be obtained from the diffusion equation, because these constant give the absolute value of neutron flux.
This is the energydependent neutron diffusion equation. Neutron diffusion equation an overview sciencedirect topics. An iteration method has been implemented to solve a neutron transport equation in a multigroup diffusion approximation. This is usually done by treating the neutron motion as a diffusion process and solving the diffusion equation numerically. In the analysis of the neutronic behaviour of a nuclear reactor, one of the most relevant parameters is the determining of the neutron flux in any region of the reactor core, as a precise assessment of this neutron flux will all. The diffusion equation can, therefore, not be exact or valid at places with strongly differing diffusion coefficients or in strongly absorbing media. This work deals with this model for nuclear reactors with hexagonal. Two step functions, properly positioned, can be summed to give a solution for finite layer placed between two semiinfinite bodies. Solution of neutrontransport multigroup equations system. Transport crosssection the effect of the scattering angular distribution on the motion of a neutron is taken into. The derivation of the diffusion equation will depend on ficks law, even though a direct derivation from the transport equation is also possible. Download file pdf nuclear reactor analysis duderstadt solutions manual nuclear reactor analysis duderstadt solutions manual 23.
This notebook is an entirely selfcontained solution to a basic neutron diffision equation for a reactor rx made up of a single fuel rod. Abstractwe study iterative methods for solving linear systems arising in the discretization of the time dependent neutron diffusion equation. This implies that the diffusion theory may show deviations from a more accurate solution of the transport equation in the proximity of external. Equation 9 is an eigenproblem, whose solution behaves in a typical way. This leads to the steady state diffusion equation which, in simple terms, states that neutron production is equal to neutron absorption plus neutron leakage for stable conditions to be maintained.
The onegroup diffusion equation that we will be stepping through time and space is. The theory of neutron transport developed rapidly during and just after the war, yet no comprehensive account of the theory has appeared until now. It is not exact, but for most of this course it is the model that we will use to describe the behavior. Unstructured grids and the multigroup neutron diffusion equation. When the diffusion equation is linear, sums of solutions are also solutions. Department of appliedmath, universidade federal do rio grande do sul, porto alegre, brazil. J,is the divergence operator x y z s is the source strength is the neutron flux neutronscm. Dont forget to sign up for the class to get really nonintrusive updates. So we will have to use some average flux and cross section that have been averaged over the property in the group energy range in question. This approximation is analogous to ficks law in species diffusion and to fouriers law in heat transfer. In nonmultiplying environment neutrons are emitted by a neutron source situated in the center of coordinate system and then they freely diffuse through media. Pdf the distribution of neutron population in nuclear reactor is described by using transport equations. Pdf numerical techniques for the neutron diffusion. Unstructured grids and the multigroup neutron diffusion.
The continuity equation for the neutron density n neutrons cm3 is. Find materials for this course in the pages linked along the left. Diffusion mse 201 callister chapter 5 introduction to materials science for engineers, ch. Fundamental concepts and language diffusion mechanisms vacancy diffusion interstitial diffusion impurities. The steadystate diffusion equation 3 substituting the source term from eq. Numerical techniques for the neutron diffusion equations in the nuclear reactors article pdf available in advanced studies in theoretical physics 6 january 2012 with 585 reads. First kind of chebyshev polynomials, neutron transport equation, diffusion approximation birinci tip chebyshev polinomlar. The present work evaluates a new version of neutron diffusion equation which is established on the fractional derivatives. Nuclear reactor theory encyclopedia of life support. The objective of this section is to compare the results obtained using the neutron diffusion hybrid equation with the results obtained using the diffusion equation defined as a fixed source problem during the startup of the reactor, as well as with the eigenvalue problem near criticality. Numerical techniques for the neutron di usion equations in. Understand origin, limitations of neutron diffusion from.
Equation 1 is known as a onedimensional diffusion equation, also often referred to as a heat equation. Chapter 2 the diffusion equation and the steady state weshallnowstudy the equations which govern the neutron field in a reactor. Here is an example that uses superposition of errorfunction solutions. This article provides some background to the mathematics of the neutron diffusion process and critical mass. The neutron fractional diffusion equation nfde can simulate a reactor core using nonlocal gradient. It is one of the computer codes maintained or developed by the nuclear engineering division. The hideous neutron transport equation has been reduced to a simple oneliner neutron diffusion equation. Solution of the multigroup neutron diffusion equations by. Heat or diffusion equation in 1d university of oxford.
Diffusion equation laboratory for reactor physics and systems behaviour neutronics comments 1 domain of application of the diffusion equation, very wide describes behaviour of the scalar flux not just the attenuation of a beam equation mathematically similar to. View notes lecture 5 neutron diffusion equation3 from mie 407 at university of toronto. With only a firstorder derivative in time, only one initial condition is needed, while the secondorder derivative in space leads to a demand for two boundary conditions. Pdf averaging the neutron diffusion equation researchgate. The diffusion equation 27 is a partial differential equation of the parabolic type. The derivation of diffusion equation is based on ficks law which is derived under many assumptions. The neutron transport equation is a balance statement that conserves neutrons. In fact the neutron flux can have any value and the critical reactor can operate at any power level. During the diffusion process some neutrons ar e absorbed and any which eventually cross the reactor boundary or surface are lost by leakag e.
Chapter 7 the diffusion equation the diffusionequation is a partial differentialequationwhich describes density. Neutron diffusion 90 if we insert the diffusion approximation 23 into our balance equation 4, we obtain. Advanced monte carlo for radiation physics, particle transport simulation and applications. Physical assumptions we consider temperature in a long thin wire of constant cross section and homogeneous material. It must be added the constant a cannot be obtained from this diffusion equation, because this constant gives the absolute value of neutron flux. Multigroup diffusion 7 recall that the cross sections and flux can vary greatly as a function of neutron energy, e. On the buildup factor from the multigroup neutron diffusion. The general mathematical model of neutron transport is provided by the linear boltzmanns transport equation and the thesis begins with its precise mathematical formulation and presentation of known con. One of possible approximations of neutron transport equation is.
The solution of twodimensional neutron diffusion equation with. The parameter \\alpha\ must be given and is referred to as the diffusion coefficient. The neutron diffusion equation is often used to perform corelevel neutronic calculations. The nfde uses a tensor of diffusion coefficients in different directions. Finite element method applied to neutron diffusion. The subject of this work is computational modeling of neutron transport relevant to economical and safe operation of nuclear facilities. Lecture 5 neutron diffusion equation3 mie 407h1f\1129h. It turns out that this set can be created by convolving the image with gaussian functions of dif ferent scales. In fact, in deriving the neutron diffusion equation it is not essential to assume energy independence of the cross sections, and the. The derivation of the diffusion equation will depend on ficks law, even though a direct derivation from the. Proceedings of the monte carlo 2000 conference, lisbon, 2326 october 2000.
Neutron balance and the diffusion equation by integrating the transport equation over all angles, we obtain an equation for the scalar flux density which can be solved over the entire domain. The neutron rate rate of formation constant is equal to v neut 1 f, where represents secondary neutrons created by ssion the 1 accounts for the neutron causing ssion being consumed in the reaction, and f represents neutron ssion. Everyone breathes a sigh of relief as it is shown to be very solvable, and a criticality relation a balance between neutrons created and destroyed links the geometry of a reactor to its material of construction. Dif3ds nodal option solves the multigroup steadystate neutron diffusion and for cartesian geometry only transport equations in two and threedimensional hexagonal and cartesian geometries. A characteristic for the element interpolation function is m. The solution of twodimensional neutron diffusion equation with delayed neutrons 341 table 3 dependence of solution on physical properties. Approximation of the neutron diffusion equation on. Iterative schemes for the neutron diffusion equation upv. This work introduces the alternatives that unstructured grids can provide. In previous section it has been considered that the environment is nonmultiplying.
The distribution of neutron population in nuclear reactor is described by using transport equations. The solution of twodimensional neutron diffusion equation. Its part of my homework to find an analytical solution to a slab reactor with a reflector. Chapter 2 the diffusion equation and the steady state. We consider the time dependent neutron diffusion equation for one energy group in cylinder coordinates, assuming translational symmetry along the cylinder axis. The multigroup neutron diffusion to equations1 space. Laboratory for reactor physics and systems behaviour neutronics. Diffusion equation and neutron diffusion theory physics forums. A description is given of a program for the ferranti mercury computer which solves the onedimensional multigroup diffusion equations in plane, cylindrical or spherical geometry, and also approximates automatically a twodimensional solution by separating the space.
We are now prepared to consider neutron diffusion in multiplying system, which contains fissionable nuclei i. The gaussian function is the greens function of the linear diffusion equation. In general, the reactor problem in the presence newtonian temperature feedback e. The helmholtz equation is derived, and the limitations on diffusion equation as well as the boundary conditions used in its application to realistic engineering and physics problems are discussed. Reconstruction of the neutron flux in a slab reactor. On the buildup factor from the multigroup neutron diffusion equation with cylindrical symmetry. Pdf this paper presents a general theoretical analysis of the neutron motion problem in a nuclear reactor, where large variations on neutron. Finite element method applied to neutron diffusion problems. The multigroup neutron diffusion equations1 space dimension sven linde summary. In general, the reactor problem in the presence newtonian temperature feedback e ects comprises a very large and.
Solving the neutron diffusion equation, and criticality relations mit 22. These equations are based ontheconceptoflocal neutron balance, which takes int file. It consists of a set of secondorder partial differential equations over the spatial coordinates that are, both in the academia and in the industry, usually solved by discretizing the neutron leakage term using a structured grid. Calculation of the power factor using the neutron diffusion. This equation is an approximation of the neutron transport equation relying on the assumption that the neutron current is proportional to the gradient of the neutron flux by means of a diffusion coefficient. Ficks laws of diffusion describe diffusion and were derived by adolf fick in 1855.
1113 319 375 43 234 1269 77 1434 556 195 980 1473 1409 363 281 741 1533 504 1305 805 1261 669 1574 848 1179 1227 571 737 495 1192 838 448 1088 1314 235 144 925