diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/materials/2d-steady-wall-balance.i b/modules/navier_stokes/test/tests/finite_volume/wcns/materials/2d-steady-wall-balance.i
new file mode 100644
index 000000000000..d0d5499fc130
--- /dev/null
+++ b/modules/navier_stokes/test/tests/finite_volume/wcns/materials/2d-steady-wall-balance.i
@@ -0,0 +1,261 @@
+L = 30
+bulk_u = 0.01
+p_ref = 101325.0
+T_in = 860
+q_source = 50000
+q2_wall = 10000
+
+A_cp = 976.78
+B_cp = 1.0634
+rho = 2000
+
+advected_interp_method = 'upwind'
+
+[Mesh]
+  [gmg]
+    type = GeneratedMeshGenerator
+    dim = 2
+    xmin = 0
+    xmax = ${L}
+    ymin = 1
+    ymax = 2.5
+    nx = 10
+  []
+[]
+
+[GlobalParams]
+  rhie_chow_user_object = 'rc'
+  advected_interp_method = ${advected_interp_method}
+  velocity_interp_method = 'rc'
+[]
+
+[UserObjects]
+  [rc]
+    type = INSFVRhieChowInterpolator
+    u = vel_x
+    v = vel_y
+    pressure = pressure
+  []
+[]
+
+[Variables]
+  [vel_x]
+    type = INSFVVelocityVariable
+    initial_condition = ${bulk_u}
+    two_term_boundary_expansion = false
+  []
+  [vel_y]
+    type = INSFVVelocityVariable
+    initial_condition = 0
+    two_term_boundary_expansion = false
+  []
+  [pressure]
+    type = INSFVPressureVariable
+    initial_condition = ${p_ref}
+    two_term_boundary_expansion = false
+  []
+  [T]
+    type = INSFVEnergyVariable
+    two_term_boundary_expansion = false
+    initial_condition = ${T_in}
+  []
+[]
+
+[FVKernels]
+  [mass]
+    type = WCNSFVMassAdvection
+    variable = pressure
+    rho = 'rho'
+  []
+
+  [u_advection]
+    type = INSFVMomentumAdvection
+    variable = vel_x
+    rho = 'rho'
+    momentum_component = 'x'
+  []
+  [u_viscosity]
+    type = INSFVMomentumDiffusion
+    variable = vel_x
+    mu = 'mu'
+    momentum_component = 'x'
+  []
+  [u_pressure]
+    type = INSFVMomentumPressure
+    variable = vel_x
+    momentum_component = 'x'
+    pressure = pressure
+  []
+  [v_advection]
+    type = INSFVMomentumAdvection
+    variable = vel_y
+    rho = 'rho'
+    momentum_component = 'y'
+  []
+  [v_viscosity]
+    type = INSFVMomentumDiffusion
+    variable = vel_y
+    mu = 'mu'
+    momentum_component = 'y'
+  []
+  [v_pressure]
+    type = INSFVMomentumPressure
+    variable = vel_y
+    momentum_component = 'y'
+    pressure = pressure
+  []
+
+  [temp_conduction]
+    type = FVDiffusion
+    coeff = 'k'
+    variable = T
+  []
+  [temp_advection]
+    type = INSFVEnergyAdvection
+    variable = T
+  []
+  [source]
+    type = FVBodyForce
+    variable = T
+    function = source_func
+  []
+[]
+
+[FVBCs]
+  [inlet-u]
+    type = INSFVInletVelocityBC
+    boundary = 'left'
+    variable = vel_x
+    functor = ${bulk_u}
+  []
+  [inlet-v]
+    type = INSFVInletVelocityBC
+    boundary = 'left'
+    variable = vel_y
+    functor = 0
+  []
+  [inlet_T]
+    type = FVDirichletBC
+    variable = T
+    boundary = 'left'
+    value = ${T_in}
+  []
+
+  [incoming_heat]
+    type = FVNeumannBC
+    variable = T
+    value = ${q2_wall}
+    boundary = 'top'
+  []
+
+  [outlet_p]
+    type = INSFVOutletPressureBC
+    boundary = 'right'
+    variable = pressure
+    function = ${p_ref}
+  []
+[]
+
+[Functions]
+  [source_func]
+    type = ParsedFunction
+    expression = '${q_source}'
+  []
+[]
+
+[FunctorMaterials]
+  [converter_to_regular_T]
+    type = FunctorADConverter
+    ad_props_in = 'T'
+    reg_props_out = 'T_nAD'
+  []
+  [ins_fv]
+    type = INSFVEnthalpyFunctorMaterial
+    temperature = 'T'
+    rho = 'rho'
+    cp = 'cp'
+
+    assume_constant_cp = false
+    h_in = 'h'
+    # fp = 'fp'
+    # pressure = 'pressure'
+  []
+
+  [rho]
+    type = ADParsedFunctorMaterial
+    property_name = 'rho'
+    expression = '${rho}'
+  []
+  [mu]
+    type = ADParsedFunctorMaterial
+    property_name = 'mu'
+    expression = '4.5e-3'
+  []
+  [k]
+    type = ADParsedFunctorMaterial
+    property_name = 'k'
+    expression = '0.7'
+  []
+  [h]
+    type = ADParsedFunctorMaterial
+    property_name = 'h'
+    functor_names = 'T ${A_cp} ${B_cp}'
+    functor_symbols = 'T A_cp B_cp'
+    expression = 'A_cp * T + B_cp * T * T / 2'
+  []
+  [cp]
+    type = ADParsedFunctorMaterial
+    property_name = 'cp'
+    functor_names = 'T ${A_cp} ${B_cp}'
+    functor_symbols = 'T A_cp B_cp'
+    expression = 'A_cp+B_cp*T'
+  []
+[]
+
+[Executioner]
+  type = Steady
+  solve_type = 'NEWTON'
+  petsc_options_iname = '-pc_type -pc_factor_shift_type'
+  petsc_options_value = 'lu       NONZERO'
+
+  nl_abs_tol = 1e-9
+  nl_max_its = 50
+  line_search = 'none'
+
+  automatic_scaling = true
+  off_diagonals_in_auto_scaling = true
+[]
+
+[Postprocessors]
+  [H_in]
+    type = VolumetricFlowRate
+    vel_x = 'vel_x'
+    advected_quantity = 'rho_h'
+    boundary = 'left'
+  []
+  [H_out]
+    type = VolumetricFlowRate
+    vel_x = 'vel_x'
+    advected_quantity = 'rho_h'
+    boundary = 'right'
+  []
+  [Q]
+    type = FunctionElementIntegral
+    function = 'source_func'
+    execute_on = 'initial'
+  []
+  [Q_wall]
+    type = FunctionSideIntegral
+    function = ${q2_wall}
+    boundary = 'top'
+  []
+  [balance_in_percent]
+    type = ParsedPostprocessor
+    expression = '(H_out + H_in - Q - Q_wall) / H_in * 100'
+    pp_names = 'H_in H_out Q Q_wall'
+  []
+[]
+
+[Outputs]
+  csv = true
+[]
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/materials/gold/2d-steady-wall-balance_out.csv b/modules/navier_stokes/test/tests/finite_volume/wcns/materials/gold/2d-steady-wall-balance_out.csv
new file mode 100644
index 000000000000..417c295edd4a
--- /dev/null
+++ b/modules/navier_stokes/test/tests/finite_volume/wcns/materials/gold/2d-steady-wall-balance_out.csv
@@ -0,0 +1,3 @@
+time,H_in,H_out,Q,Q_wall,balance_in_percent
+0,0,0,2250000,0,0
+1,-36998283.6,39548285.350731,2250000,300000,-4.7319255567647e-06
diff --git a/modules/navier_stokes/test/tests/finite_volume/wcns/materials/tests b/modules/navier_stokes/test/tests/finite_volume/wcns/materials/tests
index 66c7292a6f5c..526f8906bbcb 100644
--- a/modules/navier_stokes/test/tests/finite_volume/wcns/materials/tests
+++ b/modules/navier_stokes/test/tests/finite_volume/wcns/materials/tests
@@ -57,11 +57,18 @@
     requirement = 'The system shall be able to compute the specific enthalpy and the enthalpy with a non constant specific heat.'
     recover = false
   []
-  [validation_nonconstant_cp]
+  [verification_nonconstant_cp]
     type = 'CSVDiff'
     input = 1d_test_cpT.i
     csvdiff = 1d_test_cpT_out.csv
-    requirement = 'The system shall be able to compute the specific enthalpy and the enthalpy with a non constant specific heat, and show near perfect energy conservation.'
+    requirement = 'The system shall be able to compute the specific enthalpy and the enthalpy with a non constant specific heat, and show near perfect energy conservation in a 1D case with a volumetric source.'
+    recover = false
+  []
+  [verification_nonconstant_cp]
+    type = 'CSVDiff'
+    input = 2d-steady-wall-balance.i
+    csvdiff = 2d-steady-wall-balance_out.csv
+    requirement = 'The system shall be able to compute the specific enthalpy and the enthalpy with a non constant specific heat, and show near perfect energy conservation in a 2D case with wall heating.'
     recover = false
   []