diff --git a/framework/include/base/MooseFunctor.h b/framework/include/base/MooseFunctor.h
index 65d9e7d1746a..c9759abd1e75 100644
--- a/framework/include/base/MooseFunctor.h
+++ b/framework/include/base/MooseFunctor.h
@@ -205,6 +205,25 @@ struct SingleSidedFaceArg
   /// The subdomain ID which denotes the side of the face information we are evaluating on
   SubdomainID sub_id;
 
+  /**
+   * Make an element argument corresponding to the subdomain ID on which we are defined
+   */
+  ElemArg makeSidedElem() const
+  {
+    mooseAssert(fi, "The face must be non-null");
+#ifndef NDEBUG
+    const unsigned short matches = (sub_id == fi->elem().subdomain_id()) +
+                                   (fi->neighborPtr() && (sub_id == fi->neighbor().subdomain_id()));
+    mooseAssert(matches == 1,
+                "We should only have on match. If we have more or less, then this is not an "
+                "appropriate use of SingleSidedFaceArg");
+#endif
+
+    const Elem * const ret_elem =
+        sub_id == fi->elem().subdomain_id() ? &fi->elem() : fi->neighborPtr();
+    return {ret_elem, correct_skewness, apply_gradient_to_skewness};
+  }
+
   /**
    * Make a \p ElemArg from our data using the face information element
    */
diff --git a/framework/include/fvkernels/FVFluxKernel.h b/framework/include/fvkernels/FVFluxKernel.h
index 69dc0f97b183..cb7d40df7356 100644
--- a/framework/include/fvkernels/FVFluxKernel.h
+++ b/framework/include/fvkernels/FVFluxKernel.h
@@ -110,6 +110,20 @@ class FVFluxKernel : public FVKernel,
    */
   std::pair<SubdomainID, SubdomainID> faceArgSubdomains(const FaceInfo * face_info = nullptr) const;
 
+  /**
+   * Determine the single sided face argument when evaluating a functor on a face.
+   * This is used to perform evluations of material properties with the actual face values of
+   * their dependences, rather than interpolate the material property to the boundary.
+   * @param fi the FaceInfo for this face
+   * @param limiter_type the limiter type, to be specified if more than the default average
+   *        interpolation is required for the parameters of the functor
+   * @param correct_skewness whether to perform skew correction at the face
+   */
+  Moose::SingleSidedFaceArg singleSidedFaceArg(
+      const FaceInfo * fi = nullptr,
+      Moose::FV::LimiterType limiter_type = Moose::FV::LimiterType::CentralDifference,
+      bool correct_skewness = false) const;
+
   const bool _force_boundary_execution;
 
   std::unordered_set<BoundaryID> _boundaries_to_force;
diff --git a/framework/src/fvbcs/FVBoundaryCondition.C b/framework/src/fvbcs/FVBoundaryCondition.C
index df338bb378c4..7497eee2ffd4 100644
--- a/framework/src/fvbcs/FVBoundaryCondition.C
+++ b/framework/src/fvbcs/FVBoundaryCondition.C
@@ -92,14 +92,7 @@ FVBoundaryCondition::singleSidedFaceArg(const FaceInfo * fi,
   if (!fi)
     fi = _face_info;
   const bool use_elem = fi->faceType(_var.name()) == FaceInfo::VarFaceNeighbors::ELEM;
+  const auto sub_id = use_elem ? fi->elem().subdomain_id() : fi->neighborPtr()->subdomain_id();
 
-  if (use_elem)
-    return {fi, limiter_type, true, correct_skewness, correct_skewness, fi->elem().subdomain_id()};
-  else
-    return {fi,
-            limiter_type,
-            true,
-            correct_skewness,
-            correct_skewness,
-            fi->neighborPtr()->subdomain_id()};
+  return {fi, limiter_type, true, correct_skewness, correct_skewness, sub_id};
 }
diff --git a/framework/src/fvkernels/FVFluxKernel.C b/framework/src/fvkernels/FVFluxKernel.C
index f649842d8273..06bc2548c6d3 100644
--- a/framework/src/fvkernels/FVFluxKernel.C
+++ b/framework/src/fvkernels/FVFluxKernel.C
@@ -338,3 +338,16 @@ FVFluxKernel::faceArgSubdomains(const FaceInfo * face_info) const
 
   return Moose::FV::faceArgSubdomains(*this, *face_info);
 }
