Added passives support to InjectionProductionPhysics

Closes #31847

Author: joshuahansel

modules/porous_flow/doc/content/source/physics/CoupledInjectionProductionPhysics.md

@@ -1,19 +1,23 @@
 # CoupledInjectionProduction
 
-This [Physics](syntax/Physics/index.md) object adds dirac kernels for mass and energy
+This [Physics](syntax/Physics/index.md) object adds dirac kernels for mass, energy, and passive concentration
 sources at injection and production points within a porous matrix, where the sources
 are provided by post-processor values transferred from another application.
 
 This object creates the following post-processors:
 
 - `p_inj1`, ..., `p_inj<n_inj>`: Pressure at each injection point, via [PointValue.md].
 - `T_inj1`, ..., `T_inj<n_inj>`: Temperature at each injection point, via [PointValue.md].
+- `<passive>_inj1`, ..., `<passive>_inj<n_inj>`: Passive concentration `<passive>` at each injection point, via [PointValue.md].
 - `p_pro1`, ..., `p_pro<n_pro>`: Pressure at each production point, via [PointValue.md].
 - `T_pro1`, ..., `T_pro<n_pro>`: Temperature at each production point, via [PointValue.md].
+- `<passive>_pro1`, ..., `<passive>_pro<n_pro>`: Passive concentration `<passive>` at each production point, via [PointValue.md].
 - `mass_rate_inj1`, ..., `mass_rate_inj<n_inj>`: Mass source rate at each injection point, to be transferred from the other application using [Receiver.md].
 - `energy_rate_inj1`, ..., `energy_rate_inj<n_inj>`: Energy source rate at each injection point, to be transferred from the other application using [Receiver.md].
+- `<passive>_rate_inj1`, ..., `<passive>_rate_inj<n_inj>`: Passive concentration `<passive>` source rate at each injection point, to be transferred from the other application using [Receiver.md].
 - `mass_rate_pro1`, ..., `mass_rate_pro<n_pro>`: Mass source rate at each production point, to be transferred from the other application using [Receiver.md].
 - `energy_rate_pro1`, ..., `energy_rate_pro<n_pro>`: Energy source rate at each production point, to be transferred from the other application using [Receiver.md].
+- `<passive>_rate_pro1`, ..., `<passive>_rate_pro<n_pro>`: Passive concentration `<passive>` source rate at each production point, to be transferred from the other application using [Receiver.md].
 
 `n_inj` is the size of the [!param](/Physics/CoupledInjectionProduction/CoupledInjectionProductionPhysics/injection_points) parameter, and
 `n_pro` is the size of the [!param](/Physics/CoupledInjectionProduction/CoupledInjectionProductionPhysics/production_points) parameter.
@@ -24,7 +28,7 @@ objects are created by this object to transfer the post-processors listed above
 if the parameter is not provided, then the other application is responsible for creating the transfers.
 
 [PorousFlowPointSourceFromPostprocessor.md] dirac kernels are added for each injection and production point,
-using the `mass_rate_*` and `energy_rate_*` post-processors listed above.
+using the `mass_rate_*`, `energy_rate_*`, and `<passive>_rate_*` post-processors listed above.
 
 !syntax parameters /Physics/CoupledInjectionProduction/CoupledInjectionProductionPhysics
 

modules/porous_flow/include/physics/CoupledInjectionProductionPhysics.h

@@ -40,6 +40,9 @@ class CoupledInjectionProductionPhysics : public PhysicsBase
   /// Production points
   const std::vector<Point> & _production_points;
 
+  /// Mass fraction variable names, if any
+  const std::vector<VariableName> & _mass_fraction_vars;
+
   /// Injection and production points
   std::vector<Point> _points;
   /// Label for each point

modules/porous_flow/src/physics/CoupledInjectionProductionPhysics.C

@@ -21,6 +21,8 @@ CoupledInjectionProductionPhysics::validParams()
 
   params.addRequiredParam<std::vector<Point>>("injection_points", "List of injection points [m]");
   params.addRequiredParam<std::vector<Point>>("production_points", "List of production points [m]");
