# SUBCHAN Proproc input

# the geometry data
geometry:
  # core data
  core:
    # the fuel assembly map in the reactor core. 1 stands for the fuel assembly while 0 is moderator assembly
    # Note: Do not pad this map with zeros around the borders. The 1's should go right to the edge
    core_map:
     - [0,   0,   0,   0,   1,   2,   3,   4,   5,   6,   7,   0,   0,   0,   0  ]
     - [0,   0,   8,   9,   10,  11,  12,  13,  14,  15,  16,  17,  18,  0,   0  ]
     - [0,   19,  20,  21,  22,  23,  24,  25,  26,  27,  28,  29,  30,  31,  0  ]
     - [0,   32,  33,  34,  35,  36,  37,  38,  39,  40,  41,  42,  43,  44,  0  ]
     - [45,  46,  47,  48,  49,  50,  51,  52,  53,  54,  55,  56,  57,  58,  59 ]
     - [60,  61,  62,  63,  64,  65,  66,  67,  68,  69,  70,  71,  72,  73,  74 ]
     - [75,  76,  77,  78,  79,  80,  81,  82,  83,  84,  85,  86,  87,  88,  89 ]
     - [90,  91,  92,  93,  94,  95,  96,  97,  98,  99,  100, 101, 102, 103, 104]
     - [105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119]
     - [120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134]
     - [135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149]
     - [0,   150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 0  ]
     - [0,   163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 0  ]
     - [0,   0,   176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 0,   0  ]
     - [0,   0,   0,   0,   187, 188, 189, 190, 191, 192, 193, 0,   0,   0,   0  ]
    
    # the centroid-to-centroid distance of the connected assembly in units of m.
    # User can also input a digit with unit, such as 0.215m, 21.5cm, 215mm, .etc.
    assembly_pitch: 215mm
    # axial mesh information

    # # the length of the active region of the core
    # active_region_length: 3.6576m

    axial_mesh_number: 10
    axial_mesh_grids: [11.951, 48.527, 85.103, 121.679, 158.255, 194.831, 231.407, 267.983, 304.559, 341.135, 377.711]
    axial_mesh_grids_units: cm

  # assembly data
  assembly:
    # the assembly map where 0 stands for fuel rod and 1 stands for guide tube or instrument tube
    assembly_map:
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
      - [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    # the centroid-to-centroid distance of the connected fuel rod.
    pin_pitch: 12.6mm
    walls_around_bundle: no
    # the type of the heated elements, nuclear_fuel_rod or electric heater tubes
    heated_elements_type: nuclear_fuel_rods
    # conduction model flag, no_conduction, radial_conduction, radial_axial_conduction, radial_axial_azimuthal_conduction
    conduction_model: radial_conduction

  fuel_rods:
    # the radial nodes in the fuel pellet
    pellet_radial_nodes: 10
  # the fuel rod data
    fuel_rod:
      - id: 1
        # the fuel rod diameter which is the sum of the thickness of the fuel pellet, gap, and the clad
        outer_diameter: 9.5mm
        # the fuel pellet properties
    
        # fuel pellet diameter in units of meter.
        # User can also input a digit with unit, such as 0.215m, 21.5cm, 215mm, .etc.
        inner_diameter: 8.192mm
  
        # the fuel density in units of kg/m3.
        # User can also define a density like 10.97g/cm3,
        fuel_density: 10.97g/cm3
        # fuel thermal conductivity in units of W/m/K
        fuel_thermal_conductivity: -1
        # the fuel specific heat in units of J/kg/K
        fuel_specific_heat: -1
        # fitting factor for heat source distribution in the fuel pellet
        fitting_factor: 0.
    
        # the cladding properties
    
        # the cladding thickness is needed in SUBCHAN input.
        # User can define a cladding_thickness directly or two parameters: clad_inner_diameter and clad_outer_diameter
        clad_inner_diameter: 8.36mm
        clad_outer_diameter: 9.5mm
        #  clad_thickness: 0.57mm
        clad_density: 6.552g/cm3
        #  clad_material: zircaloy
        # the clad thermal conductivity in units of W/m/K
        clad_thermal_conductivity: -1
        # the clad specific heat in units of J/kg/K
        clad_specific_heat: -1
    
        # fuel to clad heat transfer coefficient in units of W/m2/K
        gap_conductance: 5678.3
        clad_direct_heating: 0.
    
        # gap properties
    
        # the cold fill gas pressure in units of W/m2/K
        gas_pressure: 1.0292E+04
        gas_plenum_volume: 6.488E+03
        he_molar_fraction: 1.204E-02
        xe_molar_fraction: -3.46
        ar_molar_fraction: 2.428E+02
        kr_molar_fraction: 1.5224E+01
        h_molar_fraction: 2.428E+02
        n_molar_fraction: 5.674E+03

  # the guide tube or instrument tube data
  guide_tube:
    inner_diameter: 11.22mm
    outer_diameter: 12.04mm
    clad_thickness: 0.41mm
    guide_tube_material: zircaloy

  # the spacer grids data
  spacer_grids:
    number: 7
    # each spacer grid initial height and final height in units of mm
    position:
      - [ 0, 158.17 ]
      - [ 613.44, 771.05 ]
      - [ 1135.44, 1293.05 ]
      - [ 1657.44, 1815.05 ]
      - [ 2179.44, 2337.05 ]
      - [ 2701.44, 2859.05 ]
      - [ 3223.44, 3381.05 ]
    position_units: mm
    grid_type: [ 1, 1, 1, 1, 1, 1, 1 ]
    grid_minor_loss_coefficient: [0.9065]

# control data
control:
  # the running identifier for SUBHCAN output
  running_identifier: VERA

  # the type of analysis, subchannel(default) or core
  analysis_type: subchannel

  # the units option, SI(default), AE
  units_option: SI

  # thermal hydraulic model,
  # COBRA_THREE: COBRA three-equation model (mixture mass, energy and momentum) (default)
  # COBRA_FOUR: COBRA four-equation model (vapor mass is added)
  # TWIGL_FOUR: TWIGL four-equation model
  thermal_hydraulic_model: COBRA_THREE

  # the power source data
  # SUBCHAN.IN: power source data is in SUBCHAN.IN file
  # POWRSRC.CPL: power source data is in POWRSRC.CPL file for coupling calculations
  power_source: POWRSRC.CPL

  # the solver type: gauss_elimination,  jocobi_iteration (default)
  solver: jocobi_iteration

  # the hydraulic models
  hydraulic_models:
    turbulent_mixing: 0
    single_phase_friction: 0
    two_phase_friction: 1
    void_friction: 1
    inlet_flow_boron: 0
    parameter: 1
    convergence: 1
    physical_properties: 0
    heat_transfer: 1
    water_properties: 0

  # the specific hydraulic model parameters
  hydraulic_models_parameters:
    # the turbulent mixing model parameters
    turbulent_mixing:
      # Coefficient a and b of the subcooled turbulent mixing model
      coefficients: [ 0., 0. ]
      nbbc: 1
      # Thermal conduction geometry factor (if 0.0, thermal conduction is neglected)
      thermal_conduction_geometry_factor: 0.
      # present only if NBBC > 1
      two_phase_mixing_coefficients: null

    # single_phase friction model
    single_phase_friction:
      wall_viscosity_correction: 0
      turbulent_constants: null
      laminar_flow_constants: null

    # the two phase friction model
    two_phase_friction:
      # the two phase friction multiplier option, 1: EPRI correlation, 2: Homogeneous model(default),
      # 3: Armand correlation, 4: not used, 5: Polynomial in quality
      multiplier_option: 2
      # present only if multiplier_option = 5 on the previous card
      polynomial: null

    # the void fraction model
    void_friction:
      # the subcooled boiling option, 1: EPRI correlation, 2: Levy correlation,
      #　3: Homogeneous model i.e., no subcooled boiling(default)
      subcooled_boiling: 3
      # Indicator of bulk void/quality relationship, Homogeneous model(default). See SUBCHAN document
      bulk_void_quality_relationship: 3
      # Flag for slip ratio correlation (used only if the COBRA or TWIGL four-equation model has been activated in col. 19-24 of card 2)
      slip_ratio_correlation: 1
      # slip ratio or void fraction polynomial
      bulk_void_polynomial: null
      slip_ratio_polynomial: null

    # the heat transfer model
    heat_transfer:
      # Correlation for the heat transfer coefficient in single phase
      single_phase: 0
      # Correlation for the heat transfer coefficient in the subcooled nucleate boiling region
      subcooled_boiling: 0
      # Correlation for the heat transfer coefficient in the saturated nucleate boiling region
      saturated_boiling: 0
      # Correlation to determine the CHF point on the boiling curve
      chf_point: 0
      # Correlation for the heat transfer coefficient in the transition boiling region
      transition_boiling: 0
      # Correlation for the heat transfer coefficient in the film boiling region
      film_boiling: 0
      # One card required only if single_phase = 2
      coefficients: null

    # inlet flow and boron model
    inlet_flow_boron:
      # Inlet flow option:
      flow: 0
      flow_distribution: null

      # inlet boron option
      boron: 0
      boron_distribution: null

    # parameter model
    parameter:
      cross_flow_resistance_coefficient: 0.5
      turbulent_momentum_factor: 0.0
      transverse_momentum_parameter: 0.5
      theta: 0.0

      chf_correlation: 3
      cold_wall_correction: 0
      nonuniform_axial_flux_correction: 0
      grid_loss_coefficient: 0.0
      critical_power_correlation: 0
      cpr_iterations_maximum_number: 0
      gexl_correlation_ppf: 0.
      cise_correlation_ppf: 0.
      ctc_spacing: 0.

    # convergence parameters:
    convergence:
      external_iterations_maximum: 0
      external_iterations_minimum: 0
      axial_flow_convergence: 0.
      cross_flow_convergence: 0.
      damping_factor_pressure: 0.
      damping_factor_flow: 0.
      fuel_htc_convergence: -0.01
      fuel_temperature_convergence: 0.
      internal_iterations_maximum: 0
      maximum_relative_continuity_error: 0.
      pressure_iterations_convergence: 0.
      pressure_gradient_iteration_scheme: 1

    physical_properties:
      subcooled_liquid: 0
      superheated_stream: 0
      subcooled_liquid_minimum_enthalpy: 0.
      superheated_steam_maximum_enthalpy: 0.

  # the operating condictions that include the enthalpy, mass flux, etc.
  operating_conditions:
    # the inlet enthalpy option,
    # 0: the same inlet enthalpy for all channels
    # 1: the same inlet temperature for all channels
    # 2: inlet enthalpy specified by channel
    # 3: inlet temperature specified by channel
    inlet_enthalpy_type: 0
    # the inlet enthalpy in units C, K, F, J/kg
    inlet_enthalpy: 1.282e+06
    # the inlet mass flow rate in units (kg/s).
    # There are 3 options for user to input the inlet mass flow rate:
    # 1. inlet_mass_flux: the total inlet mass flux in units of kg/s
    # 2. assembly_mass_flux: the inelt mass flux of single assembly in units of kg/s.
    # 3. average_inlet_mass_flux: the average inlet mass flux for each channels in units of kg/s/m2
    # inlet_mass_flux: 16591.67925 kg/s
    assembly_mass_flux: 85.9683 kg/s
    # average_inlet_mass_flux: 3487 kg/s/m2
    
    # the exit pressure in units of MPa
    exit_pressure: 15.51 Mpa
    # dps
    dps: 0.
    # ips
    ips: 0.
    # power normalization factor multiplying the power distribution supplied in Card 5 Linear Rod Powers for each channel and axial interval
    power_normalization_factor: 1.
    # fission power appearing as direct heating of the coolant
    moderator_direct_heating: 0.026
    # Average inlet mass flux for bypass channels in units of kg/m2/s
    inlet_mass_flux_bypass: 0.
    # inlet boron concentration for coolant channels in units of ppm
    inlet_boron_concentration: 0.
    # Fraction of the heat directly released outside the fuel rod which appears inside the channel
    channel_direct_heating: 0.
    # the Exit enthalpy (Btu/lb or J/kg), used only in the case of reverse flow at core outlet
    exit_enthalpy: 0.
    # inlet enthalpy or temperature distribution for all channels in sequence
    enthalpy_temperature_distribution: null

  # time controller
  time_controller:
    # Transient end time (s). If zero or blank, the program stops after steady state
    transient_end_time: 0.
    # number of time domains
    time_domains_number: 3
    # number of slugs of borated water into which the transient time is divided
    borated_water_slugs_number: 0
    # Starting time (s) for each slug N = 1, borated_water_slugs_number
    each_slug_starting_time: null
    # time_parameters: time step size, long output frequency, end time
    time_steps:
      - [ 0.01, 10, 0.5 ]
      - [ 0.01, 10, 1. ]
      - [ 0.01, 10, 2. ]

  # transient parameters
  transient_controller:
    # Number of time points at which pressure transient forcing function will be given
    pressure: 0
    pressure_forcing_functions: null
    # Number of time points at which power distribution forcing function will be given
    power_distribution: 0
    power_forcing_functions: null
    # Number of transient forcing functions for inlet enthalpy (or temperature), inlet flow (or specified pressure drop if DPS > 0)
    # and inlet Boron concentration (max = 10)
    inlet_enthalpy: 0
    inlet_enthalpy_forcing_functions: null
    inlet_flow: 0
    inlet_flow_forcing_functions: null
    inlet_boron_concentration: 0
    inlet_boron_forcing_functions: null

  # output controller
  output_controller:
    # Axial print option,
    # 1: every axial step printed
    # >1: only axial steps which are multiple of this numbe
    axial_print: 1
    # Time print option for short output
    time_print: 0
    # channels print options
    # 0: channels results only
    # 1: channels results + cross flow table
    # 2: channels and fuel rods results
    # 3: channel and fuel rod results + cross flow table (default)
    results_print: 3
    # channels to be printed, one or more cards specifying the identification numbers of the channels to be printed.
    channels_print: null
    # rods to be printed
    rods_print: null
    # nodes to be printed
    fuel_radial_nodes_print: null
    # plot print 0, the thermal-hydraulic plotfile is required
    plot_print: 0

# power data
power:
  # the total power in units of MW, User can also define the total power by 3.E9W
  total_power: 3411

#  2d_distributions:
#    axial_power_factor: []
#    assembly_power_distributions:
#      - []
#      - []
#      - []

  # the 3-dimensional power distribution input
  3d_distributions: null
#    - []
#    - []
#    - []


























