|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| '''Neutron transport''' is the study of the motions and interactions of [[neutron]]s with materials. Nuclear scientists and [[nuclear engineer|engineers]] often need to know where neutrons are in an apparatus, what direction they are going, and how quickly they are moving. It is commonly used to determine the behavior of [[nuclear reactor]] cores and experimental or industrial neutron [[Particle beam|beams]]. Neutron transport is a type of [[radiative transport]].
| | For your very own offense, you might contain Gunboats which can simply shoot at enemy defense coming from a long range and Landing Products which you must fill up when you train equipment for example Rifleman, Heavy, Zooka, Warrior and Takes a dive. To your village defenses, you might offer structures like Mortar, Hardware Gun, Sniper Tower, Cannon, Flamethrower, Mine, Tank Mine, Boom Cannon and Bomb Launcher to assist a person eradicate enemies.<br><br>Game titles are fun to compete against your kids. Aid you learn much more info on your kid's interests. Sharing interests with your kids like this can likewise create great conversations. It also gives an opportunity to monitor expansion of their skills.<br><br>Nevertheless, if you want avoid at the top of the competitors, there are a few simple points you want to keep in mind. Realize your foe, recognize the game and the win will be yours. It is possible to look at the aid of [http://prometeu.net clash of clans hack] tools and the other rights if you just like your course. Terribly for your convenience, beneath are the general details in this particular sport that you should try to remember of. Read all of them meticulously!<br><br>Our own acceptable abatement for not enough best stretches of green energy is essential. Without them prices would bound grow into prohibitive and cipher would likely purchase them.<br><br>We each can use this route to acquisition the discount of any time in the course of 1hr and one day. For archetype to selection the majority of send up 4 a long time, acting x = 15, 400 abnormal and / or you receive y equals 51 gems.<br><br>Need different [http://www.reddit.com/r/howto/search?q=question question] the extent which always it''s a 'strategy'" games. A good moron without strategy in nearly any respect will advance amongst most of the gamers over time. So long as the individual sign in occasionally to be sure your primary 'builders'" are building something, your game power will [http://dict.leo.org/?search=increase increase]. That''s every bit of there's going without going for walks shoes. Individuals which have been the most effective every single in the game are, typically, those who may very well be actually playing a long, plus those who handed over real cash to identify extra builders. (Applying two builders, an special one can possibly can also be obtained for 700 gems which cost $4.99 and the next i costs 1000 gems.) When it comes to four builders, you will advance amongst people about doubly as fast just like a guy with two builders.<br><br>You don''t necessarily need one of the reliable troops to win wins. A mass volume of barbarians, your first-level troop, could well totally destroy an assailant village, and strangely it''s quite enjoyable to examine the virtual carnage. |
| | |
| ==Background==
| |
| | |
| Neutron transport has roots in the [[Boltzmann equation]], which was used in the 1800s to study the kinetic theory of gases. It did not receive large-scale development until the invention of chain-reacting nuclear reactors in the 1940s. As neutron distributions came under detailed scrutiny, elegant approximations and analytic solutions were found in simple geometries. However, as computational power has increased, numerical approaches to neutron transport have become prevalent. Today, with massively parallel computers, neutron transport is still under very active development in academia and research institutions throughout the world.
| |
| | |
| ==Neutron Transport Equation==
| |
| | |
| The neutron transport equation is a balance statement that conserves neutrons. Each term represents a gain or a loss of a neutron, and the balance, in essence, claims that neutrons gained equals neutrons lost. It is formulated as follows:<ref name=Adams>{{cite book |last=Adams |first=Marvin L. |title=Introduction to Nuclear Reactor Theory |year=2009 |publisher=Texas A&M University}}</ref>
| |
| :<math>\frac{1}{v(E)}\frac{\partial\psi(\mathbf{r},E,\mathbf{\hat{\Omega}},t)}{\partial t}+\mathbf{\hat{\Omega}}\cdot\nabla\psi(\mathbf{r},E,\mathbf{\hat{\Omega}},t)+\Sigma_t(\mathbf{r},E,t)\,\psi(\mathbf{r},E,\mathbf{\hat{\Omega}},t)=\quad</math>
| |
| | |
| :<math>\quad\frac{\chi_p \left( E \right)}{4\pi}\int_0^{\infty} dE^{\prime}\nu_p \left( E^{\prime} \right) \Sigma_f \left(\mathbf{r}, E^{\prime}, t \right) \phi \left( \mathbf{r}, E^{\prime}, t \right) + \sum_{i=1}^N \frac{\chi_{di}\left( E \right)}{4\pi} \lambda_i C_i \left( \mathbf{r}, t \right)+\quad</math>
| |
| | |
| :<math>\quad\int_{4\pi}d\Omega^\prime\int^{\infty}_{0}dE^\prime\,\Sigma_s(\mathbf{r},E^\prime\rightarrow E,\mathbf{\hat{\Omega}}^\prime\rightarrow \mathbf{\hat{\Omega}},t)\psi(\mathbf{r},E^\prime,\mathbf{\hat{\Omega}^\prime},t)+s(\mathbf{r},E,\mathbf{\hat{\Omega}},t)</math>
| |
| | |
| Where:
| |
| {| class="wikitable" border="1"
| |
| |-
| |
| ! Symbol
| |
| ! Meaning
| |
| ! Comments
| |
| |-
| |
| | <math>\mathbf{r}</math>
| |
| | Position vector (i.e. x,y,z)
| |
| |
| |
| |-
| |
| |<math>E</math>
| |
| |Energy
| |
| |
| |
| |-
| |
| |<math>\mathbf{\hat{\Omega}}=\frac{\mathbf{v}(E)}{|\mathbf{v}(E)|}=\frac{\mathbf{v}(E)}{{v(E)}}</math>
| |
| |Unit vector ([[solid angle]]) in direction of motion
| |
| |
| |
| |-
| |
| |<math>t</math>
| |
| |Time
| |
| |
| |
| |-
| |
| |<math> \mathbf{v}(E) </math>
| |
| |Neutron velocity vector
| |
| |
| |
| |-
| |
| |<math>\psi(\mathbf{r},E,\mathbf{\hat{\Omega}},t)dr\,dE\,d\Omega</math>
| |
| |Angular neutron flux <br /> Amount of neutron track length in a differential volume <math>dr</math> about <math>r</math>, associated with particles of a differential energy in <math>dE</math> about <math>E</math>, moving in a differential solid angle in <math>d\Omega</math> about <math>\mathbf{\hat{\Omega}}</math>, at time <math>t</math>.
| |
| |Note integrating over all angles yields [[Neutron flux|scalar neutron flux]]<br /> <math>\phi \ = \ \int_{4\pi}d\Omega\varphi</math>
| |
| |-
| |
| |<math>\phi(\mathbf{r},E,t)dr\,dE</math>
| |
| |[[Neutron flux|Scalar neutron flux]] <br /> Amount of neutron track length in a differential volume <math>dr</math> about <math>r</math>, associated with particles of a differential energy in <math>dE</math> about <math>E</math>, at time <math>t</math>.
| |
| |
| |
| |-
| |
| |<math>\nu_p</math>
| |
| |Average number of neutrons produced per fission (e.g., 2.43 for U-235).<ref>{{cite web|title=ENDF Libraries|url=https://www-nds.iaea.org/exfor/endf.htm}}</ref>
| |
| |
| |
| |-
| |
| |<math>\chi_p(E)</math>
| |
| |Probability density function for neutrons of exit energy <math>E</math> from all neutrons produced by fission
| |
| |
| |
| |-
| |
| |<math>\chi_{di}(E)</math>
| |
| |Probability density function for neutrons of exit energy <math>E</math> from all neutrons produced by delayed neutron precursors
| |
| |
| |
| |-
| |
| |<math>\Sigma_t(\mathbf{r},E,t)</math>
| |
| |Macroscopic total [[Nuclear cross section|cross section]], which includes all possible interactions
| |
| |
| |
| |-
| |
| |<math>\Sigma_f(\mathbf{r},E^{\prime},t)</math>
| |
| |Macroscopic fission [[Nuclear cross section|cross section]], which includes all fission interactions in <math>dE^{\prime}</math> about <math>E^{\prime}</math>
| |
| |
| |
| |-
| |
| |<math>\Sigma_s(\mathbf{r},E'\rightarrow E,\mathbf{\hat{\Omega}}'\rightarrow \mathbf{\hat{\Omega}},t)dE^\prime d\Omega^\prime</math>
| |
| |Double differential scattering cross section<br />Characterizes scattering of a neutron from an incident energy <math>E</math> and direction <math>\mathbf{\hat{\Omega}}</math> to a final energy <math>E^\prime</math> in <math>dE^\prime</math> and direction <math>\mathbf{\hat{\Omega}}^\prime</math> in <math>d\Omega^\prime</math>.
| |
| |
| |
| |-
| |
| |<math>N</math>
| |
| |Number of delayed neutron precursors
| |
| |
| |
| |-
| |
| |<math>\lambda_i</math>
| |
| |Decay constant for precursor i
| |
| |
| |
| |-
| |
| |<math>C_i \left( \mathbf{r}, t \right)</math>
| |
| |Total number of precursor i in <math>\mathbf{r}</math> at time <math>t</math>
| |
| |
| |
| |-
| |
| |<math>s(\mathbf{r},E,\mathbf{\hat{\Omega}},t)</math>
| |
| |Source term
| |
| |
| |
| |-
| |
| |}
| |
| | |
| The transport equation can be applied to a given part of phase space (time t, energy E, location <math>\mathbf{r}</math>, and direction of travel <math>\mathbf{\hat{\Omega}}</math>). The first term represents the time rate of change of neutrons in the system. The second terms describes the movement of neutrons into or out of the volume of space of interest. The third term accounts for all neutrons that have a collision in that phase space. The first term on the right hand side is the production of neutrons in this phase space due to fission, while the second term on the right hand side is the production of neutrons in this phase space due to delayed neutron precursors (i.e., unstable nuclei which undergo neutron decay). The third term on the right hand side is in-scattering, these are neutrons that enter this area of phase space as a result of scattering interactions in another. The fourth term on the right is a generic source. The equation is usually solved to find <math>\phi(\mathbf{r},E)</math>, since that will allow for the calculation of reaction rates, which are of primary interest in shielding and dosimetry studies.
| |
| | |
| ==Types of neutron transport calculations==
| |
| Several basic types of neutron transport problems exist, depending on the type of problem being solved.
| |
| | |
| ===Fixed Source===
| |
| A fixed source calculation involves imposing a known neutron source on a medium and determining the resulting neutron distribution throughout the problem. This type of problem is particularly useful for shielding calculations, where a designer would like to minimize the neutron dose outside of a shield while using the least amount of shielding material. For instance, a spent nuclear fuel cask requires shielding calculations to determine how much concrete and steel is needed to safely protect the truck driver who is shipping it.
| |
| | |
| ===Criticality===
| |
| [[Nuclear fission|Fission]] is the process through which a nucleus splits into (typically two) smaller atoms. If fission is occurring, it is often of interest to know the asymptotic behavior of the system. A reactor is called “critical” if the chain reaction is self-sustaining and time-independent. If the system is not in equilibrium the asymptotic neutron distribution, or the fundamental mode, will grow or decay exponentially over time.
| |
| | |
| Criticality calculations are used to analyze steady-state multiplying media (multiplying media can undergo fission), such as a critical nuclear reactor. The loss terms (absorption, out-scattering, and leakage) and the source terms (in-scatter and fission) are proportional to the neutron flux, contrasting with fixed-source problems where the source is independent of the flux. In these calculations, the presumption of time invariance requires that neutron production exactly equals neutron loss.
| |
| | |
| Since this criticality can only be achieved by very fine manipulations of the geometry (typically via control rods in a reactor), it is unlikely that the modeled geometry will be truly critical. To allow some flexibility in the way models are set up, these problems are formulated as eigenvalue problems, where one parameter is artificially modified until criticality is reached. The most common formulations are the time-absorption and the multiplication eigenvalues, also known as the alpha and k eigenvalues. The alpha and k are the tunable quanitites.
| |
| | |
| K-eigenvalue problems are the most common in nuclear reactor analysis. The number of neutrons produced per fission is multiplicatively modified by the dominant eigenvalue. The resulting value of this eigenvalue reflects the time dependence of the neutron density in a multiplying medium.
| |
| *''k<sub>eff</sub>'' < 1, subcritical: the neutron density is decreasing as time passes;
| |
| *''k<sub>eff</sub>'' = 1, critical: the neutron density remains unchanged; and
| |
| *''k<sub>eff</sub>'' > 1, supercritical: the neutron density is increasing with time.
| |
| | |
| In the case of a [[nuclear reactor]], neutron flux and power density are proportional, hence during reactor start-up ''k<sub>eff</sub>'' > 1, during reactor operation ''k<sub>eff</sub>'' = 1 and ''k<sub>eff</sub>'' < 1 at reactor shutdown.
| |
| | |
| ==Computational Methods==
| |
| Both fixed-source and criticality calculations can be solved using [[Deterministic system (mathematics)|deterministic methods]] or [[stochastic process|stochastic methods]]. In deterministic methods the transport equation (or an approximation of it, such as [[Diffusion process|diffusion theory]]) is solved as a differential equation. In stochastic methods such as [[Monte Carlo method|Monte Carlo]] discrete particle histories are tracked and averaged in a random walk directed by measured interaction probabilities. Deterministic methods usually involve multi-group approaches while Monte Carlo can work with multi-group and continuous energy cross-section libraries. Multi-group calculations are usually iterative, because the group constants are calculated using flux-energy profiles, which are determined as the result of the neutron transport calculation.
| |
| | |
| ===Discretization in Deterministic Methods===
| |
| To numerically solve the transport equation using algebraic equations on a computer, the spatial, angular, energy, and time variables must be [[Discretization|discretized]].
| |
| * Spatial variables are typically discretized by simply breaking the geometry into many small regions on a mesh. The balance can then be solved at each mesh point using [[finite difference]] or by nodal methods.
| |
| * Angular variables can be discretized by discrete ordinates and weighting [[Numerical integration|quadrature sets]] (giving rise to the S<sub>N</sub> methods), or by functional expansion methods with the [[spherical harmonics]] (leading to the P<sub>N</sub> methods).
| |
| * Energy variables are typically discretized by the multi-group method, where each energy group represents one constant energy. As few as 2 groups can be sufficient for some [[Thermal-neutron reactor|thermal reactor]] problems, but [[Fast-neutron reactor|fast reactor]] calculations may require many more.
| |
| * The time variable is broken into discrete time steps, with time derivatives replaced with difference formulas.
| |
| | |
| ===Computer Codes Used In Neutron Transport===
| |
| | |
| ====Probabilistic codes====
| |
| *'''OpenMC''' - An [[Massachusetts Institute of Technology|MIT]] developed open source Monte Carlo code <ref>{{cite web|title=OpenMC|url=http://mit-crpg.github.io/openmc/}}</ref>
| |
| *'''[[MCNP]]''' - A [[Los Alamos National Laboratory|LANL]] developed Monte Carlo code for general radiation transport
| |
| *'''KENO''' - An [[Oak Ridge National Laboratory|ORNL]] developed Monte Carlo code for criticality analysis
| |
| *'''MCBEND''' - An ANSWERS Software Service developed Monte Carlo code for general radiation transport
| |
| *'''Serpent''' - A [[VTT Technical Research Centre of Finland|Finnish]] developed Monte Carlo neutron transport code <ref>{{cite web|title=PSG2 Serpent|url=http://montecarlo.vtt.fi/}}</ref>
| |
| | |
| ====Deterministic codes====
| |
| *'''Attila''' - A commercial transport code
| |
| *'''DRAGON''' - An open-source lattice physics code
| |
| *'''PHOENIX/ANC''' - A proprietary lattice-physics and global diffusion code suite from [[Westinghouse Electric Corporation|Westinghouse Electric]]
| |
| *'''PARTISN''' - A [[Los Alamos National Laboratory|LANL]] developed transport code based on the discrete ordinates method
| |
| *'''NEWT''' - An [[Oak Ridge National Laboratory|ORNL]] developed 2-D S<sub>N</sub> code
| |
| *'''DIF3D/VARIANT''' - An Argonne National Laboratory developed 3-D code originally developed for fast reactors
| |
| *'''DENOVO''' - A massively parallel transport code under development by [[Oak Ridge National Laboratory|ORNL]]
| |
| *'''DANTSYS'''
| |
| *'''RAPTOR-M3G''' - A proprietary parallel radiation transport code developed by [[Westinghouse Electric Company]]
| |
| *'''MPACT''' - A parallel 3D [[method of characteristics]] code under development by the [[University of Michigan]]
| |
| | |
| ==See also==
| |
| * [[Nuclear Reactor]]
| |
| * [[Boltzmann equation]]
| |
| * [[Time-dependent Neutronics and Temperatures|TINTE]]
| |
| * [[Neutron scattering]]
| |
| | |
| ==References==
| |
| | |
| {{reflist}}
| |
| * Lewis, E., & Miller, W. (1993). Computational Methods of Neutron Transport. American Nuclear Society. ISBN 0-89448-452-4.
| |
| * Duderstadt, J., & Hamilton, L. (1976). Nuclear Reactor Analysis. New York: Wiley. ISBN 0-471-22363-8.
| |
| * Marchuk, G. I., & V. I. Lebedev (1986). Numerical Methods in the Theory of Neutron Transport. Taylor & Francis. p. 123. ISBN 978-3-7186-0182-0.
| |
| | |
| ==External links==
| |
| * [http://mcnp.lanl.gov/ LANL MCNP6 website]
| |
| * [http://mcnpx.lanl.gov/ LANL MCNPX website]
| |
| * [http://montecarlo.vtt.fi/ VTT Serpent website]
| |
| | |
| [[Category:Neutron|Transport]]
| |
| [[Category:Nuclear physics]]
| |
For your very own offense, you might contain Gunboats which can simply shoot at enemy defense coming from a long range and Landing Products which you must fill up when you train equipment for example Rifleman, Heavy, Zooka, Warrior and Takes a dive. To your village defenses, you might offer structures like Mortar, Hardware Gun, Sniper Tower, Cannon, Flamethrower, Mine, Tank Mine, Boom Cannon and Bomb Launcher to assist a person eradicate enemies.
Game titles are fun to compete against your kids. Aid you learn much more info on your kid's interests. Sharing interests with your kids like this can likewise create great conversations. It also gives an opportunity to monitor expansion of their skills.
Nevertheless, if you want avoid at the top of the competitors, there are a few simple points you want to keep in mind. Realize your foe, recognize the game and the win will be yours. It is possible to look at the aid of clash of clans hack tools and the other rights if you just like your course. Terribly for your convenience, beneath are the general details in this particular sport that you should try to remember of. Read all of them meticulously!
Our own acceptable abatement for not enough best stretches of green energy is essential. Without them prices would bound grow into prohibitive and cipher would likely purchase them.
We each can use this route to acquisition the discount of any time in the course of 1hr and one day. For archetype to selection the majority of send up 4 a long time, acting x = 15, 400 abnormal and / or you receive y equals 51 gems.
Need different question the extent which always its a 'strategy'" games. A good moron without strategy in nearly any respect will advance amongst most of the gamers over time. So long as the individual sign in occasionally to be sure your primary 'builders'" are building something, your game power will increase. Thats every bit of there's going without going for walks shoes. Individuals which have been the most effective every single in the game are, typically, those who may very well be actually playing a long, plus those who handed over real cash to identify extra builders. (Applying two builders, an special one can possibly can also be obtained for 700 gems which cost $4.99 and the next i costs 1000 gems.) When it comes to four builders, you will advance amongst people about doubly as fast just like a guy with two builders.
You dont necessarily need one of the reliable troops to win wins. A mass volume of barbarians, your first-level troop, could well totally destroy an assailant village, and strangely its quite enjoyable to examine the virtual carnage.