+
+Moose::SingleSidedFaceArg
+FVFluxKernel::singleSidedFaceArg(const FaceInfo * fi,
+                                 const Moose::FV::LimiterType limiter_type,
+                                 const bool correct_skewness) const
+{
+  if (!fi)
+    fi = _face_info;
+  const bool use_elem = fi->faceType(_var.name()) == FaceInfo::VarFaceNeighbors::ELEM;
+  const auto sub_id = use_elem ? fi->elem().subdomain_id() : fi->neighborPtr()->subdomain_id();
+
+  return {fi, limiter_type, true, correct_skewness, correct_skewness, sub_id};
+}
diff --git a/modules/navier_stokes/include/fvkernels/INSFVAdvectionKernel.h b/modules/navier_stokes/include/fvkernels/INSFVAdvectionKernel.h
index 87fb6bc1ea23..988904f13147 100644
--- a/modules/navier_stokes/include/fvkernels/INSFVAdvectionKernel.h
+++ b/modules/navier_stokes/include/fvkernels/INSFVAdvectionKernel.h
@@ -23,6 +23,18 @@ class INSFVAdvectionKernel : public FVFluxKernel, public INSFVBCInterface
   INSFVAdvectionKernel(const InputParameters & params);
   void initialSetup() override;
 
+  /**
+   * Sets the advection and velocity interpolation methods
+   * @param obj The \p MooseObject with input parameters to query
+   * @param advected_interp_method The advected interpolation method we will set
+   * @param velocity_interp_method The velocity interpolation method we will set
+   * @return Whether the interpolation methods have indicated that we will need more than our base
+   * level of ghosting
+   */
+  static bool setInterpolationMethods(const MooseObject & obj,
+                                      Moose::FV::InterpMethod & advected_interp_method,
+                                      Moose::FV::InterpMethod & velocity_interp_method);
+
 protected:
   bool skipForBoundary(const FaceInfo & fi) const override;
 
diff --git a/modules/navier_stokes/src/fvkernels/INSFVAdvectionKernel.C b/modules/navier_stokes/src/fvkernels/INSFVAdvectionKernel.C
index 7f49547664f4..aaa18b55c51b 100644
--- a/modules/navier_stokes/src/fvkernels/INSFVAdvectionKernel.C
+++ b/modules/navier_stokes/src/fvkernels/INSFVAdvectionKernel.C
@@ -36,6 +36,51 @@ INSFVAdvectionKernel::validParams()
   return params;
 }
 
+bool
+INSFVAdvectionKernel::setInterpolationMethods(const MooseObject & obj,
+                                              Moose::FV::InterpMethod & advected_interp_method,
+                                              Moose::FV::InterpMethod & velocity_interp_method)
+{
+  using namespace Moose::FV;
+
+  bool need_more_ghosting = false;
+
+  const auto & advected_interp_method_in = obj.getParam<MooseEnum>("advected_interp_method");
+  if (advected_interp_method_in == "average")
+    advected_interp_method = InterpMethod::Average;
+  else if (advected_interp_method_in == "skewness-corrected")
+    advected_interp_method = Moose::FV::InterpMethod::SkewCorrectedAverage;
+  else if (advected_interp_method_in == "upwind")
+    advected_interp_method = InterpMethod::Upwind;
+  else
+  {
+    if (advected_interp_method_in == "sou")
+      advected_interp_method = InterpMethod::SOU;
+    else if (advected_interp_method_in == "min_mod")
+      advected_interp_method = InterpMethod::MinMod;
+    else if (advected_interp_method_in == "vanLeer")
+      advected_interp_method = InterpMethod::VanLeer;
+    else if (advected_interp_method_in == "quick")
+      advected_interp_method = InterpMethod::QUICK;
+    else
+      obj.mooseError("Unrecognized interpolation type ",
+                     static_cast<std::string>(advected_interp_method_in));
+
+    need_more_ghosting = true;
+  }
+
+  const auto & velocity_interp_method_in = obj.getParam<MooseEnum>("velocity_interp_method");
+  if (velocity_interp_method_in == "average")
+    velocity_interp_method = InterpMethod::Average;
+  else if (velocity_interp_method_in == "rc")
+    velocity_interp_method = InterpMethod::RhieChow;
+  else
+    obj.mooseError("Unrecognized interpolation type ",
+                   static_cast<std::string>(velocity_interp_method_in));
+
+  return need_more_ghosting;
+}
+
 INSFVAdvectionKernel::INSFVAdvectionKernel(const InputParameters & params)
   : FVFluxKernel(params),
     _rc_vel_provider(getUserObject<INSFVRhieChowInterpolator>("rhie_chow_user_object"))
@@ -45,66 +90,33 @@ INSFVAdvectionKernel::INSFVAdvectionKernel(const InputParameters & params)
              "configure script in the root MOOSE directory with the configure option "
              "'--with-ad-indexing-type=global'");
 #endif
