Use one-term boundary expansion by default for INSFV

This fixes idaholab#16822 although perhaps not 100% satisfactorily. If a user
knows they only have cells with a maximum of one extrapolated boundary
face or if they know that if the cell has multiple extrapolated boundary
faces that it is orthogonal (e.g. dCf is parallel to Sf), then they can
safely set `two_term_boundary_expansion = true` in their input. But in
the interest of allowing more things to just work out of the box (like
any old tri or tet mesh), I think it makes sense to make the one term
expansion the default. I prefer doing this over something like
projecting dCf onto the surface normal because this would cost a user
accuracy in the case where they have intelligently designed their mesh
such that there is only a maximum of one extrapolated boundary face per
cell (gmsh meshes for example avoid this).

modules/navier_stokes/src/variables/INSFVVariable.C

@@ -20,11 +20,7 @@
 InputParameters
 INSFVVariable::validParams()
 {
-  auto params = MooseVariableFVReal::validParams();
-#ifdef MOOSE_GLOBAL_AD_INDEXING
-  params.set<bool>("two_term_boundary_expansion") = true;
-#endif
-  return params;
+  return MooseVariableFVReal::validParams();
 }
 
 INSFVVariable::INSFVVariable(const InputParameters & params) : MooseVariableFVReal(params)

modules/navier_stokes/test/tests/finite_volume/ins/block_restriction/2d-rc.i

@@ -5,6 +5,10 @@ velocity_interp_method='rc'
 
 restricted_blocks = '1'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   parallel_type = 'replicated'
   [mesh]

modules/navier_stokes/test/tests/finite_volume/ins/boussinesq/boussinesq.i

@@ -10,6 +10,10 @@ rayleigh=1e3
 hot_temp=${rayleigh}
 temp_ref=${fparse hot_temp / 2.}
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/2d-mixing-length.i

@@ -9,6 +9,10 @@ mu = ${fparse rho * bulk_u * D / Re}
 advected_interp_method='upwind'
 velocity_interp_method='rc'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/2d-rc-ambient-convection.i

@@ -5,6 +5,10 @@ cp=1
 advected_interp_method='average'
 velocity_interp_method='rc'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/2d-rc-friction.i

@@ -3,6 +3,10 @@ rho=1.1
 advected_interp_method='average'
 velocity_interp_method='rc'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/2d-rc-no-slip-extrapolated-outlet-pressure.i

@@ -3,6 +3,10 @@ rho=1.1
 advected_interp_method='average'
 velocity_interp_method='rc'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/2d-rc-no-slip-outflow-bcs.i

@@ -3,6 +3,10 @@ rho=1.1
 advected_interp_method='average'
 velocity_interp_method='rc'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/2d-rc-no-slip.i

@@ -3,6 +3,10 @@ rho=1.1
 advected_interp_method='average'
 velocity_interp_method='rc'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/2d-scalar-transport.i

@@ -6,6 +6,10 @@ cp=1
 advected_interp_method='average'
 velocity_interp_method='rc'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/2d-average-no-slip.i

@@ -3,6 +3,10 @@ rho=1
 advected_interp_method='average'
 velocity_interp_method='average'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/channel-flow/cylindrical/diverging.i

@@ -3,6 +3,10 @@ rho=1
 advected_interp_method='average'
 velocity_interp_method='rc'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   file = diverging.msh
   uniform_refine = 2

modules/navier_stokes/test/tests/finite_volume/ins/lid-driven/lid-driven-with-energy.i

@@ -6,6 +6,10 @@ vel = 'velocity'
 velocity_interp_method = 'rc'
 advected_interp_method = 'average'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator

modules/navier_stokes/test/tests/finite_volume/ins/lid-driven/lid-driven.i

@@ -5,6 +5,7 @@ rho=1
   vel = 'velocity'
   velocity_interp_method = 'rc'
   advected_interp_method = 'average'
+  two_term_boundary_expansion = true
 []
 
 [Mesh]

modules/navier_stokes/test/tests/finite_volume/ins/lid-driven/transient-lid-driven-with-energy.i

@@ -6,6 +6,10 @@ vel = 'velocity'
 velocity_interp_method = 'rc'
 advected_interp_method = 'average'
 
+[GlobalParams]
+  two_term_boundary_expansion = true
+[]
+
 [Mesh]
   [gen]
     type = GeneratedMeshGenerator