MultiAppPostprocessorInterpolationTransfer

Transfer postprocessor data from sub-application into field data on the master application.

Overview

Performs a transfer of a PostProcessor value from sub-applications to a field variable on the master application using interpolation based on the location of each sub-application.

Example Input File Syntax

The following input file snippet demonstrates the use of the MultiAppPostprocessorInterpolationTransfer to interpolate the average value from two sub-applications to a field variable in the master application.

[Transfers]
  [pp_transfer]
    direction = from_multiapp
    postprocessor = average
    variable = from_sub
    type = MultiAppPostprocessorInterpolationTransfer
    multi_app = sub
  []
[]

Input Parameters

directionWhether this Transfer will be 'to' or 'from' a MultiApp, or bidirectional, by providing both FROM_MULTIAPP and TO_MULTIAPP.
C++ Type:MultiMooseEnum
Options:to_multiapp, from_multiapp
Description:Whether this Transfer will be 'to' or 'from' a MultiApp, or bidirectional, by providing both FROM_MULTIAPP and TO_MULTIAPP.
multi_appThe name of the MultiApp to use.
C++ Type:MultiAppName
Options:
Description:The name of the MultiApp to use.
postprocessorThe Postprocessor to interpolate.
C++ Type:PostprocessorName
Options:
Description:The Postprocessor to interpolate.
variableThe auxiliary variable to store the transferred values in.
C++ Type:AuxVariableName
Options:
Description:The auxiliary variable to store the transferred values in.

Required Parameters

check_multiapp_execute_onTrueWhen false the check between the multiapp and transfer execute on flags is not preformed.
Default:True
C++ Type:bool
Options:
Description:When false the check between the multiapp and transfer execute on flags is not preformed.
displaced_source_meshFalseWhether or not to use the displaced mesh for the source mesh.
Default:False
C++ Type:bool
Options:
Description:Whether or not to use the displaced mesh for the source mesh.
displaced_target_meshFalseWhether or not to use the displaced mesh for the target mesh.
Default:False
C++ Type:bool
Options:
Description:Whether or not to use the displaced mesh for the target mesh.
execute_onSAME_AS_MULTIAPPThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM, SAME_AS_MULTIAPP.
Default:SAME_AS_MULTIAPP
C++ Type:ExecFlagEnum
Options:NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM, SAME_AS_MULTIAPP
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM, SAME_AS_MULTIAPP.
interp_typeinverse_distanceThe algorithm to use for interpolation.
Default:inverse_distance
C++ Type:MooseEnum
Options:inverse_distance, radial_basis
Description:The algorithm to use for interpolation.
num_points3The number of nearest points to use for interpolation.
Default:3
C++ Type:unsigned int
Options:
Description:The number of nearest points to use for interpolation.
power2The polynomial power to use for calculation of the decay in the interpolation.
Default:2
C++ Type:double
Options:
Description:The polynomial power to use for calculation of the decay in the interpolation.
radius-1Radius to use for radial_basis interpolation. If negative then the radius is taken as the max distance between points.
Default:-1
C++ Type:double
Options:
Description:Radius to use for radial_basis interpolation. If negative then the radius is taken as the max distance between points.

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Options:
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

Input Files

(tutorials/darcy_thermo_mech/step10_multiapps/problems/step10.i)
(modules/combined/examples/geochem-porous_flow/geotes_2D/exchanger_un_quartz.i)
(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/multilevel_master.i)
(tutorials/tutorial02_multiapps/step03_coupling/02_master_picard.i)
(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/master2_quad.i)
(tutorials/tutorial02_multiapps/step02_transfers/04_master_multiscale.i)
(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/radial_master.i)
(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/master.i)
(tutorials/tutorial02_multiapps/step03_coupling/01_master.i)
(modules/combined/examples/geochem-porous_flow/geotes_2D/exchanger.i)
(tutorials/tutorial02_multiapps/step03_coupling/03_master_subcycling_picard.i)
(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/multilevel_sub.i)
(modules/combined/examples/geochem-porous_flow/geotes_weber_tensleep/exchanger.i)
(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/master_quad.i)

