Merge pull request #28060 from GiudGiud/PR_rcnotime2

Add option to not contribute to RC coefficients in the Physics
idaholab · Jul 4, 2024 · 10c4f77 · 10c4f77
2 parents c32a23e + 0c0df74
commit 10c4f77
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/modules/navier_stokes/src/physics/WCNSFVFlowPhysics.C b/modules/navier_stokes/src/physics/WCNSFVFlowPhysics.C
@@ -12,6 +12,7 @@
 #include "NSFVBase.h"
 #include "INSFVMomentumAdvection.h"
 #include "INSFVRhieChowInterpolator.h"
+#include "INSFVTimeKernel.h"
 #include "MapConversionUtils.h"
 #include "NS.h"
 
@@ -47,6 +48,10 @@ WCNSFVFlowPhysics::validParams()
 
   params += NSFVBase::commonNavierStokesFlowParams();
 
+  // Time derivative term correction
+  params.transferParam<bool>(INSFVTimeKernel::validParams(), "contribute_to_rc");
+  params.addParamNamesToGroup("contribute_to_rc", "Advanced");
+
   // Used for flow mixtures, where one phase is solid / not moving under the action of gravity
   params.addParam<MooseFunctorName>(
       "density_for_gravity_terms",
@@ -529,6 +534,7 @@ WCNSFVFlowPhysics::addINSMomentumTimeKernels()
   assignBlocks(params, _blocks);
   params.set<MooseFunctorName>(NS::density) = _density_name;
   params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
+  params.set<bool>("contribute_to_rc") = getParam<bool>("contribute_to_rc");
 
   for (const auto d : make_range(dimension()))
   {
@@ -547,6 +553,7 @@ WCNSFVFlowPhysics::addWCNSMomentumTimeKernels()
   assignBlocks(params, _blocks);
   params.set<MooseFunctorName>(NS::density) = _density_name;
   params.set<MooseFunctorName>(NS::time_deriv(NS::density)) = NS::time_deriv(_density_name);
+  params.set<bool>("contribute_to_rc") = getParam<bool>("contribute_to_rc");
 
   for (const auto d : make_range(dimension()))
   {