-  using namespace Moose::FV;
-
-  const auto & advected_interp_method = getParam<MooseEnum>("advected_interp_method");
-  if (advected_interp_method == "average")
-    _advected_interp_method = InterpMethod::Average;
-  else if (advected_interp_method == "skewness-corrected")
-    _advected_interp_method = Moose::FV::InterpMethod::SkewCorrectedAverage;
-  else if (advected_interp_method == "upwind")
-    _advected_interp_method = InterpMethod::Upwind;
-  else
+  const bool need_more_ghosting =
+      setInterpolationMethods(*this, _advected_interp_method, _velocity_interp_method);
+  if (need_more_ghosting && _tid == 0)
   {
-    if (advected_interp_method == "sou")
-      _advected_interp_method = InterpMethod::SOU;
-    else if (advected_interp_method == "min_mod")
-      _advected_interp_method = InterpMethod::MinMod;
-    else if (advected_interp_method == "vanLeer")
-      _advected_interp_method = InterpMethod::VanLeer;
-    else if (advected_interp_method == "quick")
-      _advected_interp_method = InterpMethod::QUICK;
-    else
-      mooseError("Unrecognized interpolation type ",
-                 static_cast<std::string>(advected_interp_method));
-
-    if (_tid == 0)
-    {
-      auto & factory = _app.getFactory();
-
-      auto rm_params = factory.getValidParams("ElementSideNeighborLayers");
-
-      rm_params.set<std::string>("for_whom") = name();
-      rm_params.set<MooseMesh *>("mesh") = &const_cast<MooseMesh &>(_mesh);
-      rm_params.set<Moose::RelationshipManagerType>("rm_type") =
-          Moose::RelationshipManagerType::GEOMETRIC | Moose::RelationshipManagerType::ALGEBRAIC |
-          Moose::RelationshipManagerType::COUPLING;
-      FVKernel::setRMParams(
-          _pars,
-          rm_params,
-          std::max((unsigned short)(3), _pars.get<unsigned short>("ghost_layers")));
-      mooseAssert(rm_params.areAllRequiredParamsValid(),
-                  "All relationship manager parameters should be valid.");
-
-      auto rm_obj = factory.create<RelationshipManager>(
-          "ElementSideNeighborLayers", name() + "_skew_correction", rm_params);
-
-      // Delete the resources created on behalf of the RM if it ends up not being added to the
-      // App.
-      if (!_app.addRelationshipManager(rm_obj))
-        factory.releaseSharedObjects(*rm_obj);
-    }
+    auto & factory = _app.getFactory();
+
+    auto rm_params = factory.getValidParams("ElementSideNeighborLayers");
+
+    rm_params.set<std::string>("for_whom") = name();
+    rm_params.set<MooseMesh *>("mesh") = &const_cast<MooseMesh &>(_mesh);
+    rm_params.set<Moose::RelationshipManagerType>("rm_type") =
+        Moose::RelationshipManagerType::GEOMETRIC | Moose::RelationshipManagerType::ALGEBRAIC |
+        Moose::RelationshipManagerType::COUPLING;
+    FVKernel::setRMParams(
+        _pars, rm_params, std::max((unsigned short)(3), _pars.get<unsigned short>("ghost_layers")));
+    mooseAssert(rm_params.areAllRequiredParamsValid(),
+                "All relationship manager parameters should be valid.");
+
+    auto rm_obj = factory.create<RelationshipManager>(
+        "ElementSideNeighborLayers", name() + "_skew_correction", rm_params);
+
+    // Delete the resources created on behalf of the RM if it ends up not being added to the
+    // App.
+    if (!_app.addRelationshipManager(rm_obj))
+      factory.releaseSharedObjects(*rm_obj);
   }
 
