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. In the present study, the imaginary time is incorporated into the diffusion equation for understanding of the collision problem between two micro particles. 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. By denoting the external boundary of a domain r by aer, eq. Chapter 2 the diffusion equation and the steady state weshallnowstudy the equations which govern the neutron field in a reactor. The main calculation method explored in this chapter is the neutrondiffusion equation. For the love of physics walter lewin may 16, 2011 duration. The helmholtz equation is derived, and the limitations on diffusion equation as well as the boundary conditions used in. Lecture 12 neutron diffusion equation objectives in this. Diffusion mse 201 callister chapter 5 introduction to materials science for engineers, ch. Pdf averaging the neutron diffusion equation researchgate. A comparison of some numerical methods for the advectiondi. In general, the substances of interest are mass, momentum.
A new analytical solution for the 2d advectiondispersion. Solution of diffusion equation for a groundlevel area source the first step in the solution is a transformation of equation 3. To satisfy this condition we seek for solutions in the form of an in nite series of. 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. An approximate analytical solution to the diffusion. The neutron diffusion equation is often used to perform corelevel neutronic calculations. 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. This work introduces the alternatives that unstructured grids can provide. Define i in zzo, 5 xzo, sz zol, n abk, and 5, a cu,kq, where q is the source strength cm sl. We shall then derive the neutron balance equation with diffusion. Depending on context, the same equation can be called the advectiondiffusion equation, driftdiffusion equation, or.
Analytical solutions are derived for simple neutrondiffusion problems in one neutron energy group in systems of simple geometry. Here, a finite difference method is used to discretize the neutron diffusion equation pazirandeh et al. 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. The steady state and the diffusion equation the neutron field. The neutron diffusion equation describes the neutron population in a nuclear re actor core. Derivation of diffusion equation the diffusion equation 5. Diffusion equation linear diffusion equation eqworld.
Different stages of the example should be displayed, along with prompting messages in. Lou odette american international group aig october 17, 2006 1 nonlinear drift in the continuum limit the pdf. It is revealed that the diffusivity corresponds to the angular momentum operator in quantum theory. This work deals with this model for nuclear reactors with hexagonal. Equation 1 is known as a onedimensional diffusion equation, also often referred to as a heat equation. Nonlinear reaction diffusion equation with michaelismenten kinetics and adomian decomposition method.
It is not exact, but for most of this course it is the model that we will use to describe the behavior. The parameter \\alpha\ must be given and is referred to as the diffusion coefficient. This notebook is an entirely selfcontained solution to a basic neutron diffision equation for a reactor rx made up of a single fuel rod. Twogroup diffusion theory and the approximate representation of. 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. Neutron diffusion 90 if we insert the diffusion approximation 23 into our balance equation 4, we obtain. To show the validity of these diffusion boundary conditions, we solve a sample problem. 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. 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. 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. 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. Solution of the transport equations using a moving coordinate system ole krogh jensen and bruce a.
The convectivediffusion equation is the governing equation of many important transport phenomena in building physics. Numerical solution of the neutron diffusion equation has been done by many numerical methods such as the finite difference, finite element and boundary element. These equations are based ontheconceptoflocal neutron balance, which takes int 11. 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 discuss the diffusion equation as an asymptotic limit of the neutron transport equation for large scattering cross sections. A comparison of some numerical methods for the advection. 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.
This is the energydependent neutron diffusion equation. 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 section was about conditions for a stable, selfsustained fission chain reaction and how to maintain such conditions. The derivation of diffusion equation is based on ficks law which is derived under many assumptions. Derivation of the neutron diffusion equation nuclear. Derivation of the neutron diffusion equation introduction. Chapter 7 the diffusion equation the diffusionequation is a partial differentialequationwhich describes density. Pdf solution of fixed source neutron diffusion equation. The neutron diffusion equation can be solved analytically in academic cases or using standard numerical analysis techniques such as the. In previous section we dealt with the multiplication system and we defined the infinite and finite multiplication factor.
The diffusion equation can, therefore, not be exact or valid at places with strongly differing diffusion coefficients or in strongly absorbing media. Usa received 4 march 1979 a convectiondiffusion equation arises from the conservation equations in miscible and. 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. Pdf numerical techniques for the neutron diffusion equations in. Approximation of the neutron diffusion equation on. We now add a convection term \ \boldsymbolv\cdot\nabla u \ to the diffusion equation to obtain the wellknown convectiondiffusion equation.
Fundamental concepts and language diffusion mechanisms vacancy diffusion interstitial diffusion impurities. Chapter 2 the diffusion equation and the steady state. The derivation of the diffusion equation will depend on ficks law, even though a direct derivation from the transport equation is also possible. 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. Need to solve the transport equation for an accurate estimate of local reaction rates, particularly near fuel rods and strong neutron. Before attempting to solve the equation, it is useful to understand how the analytical. Pdf this paper presents a general theoretical analysis of the neutron motion problem in a nuclear reactor, where large variations on neutron. Understand origin, limitations of neutron diffusion from. Solution of onegroup neutron diffusion equation for. Equation 19 is a nonhomogeneous ordinary differential equation that can be solved by the application of classical methods.
This problem contains no information about the spatial distribution of neutrons, because it is a point geometry problem. Then we shall discuss the common boundary conditions that are used to solve problems. Pdf the spacetime neutron diffusion equations with multigroup of delayed neutrons are a couple of the stiff nonlinear partial differential. It is commonly used to determine the behavior of nuclear reactor cores and experimental or industrial neutron beams. Unstructured grids and the multigroup neutron diffusion. This implies that the diffusion theory may show deviations from a more accurate solution of the transport equation in the proximity of external. The general solution is composed by sum of the general integral of the associated homogeneous equation and the particular solution. 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.
Solution of the transport equations using a moving. The convectiondiffusion equation convectiondiffusion without a force term. The results show that solutions of the diffusion equation with these well. The nal piece of the puzzle requires the use of an empirical physical principle of heat ow. Secondorder parabolic partial differential equations nonhomogeneous heat diffusion equation 1. Neutron transport is the study of the motions and interactions of neutrons with materials. Application of the finite element method to the three. The helmholtz equation is derived, and the limitations on diffusion equation as well as the boundary conditions used in its application to realistic.
1339 1034 353 773 276 335 993 753 118 1330 1303 1454 1256 859 1303 1214 465 843 1200 1260 1235 399 1105 593 1360 449 391 1426 725 767 129 1027 92 1318 924 1475 290 998 595