Difference between revisions of "Doppler-broadening preprocessor routine"
(→Theory) |
(→Pitfalls) |
||
(6 intermediate revisions by 3 users not shown) | |||
Line 1: | Line 1: | ||
== Introduction == | == Introduction == | ||
− | Traditionally, the effect of thermal motion on reaction rates has been taken into account by generating temperature-dependent, ''effective'' cross section libraries using dedicated processing codes such as NJOY, and using these cross sections in reactor physical calculations. This approach works well as long as there is no need to model the temperatures accurately or if the temperature of a system is very common, for instance room temperature. However, when modelling reactor conditions, the temperatures often need to be modeled at a high accuracy and, in addition, the temperatures are strongly case-dependent. Because of this, the cross section libraries often need to be tailored for each calculation case separately. | + | Traditionally, the effect of thermal motion on reaction rates has been taken into account by generating temperature-dependent, ''effective'' cross section libraries using dedicated processing codes such as NJOY <ref name="njoy">R. E. MacFarlane and D. W. Muir, H. "''NJOY99.0 Code System for Producing Pointwise and Multigroup Neutron and Photon Cross Sections from {ENDF/B} Data''", Los Alamos National Laboratory, PSR-480, (2000)</ref>, and using these cross sections in reactor physical calculations. This approach works well as long as there is no need to model the temperatures accurately or if the temperature of a system is very common, for instance room temperature. However, when modelling reactor conditions, the temperatures often need to be modeled at a high accuracy and, in addition, the temperatures are strongly case-dependent. Because of this, the cross section libraries often need to be tailored for each calculation case separately. |
To simplify the calculation chain, a Doppler-broadening preprocessor routine was implemented in Serpent 2 (and Serpent 1). The routine is able to increase the temperature of a free-atom cross section below the energy region of unresolved resonances. The Doppler-broadening is performed analytically based on the Solbrig's kernel, and the broadening is relatively fast. | To simplify the calculation chain, a Doppler-broadening preprocessor routine was implemented in Serpent 2 (and Serpent 1). The routine is able to increase the temperature of a free-atom cross section below the energy region of unresolved resonances. The Doppler-broadening is performed analytically based on the Solbrig's kernel, and the broadening is relatively fast. | ||
Line 7: | Line 7: | ||
== Use == | == Use == | ||
− | '''mat''' ''NAME'' ''DENS'' | + | '''mat''' ''NAME'' ''DENS'' ''T'' [ ''NUC<sub>1</sub>'' ''FRAC<sub>1</sub>'' ''NUC<sub>2</sub>'' ''FRAC<sub>2</sub>'' ... ] |
Sets the temperature of the material, to be used in both Doppler-broadening preprococessor and free gas treatment of elastic scattering. | Sets the temperature of the material, to be used in both Doppler-broadening preprococessor and free gas treatment of elastic scattering. | ||
Line 32: | Line 32: | ||
== Theory == | == Theory == | ||
− | <ref name="viitanen_msc">Viitanen, T. | + | The effect of thermal motion on reaction rates can be described by using so called effective cross sections in the transport calculation. The effective cross section must reproduce the true reaction rates, which can be calculated as an integral over the thermal motion of a target. The condition of constant reaction rates can be written as |
+ | |||
+ | <math> v \sigma_{\mathrm{eff}}(v,T,A) = \int |\mathbf{v'}| \sigma(|\mathbf{v'}|) P(\mathbf{V_{\mathrm{t}}},T, A) d\mathbf{v'} \Leftrightarrow </math> | ||
+ | |||
+ | <math> \sigma_{\mathrm{eff}}(v,T,A) = \frac{1}{v} \int |\mathbf{v'}| \sigma(|\mathbf{v'}|) P(\mathbf{\mathbf{V_{\mathrm{t}}}},T, A) d\mathbf{v'}, </math> | ||
+ | |||
+ | where <math>v</math> is the velocity of the neutron, <math>\sigma_{eff}</math> is the effective cross section, <math>\mathbf{v'}</math> is the relative velocity of the neutron to the target and <math>P(\mathbf{\mathbf{V_{\mathrm{t}}}},T, A)</math> is the probability density function of the target motion <math>\mathbf{\mathbf{V_{\mathrm{t}}}}</math>, which depends on temperature <math>T</math> and the atomic weight ratio of nuclide <math>A</math>. | ||
+ | |||
+ | If the target motion is assumed to obey the Maxwell-Boltzmann distribution, the equation can be developed into a widely-used Doppler-broadening formula, a.k.a. Solbrig's kernel <ref name="solbrig">A. W. Solbrig, Jr. ''"Doppler Effect in Neutron Absorption Resonances."'', Am. J. Phys., 29, 257-261, (1960).</ref> | ||
+ | |||
+ | <math> \sigma_{\mathrm{eff}}(v,T,A) = \frac{\gamma}{v^2 \sqrt{\pi}} \int_{0}^{\infty} v'^2 \sigma(v') (e^{-\gamma^2(v-v')^2} - e^{-\gamma^2(v+v')^2}) dv' </math> | ||
+ | |||
+ | <math> \gamma(T,A_{n}) = \sqrt{\frac{A M}{2 k_{\mathrm{B}}T}} </math> | ||
+ | |||
+ | in which the angular dependencies have been removed via integration over all angles. In these equations, <math>M</math> is the neutron mass and <math>k_{\mathrm{B}}</math> is the Boltzmann constant. The integral in the Doppler-broadening equation can be calculated analytically if the cross section is piece-wise linear, which is the case with ACE data. It should be emphasized that the integration methodology differs from the most widely used SIGMA1 Doppler-broadening algorithm<ref name="sigma1"> Cullen, D.E., ''"Program SIGMA1 (version 79-1): Doppler broaden evaluated cross sections in the evaluated nuclear data file/version B (ENDF/B) format"'', Lawrence Livermore National Laboratory, UCRL-50400, (1979).</ref>. | ||
+ | |||
+ | |||
+ | More information on the integration methodology and different assumptions can be found in Reference <ref name="viitanen_msc">Viitanen, T. | ||
''"Implementing a Doppler-preprocessor of cross section libraries in reactor physics code Serpent."'', M.Sc. Thesis, Helsinki University of Technology, (2009).</ref> | ''"Implementing a Doppler-preprocessor of cross section libraries in reactor physics code Serpent."'', M.Sc. Thesis, Helsinki University of Technology, (2009).</ref> | ||
Line 44: | Line 61: | ||
** the resulting cross section requires slightly more memory than what a cross section processed with NJOY would require at the same temperature and | ** the resulting cross section requires slightly more memory than what a cross section processed with NJOY would require at the same temperature and | ||
** the reconstruction tolerance of the resulting cross section is unknown. Experience has shown that the errors remain small, at least if the nearest library temperature is used as the basis for broadening. Regardless, some differences may be seen in very detailed comparisons. | ** the reconstruction tolerance of the resulting cross section is unknown. Experience has shown that the errors remain small, at least if the nearest library temperature is used as the basis for broadening. Regardless, some differences may be seen in very detailed comparisons. | ||
+ | * The preprocessor does ''NOT'' apply Doppler-broadening correctly to the KERMA coefficients used in energy deposition modes 2/3 and neutron heating detector response (MT=-4). The temperature dependence of KERMA coefficients can be accounted for (accurately) only by using ACE libraries corresponding to the exact temperature of each material. These ACE libraries can be processed using NJOY. | ||
== References == | == References == | ||
<references/> | <references/> | ||
+ | |||
+ | [[Category:Input]] | ||
+ | [[Category:Theory]] |
Latest revision as of 21:03, 28 December 2021
Contents
Introduction
Traditionally, the effect of thermal motion on reaction rates has been taken into account by generating temperature-dependent, effective cross section libraries using dedicated processing codes such as NJOY ^{[1]}, and using these cross sections in reactor physical calculations. This approach works well as long as there is no need to model the temperatures accurately or if the temperature of a system is very common, for instance room temperature. However, when modelling reactor conditions, the temperatures often need to be modeled at a high accuracy and, in addition, the temperatures are strongly case-dependent. Because of this, the cross section libraries often need to be tailored for each calculation case separately.
To simplify the calculation chain, a Doppler-broadening preprocessor routine was implemented in Serpent 2 (and Serpent 1). The routine is able to increase the temperature of a free-atom cross section below the energy region of unresolved resonances. The Doppler-broadening is performed analytically based on the Solbrig's kernel, and the broadening is relatively fast.
Use
mat NAME DENS T [ NUC_{1} FRAC_{1} NUC_{2} FRAC_{2} ... ]
Sets the temperature of the material, to be used in both Doppler-broadening preprococessor and free gas treatment of elastic scattering.
NAME | : is the material name |
DENS | : density of the material, positive value for atomic density and negative value for mass density |
T | : temperature of the material in K |
NUC_{i} | : name (ZA + library identifier) of nuclide i |
FRAC_{i} | : fraction of nuclide i, positive value for atomic fraction and negative value for mass fraction |
Theory
The effect of thermal motion on reaction rates can be described by using so called effective cross sections in the transport calculation. The effective cross section must reproduce the true reaction rates, which can be calculated as an integral over the thermal motion of a target. The condition of constant reaction rates can be written as
where is the velocity of the neutron, is the effective cross section, is the relative velocity of the neutron to the target and is the probability density function of the target motion , which depends on temperature and the atomic weight ratio of nuclide .
If the target motion is assumed to obey the Maxwell-Boltzmann distribution, the equation can be developed into a widely-used Doppler-broadening formula, a.k.a. Solbrig's kernel ^{[2]}
in which the angular dependencies have been removed via integration over all angles. In these equations, is the neutron mass and is the Boltzmann constant. The integral in the Doppler-broadening equation can be calculated analytically if the cross section is piece-wise linear, which is the case with ACE data. It should be emphasized that the integration methodology differs from the most widely used SIGMA1 Doppler-broadening algorithm^{[3]}.
More information on the integration methodology and different assumptions can be found in Reference ^{[4]}
Pitfalls
- The preprocessor does NOT adjust the temperature of thermal scattering data. For moderator materials, the temperature of thermal scattering data must be changed separately by either
- changing the library defined in the therm card or
- interpolating the thermal scattering data to the new temperature using the therm card
- The preprocessor does NOT adjust the temperature of the unresolved region probability tables. When modeling fast systems, the temperature dependence of reaction rates can be modeled (accurately) only by using ACE libraries corresponding to the exact temperature of each material. These ACE libraries can be processed using NJOY.
- The preprocessor does NOT adjust the energy grid of the cross section in the temperature adjustment process. Consequently, the resulting energy grid is suboptimal such that
- the resulting cross section requires slightly more memory than what a cross section processed with NJOY would require at the same temperature and
- the reconstruction tolerance of the resulting cross section is unknown. Experience has shown that the errors remain small, at least if the nearest library temperature is used as the basis for broadening. Regardless, some differences may be seen in very detailed comparisons.
- The preprocessor does NOT apply Doppler-broadening correctly to the KERMA coefficients used in energy deposition modes 2/3 and neutron heating detector response (MT=-4). The temperature dependence of KERMA coefficients can be accounted for (accurately) only by using ACE libraries corresponding to the exact temperature of each material. These ACE libraries can be processed using NJOY.
References
- ^ R. E. MacFarlane and D. W. Muir, H. "NJOY99.0 Code System for Producing Pointwise and Multigroup Neutron and Photon Cross Sections from {ENDF/B} Data", Los Alamos National Laboratory, PSR-480, (2000)
- ^ A. W. Solbrig, Jr. "Doppler Effect in Neutron Absorption Resonances.", Am. J. Phys., 29, 257-261, (1960).
- ^ Cullen, D.E., "Program SIGMA1 (version 79-1): Doppler broaden evaluated cross sections in the evaluated nuclear data file/version B (ENDF/B) format", Lawrence Livermore National Laboratory, UCRL-50400, (1979).
- ^ Viitanen, T. "Implementing a Doppler-preprocessor of cross section libraries in reactor physics code Serpent.", M.Sc. Thesis, Helsinki University of Technology, (2009).