Issue with Stabilization Type in MultiphasePoromechanics Solver

I am currently working with the MultiphasePoromechanics Solver in a thermal-hydro-mechanical coupling scenario. I have coupled the CompositionalMultiphaseFVM solver and the SolidMechanics_LagrangianFEM solver. However, I have encountered an issue when selecting the stabilization type.

• When using the Global stabilization type, the simulation runs but I notice an abnormal temperature rise around the injection point.
• When selecting None, I experience convergence difficulties, but the resulting temperature field is correct.

To clarify:

• The first image shows only the flow field without coupling, with a cold zone around the injection area (expected correct output).
• The second image shows the result with the Global type, where there is an abnormal temperature increase near the injection point.

Additionally, I noticed the following message when selecting the Local stabilization type:
***** LOCATION: /home/codes/GEOS/src/coreComponents/physicsSolvers/multiphysics/PoromechanicsSolver.hpp:144*****
Controlling expression (should be false): m_stabilizationType == stabilization::StabilizationType::Local
***** Rank 1: multiphasePoroelasticity/stabilizationType (1.xml, l.13): Local stabilization has been temporarily disabled

Could you please provide guidance or documentation on the settings for Global and None, and how to avoid the situation where one model converges but gives incorrect results, while the other gives correct results but cannot converge? Also, how should the Local stabilization be properly configured?

Start by reproducing the coupled MultiphasePoromechanics, CompositionalMultiphaseFVM, and SolidMechanics_LagrangianFEM cases with Global and None stabilization. Read PoromechanicsSolver.hpp:144 and the stabilizationType setting in 1.xml:13; done means the supported Local behavior and Global/None settings are clarified, and the temperature-versus-convergence discrepancy is explained.