(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/master.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [u]
  []
[]

[AuxVariables]
  [from_sub]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
  dt = 1

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [sub]
    positions = '0.2 0.2 0 0.7 0.7 0'
    type = TransientMultiApp
    app_type = MooseTestApp
    input_files = 'sub0.i sub1.i'
  []
[]

[Transfers]
  [pp_transfer]
    direction = from_multiapp
    postprocessor = average
    variable = from_sub
    type = MultiAppPostprocessorInterpolationTransfer
    multi_app = sub
  []
[]

(tutorials/darcy_thermo_mech/step10_multiapps/problems/step10.i)

[GlobalParams]
  displacements = 'disp_r disp_z'
[]

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 100
  ymax = 0.304 # Length of test chamber
  xmax = 0.0257 # Test chamber radius
[]

[Variables]
  [pressure]
  []
  [temperature]
    initial_condition = 300 # Start at room temperature
  []
[]

[AuxVariables]
  [k_eff]
    initial_condition = 15.0 # water at 20C
  []
  [velocity]
    order = CONSTANT
    family = MONOMIAL_VEC
  []
[]

[Modules/TensorMechanics/Master]
  [all]
    # This block adds all of the proper Kernels, strain calculators, and Variables
    # for TensorMechanics in the correct coordinate system (autodetected)
    add_variables = true
    strain = FINITE
    eigenstrain_names = eigenstrain
    use_automatic_differentiation = true
    generate_output = 'vonmises_stress elastic_strain_xx elastic_strain_yy strain_xx strain_yy'
  []
[]

[Kernels]
  [darcy_pressure]
    type = DarcyPressure
    variable = pressure
  []
  [heat_conduction]
    type = ADHeatConduction
    variable = temperature
  []
  [heat_conduction_time_derivative]
    type = ADHeatConductionTimeDerivative
    variable = temperature
  []
  [heat_convection]
    type = DarcyAdvection
    variable = temperature
    pressure = pressure
  []
[]

[AuxKernels]
  [velocity]
    type = DarcyVelocity
    variable = velocity
    execute_on = timestep_end
    pressure = pressure
  []
[]

[BCs]
  [inlet]
    type = DirichletBC
    variable = pressure
    boundary = bottom
    value = 4000 # (Pa) From Figure 2 from paper.  First data point for 1mm spheres.
  []
  [outlet]
    type = DirichletBC
    variable = pressure
    boundary = top
    value = 0 # (Pa) Gives the correct pressure drop from Figure 2 for 1mm spheres
  []
  [inlet_temperature]
    type = FunctionDirichletBC
    variable = temperature
    boundary = bottom
    function = 'if(t<0,350+50*t,350)'
  []
  [outlet_temperature]
    type = HeatConductionOutflow
    variable = temperature
    boundary = top
  []
  [hold_inlet]
    type = DirichletBC
    variable = disp_z
    boundary = bottom
    value = 0
  []
  [hold_center]
    type = DirichletBC
    variable = disp_r
    boundary = left
    value = 0
  []
  [hold_outside]
    type = DirichletBC
    variable = disp_r
    boundary = right
    value = 0
  []
[]

[Materials]
  viscosity_file = data/water_viscosity.csv
  density_file = data/water_density.csv
  specific_heat_file = data/water_specific_heat.csv
  thermal_expansion_file = data/water_thermal_expansion.csv
  [column]
    type = PackedColumn
    temperature = temperature
    radius = 1
    thermal_conductivity = k_eff # Use the AuxVariable instead of calculating
    fluid_viscosity_file = ${viscosity_file}
    fluid_density_file = ${density_file}
    fluid_specific_heat_file = ${specific_heat_file}
    fluid_thermal_expansion_file = ${thermal_expansion_file}
  []

  [elasticity_tensor]
    type = ADComputeIsotropicElasticityTensor
    youngs_modulus = 200e9 # (Pa) from wikipedia
    poissons_ratio = .3 # from wikipedia
  []
  [elastic_stress]
    type = ADComputeFiniteStrainElasticStress
  []
  [thermal_strain]
    type = ADComputeThermalExpansionEigenstrain
    stress_free_temperature = 300
    thermal_expansion_coeff = 1e-6
    eigenstrain_name = eigenstrain
    temperature = temperature
  []
[]

[Postprocessors]
  [average_temperature]
    type = ElementAverageValue
    variable = temperature
  []
[]

[Executioner]
  type = Transient
  start_time = -1
  end_time = 200
  steady_state_tolerance = 1e-7
  steady_state_detection = true
  dt = 0.25
  solve_type = PJFNK
  automatic_scaling = true
  compute_scaling_once = false
  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre boomeramg 500'
  line_search = none
  [TimeStepper]
    type = FunctionDT
    function = 'if(t<0,0.1,0.25)'
  []
[]


[MultiApps]
  [micro]
    type = TransientMultiApp
    app_type = DarcyThermoMechApp
    positions = '0.01285 0.0    0
                 0.01285 0.0608 0
                 0.01285 0.1216 0
                 0.01285 0.1824 0
                 0.01285 0.2432 0
                 0.01285 0.304  0'
    input_files = step10_micro.i
    execute_on = 'timestep_end'
  []
[]

[Transfers]
  [keff_from_sub]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = from_multiapp
    multi_app = micro
    variable = k_eff
    power = 1
    postprocessor = k_eff
    execute_on = 'timestep_end'
  []
  [temperature_to_sub]
    type = MultiAppVariableValueSamplePostprocessorTransfer
    direction = to_multiapp
    multi_app = micro
    source_variable = temperature
    postprocessor = temperature_in
    execute_on = 'timestep_end'
  []
[]

[Controls]
  [multiapp]
    type = TimePeriod
    disable_objects = 'MultiApps::micro Transfers::keff_from_sub Transfers::temperature_to_sub'
    start_time = '0'
    execute_on = 'initial'
  []
[]

[Outputs]
  [out]
    type = Exodus
    elemental_as_nodal = true
  []
[]

(modules/combined/examples/geochem-porous_flow/geotes_2D/exchanger_un_quartz.i)

# Model of the heat-exchanger
# The input fluid to the heat exchanger is determined by AuxVariables called production_temperature, production_rate_Na, production_rate_Cl, production_rate_SiO2 and production_rate_H2O.  These come from Postprocessors in the porous-flow simulation that measure the fluid composition at the production well.
# Given the input fluid, the exchanger cools/heats the fluid, removing any precipitates, and injects fluid back to the porous-flow simulation at temperature output_temperature and composition given by massfrac_Na, etc.
# In the absence of data concerning Quartz precipitation rates in heat exchangers, do not treat Quartz as kinetic
[GlobalParams]
  point = '0 0 0'
  reactor = reactor
[]

[TimeDependentReactionSolver]
  model_definition = definition
  include_moose_solve = false
  geochemistry_reactor_name = reactor
  charge_balance_species = "Cl-"
  swap_out_of_basis = "SiO2(aq)"
  swap_into_basis = "QuartzUnlike"
  constraint_species = "H2O              Na+                Cl-                QuartzUnlike"
  constraint_value = "  1.0E-2           0.1E-2             0.1E-2             1E-10"
  constraint_meaning = "kg_solvent_water bulk_composition bulk_composition free_mineral"
  constraint_unit = "   kg               moles              moles              moles"
  initial_temperature = 50.0
  mode = 4
  temperature = 200
  cold_temperature = 40.0
  source_species_names = 'H2O    Na+   Cl-   SiO2(aq)'
  source_species_rates = 'production_rate_H2O production_rate_Na production_rate_Cl production_rate_SiO2'
  ramp_max_ionic_strength_initial = 0 # max_ionic_strength in such a simple problem does not need ramping
  add_aux_pH = false # there is no H+ in this system
  evaluate_kinetic_rates_always = true # implicit time-marching used for stability
  execute_console_output_on = '' # only CSV output used in this example
[]

[UserObjects]
  [definition]
    type = GeochemicalModelDefinition
    database_file = "small_database.json"
    basis_species = "H2O SiO2(aq) Na+ Cl-"
    equilibrium_minerals = "QuartzUnlike"
  []
[]


[Executioner]
  type = Transient
  dt = 1E5
  end_time = 2E6 #7.76E6 # 90 days
[]

[AuxVariables]
  [production_temperature]
    initial_condition = 50 # the production_T Transfer lags one timestep behind for some reason, so give this a reasonable initial condition
  []
  [transported_H2O]
  []
  [transported_Na]
  []
  [transported_Cl]
  []
  [transported_SiO2]
  []
  [transported_mass]
  []
  [massfrac_H2O]
  []
  [massfrac_Na]
  []
  [massfrac_Cl]
  []
  [massfrac_SiO2]
  []
  [dumped_quartz]
  []
  [production_rate_H2O]
    initial_condition = 5.518533e+01 # the production_H2O Transfer lags one timestep behind for some reason (when the porous_flow simulation has finished, it correctly computes mole_rate_H2O_produced, but the Transfer gets the mole_rate_H2O_produced from the previous timestep), so give this a reasonable initial condition, otherwise this will be zero at the start of the simulation!
  []
  [production_rate_Na]
    initial_condition = 9.943302e-02
  []
  [production_rate_Cl]
    initial_condition = 9.943302e-02
  []
  [production_rate_SiO2]
    initial_condition = 2.340931e-04
  []
[]
[AuxKernels]
  [transported_H2O]
    type = GeochemistryQuantityAux
    variable = transported_H2O
    species = H2O
    quantity = transported_moles_in_original_basis
  []
  [transported_Na]
    type = GeochemistryQuantityAux
    variable = transported_Na
    species = Na+
    quantity = transported_moles_in_original_basis
  []
  [transported_Cl]
    type = GeochemistryQuantityAux
    variable = transported_Cl
    species = Cl-
    quantity = transported_moles_in_original_basis
  []
  [transported_SiO2]
    type = GeochemistryQuantityAux
    variable = transported_SiO2
    species = 'SiO2(aq)'
    quantity = transported_moles_in_original_basis
  []
  [transported_mass]
    type = ParsedAux
    args = 'transported_H2O transported_Na transported_Cl transported_SiO2'
    variable = transported_mass
    function = 'transported_H2O * 18.0152 + transported_Na * 22.9898 + transported_Cl * 35.453 + transported_SiO2 * 60.0843'
  []
  [massfrac_H2O]
    type = ParsedAux
    args = 'transported_mass transported_H2O'
    variable = massfrac_H2O
    function = '18.0152 * transported_H2O / transported_mass'
  []
  [massfrac_Na]
    type = ParsedAux
    args = 'transported_mass transported_Na'
    variable = massfrac_Na
    function = '22.9898 * transported_Na / transported_mass'
  []
  [massfrac_Cl]
    type = ParsedAux
    args = 'transported_mass transported_Cl'
    variable = massfrac_Cl
    function = '35.453 * transported_Cl / transported_mass'
  []
  [massfrac_SiO2]
    type = ParsedAux
    args = 'transported_mass transported_SiO2'
    variable = massfrac_SiO2
    function = '60.0843 * transported_SiO2 / transported_mass'
  []
  [dumped_quartz]
    type = GeochemistryQuantityAux
    variable = dumped_quartz
    species = QuartzUnlike
    quantity = moles_dumped
  []
[]
[Postprocessors]
  [cumulative_moles_precipitated_quartz]
    type = PointValue
    variable = dumped_quartz
  []
  [production_temperature]
    type = PointValue
    variable = production_temperature
  []
  [mass_heated_this_timestep]
    type = PointValue
    variable = transported_mass
  []
[]
[Outputs]
  csv = true
[]

[MultiApps]
  [porous_flow_sim]
    type = TransientMultiApp
    input_files = porous_flow.i
    cli_args = 'MultiApps/react/input_files=aquifer_un_quartz_geochemistry.i'
    execute_on = 'timestep_end'
  []
[]

[Transfers]
  [injection_T]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'solution_temperature'
    variable = 'injection_temperature'
  []
  [injection_Na]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Na'
    variable = 'injection_rate_massfrac_Na'
  []
  [injection_Cl]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Cl'
    variable = 'injection_rate_massfrac_Cl'
  []
  [injection_SiO2]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_SiO2'
    variable = 'injection_rate_massfrac_SiO2'
  []
  [injection_H2O]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_H2O'
    variable = 'injection_rate_massfrac_H2O'
  []

  [production_T]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = production_temperature
    variable = production_temperature
  []
  [production_Na]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Na_produced
    variable = production_rate_Na
  []
  [production_Cl]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Cl_produced
    variable = production_rate_Cl
  []
  [production_SiO2]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_SiO2_produced
    variable = production_rate_SiO2
  []
  [production_H2O]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_H2O_produced
    variable = production_rate_H2O
  []
[]

(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/multilevel_master.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./sub_average]
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
[]

[BCs]
  [./left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  [../]
  [./right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
  dt = 1

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [./sub]
    type = TransientMultiApp
    app_type = MooseTestApp
    positions = '0 0 0 0.5 0.5 0'
    input_files = multilevel_sub.i
  [../]
[]

[Transfers]
  [./sub_average]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = from_multiapp
    multi_app = sub
    variable = sub_average
    postprocessor = sub_average
  [../]
[]

(tutorials/tutorial02_multiapps/step03_coupling/02_master_picard.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [u]
  []
[]

[AuxVariables]
  [vt]
  []
[]

[Kernels]
  [diff]
    type = MatDiffusion
    variable = u
  []
  [force]
    type = BodyForce
    variable = u
    value = 1.
  []
  [td]
    type = TimeDerivative
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  []
[]

[Materials]
  [diff]
    type = ParsedMaterial
    f_name = D
    args = 'vt'
    function = 'vt'
  []
[]

[Executioner]
  type = Transient
  end_time = 2
  dt = 0.2

  picard_max_its = 10

  nl_abs_tol = 1e-10
  picard_rel_tol = 1e-6
  picard_abs_tol = 1e-10

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [micro]
    type = TransientMultiApp
    positions = '0.15 0.15 0  0.45 0.45 0  0.75 0.75 0'
    input_files = '02_sub_picard.i'
    execute_on = timestep_end
    output_in_position = true
  []
[]

[Transfers]
  [push_u]
    type = MultiAppVariableValueSampleTransfer
    multi_app = micro
    direction = to_multiapp
    source_variable = u
    variable = ut
  []

  [pull_v]
    type = MultiAppPostprocessorInterpolationTransfer
    multi_app = micro
    direction = from_multiapp
    variable = vt
    postprocessor = average_v
  []
[]

(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/master2_quad.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./pp_aux]
  [../]
[]

[Kernels]
  [./diff]
    type = CoefDiffusion
    variable = u
    coef = 0.1
  [../]
  [./time]
    type = TimeDerivative
    variable = u
  [../]
[]

[BCs]
  [./left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  [../]
  [./right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 20
  dt = 0.1
  solve_type = PJFNK
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [./quad]
    type = TransientMultiApp
    app_type = MooseTestApp
    positions = '0.1 0.1 0 0.9 0.1 0 0.1 0.9 0 0.9 0.9 0'
    input_files = 'quad_sub1.i'
  [../]
[]

[Transfers]
  [./sub_to_master_pp]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = from_multiapp
    multi_app = quad
    variable = pp_aux
    postprocessor = pp
  [../]
[]

(tutorials/tutorial02_multiapps/step02_transfers/04_master_multiscale.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [u]
  []
[]

[AuxVariables]
  [vt]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = u
  []
  [force]
    type = BodyForce
    variable = u
    value = 1.
  []
  [td]
    type = TimeDerivative
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  []
[]

[Executioner]
  type = Transient
  end_time = 2
  dt = 0.2

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [micro]
    type = TransientMultiApp
    positions = '0.15 0.15 0  0.45 0.45 0  0.75 0.75 0'
    input_files = '04_sub_multiscale.i'
    cli_args = 'BCs/right/value=1 BCs/right/value=2 BCs/right/value=3'
    execute_on = timestep_end
    output_in_position = true
  []
[]

[Transfers]
  [push_u]
    type = MultiAppVariableValueSampleTransfer
    multi_app = micro
    direction = to_multiapp
    source_variable = u
    variable = ut
  []

  [pull_v]
    type = MultiAppPostprocessorInterpolationTransfer
    multi_app = micro
    direction = from_multiapp
    variable = vt
    postprocessor = average_v
  []
[]

(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/radial_master.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./from_sub]
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
[]

[BCs]
  [./left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  [../]
  [./right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
  dt = 1

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [./sub]
    positions = '0.2 0.2 0 0.7 0.7 0'
    type = TransientMultiApp
    app_type = MooseTestApp
    input_files = 'sub0.i sub1.i'
  [../]
[]

[Transfers]
  [./pp_transfer]
    direction = from_multiapp
    postprocessor = average
    variable = from_sub
    type = MultiAppPostprocessorInterpolationTransfer
    multi_app = sub
    interp_type = radial_basis
    radius = 1.5
  [../]
[]

(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/master.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [u]
  []
[]

[AuxVariables]
  [from_sub]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
  dt = 1

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [sub]
    positions = '0.2 0.2 0 0.7 0.7 0'
    type = TransientMultiApp
    app_type = MooseTestApp
    input_files = 'sub0.i sub1.i'
  []
[]

[Transfers]
  [pp_transfer]
    direction = from_multiapp
    postprocessor = average
    variable = from_sub
    type = MultiAppPostprocessorInterpolationTransfer
    multi_app = sub
  []
[]

(tutorials/tutorial02_multiapps/step03_coupling/01_master.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [u]
  []
[]

[AuxVariables]
  [vt]
  []
[]

[Kernels]
  [diff]
    type = MatDiffusion
    variable = u
  []
  [force]
    type = BodyForce
    variable = u
    value = 1.
  []
  [td]
    type = TimeDerivative
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  []
[]

[Materials]
  [diff]
    type = ParsedMaterial
    f_name = D
    args = 'vt'
    function = 'vt'
  []
[]

[Executioner]
  type = Transient
  end_time = 2
  dt = 0.2

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [micro]
    type = TransientMultiApp
    positions = '0.15 0.15 0  0.45 0.45 0  0.75 0.75 0'
    input_files = '01_sub.i'
    execute_on = timestep_end
    output_in_position = true
  []
[]

[Transfers]
  [push_u]
    type = MultiAppVariableValueSampleTransfer
    multi_app = micro
    direction = to_multiapp
    source_variable = u
    variable = ut
  []

  [pull_v]
    type = MultiAppPostprocessorInterpolationTransfer
    multi_app = micro
    direction = from_multiapp
    variable = vt
    postprocessor = average_v
  []
[]

(modules/combined/examples/geochem-porous_flow/geotes_2D/exchanger.i)

# Model of the heat-exchanger
# The input fluid to the heat exchanger is determined by AuxVariables called production_temperature, production_rate_Na, production_rate_Cl, production_rate_SiO2 and production_rate_H2O.  These come from Postprocessors in the porous-flow simulation that measure the fluid composition at the production well.
# Given the input fluid, the exchanger cools/heats the fluid, removing any precipitates, and injects fluid back to the porous-flow simulation at temperature output_temperature and composition given by massfrac_Na, etc.
# In the absence of data concerning Quartz precipitation rates in heat exchangers, do not treat Quartz as kinetic
[GlobalParams]
  point = '0 0 0'
  reactor = reactor
[]

[TimeDependentReactionSolver]
  model_definition = definition
  include_moose_solve = false
  geochemistry_reactor_name = reactor
  charge_balance_species = "Cl-"
  swap_out_of_basis = "SiO2(aq)"
  swap_into_basis = "QuartzLike"
  constraint_species = "H2O              Na+              Cl-              QuartzLike"
  constraint_value = "  1.0E-2           0.1E-2           0.1E-2           1E-10"
  constraint_meaning = "kg_solvent_water bulk_composition bulk_composition free_mineral"
  constraint_unit = "   kg               moles            moles            moles"
  initial_temperature = 50.0
  mode = 4
  temperature = 200
  cold_temperature = 40.0
  source_species_names = 'H2O    Na+   Cl-   SiO2(aq)'
  source_species_rates = 'production_rate_H2O production_rate_Na production_rate_Cl production_rate_SiO2'
  ramp_max_ionic_strength_initial = 0 # max_ionic_strength in such a simple problem does not need ramping
  add_aux_pH = false # there is no H+ in this system
  evaluate_kinetic_rates_always = true # implicit time-marching used for stability
  execute_console_output_on = '' # only CSV output used in this example
[]

[UserObjects]
  [definition]
    type = GeochemicalModelDefinition
    database_file = "small_database.json"
    basis_species = "H2O SiO2(aq) Na+ Cl-"
    equilibrium_minerals = "QuartzLike"
  []
[]

[Executioner]
  type = Transient
  dt = 1E5
  end_time = 2E6 #7.76E6 # 90 days
[]

[AuxVariables]
  [production_temperature]
    initial_condition = 50 # the production_T Transfer lags one timestep behind for some reason, so give this a reasonable initial condition
  []
  [transported_H2O]
  []
  [transported_Na]
  []
  [transported_Cl]
  []
  [transported_SiO2]
  []
  [transported_mass]
  []
  [massfrac_H2O]
  []
  [massfrac_Na]
  []
  [massfrac_Cl]
  []
  [massfrac_SiO2]
  []
  [dumped_quartz]
  []
  [production_rate_H2O]
    initial_condition = 5.518533e+01 # the production_H2O Transfer lags one timestep behind for some reason (when the porous_flow simulation has finished, it correctly computes mole_rate_H2O_produced, but the Transfer gets the mole_rate_H2O_produced from the previous timestep), so give this a reasonable initial condition, otherwise this will be zero at the start of the simulation!
  []
  [production_rate_Na]
    initial_condition = 9.943302e-02
  []
  [production_rate_Cl]
    initial_condition = 9.943302e-02
  []
  [production_rate_SiO2]
    initial_condition = 2.340931e-04
  []
[]
[AuxKernels]
  [transported_H2O_auxk]
    type = GeochemistryQuantityAux
    variable = transported_H2O
    species = H2O
    quantity = transported_moles_in_original_basis
  []
  [transported_Na]
    type = GeochemistryQuantityAux
    variable = transported_Na
    species = Na+
    quantity = transported_moles_in_original_basis
  []
  [transported_Cl]
    type = GeochemistryQuantityAux
    variable = transported_Cl
    species = Cl-
    quantity = transported_moles_in_original_basis
  []
  [transported_SiO2]
    type = GeochemistryQuantityAux
    variable = transported_SiO2
    species = 'SiO2(aq)'
    quantity = transported_moles_in_original_basis
  []
  [transported_mass_auxk]
    type = ParsedAux
    args = 'transported_H2O transported_Na transported_Cl transported_SiO2'
    variable = transported_mass
    function = 'transported_H2O * 18.0152 + transported_Na * 22.9898 + transported_Cl * 35.453 + transported_SiO2 * 60.0843'
  []
  [massfrac_H2O]
    type = ParsedAux
    args = 'transported_mass transported_H2O'
    variable = massfrac_H2O
    function = '18.0152 * transported_H2O / transported_mass'
  []
  [massfrac_Na]
    type = ParsedAux
    args = 'transported_mass transported_Na'
    variable = massfrac_Na
    function = '22.9898 * transported_Na / transported_mass'
  []
  [massfrac_Cl]
    type = ParsedAux
    args = 'transported_mass transported_Cl'
    variable = massfrac_Cl
    function = '35.453 * transported_Cl / transported_mass'
  []
  [massfrac_SiO2]
    type = ParsedAux
    args = 'transported_mass transported_SiO2'
    variable = massfrac_SiO2
    function = '60.0843 * transported_SiO2 / transported_mass'
  []
  [dumped_quartz]
    type = GeochemistryQuantityAux
    variable = dumped_quartz
    species = QuartzLike
    quantity = moles_dumped
  []
[]
[Postprocessors]
  [cumulative_moles_precipitated_quartz]
    type = PointValue
    variable = dumped_quartz
  []
  [production_temperature]
    type = PointValue
    variable = production_temperature
  []
  [mass_heated_this_timestep]
    type = PointValue
    variable = transported_mass
  []
[]
[Outputs]
  csv = true
[]

[MultiApps]
  [porous_flow_sim]
    type = TransientMultiApp
    input_files = porous_flow.i
    execute_on = 'timestep_end'
  []
[]

[Transfers]
  [injection_T]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'solution_temperature'
    variable = 'injection_temperature'
  []
  [injection_Na]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Na'
    variable = 'injection_rate_massfrac_Na'
  []
  [injection_Cl]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Cl'
    variable = 'injection_rate_massfrac_Cl'
  []
  [injection_SiO2]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_SiO2'
    variable = 'injection_rate_massfrac_SiO2'
  []
  [injection_H2O]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_H2O'
    variable = 'injection_rate_massfrac_H2O'
  []

  [production_T]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = production_temperature
    variable = production_temperature
  []
  [production_Na]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Na_produced
    variable = production_rate_Na
  []
  [production_Cl]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Cl_produced
    variable = production_rate_Cl
  []
  [production_SiO2]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_SiO2_produced
    variable = production_rate_SiO2
  []
  [production_H2O]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_H2O_produced
    variable = production_rate_H2O
  []
[]

(tutorials/tutorial02_multiapps/step03_coupling/03_master_subcycling_picard.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [u]
  []
[]

[AuxVariables]
  [vt]
  []
[]

[Kernels]
  [diff]
    type = MatDiffusion
    variable = u
  []
  [force]
    type = BodyForce
    variable = u
    value = 1.
  []
  [td]
    type = TimeDerivative
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  []
[]

[Materials]
  [diff]
    type = ParsedMaterial
    f_name = D
    args = 'vt'
    function = 'vt'
  []
[]

[Executioner]
  type = Transient
  end_time = 2
  dt = 0.2

  picard_max_its = 10

  nl_abs_tol = 1e-10
  picard_rel_tol = 1e-6
  picard_abs_tol = 1e-10

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [micro]
    type = TransientMultiApp
    positions = '0.15 0.15 0  0.45 0.45 0  0.75 0.75 0'
    input_files = '03_sub_subcycling_picard.i'
    execute_on = timestep_end
    output_in_position = true
    sub_cycling = true
  []
[]

[Transfers]
  [push_u]
    type = MultiAppVariableValueSampleTransfer
    multi_app = micro
    direction = to_multiapp
    source_variable = u
    variable = ut
  []

  [pull_v]
    type = MultiAppPostprocessorInterpolationTransfer
    multi_app = micro
    direction = from_multiapp
    variable = vt
    postprocessor = average_v
  []
[]

(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/multilevel_sub.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./subsub_average]
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
  [./force]
    type = CoupledForce
    variable = u
    v = subsub_average
  [../]
[]

[BCs]
  [./left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  [../]
  [./right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  [../]
[]

[Postprocessors]
  [./sub_average]
    type = ElementAverageValue
    variable = u
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
  dt = 0.3

  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [./sub]
    type = TransientMultiApp
    app_type = MooseTestApp
    positions = '0 0 0 0.5 0.5 0'
    input_files = multilevel_subsub.i
  [../]
[]

[Transfers]
  [./subsub_average]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = from_multiapp
    multi_app = sub
    variable = subsub_average
    postprocessor = subsub_average
  [../]
[]

(modules/combined/examples/geochem-porous_flow/geotes_weber_tensleep/exchanger.i)

#########################################
#                                       #
# File written by create_input_files.py #
#                                       #
#########################################
# Model of the heat-exchanger
# The input fluid to the heat exchanger is determined by AuxVariables called production_temperature, production_rate_H, production_rate_Cl, production_rate_SO4, production_rate_HCO3, production_rate_SiO2aq, production_rate_Al, production_rate_Ca, production_rate_Mg, production_rate_Fe, production_rate_K, production_rate_Na, production_rate_Sr, production_rate_F, production_rate_BOH, production_rate_Br, production_rate_Ba, production_rate_Li, production_rate_NO3, production_rate_O2aq, production_rate_H2O.  These come from Postprocessors in the porous_flow.i simulation that measure the fluid composition at the production well.
# Given the input fluid, the exchanger cools/heats the fluid, removing any precipitates, and injects fluid back to porous_flow.i at temperature output_temperature and composition given by massfrac_H, etc.

[UserObjects]
  [definition]
    type = GeochemicalModelDefinition
    database_file = '../../../../geochemistry/database/moose_geochemdb.json'
    basis_species = 'H2O H+ Cl- SO4-- HCO3- SiO2(aq) Al+++ Ca++ Mg++ Fe++ K+ Na+ Sr++ F- B(OH)3 Br- Ba++ Li+ NO3- O2(aq)'
    equilibrium_minerals = 'Siderite Pyrrhotite Dolomite Illite Anhydrite Calcite Quartz K-feldspar Kaolinite Barite Celestite Fluorite Albite Chalcedony Goethite'
  []
[]

[TimeDependentReactionSolver]
  model_definition = definition
  include_moose_solve = false
  geochemistry_reactor_name = reactor
  swap_out_of_basis = 'NO3- O2(aq)'
  swap_into_basis = '  NH3  HS-'
  charge_balance_species = 'Cl-'
# initial conditions are unimportant because in exchanger mode all existing fluid is flushed from the system before adding the produced water
  constraint_species = 'H2O H+ Cl- SO4-- HCO3- SiO2(aq) Al+++ Ca++ Mg++ Fe++ K+ Na+ Sr++ F- B(OH)3 Br- Ba++ Li+ NH3 HS-'
  constraint_value = '1.0 1E-6 1E-6 1E-18 1E-18 1E-18    1E-18 1E-18 1E-18 1E-18 1E-18 1E-18 1E-18 1E-18 1E-18 1E-18 1E-18 1E-18 1E-18 1E-18'
  constraint_meaning = 'kg_solvent_water bulk_composition bulk_composition free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration free_concentration'
  constraint_unit = "kg moles moles molal molal molal molal molal molal molal molal molal molal molal molal molal molal molal molal molal"
  prevent_precipitation = 'Fluorite Albite Goethite'
  initial_temperature = 92
  mode = 4
  temperature = ramp_temperature # ramp up to 160degC over ~1 day so that aquifer geochemistry simulation can easily converge
  cold_temperature = 92
  heating_increments = 10
  source_species_names = ' H+ Cl- SO4-- HCO3- SiO2(aq) Al+++ Ca++ Mg++ Fe++ K+ Na+ Sr++ F- B(OH)3 Br- Ba++ Li+ NO3- O2(aq) H2O'
  source_species_rates = ' production_rate_H production_rate_Cl production_rate_SO4 production_rate_HCO3 production_rate_SiO2aq production_rate_Al production_rate_Ca production_rate_Mg production_rate_Fe production_rate_K production_rate_Na production_rate_Sr production_rate_F production_rate_BOH production_rate_Br production_rate_Ba production_rate_Li production_rate_NO3 production_rate_O2aq production_rate_H2O'
  ramp_max_ionic_strength_initial = 0 # max_ionic_strength in such a simple problem does not need ramping
[]

[GlobalParams]
  point = '0 0 0'
  reactor = reactor
[]

[AuxVariables]
  [ramp_temperature]
    initial_condition = 92
  []
  [production_temperature]
    initial_condition = 92 # the production_T Transfer lags one timestep behind for some reason, so give this a reasonable initial condition
  []
  [transported_H]
  []
  [transported_Cl]
  []
  [transported_SO4]
  []
  [transported_HCO3]
  []
  [transported_SiO2aq]
  []
  [transported_Al]
  []
  [transported_Ca]
  []
  [transported_Mg]
  []
  [transported_Fe]
  []
  [transported_K]
  []
  [transported_Na]
  []
  [transported_Sr]
  []
  [transported_F]
  []
  [transported_BOH]
  []
  [transported_Br]
  []
  [transported_Ba]
  []
  [transported_Li]
  []
  [transported_NO3]
  []
  [transported_O2aq]
  []
  [transported_H2O]
  []
  [transported_mass]
  []
  [massfrac_H]
  []
  [massfrac_Cl]
  []
  [massfrac_SO4]
  []
  [massfrac_HCO3]
  []
  [massfrac_SiO2aq]
  []
  [massfrac_Al]
  []
  [massfrac_Ca]
  []
  [massfrac_Mg]
  []
  [massfrac_Fe]
  []
  [massfrac_K]
  []
  [massfrac_Na]
  []
  [massfrac_Sr]
  []
  [massfrac_F]
  []
  [massfrac_BOH]
  []
  [massfrac_Br]
  []
  [massfrac_Ba]
  []
  [massfrac_Li]
  []
  [massfrac_NO3]
  []
  [massfrac_O2aq]
  []
  [massfrac_H2O]
  []
  [dumped_Siderite]
  []
  [dumped_Pyrrhotite]
  []
  [dumped_Dolomite]
  []
  [dumped_Illite]
  []
  [dumped_Anhydrite]
  []
  [dumped_Calcite]
  []
  [dumped_Quartz]
  []
  [dumped_K-feldspar]
  []
  [dumped_Kaolinite]
  []
  [dumped_Barite]
  []
  [dumped_Celestite]
  []
  [dumped_Fluorite]
  []
  [dumped_Albite]
  []
  [dumped_Chalcedony]
  []
  [dumped_Goethite]
  []
# The production_* Transfers lag one timestep behind for some reason (when the porous_flow simulation has finished, it correctly computes mole_rate_*_produced, but the Transfer gets the mole_rate_*_produced from the previous timestep), so give the production_rate_* reasonable initial conditions, otherwise they will be zero at the start of the simulation.
  [production_rate_H]
    initial_condition = -0.00058596786807342
  []
  [production_rate_Cl]
    initial_condition = 0.274767413291287
  []
  [production_rate_SO4]
    initial_condition = 0.012567456786868922
  []
  [production_rate_HCO3]
    initial_condition = 0.0001668295857850308
  []
  [production_rate_SiO2aq]
    initial_condition = 0.00010068057668449495
  []
  [production_rate_Al]
    initial_condition = 2.4224219572143877e-07
  []
  [production_rate_Ca]
    initial_condition = 0.0040997718654983036
  []
  [production_rate_Mg]
    initial_condition = 0.00015261342984691217
  []
  [production_rate_Fe]
    initial_condition = 0.0001550375425863269
  []
  [production_rate_K]
    initial_condition = 0.0003500651859998926
  []
  [production_rate_Na]
    initial_condition = 0.2896767602995328
  []
  [production_rate_Sr]
    initial_condition = 2.915285700108879e-05
  []
  [production_rate_F]
    initial_condition = 5.8582680830041476e-05
  []
  [production_rate_BOH]
    initial_condition = 0.0012157199878760335
  []
  [production_rate_Br]
    initial_condition = 0.00022605948665165203
  []
  [production_rate_Ba]
    initial_condition = 2.2773554030672105e-07
  []
  [production_rate_Li]
    initial_condition = 0.0023920780265869763
  []
  [production_rate_NO3]
    initial_condition = 0.000353470613973057
  []
  [production_rate_O2aq]
    initial_condition = -0.00044255942331181803
  []
  [production_rate_H2O]
    initial_condition = 10.10458252764496
  []
[]

[AuxKernels]
  [ramp_temperature]
    type = FunctionAux
    variable = ramp_temperature
    function = 'min(160, max(92, 92 + (160 - 92) * t / 1E5))'
  []
  [transported_H_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_H
    species = 'H+'
  []
  [transported_Cl_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Cl
    species = 'Cl-'
  []
  [transported_SO4_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_SO4
    species = 'SO4--'
  []
  [transported_HCO3_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_HCO3
    species = 'HCO3-'
  []
  [transported_SiO2aq_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_SiO2aq
    species = 'SiO2(aq)'
  []
  [transported_Al_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Al
    species = 'Al+++'
  []
  [transported_Ca_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Ca
    species = 'Ca++'
  []
  [transported_Mg_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Mg
    species = 'Mg++'
  []
  [transported_Fe_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Fe
    species = 'Fe++'
  []
  [transported_K_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_K
    species = 'K+'
  []
  [transported_Na_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Na
    species = 'Na+'
  []
  [transported_Sr_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Sr
    species = 'Sr++'
  []
  [transported_F_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_F
    species = 'F-'
  []
  [transported_BOH_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_BOH
    species = 'B(OH)3'
  []
  [transported_Br_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Br
    species = 'Br-'
  []
  [transported_Ba_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Ba
    species = 'Ba++'
  []
  [transported_Li_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_Li
    species = 'Li+'
  []
  [transported_NO3_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_NO3
    species = 'NO3-'
  []
  [transported_O2aq_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_O2aq
    species = 'O2(aq)'
  []
  [transported_H2O_auxk]
    type = GeochemistryQuantityAux
    quantity = transported_moles_in_original_basis
    variable = transported_H2O
    species = 'H2O'
  []
  [transported_mass_auxk]
    type = ParsedAux
    args = ' transported_H transported_Cl transported_SO4 transported_HCO3 transported_SiO2aq transported_Al transported_Ca transported_Mg transported_Fe transported_K transported_Na transported_Sr transported_F transported_BOH transported_Br transported_Ba transported_Li transported_NO3 transported_O2aq transported_H2O'
    variable = transported_mass
    function = ' transported_H * 1.0079 + transported_Cl * 35.453 + transported_SO4 * 96.0576 + transported_HCO3 * 61.0171 + transported_SiO2aq * 60.0843 + transported_Al * 26.9815 + transported_Ca * 40.08 + transported_Mg * 24.305 + transported_Fe * 55.847 + transported_K * 39.0983 + transported_Na * 22.9898 + transported_Sr * 87.62 + transported_F * 18.9984 + transported_BOH * 61.8329 + transported_Br * 79.904 + transported_Ba * 137.33 + transported_Li * 6.941 + transported_NO3 * 62.0049 + transported_O2aq * 31.9988 + transported_H2O * 18.01801802'
  []
  [massfrac_H_auxk]
    type = ParsedAux
    args = 'transported_mass transported_H'
    variable = massfrac_H
    function = '1.0079 * transported_H / transported_mass'
  []
  [massfrac_Cl_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Cl'
    variable = massfrac_Cl
    function = '35.453 * transported_Cl / transported_mass'
  []
  [massfrac_SO4_auxk]
    type = ParsedAux
    args = 'transported_mass transported_SO4'
    variable = massfrac_SO4
    function = '96.0576 * transported_SO4 / transported_mass'
  []
  [massfrac_HCO3_auxk]
    type = ParsedAux
    args = 'transported_mass transported_HCO3'
    variable = massfrac_HCO3
    function = '61.0171 * transported_HCO3 / transported_mass'
  []
  [massfrac_SiO2aq_auxk]
    type = ParsedAux
    args = 'transported_mass transported_SiO2aq'
    variable = massfrac_SiO2aq
    function = '60.0843 * transported_SiO2aq / transported_mass'
  []
  [massfrac_Al_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Al'
    variable = massfrac_Al
    function = '26.9815 * transported_Al / transported_mass'
  []
  [massfrac_Ca_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Ca'
    variable = massfrac_Ca
    function = '40.08 * transported_Ca / transported_mass'
  []
  [massfrac_Mg_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Mg'
    variable = massfrac_Mg
    function = '24.305 * transported_Mg / transported_mass'
  []
  [massfrac_Fe_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Fe'
    variable = massfrac_Fe
    function = '55.847 * transported_Fe / transported_mass'
  []
  [massfrac_K_auxk]
    type = ParsedAux
    args = 'transported_mass transported_K'
    variable = massfrac_K
    function = '39.0983 * transported_K / transported_mass'
  []
  [massfrac_Na_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Na'
    variable = massfrac_Na
    function = '22.9898 * transported_Na / transported_mass'
  []
  [massfrac_Sr_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Sr'
    variable = massfrac_Sr
    function = '87.62 * transported_Sr / transported_mass'
  []
  [massfrac_F_auxk]
    type = ParsedAux
    args = 'transported_mass transported_F'
    variable = massfrac_F
    function = '18.9984 * transported_F / transported_mass'
  []
  [massfrac_BOH_auxk]
    type = ParsedAux
    args = 'transported_mass transported_BOH'
    variable = massfrac_BOH
    function = '61.8329 * transported_BOH / transported_mass'
  []
  [massfrac_Br_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Br'
    variable = massfrac_Br
    function = '79.904 * transported_Br / transported_mass'
  []
  [massfrac_Ba_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Ba'
    variable = massfrac_Ba
    function = '137.33 * transported_Ba / transported_mass'
  []
  [massfrac_Li_auxk]
    type = ParsedAux
    args = 'transported_mass transported_Li'
    variable = massfrac_Li
    function = '6.941 * transported_Li / transported_mass'
  []
  [massfrac_NO3_auxk]
    type = ParsedAux
    args = 'transported_mass transported_NO3'
    variable = massfrac_NO3
    function = '62.0049 * transported_NO3 / transported_mass'
  []
  [massfrac_O2aq_auxk]
    type = ParsedAux
    args = 'transported_mass transported_O2aq'
    variable = massfrac_O2aq
    function = '31.9988 * transported_O2aq / transported_mass'
  []
  [massfrac_H2O_auxk]
    type = ParsedAux
    args = 'transported_mass transported_H2O'
    variable = massfrac_H2O
    function = '18.01801802 * transported_H2O / transported_mass'
  []
  [dumped_Siderite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Siderite
    species = Siderite
    quantity = moles_dumped
  []
  [dumped_Pyrrhotite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Pyrrhotite
    species = Pyrrhotite
    quantity = moles_dumped
  []
  [dumped_Dolomite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Dolomite
    species = Dolomite
    quantity = moles_dumped
  []
  [dumped_Illite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Illite
    species = Illite
    quantity = moles_dumped
  []
  [dumped_Anhydrite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Anhydrite
    species = Anhydrite
    quantity = moles_dumped
  []
  [dumped_Calcite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Calcite
    species = Calcite
    quantity = moles_dumped
  []
  [dumped_Quartz_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Quartz
    species = Quartz
    quantity = moles_dumped
  []
  [dumped_K-feldspar_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_K-feldspar
    species = K-feldspar
    quantity = moles_dumped
  []
  [dumped_Kaolinite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Kaolinite
    species = Kaolinite
    quantity = moles_dumped
  []
  [dumped_Barite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Barite
    species = Barite
    quantity = moles_dumped
  []
  [dumped_Celestite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Celestite
    species = Celestite
    quantity = moles_dumped
  []
  [dumped_Fluorite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Fluorite
    species = Fluorite
    quantity = moles_dumped
  []
  [dumped_Albite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Albite
    species = Albite
    quantity = moles_dumped
  []
  [dumped_Chalcedony_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Chalcedony
    species = Chalcedony
    quantity = moles_dumped
  []
  [dumped_Goethite_auxk]
    type = GeochemistryQuantityAux
    variable = dumped_Goethite
    species = Goethite
    quantity = moles_dumped
  []
[]

[Postprocessors]
  [cumulative_moles_precipitated_Siderite]
    type = PointValue
    variable = dumped_Siderite
  []
  [cumulative_moles_precipitated_Pyrrhotite]
    type = PointValue
    variable = dumped_Pyrrhotite
  []
  [cumulative_moles_precipitated_Dolomite]
    type = PointValue
    variable = dumped_Dolomite
  []
  [cumulative_moles_precipitated_Illite]
    type = PointValue
    variable = dumped_Illite
  []
  [cumulative_moles_precipitated_Anhydrite]
    type = PointValue
    variable = dumped_Anhydrite
  []
  [cumulative_moles_precipitated_Calcite]
    type = PointValue
    variable = dumped_Calcite
  []
  [cumulative_moles_precipitated_Quartz]
    type = PointValue
    variable = dumped_Quartz
  []
  [cumulative_moles_precipitated_K-feldspar]
    type = PointValue
    variable = dumped_K-feldspar
  []
  [cumulative_moles_precipitated_Kaolinite]
    type = PointValue
    variable = dumped_Kaolinite
  []
  [cumulative_moles_precipitated_Barite]
    type = PointValue
    variable = dumped_Barite
  []
  [cumulative_moles_precipitated_Celestite]
    type = PointValue
    variable = dumped_Celestite
  []
  [cumulative_moles_precipitated_Fluorite]
    type = PointValue
    variable = dumped_Fluorite
  []
  [cumulative_moles_precipitated_Albite]
    type = PointValue
    variable = dumped_Albite
  []
  [cumulative_moles_precipitated_Chalcedony]
    type = PointValue
    variable = dumped_Chalcedony
  []
  [cumulative_moles_precipitated_Goethite]
    type = PointValue
    variable = dumped_Goethite
  []
  [production_temperature]
    type = PointValue
    variable = production_temperature
  []
  [mass_heated_this_timestep]
    type = PointValue
    variable = transported_mass
  []
[]

[Outputs]
  csv = true
[]

[Executioner]
  type = Transient
  solve_type = Newton
  end_time = 7.76E6 # 90 days
  [TimeStepper]
    type = FunctionDT
    function = 'min(3E4, max(1E4, 0.2 * t))'
  []
[]

[MultiApps]
  [porous_flow_sim]
    type = TransientMultiApp
    input_files = porous_flow.i
    execute_on = 'timestep_end'
  []
[]

[Transfers]
  [injection_T]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'solution_temperature'
    variable = 'injection_temperature'
  []
  [injection_H]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_H'
    variable = 'injection_rate_massfrac_H'
  []
  [injection_Cl]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Cl'
    variable = 'injection_rate_massfrac_Cl'
  []
  [injection_SO4]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_SO4'
    variable = 'injection_rate_massfrac_SO4'
  []
  [injection_HCO3]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_HCO3'
    variable = 'injection_rate_massfrac_HCO3'
  []
  [injection_SiO2aq]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_SiO2aq'
    variable = 'injection_rate_massfrac_SiO2aq'
  []
  [injection_Al]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Al'
    variable = 'injection_rate_massfrac_Al'
  []
  [injection_Ca]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Ca'
    variable = 'injection_rate_massfrac_Ca'
  []
  [injection_Mg]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Mg'
    variable = 'injection_rate_massfrac_Mg'
  []
  [injection_Fe]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Fe'
    variable = 'injection_rate_massfrac_Fe'
  []
  [injection_K]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_K'
    variable = 'injection_rate_massfrac_K'
  []
  [injection_Na]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Na'
    variable = 'injection_rate_massfrac_Na'
  []
  [injection_Sr]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Sr'
    variable = 'injection_rate_massfrac_Sr'
  []
  [injection_F]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_F'
    variable = 'injection_rate_massfrac_F'
  []
  [injection_BOH]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_BOH'
    variable = 'injection_rate_massfrac_BOH'
  []
  [injection_Br]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Br'
    variable = 'injection_rate_massfrac_Br'
  []
  [injection_Ba]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Ba'
    variable = 'injection_rate_massfrac_Ba'
  []
  [injection_Li]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_Li'
    variable = 'injection_rate_massfrac_Li'
  []
  [injection_NO3]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_NO3'
    variable = 'injection_rate_massfrac_NO3'
  []
  [injection_O2aq]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_O2aq'
    variable = 'injection_rate_massfrac_O2aq'
  []
  [injection_H2O]
    type = MultiAppNearestNodeTransfer
    direction = TO_MULTIAPP
    multi_app = porous_flow_sim
    fixed_meshes = true
    source_variable = 'massfrac_H2O'
    variable = 'injection_rate_massfrac_H2O'
  []

  [production_T]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = production_temperature
    variable = production_temperature
  []
  [production_H]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_H_produced
    variable = production_rate_H
  []
  [production_Cl]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Cl_produced
    variable = production_rate_Cl
  []
  [production_SO4]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_SO4_produced
    variable = production_rate_SO4
  []
  [production_HCO3]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_HCO3_produced
    variable = production_rate_HCO3
  []
  [production_SiO2aq]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_SiO2aq_produced
    variable = production_rate_SiO2aq
  []
  [production_Al]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Al_produced
    variable = production_rate_Al
  []
  [production_Ca]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Ca_produced
    variable = production_rate_Ca
  []
  [production_Mg]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Mg_produced
    variable = production_rate_Mg
  []
  [production_Fe]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Fe_produced
    variable = production_rate_Fe
  []
  [production_K]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_K_produced
    variable = production_rate_K
  []
  [production_Na]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Na_produced
    variable = production_rate_Na
  []
  [production_Sr]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Sr_produced
    variable = production_rate_Sr
  []
  [production_F]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_F_produced
    variable = production_rate_F
  []
  [production_BOH]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_BOH_produced
    variable = production_rate_BOH
  []
  [production_Br]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Br_produced
    variable = production_rate_Br
  []
  [production_Ba]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Ba_produced
    variable = production_rate_Ba
  []
  [production_Li]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_Li_produced
    variable = production_rate_Li
  []
  [production_NO3]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_NO3_produced
    variable = production_rate_NO3
  []
  [production_O2aq]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_O2aq_produced
    variable = production_rate_O2aq
  []
  [production_H2O]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = FROM_MULTIAPP
    multi_app = porous_flow_sim
    postprocessor = mole_rate_H2O_produced
    variable = production_rate_H2O
  []
[]

(test/tests/transfers/multiapp_postprocessor_interpolation_transfer/master_quad.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./pp_aux]
  [../]
[]

[Kernels]
  [./diff]
    type = CoefDiffusion
    variable = u
    coef = 0.1
  [../]
  [./time]
    type = TimeDerivative
    variable = u
  [../]
[]

[BCs]
  [./left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  [../]
  [./right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 20
  dt = 0.1
  solve_type = PJFNK
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

[MultiApps]
  [./quad]
    type = TransientMultiApp
    app_type = MooseTestApp
    positions = '0.1 0.1 0 0.9 0.1 0 0.1 0.9 0 0.9 0.9 0'
    input_files = 'quad_sub1.i quad_sub1.i quad_sub2.i quad_sub2.i'
  [../]
[]

[Transfers]
  [./sub_to_master_pp]
    type = MultiAppPostprocessorInterpolationTransfer
    direction = from_multiapp
    multi_app = quad
    variable = pp_aux
    postprocessor = pp
  [../]
[]

Overview
Example Input File Syntax
Input Parameters
Input Files

Install MOOSE

New Users

Examples and Tutorials

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools

INL Applications and Remote Access