+  params.addRequiredParam<std::vector<VariableName>>("mass_fraction_vars",
+                                                     "List of mass fraction variables");
   params.addParam<MultiAppName>("multi_app", "MultiApp to transfer to and from");
 
   params.addClassDescription(
@@ -33,7 +35,8 @@ CoupledInjectionProductionPhysics::CoupledInjectionProductionPhysics(
     const InputParameters & parameters)
   : PhysicsBase(parameters),
     _injection_points(getParam<std::vector<Point>>("injection_points")),
-    _production_points(getParam<std::vector<Point>>("production_points"))
+    _production_points(getParam<std::vector<Point>>("production_points")),
+    _mass_fraction_vars(getParam<std::vector<VariableName>>("mass_fraction_vars"))
 {
   _points = _injection_points;
   _points.insert(_points.end(), _production_points.begin(), _production_points.end());
@@ -51,6 +54,8 @@ CoupledInjectionProductionPhysics::addDiracKernels()
   {
     addPPSourceDiracKernel(_points[i], "porepressure", "mass_rate_" + _labels[i]);
     addPPSourceDiracKernel(_points[i], "temperature", "energy_rate_" + _labels[i]);
+    for (const auto & var : _mass_fraction_vars)
+      addPPSourceDiracKernel(_points[i], var, var + "_rate_" + _labels[i]);
   }
 }
 
@@ -63,6 +68,11 @@ CoupledInjectionProductionPhysics::addPostprocessors()
     addPointValuePostprocessor("temperature", _points[i], "T_" + _labels[i]);
     addReceiverPostprocessor("mass_rate_" + _labels[i]);
     addReceiverPostprocessor("energy_rate_" + _labels[i]);
+    for (const auto & var : _mass_fraction_vars)
+    {
+      addPointValuePostprocessor(var, _points[i], var + "_" + _labels[i]);
+      addReceiverPostprocessor(var + "_rate_" + _labels[i]);
+    }
   }
 }
 
@@ -76,6 +86,11 @@ CoupledInjectionProductionPhysics::addTransfers()
       addPostprocessorTransfer("T_" + _labels[i], false);
       addPostprocessorTransfer("mass_rate_" + _labels[i], true);
       addPostprocessorTransfer("energy_rate_" + _labels[i], true);
+      for (const auto & var : _mass_fraction_vars)
+      {
+        addPostprocessorTransfer(var + "_" + _labels[i], false);
+        addPostprocessorTransfer(var + "_rate_" + _labels[i], true);
+      }
     }
 }
 

modules/porous_flow/test/tests/physics/coupled_injection_production_physics/gold/main_out.e

@@ -1,6 +1,7 @@
 gravity_vector = '0 0 0'
 initial_p = 1e5
 initial_T = 300
+initial_C1 = 1e-3
 
 porosity = 0.8
 perm = 1e-8 # permeability
@@ -34,6 +35,7 @@ pro_point1 = '3.5 1.5 1.5'
   fp = fp
   pressure_unit = Pa
   stabilization = full
+  mass_fraction_vars = 'C1'
 []
 
 [FluidProperties]
@@ -50,6 +52,9 @@ pro_point1 = '3.5 1.5 1.5'
     initial_condition = ${initial_T}
     scaling = 1e-6
   []
+  [C1]
+    initial_condition = ${initial_C1}
+  []
 []
 
 [Materials]
@@ -104,11 +109,13 @@ pro_point1 = '3.5 1.5 1.5'
     [inj_prod]
       injection_points = '${inj_point1}'
       production_points = '${pro_point1}'
+      mass_fraction_vars = 'C1'
       multi_app = sub
     []
   []
 []
 
 [Outputs]
   exodus = true
+  show = 'porepressure temperature C1'
 []

modules/porous_flow/test/tests/physics/coupled_injection_production_physics/main.i

@@ -1,8 +1,11 @@
 e_source = 1e5
 mass_rate_inj1 = 0.5
+C1_source = 2e-3
 energy_rate_inj1 = ${fparse mass_rate_inj1 * e_source}
+C1_rate_inj1 = ${fparse mass_rate_inj1 * C1_source}
 mass_rate_pro1 = -0.5
 energy_rate_pro1 = ${fparse mass_rate_pro1 * e_source}
+C1_rate_pro1 = ${fparse mass_rate_pro1 * C1_source}
 
 [Mesh]
   type = GeneratedMesh
@@ -19,12 +22,19 @@ energy_rate_pro1 = ${fparse mass_rate_pro1 * e_source}
     type = ConstantPostprocessor
     value = ${energy_rate_inj1}
   []
+  [C1_rate_inj1]
+    type = ConstantPostprocessor
+    value = ${C1_rate_inj1}
+  []
   [p_inj1]
     type = Receiver
   []
   [T_inj1]
     type = Receiver
   []
+  [C1_inj1]
+    type = Receiver
+  []
 
   [mass_rate_pro1]
     type = ConstantPostprocessor
@@ -34,12 +44,19 @@ energy_rate_pro1 = ${fparse mass_rate_pro1 * e_source}
     type = ConstantPostprocessor
     value = ${energy_rate_pro1}
   []
+  [C1_rate_pro1]
+    type = ConstantPostprocessor
+    value = ${C1_rate_pro1}
+  []
   [p_pro1]
     type = Receiver
   []
   [T_pro1]
     type = Receiver
   []
+  [C1_pro1]
+    type = Receiver
+  []
 []
 
 [Problem]