-  const auto & velocity_interp_method = getParam<MooseEnum>("velocity_interp_method");
-  if (velocity_interp_method == "average")
-    _velocity_interp_method = InterpMethod::Average;
-  else if (velocity_interp_method == "rc")
-    _velocity_interp_method = InterpMethod::RhieChow;
-  else
-    mooseError("Unrecognized interpolation type ",
-               static_cast<std::string>(velocity_interp_method));
-
   auto param_check = [&params, this](const auto & param_name)
   {
     if (params.isParamSetByUser(param_name))
diff --git a/modules/navier_stokes/src/fvkernels/INSFVMassAdvection.C b/modules/navier_stokes/src/fvkernels/INSFVMassAdvection.C
index 1af153d6a0c4..a1c2203571bf 100644
--- a/modules/navier_stokes/src/fvkernels/INSFVMassAdvection.C
+++ b/modules/navier_stokes/src/fvkernels/INSFVMassAdvection.C
@@ -35,11 +35,15 @@ ADReal
 INSFVMassAdvection::computeQpResidual()
 {
   const auto v = _rc_vel_provider.getVelocity(_velocity_interp_method, *_face_info, _tid);
-  const auto rho_interface =
-      Moose::FV::interpolate(_rho,
-                             Moose::FV::makeFace(*_face_info,
-                                                 limiterType(_advected_interp_method),
-                                                 MetaPhysicL::raw_value(v) * _normal > 0,
-                                                 faceArgSubdomains()));
-  return _normal * v * rho_interface;
+  ADReal rho_face{};
+
+  if (onBoundary(*_face_info))
+    rho_face = _rho(singleSidedFaceArg());
+  else
+    rho_face = Moose::FV::interpolate(_rho,
+                                      Moose::FV::makeFace(*_face_info,
+                                                          limiterType(_advected_interp_method),
+                                                          MetaPhysicL::raw_value(v) * _normal > 0,
+                                                          faceArgSubdomains()));
+  return _normal * v * rho_face;
 }
diff --git a/modules/navier_stokes/src/fvkernels/INSFVMomentumAdvection.C b/modules/navier_stokes/src/fvkernels/INSFVMomentumAdvection.C
index a296c87cd46c..e4939a6e1596 100644
--- a/modules/navier_stokes/src/fvkernels/INSFVMomentumAdvection.C
+++ b/modules/navier_stokes/src/fvkernels/INSFVMomentumAdvection.C
@@ -11,6 +11,7 @@
 #include "NS.h"
 #include "FVUtils.h"
 #include "INSFVRhieChowInterpolator.h"
+#include "SystemBase.h"
 
 registerMooseObject("NavierStokesApp", INSFVMomentumAdvection);
 
@@ -43,9 +44,8 @@ INSFVMomentumAdvection::initialSetup()
 
   _rho_u = std::make_unique<PiecewiseByBlockLambdaFunctor<ADReal>>(
       "rho_u",
-      [this](const auto & r, const auto & t) -> ADReal {
-        return _rho(r, t) * _var(r, t) / epsilon()(r, t);
-      },
+      [this](const auto & r, const auto & t) -> ADReal
+      { return _rho(r, t) * _var(r, t) / epsilon()(r, t); },
       std::set<ExecFlagType>({EXEC_ALWAYS}),
       _mesh,
       this->blockIDs());
@@ -59,28 +59,48 @@ INSFVMomentumAdvection::computeQpResidual()
   const bool correct_skewness = _advected_interp_method == InterpMethod::SkewCorrectedAverage;
   const auto v = _rc_vel_provider.getVelocity(_velocity_interp_method, *_face_info, _tid);
 
-  const auto [interp_coeffs, advected] =
-      interpCoeffsAndAdvected(*_rho_u,
-                              makeFace(*_face_info,
-                                       limiterType(_advected_interp_method),
-                                       MetaPhysicL::raw_value(v) * _normal > 0,
-                                       faceArgSubdomains(),
-                                       correct_skewness,
-                                       correct_skewness));
-
-  const auto elem_face = elemFromFace();
-  const auto neighbor_face = neighborFromFace();
-
-  const auto rho_elem = _rho(elem_face), rho_neighbor = _rho(neighbor_face);
-  const auto eps_elem = epsilon()(elem_face), eps_neighbor = epsilon()(neighbor_face);
-  const auto var_elem = advected.first / rho_elem * eps_elem,
-             var_neighbor = advected.second / rho_neighbor * eps_neighbor;
-
-  _ae = _normal * v * rho_elem / eps_elem * interp_coeffs.first;
-  // Minus sign because we apply a minus sign to the residual in computeResidual
-  _an = -_normal * v * rho_neighbor / eps_neighbor * interp_coeffs.second;
-
-  return _ae * var_elem - _an * var_neighbor;
+  if (onBoundary(*_face_info))
+  {
+    const auto ssf = singleSidedFaceArg();
+    const Elem * const sided_elem = ssf.makeSidedElem().elem;
+    const auto dof_number = sided_elem->dof_number(_sys.number(), _var.number(), 0);
+    const auto rhof = _rho(ssf);
+    const auto epsf = epsilon()(ssf);
+    const auto uf = _var(ssf);
+    const Real duf_du = uf.derivatives()[dof_number];
+    const auto coeff = _normal * v * rhof / epsf;
+    if (sided_elem == &_face_info->elem())
+      _ae = coeff * duf_du;
+    else
+      _an = -coeff * duf_du;
+
+    return coeff * uf;
+  }
+  else
+  {
+    const auto [interp_coeffs, advected] =
+        interpCoeffsAndAdvected(*_rho_u,
+                                makeFace(*_face_info,
+                                         limiterType(_advected_interp_method),
+                                         MetaPhysicL::raw_value(v) * _normal > 0,
+                                         faceArgSubdomains(),
+                                         correct_skewness,
+                                         correct_skewness));
+
+    const auto elem_face = elemFromFace();
+    const auto neighbor_face = neighborFromFace();
+
+    const auto rho_elem = _rho(elem_face), rho_neighbor = _rho(neighbor_face);
+    const auto eps_elem = epsilon()(elem_face), eps_neighbor = epsilon()(neighbor_face);
+    const auto var_elem = advected.first / rho_elem * eps_elem,
+               var_neighbor = advected.second / rho_neighbor * eps_neighbor;
+
+    _ae = _normal * v * rho_elem / eps_elem * interp_coeffs.first;
+    // Minus sign because we apply a minus sign to the residual in computeResidual
+    _an = -_normal * v * rho_neighbor / eps_neighbor * interp_coeffs.second;
+
+    return _ae * var_elem - _an * var_neighbor;
+  }
 }
 
 void
