From 5c9d9c514cdb1adbc1b550c7a70b80ff02f0fa4d Mon Sep 17 00:00:00 2001 From: npillardou Date: Tue, 25 Feb 2025 12:05:43 +0100 Subject: [PATCH 1/4] Add Thermal Stress calculations fix --- .../constitutive/solid/PorousSolid.hpp | 136 ++++++++++++++++-- .../ThermalMultiphasePoromechanics.hpp | 3 + .../ThermalMultiphasePoromechanics_impl.hpp | 38 ++--- .../ThermalSinglePhasePoromechanics.hpp | 3 + .../ThermalSinglePhasePoromechanics_impl.hpp | 38 ++--- .../FixedStressThermoPoromechanics_impl.hpp | 1 + 6 files changed, 173 insertions(+), 46 deletions(-) diff --git a/src/coreComponents/constitutive/solid/PorousSolid.hpp b/src/coreComponents/constitutive/solid/PorousSolid.hpp index b4a9429eb43..1bed34ff9b2 100644 --- a/src/coreComponents/constitutive/solid/PorousSolid.hpp +++ b/src/coreComponents/constitutive/solid/PorousSolid.hpp @@ -131,6 +131,69 @@ class PorousSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPorosity, dSolidDensity_dPressure = m_porosityUpdate.dGrainDensity_dPressure( k ); } + GEOS_HOST_DEVICE + void smallStrainUpdateThermoPoromechanics( localIndex const k, + localIndex const q, + real64 const & timeIncrement, + real64 const & pressure, + real64 const & pressure_n, + real64 const & temperature, + real64 const & deltaTemperatureFromLastStep, + real64 const & initialTemperature, + real64 const ( &strainIncrement )[6], + real64 ( & totalStress )[6], + real64 ( & dTotalStress_dPressure )[6], + real64 ( & dTotalStress_dTemperature )[6], + DiscretizationOps & stiffness, + integer const performStressInitialization, + real64 & porosity, + real64 & porosity_n, + real64 & dPorosity_dVolStrain, + real64 & dPorosity_dPressure, + real64 & dPorosity_dTemperature, + real64 & dSolidDensity_dPressure ) const + { + // Compute total stress increment and its derivative + computeTotalStressWithInitialTemperature( k, + q, + timeIncrement, + pressure, + temperature, + initialTemperature, + strainIncrement, + totalStress, + dTotalStress_dPressure, + dTotalStress_dTemperature, + stiffness ); + + // Compute porosity and its derivatives + real64 const deltaPressure = pressure - pressure_n; + real64 porosityInit; + computePorosity( k, + q, + deltaPressure, + deltaTemperatureFromLastStep, + strainIncrement, + porosity, + porosity_n, + porosityInit, + dPorosity_dVolStrain, + dPorosity_dPressure, + dPorosity_dTemperature ); + + // skip porosity update when doing poromechanics initialization + if( performStressInitialization ) + { + porosity = porosityInit; + dPorosity_dVolStrain = 0.0; + dPorosity_dPressure = 0.0; + dPorosity_dTemperature = 0.0; + } + + // Save the derivative of solid density wrt pressure for the computation of the body force + dSolidDensity_dPressure = m_porosityUpdate.dGrainDensity_dPressure( k ); + } + GEOS_HOST_DEVICE void smallStrainUpdatePoromechanicsFixedStress( localIndex const k, localIndex const q, @@ -139,6 +202,7 @@ class PorousSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPorosity, real64 const & pressure_n, real64 const & temperature, real64 const & temperature_n, + real64 const & initialTemperature, real64 const ( &strainIncrement )[6], real64 ( & totalStress )[6], DiscretizationOps & stiffness ) const @@ -147,16 +211,17 @@ class PorousSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPorosity, real64 dTotalStress_dTemperature[6]{}; // Compute total stress increment and its derivative - computeTotalStress( k, - q, - timeIncrement, - pressure, - temperature, - strainIncrement, - totalStress, - dTotalStress_dPressure, // To pass something here - dTotalStress_dTemperature, // To pass something here - stiffness ); + computeTotalStressWithInitialTemperature( k, + q, + timeIncrement, + pressure, + temperature, + initialTemperature, + strainIncrement, + totalStress, + dTotalStress_dPressure, // To pass something here + dTotalStress_dTemperature, // To pass something here + stiffness ); // Compute total stress increment for the porosity update real64 const bulkModulus = m_solidUpdate.getBulkModulus( k ); @@ -312,6 +377,57 @@ class PorousSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPorosity, dTotalStress_dTemperature[5] = 0; } + + GEOS_HOST_DEVICE + inline + void computeTotalStressWithInitialTemperature( localIndex const k, + localIndex const q, + real64 const & timeIncrement, + real64 const & pressure, + real64 const & temperature, + real64 const & initialTemperature, + real64 const ( &strainIncrement )[6], + real64 ( & totalStress )[6], + real64 ( & dTotalStress_dPressure )[6], + real64 ( & dTotalStress_dTemperature )[6], + DiscretizationOps & stiffness ) const + { + updateBiotCoefficientAndAssignModuli( k ); + + // Compute total stress increment and its derivative w.r.t. pressure + m_solidUpdate.smallStrainUpdate( k, + q, + timeIncrement, + strainIncrement, + totalStress, // first effective stress increment accumulated + stiffness ); + + // Add the contributions of pressure and temperature to the total stress + real64 const biotCoefficient = m_porosityUpdate.getBiotCoefficient( k ); + real64 const thermalExpansionCoefficient = m_solidUpdate.getThermalExpansionCoefficient( k ); + real64 const bulkModulus = m_solidUpdate.getBulkModulus( k ); + real64 const thermalExpansionCoefficientTimesBulkModulus = thermalExpansionCoefficient * bulkModulus; + real64 const deltaTemperature = temperature - initialTemperature; + + LvArray::tensorOps::symAddIdentity< 3 >( totalStress, -biotCoefficient * pressure - 3 * thermalExpansionCoefficientTimesBulkModulus * deltaTemperature ); + + // Compute derivatives of total stress + dTotalStress_dPressure[0] = -biotCoefficient; + dTotalStress_dPressure[1] = -biotCoefficient; + dTotalStress_dPressure[2] = -biotCoefficient; + dTotalStress_dPressure[3] = 0; + dTotalStress_dPressure[4] = 0; + dTotalStress_dPressure[5] = 0; + + dTotalStress_dTemperature[0] = -3 * thermalExpansionCoefficientTimesBulkModulus; + dTotalStress_dTemperature[1] = -3 * thermalExpansionCoefficientTimesBulkModulus; + dTotalStress_dTemperature[2] = -3 * thermalExpansionCoefficientTimesBulkModulus; + dTotalStress_dTemperature[3] = 0; + dTotalStress_dTemperature[4] = 0; + dTotalStress_dTemperature[5] = 0; + + } + }; /** diff --git a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalMultiphasePoromechanics.hpp b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalMultiphasePoromechanics.hpp index 110376c9d95..83e40f68577 100644 --- a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalMultiphasePoromechanics.hpp +++ b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalMultiphasePoromechanics.hpp @@ -326,6 +326,9 @@ class ThermalMultiphasePoromechanics : arrayView2d< real64 const > const m_rockInternalEnergy; arrayView2d< real64 const > const m_dRockInternalEnergy_dTemperature; + /// The rank-global initial temperature array + arrayView1d< real64 const > const m_initialTemperature; + /// Views on temperature arrayView1d< real64 const > const m_temperature_n; arrayView1d< real64 const > const m_temperature; diff --git a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalMultiphasePoromechanics_impl.hpp b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalMultiphasePoromechanics_impl.hpp index 98e5a37075f..3bbc7e07ddb 100644 --- a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalMultiphasePoromechanics_impl.hpp +++ b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalMultiphasePoromechanics_impl.hpp @@ -78,6 +78,7 @@ ThermalMultiphasePoromechanics( NodeManager const & nodeManager, m_rockInternalEnergy_n( inputConstitutiveType.getInternalEnergy_n() ), m_rockInternalEnergy( inputConstitutiveType.getInternalEnergy() ), m_dRockInternalEnergy_dTemperature( inputConstitutiveType.getDinternalEnergy_dTemperature() ), + m_initialTemperature( elementSubRegion.template getField< fields::flow::initialTemperature >() ), m_temperature_n( elementSubRegion.template getField< fields::flow::temperature_n >() ), m_temperature( elementSubRegion.template getField< fields::flow::temperature >() ) { @@ -129,24 +130,25 @@ smallStrainUpdate( localIndex const k, real64 dSolidDensity_dPressure = 0.0; // Step 1: call the constitutive model to update the total stress, the porosity and their derivatives - m_constitutiveUpdate.smallStrainUpdatePoromechanics( k, q, - m_dt, - m_pressure[k], - m_pressure_n[k], - stack.temperature, - stack.deltaTemperatureFromLastStep, - stack.strainIncrement, - stack.totalStress, - stack.dTotalStress_dPressure, - stack.dTotalStress_dTemperature, - stack.stiffness, - m_performStressInitialization, - porosity, - porosity_n, - dPorosity_dVolStrain, - dPorosity_dPressure, - dPorosity_dTemperature, - dSolidDensity_dPressure ); + m_constitutiveUpdate.smallStrainUpdateThermoPoromechanics( k, q, + m_dt, + m_pressure[k], + m_pressure_n[k], + stack.temperature, + stack.deltaTemperatureFromLastStep, + m_initialTemperature[k], + stack.strainIncrement, + stack.totalStress, + stack.dTotalStress_dPressure, + stack.dTotalStress_dTemperature, + stack.stiffness, + m_performStressInitialization, + porosity, + porosity_n, + dPorosity_dVolStrain, + dPorosity_dPressure, + dPorosity_dTemperature, + dSolidDensity_dPressure ); // Step 2: compute the body force computeBodyForce( k, q, diff --git a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanics.hpp b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanics.hpp index a84c66dbaec..a42f90994cd 100644 --- a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanics.hpp +++ b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanics.hpp @@ -293,6 +293,9 @@ class ThermalSinglePhasePoromechanics : arrayView2d< real64 const > const m_rockInternalEnergy; arrayView2d< real64 const > const m_dRockInternalEnergy_dTemperature; + /// The rank-global initial temperature array + arrayView1d< real64 const > const m_initialTemperature; + /// Views on temperature arrayView1d< real64 const > const m_temperature_n; arrayView1d< real64 const > const m_temperature; diff --git a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanics_impl.hpp b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanics_impl.hpp index 786ad448381..6cc03970920 100644 --- a/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanics_impl.hpp +++ b/src/coreComponents/physicsSolvers/multiphysics/poromechanicsKernels/ThermalSinglePhasePoromechanics_impl.hpp @@ -76,6 +76,7 @@ ThermalSinglePhasePoromechanics( NodeManager const & nodeManager, m_rockInternalEnergy_n( inputConstitutiveType.getInternalEnergy_n() ), m_rockInternalEnergy( inputConstitutiveType.getInternalEnergy() ), m_dRockInternalEnergy_dTemperature( inputConstitutiveType.getDinternalEnergy_dTemperature() ), + m_initialTemperature( elementSubRegion.template getField< fields::flow::initialTemperature >() ), m_temperature_n( elementSubRegion.template getField< fields::flow::temperature_n >() ), m_temperature( elementSubRegion.template getField< fields::flow::temperature >() ) {} @@ -113,24 +114,25 @@ smallStrainUpdate( localIndex const k, real64 dSolidDensity_dPressure = 0.0; // Step 1: call the constitutive model to evaluate the total stress and compute porosity - m_constitutiveUpdate.smallStrainUpdatePoromechanics( k, q, - m_dt, - m_pressure[k], - m_pressure_n[k], - stack.temperature, - stack.deltaTemperatureFromLastStep, - stack.strainIncrement, - stack.totalStress, - stack.dTotalStress_dPressure, - stack.dTotalStress_dTemperature, - stack.stiffness, - m_performStressInitialization, - porosity, - porosity_n, - dPorosity_dVolStrain, - dPorosity_dPressure, - dPorosity_dTemperature, - dSolidDensity_dPressure ); + m_constitutiveUpdate.smallStrainUpdateThermoPoromechanics( k, q, + m_dt, + m_pressure[k], + m_pressure_n[k], + stack.temperature, + stack.deltaTemperatureFromLastStep, + m_initialTemperature[k], + stack.strainIncrement, + stack.totalStress, + stack.dTotalStress_dPressure, + stack.dTotalStress_dTemperature, + stack.stiffness, + m_performStressInitialization, + porosity, + porosity_n, + dPorosity_dVolStrain, + dPorosity_dPressure, + dPorosity_dTemperature, + dSolidDensity_dPressure ); // Step 2: compute the body force computeBodyForce( k, q, diff --git a/src/coreComponents/physicsSolvers/solidMechanics/kernels/FixedStressThermoPoromechanics_impl.hpp b/src/coreComponents/physicsSolvers/solidMechanics/kernels/FixedStressThermoPoromechanics_impl.hpp index 36d70558513..8d4983a6dfb 100644 --- a/src/coreComponents/physicsSolvers/solidMechanics/kernels/FixedStressThermoPoromechanics_impl.hpp +++ b/src/coreComponents/physicsSolvers/solidMechanics/kernels/FixedStressThermoPoromechanics_impl.hpp @@ -143,6 +143,7 @@ quadraturePointKernel( localIndex const k, m_pressure_n[k], m_temperature[k], m_temperature_n[k], + m_initialTemperature[k], strainInc, totalStress, stiffness ); From 56ed6754d9904f5215cd2fc5ab78c46d48e1d58c Mon Sep 17 00:00:00 2001 From: npillardou Date: Thu, 27 Feb 2025 15:16:17 +0100 Subject: [PATCH 2/4] Add Smoke tests for THM --- .../ThermoPoro_SinglePhase_fim.xml | 425 ++++++++++++++++ .../ThermoPoro_SinglePhase_sequential.xml | 455 ++++++++++++++++++ 2 files changed, 880 insertions(+) create mode 100755 inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_fim.xml create mode 100755 inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_sequential.xml diff --git a/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_fim.xml b/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_fim.xml new file mode 100755 index 00000000000..a66a1801f44 --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_fim.xml @@ -0,0 +1,425 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_sequential.xml b/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_sequential.xml new file mode 100755 index 00000000000..fb57aa56b86 --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_sequential.xml @@ -0,0 +1,455 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + From f001d8d648b3c9e8d547d1db730f7797eaf79cb0 Mon Sep 17 00:00:00 2001 From: Jian HUANG Date: Fri, 28 Feb 2025 16:35:24 -0600 Subject: [PATCH 3/4] prepare smoke and benchmark examples --- ...glePhase_ThermalDradient_benchmark_fim.xml | 176 +++++++ ...e_ThermalDradient_benchmark_sequential.xml | 206 ++++++++ ..._SinglePhase_ThermalDradient_smoke_fim.xml | 85 ++++ ...Phase_ThermalDradient_smoke_sequential.xml | 115 +++++ ...stic_SinglePhase_ThermalGradient_base.xml} | 234 +-------- .../ThermoPoro_SinglePhase_sequential.xml | 455 ------------------ 6 files changed, 608 insertions(+), 663 deletions(-) create mode 100755 inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_benchmark_fim.xml create mode 100755 inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_benchmark_sequential.xml create mode 100755 inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_smoke_fim.xml create mode 100755 inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_smoke_sequential.xml rename inputFiles/thermoPoromechanics/{ThermoPoro_SinglePhase_fim.xml => ThermoPoroElastic_SinglePhase_ThermalGradient_base.xml} (52%) delete mode 100755 inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_sequential.xml diff --git a/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_benchmark_fim.xml b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_benchmark_fim.xml new file mode 100755 index 00000000000..a738fcaf52a --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_benchmark_fim.xml @@ -0,0 +1,176 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_benchmark_sequential.xml b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_benchmark_sequential.xml new file mode 100755 index 00000000000..89e4023ed75 --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_benchmark_sequential.xml @@ -0,0 +1,206 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_smoke_fim.xml b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_smoke_fim.xml new file mode 100755 index 00000000000..fde93b3744a --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_smoke_fim.xml @@ -0,0 +1,85 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_smoke_sequential.xml b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_smoke_sequential.xml new file mode 100755 index 00000000000..a13c4e36912 --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalDradient_smoke_sequential.xml @@ -0,0 +1,115 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_fim.xml b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalGradient_base.xml similarity index 52% rename from inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_fim.xml rename to inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalGradient_base.xml index a66a1801f44..3f2e79211b2 100755 --- a/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_fim.xml +++ b/inputFiles/thermoPoromechanics/ThermoPoroElastic_SinglePhase_ThermalGradient_base.xml @@ -1,91 +1,37 @@ - - - - - - - - - - - - - - - - - - - - - + + + + + - - @@ -138,7 +83,6 @@ - - - - + + + name="sourceTerm" + objectPath="ElementRegions/Reservoir" + component="0" + scale="-8.87" + beginTime="0" + endTime="6.30721e8" + setNames="{ source }"/> - + name="sourceTemperature" + setNames="{ source }" + objectPath="ElementRegions/Reservoir" + fieldName="temperature" + scale="293.15" + beginTime="0" + endTime="6.30721e8" /> - - @@ -289,137 +229,15 @@ name="singlePhasePoroelasticityEquilibrationStep" poromechanicsSolverName="singlePhasePoroelasticity" solidMechanicsStatisticsName="linearElasticityStatistics"/> - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + diff --git a/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_sequential.xml b/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_sequential.xml deleted file mode 100755 index fb57aa56b86..00000000000 --- a/inputFiles/thermoPoromechanics/ThermoPoro_SinglePhase_sequential.xml +++ /dev/null @@ -1,455 +0,0 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - From 3f241776574fcef05c52f4186091b1fc528f76e5 Mon Sep 17 00:00:00 2001 From: Jian HUANG Date: Fri, 28 Feb 2025 19:04:15 -0600 Subject: [PATCH 4/4] update tutorial examples accordingly --- .../ThermoPoroElastic_consolidation_base.xml | 6 +++--- ...lasticWellbore_ImperfectInterfaces_base.xml | 18 +++++++++--------- .../CasedThermoElasticWellbore_base.xml | 18 +++++++++--------- ...hermoElastic_casedContactWellbore_stress.py | 6 +++--- .../thermoElastic_casedWellbore_stress.py | 6 +++--- 5 files changed, 27 insertions(+), 27 deletions(-) diff --git a/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_base.xml b/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_base.xml index 2027c17639f..2e5c714cb81 100644 --- a/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_base.xml +++ b/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_base.xml @@ -106,7 +106,7 @@ objectPath="ElementRegions/Domain/cb1" fieldName="rockSolid_stress" component="0" - scale="2.457"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> diff --git a/inputFiles/wellbore/CasedThermoElasticWellbore_ImperfectInterfaces_base.xml b/inputFiles/wellbore/CasedThermoElasticWellbore_ImperfectInterfaces_base.xml index ab4ccf65952..4fa51f0f0d4 100644 --- a/inputFiles/wellbore/CasedThermoElasticWellbore_ImperfectInterfaces_base.xml +++ b/inputFiles/wellbore/CasedThermoElasticWellbore_ImperfectInterfaces_base.xml @@ -330,7 +330,7 @@ objectPath="ElementRegions/casing/casing" fieldName="casingSolid_stress" component="0" - scale="573913043.5"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/> + scale="0.0"/>