Divertor Monoblock During Pulsed Operation

TMAP8 is used to model tritium transport in a divertor monoblock to elucidate the effects of pulsed operation (up to fifty 1600-second plasma discharge and cool-down cycles) on the tritium in-vessel inventory source term and ex-vessel release term (i.e., tritium retention and permeation) for safety analysis. This example reproduces the results presented in Shimada et al. (2024).

General description of the simulation case and corresponding input file

Introduction

In a magnetic confinement fusion system, the divertor components will be subjected to intense particle and heat fluxes from plasma, as well as 14-MeV neutrons stemming from deuterium-tritium (D-T) reactions, thus creating extremely high temperature and hydrogen concentration gradients across tens of millimeters. Furthermore, the thermomechanical properties of the divertor materials will evolve over time in response to displacement damage and gas/solid transmutations. Thus, predictive modeling of tritium transport in the divertor poses enormous challenges, as it requires simulation capabilities capable of (1) multi-material configurations, (2) high thermal cycles, (3) complex 2D and 3D geometries, (4) and microstructural evolution resulting from displacement damage and gas/solid transmutations. For example, the ITER divertor, which consists of tungsten (W) monoblocks bonded to cooling tubes made of copper-chromium-zirconium (CuCrZr) alloy, is designed to withstand high heat fluxes (∼10 MW $var element = document.getElementById("moose-equation-ef9cd234-f13a-4def-bb37-b889accbedca");katex.render("m^{\u22122}", element, {displayMode:false,throwOnError:false});$ ) during steady-state plasma operation from intense D-T plasma flux (∼1024 $var element = document.getElementById("moose-equation-e5f18615-437c-43ef-96d0-670259c0116a");katex.render("m^{\u22122}s^{\u22121}", element, {displayMode:false,throwOnError:false});$ ). In this example, we simulate tritium retention and permeation as well as thermal transport in an ITER-like 2D W-Cu-CuCrZr monoblock and demonstrate the first three (out of four) aforementioned required simulation capabilities for tritium transport in the divertor using TMAP8, as presented in Shimada et al. (2024). This simulation was also based on the cases published in Hodille et al. (2021).

The sections below describe the simulation details and explain how they translate into the example TMAP8 input file. Not every part of the input file is explained here. If the interested reader has more questions about this example case and how to modify this input file to adapt it to a different case, feel free to reach out to the TMAP8 development team on the TMAP8 GitHub discussion page.

Divertor monoblock geometry and mesh generation

Figure 1 shows the geometry, mesh, and temperature distribution of the 2D monoblock employed in Shimada et al. (2024) and in the present example. The 2D monoblock consisted of three materials:

W monoblock ( $var element = document.getElementById("moose-equation-6ded2834-34d3-449c-a835-fbe5eb2d98bc");katex.render("r", element, {displayMode:false,throwOnError:false});$ (mm) > 8.5),
Cu interlayer (7.5 < $var element = document.getElementById("moose-equation-c645a5dc-0bfb-4080-8e97-913562fd9292");katex.render("r", element, {displayMode:false,throwOnError:false});$ (mm) < 8.5),
CuCrZr tube (6.0 < $var element = document.getElementById("moose-equation-6cadaa23-c345-40f4-b660-5020b73711d9");katex.render("r", element, {displayMode:false,throwOnError:false});$ (mm) < 7.5), and
H $var element = document.getElementById("moose-equation-7d4ce290-9c25-4396-8a2a-97ed2dfa66f5");katex.render("_2", element, {displayMode:false,throwOnError:false});$ O cooling ( $var element = document.getElementById("moose-equation-39cc0534-34c7-4c54-8303-208f932075ee");katex.render("r", element, {displayMode:false,throwOnError:false});$ (mm) < 6.0),

where $var element = document.getElementById("moose-equation-461c6e20-0f39-41d2-bfdc-8f7c1791f71d");katex.render("r=\\sqrt{x^2+y^2}", element, {displayMode:false,throwOnError:false});$ in a Cartesian coordinate system, and the center point ( $var element = document.getElementById("moose-equation-b1ec4430-7ef4-45ea-8716-8b458aecf376");katex.render("x", element, {displayMode:false,throwOnError:false});$ , $var element = document.getElementById("moose-equation-da12c010-9a24-4187-b2e6-a3b5e2d08a90");katex.render("y", element, {displayMode:false,throwOnError:false});$ ) = (0, 0) is defined at the center of the CuCrZr tube.

Figure 1: 2D monoblock: (left) geometry and mesh; (right) temperature distribution. This corresponds to Fig. 1 in Ref. Shimada et al. (2024).

MOOSE is equipped with a set of user-friendly, built-in mesh generators for creating meshes based on simple geometries (e.g., a monoblock). The built-in mesh generator ConcentricCircleMeshGenerator is used to create meshes for the 2D monoblock geometry. Only the left-hand (−14.0 < $var element = document.getElementById("moose-equation-3a487c51-c350-41fc-99a3-a8fa3f540913");katex.render("x", element, {displayMode:false,throwOnError:false});$ (mm) < 0 & −14.0 < $var element = document.getElementById("moose-equation-4ec60aa8-8e2b-4721-9ce7-760759db0fdb");katex.render("y", element, {displayMode:false,throwOnError:false});$ (mm) < 14.0) portion of the monoblock is used, such that the mesh size is halved by assuming symmetry at the $var element = document.getElementById("moose-equation-135b1169-f0e4-4d70-a0d1-82f407aceca2");katex.render("x = 0", element, {displayMode:false,throwOnError:false});$ plane. This example uses a total of 9,418 nodes, creating 28,402 degrees of freedom in the non-linear system. The finer mesh size and time step are required to compute the trapping-limited diffusion regime, as described in TMAP8 verification case ver-1d. However, a set of relatively low detrapping energies (< 0.85 eV) is used here to reduce the computing time via the relatively coarse mesh size.

In the input file, the mesh is defined as:

