diff --git a/test/tests/TEX_pipe_flow_navier_stokes_coupling/TEX_modeling_tritium_steady_state.i b/test/tests/TEX_pipe_flow_navier_stokes_coupling/TEX_modeling_tritium_steady_state.i
new file mode 100644
index 00000000..28d19625
--- /dev/null
+++ b/test/tests/TEX_pipe_flow_navier_stokes_coupling/TEX_modeling_tritium_steady_state.i
@@ -0,0 +1,203 @@
+# TMAP8 input file
+# Written by Pierre-Clément Simon - Idaho National Laboratory
+#
+# Published with:
+# P.-C. A. Simon,
+# "TBD"
+#
+# Info:
+# - This input file couples solid transport and liquid transport with Navier-Stokes
+# - It models the tritium extraction experiment (TEX) PbLi loop
+# - block 0 on the left corresponds to the solid
+# - block 1 on the right corresponds to the liquid
+# - Assumes steady state and the liquid velocity is null
+#
+#
+# Once TMAP8 is installed and built (see instructions on the TMAP8 website), run with
+# cd ~/projects/TMAP8/test/tests/TEX_pipe_flow_navier_stokes_coupling/
+# mpirun -np 1 ~/projects/TMAP8/tmap8-opt -i TEX_modeling_tritium_steady_state.i
+
+# Boundary conditions
+concentration_solid_external = 0.0
+concentration_liquid_internal = 1.0
+
+# Material properties
+diffusion_coefficient_PbLi = 1.0
+diffusion_coefficient_V = 1.0
+
+# Convection
+convection_transfer_coefficient = 1.0e8
+
+# file names
+output_name = TEX_modeling_tritium_out/TEX_modeling_tritium_out
+
+[Mesh]
+  [cmg]
+    type = CartesianMeshGenerator
+    dim = 2
+    dx = '0.5 0.5'
+    dy = '0.5 0.5'
+    ix = '5 10'
+    iy = '5 5'
+    subdomain_id = '0 1
+                    0 1'
+  []
+  [add_sideset0] # from 0 (solid) to 1 (liquid)
+    type = SideSetsBetweenSubdomainsGenerator
+    input = cmg
+    new_boundary = interface_0_1
+    primary_block = 0
+    paired_block = 1
+  []
+  [add_sideset1] # from 1 (liquid) to 0 (solid)
+    type = SideSetsBetweenSubdomainsGenerator
+    input = add_sideset0
+    new_boundary = interface_1_0
+    primary_block = 1
+    paired_block = 0
+  []
+[]
+
+[Variables]
+  [tritium_solid_fe] # finite element variable for tritium concentration in solid
+    block = 0
+  []
+  [tritium_liquid_fv] # finite volume variable for tritium concentration in liquid
+    type = MooseVariableFVReal
+    block = 1
+  []
+[]
+
+[Kernels]
+  [tritium_solid_fe_diffusion]
+    type = ADMatDiffusion
+    variable = tritium_solid_fe
+    diffusivity = ${diffusion_coefficient_V}
+  []
+[]
+
+[BCs]
+  [tritium_solid_fe_left]
+    type = ADDirichletBC
+    boundary = left
+    variable = tritium_solid_fe
+    value = ${concentration_solid_external}
+  []
+  [tritium_solid_fe_interface]
+    type = ADConvectiveHeatFluxBC
+    boundary = interface_0_1
+    variable = tritium_solid_fe
+    T_infinity_functor = tritium_liquid_fv
+    heat_transfer_coefficient_functor = ${convection_transfer_coefficient}
+  []
+[]
+
+[FVKernels]
+  [tritium_liquid_fv_diff]
+    type = FVDiffusion
+    variable = tritium_liquid_fv
+    coeff = ${diffusion_coefficient_PbLi}
+  []
+[]
+
+[FVBCs]
+  [tritium_liquid_fv_right]
+    type = FVDirichletBC
+    boundary = right
+    variable = tritium_liquid_fv
+    value = ${concentration_liquid_internal}
+  []
+  [tritium_liquid_fv_interface]
+    type = FVFunctorConvectiveHeatFluxBC
+    boundary = interface_1_0
+    variable = tritium_liquid_fv
+    T_bulk = tritium_liquid_fv
+    T_solid = tritium_solid_fe
+    heat_transfer_coefficient = ${convection_transfer_coefficient}
+    is_solid = false
+  []
+[]
+
+[Postprocessors]
+  [Tritium_amount_solid]
+    type = ElementIntegralVariablePostprocessor
+    variable = tritium_solid_fe
+    block = 0
+    execute_on = 'INITIAL TIMESTEP_END'
+  []
+  [Tritium_concentration_interface_solid]
+    type = PointValue
+    variable = tritium_solid_fe
+    point = '0.5 0 0'
+    execute_on = 'INITIAL TIMESTEP_BEGIN TIMESTEP_END LINEAR NONLINEAR'
+  []
+  [Tritium_amount_liquid]
+    type = ElementIntegralVariablePostprocessor
+    variable = tritium_liquid_fv
+    block = 1
+    execute_on = 'INITIAL TIMESTEP_END'
+  []
+  [Tritium_amount_total]
+    type = LinearCombinationPostprocessor
+    pp_names = 'Tritium_amount_solid Tritium_amount_liquid'
+    pp_coefs = '1   1'
+    execute_on = 'INITIAL TIMESTEP_END'
+  []
+  [dt]
+    type = TimestepSize
+  []
+[]
+
+[Positions]
+  [pos]
+    type = InputPositions
+    positions = '0.51 0.55 0
+                 0.61 0.55 0
+                 0.71 0.55 0
+                 0.81 0.55 0
+                 0.91 0.55 0
+                 0.99 0.55 0'
+  []
+[]
+
+[VectorPostprocessors]
+  [radial_profile_solid]
+    type = LineValueSampler
+    variable = 'tritium_solid_fe'
+    start_point = '0 0.5 0'
+    end_point = '0.5 0.5 0'
+    num_points = 10
+    sort_by = 'x'
+    execute_on = 'TIMESTEP_END FINAL'
+  []
+  [radial_profile_liquid]
+    type = PositionsFunctorValueSampler
+    functors = 'tritium_liquid_fv'
+    positions = 'pos'
+    sort_by = 'x'
+    execute_on = 'TIMESTEP_END FINAL'
+    # the finite volume variable is discontinuous at faces
+    # Note that we are not sampling on faces so it does not matter,
+    # we could set it to 'false' to have less checks when sampling
+    discontinuous = true
+  []
+[]
+
+[Preconditioning]
+  [Newtonlu]
+    type = SMP
+    full = true
+    solve_type = 'NEWTON'
+    petsc_options_iname = '-pc_type -sub_pc_type'
+    petsc_options_value = 'asm      lu          '
+  []
+[]
+
+[Executioner]
+  type = Steady
+[]
+
+[Outputs]
+  file_base = ${output_name}
+  csv = true
+[]
diff --git a/test/tests/TEX_pipe_flow_navier_stokes_coupling/linear_pipe_explicit.i b/test/tests/TEX_pipe_flow_navier_stokes_coupling/linear_pipe_explicit.i
new file mode 100644
index 00000000..3ac3cbe2
--- /dev/null
+++ b/test/tests/TEX_pipe_flow_navier_stokes_coupling/linear_pipe_explicit.i
@@ -0,0 +1,619 @@
+#######################################################################
+# Model Description
+# Vanadium Pipe Benchmark Model
+# Created by Samuel Walker on Apr23-2025
+#######################################################################
+
+
+#######################################################################
+# GEOMETRY and MESH
+#######################################################################
+
+
+#p_outlet = 1e5
+radius_inner = ${units 0.0117 m}
+radius_outer = ${units 0.0127 m}
+cylinder_height = ${units 1 m}
+#velocity_inlet = 0.465926 #2Kg/s
+#velocity_inlet = 0.698889 #3Kg/s
+velocity_inlet = 0.931852 #4Kg/s # units?
+Temperature_fluid = ${units 573.15 K}
+tritium_liquid_concentration_inlet = 1e-3 # units?
+R_gas = ${units 8.31446261815324 J/mol/K} # from PhysicalConstants.h
+#mu = 2.6
+#rho = 1.0
+advected_interp_method = 'upwind'
+
+#T_initial = 450
+
+#vanadium = '1'
+#fluid  = '2'
+
+[GlobalParams]
+  rhie_chow_user_object = 'rc'
+[]
+
+[Mesh]
+
+  # create a circle first with concentric layers (PbLi, Vanadium, Vacuum [not meshed])
+  [circle]
+    type = ConcentricCircleMeshGenerator
+    radii            = '${radius_inner} ${radius_outer}'
+    rings            = '4 4'
+    num_sectors      = 6
+    has_outer_square = false
+    preserve_volumes = false
+  []
+
+  # turn the circle into a cylinder with the corresponding heights for each block
+  # then assign new subdomains to be renamed later
+  [extrude]
+    type = AdvancedExtruderGenerator
+    input = circle
+    heights = '${cylinder_height}'
+    num_layers = '100'
+    direction = '0 0 1'
+    top_boundary    = 'outlet'
+    bottom_boundary = 'inlet'
+  []
+
+  # rename each of the blocks to make them the correct block for material assignment
+  [rename]
+    type = RenameBlockGenerator
+    input = extrude
+    old_block = '1 2'
+    new_block = 'fluid pipe'
+  []
+  [add_wall]
+    type = SideSetsBetweenSubdomainsGenerator
+    input = rename
+    primary_block = 'fluid'
+    paired_block = 'pipe'
+    new_boundary = 'wall'
+  []
+
+  [fix_sideset]
+    type = SideSetsAroundSubdomainGenerator
+    input = add_wall
+    new_boundary = 'fluid_outlet'
+    block = 'fluid'
+    #include_only_external_sides = 'true'
+    normal = '0 0 1'
+    #normals = '0.00898026 -0.00898026 1'
+    #new_boundary = 'wall_outlet'
+  []
+  [fix_sideset_2]
+    type = SideSetsAroundSubdomainGenerator
+    input = fix_sideset
+    new_boundary = 'fluid_inlet'
+    block = 'fluid'
+    include_only_external_sides = 'true'
+    normal = '0 0 -0.01'
+    #normals = '0.00898026 -0.00898026 1'
+    #new_boundary = 'wall_outlet'
+  []
+  [fix_sideset_3]
+    type = SideSetsAroundSubdomainGenerator
+    input = fix_sideset_2
+    new_boundary = 'solid_inlet'
+    block = 'pipe'
+    include_only_external_sides = 'true'
+    normal = '0 0 -0.01'
+    #normals = '0.00898026 -0.00898026 1'
+    #new_boundary = 'wall_outlet'
+  []
+  [fix_sideset_4]
+    type = SideSetsAroundSubdomainGenerator
+    input = fix_sideset_3
+    new_boundary = 'solid_outlet'
+    block = 'pipe'
+    #include_only_external_sides = 'true'
+    normal = '0 0 1'
+    #normals = '0.00898026 -0.00898026 1'
+    #new_boundary = 'wall_outlet'
+  []
+  [delete_side_sets]
+    type = BoundaryDeletionGenerator
+    input = fix_sideset_4
+    boundary_names = 'inlet outlet'
+    #normals = '0.00898026 -0.00898026 1'
+    #new_boundary = 'wall_outlet'
+  []
+[]
+
+###################################
+# TH
+###################################
+
+fluid_blocks = 'fluid'
+solid_blocks = 'pipe'
+
+[Problem]
+  linear_sys_names = 'u_system v_system w_system pressure_system tritium_liquid_system tritium_solid_system'
+  previous_nl_solution_required = true
+  kernel_coverage_check = false
+  material_coverage_check = false
+  fv_bcs_integrity_check = false
+  boundary_restricted_elem_integrity_check = false
+  #verbose_setup = true
+[]
+
+[UserObjects]
+  [rc]
+    type = RhieChowMassFlux
+    u = vel_x
+    v = vel_y
+    w = vel_z
+    pressure = pressure
+    rho = 'rho'
+    p_diffusion_kernel = p_diffusion
+    block = ${fluid_blocks}
+  []
+[]
+
+[Variables]
+  #TH Variables
+  [vel_x] # units ?
+    type = MooseLinearVariableFVReal
+    initial_condition = 0
+    block = ${fluid_blocks}
+    solver_sys = u_system
+  []
+  [vel_y] # units ?
+    type = MooseLinearVariableFVReal
+    initial_condition = 0
+    block = ${fluid_blocks}
+    solver_sys = v_system
+  []
+  [vel_z] # units ?
+    type = MooseLinearVariableFVReal
+    initial_condition = ${velocity_inlet}
+    block = ${fluid_blocks}
+    solver_sys = w_system
+  []
+  [pressure] # units ?
+    type = MooseLinearVariableFVReal
+    initial_condition = 0.2  # ?
+    block = ${fluid_blocks}
+    solver_sys = pressure_system
+  []
+
+  #Tritium Variables
+  [tritium_liquid]
+    type = MooseLinearVariableFVReal
+    initial_condition = 0.0
+    block = ${fluid_blocks}
+    solver_sys = 'tritium_liquid_system'
+  []
+  [tritium_solid]
+    type = MooseLinearVariableFVReal
+    initial_condition = 0.0
+    block = ${solid_blocks}
+    solver_sys = 'tritium_solid_system'
+  []
+[]
+
+# [FluidProperties]
+#   [fluid_properties_obj]
+#     type = LeadBismuthFluidProperties
+#   []
+# []
+
+[LinearFVKernels]
+  #Advection Kernels
+  [u_advection_stress]
+    type = LinearWCNSFVMomentumFlux
+    variable = vel_x
+    advected_interp_method = ${advected_interp_method}
+    mu = 'mu'
+    u = vel_x
+    v = vel_y
+    w = vel_z
+    momentum_component = 'x'
+    rhie_chow_user_object = 'rc'
+    use_nonorthogonal_correction = false
+  []
+  [v_advection_stress]
+    type = LinearWCNSFVMomentumFlux
+    variable = vel_y
+    advected_interp_method = ${advected_interp_method}
+    mu = 'mu'
+    u = vel_x
+    v = vel_y
+    w = vel_z
+    momentum_component = 'y'
+    rhie_chow_user_object = 'rc'
+    use_nonorthogonal_correction = false
+  []
+  [w_advection_stress]
+    type = LinearWCNSFVMomentumFlux
+    variable = vel_z
+    advected_interp_method = ${advected_interp_method}
+    mu = 'mu'
+    u = vel_x
+    v = vel_y
+    w = vel_z
+    momentum_component = 'z'
+    rhie_chow_user_object = 'rc'
+    use_nonorthogonal_correction = false
+  []
+
+  #Cross Terms
+  [u_pressure]
+    type = LinearFVMomentumPressure
+    variable = vel_x
+    pressure = pressure
+    momentum_component = 'x'
+  []
+  [v_pressure]
+    type = LinearFVMomentumPressure
+    variable = vel_y
+    pressure = pressure
+    momentum_component = 'y'
+  []
+  [w_pressure]
+    type = LinearFVMomentumPressure
+    variable = vel_z
+    pressure = pressure
+    momentum_component = 'z'
+  []
+
+  #Pressure Kernels
+  [p_diffusion]
+    type = LinearFVAnisotropicDiffusion
+    variable = pressure
+    diffusion_tensor = Ainv
+    use_nonorthogonal_correction = false
+  []
+  [HbyA_divergence]
+    type = LinearFVDivergence
+    variable = pressure
+    face_flux = HbyA
+    force_boundary_execution = true
+  []
+
+  #Tritium Kernels
+  #tritium_liquid Kernels
+  [tritium_liquid_advection]
+    type = LinearFVScalarAdvection
+    variable = tritium_liquid
+    #rhie_chow_user_object = rhie_chow_user_object
+    #vel = 'vel_x_mat vel_y_mat'
+    #vel = 'vel_x vel_y'
+    #rho = ${rho}
+    block = ${fluid_blocks}
+  []
+  [tritium_liquid_diffusion]
+    type = LinearFVDiffusion
+    diffusion_coeff = 'diffusivity_tritium_liquid'
+    variable = tritium_liquid
+    block = ${fluid_blocks}
+    use_nonorthogonal_correction = false
+  []
+  #tritium_solid Kernels - TMAP8 Here
+  [tritium_solid_diffusion]
+    type = LinearFVDiffusion
+    diffusion_coeff = 'diffusivity_tritium_solid'
+    variable = tritium_solid
+    #block = ${fluid_blocks}
+    use_nonorthogonal_correction = false
+  []
+[]
+
+[LinearFVBCs]
+  #Inlet TH BCs
+   [inlet-u]
+    type = LinearFVAdvectionDiffusionFunctorDirichletBC
+    boundary = 'fluid_inlet'
+    variable = vel_x
+    functor = '0.0'
+  []
+  [inlet-v]
+    type = LinearFVAdvectionDiffusionFunctorDirichletBC
+    boundary = 'fluid_inlet'
+    variable = vel_y
+    functor = '0.0'
+  []
+  [inlet-w]
+    type = LinearFVAdvectionDiffusionFunctorDirichletBC
+    boundary = 'fluid_inlet'
+    variable = vel_z
+    functor = 'InletLaminar'
+  []
+  [walls-u]
+    type = LinearFVAdvectionDiffusionFunctorDirichletBC
+    boundary = 'wall'
+    variable = vel_x
+    functor = 0.0
+  []
+  [walls-v]
+    type = LinearFVAdvectionDiffusionFunctorDirichletBC
+    boundary = 'wall'
+    variable = vel_y
+    functor = 0.0
+  []
+  [walls-w]
+    type = LinearFVAdvectionDiffusionFunctorDirichletBC
+    boundary = 'wall'
+    variable = vel_z
+    functor = 0.0
+  []
+
+  #Outlet TH BCs
+  [outlet_p]
+    type = LinearFVAdvectionDiffusionFunctorDirichletBC
+    boundary = 'fluid_outlet'
+    variable = pressure
+    functor = 1.2
+  []
+  [outlet_u]
+    type = LinearFVAdvectionDiffusionOutflowBC
+    variable = vel_x
+    use_two_term_expansion = false
+    boundary = 'fluid_outlet'
+  []
+  [outlet_v]
+    type = LinearFVAdvectionDiffusionOutflowBC
+    variable = vel_y
+    use_two_term_expansion = false
+    boundary = 'fluid_outlet'
+  []
+  [outlet_w]
+    type = LinearFVAdvectionDiffusionOutflowBC
+    variable = vel_z
+    use_two_term_expansion = false
+    boundary = 'fluid_outlet'
+  []
+
+  #Tritium BCs
+  #tritium_liquid inlet and outlets
+  [tritium_liquid_inlet]
+    type = LinearFVAdvectionDiffusionFunctorDirichletBC
+    boundary = 'fluid_inlet'
+    functor = ${tritium_liquid_concentration_inlet}
+    variable = tritium_liquid
+  []
+  [tritium_liquid_outlet]
+    type = LinearFVAdvectionDiffusionOutflowBC
+    boundary = 'fluid_outlet'
+    variable = tritium_liquid
+    use_two_term_expansion = false
+  []
+
+  #Mass Transfer at Interface. -TMPAP8 Here
+  [tritium_solid_source]
+    type = LinearFVConvectiveHeatTransferBC
+    boundary = 'wall'
+    T_fluid = tritium_liquid
+    T_solid = tritium_solid
+    h = 1 # mass transfer exchange rate at interface - tune
+    #functor = 1e5
+    variable = tritium_solid
+    #is_solid = true
+    #block = ${fluid_blocks}
+  []
+  [tritium_liquid_loss]
+    type = LinearFVConvectiveHeatTransferBC
+    boundary = 'wall'
+    T_fluid = tritium_liquid
+    T_solid = tritium_solid
+    h = 1
+    #functor = 1e5
+    variable = tritium_liquid
+    #is_solid = false
+    #block = ${fluid_blocks}
+  []
+
+  #Solid Outlets from diffusion
+  [tritium_solid_outer]
+    type = LinearFVAdvectionDiffusionOutflowBC
+    boundary = 'outer'
+    #value = 0
+    variable = tritium_solid
+    use_two_term_expansion = false
+    #block = ${fluid_blocks}
+  []
+  [tritium_solid_outer_2]
+    type = LinearFVAdvectionDiffusionOutflowBC
+    boundary = 'solid_outlet'
+    #value = 0
+    variable = tritium_solid
+    use_two_term_expansion = false
+    #block = ${fluid_blocks}
+  []
+  [tritium_solid_outer_3]
+    type = LinearFVAdvectionDiffusionOutflowBC
+    boundary = 'solid_inlet'
+    #value = 0
+    variable = tritium_solid
+    use_two_term_expansion = false
+    #block = ${fluid_blocks}
+  []
+[]
+
+[FunctorMaterials]
+  # [fluid_props_to_mat_props]
+  #   type = GeneralFunctorFluidProps
+  #   pressure = 'pressure'
+  #   Temperature_fluid = ${Temperature_fluid}
+  #   #Temperature_fluid = 450
+  #   speed = 'speed'
+  #   characteristic_length = 1.0
+  #   fp = fluid_properties_obj
+  #   porosity = 1.0
+  #   block = ${fluid_blocks}
+  # []
+  # [misc_mat_props]
+  #   type = ADGenericFunctorMaterial
+  #   prop_names = 'pipe_mat'
+  #   prop_value = 1
+  #   block = 'pipe'
+  #   execute_on = 'INITIAL'
+  # []
+  [rho]
+    type = ParsedFunctorMaterial
+    property_name = 'rho'
+    expression = '10520-1.19051*${Temperature_fluid}'
+    block = 'fluid'
+    #execute_on = 'Always'
+  []
+  [mu]
+    type = ParsedFunctorMaterial
+    property_name = 'mu'
+    expression = '0.187*10^-3*exp(1400/${Temperature_fluid})+1*exp(-t/10)'
+    block = 'fluid'
+    #exe
+  []
+  [diffusivity_tritium_liquid] # PbLi # ?
+    type = ParsedFunctorMaterial
+    property_name = 'diffusivity_tritium_liquid'
+    expression = '4.03*10^-8*exp(-19500/(${R_gas}*${Temperature_fluid}))'
+    block = 'fluid'
+    #exe
+  []
+  [diffusivity_tritium_solid] # Vanadium # ?
+    type = ParsedFunctorMaterial
+    property_name = 'diffusivity_tritium_solid'
+    expression = '2.0*10^-8*exp(-4200/(${R_gas}*${Temperature_fluid}))'
+    block = 'pipe'
+    #exe
+  []
+  # [rho]
+  #   type = GenericFunctorMaterial
+  #   prop_names = 'rho'
+  #   prop_value = 'rho_calc'
+  #   block = 'fluid'
+  #   #execute_on = 'INITIAL'
+  # []
+  # [mu]
+  #   type = GenericFunctorMaterial
+  #   prop_names = 'mu'
+  #   prop_value = 'mu_calc'
+  #   block = 'fluid'
+  #   #execute_on = 'INITIAL'
+  # []
+[]
+
+[Functions]
+  [InletLaminar]
+    type = ParsedFunction
+    expression = 2.*${velocity_inlet}*(1.0-(y/${radius_inner})^2.0-(x/${radius_inner})^2)
+  []
+  [rho_calc]
+    type = ParsedFunction
+    symbol_names = 'T'
+    symbol_values = '${Temperature_fluid}'
+    expression = 10520.35-1.19051*T
+  []
+  [mu_calc]
+    type = ParsedFunction
+    symbol_names = 'T'
+    symbol_values = '${Temperature_fluid}'
+    expression = 0.187*10^-3*exp(1400/T)
+  []
+[]
+
+[Executioner]
+  type = PIMPLE
+  momentum_l_abs_tol = 1e-13
+  pressure_l_abs_tol = 1e-13
+  #energy_l_abs_tol = 1e-13
+  passive_scalar_l_abs_tol = 1e-13
+  momentum_l_tol = 1e-13
+  pressure_l_tol = 1e-13
+  #energy_l_tol = 1e-12
+  passive_scalar_l_tol = 1e-13
+  rhie_chow_user_object = 'rc'
+  momentum_systems = 'u_system v_system w_system'
+  pressure_system = 'pressure_system'
+  #energy_system = 'energy_system'
+  passive_scalar_systems = 'tritium_liquid_system tritium_solid_system'
+  momentum_equation_relaxation = 0.8
+  pressure_variable_relaxation = 0.3
+  #energy_equation_relaxation = 0.9
+  passive_scalar_equation_relaxation = '0.8 0.8'
+  num_iterations = 100
+  pressure_absolute_tolerance = 1e-11
+  momentum_absolute_tolerance = 1e-11
+  #energy_absolute_tolerance = 1e-11
+  passive_scalar_absolute_tolerance = '1e-11 1e-11'
+  momentum_petsc_options_iname = '-pc_type -pc_hypre_type'
+  momentum_petsc_options_value = 'hypre boomeramg'
+  pressure_petsc_options_iname = '-pc_type -pc_hypre_type'
+  pressure_petsc_options_value = 'hypre boomeramg'
+  #energy_petsc_options_iname = '-pc_type -pc_hypre_type'
+  #energy_petsc_options_value = 'hypre boomeramg'
+  passive_scalar_petsc_options_iname = '-pc_type -pc_hypre_type'
+  passive_scalar_petsc_options_value = 'hypre boomeramg'
+  print_fields = false
+  continue_on_max_its = true
+
+ [TimeStepper]
+    type = IterationAdaptiveDT
+    dt = 0.01
+    optimal_iterations = 10
+    iteration_window = 2
+    growth_factor = 2
+    cutback_factor = 0.5
+  []
+
+  end_time = 1e10
+  steady_state_detection = true
+  steady_state_tolerance = 1e-12
+[]
+
+[Postprocessors]
+  [flow_outlet]
+    type = VolumetricFlowRate
+    boundary = 'fluid_outlet'
+    vel_x = vel_x
+    vel_y = vel_y
+    vel_z = vel_z
+    advected_quantity = 1
+  []
+  [mass_flow_outlet]
+    type = MassFluxWeightedFlowRate
+    boundary = 'fluid_outlet'
+    vel_x = vel_x
+    vel_y = vel_y
+    vel_z = vel_z
+    density = 'rho'
+    advected_quantity = 'rho'
+  []
+  [mdot_outlet]
+    type = VolumetricFlowRate
+    boundary = 'fluid_outlet'
+    vel_x = vel_x
+    vel_y = vel_y
+    vel_z = vel_z
+    #density = 'rho'
+    advected_quantity = 'rho'
+  []
+  # [rho_check]
+  #   type = ElementIntegralFunctorPostprocessor
+  #   functor = 'rho'
+  #   block = 'fluid'
+  #   #pp_names = 'rho'
+  #   #vel_x = vel_x
+  #   #vel_y = vel_y
+  #   #vel_z = vel_z
+  #   #density = 'rho'
+  #   #advected_quantity = 'rho'
+  # []
+  [tritium_liquid_avg]
+    type = ElementAverageValue
+    variable = tritium_liquid
+    block = 'fluid'
+    #pp_names = 'rho'
+    #vel_x = vel_x
+    #vel_y = vel_y
+    #vel_z = vel_z
+    #density = 'rho'
+    #advected_quantity = 'rho'
+  []
+[]
+
+[Outputs]
+  csv = true
+  exodus = true
+[]