crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 21:55:41 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

do not add 0 K energy points to union grid

Browse source
This commit is contained in:
walshjon 2014-09-10 12:20:58 -07:00
parent 4d10704a87
commit 3a9f2268de
4 changed files with 7 additions and 11 deletions
Download patch file Download diff file Expand all files Collapse all files

6
docs/source/usersguide/input.rst
View file

@ -274,7 +274,7 @@ attributes or sub-elements:
of the nuclide to which a resonance scattering treatment is to be applied.
The ``method`` attribute gives the type of resonance scattering treatment
that is to be applied to the ``nuclide``. Acceptable inputs - none of
which are case-sensitive - for the ``method`` attribute are ``ARTS``,
which are case-sensitive - for the ``method`` attribute are ``ARES``,
``CXS``, ``WCM``, and ``DBRC``. Descriptions of each of these methods
are documented here_. The ``xs_label`` attribute gives the label for the
cross section data of the ``nuclide`` at a given temperature. The
@ -291,7 +291,7 @@ attributes or sub-elements:
<resonance_scattering>
<scatterer>
<nuclide>U-238</nuclide>
<method>ARTS</method>
<method>ARES</method>
<xs_label>92238.72c</xs_label>
<xs_label_0K>92238.00c</xs_label_0K>
<E_min>5.0e-6</E_min>
@ -317,7 +317,7 @@ attributes or sub-elements:
be specified, each corresponding to a single nuclide at a single
temperature.
*Defaults*: None (scatterer), ARTS (method), 0.01 eV (E_min), 1.0 keV (E_max)
*Defaults*: None (scatterer), ARES (method), 0.01 eV (E_min), 1.0 keV (E_max)
``<run_cmfd>`` Element
----------------------

4
src/ace_header.F90
View file

@ -159,14 +159,14 @@ module ace_header
end type Nuclide
!===============================================================================
! NUCLIDE_0K temporarily contains all 0K cross section data and other parameters
! NUCLIDE0K temporarily contains all 0K cross section data and other parameters
! needed to treat resonance scattering before transferring them to NUCLIDE
!===============================================================================
type Nuclide0K
character(10) :: nuclide ! name of nuclide, e.g. U-238
character(16) :: scheme = 'arts' ! target velocity sampling scheme
character(16) :: scheme = 'ares' ! target velocity sampling scheme
character(10) :: name ! name of nuclide, e.g. 92235.03c
character(10) :: name_0K ! name of 0K nuclide, e.g. 92235.00c
real(8) :: E_min = 0.01e-6 ! lower cutoff energy for res scattering

3
src/energy_grid.F90
View file

@ -28,9 +28,6 @@ contains
do i = 1, n_nuclides_total
nuc => nuclides(i)
call add_grid_points(list, nuc % energy)
! add 0K points to the grid
if (nuc % resonant) call add_grid_points(list, nuc % energy_0K)
end do
! Set size of unionized energy grid

5
src/physics.F90
View file

@ -2,8 +2,7 @@ module physics
use ace_header, only: Nuclide, Reaction, DistEnergy
use constants
use cross_section, only: elastic_xs_0K, find_energy_index, &
union_grid_index
use cross_section, only: elastic_xs_0K
use endf, only: reaction_name
use error, only: fatal_error, warning
use fission, only: nu_total, nu_delayed
@ -885,7 +884,7 @@ contains
if (.not. reject) exit
end do
case ('arts')
case ('ares')
E_red = sqrt((awr * E) / kT)
E_low = (((E_red - 4.0_8)**2) * kT) / awr
E_up = (((E_red + 4.0_8)**2) * kT) / awr