[Mesh<<<{"href": "../../syntax/Mesh/index.html"}>>>]
  [ccmg]
    type = ConcentricCircleMeshGenerator<<<{"description": "This ConcentricCircleMeshGenerator source code is to generate concentric circle meshes.", "href": "../../source/meshgenerators/ConcentricCircleMeshGenerator.html"}>>>
    num_sectors<<<{"description": "num_sectors % 2 = 0, num_sectors > 0Number of azimuthal sectors in each quadrant'num_sectors' must be an even number."}>>> = 36
    rings<<<{"description": "Number of rings in each circle or in the enclosing square"}>>> = '1 30 20 110'
    radii<<<{"description": "Radii of major concentric circles"}>>> = '${units 6 mm -> m} ${units 7.5 mm -> m} ${units 8.5 mm -> m}'
    has_outer_square<<<{"description": "It determines if meshes for a outer square are added to concentric circle meshes."}>>> = on
    pitch<<<{"description": "The enclosing square can be added to the completed concentric circle mesh.Elements are quad meshes."}>>> = '${units 28 mm -> m}'
    portion<<<{"description": "Control of which part of mesh is created"}>>> = left_half
    preserve_volumes<<<{"description": "Volume of concentric circles can be preserved using this function."}>>> = false
    smoothing_max_it<<<{"description": "Number of Laplacian smoothing iterations"}>>> = 3
  []
  [ssbsg1]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = ccmg
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '4' # W
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '3' # Cu
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '4to3'
  []
  [ssbsg2]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = ssbsg1
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '3' # Cu
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '4' # W
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '3to4'
  []
  [ssbsg3]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = ssbsg2
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '3' # Cu
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '2' # CuCrZr
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '3to2'
  []
  [ssbsg4]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = ssbsg3
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '2' # CuCrZr
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '3' # Cu
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '2to3'
  []
  [ssbsg5]
    type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = ssbsg4
    primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '2' # CuCrZr
    paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '1' # H2O
    new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '2to1'
  []
  [bdg]
    type = BlockDeletionGenerator<<<{"description": "Mesh generator which removes elements from the specified subdomains", "href": "../../source/meshgenerators/BlockDeletionGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = ssbsg5
    block<<<{"description": "The list of blocks to be deleted"}>>> = '1' # H2O
  []
[]

Nomenclature of variables and physical parameters

Table 1 lists the variables and physical parameters used in this example with their units.

Table 1: Nomenclature of variables and physical parameters used in this example.

Variables and governing equations

To simulate tritium and thermal transport, we define two sets of partial differential equations (PDEs). First, the strong form of the mass conservation equation for solute (mobile) T atoms, $var element = document.getElementById("moose-equation-168f1f7e-9892-47a7-86e3-1c12096c56cd");katex.render("C_s", element, {displayMode:false,throwOnError:false});$ , is written as:

$(1)var element = document.getElementById("moose-equation-0ae3f778-2a11-4602-b4a3-1b8c730447ef");katex.render(" \\frac{\\partial C_s}{\\partial t} + \\nabla \\cdot (-D \\nabla C_s) + \\alpha_t \\frac{C_t^e}{N} C_s - \\alpha_r C_t = 0.", element, {displayMode:true,throwOnError:false});$

We use three sets of mass conservation equations to calculate the behaviors of solute T atoms in three different materials (i.e., W, Cu, CuCrZr). Second, the strong form of the conservation of energy equation is written as:

$(2)var element = document.getElementById("moose-equation-e6534ced-5136-4c04-a4de-08993a97b2c0");katex.render(" \\rho c_p \\frac{\\partial T}{\\partial t} - \\nabla \\cdot (k_T \\nabla T) = 0.", element, {displayMode:true,throwOnError:false});$

Table 1 lists all the symbols and definitions used in Eq. (1) and Eq. (2).

The next step is to convert these two strong-form PDEs into their weak forms by multiplying with a test function, $var element = document.getElementById("moose-equation-6436e8fb-309f-4658-b621-edeb87a464f6");katex.render("\\psi", element, {displayMode:false,throwOnError:false});$ , and integrating over a domain, $var element = document.getElementById("moose-equation-add840d0-037f-4c89-af55-1ba9c4555b65");katex.render("\\Omega", element, {displayMode:false,throwOnError:false});$ with surface $var element = document.getElementById("moose-equation-ffd14630-ae9f-4690-b53f-73110750e99f");katex.render("\\delta \\Omega", element, {displayMode:false,throwOnError:false});$ and outward-facing normal vector $var element = document.getElementById("moose-equation-4a13f6f0-7af5-4289-80d6-29a06510c713");katex.render("\\hat{n}", element, {displayMode:false,throwOnError:false});$ . Using the divergence theorem, one can obtain the weak form of the mass conservation equation (Eq. (1)) as follows:

$(3)var element = document.getElementById("moose-equation-d99f73ac-2a72-427e-a102-4795826abaf1");katex.render(" \\int_{\\Omega}\\psi \\frac{\\partial C_s}{\\partial t} -\\int_{\\delta \\Omega} \\psi D \\nabla C_s \\cdot \\hat{n} +\\int_{\\Omega} \\nabla \\psi \\cdot D \\nabla C_s +\\int_{\\Omega} \\psi \\alpha_t \\frac{C_t^e}{N} C_s -\\int_{\\Omega} \\psi \\alpha_r C_t = 0.", element, {displayMode:true,throwOnError:false});$

Similarly, the weak form of the conservation of energy equation can be written by multiplying Eq. (2) with a test function, $var element = document.getElementById("moose-equation-60cfef30-6fbc-4136-84b6-c460096d27d6");katex.render("\\psi", element, {displayMode:false,throwOnError:false});$ , and integrating over a domain, $var element = document.getElementById("moose-equation-daae29b2-8a63-4eb6-8b62-fe0b1ec16475");katex.render("\\Omega", element, {displayMode:false,throwOnError:false});$ with surface $var element = document.getElementById("moose-equation-c29eba6c-7269-4a73-99d5-4db20da480de");katex.render("\\delta \\Omega", element, {displayMode:false,throwOnError:false});$ and outward-facing normal vector $var element = document.getElementById("moose-equation-862aef57-138c-4ba7-a95f-bf28a332c131");katex.render("\\hat{n}", element, {displayMode:false,throwOnError:false});$ . Using the divergence theorem, one can obtain the weak form of the conservation of energy equation as follows:

$(4)var element = document.getElementById("moose-equation-664c10c9-b705-44ac-9037-f11a608597f5");katex.render(" \\int_{\\Omega}\\psi c_p \\frac{\\partial T}{\\partial t} -\\int_{\\delta \\Omega} \\psi k_T \\nabla T \\cdot \\hat{n} +\\int_{\\Omega} \\nabla \\psi \\cdot k_T \\nabla T = 0.", element, {displayMode:true,throwOnError:false});$

Then, to solve for the PDEs and physical phenomena, we can select appropriate kernels and boundary conditions (BCs) from MOOSE’s extensive library. The following three subsections describe each kernel, BC, and numerical method utilized in the present work.

In the input file, the variables are defined as:

[Variables<<<{"href": "../../syntax/Variables/index.html"}>>>]
  [temperature]
    order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
    family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = '${units 300 K}'
  []
  ######################### Variables for W (block = 4)
  [C_mobile_W]
    order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
    family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = '${units 1.0e-20 m^-3}'
    block = 4
  []
  [C_trapped_W]
    order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
    family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = '${units 1.0e-15 m^-3}'
    block = 4
  []
  ######################### Variables for Cu (block = 3)
  [C_mobile_Cu]
    order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
    family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = '${units 5.0e-17 m^-3}'
    block = 3
  []
  [C_trapped_Cu]
    order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
    family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = '${units 1.0e-15 m^-3}'
    block = 3
  []
  ######################### Variables for CuCrZr (block = 2)
  [C_mobile_CuCrZr]
    order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
    family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = '${units 1.0e-15 m^-3}'
    block = 2
  []
  [C_trapped_CuCrZr]
    order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
    family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
    initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = '${units 1.0e-15 m^-3}'
    block = 2
  []
[]

Note the usage of initial_condition and block parameters in order to set initial conditions and material block applicability for each variable.

Kernels

For the mass conservation equation, we use ADTimeDerivative and ADMatDiffusion to represent the 1 $var element = document.getElementById("moose-equation-7a970a4b-7b99-4064-810e-7644db8c9425");katex.render("^{\\text{st}}", element, {displayMode:false,throwOnError:false});$ and 3 $var element = document.getElementById("moose-equation-be417e8f-6d8c-4243-b737-7d1b96d999f9");katex.render("^{\\text{rd}}", element, {displayMode:false,throwOnError:false});$ terms of Eq. (3). The TMAP8 kernels TrappingNodalKernel and ReleasingNodalKernel represent the 4 $var element = document.getElementById("moose-equation-77f38918-02ae-493c-a574-d3c038a052c0");katex.render("^{\\text{th}}", element, {displayMode:false,throwOnError:false});$ and 5 $var element = document.getElementById("moose-equation-f07d6082-2c8f-4b7a-acae-7983b80c7793");katex.render("^{\\text{th}}", element, {displayMode:false,throwOnError:false});$ terms of Eq. (3), respectively, and simulate the trapping/release behavior of hydrogen isotopes in/from trap sites. For the conservation of energy equation, SpecificHeatConductionTimeDerivative and HeatConduction solve the 1 $var element = document.getElementById("moose-equation-779007af-5d7e-44cf-942f-a6d6330c67c6");katex.render("^{\\text{st}}", element, {displayMode:false,throwOnError:false});$ and 3 $var element = document.getElementById("moose-equation-fbc41932-1053-4d27-a86e-aea750496458");katex.render("^{\\text{rd}}", element, {displayMode:false,throwOnError:false});$ terms of Eq. (4). ADTimeDerivative and ADMatDiffusion are kernels from the core MOOSE framework, and SpecificHeatConductionTimeDerivative and HeatConduction are kernels from the MOOSE Heat Transfer Module. These pre-made kernels are commonly re-used to represent time derivative, diffusion, and heat conduction terms in a given material within a variety of MOOSE-based simulations.

[Kernels<<<{"href": "../../syntax/Kernels/index.html"}>>>]
  ############################## Kernels for W (block = 4)
  [diff_W]
    type = ADMatDiffusion<<<{"description": "Diffusion equation kernel that takes an isotropic diffusivity from a material property", "href": "../../source/kernels/ADMatDiffusion.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
    diffusivity<<<{"description": "The diffusivity value or material property"}>>> = diffusivity_W
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [time_diff_W]
    type = ADTimeDerivative<<<{"description": "The time derivative operator with the weak form of $(\\psi_i, \\frac{\\partial u_h}{\\partial t})$.", "href": "../../source/kernels/ADTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [coupled_time_W]
    type = ScaledCoupledTimeDerivative<<<{"description": "Time derivative Kernel that acts on a coupled variable. Weak form: $(\\psi_i, \\frac{\\partial v_h}{\\partial t})$.", "href": "../../source/kernels/ScaledCoupledTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
    v<<<{"description": "Coupled variable"}>>> = C_trapped_W
    factor<<<{"description": "The factor by which to scale"}>>> = 1e0
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [heat_conduction_W]
    type = HeatConduction<<<{"description": "Diffusive heat conduction term $-\\nabla\\cdot(k\\nabla T)$ of the thermal energy conservation equation", "href": "../../source/kernels/HeatConduction.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    diffusion_coefficient<<<{"description": "Property name of the diffusion coefficient"}>>> = thermal_conductivity_W
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [time_heat_conduction_W]
    type = SpecificHeatConductionTimeDerivative<<<{"description": "Time derivative term $\\rho c_p \\frac{\\partial T}{\\partial t}$ of the heat equation with the specific heat $c_p$ and the density $\\rho$ as arguments.", "href": "../../source/kernels/SpecificHeatConductionTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    specific_heat<<<{"description": "Property name of the specific heat material property"}>>> = specific_heat_W
    density<<<{"description": "Property name of the density material property"}>>> = density_W
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  ############################## Kernels for Cu (block = 3)
  [diff_Cu]
    type = ADMatDiffusion<<<{"description": "Diffusion equation kernel that takes an isotropic diffusivity from a material property", "href": "../../source/kernels/ADMatDiffusion.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_Cu
    diffusivity<<<{"description": "The diffusivity value or material property"}>>> = diffusivity_Cu
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [time_diff_Cu]
    type = ADTimeDerivative<<<{"description": "The time derivative operator with the weak form of $(\\psi_i, \\frac{\\partial u_h}{\\partial t})$.", "href": "../../source/kernels/ADTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_Cu
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [coupled_time_Cu]
    type = ScaledCoupledTimeDerivative<<<{"description": "Time derivative Kernel that acts on a coupled variable. Weak form: $(\\psi_i, \\frac{\\partial v_h}{\\partial t})$.", "href": "../../source/kernels/ScaledCoupledTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_Cu
    v<<<{"description": "Coupled variable"}>>> = C_trapped_Cu
    factor<<<{"description": "The factor by which to scale"}>>> = 1e0
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [heat_conduction_Cu]
    type = HeatConduction<<<{"description": "Diffusive heat conduction term $-\\nabla\\cdot(k\\nabla T)$ of the thermal energy conservation equation", "href": "../../source/kernels/HeatConduction.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    diffusion_coefficient<<<{"description": "Property name of the diffusion coefficient"}>>> = thermal_conductivity_Cu
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [time_heat_conduction_Cu]
    type = SpecificHeatConductionTimeDerivative<<<{"description": "Time derivative term $\\rho c_p \\frac{\\partial T}{\\partial t}$ of the heat equation with the specific heat $c_p$ and the density $\\rho$ as arguments.", "href": "../../source/kernels/SpecificHeatConductionTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    specific_heat<<<{"description": "Property name of the specific heat material property"}>>> = specific_heat_Cu
    density<<<{"description": "Property name of the density material property"}>>> = density_Cu
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  ############################## Kernels for CuCrZr (block = 2)
  [diff_CuCrZr]
    type = ADMatDiffusion<<<{"description": "Diffusion equation kernel that takes an isotropic diffusivity from a material property", "href": "../../source/kernels/ADMatDiffusion.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_CuCrZr
    diffusivity<<<{"description": "The diffusivity value or material property"}>>> = diffusivity_CuCrZr
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [time_diff_CuCrZr]
    type = ADTimeDerivative<<<{"description": "The time derivative operator with the weak form of $(\\psi_i, \\frac{\\partial u_h}{\\partial t})$.", "href": "../../source/kernels/ADTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_CuCrZr
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [coupled_time_CuCrZr]
    type = ScaledCoupledTimeDerivative<<<{"description": "Time derivative Kernel that acts on a coupled variable. Weak form: $(\\psi_i, \\frac{\\partial v_h}{\\partial t})$.", "href": "../../source/kernels/ScaledCoupledTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_CuCrZr
    v<<<{"description": "Coupled variable"}>>> = C_trapped_CuCrZr
    factor<<<{"description": "The factor by which to scale"}>>> = 1e0
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [heat_conduction_CuCrZr]
    type = HeatConduction<<<{"description": "Diffusive heat conduction term $-\\nabla\\cdot(k\\nabla T)$ of the thermal energy conservation equation", "href": "../../source/kernels/HeatConduction.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    diffusion_coefficient<<<{"description": "Property name of the diffusion coefficient"}>>> = thermal_conductivity_CuCrZr
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [time_heat_conduction_CuCrZr]
    type = SpecificHeatConductionTimeDerivative<<<{"description": "Time derivative term $\\rho c_p \\frac{\\partial T}{\\partial t}$ of the heat equation with the specific heat $c_p$ and the density $\\rho$ as arguments.", "href": "../../source/kernels/SpecificHeatConductionTimeDerivative.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    specific_heat<<<{"description": "Property name of the specific heat material property"}>>> = specific_heat_CuCrZr
    density<<<{"description": "Property name of the density material property"}>>> = density_CuCrZr
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
[]

NodalKernels are used for non-diffusive variables (trapped species).

[NodalKernels<<<{"href": "../../syntax/NodalKernels/index.html"}>>>]
  ############################## NodalKernels for W (block = 4)
  [time_W]
    type = TimeDerivativeNodalKernel<<<{"description": "Forms the contribution to the residual and jacobian of the time derivative term from an ODE being solved at all nodes.", "href": "../../source/nodalkernels/TimeDerivativeNodalKernel.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_W
  []
  [trapping_W]
    type = TrappingNodalKernel<<<{"description": "Implements a residual describing the trapping of a species in a material.", "href": "../../source/nodal_kernels/TrappingNodalKernel.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_W
    temperature<<<{"description": "The temperature (K)"}>>> = temperature
    alpha_t<<<{"description": "The trapping rate coefficient. This has units of 1/s (e.g. no number densities are involved)"}>>> = 2.75e11 # 1e15
    N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = 1.0e0 # = (1e0) x (${tungsten_atomic_density} #/m^3)
    # Ct0 = 1.0e-4                # E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
    Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = 1.0e-4 # E.A. Hodille et al 2021 Nucl. Fusion 61 1260033, trap 2
    trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix (-)"}>>> = 1.0e0 # 1.0e1
    mobile_concentration<<<{"description": "The variable representing the mobile concentration of solute particles (1/m^3)"}>>> = 'C_mobile_W'
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [release_W]
    type = ReleasingNodalKernel<<<{"description": "Implements a residual describing the release of trapped species in a material.", "href": "../../source/nodal_kernels/ReleasingNodalKernel.html"}>>>
    alpha_r<<<{"description": "The release rate coefficient (1/s)"}>>> = 8.4e12 # 1.0e13
    temperature<<<{"description": "The temperature (K)"}>>> = temperature
    # detrapping_energy = 9863.9    # = 0.85 eV    E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
    detrapping_energy<<<{"description": "The detrapping energy (K)"}>>> = 11604.6 # = 1.00 eV    E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 2
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_W
  []
  ############################## NodalKernels for Cu (block = 3)
  [time_Cu]
    type = TimeDerivativeNodalKernel<<<{"description": "Forms the contribution to the residual and jacobian of the time derivative term from an ODE being solved at all nodes.", "href": "../../source/nodalkernels/TimeDerivativeNodalKernel.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_Cu
  []
  [trapping_Cu]
    type = TrappingNodalKernel<<<{"description": "Implements a residual describing the trapping of a species in a material.", "href": "../../source/nodal_kernels/TrappingNodalKernel.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_Cu
    temperature<<<{"description": "The temperature (K)"}>>> = temperature
    alpha_t<<<{"description": "The trapping rate coefficient. This has units of 1/s (e.g. no number densities are involved)"}>>> = 2.75e11 # 1e15
    N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = 1.0e0 # = ${tungsten_atomic_density} #/m^3 (W lattice density)
    Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = 5.0e-5 # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
    trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix (-)"}>>> = 1.0e0 # 1.0e1
    mobile_concentration<<<{"description": "The variable representing the mobile concentration of solute particles (1/m^3)"}>>> = 'C_mobile_Cu'
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [release_Cu]
    type = ReleasingNodalKernel<<<{"description": "Implements a residual describing the release of trapped species in a material.", "href": "../../source/nodal_kernels/ReleasingNodalKernel.html"}>>>
    alpha_r<<<{"description": "The release rate coefficient (1/s)"}>>> = 8.4e12 # 1.0e13
    temperature<<<{"description": "The temperature (K)"}>>> = temperature
    detrapping_energy<<<{"description": "The detrapping energy (K)"}>>> = 5802.3 # = 0.50eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_Cu
  []
  ############################## NodalKernels for CuCrZr (block = 2)
  [time_CuCrZr]
    type = TimeDerivativeNodalKernel<<<{"description": "Forms the contribution to the residual and jacobian of the time derivative term from an ODE being solved at all nodes.", "href": "../../source/nodalkernels/TimeDerivativeNodalKernel.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_CuCrZr
  []
  [trapping_CuCrZr]
    type = TrappingNodalKernel<<<{"description": "Implements a residual describing the trapping of a species in a material.", "href": "../../source/nodal_kernels/TrappingNodalKernel.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_CuCrZr
    temperature<<<{"description": "The temperature (K)"}>>> = temperature
    alpha_t<<<{"description": "The trapping rate coefficient. This has units of 1/s (e.g. no number densities are involved)"}>>> = 2.75e11 # 1e15
    N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = 1.0e0 # = ${tungsten_atomic_density} #/m^3 (W lattice density)
    Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = 5.0e-5 # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
    # Ct0 = 4.0e-2                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
    trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix (-)"}>>> = 1.0e0 # 1.0e1
    mobile_concentration<<<{"description": "The variable representing the mobile concentration of solute particles (1/m^3)"}>>> = 'C_mobile_CuCrZr'
    extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
  []
  [release_CuCrZr]
    type = ReleasingNodalKernel<<<{"description": "Implements a residual describing the release of trapped species in a material.", "href": "../../source/nodal_kernels/ReleasingNodalKernel.html"}>>>
    alpha_r<<<{"description": "The release rate coefficient (1/s)"}>>> = 8.4e12 # 1.0e13
    temperature<<<{"description": "The temperature (K)"}>>> = temperature
    detrapping_energy<<<{"description": "The detrapping energy (K)"}>>> = 5802.3 # = 0.50eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
    # detrapping_energy = 9631.8   # = 0.83 eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_CuCrZr
  []
[]

AuxVariables and AuxKernels

AuxVariables are used to track quantities that are not solved for by the PDEs, but are needed to compute materials properties, or are desirable to obtain as outputs, such as the total concentration of tritium or flux values. AuxKernels provide the expressions that define the AuxVariables.

[AuxVariables<<<{"href": "../../syntax/AuxVariables/index.html"}>>>]
  [flux_y]
    order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
    family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = MONOMIAL
  []
  ############################## AuxVariables for W (block = 4)
  [Sc_C_mobile_W]
    block = 4
  []
  [Sc_C_trapped_W]
    block = 4
  []
  [C_total_W]
    block = 4
  []
  [Sc_C_total_W]
    block = 4
  []
  [S_empty_W]
    block = 4
  []
  [Sc_S_empty_W]
    block = 4
  []
  [S_trapped_W]
    block = 4
  []
  [Sc_S_trapped_W]
    block = 4
  []
  [S_total_W]
    block = 4
  []
  [Sc_S_total_W]
    block = 4
  []
  ############################## AuxVariables for Cu (block = 3)
  [Sc_C_mobile_Cu]
    block = 3
  []
  [Sc_C_trapped_Cu]
    block = 3
  []
  [C_total_Cu]
    block = 3
  []
  [Sc_C_total_Cu]
    block = 3
  []
  [S_empty_Cu]
    block = 3
  []
  [Sc_S_empty_Cu]
    block = 3
  []
  [S_trapped_Cu]
    block = 3
  []
  [Sc_S_trapped_Cu]
    block = 3
  []
  [S_total_Cu]
    block = 3
  []
  [Sc_S_total_Cu]
    block = 3
  []
  ############################## AuxVariables for CuCrZr (block = 2)
  [Sc_C_mobile_CuCrZr]
    block = 2
  []
  [Sc_C_trapped_CuCrZr]
    block = 2
  []
  [C_total_CuCrZr]
    block = 2
  []
  [Sc_C_total_CuCrZr]
    block = 2
  []
  [S_empty_CuCrZr]
    block = 2
  []
  [Sc_S_empty_CuCrZr]
    block = 2
  []
  [S_trapped_CuCrZr]
    block = 2
  []
  [Sc_S_trapped_CuCrZr]
    block = 2
  []
  [S_total_CuCrZr]
    block = 2
  []
  [Sc_S_total_CuCrZr]
    block = 2
  []
[]

[AuxKernels<<<{"href": "../../syntax/AuxKernels/index.html"}>>>]
  ############################## AuxKernels for W (block = 4)
  [Scaled_mobile_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_mobile_W
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_mobile_W
  []
  [Scaled_trapped_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_trapped_W
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_W
  []
  [total_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = C_total_W
    type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
    expression<<<{"description": "Parsed function expression to compute"}>>> = 'C_mobile_W + C_trapped_W'
    coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'C_mobile_W C_trapped_W'
  []
  [Scaled_total_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_total_W
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_total_W
  []
  [empty_sites_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_empty_W
    type = EmptySitesAux<<<{"description": "Calculates the concentration of empty trapping sites.", "href": "../../source/auxkernels/EmptySitesAux.html"}>>>
    N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = '${units 1.0e0 m^-3}' # = ${tungsten_atomic_density} #/m^3 (W lattice density)
    # Ct0 = ${units 1.0e-4 m^-3}   # E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
    Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = '${units 1.0e-4 m^-3}' # E.A. Hodille et al 2021 Nucl. Fusion 61 1260033, trap 2
    trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix"}>>> = 1.0e0 # 1.0e1
    trapped_concentration_variables<<<{"description": "Variables representing trapped particle concentrations."}>>> = C_trapped_W
  []
  [scaled_empty_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_empty_W
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_empty_W
  []
  [trapped_sites_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_trapped_W
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = 1e0
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_W
  []
  [scaled_trapped_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_trapped_W
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_trapped_W
  []
  [total_sites_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_total_W
    type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
    expression<<<{"description": "Parsed function expression to compute"}>>> = 'S_trapped_W + S_empty_W'
    coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'S_trapped_W S_empty_W'
  []
  [scaled_total_W]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_total_W
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_total_W
  []
  ############################## AuxKernels for Cu (block = 3)
  [Scaled_mobile_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_mobile_Cu
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_mobile_Cu
  []
  [Scaled_trapped_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_trapped_Cu
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_Cu
  []
  [total_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = C_total_Cu
    type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
    expression<<<{"description": "Parsed function expression to compute"}>>> = 'C_mobile_Cu + C_trapped_Cu'
    coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'C_mobile_Cu C_trapped_Cu'
  []
  [Scaled_total_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_total_Cu
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_total_Cu
  []
  [empty_sites_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_empty_Cu
    type = EmptySitesAux<<<{"description": "Calculates the concentration of empty trapping sites.", "href": "../../source/auxkernels/EmptySitesAux.html"}>>>
    N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = '${units 1.0e0 m^-3}' # = ${tungsten_atomic_density} #/m^3 (W lattice density)
    Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = '${units 5.0e-5 m^-3}' # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
    trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix"}>>> = 1.0e0 # 1.0e1
    trapped_concentration_variables<<<{"description": "Variables representing trapped particle concentrations."}>>> = C_trapped_Cu
  []
  [scaled_empty_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_empty_Cu
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_empty_Cu
  []
  [trapped_sites_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_trapped_Cu
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = 1e0
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_Cu
  []
  [scaled_trapped_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_trapped_Cu
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_trapped_Cu
  []
  [total_sites_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_total_Cu
    type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
    expression<<<{"description": "Parsed function expression to compute"}>>> = 'S_trapped_Cu + S_empty_Cu'
    coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'S_trapped_Cu S_empty_Cu'
  []
  [scaled_total_Cu]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_total_Cu
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_total_Cu
  []
  ############################## AuxKernels for CuCrZr (block = 2)
  [Scaled_mobile_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_mobile_CuCrZr
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_mobile_CuCrZr
  []
  [Scaled_trapped_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_trapped_CuCrZr
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_CuCrZr
  []
  [total_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = C_total_CuCrZr
    type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
    expression<<<{"description": "Parsed function expression to compute"}>>> = 'C_mobile_CuCrZr + C_trapped_CuCrZr'
    coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'C_mobile_CuCrZr C_trapped_CuCrZr'
  []
  [Scaled_total_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_total_CuCrZr
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_total_CuCrZr
  []
  [empty_sites_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_empty_CuCrZr
    type = EmptySitesAux<<<{"description": "Calculates the concentration of empty trapping sites.", "href": "../../source/auxkernels/EmptySitesAux.html"}>>>
    N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = '${units 1.0e0 m^-3}' # = ${tungsten_atomic_density} #/m^3 (W lattice density)
    Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = '${units 5.0e-5 m^-3}' # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
    # Ct0 = ${units 4.0e-2 m^-3} # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
    trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix"}>>> = 1.0e0 # 1.0e1
    trapped_concentration_variables<<<{"description": "Variables representing trapped particle concentrations."}>>> = C_trapped_CuCrZr
  []
  [scaled_empty_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_empty_CuCrZr
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_empty_CuCrZr
  []
  [trapped_sites_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_trapped_CuCrZr
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = 1e0
    source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_CuCrZr
  []
  [scaled_trapped_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_trapped_CuCrZr
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_trapped_CuCrZr
  []
  [total_sites_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_total_CuCrZr
    type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
    expression<<<{"description": "Parsed function expression to compute"}>>> = 'S_trapped_CuCrZr + S_empty_CuCrZr'
    coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'S_trapped_CuCrZr S_empty_CuCrZr'
  []
  [scaled_total_CuCrZr]
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_total_CuCrZr
    type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
    normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
    source_variable<<<{"description": "The variable to be normalized"}>>> = S_total_CuCrZr
  []
  [flux_y_W]
    type = DiffusionFluxAux<<<{"description": "Compute components of flux vector for diffusion problems $(\\vec{J} = -D \\nabla C)$.", "href": "../../source/auxkernels/DiffusionFluxAux.html"}>>>
    diffusivity<<<{"description": "The name of the diffusivity material property that will be used in the flux computation."}>>> = diffusivity_W
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = flux_y
    diffusion_variable<<<{"description": "The name of the variable"}>>> = C_mobile_W
    component<<<{"description": "The desired component of flux."}>>> = y
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
  []
  [flux_y_Cu]
    type = DiffusionFluxAux<<<{"description": "Compute components of flux vector for diffusion problems $(\\vec{J} = -D \\nabla C)$.", "href": "../../source/auxkernels/DiffusionFluxAux.html"}>>>
    diffusivity<<<{"description": "The name of the diffusivity material property that will be used in the flux computation."}>>> = diffusivity_Cu
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = flux_y
    diffusion_variable<<<{"description": "The name of the variable"}>>> = C_mobile_Cu
    component<<<{"description": "The desired component of flux."}>>> = y
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
  []
  [flux_y_CuCrZr]
    type = DiffusionFluxAux<<<{"description": "Compute components of flux vector for diffusion problems $(\\vec{J} = -D \\nabla C)$.", "href": "../../source/auxkernels/DiffusionFluxAux.html"}>>>
    diffusivity<<<{"description": "The name of the diffusivity material property that will be used in the flux computation."}>>> = diffusivity_CuCrZr
    variable<<<{"description": "The name of the variable that this object applies to"}>>> = flux_y
    diffusion_variable<<<{"description": "The name of the variable"}>>> = C_mobile_CuCrZr
    component<<<{"description": "The desired component of flux."}>>> = y
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
  []
[]

Boundary conditions and history

For the mass conservation equation, we use FunctionNeumannBC and DirichletBC to solve the 2 $var element = document.getElementById("moose-equation-6c04065a-741a-496d-86eb-2f0f081816cd");katex.render("^{\\text{nd}}", element, {displayMode:false,throwOnError:false});$ term of Eq. (3). FunctionNeumannBC is used to treat the time-dependent plasma exposure at the top (plasma-exposed) surface, and DirichletBC is used to set the BC of the solute T atom concentration to zero at the inner CuCrZr tube (at $var element = document.getElementById("moose-equation-786763b3-a871-4416-9449-b2d79b2ea7e2");katex.render("r = 6", element, {displayMode:false,throwOnError:false});$ mm). For the energy conservation equation, we use FunctionNeumannBC and DirichletBC to solve the 2 $var element = document.getElementById("moose-equation-89e5c816-906d-4151-9c48-85a3f32aec5b");katex.render("^{\\text{nd}}", element, {displayMode:false,throwOnError:false});$ term of Eq. (4). FunctionNeumannBC is used to treat the time-dependent heat flux at the top (plasma-exposed) surface, and FunctionDirichletBC is used to treat the temperature increase in the cooling tube. FunctionNeumannBC, FunctionDirichletBC, DirichletBC, and FunctionDirichletBC are MOOSE objects commonly used to represent the BCs of the variables to be solved.

We simulate a 20,000-second plasma discharge, with each 1,600-second cycle consisting of a 100-second plasma ramp-up, a 400-second steady-state plasma discharge, a 100-second plasma ramp-down, and a 1,000-second waiting phase. Up to 50 cycles are simulated to achieve the total discharge. Figure 2 shows the integrated (solute, total and trapped) tritium concentration profiles in the monoblock. It shows that implantation fluxes and temperatures vary linearly up to (from) their steady-state values from (up to) their initial values during ramp-up (ramp-down).

Figure 2: Temperature profiles (orange) and integrated tritium concentration profiles (blue) during two 1,600-second-cycle plasma discharges. This corresponds to Fig. 2 in Shimada et al. (2024).

During the steady-state plasma discharge, we set a heat flux of 10 MW $var element = document.getElementById("moose-equation-5e8506ee-3642-4acd-a499-a9078d770879");katex.render("\\cdot", element, {displayMode:false,throwOnError:false});$ m $var element = document.getElementById("moose-equation-6d6ff5a0-5e28-4941-9aaf-3ef4f2313abe");katex.render("^{-2}", element, {displayMode:false,throwOnError:false});$ at the top of the 2D monoblock (at $var element = document.getElementById("moose-equation-73d10501-1a72-4214-b623-969e6a28fb9b");katex.render("y = 14.0", element, {displayMode:false,throwOnError:false});$ mm) and a cooling temperature of 552 K at the inner CuCrZr tube (at $var element = document.getElementById("moose-equation-5222d216-134f-476c-bed0-dfc146b3b3a5");katex.render("r = 6.0", element, {displayMode:false,throwOnError:false});$ mm). We assume a 100% T plasma with a 5.0 $var element = document.getElementById("moose-equation-bbb2df48-2f80-4ec4-b71d-919f0ca88a6d");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-dbad83af-efe8-4eee-9a06-3c9ecd5adcb0");katex.render("^{23}", element, {displayMode:false,throwOnError:false});$ m $var element = document.getElementById("moose-equation-218c99b1-ea49-4ad3-b62d-258706ee075d");katex.render("^{-2}", element, {displayMode:false,throwOnError:false});$ $var element = document.getElementById("moose-equation-deb9c97e-7f53-41fa-86c4-47c548f84cad");katex.render("\\cdot", element, {displayMode:false,throwOnError:false});$ s $var element = document.getElementById("moose-equation-834d8781-dd0f-4701-b21b-2b26f9e89289");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ plasma particle flux (which is half of the full 1.0 $var element = document.getElementById("moose-equation-10f6dd4f-5370-4fa9-9db1-5f13ca3a56cc");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-ebf5ce72-76dd-4c89-a182-1d1d7d7b200a");katex.render("^{24}", element, {displayMode:false,throwOnError:false});$ m $var element = document.getElementById("moose-equation-6142bcc5-e3b5-4de0-9bb3-a1f31c82c11b");katex.render("^{-2}", element, {displayMode:false,throwOnError:false});$ $var element = document.getElementById("moose-equation-1a768355-63e8-41c7-b5ab-990a9352f537");katex.render("\\cdot", element, {displayMode:false,throwOnError:false});$ s $var element = document.getElementById("moose-equation-8ec68c25-1876-478f-a170-08f65e0d3cd4");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ DT plasma particle flux), and only 0.1% of the incident plasma particle flux (5.0 $var element = document.getElementById("moose-equation-af7ea11f-f339-4636-a75e-7dde6f9810c8");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-52e41399-92d8-4c91-86c7-5aaae2091a69");katex.render("^{20}", element, {displayMode:false,throwOnError:false});$ m $var element = document.getElementById("moose-equation-224d5cb3-e529-4b1d-918a-e404e065b2c8");katex.render("^{-2}", element, {displayMode:false,throwOnError:false});$ $var element = document.getElementById("moose-equation-a91bada7-80e6-4f6a-a5d6-0771e9cc2ece");katex.render("\\cdot", element, {displayMode:false,throwOnError:false});$ s $var element = document.getElementById("moose-equation-32e34798-25a0-4ba2-abae-26565a33c47c");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ ) entered the first layer of mesh at the exposed surface ( $var element = document.getElementById("moose-equation-a6d91683-099d-4607-8e1f-2b1a54139b9d");katex.render("y = 14.0", element, {displayMode:false,throwOnError:false});$ mm) as in Shimada et al. (2024). The solute T atom concentration is set to zero at the inner CuCrZr tube (at $var element = document.getElementById("moose-equation-f4030573-3744-4084-bd68-bdf276fe165b");katex.render("r = 6.0", element, {displayMode:false,throwOnError:false});$ mm).

We treated this plasma exposure by setting the flux BC of the solute T atom concentration at the exposed surface ( $var element = document.getElementById("moose-equation-52495c0b-6c5f-4180-8431-b667fb3e9a06");katex.render("y = 14.0", element, {displayMode:false,throwOnError:false});$ mm) as a simplification of the complex plasma implantation and recombination phenomena, which would require a very fine mesh and increase computational costs.

[BCs<<<{"href": "../../syntax/BCs/index.html"}>>>]
  [C_mob_W_top_flux]
    type = FunctionNeumannBC<<<{"description": "Imposes the integrated boundary condition $\\frac{\\partial u}{\\partial n}=h(t,\\vec{x})$, where $h$ is a (possibly) time and space-dependent MOOSE Function.", "href": "../../source/bcs/FunctionNeumannBC.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = 'top'
    function<<<{"description": "The function."}>>> = mobile_flux_bc_func
  []
  [mobile_tube]
    type = DirichletBC<<<{"description": "Imposes the essential boundary condition $u=g$, where $g$ is a constant, controllable value.", "href": "../../source/bcs/DirichletBC.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_CuCrZr
    value<<<{"description": "Value of the BC"}>>> = 1.0e-18
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = '2to1'
  []
  [temp_top]
    type = FunctionNeumannBC<<<{"description": "Imposes the integrated boundary condition $\\frac{\\partial u}{\\partial n}=h(t,\\vec{x})$, where $h$ is a (possibly) time and space-dependent MOOSE Function.", "href": "../../source/bcs/FunctionNeumannBC.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = 'top'
    function<<<{"description": "The function."}>>> = temp_flux_bc_func
  []
  [temp_tube]
    type = FunctionDirichletBC<<<{"description": "Imposes the essential boundary condition $u=g(t,\\vec{x})$, where $g$ is a (possibly) time and space-dependent MOOSE Function.", "href": "../../source/bcs/FunctionDirichletBC.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
    boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = '2to1'
    function<<<{"description": "The forcing function."}>>> = temp_inner_func
  []
[]

Functions

In this example, Functions are used to define the time-dependent tritium and heat flux at the exposed surface, as well as the coolant temperature. Functions can also be spatially-dependent. The functions are then provided to BCs, as described above.

[Functions<<<{"href": "../../syntax/Functions/index.html"}>>>]
  ### Maximum mobile flux of 7.90e-13 at the top surface (1.0e-4 [m])
  ### 1.80e23/m^2-s = (5.0e23/m^2-s) *(1-0.999) = (7.90e-13)*(${tungsten_atomic_density})/(1.0e-4)  at steady state
  [mobile_flux_bc_func]
    type = ParsedFunction<<<{"description": "Function created by parsing a string", "href": "../../source/functions/MooseParsedFunction.html"}>>>
    expression<<<{"description": "The user defined function."}>>> = 'if((t % 1600) < 100.0, 0.0      + 7.90e-13*(t % 1600)/100,
                        if((t % 1600) < 500.0, 7.90e-13,
                        if((t % 1600) < 600.0, 7.90e-13 - 7.90e-13*((t % 1600)-500)/100, 0.0)))'
  []
  ### Heat flux of 10MW/m^2 at steady state
  [temp_flux_bc_func]
    type = ParsedFunction<<<{"description": "Function created by parsing a string", "href": "../../source/functions/MooseParsedFunction.html"}>>>
    expression<<<{"description": "The user defined function."}>>> = 'if((t % 1600) < 100.0, 0.0   + 1.0e7*(t % 1600)/100,
                        if((t % 1600) < 500.0, 1.0e7,
                        if((t % 1600) < 600.0, 1.0e7 - 1.0e7*((t % 1600)-500)/100, 300)))'
  []
  ### Maximum coolant temperature of 552K at steady state
  [temp_inner_func]
    type = ParsedFunction<<<{"description": "Function created by parsing a string", "href": "../../source/functions/MooseParsedFunction.html"}>>>
    expression<<<{"description": "The user defined function."}>>> = 'if((t % 1600) < 100.0, 300.0 + (552-300)*(t % 1600)/100,
                        if((t % 1600) < 500.0, 552,
                        if((t % 1600) < 600.0, 552.0 - (552-300)*((t % 1600)-500)/100, 300)))'
  []
  [timestep_function]
    type = ParsedFunction<<<{"description": "Function created by parsing a string", "href": "../../source/functions/MooseParsedFunction.html"}>>>
    expression<<<{"description": "The user defined function."}>>> = 'if((t % 1600) <   10.0,  20,
                      if((t % 1600) <   90.0,  40,
                      if((t % 1600) <  110.0,  20,
                      if((t % 1600) <  480.0,  40,
                      if((t % 1600) <  500.0,  20,
                      if((t % 1600) <  590.0,   4,
                      if((t % 1600) <  610.0,  20,
                      if((t % 1600) < 1500.0, 200,
                      if((t % 1600) < 1600.0,  40,  2)))))))))'
  []
[]

Material properties

We use the tritium mass transport properties listed in Table 2 to solve the mass conservation equation (Eq. (3)) for two variables: $var element = document.getElementById("moose-equation-9a2c56fb-ff2b-4bfa-b83e-7fe442d8b15d");katex.render("C_s", element, {displayMode:false,throwOnError:false});$ and $var element = document.getElementById("moose-equation-5e55de5d-b3ca-4c23-9386-7ed369668473");katex.render("C_t", element, {displayMode:false,throwOnError:false});$ (i.e., solute and trapped T atom concentrations) in three materials and the thermal properties (i.e., density, temperature-dependent specific heat, and thermal conductivity) listed in Table 3 to solve the energy conservation equation (Eq. (4)) for one variable: the temperature $var element = document.getElementById("moose-equation-a48fe2c2-e664-4262-8a32-31a6d941bb8e");katex.render("T", element, {displayMode:false,throwOnError:false});$ . The diffusivity is defined as $var element = document.getElementById("moose-equation-2d8a3524-6dc0-4485-8ce9-23f6142f00ff");katex.render("D=D_0 \\exp\u2061(-E_D/k_B/T)", element, {displayMode:false,throwOnError:false});$ and the solubility is defined as $var element = document.getElementById("moose-equation-b288c5e9-6231-4d31-ae03-b2b9cf94fa47");katex.render("K_s=K_{s,0} \\exp\u2061(-E_s/k_B/T)", element, {displayMode:false,throwOnError:false});$ .

Table 2: Tritium mass transport and trapping properties used in the W, Cu, and CuCrZr layers. $var element = document.getElementById("moose-equation-a4e21223-d9f4-4db2-9110-abd59c8994ec");katex.render("^{**}", element, {displayMode:false,throwOnError:false});$ NOTE: Two-component solubility in W kept the maximum solubility ratio between W and Cu to 104 at low temperature. The references are provided in Shimada et al. (2024).

Material	$var element = document.getElementById("moose-equation-5f8b582d-79d8-4d86-bd70-7f9603a96f09");katex.render("D_0", element, {displayMode:false,throwOnError:false});$ (m $var element = document.getElementById("moose-equation-fc94306f-41fc-47a4-a76e-7de1601132f8");katex.render("^2", element, {displayMode:false,throwOnError:false});$ /s)	$var element = document.getElementById("moose-equation-fbefc6b9-871b-4af4-a4bf-66d7f22e0951");katex.render("E_D", element, {displayMode:false,throwOnError:false});$ (eV)	$var element = document.getElementById("moose-equation-ecef136f-b1f6-4032-bec6-139de1dbafd4");katex.render("K_{s,0}", element, {displayMode:false,throwOnError:false});$ (Pa $var element = document.getElementById("moose-equation-bee71dc7-7248-4f78-ac61-b4045c6d5db9");katex.render("^{1/2}", element, {displayMode:false,throwOnError:false});$ )	$var element = document.getElementById("moose-equation-b3f593b1-c4a5-44c4-b695-a634a5ecffff");katex.render("E_s", element, {displayMode:false,throwOnError:false});$ (eV)	Detrapping energy: $var element = document.getElementById("moose-equation-8ecb2d61-e21a-4f64-ae34-3370ab289e4b");katex.render("E_{dt}", element, {displayMode:false,throwOnError:false});$ (eV)	Trap density: $var element = document.getElementById("moose-equation-4f61512e-a5e3-449d-9d2c-52d9eabbd181");katex.render("n_{trap}", element, {displayMode:false,throwOnError:false});$ (at.fr.)
W	2.4 $var element = document.getElementById("moose-equation-74772804-33db-4564-8ba2-07e704df608c");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-531376db-5a11-4096-80a1-6682ef5c1db1");katex.render("^{-7}", element, {displayMode:false,throwOnError:false});$	0.39	1.87 $var element = document.getElementById("moose-equation-28a7179d-a6ed-4107-a3af-a2ef1c8b3d24");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-35dacdd6-9dd5-4135-9d15-3879922ff013");katex.render("^{24}", element, {displayMode:false,throwOnError:false});$	1.04	0.85	1.0 $var element = document.getElementById("moose-equation-91d78e2b-6fac-4016-be80-a7212995cb13");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-95ab544b-36d3-4a01-8d90-04c11ee22bc3");katex.render("^{-4}", element, {displayMode:false,throwOnError:false});$
W			3.14 $var element = document.getElementById("moose-equation-5cf899fb-5df7-48cc-8c84-3addd8c993c5");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-6adca820-3c7f-4d39-9eb6-9e5720db306e");katex.render("^{20}", element, {displayMode:false,throwOnError:false});$ $var element = document.getElementById("moose-equation-bd264235-3eae-426a-a013-468ffcdd2cd2");katex.render("^{**}", element, {displayMode:false,throwOnError:false});$	0.57 $var element = document.getElementById("moose-equation-b687635f-887d-4ea7-867e-9e716cddd2de");katex.render("^{**}", element, {displayMode:false,throwOnError:false});$
Cu	6.6 $var element = document.getElementById("moose-equation-f0922e5f-cae2-414d-ae09-2952800fafd0");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-96e262c5-e9a8-4341-be64-5e00bc6f4293");katex.render("^{-7}", element, {displayMode:false,throwOnError:false});$	0.39	3.14 $var element = document.getElementById("moose-equation-5e3a52bb-2615-4715-b477-8abe84bb4011");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-35397739-bc35-4bde-a579-20e431483bf8");katex.render("^{24}", element, {displayMode:false,throwOnError:false});$	0.57	0.50	5.0 $var element = document.getElementById("moose-equation-8623ccd9-5c55-463b-b322-565fa6b533c1");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-b7a78ea7-075a-4a55-b3f9-4f73dc4041dd");katex.render("^{-5}", element, {displayMode:false,throwOnError:false});$
CuCrZr	3.9 $var element = document.getElementById("moose-equation-67112c6f-be2f-4576-8170-6c327347a384");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-c5b4c246-f756-4430-94b9-522171a4cf37");katex.render("^{-7}", element, {displayMode:false,throwOnError:false});$	0.42	4.28 $var element = document.getElementById("moose-equation-71d4915b-ae86-404c-b42e-8662c4a9e684");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-baee0d76-3d72-453c-9481-f77e72eea4fd");katex.render("^{23}", element, {displayMode:false,throwOnError:false});$	0.39	0.83	5.0 $var element = document.getElementById("moose-equation-1f013347-7f36-4ec3-98be-cf996c7681b8");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-0e133723-e8ad-451d-9bba-e90b43f5c72c");katex.render("^{-5}", element, {displayMode:false,throwOnError:false});$

Table 3: Thermal properties used in the W, Cu, and CuCrZr layers. The references are provided in Shimada et al. (2024).

Material	Density: $var element = document.getElementById("moose-equation-f43a4c87-759d-49dc-a0aa-450db7e673e5");katex.render("\\rho", element, {displayMode:false,throwOnError:false});$ (g m $var element = document.getElementById("moose-equation-07c04165-3710-4e31-a4fa-972186f41aba");katex.render("^{-3}", element, {displayMode:false,throwOnError:false});$ )	Specific heat: $var element = document.getElementById("moose-equation-3a6a9bf7-28eb-4f24-89cb-a96703ce619e");katex.render("c_p", element, {displayMode:false,throwOnError:false});$ (J kg $var element = document.getElementById("moose-equation-8a7c0f0f-31d1-462a-bb5b-b5f3fff77999");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ K $var element = document.getElementById("moose-equation-3a054376-28a0-4b95-a290-6fb909263cfb");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ )	Thermal conductivity: $var element = document.getElementById("moose-equation-fb264a49-9b43-4a51-893d-70845eacceab");katex.render("k_T", element, {displayMode:false,throwOnError:false});$ (W m $var element = document.getElementById("moose-equation-ff26168c-06bc-42cd-81d9-22a5ccd2e8fe");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ K $var element = document.getElementById("moose-equation-18f1040c-164f-4a0f-93fc-64010db17567");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ )
W	19,300	1.16 $var element = document.getElementById("moose-equation-64f9712a-b781-4883-990a-4cc2f052d8fa");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-1493a743-e3f3-483a-8951-3a6f090ff21c");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ +7.11 $var element = document.getElementById("moose-equation-6788af43-a50b-4a3a-badd-abd6d94fe931");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-c81fd956-4625-4cd2-bcb0-6245b0aa93e6");katex.render("^{-2}", element, {displayMode:false,throwOnError:false});$ T –6.58 $var element = document.getElementById("moose-equation-6883971e-ea55-4a3a-b684-b28c81adcf0e");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-1fbeee9b-a9a3-4658-8ac1-60b596ff17de");katex.render("^{-5}", element, {displayMode:false,throwOnError:false});$ T $var element = document.getElementById("moose-equation-325281d1-34b6-47fb-8d7a-2efb41342cf8");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ +3.24 $var element = document.getElementById("moose-equation-08d42983-fd05-4e66-b2f5-ee18ff205f18");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-6d75d502-e90d-4863-aaa8-a45298fbe0e8");katex.render("^{-8}", element, {displayMode:false,throwOnError:false});$ T $var element = document.getElementById("moose-equation-c5d944c3-b242-4ae9-ae4a-2b960727eb08");katex.render("^{3}", element, {displayMode:false,throwOnError:false});$ –5.45 $var element = document.getElementById("moose-equation-c388dad6-3be3-460b-97eb-70e457cb2b23");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-08551a10-9cd1-47de-92f7-18cf9c606ba0");katex.render("^{-12}", element, {displayMode:false,throwOnError:false});$ T $var element = document.getElementById("moose-equation-09d95500-ceb3-498c-b9d3-6b471863a8f9");katex.render("^{4}", element, {displayMode:false,throwOnError:false});$ (293 < T (K) < 2500)	2.41 $var element = document.getElementById("moose-equation-7492c03b-a77a-4252-b75e-591aa6955f9f");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-922c66c8-630b-4d4e-ac3e-818960c73ae4");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ –2.90 $var element = document.getElementById("moose-equation-12db8a46-0881-4e74-9d2f-8cbb5f412995");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-6d8ae9cd-11dc-4077-a24e-14c65d954d6f");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ T + 2.54 $var element = document.getElementById("moose-equation-12eb0ba4-bc9f-42ce-94bc-c497188c57fd");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-7c1c0fb0-67b9-43a5-b2f9-e587f6856416");katex.render("^{-4}", element, {displayMode:false,throwOnError:false});$ T $var element = document.getElementById("moose-equation-947ca239-d900-4ba3-96f5-d597aaed80fb");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ –1.03 $var element = document.getElementById("moose-equation-6f347040-7e24-4e01-9596-73b0105ef230");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-cf113fdb-7e1c-4bfb-a941-d300c1f35ad1");katex.render("^{-7}", element, {displayMode:false,throwOnError:false});$ T $var element = document.getElementById("moose-equation-bf065886-576d-4ce4-8490-70ec4299b1be");katex.render("^{3}", element, {displayMode:false,throwOnError:false});$ +1.52 $var element = document.getElementById("moose-equation-966217c8-1ecc-49e7-ae37-62b3edfb35a5");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-fa74a887-dcd4-477c-bc30-1501bcd89840");katex.render("^{-11}", element, {displayMode:false,throwOnError:false});$ T $var element = document.getElementById("moose-equation-13735bf6-b1bc-4b00-afeb-b60adc0970c0");katex.render("^{4}", element, {displayMode:false,throwOnError:false});$ (293 < T (K) < 2500)
Cu	8,960	4.21 $var element = document.getElementById("moose-equation-de9f1ca1-bce3-4ac6-b1a8-3a64dcb115da");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-db7a2fec-6fd4-48cc-be76-d7469190634c");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ –6.85 $var element = document.getElementById("moose-equation-52764646-b1ff-4d41-a1d8-ee997fd98a23");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-51d8bd42-4acc-4fb8-9a42-09dac8e5bf31");katex.render("^{-2}", element, {displayMode:false,throwOnError:false});$ T (293 < T (K) < 873)	3.16 $var element = document.getElementById("moose-equation-0a9993cc-8298-416c-ad8b-fc1ad0b30659");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-09a68756-0211-4567-9ca7-09714356c480");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ +3.18 $var element = document.getElementById("moose-equation-26711bf8-bc7d-4109-b05d-630f37935888");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-dc54078c-9013-4d21-90a0-95315eaf2012");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ T –3.49 $var element = document.getElementById("moose-equation-79571363-624c-47a0-b26b-f548fa305bb6");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-0bfcfff5-b40c-47c7-a2a4-c46910721622");katex.render("^{-4}", element, {displayMode:false,throwOnError:false});$ T $var element = document.getElementById("moose-equation-3fd30e70-714e-4cdb-954a-66605af52db7");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ +1.66 $var element = document.getElementById("moose-equation-f07aedaf-98b0-495b-91ef-7c2f0e467c23");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-cafa6a67-ae08-4509-94d2-d5f4a7dab214");katex.render("^{-7}", element, {displayMode:false,throwOnError:false});$ T $var element = document.getElementById("moose-equation-7e3fb093-0e5f-40b8-9cda-3c2676e2051c");katex.render("^{3}", element, {displayMode:false,throwOnError:false});$ (293 < T (K) < 873)
CuCrZr	8,900	390	3.87 $var element = document.getElementById("moose-equation-c2a4e3ea-fe05-4003-a04d-baddd9a07175");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-ddb68b9e-68d4-4f52-9b06-e8565c33ec4a");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ –1.28 $var element = document.getElementById("moose-equation-b5de72eb-003c-42c0-b43b-72d7ef42bb2b");katex.render("\\times", element, {displayMode:false,throwOnError:false});$ 10 $var element = document.getElementById("moose-equation-b6aa3d62-3b87-4e41-ae90-7c5f92a97357");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ T (293 < T (K) < 927)

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
  ############################## Materials for W (block = 4)
  [diffusivity_W]
    type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = diffusivity_W
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '2.4e-7*exp(-4525.8/temperature)' # H diffusivity in W
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  [solubility_W]
    type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = solubility_W
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    # expression = '2.95e-5 *exp(-12069.0/temperature)'              # H solubility in W = (1.87e24)/(${tungsten_atomic_density}) [#/m^3]
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '2.95e-5 *exp(-12069.0/temperature) + 4.95e-8 * exp(-6614.6/temperature)' # H solubility in W = (1.87e24)/(${tungsten_atomic_density}) [#/m^3]
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  [converter_to_regular_W]
    type = MaterialADConverter<<<{"description": "Converts regular material properties to AD properties and vice versa", "href": "../../source/materials/MaterialADConverter.html"}>>>
    ad_props_in<<<{"description": "The names of the AD material properties to convert to regular properties"}>>> = 'diffusivity_W'
    reg_props_out<<<{"description": "The names of the output regular properties"}>>> = 'diffusivity_W_nonAD'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
  []
  [heat_transfer_W]
    type = GenericConstantMaterial<<<{"description": "Declares material properties based on names and values prescribed by input parameters.", "href": "../../source/materials/GenericConstantMaterial.html"}>>>
    prop_names<<<{"description": "The names of the properties this material will have"}>>> = 'density_W'
    prop_values<<<{"description": "The values associated with the named properties"}>>> = '19300' # [g/m^3]
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
  []
  [specific_heat_W]
    type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = specific_heat_W
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '1.16e2 + 7.11e-2 * temperature - 6.58e-5 * temperature^2 + 3.24e-8 * temperature^3 -5.45e-12 * temperature^4' # ~ 132[J/kg-K]
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  [thermal_conductivity_W]
    type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = thermal_conductivity_W
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    # expression = '-7.8e-9 * temperature^3 + 5.0e-5 * temperature^2 - 1.1e-1 * temperature + 1.8e2'    # ~ 173.0 [ W/m-K]   from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '2.41e2 - 2.90e-1 * temperature + 2.54e-4 * temperature^2 - 1.03e-7 * temperature^3 + 1.52e-11 * temperature^4' # ~ 173.0 [ W/m-K]
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  ############################## Materials for Cu (block = 3)
  [diffusivity_Cu]
    type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = diffusivity_Cu
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '6.60e-7*exp(-4525.8/temperature)' # H diffusivity in Cu
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  [solubility_Cu]
    type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = solubility_Cu
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '4.95e-5*exp(-6614.6/temperature)' # H solubility in Cu = (3.14e24)/(${tungsten_atomic_density}) [#/m^3]
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  [converter_to_regular_Cu]
    type = MaterialADConverter<<<{"description": "Converts regular material properties to AD properties and vice versa", "href": "../../source/materials/MaterialADConverter.html"}>>>
    ad_props_in<<<{"description": "The names of the AD material properties to convert to regular properties"}>>> = 'diffusivity_Cu'
    reg_props_out<<<{"description": "The names of the output regular properties"}>>> = 'diffusivity_Cu_nonAD'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
  []
  [heat_transfer_Cu]
    type = GenericConstantMaterial<<<{"description": "Declares material properties based on names and values prescribed by input parameters.", "href": "../../source/materials/GenericConstantMaterial.html"}>>>
    prop_names<<<{"description": "The names of the properties this material will have"}>>> = 'density_Cu'
    prop_values<<<{"description": "The values associated with the named properties"}>>> = '8960.0' # [g/m^3]
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
  []
  [specific_heat_Cu]
    type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = specific_heat_Cu
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '3.16e2 + 3.18e-1 * temperature - 3.49e-4 * temperature^2 + 1.66e-7 * temperature^3' # ~ 384 [J/kg-K]
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  [thermal_conductivity_Cu]
    type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = thermal_conductivity_Cu
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    # expression = '-3.9e-8 * temperature^3 + 3.8e-5 * temperature^2 - 7.9e-2 * temperature + 4.0e2'    # ~ 401.0  [ W/m-K] from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '4.21e2 - 6.85e-2 * temperature' # ~ 400.0 [ W/m-K]
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  ############################## Materials for CuCrZr (block = 2)
  [diffusivity_CuCrZr]
    type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = diffusivity_CuCrZr
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '3.90e-7*exp(-4873.9/temperature)' # H diffusivity in CuCrZr
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  [solubility_CuCrZr]
    type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = solubility_CuCrZr
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '6.75e-6*exp(-4525.8/temperature)' # H solubility in CuCrZr = (4.28e23)/(${tungsten_atomic_density}) [#/m^3]
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  [converter_to_regular_CuCrZr]
    type = MaterialADConverter<<<{"description": "Converts regular material properties to AD properties and vice versa", "href": "../../source/materials/MaterialADConverter.html"}>>>
    ad_props_in<<<{"description": "The names of the AD material properties to convert to regular properties"}>>> = 'diffusivity_CuCrZr'
    reg_props_out<<<{"description": "The names of the output regular properties"}>>> = 'diffusivity_CuCrZr_nonAD'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
  []
  [heat_transfer_CuCrZr]
    type = GenericConstantMaterial<<<{"description": "Declares material properties based on names and values prescribed by input parameters.", "href": "../../source/materials/GenericConstantMaterial.html"}>>>
    prop_names<<<{"description": "The names of the properties this material will have"}>>> = 'density_CuCrZr specific_heat_CuCrZr'
    prop_values<<<{"description": "The values associated with the named properties"}>>> = '8900.0 390.0' # [g/m^3], [ W/m-K], [J/kg-K]
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
  []
  [thermal_conductivity_CuCrZr]
    type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
    property_name<<<{"description": "Name of the parsed material property"}>>> = thermal_conductivity_CuCrZr
    coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    # expression = '5.3e-7 * temperature^3 - 6.5e-4 * temperature^2 + 2.6e-1 * temperature + 3.1e2'    # ~ 320.0  [ W/m-K] from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
    expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '3.87e2 - 1.28e-1 * temperature' # ~ 349 [ W/m-K]
    outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
  []
  ############################## Materials for others
  [interface_jump_4to3]
    type = SolubilityRatioMaterial<<<{"description": "Calculates the jump in concentration across an interface.", "href": "../../source/materials/SolubilityRatioMaterial.html"}>>>
    solubility_primary<<<{"description": "The material property on the primary side of the interface"}>>> = solubility_W
    solubility_secondary<<<{"description": "The material property on the secondary side of the interface"}>>> = solubility_Cu
    boundary<<<{"description": "The list of boundaries (ids or names) from the mesh where this object applies"}>>> = '4to3'
    concentration_primary<<<{"description": "Primary side non-linear variable for jump computation"}>>> = C_mobile_W
    concentration_secondary<<<{"description": "Secondary side non-linear variable for jump computation"}>>> = C_mobile_Cu
  []
  [interface_jump_3to2]
    type = SolubilityRatioMaterial<<<{"description": "Calculates the jump in concentration across an interface.", "href": "../../source/materials/SolubilityRatioMaterial.html"}>>>
    solubility_primary<<<{"description": "The material property on the primary side of the interface"}>>> = solubility_Cu
    solubility_secondary<<<{"description": "The material property on the secondary side of the interface"}>>> = solubility_CuCrZr
    boundary<<<{"description": "The list of boundaries (ids or names) from the mesh where this object applies"}>>> = '3to2'
    concentration_primary<<<{"description": "Primary side non-linear variable for jump computation"}>>> = C_mobile_Cu
    concentration_secondary<<<{"description": "Secondary side non-linear variable for jump computation"}>>> = C_mobile_CuCrZr
  []
[]

Interface Kernels

TMAP8 assumes equilibrium between the chemical potentials of the diffusing species similar to how TMAP4 and TMAP7 treated the diffusing species across two different materials with different chemical potentials. SolubilityRatioMaterial is used to treat solute concentration differences stemming from a difference in T solubilities across the interface. The solubility ratio jump is calculated via the following: $(5)var element = document.getElementById("moose-equation-2c8ddebe-7421-4d57-ad3a-e1c764acf523");katex.render(" J = \\frac{C_{s,1}}{K_{s,1}} - \\frac{C_{s,2}}{K_{s,2}}.", element, {displayMode:true,throwOnError:false});$

ADPenaltyInterfaceDiffusion is used to conserve the particle flux at this interface between two different solubilities. The extremely low tritium solubility in W at low temperature leads to an extremely low solute concentration in W, creating a large solute concentration difference between W and Cu. Two-component solubility in W is used to keep the maximum solubility ratio between W and Cu to 104 at low temperature to avoid the convergence issue associated with calculating two significantly different solute concentration in W and Cu.

[InterfaceKernels<<<{"href": "../../syntax/InterfaceKernels/index.html"}>>>]
  [tied_4to3]
    type = ADPenaltyInterfaceDiffusion<<<{"description": "A penalty-based interface condition that forcesthe continuity of variables and the flux equivalence across an interface.", "href": "../../source/interfacekernels/PenaltyInterfaceDiffusion.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
    neighbor_var<<<{"description": "The variable on the other side of the interface."}>>> = C_mobile_Cu
    penalty<<<{"description": "The penalty that penalizes jump between primary and neighbor variables."}>>> = 0.05 #  it will not converge with > 0.1, but it creates negative C_mobile _Cu with << 0.1
    # jump_prop_name = solubility_ratio_4to3
    jump_prop_name<<<{"description": "the name of the material property that calculates the jump."}>>> = solubility_ratio
    boundary<<<{"description": "The list of boundaries (ids or names) from the mesh where this object applies"}>>> = '4to3'
  []
  [tied_3to2]
    type = ADPenaltyInterfaceDiffusion<<<{"description": "A penalty-based interface condition that forcesthe continuity of variables and the flux equivalence across an interface.", "href": "../../source/interfacekernels/PenaltyInterfaceDiffusion.html"}>>>
    variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_Cu
    neighbor_var<<<{"description": "The variable on the other side of the interface."}>>> = C_mobile_CuCrZr
    penalty<<<{"description": "The penalty that penalizes jump between primary and neighbor variables."}>>> = 0.05 #  it will not converge with > 0.1, but it creates negative C_mobile _Cu with << 0.1
    # jump_prop_name = solubility_ratio_3to2
    jump_prop_name<<<{"description": "the name of the material property that calculates the jump."}>>> = solubility_ratio
    boundary<<<{"description": "The list of boundaries (ids or names) from the mesh where this object applies"}>>> = '3to2'
  []
[]

Numerical method

We use a standard preconditioner: the “single matrix preconditioner”. The Newton method is used to model the transient tritium and thermal transport in a 2D monoblock. It is important to note that MOOSE is equipped with the built-in Message Passing Interface (MPI) protocol, as tritium and thermal transport analysis of fifty 1,600-second cycle plasma discharges in the 2D monoblock is performed in under 2 hours using a single device/computer (3.5 GHz Apple M2 Pro, 10-Core CPU/16-Core GPU) with this MPI feature.

[Executioner<<<{"href": "../../syntax/Executioner/index.html"}>>>]
  type = Transient
  scheme = bdf2
  solve_type = NEWTON
  petsc_options_iname = '-pc_type'
  petsc_options_value = 'lu'
  nl_rel_tol = 1e-6 # 1e-8 works for 1 cycle
  nl_abs_tol = 1e-7 # 1e-11 works for 1 cycle
  end_time = 8.0e4 # 50 ITER shots (3.0e4 s plasma, 2.0e4 SSP)
  automatic_scaling = true
  line_search = 'none'
  dtmin = 1e-4
  nl_max_its = 18
  [TimeStepper<<<{"href": "../../syntax/Executioner/TimeStepper/index.html"}>>>]
    type = IterationAdaptiveDT
    dt = 20
    optimal_iterations = 15
    iteration_window = 1
    growth_factor = 1.2
    cutback_factor = 0.8
    timestep_limiting_postprocessor = timestep_max_pp
  []
[]

Results

The simulation results from this example are shown in Figure 3 and Figure 4. For more results, information, and discussion about the results for this example case and their significance, the reader is referred to Shimada et al. (2024).

Figure 3: Tritium concentration profile in W (left), Cu (center), and CuCrZr (right) after ten 1,600-second cycles ( $var element = document.getElementById("moose-equation-0b1242ad-334c-47ce-9eb8-5ce8b64e771d");katex.render("t = 14912", element, {displayMode:false,throwOnError:false});$ s). This corresponds to Fig. 4A in Shimada et al. (2024).

Figure 4: Tritium concentration profile in W (left), Cu (center), and CuCrZr (right) after fifty 1,600-second cycles ( $var element = document.getElementById("moose-equation-d646b961-798f-41bf-82df-1d0a8156ee88");katex.render("t = 78912", element, {displayMode:false,throwOnError:false});$ s). This corresponds to Fig. 4B in Shimada et al. (2024).

warning:The exodus file in gold is a smaller version of the output

The input file (test/tests/divertor_monoblock/divertor_monoblock.i) returns the outputs that were used in Shimada et al. (2024). However, a slightly modified version of this input is run in (test/tests/divertor_monoblock/tests) as part of TMAP8's Software Quality Assurance process: It simulates only one pulse cycle, has a coarser mesh, and outputs the results less regularly to limit the file size. As a result, the exodus file in the test gold directory is a smaller version of the output generated when running the full input file.

Complete input file

Below is the complete input file, which can be run reliably with approximately 4 processor cores. Note that this input file has not been optimized for computational costs.

### Nomenclatures
### C_mobile_j      mobile H concentration in "j" material, where j = CuCrZr, Cu, W
### C_trapped_j     trapped H concentration in "j" material, where j = CuCrZr, Cu, W
### C_total_j       total H concentration in "j" material, where j = CuCrZr, Cu, W
###
### S_empty_j       empty site concentration in "j" material, where j = CuCrZr, Cu, W
### S_trapped_j     trapped site concentration in "j" material, where j = CuCrZr, Cu, W
### S_total_j       total site H concentration in "j" material, where j = CuCrZr, Cu, W
###
### F_permeation    permeation flux
### F_recombination recombination flux
###
### Sc_             Scaled
### Int_            Integrated
### ScInt_          Scaled and integrated

tungsten_atomic_density = ${units 6.338e28 m^-3}

## Mesh generated by ConcentricCircleMeshGenerator and saved as "2DMonoblock.e"
[Mesh<<<{"href": "../../syntax/Mesh/index.html"}>>>]
    [ccmg]
        type = ConcentricCircleMeshGenerator<<<{"description": "This ConcentricCircleMeshGenerator source code is to generate concentric circle meshes.", "href": "../../source/meshgenerators/ConcentricCircleMeshGenerator.html"}>>>
        num_sectors<<<{"description": "num_sectors % 2 = 0, num_sectors > 0Number of azimuthal sectors in each quadrant'num_sectors' must be an even number."}>>> = 36
        rings<<<{"description": "Number of rings in each circle or in the enclosing square"}>>> = '1 30 20 110'
        radii<<<{"description": "Radii of major concentric circles"}>>> = '${units 6 mm -> m} ${units 7.5 mm -> m} ${units 8.5 mm -> m}'
        has_outer_square<<<{"description": "It determines if meshes for a outer square are added to concentric circle meshes."}>>> = on
        pitch<<<{"description": "The enclosing square can be added to the completed concentric circle mesh.Elements are quad meshes."}>>> = ${units 28 mm -> m}
        portion<<<{"description": "Control of which part of mesh is created"}>>> = left_half
        preserve_volumes<<<{"description": "Volume of concentric circles can be preserved using this function."}>>> = false
        smoothing_max_it<<<{"description": "Number of Laplacian smoothing iterations"}>>> = 3
    []
    [ssbsg1]
        type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
        input<<<{"description": "The mesh we want to modify"}>>> = ccmg
        primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '4'     # W
        paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '3'      # Cu
        new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '4to3'
    []
    [ssbsg2]
        type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
        input<<<{"description": "The mesh we want to modify"}>>> = ssbsg1
        primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '3'     # Cu
        paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '4'      # W
        new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '3to4'
    []
    [ssbsg3]
        type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
        input<<<{"description": "The mesh we want to modify"}>>> = ssbsg2
        primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '3'     # Cu
        paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '2'      # CuCrZr
        new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '3to2'
    []
    [ssbsg4]
        type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
        input<<<{"description": "The mesh we want to modify"}>>> = ssbsg3
        primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '2'     # CuCrZr
        paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '3'      # Cu
        new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '2to3'
    []
    [ssbsg5]
        type = SideSetsBetweenSubdomainsGenerator<<<{"description": "MeshGenerator that creates a sideset composed of the nodes located between two or more subdomains.", "href": "../../source/meshgenerators/SideSetsBetweenSubdomainsGenerator.html"}>>>
        input<<<{"description": "The mesh we want to modify"}>>> = ssbsg4
        primary_block<<<{"description": "The primary set of blocks for which to draw a sideset between"}>>> = '2'     # CuCrZr
        paired_block<<<{"description": "The paired set of blocks for which to draw a sideset between"}>>> = '1'      # H2O
        new_boundary<<<{"description": "The list of boundary names to create on the supplied subdomain"}>>> = '2to1'
    []
    [bdg]
        type = BlockDeletionGenerator<<<{"description": "Mesh generator which removes elements from the specified subdomains", "href": "../../source/meshgenerators/BlockDeletionGenerator.html"}>>>
        input<<<{"description": "The mesh we want to modify"}>>> = ssbsg5
        block<<<{"description": "The list of blocks to be deleted"}>>> = '1'             # H2O
    []
[]

[Problem<<<{"href": "../../syntax/Problem/index.html"}>>>]
    type = ReferenceResidualProblem
    extra_tag_vectors = 'ref'
    reference_vector = 'ref'
[]

[Variables<<<{"href": "../../syntax/Variables/index.html"}>>>]
    [temperature]
        order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
        family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
        initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = ${units 300 K}
    []
    ######################### Variables for W (block = 4)
    [C_mobile_W]
        order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
        family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
        initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = ${units 1.0e-20 m^-3}
        block = 4
    []
    [C_trapped_W]
        order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
        family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
        initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = ${units 1.0e-15 m^-3}
        block = 4
    []
    ######################### Variables for Cu (block = 3)
    [C_mobile_Cu]
        order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
        family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
        initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = ${units 5.0e-17 m^-3}
        block = 3
    []
    [C_trapped_Cu]
        order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
        family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
        initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = ${units 1.0e-15 m^-3}
        block = 3
    []
    ######################### Variables for CuCrZr (block = 2)
    [C_mobile_CuCrZr]
        order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
        family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
        initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = ${units 1.0e-15 m^-3}
        block = 2
    []
    [C_trapped_CuCrZr]
        order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
        family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = LAGRANGE
        initial_condition<<<{"description": "Specifies a constant initial condition for this variable"}>>> = ${units 1.0e-15 m^-3}
        block = 2
    []
[]

[AuxVariables<<<{"href": "../../syntax/AuxVariables/index.html"}>>>]
    [flux_y]
        order<<<{"description": "Specifies the order of the FE shape function to use for this variable (additional orders not listed are allowed)"}>>> = FIRST
        family<<<{"description": "Specifies the family of FE shape functions to use for this variable"}>>> = MONOMIAL
    []
    ############################## AuxVariables for W (block = 4)
    [Sc_C_mobile_W]
        block = 4
    []
    [Sc_C_trapped_W]
        block = 4
    []
    [C_total_W]
        block = 4
    []
    [Sc_C_total_W]
        block = 4
    []
    [S_empty_W]
        block = 4
    []
    [Sc_S_empty_W]
        block = 4
    []
    [S_trapped_W]
        block = 4
    []
    [Sc_S_trapped_W]
        block = 4
    []
    [S_total_W]
        block = 4
    []
    [Sc_S_total_W]
        block = 4
    []
    ############################## AuxVariables for Cu (block = 3)
    [Sc_C_mobile_Cu]
        block = 3
    []
    [Sc_C_trapped_Cu]
        block = 3
    []
    [C_total_Cu]
        block = 3
    []
    [Sc_C_total_Cu]
        block = 3
    []
    [S_empty_Cu]
        block = 3
    []
    [Sc_S_empty_Cu]
        block = 3
    []
    [S_trapped_Cu]
        block = 3
    []
    [Sc_S_trapped_Cu]
        block = 3
    []
    [S_total_Cu]
        block = 3
    []
    [Sc_S_total_Cu]
        block = 3
    []
    ############################## AuxVariables for CuCrZr (block = 2)
    [Sc_C_mobile_CuCrZr]
        block = 2
    []
    [Sc_C_trapped_CuCrZr]
        block = 2
    []
    [C_total_CuCrZr]
        block = 2
    []
    [Sc_C_total_CuCrZr]
        block = 2
    []
    [S_empty_CuCrZr]
        block = 2
    []
    [Sc_S_empty_CuCrZr]
        block = 2
    []
    [S_trapped_CuCrZr]
        block = 2
    []
    [Sc_S_trapped_CuCrZr]
        block = 2
    []
    [S_total_CuCrZr]
        block = 2
    []
    [Sc_S_total_CuCrZr]
        block = 2
    []
[]

[Kernels<<<{"href": "../../syntax/Kernels/index.html"}>>>]
    ############################## Kernels for W (block = 4)
    [diff_W]
        type = ADMatDiffusion<<<{"description": "Diffusion equation kernel that takes an isotropic diffusivity from a material property", "href": "../../source/kernels/ADMatDiffusion.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
        diffusivity<<<{"description": "The diffusivity value or material property"}>>> = diffusivity_W
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [time_diff_W]
        type = ADTimeDerivative<<<{"description": "The time derivative operator with the weak form of $(\\psi_i, \\frac{\\partial u_h}{\\partial t})$.", "href": "../../source/kernels/ADTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [coupled_time_W]
        type = ScaledCoupledTimeDerivative<<<{"description": "Time derivative Kernel that acts on a coupled variable. Weak form: $(\\psi_i, \\frac{\\partial v_h}{\\partial t})$.", "href": "../../source/kernels/ScaledCoupledTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
        v<<<{"description": "Coupled variable"}>>> = C_trapped_W
        factor<<<{"description": "The factor by which to scale"}>>> = 1e0
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [heat_conduction_W]
        type = HeatConduction<<<{"description": "Diffusive heat conduction term $-\\nabla\\cdot(k\\nabla T)$ of the thermal energy conservation equation", "href": "../../source/kernels/HeatConduction.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
        diffusion_coefficient<<<{"description": "Property name of the diffusion coefficient"}>>> = thermal_conductivity_W
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [time_heat_conduction_W]
        type = SpecificHeatConductionTimeDerivative<<<{"description": "Time derivative term $\\rho c_p \\frac{\\partial T}{\\partial t}$ of the heat equation with the specific heat $c_p$ and the density $\\rho$ as arguments.", "href": "../../source/kernels/SpecificHeatConductionTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
        specific_heat<<<{"description": "Property name of the specific heat material property"}>>> = specific_heat_W
        density<<<{"description": "Property name of the density material property"}>>> = density_W
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    ############################## Kernels for Cu (block = 3)
    [diff_Cu]
        type = ADMatDiffusion<<<{"description": "Diffusion equation kernel that takes an isotropic diffusivity from a material property", "href": "../../source/kernels/ADMatDiffusion.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_Cu
        diffusivity<<<{"description": "The diffusivity value or material property"}>>> = diffusivity_Cu
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [time_diff_Cu]
        type = ADTimeDerivative<<<{"description": "The time derivative operator with the weak form of $(\\psi_i, \\frac{\\partial u_h}{\\partial t})$.", "href": "../../source/kernels/ADTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_Cu
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [coupled_time_Cu]
        type = ScaledCoupledTimeDerivative<<<{"description": "Time derivative Kernel that acts on a coupled variable. Weak form: $(\\psi_i, \\frac{\\partial v_h}{\\partial t})$.", "href": "../../source/kernels/ScaledCoupledTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_Cu
        v<<<{"description": "Coupled variable"}>>> = C_trapped_Cu
        factor<<<{"description": "The factor by which to scale"}>>> = 1e0
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [heat_conduction_Cu]
        type = HeatConduction<<<{"description": "Diffusive heat conduction term $-\\nabla\\cdot(k\\nabla T)$ of the thermal energy conservation equation", "href": "../../source/kernels/HeatConduction.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
        diffusion_coefficient<<<{"description": "Property name of the diffusion coefficient"}>>> = thermal_conductivity_Cu
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [time_heat_conduction_Cu]
        type = SpecificHeatConductionTimeDerivative<<<{"description": "Time derivative term $\\rho c_p \\frac{\\partial T}{\\partial t}$ of the heat equation with the specific heat $c_p$ and the density $\\rho$ as arguments.", "href": "../../source/kernels/SpecificHeatConductionTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
        specific_heat<<<{"description": "Property name of the specific heat material property"}>>> = specific_heat_Cu
        density<<<{"description": "Property name of the density material property"}>>> = density_Cu
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    ############################## Kernels for CuCrZr (block = 2)
    [diff_CuCrZr]
        type = ADMatDiffusion<<<{"description": "Diffusion equation kernel that takes an isotropic diffusivity from a material property", "href": "../../source/kernels/ADMatDiffusion.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_CuCrZr
        diffusivity<<<{"description": "The diffusivity value or material property"}>>> = diffusivity_CuCrZr
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [time_diff_CuCrZr]
        type = ADTimeDerivative<<<{"description": "The time derivative operator with the weak form of $(\\psi_i, \\frac{\\partial u_h}{\\partial t})$.", "href": "../../source/kernels/ADTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_CuCrZr
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [coupled_time_CuCrZr]
        type = ScaledCoupledTimeDerivative<<<{"description": "Time derivative Kernel that acts on a coupled variable. Weak form: $(\\psi_i, \\frac{\\partial v_h}{\\partial t})$.", "href": "../../source/kernels/ScaledCoupledTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_CuCrZr
        v<<<{"description": "Coupled variable"}>>> = C_trapped_CuCrZr
        factor<<<{"description": "The factor by which to scale"}>>> = 1e0
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [heat_conduction_CuCrZr]
        type = HeatConduction<<<{"description": "Diffusive heat conduction term $-\\nabla\\cdot(k\\nabla T)$ of the thermal energy conservation equation", "href": "../../source/kernels/HeatConduction.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
        diffusion_coefficient<<<{"description": "Property name of the diffusion coefficient"}>>> = thermal_conductivity_CuCrZr
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [time_heat_conduction_CuCrZr]
        type = SpecificHeatConductionTimeDerivative<<<{"description": "Time derivative term $\\rho c_p \\frac{\\partial T}{\\partial t}$ of the heat equation with the specific heat $c_p$ and the density $\\rho$ as arguments.", "href": "../../source/kernels/SpecificHeatConductionTimeDerivative.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
        specific_heat<<<{"description": "Property name of the specific heat material property"}>>> = specific_heat_CuCrZr
        density<<<{"description": "Property name of the density material property"}>>> = density_CuCrZr
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
[]

[AuxKernels<<<{"href": "../../syntax/AuxKernels/index.html"}>>>]
    ############################## AuxKernels for W (block = 4)
    [Scaled_mobile_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_mobile_W
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_mobile_W
    []
    [Scaled_trapped_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_trapped_W
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_W
    []
    [total_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = C_total_W
        type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
        expression<<<{"description": "Parsed function expression to compute"}>>> = 'C_mobile_W + C_trapped_W'
        coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'C_mobile_W C_trapped_W'
    []
    [Scaled_total_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_total_W
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_total_W
    []
    [empty_sites_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_empty_W
        type = EmptySitesAux<<<{"description": "Calculates the concentration of empty trapping sites.", "href": "../../source/auxkernels/EmptySitesAux.html"}>>>
        N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = ${units 1.0e0 m^-3}       # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        # Ct0 = ${units 1.0e-4 m^-3}   # E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = ${units 1.0e-4 m^-3}    # E.A. Hodille et al 2021 Nucl. Fusion 61 1260033, trap 2
        trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix"}>>> = 1.0e0         # 1.0e1
        trapped_concentration_variables<<<{"description": "Variables representing trapped particle concentrations."}>>> = C_trapped_W
    []
    [scaled_empty_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_empty_W
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_empty_W
    []
    [trapped_sites_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_trapped_W
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = 1e0
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_W
    []
    [scaled_trapped_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_trapped_W
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_trapped_W
    []
    [total_sites_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_total_W
        type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
        expression<<<{"description": "Parsed function expression to compute"}>>> = 'S_trapped_W + S_empty_W'
        coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'S_trapped_W S_empty_W'
    []
    [scaled_total_W]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_total_W
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_total_W
    []
    ############################## AuxKernels for Cu (block = 3)
    [Scaled_mobile_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_mobile_Cu
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_mobile_Cu
    []
    [Scaled_trapped_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_trapped_Cu
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_Cu
    []
    [total_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = C_total_Cu
        type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
        expression<<<{"description": "Parsed function expression to compute"}>>> = 'C_mobile_Cu + C_trapped_Cu'
        coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'C_mobile_Cu C_trapped_Cu'
    []
    [Scaled_total_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_total_Cu
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_total_Cu
    []
    [empty_sites_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_empty_Cu
        type = EmptySitesAux<<<{"description": "Calculates the concentration of empty trapping sites.", "href": "../../source/auxkernels/EmptySitesAux.html"}>>>
        N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = ${units 1.0e0 m^-3}     # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = ${units 5.0e-5 m^-3}  # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix"}>>> = 1.0e0       # 1.0e1
        trapped_concentration_variables<<<{"description": "Variables representing trapped particle concentrations."}>>> = C_trapped_Cu
    []
    [scaled_empty_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_empty_Cu
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_empty_Cu
    []
    [trapped_sites_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_trapped_Cu
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = 1e0
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_Cu
    []
    [scaled_trapped_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_trapped_Cu
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_trapped_Cu
    []
    [total_sites_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_total_Cu
        type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
        expression<<<{"description": "Parsed function expression to compute"}>>> = 'S_trapped_Cu + S_empty_Cu'
        coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'S_trapped_Cu S_empty_Cu'
    []
    [scaled_total_Cu]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_total_Cu
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_total_Cu
    []
    ############################## AuxKernels for CuCrZr (block = 2)
    [Scaled_mobile_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_mobile_CuCrZr
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_mobile_CuCrZr
    []
    [Scaled_trapped_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_trapped_CuCrZr
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_CuCrZr
    []
    [total_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = C_total_CuCrZr
        type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
        expression<<<{"description": "Parsed function expression to compute"}>>> = 'C_mobile_CuCrZr + C_trapped_CuCrZr'
        coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'C_mobile_CuCrZr C_trapped_CuCrZr'
    []
    [Scaled_total_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_C_total_CuCrZr
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_total_CuCrZr
    []
    [empty_sites_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_empty_CuCrZr
        type = EmptySitesAux<<<{"description": "Calculates the concentration of empty trapping sites.", "href": "../../source/auxkernels/EmptySitesAux.html"}>>>
        N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = ${units 1.0e0 m^-3}     # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = ${units 5.0e-5 m^-3}  # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # Ct0 = ${units 4.0e-2 m^-3} # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix"}>>> = 1.0e0       # 1.0e1
        trapped_concentration_variables<<<{"description": "Variables representing trapped particle concentrations."}>>> = C_trapped_CuCrZr
    []
    [scaled_empty_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_empty_CuCrZr
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_empty_CuCrZr
    []
    [trapped_sites_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_trapped_CuCrZr
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = 1e0
        source_variable<<<{"description": "The variable to be normalized"}>>> = C_trapped_CuCrZr
    []
    [scaled_trapped_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_trapped_CuCrZr
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_trapped_CuCrZr
    []
    [total_sites_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = S_total_CuCrZr
        type = ParsedAux<<<{"description": "Sets a field variable value to the evaluation of a parsed expression.", "href": "../../source/auxkernels/ParsedAux.html"}>>>
        expression<<<{"description": "Parsed function expression to compute"}>>> = 'S_trapped_CuCrZr + S_empty_CuCrZr'
        coupled_variables<<<{"description": "Vector of coupled variable names"}>>> = 'S_trapped_CuCrZr S_empty_CuCrZr'
    []
    [scaled_total_CuCrZr]
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = Sc_S_total_CuCrZr
        type = NormalizationAux<<<{"description": "Normalizes a variable based on a Postprocessor value.", "href": "../../source/auxkernels/NormalizationAux.html"}>>>
        normal_factor<<<{"description": "The normalization factor"}>>> = ${tungsten_atomic_density}
        source_variable<<<{"description": "The variable to be normalized"}>>> = S_total_CuCrZr
    []
    [flux_y_W]
        type = DiffusionFluxAux<<<{"description": "Compute components of flux vector for diffusion problems $(\\vec{J} = -D \\nabla C)$.", "href": "../../source/auxkernels/DiffusionFluxAux.html"}>>>
        diffusivity<<<{"description": "The name of the diffusivity material property that will be used in the flux computation."}>>> = diffusivity_W
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = flux_y
        diffusion_variable<<<{"description": "The name of the variable"}>>> = C_mobile_W
        component<<<{"description": "The desired component of flux."}>>> = y
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    []
    [flux_y_Cu]
        type = DiffusionFluxAux<<<{"description": "Compute components of flux vector for diffusion problems $(\\vec{J} = -D \\nabla C)$.", "href": "../../source/auxkernels/DiffusionFluxAux.html"}>>>
        diffusivity<<<{"description": "The name of the diffusivity material property that will be used in the flux computation."}>>> = diffusivity_Cu
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = flux_y
        diffusion_variable<<<{"description": "The name of the variable"}>>> = C_mobile_Cu
        component<<<{"description": "The desired component of flux."}>>> = y
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    []
    [flux_y_CuCrZr]
        type = DiffusionFluxAux<<<{"description": "Compute components of flux vector for diffusion problems $(\\vec{J} = -D \\nabla C)$.", "href": "../../source/auxkernels/DiffusionFluxAux.html"}>>>
        diffusivity<<<{"description": "The name of the diffusivity material property that will be used in the flux computation."}>>> = diffusivity_CuCrZr
        variable<<<{"description": "The name of the variable that this object applies to"}>>> = flux_y
        diffusion_variable<<<{"description": "The name of the variable"}>>> = C_mobile_CuCrZr
        component<<<{"description": "The desired component of flux."}>>> = y
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    []
[]

[InterfaceKernels<<<{"href": "../../syntax/InterfaceKernels/index.html"}>>>]
    [tied_4to3]
        type = ADPenaltyInterfaceDiffusion<<<{"description": "A penalty-based interface condition that forcesthe continuity of variables and the flux equivalence across an interface.", "href": "../../source/interfacekernels/PenaltyInterfaceDiffusion.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
        neighbor_var<<<{"description": "The variable on the other side of the interface."}>>> = C_mobile_Cu
        penalty<<<{"description": "The penalty that penalizes jump between primary and neighbor variables."}>>> = 0.05            #  it will not converge with > 0.1, but it creates negative C_mobile _Cu with << 0.1
        # jump_prop_name = solubility_ratio_4to3
        jump_prop_name<<<{"description": "the name of the material property that calculates the jump."}>>> = solubility_ratio
        boundary<<<{"description": "The list of boundaries (ids or names) from the mesh where this object applies"}>>> = '4to3'
    []
    [tied_3to2]
        type = ADPenaltyInterfaceDiffusion<<<{"description": "A penalty-based interface condition that forcesthe continuity of variables and the flux equivalence across an interface.", "href": "../../source/interfacekernels/PenaltyInterfaceDiffusion.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_Cu
        neighbor_var<<<{"description": "The variable on the other side of the interface."}>>> = C_mobile_CuCrZr
        penalty<<<{"description": "The penalty that penalizes jump between primary and neighbor variables."}>>> = 0.05            #  it will not converge with > 0.1, but it creates negative C_mobile _Cu with << 0.1
        # jump_prop_name = solubility_ratio_3to2
        jump_prop_name<<<{"description": "the name of the material property that calculates the jump."}>>> = solubility_ratio
        boundary<<<{"description": "The list of boundaries (ids or names) from the mesh where this object applies"}>>> = '3to2'
    []
[]

[NodalKernels<<<{"href": "../../syntax/NodalKernels/index.html"}>>>]
    ############################## NodalKernels for W (block = 4)
    [time_W]
        type = TimeDerivativeNodalKernel<<<{"description": "Forms the contribution to the residual and jacobian of the time derivative term from an ODE being solved at all nodes.", "href": "../../source/nodalkernels/TimeDerivativeNodalKernel.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_W
    []
    [trapping_W]
        type = TrappingNodalKernel<<<{"description": "Implements a residual describing the trapping of a species in a material.", "href": "../../source/nodal_kernels/TrappingNodalKernel.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_W
        temperature<<<{"description": "The temperature (K)"}>>> = temperature
        alpha_t<<<{"description": "The trapping rate coefficient. This has units of 1/s (e.g. no number densities are involved)"}>>> = 2.75e11      # 1e15
        N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = 1.0e0  # = (1e0) x (${tungsten_atomic_density} #/m^3)
        # Ct0 = 1.0e-4                # E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = 1.0e-4                # E.A. Hodille et al 2021 Nucl. Fusion 61 1260033, trap 2
        trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix (-)"}>>> = 1.0e0       # 1.0e1
        mobile_concentration<<<{"description": "The variable representing the mobile concentration of solute particles (1/m^3)"}>>> = 'C_mobile_W'
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [release_W]
        type = ReleasingNodalKernel<<<{"description": "Implements a residual describing the release of trapped species in a material.", "href": "../../source/nodal_kernels/ReleasingNodalKernel.html"}>>>
        alpha_r<<<{"description": "The release rate coefficient (1/s)"}>>> = 8.4e12    # 1.0e13
        temperature<<<{"description": "The temperature (K)"}>>> = temperature
        # detrapping_energy = 9863.9    # = 0.85 eV    E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        detrapping_energy<<<{"description": "The detrapping energy (K)"}>>> = 11604.6   # = 1.00 eV    E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 2
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_W
    []
    ############################## NodalKernels for Cu (block = 3)
    [time_Cu]
        type = TimeDerivativeNodalKernel<<<{"description": "Forms the contribution to the residual and jacobian of the time derivative term from an ODE being solved at all nodes.", "href": "../../source/nodalkernels/TimeDerivativeNodalKernel.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_Cu
    []
    [trapping_Cu]
        type = TrappingNodalKernel<<<{"description": "Implements a residual describing the trapping of a species in a material.", "href": "../../source/nodal_kernels/TrappingNodalKernel.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_Cu
        temperature<<<{"description": "The temperature (K)"}>>> = temperature
        alpha_t<<<{"description": "The trapping rate coefficient. This has units of 1/s (e.g. no number densities are involved)"}>>> = 2.75e11      # 1e15
        N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = 1.0e0  # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = 5.0e-5                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix (-)"}>>> = 1.0e0       # 1.0e1
        mobile_concentration<<<{"description": "The variable representing the mobile concentration of solute particles (1/m^3)"}>>> = 'C_mobile_Cu'
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [release_Cu]
        type = ReleasingNodalKernel<<<{"description": "Implements a residual describing the release of trapped species in a material.", "href": "../../source/nodal_kernels/ReleasingNodalKernel.html"}>>>
        alpha_r<<<{"description": "The release rate coefficient (1/s)"}>>> = 8.4e12    # 1.0e13
        temperature<<<{"description": "The temperature (K)"}>>> = temperature
        detrapping_energy<<<{"description": "The detrapping energy (K)"}>>> = 5802.3    # = 0.50eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_Cu
    []
    ############################## NodalKernels for CuCrZr (block = 2)
    [time_CuCrZr]
        type = TimeDerivativeNodalKernel<<<{"description": "Forms the contribution to the residual and jacobian of the time derivative term from an ODE being solved at all nodes.", "href": "../../source/nodalkernels/TimeDerivativeNodalKernel.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_CuCrZr
    []
    [trapping_CuCrZr]
        type = TrappingNodalKernel<<<{"description": "Implements a residual describing the trapping of a species in a material.", "href": "../../source/nodal_kernels/TrappingNodalKernel.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_CuCrZr
        temperature<<<{"description": "The temperature (K)"}>>> = temperature
        alpha_t<<<{"description": "The trapping rate coefficient. This has units of 1/s (e.g. no number densities are involved)"}>>> = 2.75e11      # 1e15
        N<<<{"description": "The atomic number density of the host material (1/m^3)"}>>> = 1.0e0  # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0<<<{"description": "The fraction of host sites that can contribute to trapping as a function (-)"}>>> = 5.0e-5                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # Ct0 = 4.0e-2                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        trap_per_free<<<{"description": "An estimate for the ratio of the concentration magnitude of trapped species to free species. Setting a value for this can be helpful in producing a well-scaled matrix (-)"}>>> = 1.0e0       # 1.0e1
        mobile_concentration<<<{"description": "The variable representing the mobile concentration of solute particles (1/m^3)"}>>> = 'C_mobile_CuCrZr'
        extra_vector_tags<<<{"description": "The extra tags for the vectors this Kernel should fill"}>>> = ref
    []
    [release_CuCrZr]
        type = ReleasingNodalKernel<<<{"description": "Implements a residual describing the release of trapped species in a material.", "href": "../../source/nodal_kernels/ReleasingNodalKernel.html"}>>>
        alpha_r<<<{"description": "The release rate coefficient (1/s)"}>>> = 8.4e12    # 1.0e13
        temperature<<<{"description": "The temperature (K)"}>>> = temperature
        detrapping_energy<<<{"description": "The detrapping energy (K)"}>>> = 5802.3    # = 0.50eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # detrapping_energy = 9631.8   # = 0.83 eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_trapped_CuCrZr
    []
[]

[BCs<<<{"href": "../../syntax/BCs/index.html"}>>>]
    [C_mob_W_top_flux]
        type = FunctionNeumannBC<<<{"description": "Imposes the integrated boundary condition $\\frac{\\partial u}{\\partial n}=h(t,\\vec{x})$, where $h$ is a (possibly) time and space-dependent MOOSE Function.", "href": "../../source/bcs/FunctionNeumannBC.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_W
        boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = 'top'
        function<<<{"description": "The function."}>>> = mobile_flux_bc_func
    []
    [mobile_tube]
        type = DirichletBC<<<{"description": "Imposes the essential boundary condition $u=g$, where $g$ is a constant, controllable value.", "href": "../../source/bcs/DirichletBC.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = C_mobile_CuCrZr
        value<<<{"description": "Value of the BC"}>>> = 1.0e-18
        boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = '2to1'
    []
    [temp_top]
        type = FunctionNeumannBC<<<{"description": "Imposes the integrated boundary condition $\\frac{\\partial u}{\\partial n}=h(t,\\vec{x})$, where $h$ is a (possibly) time and space-dependent MOOSE Function.", "href": "../../source/bcs/FunctionNeumannBC.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
        boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = 'top'
        function<<<{"description": "The function."}>>> = temp_flux_bc_func
    []
    [temp_tube]
        type = FunctionDirichletBC<<<{"description": "Imposes the essential boundary condition $u=g(t,\\vec{x})$, where $g$ is a (possibly) time and space-dependent MOOSE Function.", "href": "../../source/bcs/FunctionDirichletBC.html"}>>>
        variable<<<{"description": "The name of the variable that this residual object operates on"}>>> = temperature
        boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = '2to1'
        function<<<{"description": "The forcing function."}>>> = temp_inner_func
    []
[]

[Functions<<<{"href": "../../syntax/Functions/index.html"}>>>]
    ### Maximum mobile flux of 7.90e-13 at the top surface (1.0e-4 [m])
    ### 1.80e23/m^2-s = (5.0e23/m^2-s) *(1-0.999) = (7.90e-13)*(${tungsten_atomic_density})/(1.0e-4)  at steady state
    [mobile_flux_bc_func]
        type = ParsedFunction<<<{"description": "Function created by parsing a string", "href": "../../source/functions/MooseParsedFunction.html"}>>>
        expression<<<{"description": "The user defined function."}>>> =   'if((t % 1600) < 100.0, 0.0      + 7.90e-13*(t % 1600)/100,
                        if((t % 1600) < 500.0, 7.90e-13,
                        if((t % 1600) < 600.0, 7.90e-13 - 7.90e-13*((t % 1600)-500)/100, 0.0)))'
    []
    ### Heat flux of 10MW/m^2 at steady state
    [temp_flux_bc_func]
        type = ParsedFunction<<<{"description": "Function created by parsing a string", "href": "../../source/functions/MooseParsedFunction.html"}>>>
        expression<<<{"description": "The user defined function."}>>> =   'if((t % 1600) < 100.0, 0.0   + 1.0e7*(t % 1600)/100,
                        if((t % 1600) < 500.0, 1.0e7,
                        if((t % 1600) < 600.0, 1.0e7 - 1.0e7*((t % 1600)-500)/100, 300)))'
    []
    ### Maximum coolant temperature of 552K at steady state
    [temp_inner_func]
        type = ParsedFunction<<<{"description": "Function created by parsing a string", "href": "../../source/functions/MooseParsedFunction.html"}>>>
        expression<<<{"description": "The user defined function."}>>> =   'if((t % 1600) < 100.0, 300.0 + (552-300)*(t % 1600)/100,
                        if((t % 1600) < 500.0, 552,
                        if((t % 1600) < 600.0, 552.0 - (552-300)*((t % 1600)-500)/100, 300)))'
    []
    [timestep_function]
        type = ParsedFunction<<<{"description": "Function created by parsing a string", "href": "../../source/functions/MooseParsedFunction.html"}>>>
        expression<<<{"description": "The user defined function."}>>> = 'if((t % 1600) <   10.0,  20,
                      if((t % 1600) <   90.0,  40,
                      if((t % 1600) <  110.0,  20,
                      if((t % 1600) <  480.0,  40,
                      if((t % 1600) <  500.0,  20,
                      if((t % 1600) <  590.0,   4,
                      if((t % 1600) <  610.0,  20,
                      if((t % 1600) < 1500.0, 200,
                      if((t % 1600) < 1600.0,  40,  2)))))))))'
    []
[]

[Materials<<<{"href": "../../syntax/Materials/index.html"}>>>]
    ############################## Materials for W (block = 4)
    [diffusivity_W]
        type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = diffusivity_W
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '2.4e-7*exp(-4525.8/temperature)'    # H diffusivity in W
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    [solubility_W]
        type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = solubility_W
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        # expression = '2.95e-5 *exp(-12069.0/temperature)'              # H solubility in W = (1.87e24)/(${tungsten_atomic_density}) [#/m^3]
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '2.95e-5 *exp(-12069.0/temperature) + 4.95e-8 * exp(-6614.6/temperature)'    # H solubility in W = (1.87e24)/(${tungsten_atomic_density}) [#/m^3]
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    [converter_to_regular_W]
        type = MaterialADConverter<<<{"description": "Converts regular material properties to AD properties and vice versa", "href": "../../source/materials/MaterialADConverter.html"}>>>
        ad_props_in<<<{"description": "The names of the AD material properties to convert to regular properties"}>>> = 'diffusivity_W'
        reg_props_out<<<{"description": "The names of the output regular properties"}>>> = 'diffusivity_W_nonAD'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    []
    [heat_transfer_W]
        type = GenericConstantMaterial<<<{"description": "Declares material properties based on names and values prescribed by input parameters.", "href": "../../source/materials/GenericConstantMaterial.html"}>>>
        prop_names<<<{"description": "The names of the properties this material will have"}>>> = 'density_W'
        prop_values<<<{"description": "The values associated with the named properties"}>>> = '19300'                # [g/m^3]
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    []
    [specific_heat_W]
        type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = specific_heat_W
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '1.16e2 + 7.11e-2 * temperature - 6.58e-5 * temperature^2 + 3.24e-8 * temperature^3 -5.45e-12 * temperature^4'    # ~ 132[J/kg-K]
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    [thermal_conductivity_W]
        type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = thermal_conductivity_W
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
        # expression = '-7.8e-9 * temperature^3 + 5.0e-5 * temperature^2 - 1.1e-1 * temperature + 1.8e2'    # ~ 173.0 [ W/m-K]   from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '2.41e2 - 2.90e-1 * temperature + 2.54e-4 * temperature^2 - 1.03e-7 * temperature^3 + 1.52e-11 * temperature^4'    # ~ 173.0 [ W/m-K]
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    ############################## Materials for Cu (block = 3)
    [diffusivity_Cu]
        type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = diffusivity_Cu
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '6.60e-7*exp(-4525.8/temperature)'    # H diffusivity in Cu
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    [solubility_Cu]
        type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = solubility_Cu
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '4.95e-5*exp(-6614.6/temperature)'    # H solubility in Cu = (3.14e24)/(${tungsten_atomic_density}) [#/m^3]
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    [converter_to_regular_Cu]
        type = MaterialADConverter<<<{"description": "Converts regular material properties to AD properties and vice versa", "href": "../../source/materials/MaterialADConverter.html"}>>>
        ad_props_in<<<{"description": "The names of the AD material properties to convert to regular properties"}>>> = 'diffusivity_Cu'
        reg_props_out<<<{"description": "The names of the output regular properties"}>>> = 'diffusivity_Cu_nonAD'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    []
    [heat_transfer_Cu]
        type = GenericConstantMaterial<<<{"description": "Declares material properties based on names and values prescribed by input parameters.", "href": "../../source/materials/GenericConstantMaterial.html"}>>>
        prop_names<<<{"description": "The names of the properties this material will have"}>>> = 'density_Cu'
        prop_values<<<{"description": "The values associated with the named properties"}>>> = '8960.0'                # [g/m^3]
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    []
    [specific_heat_Cu]
        type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = specific_heat_Cu
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '3.16e2 + 3.18e-1 * temperature - 3.49e-4 * temperature^2 + 1.66e-7 * temperature^3'    # ~ 384 [J/kg-K]
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    [thermal_conductivity_Cu]
        type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = thermal_conductivity_Cu
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
        # expression = '-3.9e-8 * temperature^3 + 3.8e-5 * temperature^2 - 7.9e-2 * temperature + 4.0e2'    # ~ 401.0  [ W/m-K] from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '4.21e2 - 6.85e-2 * temperature'    # ~ 400.0 [ W/m-K]
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    ############################## Materials for CuCrZr (block = 2)
    [diffusivity_CuCrZr]
        type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = diffusivity_CuCrZr
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '3.90e-7*exp(-4873.9/temperature)'    # H diffusivity in CuCrZr
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    [solubility_CuCrZr]
        type = ADParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = solubility_CuCrZr
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '6.75e-6*exp(-4525.8/temperature)'    # H solubility in CuCrZr = (4.28e23)/(${tungsten_atomic_density}) [#/m^3]
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    [converter_to_regular_CuCrZr]
        type = MaterialADConverter<<<{"description": "Converts regular material properties to AD properties and vice versa", "href": "../../source/materials/MaterialADConverter.html"}>>>
        ad_props_in<<<{"description": "The names of the AD material properties to convert to regular properties"}>>> = 'diffusivity_CuCrZr'
        reg_props_out<<<{"description": "The names of the output regular properties"}>>> = 'diffusivity_CuCrZr_nonAD'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    []
    [heat_transfer_CuCrZr]
        type = GenericConstantMaterial<<<{"description": "Declares material properties based on names and values prescribed by input parameters.", "href": "../../source/materials/GenericConstantMaterial.html"}>>>
        prop_names<<<{"description": "The names of the properties this material will have"}>>> = 'density_CuCrZr specific_heat_CuCrZr'
        prop_values<<<{"description": "The values associated with the named properties"}>>> = '8900.0 390.0'            # [g/m^3], [ W/m-K], [J/kg-K]
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    []
    [thermal_conductivity_CuCrZr]
        type = ParsedMaterial<<<{"description": "Parsed expression Material.", "href": "../../source/materials/ParsedMaterial.html"}>>>
        property_name<<<{"description": "Name of the parsed material property"}>>> = thermal_conductivity_CuCrZr
        coupled_variables<<<{"description": "Vector of variables used in the parsed function"}>>> = 'temperature'
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
        # expression = '5.3e-7 * temperature^3 - 6.5e-4 * temperature^2 + 2.6e-1 * temperature + 3.1e2'    # ~ 320.0  [ W/m-K] from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression<<<{"description": "Parsed function (see FParser) expression for the parsed material"}>>> = '3.87e2 - 1.28e-1 * temperature'    # ~ 349 [ W/m-K]
        outputs<<<{"description": "Vector of output names where you would like to restrict the output of variables(s) associated with this object"}>>> = all
    []
    ############################## Materials for others
    [interface_jump_4to3]
        type = SolubilityRatioMaterial<<<{"description": "Calculates the jump in concentration across an interface.", "href": "../../source/materials/SolubilityRatioMaterial.html"}>>>
        solubility_primary<<<{"description": "The material property on the primary side of the interface"}>>> = solubility_W
        solubility_secondary<<<{"description": "The material property on the secondary side of the interface"}>>> = solubility_Cu
        boundary<<<{"description": "The list of boundaries (ids or names) from the mesh where this object applies"}>>> = '4to3'
        concentration_primary<<<{"description": "Primary side non-linear variable for jump computation"}>>> = C_mobile_W
        concentration_secondary<<<{"description": "Secondary side non-linear variable for jump computation"}>>> = C_mobile_Cu
    []
    [interface_jump_3to2]
        type = SolubilityRatioMaterial<<<{"description": "Calculates the jump in concentration across an interface.", "href": "../../source/materials/SolubilityRatioMaterial.html"}>>>
        solubility_primary<<<{"description": "The material property on the primary side of the interface"}>>> = solubility_Cu
        solubility_secondary<<<{"description": "The material property on the secondary side of the interface"}>>> = solubility_CuCrZr
        boundary<<<{"description": "The list of boundaries (ids or names) from the mesh where this object applies"}>>> = '3to2'
        concentration_primary<<<{"description": "Primary side non-linear variable for jump computation"}>>> = C_mobile_Cu
        concentration_secondary<<<{"description": "Secondary side non-linear variable for jump computation"}>>> = C_mobile_CuCrZr
    []
[]

[Postprocessors<<<{"href": "../../syntax/Postprocessors/index.html"}>>>]
    ############################################################ Postprocessors for W (block = 4)
    [F_recombination]
        type = SideDiffusiveFluxAverage<<<{"description": "Computes the integral of the diffusive flux over the specified boundary", "href": "../../source/postprocessors/SideDiffusiveFluxAverage.html"}>>>
        boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = 'top'
        diffusivity<<<{"description": "The name of the diffusivity material property that will be used in the flux computation. This must be provided if the variable is of finite element type"}>>> = 5.01e-24   # (3.01604928)/(6.02e23)/[gram(T)/m^2]
        # diffusivity = 5.508e-19   # (1.0e3)*(1.0e3)/(6.02e23)/(3.01604928) [gram(T)/m^2]
        variable<<<{"description": "The name of the variable which this postprocessor integrates"}>>> = Sc_C_total_W
    []
    [F_permeation]
        type = SideDiffusiveFluxAverage<<<{"description": "Computes the integral of the diffusive flux over the specified boundary", "href": "../../source/postprocessors/SideDiffusiveFluxAverage.html"}>>>
        boundary<<<{"description": "The list of boundary IDs from the mesh where this object applies"}>>> = '2to1'
        diffusivity<<<{"description": "The name of the diffusivity material property that will be used in the flux computation. This must be provided if the variable is of finite element type"}>>> = 5.01e-24   # (3.01604928)/(6.02e23)/[gram(T)/m^2]
        # diffusivity = 5.508e-19   # (1.0e3)*(1.0e3)/(6.02e23)/(3.01604928) [gram(T)/m^2]
        variable<<<{"description": "The name of the variable which this postprocessor integrates"}>>> = Sc_C_total_CuCrZr
    []

    [Int_C_mobile_W]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_mobile_W
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    []
    [ScInt_C_mobile_W]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> =  Int_C_mobile_W
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_W]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_trapped_W
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    []
    [ScInt_C_trapped_W]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> = Int_C_trapped_W
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_W]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_total_W
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 4
    []
    [ScInt_C_total_W]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> = Int_C_total_W
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    # ############################################################ Postprocessors for Cu (block = 3)
    [Int_C_mobile_Cu]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_mobile_Cu
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    []
    [ScInt_C_mobile_Cu]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> =  Int_C_mobile_Cu
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_Cu]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_trapped_Cu
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    []
    [ScInt_C_trapped_Cu]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> = Int_C_trapped_Cu
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.44e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_Cu]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_total_Cu
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 3
    []
    [ScInt_C_total_Cu]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> = Int_C_total_Cu
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    # ############################################################ Postprocessors for CuCrZr (block = 2)
    [Int_C_mobile_CuCrZr]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_mobile_CuCrZr
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    []
    [ScInt_C_mobile_CuCrZr]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> =  Int_C_mobile_CuCrZr
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_CuCrZr]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_trapped_CuCrZr
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    []
    [ScInt_C_trapped_CuCrZr]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> = Int_C_trapped_CuCrZr
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.44e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_CuCrZr]
        type = ElementIntegralVariablePostprocessor<<<{"description": "Computes a volume integral of the specified variable", "href": "../../source/postprocessors/ElementIntegralVariablePostprocessor.html"}>>>
        variable<<<{"description": "The name of the variable that this object operates on"}>>> = C_total_CuCrZr
        block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = 2
    []
    [ScInt_C_total_CuCrZr]
        type = ScalePostprocessor<<<{"description": "Scales a post-processor by a value", "href": "../../source/postprocessors/ScalePostprocessor.html"}>>>
        value<<<{"description": "The postprocessor to be scaled"}>>> = Int_C_total_CuCrZr
        scaling_factor<<<{"description": "The scaling factor"}>>> = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    ############################################################ Postprocessors for others
    [dt]
        type = TimestepSize<<<{"description": "Reports the timestep size", "href": "../../source/postprocessors/TimestepSize.html"}>>>
    []
    [temperature_top]
        type = PointValue<<<{"description": "Compute the value of a variable at a specified location", "href": "../../source/postprocessors/PointValue.html"}>>>
        variable<<<{"description": "The name of the variable that this postprocessor operates on."}>>> = temperature
        point<<<{"description": "The physical point where the solution will be evaluated."}>>> = '0 14.0e-3 0'
    []
    [temperature_tube]
        type = PointValue<<<{"description": "Compute the value of a variable at a specified location", "href": "../../source/postprocessors/PointValue.html"}>>>
        variable<<<{"description": "The name of the variable that this postprocessor operates on."}>>> = temperature
        point<<<{"description": "The physical point where the solution will be evaluated."}>>> = '0 6.0e-3 0'
    []
    # limit timestep
    [timestep_max_pp] # s
    type = FunctionValuePostprocessor<<<{"description": "Computes the value of a supplied function at a single point (scalable)", "href": "../../source/postprocessors/FunctionValuePostprocessor.html"}>>>
    function<<<{"description": "The function which supplies the postprocessor value."}>>> = timestep_function
  []
[]

[VectorPostprocessors<<<{"href": "../../syntax/VectorPostprocessors/index.html"}>>>]
    [line]
        type = LineValueSampler<<<{"description": "Samples variable(s) along a specified line", "href": "../../source/vectorpostprocessors/LineValueSampler.html"}>>>
        start_point<<<{"description": "The beginning of the line"}>>> = '0 14.0e-3 0'
        end_point<<<{"description": "The ending of the line"}>>> = '0 6.0e-3 0'
        num_points<<<{"description": "The number of points to sample along the line"}>>> = 100
        sort_by<<<{"description": "What to sort the samples by"}>>> = 'y'
        variable<<<{"description": "The names of the variables that this VectorPostprocessor operates on"}>>> = 'C_total_W C_total_Cu C_total_CuCrZr C_mobile_W C_mobile_Cu C_mobile_CuCrZr C_trapped_W C_trapped_Cu C_trapped_CuCrZr flux_y temperature'
        execute_on<<<{"description": "The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html."}>>> = timestep_end
    []
[]

[Preconditioning<<<{"href": "../../syntax/Preconditioning/index.html"}>>>]
    [smp]
        type = SMP<<<{"description": "Single matrix preconditioner (SMP) builds a preconditioner using user defined off-diagonal parts of the Jacobian.", "href": "../../source/preconditioners/SingleMatrixPreconditioner.html"}>>>
        full<<<{"description": "Set to true if you want the full set of couplings between variables simply for convenience so you don't have to set every off_diag_row and off_diag_column combination."}>>> = true
    []
[]

[Executioner<<<{"href": "../../syntax/Executioner/index.html"}>>>]
    type = Transient
    scheme = bdf2
    solve_type = NEWTON
    petsc_options_iname = '-pc_type'
    petsc_options_value = 'lu'
    nl_rel_tol  = 1e-6 # 1e-8 works for 1 cycle
    nl_abs_tol  = 1e-7 # 1e-11 works for 1 cycle
    end_time = 8.0e4   # 50 ITER shots (3.0e4 s plasma, 2.0e4 SSP)
    automatic_scaling = true
    line_search = 'none'
    dtmin = 1e-4
    nl_max_its = 18
    [TimeStepper<<<{"href": "../../syntax/Executioner/TimeStepper/index.html"}>>>]
        type = IterationAdaptiveDT
        dt = 20
        optimal_iterations = 15
        iteration_window = 1
        growth_factor = 1.2
        cutback_factor = 0.8
        timestep_limiting_postprocessor = timestep_max_pp
    []
[]

[Outputs<<<{"href": "../../syntax/Outputs/index.html"}>>>]
    [exodus]
        type = Exodus<<<{"description": "Object for output data in the Exodus format", "href": "../../source/outputs/Exodus.html"}>>>
        sync_only<<<{"description": "Only export results at sync times"}>>> = false
        # output at key moment in the first two cycles, and then at the end of the simulation
        sync_times<<<{"description": "Times at which the output and solution is forced to occur"}>>> = '110.0 480.0 590.0 1600.0 1710.0 2080.0 2190.0 3400.0 8.0e4'
    []
    csv<<<{"description": "Output the scalar variable and postprocessors to a *.csv file using the default CSV output."}>>> = true
    hide<<<{"description": "A list of the variables and postprocessors that should NOT be output to the Exodus file (may include Variables, ScalarVariables, and Postprocessor names)."}>>> = 'dt
            Int_C_mobile_W Int_C_trapped_W Int_C_total_W
            Int_C_mobile_Cu Int_C_trapped_Cu Int_C_total_Cu
            Int_C_mobile_CuCrZr Int_C_trapped_CuCrZr Int_C_total_CuCrZr'
    perf_graph<<<{"description": "Enable printing of the performance graph to the screen (Console)"}>>> = true
[]

References

E. A. Hodille, R. Delaporte-Mathurin, J. Denis, M. Pecovnik, E. Bernard, Y. Ferro, R. Sakamoto, Y. Charles, J. Mougenot, A. De Backer, C. S. Becquart, S. Markelj, and C. Grisolia. Modelling of hydrogen isotopes trapping, diffusion and permeation in divertor monoblocks under ITER-like conditions. Nuclear Fusion, 61:126003, 10 2021. URL: https://iopscience.iop.org/article/10.1088/1741-4326/ac2abc, doi:10.1088/1741-4326/AC2ABC.

@article{Hodille2021126003,
    author = "Hodille, E. A. and Delaporte-Mathurin, R. and Denis, J. and Pecovnik, M. and Bernard, E. and Ferro, Y. and Sakamoto, R. and Charles, Y. and Mougenot, J. and Backer, A. De and Becquart, C. S. and Markelj, S. and Grisolia, C.",
    doi = "10.1088/1741-4326/AC2ABC",
    issn = "0029-5515",
    issue = "12",
    journal = "Nuclear Fusion",
    month = "10",
    pages = "126003",
    publisher = "IOP Publishing",
    title = "Modelling of hydrogen isotopes trapping, diffusion and permeation in divertor monoblocks under {ITER}-like conditions",
    volume = "61",
    url = "https://iopscience.iop.org/article/10.1088/1741-4326/ac2abc",
    year = "2021"
}

Masashi Shimada, Pierre-Clément A. Simon, Casey T. Icenhour, and Gyanender Singh. Toward a high-fidelity tritium transport modeling for retention and permeation experiments. Fusion Engineering and Design, 203:114438, 6 2024. URL: https://linkinghub.elsevier.com/retrieve/pii/S0920379624002916, doi:10.1016/J.FUSENGDES.2024.114438.

@article{Shimada2024114438,
    author = "Shimada, Masashi and Simon, Pierre-Clément A. and Icenhour, Casey T. and Singh, Gyanender",
    doi = "10.1016/J.FUSENGDES.2024.114438",
    issn = "0920-3796",
    journal = "Fusion Engineering and Design",
    month = "6",
    pages = "114438",
    publisher = "North-Holland",
    title = "Toward a high-fidelity tritium transport modeling for retention and permeation experiments",
    volume = "203",
    url = "https://linkinghub.elsevier.com/retrieve/pii/S0920379624002916",
    year = "2024"
}

(test/tests/divertor_monoblock/divertor_monoblock.i)

### Nomenclatures
### C_mobile_j      mobile H concentration in "j" material, where j = CuCrZr, Cu, W
### C_trapped_j     trapped H concentration in "j" material, where j = CuCrZr, Cu, W
### C_total_j       total H concentration in "j" material, where j = CuCrZr, Cu, W
###
### S_empty_j       empty site concentration in "j" material, where j = CuCrZr, Cu, W
### S_trapped_j     trapped site concentration in "j" material, where j = CuCrZr, Cu, W
### S_total_j       total site H concentration in "j" material, where j = CuCrZr, Cu, W
###
### F_permeation    permeation flux
### F_recombination recombination flux
###
### Sc_             Scaled
### Int_            Integrated
### ScInt_          Scaled and integrated

tungsten_atomic_density = ${units 6.338e28 m^-3}

## Mesh generated by ConcentricCircleMeshGenerator and saved as "2DMonoblock.e"
[Mesh]
    [ccmg]
        type = ConcentricCircleMeshGenerator
        num_sectors = 36
        rings = '1 30 20 110'
        radii = '${units 6 mm -> m} ${units 7.5 mm -> m} ${units 8.5 mm -> m}'
        has_outer_square = on
        pitch = ${units 28 mm -> m}
        portion = left_half
        preserve_volumes = false
        smoothing_max_it = 3
    []
    [ssbsg1]
        type = SideSetsBetweenSubdomainsGenerator
        input = ccmg
        primary_block = '4'     # W
        paired_block = '3'      # Cu
        new_boundary = '4to3'
    []
    [ssbsg2]
        type = SideSetsBetweenSubdomainsGenerator
        input = ssbsg1
        primary_block = '3'     # Cu
        paired_block = '4'      # W
        new_boundary = '3to4'
    []
    [ssbsg3]
        type = SideSetsBetweenSubdomainsGenerator
        input = ssbsg2
        primary_block = '3'     # Cu
        paired_block = '2'      # CuCrZr
        new_boundary = '3to2'
    []
    [ssbsg4]
        type = SideSetsBetweenSubdomainsGenerator
        input = ssbsg3
        primary_block = '2'     # CuCrZr
        paired_block = '3'      # Cu
        new_boundary = '2to3'
    []
    [ssbsg5]
        type = SideSetsBetweenSubdomainsGenerator
        input = ssbsg4
        primary_block = '2'     # CuCrZr
        paired_block = '1'      # H2O
        new_boundary = '2to1'
    []
    [bdg]
        type = BlockDeletionGenerator
        input = ssbsg5
        block = '1'             # H2O
    []
[]

[Problem]
    type = ReferenceResidualProblem
    extra_tag_vectors = 'ref'
    reference_vector = 'ref'
[]

[Variables]
    [temperature]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 300 K}
    []
    ######################### Variables for W (block = 4)
    [C_mobile_W]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-20 m^-3}
        block = 4
    []
    [C_trapped_W]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-15 m^-3}
        block = 4
    []
    ######################### Variables for Cu (block = 3)
    [C_mobile_Cu]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 5.0e-17 m^-3}
        block = 3
    []
    [C_trapped_Cu]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-15 m^-3}
        block = 3
    []
    ######################### Variables for CuCrZr (block = 2)
    [C_mobile_CuCrZr]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-15 m^-3}
        block = 2
    []
    [C_trapped_CuCrZr]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-15 m^-3}
        block = 2
    []
[]

[AuxVariables]
    [flux_y]
        order = FIRST
        family = MONOMIAL
    []
    ############################## AuxVariables for W (block = 4)
    [Sc_C_mobile_W]
        block = 4
    []
    [Sc_C_trapped_W]
        block = 4
    []
    [C_total_W]
        block = 4
    []
    [Sc_C_total_W]
        block = 4
    []
    [S_empty_W]
        block = 4
    []
    [Sc_S_empty_W]
        block = 4
    []
    [S_trapped_W]
        block = 4
    []
    [Sc_S_trapped_W]
        block = 4
    []
    [S_total_W]
        block = 4
    []
    [Sc_S_total_W]
        block = 4
    []
    ############################## AuxVariables for Cu (block = 3)
    [Sc_C_mobile_Cu]
        block = 3
    []
    [Sc_C_trapped_Cu]
        block = 3
    []
    [C_total_Cu]
        block = 3
    []
    [Sc_C_total_Cu]
        block = 3
    []
    [S_empty_Cu]
        block = 3
    []
    [Sc_S_empty_Cu]
        block = 3
    []
    [S_trapped_Cu]
        block = 3
    []
    [Sc_S_trapped_Cu]
        block = 3
    []
    [S_total_Cu]
        block = 3
    []
    [Sc_S_total_Cu]
        block = 3
    []
    ############################## AuxVariables for CuCrZr (block = 2)
    [Sc_C_mobile_CuCrZr]
        block = 2
    []
    [Sc_C_trapped_CuCrZr]
        block = 2
    []
    [C_total_CuCrZr]
        block = 2
    []
    [Sc_C_total_CuCrZr]
        block = 2
    []
    [S_empty_CuCrZr]
        block = 2
    []
    [Sc_S_empty_CuCrZr]
        block = 2
    []
    [S_trapped_CuCrZr]
        block = 2
    []
    [Sc_S_trapped_CuCrZr]
        block = 2
    []
    [S_total_CuCrZr]
        block = 2
    []
    [Sc_S_total_CuCrZr]
        block = 2
    []
[]

[Kernels]
    ############################## Kernels for W (block = 4)
    [diff_W]
        type = ADMatDiffusion
        variable = C_mobile_W
        diffusivity = diffusivity_W
        block = 4
        extra_vector_tags = ref
    []
    [time_diff_W]
        type = ADTimeDerivative
        variable = C_mobile_W
        block = 4
        extra_vector_tags = ref
    []
    [coupled_time_W]
        type = ScaledCoupledTimeDerivative
        variable = C_mobile_W
        v = C_trapped_W
        factor = 1e0
        block = 4
        extra_vector_tags = ref
    []
    [heat_conduction_W]
        type = HeatConduction
        variable = temperature
        diffusion_coefficient = thermal_conductivity_W
        block = 4
        extra_vector_tags = ref
    []
    [time_heat_conduction_W]
        type = SpecificHeatConductionTimeDerivative
        variable = temperature
        specific_heat = specific_heat_W
        density = density_W
        block = 4
        extra_vector_tags = ref
    []
    ############################## Kernels for Cu (block = 3)
    [diff_Cu]
        type = ADMatDiffusion
        variable = C_mobile_Cu
        diffusivity = diffusivity_Cu
        block = 3
        extra_vector_tags = ref
    []
    [time_diff_Cu]
        type = ADTimeDerivative
        variable = C_mobile_Cu
        block = 3
        extra_vector_tags = ref
    []
    [coupled_time_Cu]
        type = ScaledCoupledTimeDerivative
        variable = C_mobile_Cu
        v = C_trapped_Cu
        factor = 1e0
        block = 3
        extra_vector_tags = ref
    []
    [heat_conduction_Cu]
        type = HeatConduction
        variable = temperature
        diffusion_coefficient = thermal_conductivity_Cu
        block = 3
        extra_vector_tags = ref
    []
    [time_heat_conduction_Cu]
        type = SpecificHeatConductionTimeDerivative
        variable = temperature
        specific_heat = specific_heat_Cu
        density = density_Cu
        block = 3
        extra_vector_tags = ref
    []
    ############################## Kernels for CuCrZr (block = 2)
    [diff_CuCrZr]
        type = ADMatDiffusion
        variable = C_mobile_CuCrZr
        diffusivity = diffusivity_CuCrZr
        block = 2
        extra_vector_tags = ref
    []
    [time_diff_CuCrZr]
        type = ADTimeDerivative
        variable = C_mobile_CuCrZr
        block = 2
        extra_vector_tags = ref
    []
    [coupled_time_CuCrZr]
        type = ScaledCoupledTimeDerivative
        variable = C_mobile_CuCrZr
        v = C_trapped_CuCrZr
        factor = 1e0
        block = 2
        extra_vector_tags = ref
    []
    [heat_conduction_CuCrZr]
        type = HeatConduction
        variable = temperature
        diffusion_coefficient = thermal_conductivity_CuCrZr
        block = 2
        extra_vector_tags = ref
    []
    [time_heat_conduction_CuCrZr]
        type = SpecificHeatConductionTimeDerivative
        variable = temperature
        specific_heat = specific_heat_CuCrZr
        density = density_CuCrZr
        block = 2
        extra_vector_tags = ref
    []
[]

[AuxKernels]
    ############################## AuxKernels for W (block = 4)
    [Scaled_mobile_W]
        variable = Sc_C_mobile_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_mobile_W
    []
    [Scaled_trapped_W]
        variable = Sc_C_trapped_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_trapped_W
    []
    [total_W]
        variable = C_total_W
        type = ParsedAux
        expression = 'C_mobile_W + C_trapped_W'
        coupled_variables = 'C_mobile_W C_trapped_W'
    []
    [Scaled_total_W]
        variable = Sc_C_total_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_total_W
    []
    [empty_sites_W]
        variable = S_empty_W
        type = EmptySitesAux
        N = ${units 1.0e0 m^-3}       # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        # Ct0 = ${units 1.0e-4 m^-3}   # E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        Ct0 = ${units 1.0e-4 m^-3}    # E.A. Hodille et al 2021 Nucl. Fusion 61 1260033, trap 2
        trap_per_free = 1.0e0         # 1.0e1
        trapped_concentration_variables = C_trapped_W
    []
    [scaled_empty_W]
        variable = Sc_S_empty_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_empty_W
    []
    [trapped_sites_W]
        variable = S_trapped_W
        type = NormalizationAux
        normal_factor = 1e0
        source_variable = C_trapped_W
    []
    [scaled_trapped_W]
        variable = Sc_S_trapped_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_trapped_W
    []
    [total_sites_W]
        variable = S_total_W
        type = ParsedAux
        expression = 'S_trapped_W + S_empty_W'
        coupled_variables = 'S_trapped_W S_empty_W'
    []
    [scaled_total_W]
        variable = Sc_S_total_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_total_W
    []
    ############################## AuxKernels for Cu (block = 3)
    [Scaled_mobile_Cu]
        variable = Sc_C_mobile_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_mobile_Cu
    []
    [Scaled_trapped_Cu]
        variable = Sc_C_trapped_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_trapped_Cu
    []
    [total_Cu]
        variable = C_total_Cu
        type = ParsedAux
        expression = 'C_mobile_Cu + C_trapped_Cu'
        coupled_variables = 'C_mobile_Cu C_trapped_Cu'
    []
    [Scaled_total_Cu]
        variable = Sc_C_total_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_total_Cu
    []
    [empty_sites_Cu]
        variable = S_empty_Cu
        type = EmptySitesAux
        N = ${units 1.0e0 m^-3}     # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0 = ${units 5.0e-5 m^-3}  # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        trap_per_free = 1.0e0       # 1.0e1
        trapped_concentration_variables = C_trapped_Cu
    []
    [scaled_empty_Cu]
        variable = Sc_S_empty_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_empty_Cu
    []
    [trapped_sites_Cu]
        variable = S_trapped_Cu
        type = NormalizationAux
        normal_factor = 1e0
        source_variable = C_trapped_Cu
    []
    [scaled_trapped_Cu]
        variable = Sc_S_trapped_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_trapped_Cu
    []
    [total_sites_Cu]
        variable = S_total_Cu
        type = ParsedAux
        expression = 'S_trapped_Cu + S_empty_Cu'
        coupled_variables = 'S_trapped_Cu S_empty_Cu'
    []
    [scaled_total_Cu]
        variable = Sc_S_total_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_total_Cu
    []
    ############################## AuxKernels for CuCrZr (block = 2)
    [Scaled_mobile_CuCrZr]
        variable = Sc_C_mobile_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_mobile_CuCrZr
    []
    [Scaled_trapped_CuCrZr]
        variable = Sc_C_trapped_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_trapped_CuCrZr
    []
    [total_CuCrZr]
        variable = C_total_CuCrZr
        type = ParsedAux
        expression = 'C_mobile_CuCrZr + C_trapped_CuCrZr'
        coupled_variables = 'C_mobile_CuCrZr C_trapped_CuCrZr'
    []
    [Scaled_total_CuCrZr]
        variable = Sc_C_total_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_total_CuCrZr
    []
    [empty_sites_CuCrZr]
        variable = S_empty_CuCrZr
        type = EmptySitesAux
        N = ${units 1.0e0 m^-3}     # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0 = ${units 5.0e-5 m^-3}  # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # Ct0 = ${units 4.0e-2 m^-3} # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        trap_per_free = 1.0e0       # 1.0e1
        trapped_concentration_variables = C_trapped_CuCrZr
    []
    [scaled_empty_CuCrZr]
        variable = Sc_S_empty_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_empty_CuCrZr
    []
    [trapped_sites_CuCrZr]
        variable = S_trapped_CuCrZr
        type = NormalizationAux
        normal_factor = 1e0
        source_variable = C_trapped_CuCrZr
    []
    [scaled_trapped_CuCrZr]
        variable = Sc_S_trapped_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_trapped_CuCrZr
    []
    [total_sites_CuCrZr]
        variable = S_total_CuCrZr
        type = ParsedAux
        expression = 'S_trapped_CuCrZr + S_empty_CuCrZr'
        coupled_variables = 'S_trapped_CuCrZr S_empty_CuCrZr'
    []
    [scaled_total_CuCrZr]
        variable = Sc_S_total_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_total_CuCrZr
    []
    [flux_y_W]
        type = DiffusionFluxAux
        diffusivity = diffusivity_W
        variable = flux_y
        diffusion_variable = C_mobile_W
        component = y
        block = 4
    []
    [flux_y_Cu]
        type = DiffusionFluxAux
        diffusivity = diffusivity_Cu
        variable = flux_y
        diffusion_variable = C_mobile_Cu
        component = y
        block = 3
    []
    [flux_y_CuCrZr]
        type = DiffusionFluxAux
        diffusivity = diffusivity_CuCrZr
        variable = flux_y
        diffusion_variable = C_mobile_CuCrZr
        component = y
        block = 2
    []
[]

[InterfaceKernels]
    [tied_4to3]
        type = ADPenaltyInterfaceDiffusion
        variable = C_mobile_W
        neighbor_var = C_mobile_Cu
        penalty = 0.05            #  it will not converge with > 0.1, but it creates negative C_mobile _Cu with << 0.1
        # jump_prop_name = solubility_ratio_4to3
        jump_prop_name = solubility_ratio
        boundary = '4to3'
    []
    [tied_3to2]
        type = ADPenaltyInterfaceDiffusion
        variable = C_mobile_Cu
        neighbor_var = C_mobile_CuCrZr
        penalty = 0.05            #  it will not converge with > 0.1, but it creates negative C_mobile _Cu with << 0.1
        # jump_prop_name = solubility_ratio_3to2
        jump_prop_name = solubility_ratio
        boundary = '3to2'
    []
[]

[NodalKernels]
    ############################## NodalKernels for W (block = 4)
    [time_W]
        type = TimeDerivativeNodalKernel
        variable = C_trapped_W
    []
    [trapping_W]
        type = TrappingNodalKernel
        variable = C_trapped_W
        temperature = temperature
        alpha_t = 2.75e11      # 1e15
        N = 1.0e0  # = (1e0) x (${tungsten_atomic_density} #/m^3)
        # Ct0 = 1.0e-4                # E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        Ct0 = 1.0e-4                # E.A. Hodille et al 2021 Nucl. Fusion 61 1260033, trap 2
        trap_per_free = 1.0e0       # 1.0e1
        mobile_concentration = 'C_mobile_W'
        extra_vector_tags = ref
    []
    [release_W]
        type = ReleasingNodalKernel
        alpha_r = 8.4e12    # 1.0e13
        temperature = temperature
        # detrapping_energy = 9863.9    # = 0.85 eV    E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        detrapping_energy = 11604.6   # = 1.00 eV    E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 2
        variable = C_trapped_W
    []
    ############################## NodalKernels for Cu (block = 3)
    [time_Cu]
        type = TimeDerivativeNodalKernel
        variable = C_trapped_Cu
    []
    [trapping_Cu]
        type = TrappingNodalKernel
        variable = C_trapped_Cu
        temperature = temperature
        alpha_t = 2.75e11      # 1e15
        N = 1.0e0  # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0 = 5.0e-5                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        trap_per_free = 1.0e0       # 1.0e1
        mobile_concentration = 'C_mobile_Cu'
        extra_vector_tags = ref
    []
    [release_Cu]
        type = ReleasingNodalKernel
        alpha_r = 8.4e12    # 1.0e13
        temperature = temperature
        detrapping_energy = 5802.3    # = 0.50eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        variable = C_trapped_Cu
    []
    ############################## NodalKernels for CuCrZr (block = 2)
    [time_CuCrZr]
        type = TimeDerivativeNodalKernel
        variable = C_trapped_CuCrZr
    []
    [trapping_CuCrZr]
        type = TrappingNodalKernel
        variable = C_trapped_CuCrZr
        temperature = temperature
        alpha_t = 2.75e11      # 1e15
        N = 1.0e0  # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0 = 5.0e-5                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # Ct0 = 4.0e-2                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        trap_per_free = 1.0e0       # 1.0e1
        mobile_concentration = 'C_mobile_CuCrZr'
        extra_vector_tags = ref
    []
    [release_CuCrZr]
        type = ReleasingNodalKernel
        alpha_r = 8.4e12    # 1.0e13
        temperature = temperature
        detrapping_energy = 5802.3    # = 0.50eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # detrapping_energy = 9631.8   # = 0.83 eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        variable = C_trapped_CuCrZr
    []
[]

[BCs]
    [C_mob_W_top_flux]
        type = FunctionNeumannBC
        variable = C_mobile_W
        boundary = 'top'
        function = mobile_flux_bc_func
    []
    [mobile_tube]
        type = DirichletBC
        variable = C_mobile_CuCrZr
        value = 1.0e-18
        boundary = '2to1'
    []
    [temp_top]
        type = FunctionNeumannBC
        variable = temperature
        boundary = 'top'
        function = temp_flux_bc_func
    []
    [temp_tube]
        type = FunctionDirichletBC
        variable = temperature
        boundary = '2to1'
        function = temp_inner_func
    []
[]

[Functions]
    ### Maximum mobile flux of 7.90e-13 at the top surface (1.0e-4 [m])
    ### 1.80e23/m^2-s = (5.0e23/m^2-s) *(1-0.999) = (7.90e-13)*(${tungsten_atomic_density})/(1.0e-4)  at steady state
    [mobile_flux_bc_func]
        type = ParsedFunction
        expression =   'if((t % 1600) < 100.0, 0.0      + 7.90e-13*(t % 1600)/100,
                        if((t % 1600) < 500.0, 7.90e-13,
                        if((t % 1600) < 600.0, 7.90e-13 - 7.90e-13*((t % 1600)-500)/100, 0.0)))'
    []
    ### Heat flux of 10MW/m^2 at steady state
    [temp_flux_bc_func]
        type = ParsedFunction
        expression =   'if((t % 1600) < 100.0, 0.0   + 1.0e7*(t % 1600)/100,
                        if((t % 1600) < 500.0, 1.0e7,
                        if((t % 1600) < 600.0, 1.0e7 - 1.0e7*((t % 1600)-500)/100, 300)))'
    []
    ### Maximum coolant temperature of 552K at steady state
    [temp_inner_func]
        type = ParsedFunction
        expression =   'if((t % 1600) < 100.0, 300.0 + (552-300)*(t % 1600)/100,
                        if((t % 1600) < 500.0, 552,
                        if((t % 1600) < 600.0, 552.0 - (552-300)*((t % 1600)-500)/100, 300)))'
    []
    [timestep_function]
        type = ParsedFunction
        expression = 'if((t % 1600) <   10.0,  20,
                      if((t % 1600) <   90.0,  40,
                      if((t % 1600) <  110.0,  20,
                      if((t % 1600) <  480.0,  40,
                      if((t % 1600) <  500.0,  20,
                      if((t % 1600) <  590.0,   4,
                      if((t % 1600) <  610.0,  20,
                      if((t % 1600) < 1500.0, 200,
                      if((t % 1600) < 1600.0,  40,  2)))))))))'
    []
[]

[Materials]
    ############################## Materials for W (block = 4)
    [diffusivity_W]
        type = ADParsedMaterial
        property_name = diffusivity_W
        coupled_variables = 'temperature'
        block = 4
        expression = '2.4e-7*exp(-4525.8/temperature)'    # H diffusivity in W
        outputs = all
    []
    [solubility_W]
        type = ADParsedMaterial
        property_name = solubility_W
        coupled_variables = 'temperature'
        block = 4
        # expression = '2.95e-5 *exp(-12069.0/temperature)'              # H solubility in W = (1.87e24)/(${tungsten_atomic_density}) [#/m^3]
        expression = '2.95e-5 *exp(-12069.0/temperature) + 4.95e-8 * exp(-6614.6/temperature)'    # H solubility in W = (1.87e24)/(${tungsten_atomic_density}) [#/m^3]
        outputs = all
    []
    [converter_to_regular_W]
        type = MaterialADConverter
        ad_props_in = 'diffusivity_W'
        reg_props_out = 'diffusivity_W_nonAD'
        block = 4
    []
    [heat_transfer_W]
        type = GenericConstantMaterial
        prop_names = 'density_W'
        prop_values = '19300'                # [g/m^3]
        block = 4
    []
    [specific_heat_W]
        type = ParsedMaterial
        property_name = specific_heat_W
        coupled_variables = 'temperature'
        block = 4
        expression = '1.16e2 + 7.11e-2 * temperature - 6.58e-5 * temperature^2 + 3.24e-8 * temperature^3 -5.45e-12 * temperature^4'    # ~ 132[J/kg-K]
        outputs = all
    []
    [thermal_conductivity_W]
        type = ParsedMaterial
        property_name = thermal_conductivity_W
        coupled_variables = 'temperature'
        block = 4
        # expression = '-7.8e-9 * temperature^3 + 5.0e-5 * temperature^2 - 1.1e-1 * temperature + 1.8e2'    # ~ 173.0 [ W/m-K]   from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression = '2.41e2 - 2.90e-1 * temperature + 2.54e-4 * temperature^2 - 1.03e-7 * temperature^3 + 1.52e-11 * temperature^4'    # ~ 173.0 [ W/m-K]
        outputs = all
    []
    ############################## Materials for Cu (block = 3)
    [diffusivity_Cu]
        type = ADParsedMaterial
        property_name = diffusivity_Cu
        coupled_variables = 'temperature'
        block = 3
        expression = '6.60e-7*exp(-4525.8/temperature)'    # H diffusivity in Cu
        outputs = all
    []
    [solubility_Cu]
        type = ADParsedMaterial
        property_name = solubility_Cu
        coupled_variables = 'temperature'
        block = 3
        expression = '4.95e-5*exp(-6614.6/temperature)'    # H solubility in Cu = (3.14e24)/(${tungsten_atomic_density}) [#/m^3]
        outputs = all
    []
    [converter_to_regular_Cu]
        type = MaterialADConverter
        ad_props_in = 'diffusivity_Cu'
        reg_props_out = 'diffusivity_Cu_nonAD'
        block = 3
    []
    [heat_transfer_Cu]
        type = GenericConstantMaterial
        prop_names = 'density_Cu'
        prop_values = '8960.0'                # [g/m^3]
        block = 3
    []
    [specific_heat_Cu]
        type = ParsedMaterial
        property_name = specific_heat_Cu
        coupled_variables = 'temperature'
        block = 3
        expression = '3.16e2 + 3.18e-1 * temperature - 3.49e-4 * temperature^2 + 1.66e-7 * temperature^3'    # ~ 384 [J/kg-K]
        outputs = all
    []
    [thermal_conductivity_Cu]
        type = ParsedMaterial
        property_name = thermal_conductivity_Cu
        coupled_variables = 'temperature'
        block = 3
        # expression = '-3.9e-8 * temperature^3 + 3.8e-5 * temperature^2 - 7.9e-2 * temperature + 4.0e2'    # ~ 401.0  [ W/m-K] from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression = '4.21e2 - 6.85e-2 * temperature'    # ~ 400.0 [ W/m-K]
        outputs = all
    []
    ############################## Materials for CuCrZr (block = 2)
    [diffusivity_CuCrZr]
        type = ADParsedMaterial
        property_name = diffusivity_CuCrZr
        coupled_variables = 'temperature'
        block = 2
        expression = '3.90e-7*exp(-4873.9/temperature)'    # H diffusivity in CuCrZr
        outputs = all
    []
    [solubility_CuCrZr]
        type = ADParsedMaterial
        property_name = solubility_CuCrZr
        coupled_variables = 'temperature'
        block = 2
        expression = '6.75e-6*exp(-4525.8/temperature)'    # H solubility in CuCrZr = (4.28e23)/(${tungsten_atomic_density}) [#/m^3]
        outputs = all
    []
    [converter_to_regular_CuCrZr]
        type = MaterialADConverter
        ad_props_in = 'diffusivity_CuCrZr'
        reg_props_out = 'diffusivity_CuCrZr_nonAD'
        block = 2
    []
    [heat_transfer_CuCrZr]
        type = GenericConstantMaterial
        prop_names = 'density_CuCrZr specific_heat_CuCrZr'
        prop_values = '8900.0 390.0'            # [g/m^3], [ W/m-K], [J/kg-K]
        block = 2
    []
    [thermal_conductivity_CuCrZr]
        type = ParsedMaterial
        property_name = thermal_conductivity_CuCrZr
        coupled_variables = 'temperature'
        block = 2
        # expression = '5.3e-7 * temperature^3 - 6.5e-4 * temperature^2 + 2.6e-1 * temperature + 3.1e2'    # ~ 320.0  [ W/m-K] from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression = '3.87e2 - 1.28e-1 * temperature'    # ~ 349 [ W/m-K]
        outputs = all
    []
    ############################## Materials for others
    [interface_jump_4to3]
        type = SolubilityRatioMaterial
        solubility_primary = solubility_W
        solubility_secondary = solubility_Cu
        boundary = '4to3'
        concentration_primary = C_mobile_W
        concentration_secondary = C_mobile_Cu
    []
    [interface_jump_3to2]
        type = SolubilityRatioMaterial
        solubility_primary = solubility_Cu
        solubility_secondary = solubility_CuCrZr
        boundary = '3to2'
        concentration_primary = C_mobile_Cu
        concentration_secondary = C_mobile_CuCrZr
    []
[]

[Postprocessors]
    ############################################################ Postprocessors for W (block = 4)
    [F_recombination]
        type = SideDiffusiveFluxAverage
        boundary = 'top'
        diffusivity = 5.01e-24   # (3.01604928)/(6.02e23)/[gram(T)/m^2]
        # diffusivity = 5.508e-19   # (1.0e3)*(1.0e3)/(6.02e23)/(3.01604928) [gram(T)/m^2]
        variable = Sc_C_total_W
    []
    [F_permeation]
        type = SideDiffusiveFluxAverage
        boundary = '2to1'
        diffusivity = 5.01e-24   # (3.01604928)/(6.02e23)/[gram(T)/m^2]
        # diffusivity = 5.508e-19   # (1.0e3)*(1.0e3)/(6.02e23)/(3.01604928) [gram(T)/m^2]
        variable = Sc_C_total_CuCrZr
    []

    [Int_C_mobile_W]
        type = ElementIntegralVariablePostprocessor
        variable = C_mobile_W
        block = 4
    []
    [ScInt_C_mobile_W]
        type = ScalePostprocessor
        value =  Int_C_mobile_W
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_W]
        type = ElementIntegralVariablePostprocessor
        variable = C_trapped_W
        block = 4
    []
    [ScInt_C_trapped_W]
        type = ScalePostprocessor
        value = Int_C_trapped_W
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_W]
        type = ElementIntegralVariablePostprocessor
        variable = C_total_W
        block = 4
    []
    [ScInt_C_total_W]
        type = ScalePostprocessor
        value = Int_C_total_W
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    # ############################################################ Postprocessors for Cu (block = 3)
    [Int_C_mobile_Cu]
        type = ElementIntegralVariablePostprocessor
        variable = C_mobile_Cu
        block = 3
    []
    [ScInt_C_mobile_Cu]
        type = ScalePostprocessor
        value =  Int_C_mobile_Cu
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_Cu]
        type = ElementIntegralVariablePostprocessor
        variable = C_trapped_Cu
        block = 3
    []
    [ScInt_C_trapped_Cu]
        type = ScalePostprocessor
        value = Int_C_trapped_Cu
        scaling_factor = 3.44e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_Cu]
        type = ElementIntegralVariablePostprocessor
        variable = C_total_Cu
        block = 3
    []
    [ScInt_C_total_Cu]
        type = ScalePostprocessor
        value = Int_C_total_Cu
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    # ############################################################ Postprocessors for CuCrZr (block = 2)
    [Int_C_mobile_CuCrZr]
        type = ElementIntegralVariablePostprocessor
        variable = C_mobile_CuCrZr
        block = 2
    []
    [ScInt_C_mobile_CuCrZr]
        type = ScalePostprocessor
        value =  Int_C_mobile_CuCrZr
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_CuCrZr]
        type = ElementIntegralVariablePostprocessor
        variable = C_trapped_CuCrZr
        block = 2
    []
    [ScInt_C_trapped_CuCrZr]
        type = ScalePostprocessor
        value = Int_C_trapped_CuCrZr
        scaling_factor = 3.44e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_CuCrZr]
        type = ElementIntegralVariablePostprocessor
        variable = C_total_CuCrZr
        block = 2
    []
    [ScInt_C_total_CuCrZr]
        type = ScalePostprocessor
        value = Int_C_total_CuCrZr
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    ############################################################ Postprocessors for others
    [dt]
        type = TimestepSize
    []
    [temperature_top]
        type = PointValue
        variable = temperature
        point = '0 14.0e-3 0'
    []
    [temperature_tube]
        type = PointValue
        variable = temperature
        point = '0 6.0e-3 0'
    []
    # limit timestep
    [timestep_max_pp] # s
    type = FunctionValuePostprocessor
    function = timestep_function
  []
[]

[VectorPostprocessors]
    [line]
        type = LineValueSampler
        start_point = '0 14.0e-3 0'
        end_point = '0 6.0e-3 0'
        num_points = 100
        sort_by = 'y'
        variable = 'C_total_W C_total_Cu C_total_CuCrZr C_mobile_W C_mobile_Cu C_mobile_CuCrZr C_trapped_W C_trapped_Cu C_trapped_CuCrZr flux_y temperature'
        execute_on = timestep_end
    []
[]


[Preconditioning]
    [smp]
        type = SMP
        full = true
    []
[]

[Executioner]
    type = Transient
    scheme = bdf2
    solve_type = NEWTON
    petsc_options_iname = '-pc_type'
    petsc_options_value = 'lu'
    nl_rel_tol  = 1e-6 # 1e-8 works for 1 cycle
    nl_abs_tol  = 1e-7 # 1e-11 works for 1 cycle
    end_time = 8.0e4   # 50 ITER shots (3.0e4 s plasma, 2.0e4 SSP)
    automatic_scaling = true
    line_search = 'none'
    dtmin = 1e-4
    nl_max_its = 18
    [TimeStepper]
        type = IterationAdaptiveDT
        dt = 20
        optimal_iterations = 15
        iteration_window = 1
        growth_factor = 1.2
        cutback_factor = 0.8
        timestep_limiting_postprocessor = timestep_max_pp
    []
[]

[Outputs]
    [exodus]
        type = Exodus
        sync_only = false
        # output at key moment in the first two cycles, and then at the end of the simulation
        sync_times = '110.0 480.0 590.0 1600.0 1710.0 2080.0 2190.0 3400.0 8.0e4'
    []
    csv = true
    hide = 'dt
            Int_C_mobile_W Int_C_trapped_W Int_C_total_W
            Int_C_mobile_Cu Int_C_trapped_Cu Int_C_total_Cu
            Int_C_mobile_CuCrZr Int_C_trapped_CuCrZr Int_C_total_CuCrZr'
    perf_graph = true
[]

(test/tests/divertor_monoblock/tests)

[Tests]
  issues = '#144'
  design = 'divertor_monoblock/index.md'
  [Shimada2024_input_check]
    type = RunApp
    input = 'divertor_monoblock.i'
    requirement = 'The system shall maintain a working input file to model heat and tritium transport in a divertor monoblock during pulsed operation.'
    check_input = True
    method = opt
  []
  [Shimada2024_run]
    type = Exodiff
    input = 'divertor_monoblock.i'
    exodiff = 'divertor_monoblock_exodus.e'
    cli_args = 'Outputs/exodus/sync_only=true Executioner/end_time=1600 Mesh/ccmg/rings="1 10 8 40" Executioner/nl_rel_tol=1e-8 Executioner/nl_abs_tol=1e-11'
    requirement = 'The system shall model heat and tritium transport in a divertor monoblock during pulsed operation.'
    rel_err = 8e-4 # increasing the relative error tolerance because the test is sensitive to the environment. This error is still small for physical problems.
    heavy = true
    min_parallel = 2
    max_time = 1000
  []
[]

(test/tests/divertor_monoblock/divertor_monoblock.i)

### Nomenclatures
### C_mobile_j      mobile H concentration in "j" material, where j = CuCrZr, Cu, W
### C_trapped_j     trapped H concentration in "j" material, where j = CuCrZr, Cu, W
### C_total_j       total H concentration in "j" material, where j = CuCrZr, Cu, W
###
### S_empty_j       empty site concentration in "j" material, where j = CuCrZr, Cu, W
### S_trapped_j     trapped site concentration in "j" material, where j = CuCrZr, Cu, W
### S_total_j       total site H concentration in "j" material, where j = CuCrZr, Cu, W
###
### F_permeation    permeation flux
### F_recombination recombination flux
###
### Sc_             Scaled
### Int_            Integrated
### ScInt_          Scaled and integrated

tungsten_atomic_density = ${units 6.338e28 m^-3}

## Mesh generated by ConcentricCircleMeshGenerator and saved as "2DMonoblock.e"
[Mesh]
    [ccmg]
        type = ConcentricCircleMeshGenerator
        num_sectors = 36
        rings = '1 30 20 110'
        radii = '${units 6 mm -> m} ${units 7.5 mm -> m} ${units 8.5 mm -> m}'
        has_outer_square = on
        pitch = ${units 28 mm -> m}
        portion = left_half
        preserve_volumes = false
        smoothing_max_it = 3
    []
    [ssbsg1]
        type = SideSetsBetweenSubdomainsGenerator
        input = ccmg
        primary_block = '4'     # W
        paired_block = '3'      # Cu
        new_boundary = '4to3'
    []
    [ssbsg2]
        type = SideSetsBetweenSubdomainsGenerator
        input = ssbsg1
        primary_block = '3'     # Cu
        paired_block = '4'      # W
        new_boundary = '3to4'
    []
    [ssbsg3]
        type = SideSetsBetweenSubdomainsGenerator
        input = ssbsg2
        primary_block = '3'     # Cu
        paired_block = '2'      # CuCrZr
        new_boundary = '3to2'
    []
    [ssbsg4]
        type = SideSetsBetweenSubdomainsGenerator
        input = ssbsg3
        primary_block = '2'     # CuCrZr
        paired_block = '3'      # Cu
        new_boundary = '2to3'
    []
    [ssbsg5]
        type = SideSetsBetweenSubdomainsGenerator
        input = ssbsg4
        primary_block = '2'     # CuCrZr
        paired_block = '1'      # H2O
        new_boundary = '2to1'
    []
    [bdg]
        type = BlockDeletionGenerator
        input = ssbsg5
        block = '1'             # H2O
    []
[]

[Problem]
    type = ReferenceResidualProblem
    extra_tag_vectors = 'ref'
    reference_vector = 'ref'
[]

[Variables]
    [temperature]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 300 K}
    []
    ######################### Variables for W (block = 4)
    [C_mobile_W]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-20 m^-3}
        block = 4
    []
    [C_trapped_W]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-15 m^-3}
        block = 4
    []
    ######################### Variables for Cu (block = 3)
    [C_mobile_Cu]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 5.0e-17 m^-3}
        block = 3
    []
    [C_trapped_Cu]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-15 m^-3}
        block = 3
    []
    ######################### Variables for CuCrZr (block = 2)
    [C_mobile_CuCrZr]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-15 m^-3}
        block = 2
    []
    [C_trapped_CuCrZr]
        order = FIRST
        family = LAGRANGE
        initial_condition = ${units 1.0e-15 m^-3}
        block = 2
    []
[]

[AuxVariables]
    [flux_y]
        order = FIRST
        family = MONOMIAL
    []
    ############################## AuxVariables for W (block = 4)
    [Sc_C_mobile_W]
        block = 4
    []
    [Sc_C_trapped_W]
        block = 4
    []
    [C_total_W]
        block = 4
    []
    [Sc_C_total_W]
        block = 4
    []
    [S_empty_W]
        block = 4
    []
    [Sc_S_empty_W]
        block = 4
    []
    [S_trapped_W]
        block = 4
    []
    [Sc_S_trapped_W]
        block = 4
    []
    [S_total_W]
        block = 4
    []
    [Sc_S_total_W]
        block = 4
    []
    ############################## AuxVariables for Cu (block = 3)
    [Sc_C_mobile_Cu]
        block = 3
    []
    [Sc_C_trapped_Cu]
        block = 3
    []
    [C_total_Cu]
        block = 3
    []
    [Sc_C_total_Cu]
        block = 3
    []
    [S_empty_Cu]
        block = 3
    []
    [Sc_S_empty_Cu]
        block = 3
    []
    [S_trapped_Cu]
        block = 3
    []
    [Sc_S_trapped_Cu]
        block = 3
    []
    [S_total_Cu]
        block = 3
    []
    [Sc_S_total_Cu]
        block = 3
    []
    ############################## AuxVariables for CuCrZr (block = 2)
    [Sc_C_mobile_CuCrZr]
        block = 2
    []
    [Sc_C_trapped_CuCrZr]
        block = 2
    []
    [C_total_CuCrZr]
        block = 2
    []
    [Sc_C_total_CuCrZr]
        block = 2
    []
    [S_empty_CuCrZr]
        block = 2
    []
    [Sc_S_empty_CuCrZr]
        block = 2
    []
    [S_trapped_CuCrZr]
        block = 2
    []
    [Sc_S_trapped_CuCrZr]
        block = 2
    []
    [S_total_CuCrZr]
        block = 2
    []
    [Sc_S_total_CuCrZr]
        block = 2
    []
[]

[Kernels]
    ############################## Kernels for W (block = 4)
    [diff_W]
        type = ADMatDiffusion
        variable = C_mobile_W
        diffusivity = diffusivity_W
        block = 4
        extra_vector_tags = ref
    []
    [time_diff_W]
        type = ADTimeDerivative
        variable = C_mobile_W
        block = 4
        extra_vector_tags = ref
    []
    [coupled_time_W]
        type = ScaledCoupledTimeDerivative
        variable = C_mobile_W
        v = C_trapped_W
        factor = 1e0
        block = 4
        extra_vector_tags = ref
    []
    [heat_conduction_W]
        type = HeatConduction
        variable = temperature
        diffusion_coefficient = thermal_conductivity_W
        block = 4
        extra_vector_tags = ref
    []
    [time_heat_conduction_W]
        type = SpecificHeatConductionTimeDerivative
        variable = temperature
        specific_heat = specific_heat_W
        density = density_W
        block = 4
        extra_vector_tags = ref
    []
    ############################## Kernels for Cu (block = 3)
    [diff_Cu]
        type = ADMatDiffusion
        variable = C_mobile_Cu
        diffusivity = diffusivity_Cu
        block = 3
        extra_vector_tags = ref
    []
    [time_diff_Cu]
        type = ADTimeDerivative
        variable = C_mobile_Cu
        block = 3
        extra_vector_tags = ref
    []
    [coupled_time_Cu]
        type = ScaledCoupledTimeDerivative
        variable = C_mobile_Cu
        v = C_trapped_Cu
        factor = 1e0
        block = 3
        extra_vector_tags = ref
    []
    [heat_conduction_Cu]
        type = HeatConduction
        variable = temperature
        diffusion_coefficient = thermal_conductivity_Cu
        block = 3
        extra_vector_tags = ref
    []
    [time_heat_conduction_Cu]
        type = SpecificHeatConductionTimeDerivative
        variable = temperature
        specific_heat = specific_heat_Cu
        density = density_Cu
        block = 3
        extra_vector_tags = ref
    []
    ############################## Kernels for CuCrZr (block = 2)
    [diff_CuCrZr]
        type = ADMatDiffusion
        variable = C_mobile_CuCrZr
        diffusivity = diffusivity_CuCrZr
        block = 2
        extra_vector_tags = ref
    []
    [time_diff_CuCrZr]
        type = ADTimeDerivative
        variable = C_mobile_CuCrZr
        block = 2
        extra_vector_tags = ref
    []
    [coupled_time_CuCrZr]
        type = ScaledCoupledTimeDerivative
        variable = C_mobile_CuCrZr
        v = C_trapped_CuCrZr
        factor = 1e0
        block = 2
        extra_vector_tags = ref
    []
    [heat_conduction_CuCrZr]
        type = HeatConduction
        variable = temperature
        diffusion_coefficient = thermal_conductivity_CuCrZr
        block = 2
        extra_vector_tags = ref
    []
    [time_heat_conduction_CuCrZr]
        type = SpecificHeatConductionTimeDerivative
        variable = temperature
        specific_heat = specific_heat_CuCrZr
        density = density_CuCrZr
        block = 2
        extra_vector_tags = ref
    []
[]

[AuxKernels]
    ############################## AuxKernels for W (block = 4)
    [Scaled_mobile_W]
        variable = Sc_C_mobile_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_mobile_W
    []
    [Scaled_trapped_W]
        variable = Sc_C_trapped_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_trapped_W
    []
    [total_W]
        variable = C_total_W
        type = ParsedAux
        expression = 'C_mobile_W + C_trapped_W'
        coupled_variables = 'C_mobile_W C_trapped_W'
    []
    [Scaled_total_W]
        variable = Sc_C_total_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_total_W
    []
    [empty_sites_W]
        variable = S_empty_W
        type = EmptySitesAux
        N = ${units 1.0e0 m^-3}       # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        # Ct0 = ${units 1.0e-4 m^-3}   # E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        Ct0 = ${units 1.0e-4 m^-3}    # E.A. Hodille et al 2021 Nucl. Fusion 61 1260033, trap 2
        trap_per_free = 1.0e0         # 1.0e1
        trapped_concentration_variables = C_trapped_W
    []
    [scaled_empty_W]
        variable = Sc_S_empty_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_empty_W
    []
    [trapped_sites_W]
        variable = S_trapped_W
        type = NormalizationAux
        normal_factor = 1e0
        source_variable = C_trapped_W
    []
    [scaled_trapped_W]
        variable = Sc_S_trapped_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_trapped_W
    []
    [total_sites_W]
        variable = S_total_W
        type = ParsedAux
        expression = 'S_trapped_W + S_empty_W'
        coupled_variables = 'S_trapped_W S_empty_W'
    []
    [scaled_total_W]
        variable = Sc_S_total_W
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_total_W
    []
    ############################## AuxKernels for Cu (block = 3)
    [Scaled_mobile_Cu]
        variable = Sc_C_mobile_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_mobile_Cu
    []
    [Scaled_trapped_Cu]
        variable = Sc_C_trapped_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_trapped_Cu
    []
    [total_Cu]
        variable = C_total_Cu
        type = ParsedAux
        expression = 'C_mobile_Cu + C_trapped_Cu'
        coupled_variables = 'C_mobile_Cu C_trapped_Cu'
    []
    [Scaled_total_Cu]
        variable = Sc_C_total_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_total_Cu
    []
    [empty_sites_Cu]
        variable = S_empty_Cu
        type = EmptySitesAux
        N = ${units 1.0e0 m^-3}     # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0 = ${units 5.0e-5 m^-3}  # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        trap_per_free = 1.0e0       # 1.0e1
        trapped_concentration_variables = C_trapped_Cu
    []
    [scaled_empty_Cu]
        variable = Sc_S_empty_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_empty_Cu
    []
    [trapped_sites_Cu]
        variable = S_trapped_Cu
        type = NormalizationAux
        normal_factor = 1e0
        source_variable = C_trapped_Cu
    []
    [scaled_trapped_Cu]
        variable = Sc_S_trapped_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_trapped_Cu
    []
    [total_sites_Cu]
        variable = S_total_Cu
        type = ParsedAux
        expression = 'S_trapped_Cu + S_empty_Cu'
        coupled_variables = 'S_trapped_Cu S_empty_Cu'
    []
    [scaled_total_Cu]
        variable = Sc_S_total_Cu
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_total_Cu
    []
    ############################## AuxKernels for CuCrZr (block = 2)
    [Scaled_mobile_CuCrZr]
        variable = Sc_C_mobile_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_mobile_CuCrZr
    []
    [Scaled_trapped_CuCrZr]
        variable = Sc_C_trapped_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_trapped_CuCrZr
    []
    [total_CuCrZr]
        variable = C_total_CuCrZr
        type = ParsedAux
        expression = 'C_mobile_CuCrZr + C_trapped_CuCrZr'
        coupled_variables = 'C_mobile_CuCrZr C_trapped_CuCrZr'
    []
    [Scaled_total_CuCrZr]
        variable = Sc_C_total_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = C_total_CuCrZr
    []
    [empty_sites_CuCrZr]
        variable = S_empty_CuCrZr
        type = EmptySitesAux
        N = ${units 1.0e0 m^-3}     # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0 = ${units 5.0e-5 m^-3}  # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # Ct0 = ${units 4.0e-2 m^-3} # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        trap_per_free = 1.0e0       # 1.0e1
        trapped_concentration_variables = C_trapped_CuCrZr
    []
    [scaled_empty_CuCrZr]
        variable = Sc_S_empty_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_empty_CuCrZr
    []
    [trapped_sites_CuCrZr]
        variable = S_trapped_CuCrZr
        type = NormalizationAux
        normal_factor = 1e0
        source_variable = C_trapped_CuCrZr
    []
    [scaled_trapped_CuCrZr]
        variable = Sc_S_trapped_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_trapped_CuCrZr
    []
    [total_sites_CuCrZr]
        variable = S_total_CuCrZr
        type = ParsedAux
        expression = 'S_trapped_CuCrZr + S_empty_CuCrZr'
        coupled_variables = 'S_trapped_CuCrZr S_empty_CuCrZr'
    []
    [scaled_total_CuCrZr]
        variable = Sc_S_total_CuCrZr
        type = NormalizationAux
        normal_factor = ${tungsten_atomic_density}
        source_variable = S_total_CuCrZr
    []
    [flux_y_W]
        type = DiffusionFluxAux
        diffusivity = diffusivity_W
        variable = flux_y
        diffusion_variable = C_mobile_W
        component = y
        block = 4
    []
    [flux_y_Cu]
        type = DiffusionFluxAux
        diffusivity = diffusivity_Cu
        variable = flux_y
        diffusion_variable = C_mobile_Cu
        component = y
        block = 3
    []
    [flux_y_CuCrZr]
        type = DiffusionFluxAux
        diffusivity = diffusivity_CuCrZr
        variable = flux_y
        diffusion_variable = C_mobile_CuCrZr
        component = y
        block = 2
    []
[]

[InterfaceKernels]
    [tied_4to3]
        type = ADPenaltyInterfaceDiffusion
        variable = C_mobile_W
        neighbor_var = C_mobile_Cu
        penalty = 0.05            #  it will not converge with > 0.1, but it creates negative C_mobile _Cu with << 0.1
        # jump_prop_name = solubility_ratio_4to3
        jump_prop_name = solubility_ratio
        boundary = '4to3'
    []
    [tied_3to2]
        type = ADPenaltyInterfaceDiffusion
        variable = C_mobile_Cu
        neighbor_var = C_mobile_CuCrZr
        penalty = 0.05            #  it will not converge with > 0.1, but it creates negative C_mobile _Cu with << 0.1
        # jump_prop_name = solubility_ratio_3to2
        jump_prop_name = solubility_ratio
        boundary = '3to2'
    []
[]

[NodalKernels]
    ############################## NodalKernels for W (block = 4)
    [time_W]
        type = TimeDerivativeNodalKernel
        variable = C_trapped_W
    []
    [trapping_W]
        type = TrappingNodalKernel
        variable = C_trapped_W
        temperature = temperature
        alpha_t = 2.75e11      # 1e15
        N = 1.0e0  # = (1e0) x (${tungsten_atomic_density} #/m^3)
        # Ct0 = 1.0e-4                # E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        Ct0 = 1.0e-4                # E.A. Hodille et al 2021 Nucl. Fusion 61 1260033, trap 2
        trap_per_free = 1.0e0       # 1.0e1
        mobile_concentration = 'C_mobile_W'
        extra_vector_tags = ref
    []
    [release_W]
        type = ReleasingNodalKernel
        alpha_r = 8.4e12    # 1.0e13
        temperature = temperature
        # detrapping_energy = 9863.9    # = 0.85 eV    E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 1
        detrapping_energy = 11604.6   # = 1.00 eV    E.A. Hodille et al 2021 Nucl. Fusion 61 126003, trap 2
        variable = C_trapped_W
    []
    ############################## NodalKernels for Cu (block = 3)
    [time_Cu]
        type = TimeDerivativeNodalKernel
        variable = C_trapped_Cu
    []
    [trapping_Cu]
        type = TrappingNodalKernel
        variable = C_trapped_Cu
        temperature = temperature
        alpha_t = 2.75e11      # 1e15
        N = 1.0e0  # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0 = 5.0e-5                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        trap_per_free = 1.0e0       # 1.0e1
        mobile_concentration = 'C_mobile_Cu'
        extra_vector_tags = ref
    []
    [release_Cu]
        type = ReleasingNodalKernel
        alpha_r = 8.4e12    # 1.0e13
        temperature = temperature
        detrapping_energy = 5802.3    # = 0.50eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 3
        variable = C_trapped_Cu
    []
    ############################## NodalKernels for CuCrZr (block = 2)
    [time_CuCrZr]
        type = TimeDerivativeNodalKernel
        variable = C_trapped_CuCrZr
    []
    [trapping_CuCrZr]
        type = TrappingNodalKernel
        variable = C_trapped_CuCrZr
        temperature = temperature
        alpha_t = 2.75e11      # 1e15
        N = 1.0e0  # = ${tungsten_atomic_density} #/m^3 (W lattice density)
        Ct0 = 5.0e-5                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # Ct0 = 4.0e-2                # R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        trap_per_free = 1.0e0       # 1.0e1
        mobile_concentration = 'C_mobile_CuCrZr'
        extra_vector_tags = ref
    []
    [release_CuCrZr]
        type = ReleasingNodalKernel
        alpha_r = 8.4e12    # 1.0e13
        temperature = temperature
        detrapping_energy = 5802.3    # = 0.50eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 4
        # detrapping_energy = 9631.8   # = 0.83 eV  R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038, trap 5
        variable = C_trapped_CuCrZr
    []
[]

[BCs]
    [C_mob_W_top_flux]
        type = FunctionNeumannBC
        variable = C_mobile_W
        boundary = 'top'
        function = mobile_flux_bc_func
    []
    [mobile_tube]
        type = DirichletBC
        variable = C_mobile_CuCrZr
        value = 1.0e-18
        boundary = '2to1'
    []
    [temp_top]
        type = FunctionNeumannBC
        variable = temperature
        boundary = 'top'
        function = temp_flux_bc_func
    []
    [temp_tube]
        type = FunctionDirichletBC
        variable = temperature
        boundary = '2to1'
        function = temp_inner_func
    []
[]

[Functions]
    ### Maximum mobile flux of 7.90e-13 at the top surface (1.0e-4 [m])
    ### 1.80e23/m^2-s = (5.0e23/m^2-s) *(1-0.999) = (7.90e-13)*(${tungsten_atomic_density})/(1.0e-4)  at steady state
    [mobile_flux_bc_func]
        type = ParsedFunction
        expression =   'if((t % 1600) < 100.0, 0.0      + 7.90e-13*(t % 1600)/100,
                        if((t % 1600) < 500.0, 7.90e-13,
                        if((t % 1600) < 600.0, 7.90e-13 - 7.90e-13*((t % 1600)-500)/100, 0.0)))'
    []
    ### Heat flux of 10MW/m^2 at steady state
    [temp_flux_bc_func]
        type = ParsedFunction
        expression =   'if((t % 1600) < 100.0, 0.0   + 1.0e7*(t % 1600)/100,
                        if((t % 1600) < 500.0, 1.0e7,
                        if((t % 1600) < 600.0, 1.0e7 - 1.0e7*((t % 1600)-500)/100, 300)))'
    []
    ### Maximum coolant temperature of 552K at steady state
    [temp_inner_func]
        type = ParsedFunction
        expression =   'if((t % 1600) < 100.0, 300.0 + (552-300)*(t % 1600)/100,
                        if((t % 1600) < 500.0, 552,
                        if((t % 1600) < 600.0, 552.0 - (552-300)*((t % 1600)-500)/100, 300)))'
    []
    [timestep_function]
        type = ParsedFunction
        expression = 'if((t % 1600) <   10.0,  20,
                      if((t % 1600) <   90.0,  40,
                      if((t % 1600) <  110.0,  20,
                      if((t % 1600) <  480.0,  40,
                      if((t % 1600) <  500.0,  20,
                      if((t % 1600) <  590.0,   4,
                      if((t % 1600) <  610.0,  20,
                      if((t % 1600) < 1500.0, 200,
                      if((t % 1600) < 1600.0,  40,  2)))))))))'
    []
[]

[Materials]
    ############################## Materials for W (block = 4)
    [diffusivity_W]
        type = ADParsedMaterial
        property_name = diffusivity_W
        coupled_variables = 'temperature'
        block = 4
        expression = '2.4e-7*exp(-4525.8/temperature)'    # H diffusivity in W
        outputs = all
    []
    [solubility_W]
        type = ADParsedMaterial
        property_name = solubility_W
        coupled_variables = 'temperature'
        block = 4
        # expression = '2.95e-5 *exp(-12069.0/temperature)'              # H solubility in W = (1.87e24)/(${tungsten_atomic_density}) [#/m^3]
        expression = '2.95e-5 *exp(-12069.0/temperature) + 4.95e-8 * exp(-6614.6/temperature)'    # H solubility in W = (1.87e24)/(${tungsten_atomic_density}) [#/m^3]
        outputs = all
    []
    [converter_to_regular_W]
        type = MaterialADConverter
        ad_props_in = 'diffusivity_W'
        reg_props_out = 'diffusivity_W_nonAD'
        block = 4
    []
    [heat_transfer_W]
        type = GenericConstantMaterial
        prop_names = 'density_W'
        prop_values = '19300'                # [g/m^3]
        block = 4
    []
    [specific_heat_W]
        type = ParsedMaterial
        property_name = specific_heat_W
        coupled_variables = 'temperature'
        block = 4
        expression = '1.16e2 + 7.11e-2 * temperature - 6.58e-5 * temperature^2 + 3.24e-8 * temperature^3 -5.45e-12 * temperature^4'    # ~ 132[J/kg-K]
        outputs = all
    []
    [thermal_conductivity_W]
        type = ParsedMaterial
        property_name = thermal_conductivity_W
        coupled_variables = 'temperature'
        block = 4
        # expression = '-7.8e-9 * temperature^3 + 5.0e-5 * temperature^2 - 1.1e-1 * temperature + 1.8e2'    # ~ 173.0 [ W/m-K]   from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression = '2.41e2 - 2.90e-1 * temperature + 2.54e-4 * temperature^2 - 1.03e-7 * temperature^3 + 1.52e-11 * temperature^4'    # ~ 173.0 [ W/m-K]
        outputs = all
    []
    ############################## Materials for Cu (block = 3)
    [diffusivity_Cu]
        type = ADParsedMaterial
        property_name = diffusivity_Cu
        coupled_variables = 'temperature'
        block = 3
        expression = '6.60e-7*exp(-4525.8/temperature)'    # H diffusivity in Cu
        outputs = all
    []
    [solubility_Cu]
        type = ADParsedMaterial
        property_name = solubility_Cu
        coupled_variables = 'temperature'
        block = 3
        expression = '4.95e-5*exp(-6614.6/temperature)'    # H solubility in Cu = (3.14e24)/(${tungsten_atomic_density}) [#/m^3]
        outputs = all
    []
    [converter_to_regular_Cu]
        type = MaterialADConverter
        ad_props_in = 'diffusivity_Cu'
        reg_props_out = 'diffusivity_Cu_nonAD'
        block = 3
    []
    [heat_transfer_Cu]
        type = GenericConstantMaterial
        prop_names = 'density_Cu'
        prop_values = '8960.0'                # [g/m^3]
        block = 3
    []
    [specific_heat_Cu]
        type = ParsedMaterial
        property_name = specific_heat_Cu
        coupled_variables = 'temperature'
        block = 3
        expression = '3.16e2 + 3.18e-1 * temperature - 3.49e-4 * temperature^2 + 1.66e-7 * temperature^3'    # ~ 384 [J/kg-K]
        outputs = all
    []
    [thermal_conductivity_Cu]
        type = ParsedMaterial
        property_name = thermal_conductivity_Cu
        coupled_variables = 'temperature'
        block = 3
        # expression = '-3.9e-8 * temperature^3 + 3.8e-5 * temperature^2 - 7.9e-2 * temperature + 4.0e2'    # ~ 401.0  [ W/m-K] from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression = '4.21e2 - 6.85e-2 * temperature'    # ~ 400.0 [ W/m-K]
        outputs = all
    []
    ############################## Materials for CuCrZr (block = 2)
    [diffusivity_CuCrZr]
        type = ADParsedMaterial
        property_name = diffusivity_CuCrZr
        coupled_variables = 'temperature'
        block = 2
        expression = '3.90e-7*exp(-4873.9/temperature)'    # H diffusivity in CuCrZr
        outputs = all
    []
    [solubility_CuCrZr]
        type = ADParsedMaterial
        property_name = solubility_CuCrZr
        coupled_variables = 'temperature'
        block = 2
        expression = '6.75e-6*exp(-4525.8/temperature)'    # H solubility in CuCrZr = (4.28e23)/(${tungsten_atomic_density}) [#/m^3]
        outputs = all
    []
    [converter_to_regular_CuCrZr]
        type = MaterialADConverter
        ad_props_in = 'diffusivity_CuCrZr'
        reg_props_out = 'diffusivity_CuCrZr_nonAD'
        block = 2
    []
    [heat_transfer_CuCrZr]
        type = GenericConstantMaterial
        prop_names = 'density_CuCrZr specific_heat_CuCrZr'
        prop_values = '8900.0 390.0'            # [g/m^3], [ W/m-K], [J/kg-K]
        block = 2
    []
    [thermal_conductivity_CuCrZr]
        type = ParsedMaterial
        property_name = thermal_conductivity_CuCrZr
        coupled_variables = 'temperature'
        block = 2
        # expression = '5.3e-7 * temperature^3 - 6.5e-4 * temperature^2 + 2.6e-1 * temperature + 3.1e2'    # ~ 320.0  [ W/m-K] from R. Delaporte-Mathurin et al 2021 Nucl. Fusion 61 036038,
        expression = '3.87e2 - 1.28e-1 * temperature'    # ~ 349 [ W/m-K]
        outputs = all
    []
    ############################## Materials for others
    [interface_jump_4to3]
        type = SolubilityRatioMaterial
        solubility_primary = solubility_W
        solubility_secondary = solubility_Cu
        boundary = '4to3'
        concentration_primary = C_mobile_W
        concentration_secondary = C_mobile_Cu
    []
    [interface_jump_3to2]
        type = SolubilityRatioMaterial
        solubility_primary = solubility_Cu
        solubility_secondary = solubility_CuCrZr
        boundary = '3to2'
        concentration_primary = C_mobile_Cu
        concentration_secondary = C_mobile_CuCrZr
    []
[]

[Postprocessors]
    ############################################################ Postprocessors for W (block = 4)
    [F_recombination]
        type = SideDiffusiveFluxAverage
        boundary = 'top'
        diffusivity = 5.01e-24   # (3.01604928)/(6.02e23)/[gram(T)/m^2]
        # diffusivity = 5.508e-19   # (1.0e3)*(1.0e3)/(6.02e23)/(3.01604928) [gram(T)/m^2]
        variable = Sc_C_total_W
    []
    [F_permeation]
        type = SideDiffusiveFluxAverage
        boundary = '2to1'
        diffusivity = 5.01e-24   # (3.01604928)/(6.02e23)/[gram(T)/m^2]
        # diffusivity = 5.508e-19   # (1.0e3)*(1.0e3)/(6.02e23)/(3.01604928) [gram(T)/m^2]
        variable = Sc_C_total_CuCrZr
    []

    [Int_C_mobile_W]
        type = ElementIntegralVariablePostprocessor
        variable = C_mobile_W
        block = 4
    []
    [ScInt_C_mobile_W]
        type = ScalePostprocessor
        value =  Int_C_mobile_W
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_W]
        type = ElementIntegralVariablePostprocessor
        variable = C_trapped_W
        block = 4
    []
    [ScInt_C_trapped_W]
        type = ScalePostprocessor
        value = Int_C_trapped_W
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_W]
        type = ElementIntegralVariablePostprocessor
        variable = C_total_W
        block = 4
    []
    [ScInt_C_total_W]
        type = ScalePostprocessor
        value = Int_C_total_W
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    # ############################################################ Postprocessors for Cu (block = 3)
    [Int_C_mobile_Cu]
        type = ElementIntegralVariablePostprocessor
        variable = C_mobile_Cu
        block = 3
    []
    [ScInt_C_mobile_Cu]
        type = ScalePostprocessor
        value =  Int_C_mobile_Cu
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_Cu]
        type = ElementIntegralVariablePostprocessor
        variable = C_trapped_Cu
        block = 3
    []
    [ScInt_C_trapped_Cu]
        type = ScalePostprocessor
        value = Int_C_trapped_Cu
        scaling_factor = 3.44e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_Cu]
        type = ElementIntegralVariablePostprocessor
        variable = C_total_Cu
        block = 3
    []
    [ScInt_C_total_Cu]
        type = ScalePostprocessor
        value = Int_C_total_Cu
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    # ############################################################ Postprocessors for CuCrZr (block = 2)
    [Int_C_mobile_CuCrZr]
        type = ElementIntegralVariablePostprocessor
        variable = C_mobile_CuCrZr
        block = 2
    []
    [ScInt_C_mobile_CuCrZr]
        type = ScalePostprocessor
        value =  Int_C_mobile_CuCrZr
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_trapped_CuCrZr]
        type = ElementIntegralVariablePostprocessor
        variable = C_trapped_CuCrZr
        block = 2
    []
    [ScInt_C_trapped_CuCrZr]
        type = ScalePostprocessor
        value = Int_C_trapped_CuCrZr
        scaling_factor = 3.44e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    [Int_C_total_CuCrZr]
        type = ElementIntegralVariablePostprocessor
        variable = C_total_CuCrZr
        block = 2
    []
    [ScInt_C_total_CuCrZr]
        type = ScalePostprocessor
        value = Int_C_total_CuCrZr
        scaling_factor = 3.491e10   # (1.0e3)*(1.0e3)*(${tungsten_atomic_density})/(6.02e23)/(3.01604928) [gram(T)/m^2]
    []
    ############################################################ Postprocessors for others
    [dt]
        type = TimestepSize
    []
    [temperature_top]
        type = PointValue
        variable = temperature
        point = '0 14.0e-3 0'
    []
    [temperature_tube]
        type = PointValue
        variable = temperature
        point = '0 6.0e-3 0'
    []
    # limit timestep
    [timestep_max_pp] # s
    type = FunctionValuePostprocessor
    function = timestep_function
  []
[]

[VectorPostprocessors]
    [line]
        type = LineValueSampler
        start_point = '0 14.0e-3 0'
        end_point = '0 6.0e-3 0'
        num_points = 100
        sort_by = 'y'
        variable = 'C_total_W C_total_Cu C_total_CuCrZr C_mobile_W C_mobile_Cu C_mobile_CuCrZr C_trapped_W C_trapped_Cu C_trapped_CuCrZr flux_y temperature'
        execute_on = timestep_end
    []
[]


[Preconditioning]
    [smp]
        type = SMP
        full = true
    []
[]

[Executioner]
    type = Transient
    scheme = bdf2
    solve_type = NEWTON
    petsc_options_iname = '-pc_type'
    petsc_options_value = 'lu'
    nl_rel_tol  = 1e-6 # 1e-8 works for 1 cycle
    nl_abs_tol  = 1e-7 # 1e-11 works for 1 cycle
    end_time = 8.0e4   # 50 ITER shots (3.0e4 s plasma, 2.0e4 SSP)
    automatic_scaling = true
    line_search = 'none'
    dtmin = 1e-4
    nl_max_its = 18
    [TimeStepper]
        type = IterationAdaptiveDT
        dt = 20
        optimal_iterations = 15
        iteration_window = 1
        growth_factor = 1.2
        cutback_factor = 0.8
        timestep_limiting_postprocessor = timestep_max_pp
    []
[]

[Outputs]
    [exodus]
        type = Exodus
        sync_only = false
        # output at key moment in the first two cycles, and then at the end of the simulation
        sync_times = '110.0 480.0 590.0 1600.0 1710.0 2080.0 2190.0 3400.0 8.0e4'
    []
    csv = true
    hide = 'dt
            Int_C_mobile_W Int_C_trapped_W Int_C_total_W
            Int_C_mobile_Cu Int_C_trapped_Cu Int_C_total_Cu
            Int_C_mobile_CuCrZr Int_C_trapped_CuCrZr Int_C_total_CuCrZr'
    perf_graph = true
[]

General description of the simulation case and corresponding input file
Results
Complete input file
References

Symbol	Variable or Physical Property	Unit
$var element = document.getElementById("moose-equation-ee4ffae1-c66e-4f23-b172-43ef565dd9b6");katex.render("C_s", element, {displayMode:false,throwOnError:false});$	Concentration of solute (mobile) species	m $var element = document.getElementById("moose-equation-ded701c9-b552-41cf-8287-6960296e22e1");katex.render("^{\u22123}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-c9e9a657-58d8-4e2e-b232-d308faa4b365");katex.render("C_t", element, {displayMode:false,throwOnError:false});$	Concentration of trapped species	m $var element = document.getElementById("moose-equation-9b1b8ecf-441d-4dc5-8491-bb045d1441f7");katex.render("^{\u22123}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-13097708-8fd3-4fd7-ba21-df4ab09dee09");katex.render("T", element, {displayMode:false,throwOnError:false});$	Temperature	K
$var element = document.getElementById("moose-equation-42649e70-f534-4486-80f1-6cc5829a4f21");katex.render("D", element, {displayMode:false,throwOnError:false});$	Diffusivity of solute (mobile) species	m $var element = document.getElementById("moose-equation-0878e0df-74dc-4c19-a541-7eb1e98130f0");katex.render("^{2}", element, {displayMode:false,throwOnError:false});$ s $var element = document.getElementById("moose-equation-5598b8ae-7af3-4fd9-bce4-6368d306b799");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-ec4269ca-d7d0-49de-80ba-4da73f4312c7");katex.render("K_s", element, {displayMode:false,throwOnError:false});$	Solubility of solute (mobile) species	m $var element = document.getElementById("moose-equation-2a9fbe36-6211-4ef7-a8f1-40948b8d71cf");katex.render("^{\u22123}", element, {displayMode:false,throwOnError:false});$ Pa $var element = document.getElementById("moose-equation-b322d565-3046-4952-abf5-19f39907a9e8");katex.render("^{-1/2}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-cb1f2007-a77b-4b7f-ad22-49f17f27395b");katex.render("C_{total}", element, {displayMode:false,throwOnError:false});$	Total concentration of species $var element = document.getElementById("moose-equation-1521c4e9-c20b-47ae-8aa3-000f28e95e54");katex.render("C_{total}=C_s + C_t", element, {displayMode:false,throwOnError:false});$	m $var element = document.getElementById("moose-equation-56e4a7ec-cc8b-449f-8657-f8e2c440aed0");katex.render("^{\u22123}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-ef9923e4-7f05-4439-ae46-6649f49636fd");katex.render("C_t^e", element, {displayMode:false,throwOnError:false});$	Concentration of empty trapping sites	m $var element = document.getElementById("moose-equation-5ea118af-17fe-44e5-83f1-04ef17dbf933");katex.render("^{\u22123}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-e430e47c-bcfd-451c-a986-24edb57a3ccb");katex.render("C_t^0", element, {displayMode:false,throwOnError:false});$	Concentration of trapping sites $var element = document.getElementById("moose-equation-af91a6c1-d944-4c19-8274-870eb5a6eb72");katex.render("C_t^e=C_t^0-C_t", element, {displayMode:false,throwOnError:false});$	m $var element = document.getElementById("moose-equation-d74927f3-65a0-42d0-ad78-c66e2c0442b2");katex.render("^{\u22123}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-b81585bb-1cc6-426d-9e57-b85a116f09d6");katex.render("\\alpha_t", element, {displayMode:false,throwOnError:false});$	Trapping rate coefficient, $var element = document.getElementById("moose-equation-c35a0cd9-41db-49d4-8bfa-35403f119fa0");katex.render("2.75 \\times 10^{11}", element, {displayMode:false,throwOnError:false});$	s $var element = document.getElementById("moose-equation-8a6792b3-2a1b-4cd4-98c3-851d8a27de95");katex.render("^{\u22121}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-60b93253-b3f3-4a2c-b8a5-84e177f2a928");katex.render("\\alpha_r", element, {displayMode:false,throwOnError:false});$	Release rate coefficient $var element = document.getElementById("moose-equation-703db0db-2f24-4291-b732-c510f8c080e0");katex.render("\\alpha_r=\\alpha_{r0} \\exp\u2061((-E_{dt})/(k_b T))", element, {displayMode:false,throwOnError:false});$	s $var element = document.getElementById("moose-equation-8e9049a7-aab3-4a4e-8e49-4d1c9d3f8ed9");katex.render("^{\u22121}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-4cf9492a-da82-49b0-b970-955069edf21f");katex.render("\\alpha_{r0}", element, {displayMode:false,throwOnError:false});$	Pre-exponential factor, $var element = document.getElementById("moose-equation-b24fd735-db35-4726-822c-bf02a3c8624e");katex.render("8.40 \\times 10^{12}", element, {displayMode:false,throwOnError:false});$	s $var element = document.getElementById("moose-equation-5563a7ca-3c44-4483-b6c3-d006c0baa935");katex.render("^{\u22121}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-f6d4931a-cf2a-413e-ba9e-61910200a3e6");katex.render("E_{dt}", element, {displayMode:false,throwOnError:false});$	Detrapping energy	eV
$var element = document.getElementById("moose-equation-a6e4763f-24c1-41e0-ad5c-8912471a7862");katex.render("N", element, {displayMode:false,throwOnError:false});$	Atomic number density	m $var element = document.getElementById("moose-equation-bb3761a5-91bd-4043-b3bd-59e8edd0a12a");katex.render("^{\u22123}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-cae1e864-8340-41c6-a35c-09c78cb2bb7f");katex.render("\\rho", element, {displayMode:false,throwOnError:false});$	Density	g m $var element = document.getElementById("moose-equation-755cfcb7-df96-464f-8328-56b85e2f7597");katex.render("^{-3}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-b287ced0-2854-4eb6-b80f-44f3887822a1");katex.render("c_p", element, {displayMode:false,throwOnError:false});$	Specific heat	J kg $var element = document.getElementById("moose-equation-000824a4-6cb6-4b41-befa-c88be1f05e17");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ K $var element = document.getElementById("moose-equation-8815eba9-1c45-4b68-8b2f-4e829f7eb57c");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$
$var element = document.getElementById("moose-equation-a4943ae0-9da7-42e3-a35b-d2502d00c510");katex.render("k_T", element, {displayMode:false,throwOnError:false});$	Thermal conductivity	W m $var element = document.getElementById("moose-equation-38962d7c-94a3-424d-a2bd-884a4a2beaad");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$ K $var element = document.getElementById("moose-equation-c2934a97-fe02-43d9-9202-84d474c3581d");katex.render("^{-1}", element, {displayMode:false,throwOnError:false});$