The steady state and the diffusion equation the neutron field basic field quantity in reactor physics is the neutron angular flux density distribution. Analytical solutions of the diffusion differential equation kit. Chapter 2 the diffusion equation and the steady state. Pdf solution of fixed source neutron diffusion equation. Solution of fixed source neutron diffusion equation via homotopy perturbation method conference paper pdf available january 2010 with 107 reads how we measure reads. In previous section we dealt with the multiplication system and we defined the infinite and finite multiplication factor. The hideous neutron transport equation has been reduced to a simple oneliner neutron diffusion equation. Unstructured grids and the multigroup neutron diffusion. 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. Solution of the multigroup neutron diffusion equations by the finite element method. Twogroup diffusion theory and the approximate representation of the diffusion equation using finite differences applied to a discrete spatial mesh are introduced. The convectiondiffusion partial differential equation pde solved is, where is the diffusion parameter, is the advection parameter also called the transport parameter, and is the convection parameter. Get a printable copy pdf file of the complete article 246k, or click on a page image below.
An iteration method has been implemented to solve a neutron transport equation in a multigroup diffusion approximation. Diffusion equation and stochastic processes ncbi nih. What links here related changes upload file special pages permanent link page. The diffusion equation the diffusion equation, like the wave equation, provides a way to analyse some important physical processes that require evaluation as a function of space and time. 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 since these low order models do not come close to. The functions plug and gaussian runs the case with \ix\ as a discontinuous plug or a smooth gaussian function, respectively. A generalized equa tion is presented describing diffusion. Nonlinear reaction diffusion equation with michaelismenten. We then go on to consider the applications of both di usion equations. 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. This partial differential equation is dissipative but not dispersive. Lecture 28 solution of heat equation via fourier transforms and convolution theorem relvant sections of text. Solving the convectiondiffusion equation in 1d using.
It is commonly used to determine the behavior of nuclear reactor cores and experimental or industrial neutron beams. These two routines are combined by a subroutiw crossace. Iterative schemes for the neutron diffusion equation. This section was about conditions for a stable, selfsustained fission chain reaction and how to maintain such conditions. We consider the laxwendroff scheme which is explicit, the cranknicolson scheme which is implicit, and a nonstandard finite difference scheme mickens 1991. You can specify using the initial conditions button. Solution of the transport equations using a moving coordinate. This is the energydependent neutron diffusion equation. The diffusion equation is a parabolic partial differential equation.
Steadystate diffusion when the concentration field is independent of time and d is independent of c, fick. Other posts in the series concentrate on derivative approximation, the cranknicolson implicit method and the tridiagonal matrix solverthomas algorithm. This video describes the neutron diffusion in nuclear reactors. Neutron transport is the study of the motions and interactions of neutrons with materials. This implies that the diffusion theory may show deviations from a more accurate solution of the transport equation in the proximity of external. Discretization and solution of convection diffusion problems. A comparison of some numerical methods for the advectiondi. 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 those for other physics phenomena, e. This post is part of a series of finite difference method articles. To run this example from the base fipy directory, type.
Nonlinear diffusion these notes summarize the way i present this material, for my bene. In section 3 a number of desirable properties of where b is again arbitrary, but not necessarily the same as this nonlinear smoothing process are presented. Partial differential equations, piecewise constant arguments, oscillations, and stability. Chapter 2 the diffusion equation and the steady state weshallnowstudy the equations which govern the neutron field in a reactor. 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. Analytical solution to the onedimensional advection. The diffusion equation can, therefore, not be exact or valid at places with strongly differing diffusion coefficients or in strongly absorbing media. Reactor physics tutorial classi cation of time problems classi cation of time problems timedependent neutron population i short time problems seconds tens of minutes i reactor conditions altered change in k i intermediate time problems hours 1 or 2 days. Whatelse canwesayaboutthe shape of the distribution of particles. Solution of the multigroup neutron diffusion equations by the finite element method misfeldt, i.
The derivation of diffusion equation is based on ficks law which is derived under many assumptions. 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. Different stages of the example should be displayed, along with prompting messages in. These equations are based ontheconceptoflocal neutron balance, which takes int university of washington, seattle, washington 98195. Geometric heat equation and nonlinear diffusion of shapes and. Neutron diffusion 90 if we insert the diffusion approximation 23 into our balance equation 4, we obtain.
Pdf solution of neutrontransport multigroup equations. Neutron diffusion equation 9 to integrate equation 3, we must take into account that it constitutes a system of stiff differential equations, mainly due to the elements of the diagonal. Development of a three dimensional neutron diffusion code. 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. In order to obtain that, we must then use the diffusion equation. Dif3ds nodal option solves the multigroup steadystate neutron diffusion and for cartesian geometry only transport equations in two and threedimensional hexagonal and cartesian geometries.
For elliptic diffusion equations with random coefficient and source term, the probability measure of the solution. Steadystate diffusion when the concentration field is independent of time and d is independent of c, fick 2c0 s second law is reduced to laplaces equation, for simple geometries, such as permeation through a thin membrane, laplaces equation can. Iterative solution algorithms krylov subspace methods splitting methods multigrid. But everything in here is said in more detail, and better, in weickerts paper. The derivation of the diffusion equation will depend on ficks law, even though a direct derivation from the transport equation is also possible. Chapter 8 the reactiondiffusion equations reactiondiffusion rd equations arise naturally in systems consisting of many interacting components, e. These equations are based ontheconceptoflocal neutron balance, which takes int neutron diffusion theory. Solving the convectiondiffusion equation in 1d using finite. Three numerical methods have been used to solve the onedimensional advectiondiffusion equation with constant coefficients. Sensitivity of uncertainty propagation for the elliptic diffusion equation. The nuclear group constants for the diffusion equation are expressed as functions of the thermal hydraulic condition at each block in the core. Usa received 4 march 1979 a convectiondiffusion equation arises from the conservation equations in miscible and.
It is one of the computer codes maintained or developed by the nuclear engineering division. Discretization strategies finite element methods inadequacy of galerkin methods stabilization. Solving the diffusion equation explicitly quantstart. References references are in electronic format in file c784. Heat or diffusion equation in 1d derivation of the 1d heat equation separation of variables refresher worked examples kreysig, 8th edn, sections 11. The characterization of reactionconvectiondiffusion processes. Numerical solution of the 1d advectiondiffusion equation. Physical assumptions we consider temperature in a long thin wire of constant cross section and homogeneous material. The helmholtz equation is derived, and the limitations on diffusion equation as well as the boundary conditions used in its application to. A thermoelectric generator containing plutonium dioxide, used as a source of.
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. Feb 15, 2017 this video describes the neutron diffusion in nuclear reactors. Experiments with these two functions reveal some important observations. Finite difference method for solving neutron diffusion equation in hexagonal geometry conference paper pdf available september 2009 with 1,159 reads how we measure reads. This work introduces the alternatives that unstructured grids can provide. A code to solve one, two, and threedimensional finitedifference diffusion theory problems, anl8264, argonne national laboratory, argonne, il 1984. Hence, for its integration, it is convenient to use an implicit backward difference formula bdf 7. A comparison of some numerical methods for the advection. The diffusion equation is a partial differential equation which describes density fluc tuations in a material undergoing diffusion.
If the two coefficients and are constants then they are referred to as solute dispersion coefficient and uniform velocity, respectively, and the above equation reduces to equation 1. Diffusion equation linear diffusion equation eqworld. Introduction we have seen that the transport equation is exact, but difficult to solve. Unstructured grids and the multigroup neutron diffusion equation. 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 convection diffusion partial differential equation pde solved is, where is the diffusion parameter, is the advection parameter also called the transport parameter, and is the convection parameter. Chapter 8 the reaction diffusion equations reaction diffusion rd equations arise naturally in systems consisting of many interacting components, e. A detailed account of the construction of similarity variables to obtain selfsimilar solutions are shown for both the porous medium equation, for all n, and the thin film equation, in the case of n 1, with given initial conditions. The image data is written every so many 10 iterations through the time evolution of the temperature field, thereby allowing the production of a movie postmortem. Setting source to 0 solves the diffusion equation with no source. The neutron diffusion equation is often used to perform corelevel neutronic calculations.
Equation is known as a onedimensional diffusion equation, also often referred to as a heat equation. Solution of the transport equations using a moving coordinate system ole krogh jensen and bruce a. It is not exact, but for most of this course it is the model that we will use to describe the behavior. Finlayson department of chemical engineering, university of washington, seattle, washington 98195. Tofind out, we haveto workouttheprobabilities that theparticles step different distancesto therightortotheleft. This article provides some background to the mathematics of the neutron diffusion process and critical mass. The convectiondiffusion equation introduction and examples 2. Heat or diffusion equation in 1d university of oxford. Full text is available as a scanned copy of the original print version. The convection diffusion equation introduction and examples 2. Solution of the multigroup neutron diffusion equations by.
1085 380 352 859 1333 466 435 1300 1220 23 1196 1402 1417 667 1183 1194 1596 1306 553 1565 1183 1230 373 247 1435 1063 332 1295 978 648 1137 660 1426 1336 108 402