diff --git a/modules/navier_stokes/src/postprocessors/VolumetricFlowRate.C b/modules/navier_stokes/src/postprocessors/VolumetricFlowRate.C
index eb5616423a66..45d173968d6d 100644
--- a/modules/navier_stokes/src/postprocessors/VolumetricFlowRate.C
+++ b/modules/navier_stokes/src/postprocessors/VolumetricFlowRate.C
@@ -10,6 +10,7 @@
 #include "VolumetricFlowRate.h"
 #include "MathFVUtils.h"
 #include "INSFVRhieChowInterpolator.h"
+#include "INSFVAdvectionKernel.h"
 #include <math.h>
 
 registerMooseObject("NavierStokesApp", VolumetricFlowRate);
@@ -86,25 +87,8 @@ VolumetricFlowRate::VolumetricFlowRate(const InputParameters & parameters)
                  "supplied.");
   }
 
-  using namespace Moose::FV;
-
-  const auto & advected_interp_method = getParam<MooseEnum>("advected_interp_method");
-  if (advected_interp_method == "average")
-    _advected_interp_method = InterpMethod::Average;
-  else if (advected_interp_method == "upwind")
-    _advected_interp_method = InterpMethod::Upwind;
-  else
-    mooseError("Unrecognized advected quantity interpolation type ",
-               static_cast<std::string>(advected_interp_method));
-
-  const auto & velocity_interp_method = getParam<MooseEnum>("velocity_interp_method");
-  if (velocity_interp_method == "average")
-    _velocity_interp_method = InterpMethod::Average;
-  else if (velocity_interp_method == "rc")
-    _velocity_interp_method = InterpMethod::RhieChow;
-  else
-    mooseError("Unrecognized interpolation type ",
-               static_cast<std::string>(velocity_interp_method));
+  INSFVAdvectionKernel::setInterpolationMethods(
+      *this, _advected_interp_method, _velocity_interp_method);
 }
 
 Real
@@ -120,19 +104,17 @@ VolumetricFlowRate::computeQpIntegral()
     // Get face value for velocity
     const auto vel =
         MetaPhysicL::raw_value(_rc_uo->getVelocity(_velocity_interp_method, *fi, _tid));
+    const bool correct_skewness =
+        _advected_interp_method == Moose::FV::InterpMethod::SkewCorrectedAverage;
 
-    const auto elem_face = Moose::FV::elemFromFace(*_rc_uo, *fi);
-    const auto neighbor_face = Moose::FV::neighborFromFace(*_rc_uo, *fi);
+    const auto ssf = Moose::SingleSidedFaceArg({fi,
+                                                Moose::FV::LimiterType::CentralDifference,
+                                                true,
+                                                correct_skewness,
+                                                correct_skewness,
+                                                _current_elem->subdomain_id()});
 
-    Real adv_quant_interface;
-    Moose::FV::interpolate(_advected_interp_method,
-                           adv_quant_interface,
-                           MetaPhysicL::raw_value((*_adv_quant)(elem_face)),
-                           MetaPhysicL::raw_value((*_adv_quant)(neighbor_face)),
-                           vel,
-                           *fi,
-                           true);
-    return fi->normal() * vel * adv_quant_interface;
+    return fi->normal() * MetaPhysicL::raw_value((*_adv_quant)(ssf)) * vel;
   }
   else
 #endif
