From 805e889f1b988aad7280f1aadc5e46a66f56a1b9 Mon Sep 17 00:00:00 2001 From: Pavel Tomin Date: Thu, 12 Dec 2024 10:44:27 -0600 Subject: [PATCH] Update ThermalSinglePhasePoromechanicsEFEM_impl.hpp --- .../ThermalSinglePhasePoromechanicsEFEM_impl.hpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanicsEFEM_impl.hpp b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanicsEFEM_impl.hpp index fd1a33dce8..ac4475339d 100644 --- a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanicsEFEM_impl.hpp +++ b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanicsEFEM_impl.hpp @@ -176,7 +176,7 @@ complete( localIndex const k, // Energy balance accumulation real64 const volume = m_elementVolumeFrac( embSurfIndex ) + m_deltaVolume( embSurfIndex ); real64 const fluidEnergy = m_fluidDensity( embSurfIndex, 0 ) * m_fluidInternalEnergy( embSurfIndex, 0 ) * volume; - real64 const fluidEnergy_n = m_energy_n[embSurfIndex]; // TODO where is solid energy? + real64 const fluidEnergy_n = m_energy_n[embSurfIndex]; stack.dFluidMassIncrement_dTemperature = m_dFluidDensity_dTemperature( embSurfIndex, 0 ) * volume;