Fixes PrecursorAction with velocity functions, and supports higher order precursor variable initialization

include/actions/PrecursorAction.h

@@ -51,6 +51,13 @@ class PrecursorAction : public VariableNotAMooseObjectAction
   */
   void addTimeDerivative(const std::string & var_name);
 
+  /**
+  * Adds appropriate conservative advection kernel
+  *
+  * @param var_name The name of the variable the kernel acts on
+  */
+  void addAdvection(const std::string & var_name);
+
   /**
   * Adds DGAdvection kernel
   *
@@ -72,6 +79,11 @@ class PrecursorAction : public VariableNotAMooseObjectAction
   */
   void addInflowBC(const std::string & var_name);
 
+  /**
+   * Adds a PostprocessorPenaltyDirichletBC to stabilize precursor conc at inlet
+   */
+  void addPenaltyBC(const std::string & var_name);
+
   /**
   * Adds random initial conditions for Jacobian testing
   *
@@ -127,5 +139,9 @@ class PrecursorAction : public VariableNotAMooseObjectAction
   /// optional object name suffix
   std::string _object_suffix;
 
+  /// whether input file is for the outer loop
   bool _is_loopapp;
+
+  /// velocity type
+  MooseEnum _velocity_type;
 };

include/postprocessors/SideFunctionWeightedIntegralPostprocessor.h

@@ -0,0 +1,21 @@
+#pragma once
+
+#include "SideIntegralVariablePostprocessor.h"
+
+/**
+ * This postprocessor computes the weighted integral of a
+ * variable along a boundary, weighted by the user-provided function.
+ */
+class SideFunctionWeightedIntegralPostprocessor : public SideIntegralVariablePostprocessor
+{
+public:
+  SideFunctionWeightedIntegralPostprocessor(const InputParameters & parameters);
+
+  static InputParameters validParams();
+
+protected:
+  virtual Real computeQpIntegral() override;
+
+  // Weight function
+  const Function & _weight_func;
+};

problems/2021-cnrs-benchmark/phase-0/nts.i

@@ -44,7 +44,7 @@
   [pres]
     var_name_base = pre
     outlet_boundaries = ''
-    constant_velocity_values = true
+    velocity_type = constant
     u_def = 0
     v_def = 0
     w_def = 0

problems/2021-cnrs-benchmark/phase-0/tests

@@ -1,7 +1,7 @@
 [Tests]
   [vel_field_syntax]
     type = RunApp
-    input = nts.i
+    input = vel-field.i
     check_input = True
     method = opt
   []
@@ -13,7 +13,7 @@
   []
   [temperature_syntax]
     type = RunApp
-    input = nts.i
+    input = temperature.i
     check_input = True
     method = opt
   []

problems/2021-cnrs-benchmark/phase-1/buoyancy-nts.i

@@ -42,7 +42,7 @@
   [pres]
     var_name_base = pre
     outlet_boundaries = ''
-    constant_velocity_values = false
+    velocity_type = variable
     uvel = vel_x
     vvel = vel_y
     nt_exp_form = false

problems/2021-cnrs-benchmark/phase-1/circ-fuel.i

@@ -45,7 +45,7 @@
   [pres]
     var_name_base = pre
     outlet_boundaries = ''
-    constant_velocity_values = false
+    velocity_type = variable
     uvel = ux
     vvel = uy
     nt_exp_form = false

problems/2021-cnrs-benchmark/phase-1/full-coupling-nts.i

@@ -42,7 +42,7 @@
   [pres]
     var_name_base = pre
     outlet_boundaries = ''
-    constant_velocity_values = false
+    velocity_type = variable
     uvel = vel_x
     vvel = vel_y
     nt_exp_form = false

problems/2021-cnrs-benchmark/phase-1/power-coupling-nts.i

@@ -42,7 +42,7 @@
   [pres]
     var_name_base = pre
     outlet_boundaries = ''
-    constant_velocity_values = false
+    velocity_type = variable
     uvel = ux
     vvel = uy
     nt_exp_form = false

problems/2021-cnrs-benchmark/phase-2/transient.i

@@ -41,6 +41,7 @@ dt = 0.00625        # Timestep size = 1 / freq / 200 = 0.00625 s
 
 [Problem]
   type = FEProblem
+  allow_initial_conditions_with_restart = true
 []
 
 [Variables]
@@ -69,7 +70,7 @@ dt = 0.00625        # Timestep size = 1 / freq / 200 = 0.00625 s
   [pres]
     var_name_base = pre
     outlet_boundaries = ''
-    constant_velocity_values = false
+    velocity_type = variable
     uvel = velx
     vvel = vely
     nt_exp_form = false
@@ -214,7 +215,7 @@ dt = 0.00625        # Timestep size = 1 / freq / 200 = 0.00625 s
 [Functions]
   [func_alpha]
     type = ParsedFunction
-    value = '1 + 0.1 * sin(2*pi*t*${freq})' # Perturbation frequency = 0.8Hz
+    expression = '1 + 0.1 * sin(2*pi*t*${freq})' # Perturbation frequency = 0.8Hz
   []
   [velxf]
     type = SolutionFunction