diff --git a/modules/navier_stokes/test/tests/finite_volume/fvbcs/wall_function/gold/Re_t395_out.e b/modules/navier_stokes/test/tests/finite_volume/fvbcs/wall_function/gold/Re_t395_out.e
index 316fc9826a4b..0a02620870b0 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/fvbcs/wall_function/gold/Re_t395_out.e and b/modules/navier_stokes/test/tests/finite_volume/fvbcs/wall_function/gold/Re_t395_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/fviks/convection/gold/convection_channel_out.e b/modules/navier_stokes/test/tests/finite_volume/fviks/convection/gold/convection_channel_out.e
index f9108121e800..ac58e0754f49 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/fviks/convection/gold/convection_channel_out.e and b/modules/navier_stokes/test/tests/finite_volume/fviks/convection/gold/convection_channel_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/fvkernels/flow_diode/gold/friction_out.e b/modules/navier_stokes/test/tests/finite_volume/fvkernels/flow_diode/gold/friction_out.e
index 5948270525e6..0b288d59c1fc 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/fvkernels/flow_diode/gold/friction_out.e and b/modules/navier_stokes/test/tests/finite_volume/fvkernels/flow_diode/gold/friction_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/block_restriction/gold/2d-rc_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/block_restriction/gold/2d-rc_out.e
index fd4d7ecdf75b..edf68674a5e0 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/block_restriction/gold/2d-rc_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/block_restriction/gold/2d-rc_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/2d-rc-no-slip-xyz-mass-conservation.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/2d-rc-no-slip-xyz-mass-conservation.e
index b192e9984ab6..140089f3bbfb 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/2d-rc-no-slip-xyz-mass-conservation.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/2d-rc-no-slip-xyz-mass-conservation.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/diverging-free-slip.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/diverging-free-slip.e
index dedb9316d30c..c22df7abd0cc 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/diverging-free-slip.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/diverging-free-slip.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/no-slip-tris_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/no-slip-tris_out.e
index a53bccae01c8..f9b9371c0c83 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/no-slip-tris_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/gold/no-slip-tris_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-mixing-length_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-mixing-length_out.e
index f02482e5a13f..b203dedef7cb 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-mixing-length_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-mixing-length_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-friction-quad_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-friction-quad_out.e
index f684d02b07bd..9c12c55a0c4a 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-friction-quad_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-friction-quad_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-friction_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-friction_out.e
index c01430d1af41..78fc45746532 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-friction_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-friction_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip-extrapolated-outlet-pressure_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip-extrapolated-outlet-pressure_out.e
index 017c79253cac..c58c1eaabfae 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip-extrapolated-outlet-pressure_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip-extrapolated-outlet-pressure_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip-tris.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip-tris.e
index d2e92285fbd0..802a03df3362 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip-tris.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip-tris.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip_out.e
index 46245cc46ccb..a8e0c9420094 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-rc-no-slip_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-scalar-transport_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-scalar-transport_out.e
index ce49753eb93d..0185fd80f53d 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-scalar-transport_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/2d-scalar-transport_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/3d-rc-no-slip_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/3d-rc-no-slip_out.e
index 28e9976715c7..9dabd42d2b0b 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/3d-rc-no-slip_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/3d-rc-no-slip_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/discontinuous-body-forces_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/discontinuous-body-forces_out.e
index 565282b35861..f778c5b918d1 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/discontinuous-body-forces_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/gold/discontinuous-body-forces_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/gold/mixing_length_eddy_viscosity_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/gold/mixing_length_eddy_viscosity_out.e
index 8390fcedd793..e616eca16d5b 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/gold/mixing_length_eddy_viscosity_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/gold/mixing_length_eddy_viscosity_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/mixing_length_eddy_viscosity.i b/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/mixing_length_eddy_viscosity.i
index be3cf08cbd12..499fb32a3609 100644
--- a/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/mixing_length_eddy_viscosity.i
+++ b/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/mixing_length_eddy_viscosity.i
@@ -3,7 +3,7 @@ von_karman_const = 0.41
 H = 1 #halfwidth of the channel
 L = 150
 
-Re = 13700
+Re = 100
 
 rho = 1
 bulk_u = 1
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/tests b/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/tests
index 67169ccf7ee9..115a42375c8d 100644
--- a/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/tests
+++ b/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_eddy_viscosity_aux/tests
@@ -13,5 +13,6 @@
     mesh_mode = REPLICATED
     valgrind = 'none'
     recover = false
