Analytical solutions are derived for simple neutrondiffusion problems in one neutron energy group in systems of simple geometry. The derivation of diffusion equation is based on ficks law which is derived under many assumptions. A new analytical solution for the 2d advectiondispersion. With the closedform solutions of the diffusion equation established in the present study, we also derived the mathematical conditions under which the diffusion and surface. A comparison of some numerical methods for the advection. Nuclear scientists and engineers often need to know where neutrons are in an apparatus, what direction they are going, and how quickly they are moving. Neutron diffusion 90 if we insert the diffusion approximation 23 into our balance equation 4, we obtain. The plots all use the same colour range, defined by vmin and vmax, so it doesnt matter which one we pass in the first argument to lorbar the state of the system is plotted as an image at four different stages of its evolution. Solution of onegroup neutron diffusion equation for. Pdf the spacetime neutron diffusion equations with multigroup of delayed neutrons are a couple of the stiff nonlinear partial differential. Need to solve the transport equation for an accurate estimate of local reaction rates, particularly near fuel rods and strong neutron. Pdf averaging the neutron diffusion equation researchgate. Derivation of the neutron diffusion equation nuclear.
To show the validity of these diffusion boundary conditions, we solve a sample problem. Fundamental concepts and language diffusion mechanisms vacancy diffusion interstitial diffusion impurities. The paper deals in its first part with the general formulation of the convectivediffusion equation and with the numerical solution of this equation by means of. For the love of physics walter lewin may 16, 2011 duration. We now add a convection term \ \boldsymbolv\cdot\nabla u \ to the diffusion equation to obtain the wellknown convectiondiffusion equation. Application of the finite element method to the three. Approximation of the neutron diffusion equation on. Derivation of diffusion equation the diffusion equation 5. The helmholtz equation is derived, and the limitations on diffusion equation as well as the boundary conditions used in. Neutron transport is the study of the motions and interactions of neutrons with materials. 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. Depending on context, the same equation can be called the advectiondiffusion equation, driftdiffusion equation, or.
Different stages of the example should be displayed, along with prompting messages in. Multigroup diffusion 6 this work is detailed in garland1975 but for the present discussion, the main point to note is the inadequacy of the onegroup model or even the twogroup model since the appropriate cross sections are not explicitly available and. The neutron diffusion equation can be solved analytically in academic cases or using standard numerical analysis techniques such as the. Diffusion equation linear diffusion equation eqworld. 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. Understand origin, limitations of neutron diffusion from. Derivation of the neutron diffusion equation introduction. The neutron diffusion equation is often used to perform corelevel neutronic calculations. Twogroup diffusion theory and the approximate representation of. To run this example from the base fipy directory, type. A comparison of some numerical methods for the advectiondi. Solutions of diffusion equation 2461 to a great extent, this lack of consistency is due to the limited information that can be derived from the infiniteseries solutions. An approximate analytical solution to the diffusion. It is commonly used to determine the behavior of nuclear reactor cores and experimental or industrial neutron beams.
By denoting the external boundary of a domain r by aer, eq. Pdf numerical techniques for the neutron diffusion equations in. Steadystate diffusion when the concentration field is independent of time and d is independent of c, ficks second law is reduced to laplaces equation, 2c 0 for simple geometries, such as permeation through a thin membrane, laplaces equation can be solved by integration. Lecture 12 neutron diffusion equation objectives in this. It is revealed that the diffusivity corresponds to the angular momentum operator in quantum theory. In previous section we dealt with the multiplication system and we defined the infinite and finite multiplication factor. To satisfy this condition we seek for solutions in the form of an in nite series of. Nonlinear reaction diffusion equation with michaelismenten kinetics and adomian decomposition method. This is the energydependent neutron diffusion equation. Summary the neutron diffusion equation is a basic balance equation that describes the transport of neutrons in space, energy, and time. This notebook is an entirely selfcontained solution to a basic neutron diffision equation for a reactor rx made up of a single fuel rod.
It is not exact, but for most of this course it is the model that we will use to describe the behavior. The general solution is composed by sum of the general integral of the associated homogeneous equation and the particular solution. In other words, we assume that the lateral surface of the bar is perfectly insulated so no heat can be gained or lost through it. In the present study, the imaginary time is incorporated into the diffusion equation for understanding of the collision problem between two micro particles. Before attempting to solve the equation, it is useful to understand how the analytical. The diffusion equation to derive the homogeneous heatconduction equation we assume that there are no internal sources of heat along the bar, and that the heat can only enter the bar through its ends. Unstructured grids and the multigroup neutron diffusion. The main calculation method explored in this chapter is the neutrondiffusion equation. This work introduces the alternatives that unstructured grids can provide. The convectiondiffusion equation convectiondiffusion without a force term. Solution of the transport equations using a moving coordinate system ole krogh jensen and bruce a. Equation 19 is a nonhomogeneous ordinary differential equation that can be solved by the application of classical methods.
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. This implies that the diffusion theory may show deviations from a more accurate solution of the transport equation in the proximity of external. The neutron diffusion equation describes the neutron population in a nuclear re actor core. We let cx,y,z,t be the density mass per unit volume of a diffusing substance x, and let e be any small subregion of the region where diffusion is occurring. The nal piece of the puzzle requires the use of an empirical physical principle of heat ow. Finlayson department of chemical engineering, university of washington, seattle, washington 98195. This problem contains no information about the spatial distribution of neutrons, because it is a point geometry problem.
Here, a finite difference method is used to discretize the neutron diffusion equation pazirandeh et al. The diffusion equation parabolic d is the diffusion coefficient is such that we ask for what is the value of the field wave at a later time t knowing the field at an initial time t0 and subject to some specific boundary conditions at all times. Equation 1 is known as a onedimensional diffusion equation, also often referred to as a heat equation. The convectiondiffusion equation is a combination of the diffusion and convection equations, and describes physical phenomena where particles, energy, or other physical quantities are transferred inside a physical system due to two processes. We shall then derive the neutron balance equation with diffusion. The helmholtz equation is derived, and the limitations on diffusion equation as well as the boundary conditions used in its application to realistic. Pdf this paper presents a general theoretical analysis of the neutron motion problem in a nuclear reactor, where large variations on neutron. Solution of the neutron diffusion equation by the finite element method in the general multigroup formalism, the neutron diffusion equation is represented by a coupled system of differential equations on the scalar flux 4, where the notation is conventional. Define i in zzo, 5 xzo, sz zol, n abk, and 5, a cu,kq, where q is the source strength cm sl. Pdf solution of fixed source neutron diffusion equation.
Chapter 2 the diffusion equation and the steady state. This section was about conditions for a stable, selfsustained fission chain reaction and how to maintain such conditions. Chapter 2 the diffusion equation and the steady state weshallnowstudy the equations which govern the neutron field in a reactor. The convectivediffusion equation is the governing equation of many important transport phenomena in building physics. Solution of the transport equations using a moving. Nonlinear reaction diffusion equation with michaelis. This work deals with this model for nuclear reactors with hexagonal.
Usa received 4 march 1979 a convectiondiffusion equation arises from the conservation equations in miscible and. Secondorder parabolic partial differential equations nonhomogeneous heat diffusion equation 1. We discuss the diffusion equation as an asymptotic limit of the neutron transport equation for large scattering cross sections. The steady state and the diffusion equation the neutron field. The diffusion equation can, therefore, not be exact or valid at places with strongly differing diffusion coefficients or in strongly absorbing media. Diffusion mse 201 callister chapter 5 introduction to materials science for engineers, ch. In general, the substances of interest are mass, momentum. Lou odette american international group aig october 17, 2006 1 nonlinear drift in the continuum limit the pdf. Numerical solution of the neutron diffusion equation has been done by many numerical methods such as the finite difference, finite element and boundary element. The parameter \\alpha\ must be given and is referred to as the diffusion coefficient. The derivation of the diffusion equation will depend on ficks law, even though a direct derivation from the transport equation is also possible. Solution of diffusion equation for a groundlevel area source the first step in the solution is a transformation of equation 3. The results show that solutions of the diffusion equation with these well.
1059 326 224 590 660 987 868 1134 676 1101 1086 1526 459 205 889 1484 501 1313 756 1270 300 1517 166 388 562 208 363 316 716 771 506 92 36 7 329 255