From c11a5d7da8781b0d1f120f0c5e14a6ba0dd60c72 Mon Sep 17 00:00:00 2001 From: dkachuma Date: Fri, 30 Jan 2026 16:05:16 -0600 Subject: [PATCH] fix: Fix bug in L2 norm calculation for isothermal --- .../physicsSolvers/fluidFlow/CompositionalMultiphaseFVM.cpp | 2 ++ 1 file changed, 2 insertions(+) diff --git a/src/coreComponents/physicsSolvers/fluidFlow/CompositionalMultiphaseFVM.cpp b/src/coreComponents/physicsSolvers/fluidFlow/CompositionalMultiphaseFVM.cpp index f4a8d3c3ffd..32b883624b3 100644 --- a/src/coreComponents/physicsSolvers/fluidFlow/CompositionalMultiphaseFVM.cpp +++ b/src/coreComponents/physicsSolvers/fluidFlow/CompositionalMultiphaseFVM.cpp @@ -538,6 +538,8 @@ real64 CompositionalMultiphaseFVM::calculateResidualNorm( real64 const & GEOS_UN // volume subRegionResidualNorm[1] = subRegionFlowResidualNorm[1]; subRegionResidualNormalizer[1] = subRegionFlowResidualNormalizer[1]; + // Ensure that the normalizer for the energy equation is non-zero + subRegionResidualNormalizer[2] = m_nonlinearSolverParameters.m_minNormalizer; } // step 2: first reduction across meshBodies/regions/subRegions