+    abs_zero = 3e-10
   []
 []
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_total_viscosity_material/gold/mixing_length_total_viscosity_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_total_viscosity_material/gold/mixing_length_total_viscosity_out.e
index 132c6a5b3e4d..7b2efabf7bd5 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_total_viscosity_material/gold/mixing_length_total_viscosity_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/mixing_length_total_viscosity_material/gold/mixing_length_total_viscosity_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/multiapp-scalar-transport/gold/fluid-flow_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/multiapp-scalar-transport/gold/fluid-flow_out.e
index 07247fe1b0af..646c663583e5 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/multiapp-scalar-transport/gold/fluid-flow_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/multiapp-scalar-transport/gold/fluid-flow_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/multiapp-scalar-transport/gold/fluid-flow_out_scalar0.e b/modules/navier_stokes/test/tests/finite_volume/ins/multiapp-scalar-transport/gold/fluid-flow_out_scalar0.e
index 336bf5396111..3fa73a9eb07a 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/multiapp-scalar-transport/gold/fluid-flow_out_scalar0.e and b/modules/navier_stokes/test/tests/finite_volume/ins/multiapp-scalar-transport/gold/fluid-flow_out_scalar0.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/capped_mixing_length.i b/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/capped_mixing_length.i
index f0f4fc729793..d7c099762b4a 100644
--- a/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/capped_mixing_length.i
+++ b/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/capped_mixing_length.i
@@ -3,7 +3,7 @@ von_karman_const = 0.41
 H = 1 #halfwidth of the channel
 L = 150
 
-Re = 13700
+Re = 100
 
 rho = 1
 bulk_u = 1
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/gold/capped_mixing_length_out.e b/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/gold/capped_mixing_length_out.e
index f39f95914596..124b39fbdb1a 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/gold/capped_mixing_length_out.e and b/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/gold/capped_mixing_length_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/tests b/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/tests
index 387e1e82e45f..97af23f33ab3 100644
--- a/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/tests
+++ b/modules/navier_stokes/test/tests/finite_volume/ins/wall_distance_capped_mixing_length_aux/tests
@@ -13,5 +13,6 @@
     mesh_mode = REPLICATED
     valgrind = 'none'
     recover = false
+    abs_zero = 3e-10
   []
 []
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/2d-rc-rz-by-parts_out.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/2d-rc-rz-by-parts_out.e
index 3f234bcef49d..894553856a17 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/2d-rc-rz-by-parts_out.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/2d-rc-rz-by-parts_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-friction.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-friction.e
index fe0a972c708a..76f3bf8e1575 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-friction.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-friction.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-match-insfv.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-match-insfv.e
index ca8553bbffa0..5851498b1b0f 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-match-insfv.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-match-insfv.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-mean-pressure.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-mean-pressure.e
index e439c2d74fad..a3cf18004681 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-mean-pressure.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-mean-pressure.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-rz-symmetry.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-rz-symmetry.e
index dee7079e65f2..0c864b147358 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-rz-symmetry.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-rz-symmetry.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-symmetry.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-symmetry.e
index d98c0c18c95f..b23bc4bd1e7f 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-symmetry.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip-symmetry.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip.e
index 2f29e3e753d8..32c3d0296526 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/gold/rc-no-slip.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/2d-transient_out.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/2d-transient_out.e
index 595ccce97688..85bee813c6a8 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/2d-transient_out.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/2d-transient_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-boussinesq.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-boussinesq.e
index f8f7870ed28d..484812902bdf 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-boussinesq.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-boussinesq.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-effective.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-effective.e
index 00ddd906c972..fd86287be553 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-effective.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-effective.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-rz.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-rz.e
index 6a2b10bb8a4a..80b366885198 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-rz.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-rz.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-solid-fluid.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-solid-fluid.e
index b08eb724acdd..c128f3cc5354 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-solid-fluid.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated-solid-fluid.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated.e
index f2b72fe614a7..51e9a0689086 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/gold/rc-heated.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/tests b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/tests
index 525204d89ad6..b93dc2148ec6 100644
--- a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/tests
+++ b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/heated/tests
@@ -19,7 +19,7 @@
       input = 2d-rc-heated.i
       exodiff = rc-heated-rz.e
       method = "!dbg"
