fix: use full reconstruction with int_comp to enable THINC with viscous flows

Author: sbryngelson

examples/3D_rayleigh_taylor_muscl/case.py

@@ -57,6 +57,7 @@
     "recon_type": 2,
     "muscl_order": 2,
     "muscl_lim": 4,
+    "int_comp": 1,
     "avg_state": 2,
     "riemann_solver": 2,
     "wave_speeds": 1,

src/simulation/m_rhs.fpp

@@ -655,7 +655,7 @@ contains
                 call nvtxStartRange("RHS-RECONSTRUCTION")
 
                 if (.not. surface_tension) then
-                    if (all(Re_size == 0)) then
+                    if (all(Re_size == 0) .or. int_comp > 0) then
                         ! Reconstruct densitiess
                         iv%beg = 1; iv%end = sys_size
                         call s_reconstruct_cell_boundary_values(q_prim_qp%vf(1:sys_size), qL_rsx_vf, qL_rsy_vf, qL_rsz_vf, &

toolchain/mfc/case_validator.py

@@ -391,10 +391,6 @@ def check_interface_compression(self):
         self.prohibit(int_comp not in [0, 1, 2], "int_comp must be 0 (off), 1 (THINC), or 2 (MTHINC)")
         self.prohibit(int_comp == 2 and n == 0, "int_comp = 2 (MTHINC) requires at least 2D (n > 0)")
         self.prohibit(int_comp != 0 and num_fluids != 2, "int_comp > 0 requires num_fluids = 2")
-        viscous = self.get("viscous", "F") == "T"
-        surface_tension = self.get("surface_tension", "F") == "T"
-        self.prohibit(int_comp != 0 and viscous, "int_comp > 0 is not supported with viscosity (reconstruction is split; compression would be silently skipped)")
-        self.prohibit(int_comp != 0 and surface_tension, "int_comp > 0 is not supported with surface tension (reconstruction is split; compression would be silently skipped)")
 
         recon_type = self.get("recon_type", 1)
         if recon_type == 1:  # WENO