Difference between revisions of "Kraken tutorial 2024"

title "LDR-Lite"

global geometry square neutronics diffusion

library type polynomial path "../inputs/xs/narrow_2g.xs"

core width_xy      21.5036
     rows          9
     columns       9
     axial_nodes   22
     axial_heights
      15.000000                ! Support plate
       1.748000                ! Bottom zirconium plugs
       5.715000                ! Grid 1
  8 *  5.535625                ! Fuel
       5.715000                ! Grid 2
  8 *  5.535625                ! Fuel
       4.007000                ! Top plenum
      14.712000                ! Structures above active core
     division_x      2
     division_y      2
     radial_boundary_coeff 0.031477  ! Extrapolation length
     axial_boundary_coeff  0.031477  ! Extrapolation length

include "../inputs/ants_includes/fuel_assemblies.inc"
include "../inputs/ants_includes/radial_reflector.inc"
include "../inputs/ants_includes/control_rods.inc"
include "../inputs/ants_includes/core.inc"


iteration power                    1.0
          max_outer                10000
          control_variable         keff
          xenon_state              zero
          samarium_state           zero

# ---

simulationControl
{
    nPhases            2
    nIterMax           200
    urfU               0.125
    urfP               0.125
    urfH               0.25
    gravity            -9.81
    readRestart        0
    standAlone         0
    writeFormat        ascii
    materialProperties libTable
    thermalBC          power
}

# ---

inlet
{
    type               massFlowInlet
    massFlowRate       6.36    0.0   
    pressure           6.10e5         
    temperature        367.6 367.6   
    volumeFraction     1.0     0.0
}

outlet
{
    pressure           6.20e5
    temperature        380 380
    volumeFraction     1.0 0.0
}

fuelType rr
{
    nCells             20

    ! Align bottom of active core to z = 0.0

    bottomElevation    -0.05715

    cellHeight
       0.05715000                ! Grid 1
       0.05715000                ! Grid 1
  8 *  0.05535625                ! Fuel
       0.05715000                ! Grid 2
  8 *  0.05535625                ! Fuel
       0.05535625

    hydraulicDiameter   20 * 1.153052E-02

    thermalDiameter     20 * 1.278587E-02

    porosity            20 * 5.446598E-01

}

fuelType f1
{
    nCells             20

    ! Align bottom of active core to z = 0.0

    bottomElevation    -0.05715

    cellHeight
       0.05715000                ! Grid 1
       0.05715000                ! Grid 1
  8 *  0.05535625                ! Fuel
       0.05715000                ! Grid 2
  8 *  0.05535625                ! Fuel
       0.05535625

    hydraulicDiameter   20 * 1.153052E-02

    thermalDiameter     20 * 1.278587E-02

    porosity            20 * 5.446598E-01

}

core
{
    coreLattice        square
    assemblyPitch      0.2150360
}

!
!

loadPattern
{
  0   0  rr  rr  rr  rr  rr   0  0
  0  rr  rr  f1  f1  f1  rr  rr  0
 rr  rr  f1  f1  f1  f1  f1  rr rr
 rr  f1  f1  f1  f1  f1  f1  f1 rr
 rr  f1  f1  f1  f1  f1  f1  f1 rr
 rr  f1  f1  f1  f1  f1  f1  f1 rr
 rr  rr  f1  f1  f1  f1  f1  rr rr
  0  rr  rr  f1  f1  f1  rr  rr  0
  0   0  rr  rr  rr  rr  rr   0  0
}

% --- Set library files to use

set optionsfile  "./finix_files/finix.options"
set rodfile      "./finix_files/finix.rods"
set scenariofile "./finix_files/finix.scenario"

% --- Rod types

rodtype F1 rod_B options_default scenario_default
           rgb 255 150 150
           multiplicity 264

% --- Radial meshes (input of power, output of temperature)

set inrmesh  1 0.0 0.40950
set outrmesh 1 0.0 0.40950

% --- Set axial meshing (centimetres)
%     set zmesh <nz>   <zlims>
%     <nz> = Number of axial layers (> 0 = equal spacing, < 0 = arbitrary spacing)
%     <zlims> = axial limits:
%        if <nz> > 0: <zlims> = <zmin> <zmax>
%        if <nz> < 0: <zlims> = <z1> <z2> ... <zn>

set zmesh -18
    0.00000
    5.71500
  11.250625
  16.78625
  22.321875
  27.8575
  33.393125
  38.92875
  44.464375
  50
  55.715
  61.250625
  66.78625
  72.321875
  77.8575
  83.393125
  88.928750
  94.464375
  100

% --- Minicore level case

lat SMR 1 0.0 0.0 9 9 21.50360
..  ..  ..  ..  ..  ..  ..  .. ..
..  ..  ..  F1  F1  F1  ..  .. ..
..  ..  F1  F1  F1  F1  F1  .. ..
..  F1  F1  F1  F1  F1  F1  F1 ..
..  F1  F1  F1  F1  F1  F1  F1 ..
..  F1  F1  F1  F1  F1  F1  F1 ..
..  ..  F1  F1  F1  F1  F1  .. ..
..  ..  ..  F1  F1  F1  ..  .. ..
..  ..  ..  ..  ..  ..  ..  .. ..

% --- Set the minicore lattice as the top level lattice

set toplattice SMR

% --- Horizontal plot

plot 3 500 500