-      abs_zero = 1e-9
+      abs_zero = 2e-9
       detail = 'in rz geometry,'
       ad_indexing_type = 'global'
       valgrind = HEAVY
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/rc-discontinuous-epsjump.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/rc-discontinuous-epsjump.e
index 64b1cef4ed9c..e08ffe306eac 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/rc-discontinuous-epsjump.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/rc-discontinuous-epsjump.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/rc-smooth-epsjump.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/rc-smooth-epsjump.e
index f0b955ee95ec..4f3398c0005a 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/rc-smooth-epsjump.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/rc-smooth-epsjump.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/reconstructed.e b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/reconstructed.e
index 239190f8cc8d..7d6f44508059 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/reconstructed.e and b/modules/navier_stokes/test/tests/finite_volume/pins/channel-flow/porosity_jump/gold/reconstructed.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/pwcns/channel-flow/gold/2d-transient_out.e b/modules/navier_stokes/test/tests/finite_volume/pwcns/channel-flow/gold/2d-transient_out.e
index 19569e8bb714..7cc678c3c9b6 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/pwcns/channel-flow/gold/2d-transient_out.e and b/modules/navier_stokes/test/tests/finite_volume/pwcns/channel-flow/gold/2d-transient_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/dirichlet_bcs_mdot_out.e b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/dirichlet_bcs_mdot_out.e
index ac5d5b04fb06..be3f1aa6f7ed 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/dirichlet_bcs_mdot_out.e and b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/dirichlet_bcs_mdot_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/dirichlet_bcs_velocity_out.e b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/dirichlet_bcs_velocity_out.e
index cc770cc717b5..3dcf0e796b10 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/dirichlet_bcs_velocity_out.e and b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/dirichlet_bcs_velocity_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_direct_out.e b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_direct_out.e
index f8e826c20bcd..596fe6404c9c 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_direct_out.e and b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_direct_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_mdot_out.e b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_mdot_out.e
index 34a84a0870e3..0cccfd6596b2 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_mdot_out.e and b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_mdot_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_velocity_out.e b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_velocity_out.e
index d5e0b5c8b688..376a83a02a03 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_velocity_out.e and b/modules/navier_stokes/test/tests/finite_volume/wcns/boundary_conditions/gold/flux_bcs_velocity_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/channel-flow/gold/2d-transient_out.e b/modules/navier_stokes/test/tests/finite_volume/wcns/channel-flow/gold/2d-transient_out.e
index 2d1d95065684..e60556d65398 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/wcns/channel-flow/gold/2d-transient_out.e and b/modules/navier_stokes/test/tests/finite_volume/wcns/channel-flow/gold/2d-transient_out.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/channel-flow/gold/2d-turbulent.e b/modules/navier_stokes/test/tests/finite_volume/wcns/channel-flow/gold/2d-turbulent.e
index 0cb41271b0cc..aee52c3469ad 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/wcns/channel-flow/gold/2d-turbulent.e and b/modules/navier_stokes/test/tests/finite_volume/wcns/channel-flow/gold/2d-turbulent.e differ
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/materials/gold/2d-transient_out.e b/modules/navier_stokes/test/tests/finite_volume/wcns/materials/gold/2d-transient_out.e
index e2fb9f1bc808..644df577058d 100644
Binary files a/modules/navier_stokes/test/tests/finite_volume/wcns/materials/gold/2d-transient_out.e and b/modules/navier_stokes/test/tests/finite_volume/wcns/materials/gold/2d-transient_out.e differ
diff --git a/modules/navier_stokes/test/tests/postprocessors/flow_rates/gold/diverging_INSFV_quad_xy_noslip.csv b/modules/navier_stokes/test/tests/postprocessors/flow_rates/gold/diverging_INSFV_quad_xy_noslip.csv
index afb34955369d..27620c6d58ea 100644
--- a/modules/navier_stokes/test/tests/postprocessors/flow_rates/gold/diverging_INSFV_quad_xy_noslip.csv
+++ b/modules/navier_stokes/test/tests/postprocessors/flow_rates/gold/diverging_INSFV_quad_xy_noslip.csv
@@ -1,3 +1,3 @@
 time,inlet_mass_constant,inlet_mass_matprop,inlet_mass_variable,inlet_momentum_x,inlet_momentum_y,mid1_advected_energy,mid1_mass,mid2_advected_energy,mid2_mass,outlet_advected_energy,outlet_mass
 0,0,0,0,0,0,0,0,0,0,0,0
-1,-1,-1,-1,0,-1,310,1,310,1,310,1
+1,-1,-1,-1,0,-1,310,1,310,1,309.98396702432,1
diff --git a/modules/navier_stokes/test/tests/postprocessors/flow_rates/gold/diverging_PINSFV_quad_xy_noslip.csv b/modules/navier_stokes/test/tests/postprocessors/flow_rates/gold/diverging_PINSFV_quad_xy_noslip.csv
index afb34955369d..afa6b24b88a3 100644
--- a/modules/navier_stokes/test/tests/postprocessors/flow_rates/gold/diverging_PINSFV_quad_xy_noslip.csv
+++ b/modules/navier_stokes/test/tests/postprocessors/flow_rates/gold/diverging_PINSFV_quad_xy_noslip.csv
@@ -1,3 +1,3 @@
 time,inlet_mass_constant,inlet_mass_matprop,inlet_mass_variable,inlet_momentum_x,inlet_momentum_y,mid1_advected_energy,mid1_mass,mid2_advected_energy,mid2_mass,outlet_advected_energy,outlet_mass
 0,0,0,0,0,0,0,0,0,0,0,0
-1,-1,-1,-1,0,-1,310,1,310,1,310,1
+1,-1,-1,-1,0,-1,310,1,310,1,309.98506259795,1