Model Name

Contact: Model Developer, model_developer.at.institute.com (contact email)

Example of a page tag for the automated model indexing (remove # before tag)

Markdown
HTML

#!tag name=model_name
      description=<Short blurb describing the model>
      image=<Link to the image on the online website. See other model documentation pages for examples>
      pairs=reactor_type:MSR/FHR/SFR/...
                       reactor:<reactor name>
                       geometry:<reactor part being simulated>
                       simulation_type:<type of simulation>
                       input_features:<specific features in use>
                       transient:steady/transient/startup/RIA/...
                       codes_used:BlueCrab;Griffin;Pronghorn;MOOSE_NavierStokes;...
                       open_source:partial/true
                       computing_needs:workstation/HPC
                       fiscal_year:2024
                       institution:ANL/INL/...
                       sponsor:NRC/NEAMS/NEUP/...

#!tag name=model_name description=<Short blurb describing the model> image=<Link to the image on the online website. See other model documentation pages for examples> pairs=reactor_type:MSR/FHR/SFR/... reactor:<reactor name> geometry:<reactor part being simulated> simulation_type:<type of simulation> input_features:<specific features in use> transient:steady/transient/startup/RIA/... codes_used:BlueCrab;Griffin;Pronghorn;MOOSE_NavierStokes;... open_source:partial/true computing_needs:workstation/HPC fiscal_year:2024 institution:ANL/INL/... sponsor:NRC/NEAMS/NEUP/...

High Level Summary of Model

High level introduction is expected here about the model of interest. Figures and tables can be added in the model description, of which the examples are provided below.

Example 1: Example of adding a figure in VTB documentation.

Markdown
HTML

!media msr/msre/SAM_MSRE_1D.png
       style=width:40%
       id=msre_sam
       caption=The steady-state temperature distribution in the 1-D MSRE primary loop.

Figure 1: The steady-state temperature distribution in the 1-D MSRE primary loop.

note

One should put a copy of the figure at the specific file path. The VTB media is now arranged by reactor types under the folder /media. For more detailed instructions on figure formatting, please refer to the MooseDocs media tutorial.

Example 2: Example of adding a table in VTB documentation.

Markdown
HTML

!table id=fuel_salt_properties caption=Thermophysical properties of the fuel salt.
|   |   | Unit  | LiF-BeF$_4$-ZrF$_4$-UF$_4$  |
| :- | :- | :- | :- |
| Melting temperature | $T_{melt}$ | $K$ | $722.15$  |
| Density | $\rho$ | $kg/m^3$  | $2553.3-0.562\bullet T$ |
| Dynamic viscosity | $\mu$ | $Pa\bullet s$ | $8.4\times 10^{-5} exp(4340/T)$ |
| Thermal conductivity | $k$ | $W/(m\bullet K)$ | $1.0$ |
| Specific heat capacity | $c_p$ | $J/(kg\bullet K)$ | $2009.66$ |

Table 1: Thermophysical properties of the fuel salt.

		Unit	LiF-BeF $var element = document.getElementById("moose-equation-be4c4e65-cb3d-4058-aab4-901cd4d1466b");katex.render("_4", element, {displayMode:false,throwOnError:false});$ -ZrF $var element = document.getElementById("moose-equation-3368d00c-0179-44ee-b75b-dceedb291655");katex.render("_4", element, {displayMode:false,throwOnError:false});$ -UF $var element = document.getElementById("moose-equation-701470d3-d32d-49c8-8990-033856a1f849");katex.render("_4", element, {displayMode:false,throwOnError:false});$
Melting temperature	$var element = document.getElementById("moose-equation-c5e7651e-fc6d-46e8-8cad-a1d771dee62d");katex.render("T_{melt}", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-d6883830-bd46-4f69-aa5e-6980201224f8");katex.render("K", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-e7be1676-e481-4d7b-be48-88ad44dce982");katex.render("722.15", element, {displayMode:false,throwOnError:false});$
Density	$var element = document.getElementById("moose-equation-0f5f4b48-9af2-46b5-a5f2-5e195d75fa85");katex.render("\\rho", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-8dea4255-bfa2-42ac-8a8c-1d4dcf110c2a");katex.render("kg/m^3", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-95ebd9e6-06fb-4011-ab5c-a14a21ea1afa");katex.render("2553.3-0.562\\bullet T", element, {displayMode:false,throwOnError:false});$
Dynamic viscosity	$var element = document.getElementById("moose-equation-06613f67-5cd7-477e-b469-fbd7f4eee5c1");katex.render("\\mu", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-01a66a31-06c5-4753-b765-8c6f695c7d47");katex.render("Pa\\bullet s", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-a3aaaa89-db6f-425f-9668-80ebb756a69e");katex.render("8.4\\times 10^{-5} exp(4340/T)", element, {displayMode:false,throwOnError:false});$
Thermal conductivity	$var element = document.getElementById("moose-equation-9e676b76-29d4-49fe-9ea6-1a52975086e6");katex.render("k", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-32b7eb0b-7dba-4d18-9713-d1723c134dcf");katex.render("W/(m\\bullet K)", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-a95d801c-747b-4e95-8449-bdf41b8c2144");katex.render("1.0", element, {displayMode:false,throwOnError:false});$
Specific heat capacity	$var element = document.getElementById("moose-equation-82e328c6-5b3a-4b88-b1bf-e06a9f026fdf");katex.render("c_p", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-446f0dd9-40bc-4b4a-a2e8-3f7ca139f111");katex.render("J/(kg\\bullet K)", element, {displayMode:false,throwOnError:false});$	$var element = document.getElementById("moose-equation-9fd557bb-9ed5-47e7-8448-0b0fcbf0fe32");katex.render("2009.66", element, {displayMode:false,throwOnError:false});$

note

For more detailed instructions on table formatting, please refer to the MooseDocs table tutorial.

Example 3: Example of adding equations in VTB documentation.

Markdown
HTML

NekRS solves the incompressible Navier-Stokes equations:

\begin{equation}
\label{eq:mass}
\nabla\cdot\mathbf u=0
\end{equation}

\begin{equation}
\label{eq:momentum}
\rho_f\left(\frac{\partial\mathbf u}{\partial t}+\mathbf u\cdot\nabla\mathbf u\right)=-\nabla P+\nabla\cdot\tau+\rho\ \mathbf f
\end{equation}

\begin{equation}
\label{eq:energy}
\rho_f C_{p,f}\left(\frac{\partial T_f}{\partial t}+\mathbf u\cdot\nabla T_f\right)=\nabla\cdot\left(k_f\nabla T_f\right)+\dot{q}_f
\end{equation}

NekRS solves the incompressible Navier-Stokes equations:

$(1)var element = document.getElementById("moose-equation-10661c7e-785e-4c6c-908d-c33e89423408");katex.render(" \\nabla\\cdot\\mathbf u=0", element, {displayMode:true,throwOnError:false});$

$(2)var element = document.getElementById("moose-equation-7b567b41-5736-42b5-be81-bfc0e6946d03");katex.render(" \\rho_f\\left(\\frac{\\partial\\mathbf u}{\\partial t}+\\mathbf u\\cdot\\nabla\\mathbf u\\right)=-\\nabla P+\\nabla\\cdot\\tau+\\rho\\ \\mathbf f", element, {displayMode:true,throwOnError:false});$

$(3)var element = document.getElementById("moose-equation-42d317dc-3e5d-48b1-8597-259aa4ea0672");katex.render(" \\rho_f C_{p,f}\\left(\\frac{\\partial T_f}{\\partial t}+\\mathbf u\\cdot\\nabla T_f\\right)=\\nabla\\cdot\\left(k_f\\nabla T_f\\right)+\\dot{q}_f", element, {displayMode:true,throwOnError:false});$

note

For more detailed instructions on equation formatting, please refer to the MooseDocs equation tutorial.

Adding references

If the model has been previously published, one should provide the reference to the published works together with key references used for model development. MooseDocs supports BibTex style references, and some examples about how to cite technical report, conference proceeding and journal article are given below.


@techreport{Hu2017,
        address = {Lemont, IL},
        author = {Hu, Rui},
        number = {ANL/NE-17/4},
        mendeley-tags = {ANL/NE-17/4},
        institution = {Argonne National Laboratory},
        title = {{SAM Theory Manual}},
        year = {2017}
}


@InProceedings{novak_2021c,
  title       = {{Coupled {Monte} {Carlo} and Thermal-Hydraulics Modeling of a Prismatic Gas Reactor Fuel Assembly Using {Cardinal}}},
  author     = {A.J. Novak and D. Andrs and P. Shriwise and D. Shaver and P.K. Romano and E. Merzari and P. Keutelian},
  booktitle  = {{Proceedings of Physor}},
  year       = 2021
}


@article{Fang2021,
        author = {Fang, Jun and Shaver, Dillon R. and Min, Misun and Fischer, Paul
                and Lan, Yu-Hsiang and Rahaman, Ronald and Romano, Paul
                and Benhamadouche, Sofiane and Hassan, Yassin A.
                and Kraus, Adam and Merzari, Elia},
        doi = {10.1016/j.nucengdes.2021.111143},
        journal = {Nuclear Engineering and Design},
        title = {{Feasibility of Full-Core Pin Resolved CFD Simulations of
                Small Modular Reactor with Momentum Sources}},
        volume = {378},
        year = {2021}
}

The bibtex entries must be added into the VTB bibliography file vtb.bib. A reference list will be generated automatically at the end of the VTB documentation page.

Example 4: Example of inserting citations in VTB documentation.

Markdown
HTML

Description part A [!citep](Hu2017), description part B [!citep](novak_2021c), description part C [!citep](Fang2021)

Description part A (Hu, 2017), description part B (Novak et al., 2021), description part C (Fang et al., 2021)

Other key model details to provide in the high-level summary

Reactor Type(s): e.g., Liquid Metal Cooled Fast Reactor
Type of simulation: e.g., 3D core multiphysics (neutronics-TH) transient, 1D system steady-state
NEAMS MOOSE Codes Used: Griffin, Heat Conduction, Cardinal, etc.
NEAMS non-MOOSE Codes Used: MC2-3, Nek5000/RS, etc.
Non-NEAMS Codes Used: Cubit, OpenFOAM, etc.

note

More information is available online about the list of NEAMS codes and the codes that have been used for the models available in VTB repository. Providing the related information in documentation can better help the VTB team to categorize/sort the model.

Computational Model Description

Walk through the main kernel/blocks, show snippets of example inputs when needed. The markdown source can be as simple as the following to show the EOS block from the msre_loop.i input file:

Example 5: Example of including input block in VTB documentation.

Markdown
HTML

!listing msr/msre/steady_state/msre_loop_1d.i block=EOS language=cpp

[EOS]
  [./fuel_salt_eos]
    type = PTFunctionsEOS
    rho = fuel_salt_rho_func
    mu = fuel_salt_mu_func
    enthalpy = fuel_salt_enthalpy_func
    cp = 2009.66
    k = 1.0
  [../]
  [./hx_salt_eos]
    type = SaltEquationOfState
    salt_type = Flibe
  [../]
[]

Results

Here one can document the most important results from the related simulations. Please refer to aforementioned examples about how to add figures and tables when discussing the main results and discoveries.

Run Command

Please provide the run command for model execution on INL’s HPC Cluster or otherwise. Below is a simple example how to run a SAM system modeling job.


mpirun -np 48 /path/to/sam-opt -i sam_input.i

References

Jun Fang, Dillon R. Shaver, Misun Min, Paul Fischer, Yu-Hsiang Lan, Ronald Rahaman, Paul Romano, Sofiane Benhamadouche, Yassin A. Hassan, Adam Kraus, and Elia Merzari. Feasibility of Full-Core Pin Resolved CFD Simulations of Small Modular Reactor with Momentum Sources. Nuclear Engineering and Design, 2021. doi:10.1016/j.nucengdes.2021.111143.

@article{Fang2021,
    author = "Fang, Jun and Shaver, Dillon R. and Min, Misun and Fischer, Paul and Lan, Yu-Hsiang and Rahaman, Ronald and Romano, Paul and Benhamadouche, Sofiane and Hassan, Yassin A. and Kraus, Adam and Merzari, Elia",
    doi = "10.1016/j.nucengdes.2021.111143",
    journal = "Nuclear Engineering and Design",
    title = "{Feasibility of Full-Core Pin Resolved CFD Simulations of Small Modular Reactor with Momentum Sources}",
    volume = "378",
    year = "2021"
}

Rui Hu. SAM Theory Manual. Technical Report ANL/NE-17/4, Argonne National Laboratory, Lemont, IL, 2017.

@techreport{Hu2017,
    author = "Hu, Rui",
    address = "Lemont, IL",
    number = "ANL/NE-17/4",
    mendeley-tags = "ANL/NE-17/4",
    institution = "Argonne National Laboratory",
    title = "{SAM Theory Manual}",
    year = "2017"
}

A.J. Novak, D. Andrs, P. Shriwise, D. Shaver, P.K. Romano, E. Merzari, and P. Keutelian. Coupled Monte Carlo and Thermal-Hydraulics Modeling of a Prismatic Gas Reactor Fuel Assembly Using Cardinal. In Proceedings of Physor. 2021.

@InProceedings{novak_2021c,
    author = "Novak, A.J. and Andrs, D. and Shriwise, P. and Shaver, D. and Romano, P.K. and Merzari, E. and Keutelian, P.",
    title = "{Coupled {Monte} {Carlo} and Thermal-Hydraulics Modeling of a Prismatic Gas Reactor Fuel Assembly Using {Cardinal}}",
    booktitle = "{Proceedings of Physor}",
    year = "2021"
}

(msr/msre/steady_state/msre_loop_1d.i)

# Modeling the MSRE
# Steady state simulation
# Application : SAM
# If using or referring to this model, please cite as explained in
# https://mooseframework.inl.gov/virtual_test_bed/citing.html

[GlobalParams]
  global_init_P = 1e5                        # Global initial fluid pressure
  global_init_V = 0.001                      # Global initial fluid velocity
  global_init_T = 908.15                     # Global initial temperature for fluid and solid
  Tsolid_sf     = 1e-3
  gravity       = '0 -9.8 0'
  scaling_factor_var = '1 1e-3 1e-6'  # fluid model solver parameters
  p_order             = 2
[]

[Functions]
  [./fuel_salt_rho_func] # Linear fitting used by He, 2016
    type = PiecewiseLinear
    x    = '750       1200'
    y    = '2285.31   2032.41'
  [../]
  [./fuel_salt_enthalpy_func] # Approximated by Cp*T
    type = PiecewiseLinear
    x    = '750        1200'
    y    = '1.51E+06   2.41E+06'
  [../]
  [./fuel_salt_mu_func]
    type = PiecewiseLinear
    x    = '750  760  770  780  790  800  810  820
            830  840  850  860  870  880  890  900
            910  920  930  940  950  960  970  980
            990  1000 1010 1020 1030 1040 1050 1060
            1070 1080 1090 1100 1110 1120 1130 1140
            1150 1160 1170 1180 1190 1200'
    y    = '2.7378E-02 2.5371E-02 2.3557E-02 2.1915E-02 2.0424E-02 1.9069E-02 1.7834E-02 1.6706E-02
            1.5674E-02 1.4728E-02 1.3859E-02 1.3060E-02 1.2324E-02 1.1645E-02 1.1017E-02 1.0436E-02
            9.8976E-03 9.3976E-03 8.9328E-03 8.5001E-03 8.0969E-03 7.7206E-03 7.3690E-03 7.0402E-03
            6.7322E-03 6.4434E-03 6.1724E-03 5.9178E-03 5.6783E-03 5.4529E-03 5.2404E-03 5.0400E-03
            4.8508E-03 4.6720E-03 4.5029E-03 4.3428E-03 4.1911E-03 4.0473E-03 3.9109E-03 3.7813E-03
            3.6582E-03 3.5411E-03 3.4297E-03 3.3235E-03 3.2224E-03 3.1259E-03'
  [../]
[]

[EOS]
  [./fuel_salt_eos]
    type     = PTFunctionsEOS
    rho      = fuel_salt_rho_func
    mu       = fuel_salt_mu_func
    enthalpy = fuel_salt_enthalpy_func
    cp       = 2009.66
    k        = 1.0
  [../]
  [./hx_salt_eos]
    type      = SaltEquationOfState
    salt_type = Flibe
  [../]
[]

[MaterialProperties]
  [./alloy-mat]                              # Based on Hastelloy N alloy
    type = SolidMaterialProps
    k    = 23.6                              # Thermal conductivity
    Cp   = 578                               # Specific heat
    rho  = 8.86e3                            # Density
  [../]
[]

[Components]
  [./downcomer]
    type           = PBOneDFluidComponent
    A              = 0.1589
    Dh             = 0.0508
    length         = 1.7272
    n_elems        = 18
    orientation    = '0 -1 0'
    position       = '-0.7366 1.7272 0'
    eos            = fuel_salt_eos
  [../]
  [./j_dn_pl]
    type    = PBBranch
    Area    = 0.1155
    K       = '0.0 0.0'
    eos     = fuel_salt_eos
    inputs  = 'downcomer(out)'
    outputs = 'iplnm(in)'
  [../]
  [./iplnm]        # inlet plenum is connected to core bottom
    type           = PBOneDFluidComponent
    A              = 0.3932
    Dh             = 0.6997
    length         = 0.7366
    n_elems        = 8
    orientation    = '1 0 0'
    position       = '-0.7366 0 0'
    eos            = fuel_salt_eos
  [../]

  [./j_ip_c]
    type    = PBBranch
    Area    = 0.1155
    K       = '0.0 0.0'
    eos     = fuel_salt_eos
    inputs  = 'iplnm(out)'
    outputs = 'core(in)'
  [../]

  [./core]         # 1-D representation of the core
    type           = PBOneDFluidComponent
    A              = 0.3512 # consider a porosity of 0.225
    Dh             = 0.6687
    length         = 1.7272
    n_elems        = 20
    orientation    = '0 1 0'
    position       = '0 0 0'
    eos            = fuel_salt_eos
    heat_source    = 1.65e7
  [../]

  [./j_c_up]
    type    = PBBranch
    Area    = 0.1155
    K       = '0.0 0.0'
    eos     = fuel_salt_eos
    inputs  = 'core(out)'
    outputs = 'uplnm(in)'
  [../]

  [./uplnm]        # upper plenum is connected to core exit
    type           = PBOneDFluidComponent
    A              = 0.3932
    Dh             = 0.6997
    length         = 0.4346
    n_elems        = 6
    orientation    = '0 1 0'
    position       = '0 1.7272 0'
    eos            = fuel_salt_eos
  [../]

  [./j_up_ps1]
    type    = PBBranch
    Area    = 0.1155
    K       = '0.0 0.0'
    eos     = fuel_salt_eos
    inputs  = 'uplnm(out)'
    outputs = 'pipe1_s1(in)'
  [../]

  [./pipe1_s1]  # Connecting core to pump, horizontal section
    type        = PBOneDFluidComponent
    A           = 0.01267
    Dh          = 0.127
    length      = 1.8288
    n_elems     = 19
    orientation = '1 0 0'
    position    = '0 2.1618 0'
    eos         = fuel_salt_eos
  [../]

  [./j1]
    type    = PBBranch
    Area    = 0.01292
    K       = '0.0 0.0'
    eos     = fuel_salt_eos
    inputs  = 'pipe1_s1(out)'
    outputs = 'pipe1_s2(in)'
  [../]

  [./pipe1_s2]  # vertical section
    type        = PBOneDFluidComponent
    A           = 0.01267
    Dh          = 0.127
    length      = 0.8128
    n_elems     = 9
    orientation = '0 1 0'
    position    = '1.8288 2.1618 0'
    eos         = fuel_salt_eos
  [../]

  #
  # ====== pump ======
  #

  [./pump]
    type      = PBPump
    Area      = 0.01292
    K         = '0.15 0.1'
    eos       = fuel_salt_eos
    inputs    = 'pipe1_s2(out)'
    outputs   = 'pipe2(in)'
    initial_P = 1.1e5
  # Head_fn   = f_pump_head
    Head      = 43909.58
  []

  [./pipe2]     # Connecting the pump to HX
    type        = PBOneDFluidComponent
    A           = 0.01267
    Dh          = 0.127
    length      = 1.0668
    n_elems     = 11
    orientation = '-1 0 0'
    position    = '1.8288 2.9746 0'
    eos         = fuel_salt_eos
  [../]
  [./j2]    # junction connect to heat exchanger
    type    = PBBranch
    Area    = 0.01267
    K       = '0.0 0.0'
    eos     = fuel_salt_eos
    inputs  = 'pipe2(out)'
    outputs = 'hx_shell(in)'
  [../]

  #
  # ====== Customized U-tube heat exchanger ======
  #
  [./hx_shell]
    type        = PBOneDFluidComponent
    position    = '0.762 2.9746 0'
    orientation = '-1 0 0'
    length      = 2.5298
    n_elems     = 26
    eos         = fuel_salt_eos
    heat_source = 0
    A           = 1.0183E-01
    Dh          = 2.0945E-02
  [../]

  [./hx_tube1]
    type        = PBOneDFluidComponent
    position    = '-1.7678 3.0762 0'
    orientation = '1 0 0'
    length      = 2.5298
    n_elems     = 26
    eos         = hx_salt_eos
    heat_source = 0
    A           = 2.7885E-02
    Dh          = 1.0566E-02
    initial_T   = 824.8167
  [../]
  [./hx_j1]
    type      = PBBranch
    eos       = hx_salt_eos
    inputs    = 'hx_tube1(out)'
    outputs   = 'hx_tube2(in)'
    K         = '0 0'
    Area      = 2.7885E-02
    initial_T = 824.8167
  [../]
  [./hx_tube2]
    type        = PBOneDFluidComponent
    position    = '0.762 3.0762 0'
    orientation = '0 -1 0'
    length      = 0.2032
    n_elems     = 2
    eos         = hx_salt_eos
    heat_source = 0
    A           = 2.7885E-02
    Dh          = 1.0566E-02
    initial_T   = 824.8167
  [../]
  [./hx_j2]
    type      = PBBranch
    eos       = hx_salt_eos
    inputs    = 'hx_tube2(out)'
    outputs   = 'hx_tube3(in)'
    K         = '0 0'
    Area      = 2.7885E-02
    initial_T = 824.8167
  [../]
  [./hx_tube3]
    type        = PBOneDFluidComponent
    position    = '0.762 2.873 0'
    orientation = '-1 0 0'
    length      = 2.5298
    n_elems     = 26
    eos         = hx_salt_eos
    heat_source = 0
    A           = 2.7885E-02
    Dh          = 1.0566E-02
    initial_T   = 824.8167
  [../]

  [./hx_s_in]
    type  = PBTDJ
    input = 'hx_tube1(in)'
    eos   = hx_salt_eos
    v_bc  = 1.6
    T_bc  = 824.8167
  [../]
  [./hx_s_out]
    type  = PBTDV
    input = 'hx_tube3(out)'
    eos   = hx_salt_eos
    p_bc  = 1.0e5
    T_bc  = 866.4833
  [../]

  [./hx_wall1]
    type               = PBCoupledHeatStructure
    position           = '-1.7678 3.0762 0'
    orientation        = '1 0 0'
    length             = 2.5298
    hs_type            = cylinder
    radius_i           = 5.2832E-03
    width_of_hs        = 1.0668E-03
    elem_number_radial = 2
    elem_number_axial  = 26
    dim_hs             = 2
    material_hs        = 'alloy-mat'
    Ts_init            = 824.8167

    HS_BC_type                    = 'Coupled Coupled'
    name_comp_left                = hx_tube1
    HT_surface_area_density_left  = 8.6290E+02
    name_comp_right               = hx_shell
    HT_surface_area_density_right = 2.3629E+02
  [../]

  [./hx_wall2]
    type               = PBCoupledHeatStructure
    position           = '0.762 2.873 0'
    orientation        = '-1 0 0'
    length             = 2.5298
    hs_type            = cylinder
    radius_i           = 5.2832E-03
    width_of_hs        = 1.0668E-03
    elem_number_radial = 2
    elem_number_axial  = 26
    dim_hs             = 2
    material_hs        = 'alloy-mat'
    Ts_init            = 824.8167

    HS_BC_type                    = 'Coupled Coupled'
    name_comp_left                = hx_tube3
    HT_surface_area_density_left  = 8.6290E+02
    name_comp_right               = hx_shell
    HT_surface_area_density_right = 2.3629E+02
  [../]

  [./j3]
    type    = PBBranch
    Area    = 0.01267
    K       = '0.0 1e3 0.0'
    eos     = fuel_salt_eos
    inputs  = 'hx_shell(out) pipe_ref(out)'
    outputs = 'pipe3_s1(in)'
  [../]

  [./pipe_ref]
    type        = PBOneDFluidComponent
    A           = 0.01267
    Dh          = 0.127
    length      = 0.1
    n_elems     = 2
    orientation = '1 0 0'
    position    = '-1.8678 2.9746 0'
    eos         = fuel_salt_eos
  [../]
  [./ref_p]
    type  = PBTDV
    eos   = fuel_salt_eos
    T_bc  = 908.15
    p_bc  = 1.233351e+05
    input = 'pipe_ref(in)'
  [../]

  [./pipe3_s1]  # Connecting hx to downcomer
    type        = PBOneDFluidComponent
    A           = 0.01267
    Dh          = 0.127
    length      = 1.2474
    n_elems     = 13
    orientation = '0 -1 0'
    position    = '-1.7678 2.9746 0'
    eos         = fuel_salt_eos
  [../]

  [./j4]
    type    = PBBranch
    Area    = 0.01267
    K       = '0.0 0.0'
    eos     = fuel_salt_eos
    inputs  = 'pipe3_s1(out)'
    outputs = 'pipe3_s2(in)'
  [../]

  [./pipe3_s2]
    type        = PBOneDFluidComponent
    A           = 0.01267
    Dh          = 0.127
    length      = 1.0312
    n_elems     = 11
    orientation = '1 0 0'
    position    = '-1.7678 1.7272 0'
    eos         = fuel_salt_eos
  [../]

  [./j5]
    type    = PBBranch
    Area    = 0.01267
    K       = '0.0 0.0'
    eos     = fuel_salt_eos
    inputs  = 'pipe3_s2(out)'
    outputs = 'downcomer(in)'
  [../]
[]

[Postprocessors]
  [./Core_P_out]                               # Pressure at Core outlet/Loop inlet
    type     = ComponentBoundaryVariableValue
    variable = pressure
    input    = core(out)
  [../]
  [./Core_vel_in]
    type     = ComponentBoundaryVariableValue
    variable = velocity
    input    = core(in)
  [../]
  [./Core_T_out]
    type     = ComponentBoundaryVariableValue
    variable = temperature
    input    = core(out)
  [../]
  [./Core_T_in]
    type     = ComponentBoundaryVariableValue
    variable = temperature
    input    = core(in)
  [../]
  [./Core_P_in]
    type     = ComponentBoundaryVariableValue
    variable = pressure
    input    = core(in)
  [../]
  [./HX_Tin_p]
    type     = ComponentBoundaryVariableValue
    variable = temperature
    input    = hx_shell(in)
  [../]
  [./HX_Tout_p]
    type     = ComponentBoundaryVariableValue
    variable = temperature
    input    = hx_shell(out)
  [../]
[]

[Preconditioning]
  [./SMP_PJFNK]
    type                = SMP
    full                = true
    solve_type          = 'PJFNK'
    # pc_factor_shift are added automatically by SAM, they are added here for BlueCRAB
    petsc_options_iname = '-pc_type -ksp_gmres_restart -pc_factor_shift_type -pc_factor_shift_amount'
    petsc_options_value = 'lu 101 NONZERO 1e-9'
  [../]
[]

#[Problem]
#  restart_file_base = loop_1d_checkpoint_cp/0254
#[]

[Executioner]
  type       = Transient
  dt         = 0.2
  dtmin      = 1.e-3
  dtmax      = 10.0
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-6
  nl_max_its = 10
  l_tol      = 1e-6
  l_max_its  = 200
  start_time = 0
  end_time   = 300
  num_steps  = 100000
  [./Quadrature]
    type  = SIMPSON
    order = SECOND
  [../]
[]

[Outputs]
  print_linear_residuals = false
  perf_graph             = true
  [./out_displaced]
    type          = Exodus
    use_displaced = true
    execute_on    = 'initial timestep_end'
    sequence      = false
  [../]
  [./csv]
    type = CSV
  [../]
  [./checkpoint]
    type      = Checkpoint
    num_files = 1
  [../]
  [./console]
    type               = Console
    execute_scalars_on = 'none'
  [../]
[]

High Level Summary of Model
Computational Model Description
Results
Run Command
References