diff --git a/.integrated_tests.yaml b/.integrated_tests.yaml index b24f3a84cfc..78819e860b8 100644 --- a/.integrated_tests.yaml +++ b/.integrated_tests.yaml @@ -1,6 +1,6 @@ baselines: bucket: geosx - baseline: integratedTests/baseline_integratedTests-pr4114-17283-3bb33cb + baseline: integratedTests/baseline_integratedTests-pr3884-17320-39debdc allow_fail: all: '' diff --git a/BASELINE_NOTES.md b/BASELINE_NOTES.md index bdc04d96954..f01cc4a708d 100644 --- a/BASELINE_NOTES.md +++ b/BASELINE_NOTES.md @@ -5,6 +5,9 @@ This file is designed to track changes to the integrated test baselines. Any developer who updates the baseline ID in the .integrated_tests.yaml file is expected to create an entry in this file with the pull request number, date, and their justification for rebaselining. These notes should be in reverse-chronological order, and use the following time format: (YYYY-MM-DD). +PR #3884 (2026-08-16) +Total stress fix in the thermo-poromechanics model + PR #4114 (2026-08-14) Stop dumping linear systems from ATS decks (`writeLinearSystem` no longer set). That flag is stored in restart files, so `perf_status_test` restartchecks need a new baseline. diff --git a/inputFiles/thermoPoromechanics/ThermoDruckerPrager_1DCooling_fim_smoke.xml b/inputFiles/thermoPoromechanics/ThermoDruckerPrager_1DCooling_fim_smoke.xml new file mode 100644 index 00000000000..fa6d45549b0 --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoDruckerPrager_1DCooling_fim_smoke.xml @@ -0,0 +1,159 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoElastic_1DCooling_fim_smoke.xml b/inputFiles/thermoPoromechanics/ThermoElastic_1DCooling_fim_smoke.xml new file mode 100644 index 00000000000..d03c8003a39 --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoElastic_1DCooling_fim_smoke.xml @@ -0,0 +1,155 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoElastic_1DCooling_sq_smoke.xml b/inputFiles/thermoPoromechanics/ThermoElastic_1DCooling_sq_smoke.xml new file mode 100644 index 00000000000..6c56d1f51cf --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoElastic_1DCooling_sq_smoke.xml @@ -0,0 +1,163 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoMech_1DCooling_base.xml b/inputFiles/thermoPoromechanics/ThermoMech_1DCooling_base.xml new file mode 100644 index 00000000000..858d025e01b --- /dev/null +++ b/inputFiles/thermoPoromechanics/ThermoMech_1DCooling_base.xml @@ -0,0 +1,116 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_base.xml b/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_base.xml index 2e5c714cb81..af59a4c1efd 100644 --- a/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_base.xml +++ b/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_base.xml @@ -196,7 +196,12 @@ name="displacementCollection" objectPath="nodeManager" fieldName="totalDisplacement"/> - + + + @@ -216,7 +221,12 @@ name="displacementHistoryOutput" sources="{ /Tasks/displacementCollection }" filename="displacementHistory"/> - + + + diff --git a/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_benchmark_base.xml b/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_benchmark_base.xml index a938018d6d9..d631e13dbc6 100644 --- a/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_benchmark_base.xml +++ b/inputFiles/thermoPoromechanics/ThermoPoroElastic_consolidation_benchmark_base.xml @@ -85,6 +85,16 @@ endTime="1e-2" forceDt="1e-2" target="/Outputs/displacementHistoryOutput"/> + + - + + + + + + + + + + + + + + + - + + + - + + + smallStrainUpdate( k, q, timeIncrement, strainIncrement, stress, stiffness.m_c ); } - GEOS_HOST_DEVICE virtual real64 getStrainEnergyDensity( localIndex const k, localIndex const q ) const override final diff --git a/src/coreComponents/constitutive/solid/DamageVolDev.hpp b/src/coreComponents/constitutive/solid/DamageVolDev.hpp index 951eded451f..cad75993cbd 100644 --- a/src/coreComponents/constitutive/solid/DamageVolDev.hpp +++ b/src/coreComponents/constitutive/solid/DamageVolDev.hpp @@ -153,7 +153,6 @@ class DamageVolDevUpdates : public DamageUpdates< UPDATE_BASE > } } - GEOS_HOST_DEVICE virtual real64 getStrainEnergyDensity( localIndex const k, localIndex const q ) const override final diff --git a/src/coreComponents/constitutive/solid/DelftEgg.hpp b/src/coreComponents/constitutive/solid/DelftEgg.hpp index 167301468f0..5d9afab89f4 100644 --- a/src/coreComponents/constitutive/solid/DelftEgg.hpp +++ b/src/coreComponents/constitutive/solid/DelftEgg.hpp @@ -449,8 +449,6 @@ void DelftEggUpdates::smallStrainUpdate( localIndex const k, smallStrainUpdate( k, q, timeIncrement, strainIncrement, stress, stiffness.m_c ); } - - /** * @class DelftEgg * diff --git a/src/coreComponents/constitutive/solid/DruckerPrager.hpp b/src/coreComponents/constitutive/solid/DruckerPrager.hpp index e59494d7b06..02e2a44fcf7 100644 --- a/src/coreComponents/constitutive/solid/DruckerPrager.hpp +++ b/src/coreComponents/constitutive/solid/DruckerPrager.hpp @@ -339,8 +339,6 @@ void DruckerPragerUpdates::smallStrainUpdate( localIndex const k, smallStrainUpdate( k, q, timeIncrement, strainIncrement, stress, stiffness.m_c ); } - - /** * @class DruckerPrager * diff --git a/src/coreComponents/constitutive/solid/DruckerPragerExtended.hpp b/src/coreComponents/constitutive/solid/DruckerPragerExtended.hpp index dc79394677f..ebe6e8d14b9 100644 --- a/src/coreComponents/constitutive/solid/DruckerPragerExtended.hpp +++ b/src/coreComponents/constitutive/solid/DruckerPragerExtended.hpp @@ -116,7 +116,6 @@ class DruckerPragerExtendedUpdates : public ElasticIsotropicUpdates real64 ( &stress )[6], real64 ( &stiffness )[6][6] ) const override; - GEOS_HOST_DEVICE inline virtual void saveConvergedState( localIndex const k, @@ -365,8 +364,6 @@ void DruckerPragerExtendedUpdates::smallStrainUpdate( localIndex const k, smallStrainUpdate( k, q, timeIncrement, strainIncrement, stress, stiffness.m_c ); } - - /** * @class DruckerPragerExtended * diff --git a/src/coreComponents/constitutive/solid/PorousDamageSolid.hpp b/src/coreComponents/constitutive/solid/PorousDamageSolid.hpp index d32923f83de..99e85bc4223 100644 --- a/src/coreComponents/constitutive/solid/PorousDamageSolid.hpp +++ b/src/coreComponents/constitutive/solid/PorousDamageSolid.hpp @@ -59,7 +59,7 @@ class PorousDamageSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPor real64 const & timeIncrement, real64 const & pressure_n, real64 const & pressure, - real64 const & temperature, + real64 const & deltaTemperature, real64 const & deltaTemperatureFromLastStep, real64 const ( &strainIncrement )[6], real64 ( & totalStress )[6], @@ -80,7 +80,8 @@ class PorousDamageSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPor timeIncrement, pressure_n, pressure, - temperature, + deltaTemperature, + deltaTemperatureFromLastStep, strainIncrement, totalStress, dTotalStress_dPressure, @@ -266,33 +267,45 @@ class PorousDamageSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPor real64 const & timeIncrement, real64 const & pressure_n, real64 const & pressure, - real64 const & temperature, + real64 const & deltaTemperature, + real64 const & deltaTemperatureFromLastStep, real64 const ( &strainIncrement )[6], real64 ( & totalStress )[6], real64 ( & dTotalStress_dPressure )[6], real64 ( & dTotalStress_dTemperature )[6], DiscretizationOps & stiffness ) const { + GEOS_UNUSED_VAR( deltaTemperature ); + updateBiotCoefficientAndAssignModuli( k ); + real64 const thermalExpansionCoefficient = m_solidUpdate.getThermalExpansionCoefficient( k ); + + real64 strainIncrementNoThermalStrain[6]{}; + strainIncrementNoThermalStrain[0] = strainIncrement[0] - thermalExpansionCoefficient * deltaTemperatureFromLastStep; + strainIncrementNoThermalStrain[1] = strainIncrement[1] - thermalExpansionCoefficient * deltaTemperatureFromLastStep; + strainIncrementNoThermalStrain[2] = strainIncrement[2] - thermalExpansionCoefficient * deltaTemperatureFromLastStep; + strainIncrementNoThermalStrain[3] = strainIncrement[3]; + strainIncrementNoThermalStrain[4] = strainIncrement[4]; + strainIncrementNoThermalStrain[5] = strainIncrement[5]; + // Compute total stress increment and its derivative w.r.t. pressure m_solidUpdate.smallStrainUpdate( k, q, timeIncrement, - strainIncrement, - totalStress, // first effective stress increment accumulated + strainIncrementNoThermalStrain, + totalStress, stiffness ); + stiffness.computeTemperatureDerivative( thermalExpansionCoefficient, dTotalStress_dTemperature ); + // 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 pressureDamage = m_solidUpdate.pressureDamageFunction( k, q ); real64 const damagedBiotCoefficient = pressureDamage * biotCoefficient; - LvArray::tensorOps::symAddIdentity< 3 >( totalStress, pressureDamage * pressure_n - damagedBiotCoefficient * pressure - 3 * thermalExpansionCoefficientTimesBulkModulus * temperature ); + LvArray::tensorOps::symAddIdentity< 3 >( totalStress, pressureDamage * pressure_n - damagedBiotCoefficient * pressure ); dTotalStress_dPressure[0] = -damagedBiotCoefficient; dTotalStress_dPressure[1] = -damagedBiotCoefficient; @@ -300,13 +313,6 @@ class PorousDamageSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPor 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/constitutive/solid/PorousSolid.hpp b/src/coreComponents/constitutive/solid/PorousSolid.hpp index d7f27433787..b586181b420 100644 --- a/src/coreComponents/constitutive/solid/PorousSolid.hpp +++ b/src/coreComponents/constitutive/solid/PorousSolid.hpp @@ -98,6 +98,7 @@ class PorousSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPorosity, timeIncrement, pressure, deltaTemperature, + deltaTemperatureFromLastStep, strainIncrement, totalStress, dTotalStress_dPressure, @@ -150,11 +151,13 @@ class PorousSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPorosity, // Compute total stress increment and its derivative real64 const deltaTemperature = temperature - referenceTemperature; + real64 const deltaTemperatureFromLastStep = temperature - temperature_n; computeTotalStress( k, q, timeIncrement, pressure, deltaTemperature, + deltaTemperatureFromLastStep, strainIncrement, totalStress, dTotalStress_dPressure, // To pass something here @@ -275,29 +278,41 @@ class PorousSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPorosity, real64 const & timeIncrement, real64 const & pressure, real64 const & deltaTemperature, + real64 const & deltaTemperatureFromLastStep, real64 const ( &strainIncrement )[6], real64 ( & totalStress )[6], real64 ( & dTotalStress_dPressure )[6], real64 ( & dTotalStress_dTemperature )[6], DiscretizationOps & stiffness ) const { + GEOS_UNUSED_VAR( deltaTemperature ); + updateBiotCoefficientAndAssignModuli( k ); - // Compute total stress increment and its derivative w.r.t. pressure + real64 const thermalExpansionCoefficient = m_solidUpdate.getThermalExpansionCoefficient( k ); + + real64 strainIncrementNoThermalStrain[6]{}; + strainIncrementNoThermalStrain[0] = strainIncrement[0] - thermalExpansionCoefficient * deltaTemperatureFromLastStep; + strainIncrementNoThermalStrain[1] = strainIncrement[1] - thermalExpansionCoefficient * deltaTemperatureFromLastStep; + strainIncrementNoThermalStrain[2] = strainIncrement[2] - thermalExpansionCoefficient * deltaTemperatureFromLastStep; + strainIncrementNoThermalStrain[3] = strainIncrement[3]; + strainIncrementNoThermalStrain[4] = strainIncrement[4]; + strainIncrementNoThermalStrain[5] = strainIncrement[5]; + + // Compute total stress increment and its derivative w.r.t. pressure and temperature m_solidUpdate.smallStrainUpdate( k, q, timeIncrement, - strainIncrement, - totalStress, // first effective stress increment accumulated + strainIncrementNoThermalStrain, + totalStress, stiffness ); + stiffness.computeTemperatureDerivative( thermalExpansionCoefficient, dTotalStress_dTemperature ); + // 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; - LvArray::tensorOps::symAddIdentity< 3 >( totalStress, -biotCoefficient * pressure - 3 * thermalExpansionCoefficientTimesBulkModulus * deltaTemperature ); + LvArray::tensorOps::symAddIdentity< 3 >( totalStress, -biotCoefficient * pressure ); // Compute derivatives of total stress dTotalStress_dPressure[0] = -biotCoefficient; @@ -307,13 +322,6 @@ class PorousSolidUpdates : public CoupledSolidUpdates< SOLID_TYPE, BiotPorosity, 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/constitutive/solid/SolidModelDiscretizationOpsFullyAnisotropic.hpp b/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsFullyAnisotropic.hpp index 19a7db1bea3..101dd175048 100644 --- a/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsFullyAnisotropic.hpp +++ b/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsFullyAnisotropic.hpp @@ -65,6 +65,15 @@ struct SolidModelDiscretizationOpsFullyAnisotropic : public SolidModelDiscretiza LvArray::tensorOps::scale< 6, 6 >( m_c, scale ); } + GEOS_HOST_DEVICE + inline + void computeTemperatureDerivative( real64 const thermalExpansionCoefficient, real64 ( & dStress_dT )[6] ) const + { + for( int i = 0; i < 6; ++i ) + { + dStress_dT[i] = -thermalExpansionCoefficient * ( m_c[i][0] + m_c[i][1] + m_c[i][2] ); + } + } real64 m_c[6][6]; }; diff --git a/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsIsotropic.hpp b/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsIsotropic.hpp index 78172ee75f9..68a4922e1d7 100644 --- a/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsIsotropic.hpp +++ b/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsIsotropic.hpp @@ -103,6 +103,19 @@ struct SolidModelDiscretizationOpsIsotropic : public SolidModelDiscretizationOps m_shearModulus *= scale; } + GEOS_HOST_DEVICE + inline + void computeTemperatureDerivative( real64 const thermalExpansionCoefficient, real64 ( & dStress_dT )[6] ) const + { + real64 const val = -3.0 * thermalExpansionCoefficient * m_bulkModulus; + dStress_dT[0] = val; + dStress_dT[1] = val; + dStress_dT[2] = val; + dStress_dT[3] = 0; + dStress_dT[4] = 0; + dStress_dT[5] = 0; + } + real64 m_bulkModulus; ///< Bulk modulus real64 m_shearModulus; ///< Shear modulus }; diff --git a/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsOrthotropic.hpp b/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsOrthotropic.hpp index a16fdb8fb23..859fc8fc3e9 100644 --- a/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsOrthotropic.hpp +++ b/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsOrthotropic.hpp @@ -72,6 +72,18 @@ struct SolidModelDiscretizationOpsOrthotropic : public SolidModelDiscretizationO m_c66 *= scale; } + GEOS_HOST_DEVICE + inline + void computeTemperatureDerivative( real64 const thermalExpansionCoefficient, real64 ( & dStress_dT )[6] ) const + { + dStress_dT[0] = -thermalExpansionCoefficient * ( m_c11 + m_c12 + m_c13 ); + dStress_dT[1] = -thermalExpansionCoefficient * ( m_c12 + m_c22 + m_c23 ); + dStress_dT[2] = -thermalExpansionCoefficient * ( m_c13 + m_c23 + m_c33 ); + dStress_dT[3] = 0; + dStress_dT[4] = 0; + dStress_dT[5] = 0; + } + real64 m_c11; real64 m_c12; real64 m_c13; diff --git a/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsTransverseIsotropic.hpp b/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsTransverseIsotropic.hpp index 9cf78b133de..bfcb13b2202 100644 --- a/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsTransverseIsotropic.hpp +++ b/src/coreComponents/constitutive/solid/SolidModelDiscretizationOpsTransverseIsotropic.hpp @@ -74,6 +74,20 @@ struct SolidModelDiscretizationOpsTransverseIsotropic : public SolidModelDiscret m_c66 *= scale; } + // dStress_dT[i] = -alpha * sum_{j=0}^{2} C[i][j], where C12 = C11 - 2*C66 + // Shear rows are zero because TI symmetry decouples normal and shear components. + GEOS_HOST_DEVICE + inline + void computeTemperatureDerivative( real64 const thermalExpansionCoefficient, real64 ( & dStress_dT )[6] ) const + { + dStress_dT[0] = -thermalExpansionCoefficient * ( 2 * m_c11 - 2 * m_c66 + m_c13 ); + dStress_dT[1] = dStress_dT[0]; + dStress_dT[2] = -thermalExpansionCoefficient * ( 2 * m_c13 + m_c33 ); + dStress_dT[3] = 0; + dStress_dT[4] = 0; + dStress_dT[5] = 0; + } + real64 m_c11; ///< (1,1) element of TI stiffness matrix real64 m_c13; ///< (1,3) element of TI stiffness matrix real64 m_c33; ///< (3,3) element of TI stiffness matrix diff --git a/src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsLagrangianFEM.cpp b/src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsLagrangianFEM.cpp index 2e2ab2d365d..4072b5683a9 100644 --- a/src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsLagrangianFEM.cpp +++ b/src/coreComponents/physicsSolvers/solidMechanics/SolidMechanicsLagrangianFEM.cpp @@ -50,6 +50,7 @@ #include "physicsSolvers/solidMechanics/kernels/SolidMechanicsKernelsDispatchTypeList.hpp" #include "physicsSolvers/solidMechanics/kernels/SolidMechanicsFixedStressThermoPoromechanicsKernelsDispatchTypeList.hpp" #include "physicsSolvers/fluidFlow/FlowSolverBase.hpp" +#include "physicsSolvers/fluidFlow/FlowSolverBaseFields.hpp" namespace geos { @@ -1002,6 +1003,16 @@ void SolidMechanicsLagrangianFEM::implicitStepComplete( real64 const & GEOS_UNUS solidMechanics::arrayView2dLayoutStrain avgPlasticStrain = subRegion.getField< solidMechanics::averagePlasticStrain >(); solidMechanics::arrayView2dLayoutAvgStress avgStress = subRegion.getField< solidMechanics::averageStress >(); + arrayView1d< real64 const > const temperature = + subRegion.hasField< fields::flow::temperature >() + ? subRegion.getField< fields::flow::temperature >().toViewConst() + : arrayView1d< real64 const >{}; + + arrayView1d< real64 const > const temperature_n = + subRegion.hasField< fields::flow::temperature_n >() + ? subRegion.getField< fields::flow::temperature_n >().toViewConst() + : arrayView1d< real64 const >{}; + constitutive::ConstitutivePassThru< SolidBase >::execute( constitutiveRelation, [&] ( auto & solidModel ) { using SOLID_TYPE = TYPEOFREF( solidModel ); @@ -1021,7 +1032,9 @@ void SolidMechanicsLagrangianFEM::implicitStepComplete( real64 const & GEOS_UNUS avgStrain, avgPlasticStrain, stress, - avgStress ); + avgStress, + temperature, + temperature_n ); } ); diff --git a/src/coreComponents/physicsSolvers/solidMechanics/kernels/StressStrainAverageKernels.hpp b/src/coreComponents/physicsSolvers/solidMechanics/kernels/StressStrainAverageKernels.hpp index 9765e7a262d..216ee15ee7f 100644 --- a/src/coreComponents/physicsSolvers/solidMechanics/kernels/StressStrainAverageKernels.hpp +++ b/src/coreComponents/physicsSolvers/solidMechanics/kernels/StressStrainAverageKernels.hpp @@ -28,6 +28,7 @@ #include "mesh/CellElementSubRegion.hpp" #include "mesh/utilities/AverageOverQuadraturePointsKernel.hpp" #include "physicsSolvers/solidMechanics/SolidMechanicsFields.hpp" +#include "physicsSolvers/fluidFlow/FlowSolverBaseFields.hpp" namespace geos { @@ -66,6 +67,8 @@ class AverageStressStrainOverQuadraturePoints : * @param avgStrain the strain averaged over quadrature points * @param stress the stress solution field * @param avgStress the stress averaged over quadrature points + * @param temperature the temperature field (may be empty for non-thermal simulations) + * @param temperature_n the temperature at the previous time step (may be empty for non-thermal simulations) */ AverageStressStrainOverQuadraturePoints( NodeManager & nodeManager, EdgeManager const & edgeManager, @@ -78,7 +81,9 @@ class AverageStressStrainOverQuadraturePoints : fields::solidMechanics::arrayView2dLayoutStrain const avgStrain, fields::solidMechanics::arrayView2dLayoutStrain const avgPlasticStrain, arrayView3d< real64 const, solid::STRESS_USD > const stress, - fields::solidMechanics::arrayView2dLayoutAvgStress const avgStress ): + fields::solidMechanics::arrayView2dLayoutAvgStress const avgStress, + arrayView1d< real64 const > const temperature, + arrayView1d< real64 const > const temperature_n ): Base( nodeManager, edgeManager, faceManager, @@ -90,7 +95,9 @@ class AverageStressStrainOverQuadraturePoints : m_avgStrain( avgStrain ), m_avgPlasticStrain( avgPlasticStrain ), m_stress( stress ), - m_avgStress( avgStress ) + m_avgStress( avgStress ), + m_temperature( temperature ), + m_temperature_n( temperature_n ) {} /** @@ -151,6 +158,11 @@ class AverageStressStrainOverQuadraturePoints : real64 elasticStrainInc[6] = {0.0}; m_solidUpdate.getElasticStrainInc( k, q, elasticStrainInc ); + real64 const thermalExpansionCoefficient = m_solidUpdate.getThermalExpansionCoefficient( k ); + real64 const deltaTemperature = ( m_temperature.size() > 0 ) + ? ( m_temperature[k] - m_temperature_n[k] ) + : 0.0; + real64 conversionFactor[6] = {1.0, 1.0, 1.0, 0.5, 0.5, 0.5}; // used for converting from engineering shear to tensor shear for( int icomp = 0; icomp < 6; ++icomp ) @@ -158,12 +170,17 @@ class AverageStressStrainOverQuadraturePoints : m_avgStrain[k][icomp] += conversionFactor[icomp]*detJxW*strain[icomp]/m_elementVolume[k]; m_avgStress[k][icomp] += detJxW*m_stress[k][q][icomp]/m_elementVolume[k]; + // Thermal strain is purely volumetric: subtract thermal strain from normal components so that + // only the mechanical (plastic) part is accumulated. + real64 const thermalStrainInc = ( icomp < 3 ) ? thermalExpansionCoefficient * deltaTemperature : 0.0; + real64 const mechanicalStrainInc = strainInc[icomp] - thermalStrainInc; + // This is a hack to handle boundary conditions such as those seen in plane-strain wellbore problems // Essentially, if bcs are constraining the strain (and thus total displacement), we do not accumulate any plastic strain (regardless // of stresses in material law) - if( std::abs( strainInc[icomp] ) > 1.0e-8 ) + if( std::abs( mechanicalStrainInc ) > 1.0e-8 ) { - m_avgPlasticStrain[k][icomp] += conversionFactor[icomp]*detJxW*(strainInc[icomp] - elasticStrainInc[icomp])/m_elementVolume[k]; + m_avgPlasticStrain[k][icomp] += conversionFactor[icomp]*detJxW*(mechanicalStrainInc - elasticStrainInc[icomp])/m_elementVolume[k]; } } } @@ -217,6 +234,12 @@ class AverageStressStrainOverQuadraturePoints : /// The average stress fields::solidMechanics::arrayView2dLayoutAvgStress const m_avgStress; + /// The temperature at the current time step (empty for non-thermal simulations) + arrayView1d< real64 const > const m_temperature; + + /// The temperature at the previous time step (empty for non-thermal simulations) + arrayView1d< real64 const > const m_temperature_n; + }; @@ -258,11 +281,14 @@ class AverageStressStrainOverQuadraturePointsKernelFactory fields::solidMechanics::arrayView2dLayoutStrain const avgStrain, fields::solidMechanics::arrayView2dLayoutStrain const avgPlasticStrain, arrayView3d< real64 const, solid::STRESS_USD > const stress, - fields::solidMechanics::arrayView2dLayoutAvgStress const avgStress ) + fields::solidMechanics::arrayView2dLayoutAvgStress const avgStress, + arrayView1d< real64 const > const temperature = {}, + arrayView1d< real64 const > const temperature_n = {} ) { AverageStressStrainOverQuadraturePoints< FE_TYPE, SOLID_TYPE > kernel( nodeManager, edgeManager, faceManager, elementSubRegion, finiteElementSpace, - solidModel, displacement, displacementInc, avgStrain, avgPlasticStrain, stress, avgStress ); + solidModel, displacement, displacementInc, avgStrain, avgPlasticStrain, stress, avgStress, + temperature, temperature_n ); AverageStressStrainOverQuadraturePoints< FE_TYPE, SOLID_TYPE >::template kernelLaunch< POLICY >( elementSubRegion.size(), kernel ); diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/Example.rst b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/Example.rst index 32461b48560..554958a7160 100644 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/Example.rst +++ b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/Example.rst @@ -187,19 +187,30 @@ On the top surface, we impose the traction boundary condition and the non-isothe Inspecting results --------------------------------- -We request an output of the displacements, pressure, and temperature using the **TimeHistory** feature of GEOS. -The figures below compare the results from GEOS (dashed line) and the corresponding analytical solution (solid line) as a function of time at different locations of the slab. -We obtain a very good match, confirming that GEOS can accurately capture the thermo-poromechanical coupling on this example. The first figure illustrates this good agreement for the pressure evolution. - -.. plot:: docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationPressureFigure.py - -The second figure confirms the good match with the analytical solution for the temperature. - -.. plot:: docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationTemperatureFigure.py - -The third figure shows that GEOS is also able to match the vertical displacement (settlement) analytical solution. - -.. plot:: docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationDisplacementFigure.py + We request an output of the displacements, pressure, temperature, and stress using the + **TimeHistory** feature of GEOS. + The figure below compares the results from GEOS (open circles) with the corresponding + analytical solution (solid lines) as a function of time at different locations along the slab. + We obtain a very good match for all four quantities, confirming that GEOS accurately captures + the thermo-poromechanical coupling on this example. + + - The top-left panel shows the excellent agreement for the **pore pressure** evolution: + the initial thermally- and mechanically-induced overpressure gradually dissipates as the + fluid drains through the top surface. + - The top-right panel confirms the good match for the **temperature**, which diffuses from + the heated top boundary into the column. + - The bottom-left panel shows that GEOS also reproduces the **vertical displacement** + (settlement), including the characteristic non-monotonic response (thermal expansion + followed by consolidation). + - The bottom-right panel compares the **total stress**. By equilibrium, the total vertical + stress remains equal to the applied surface load at all times and depths + (:math:`\sigma_{yy} = -F`), which GEOS recovers exactly. The total horizontal stress + :math:`\sigma_{xx} = \sigma_{zz}` follows the uniaxial-strain thermo-poroelastic relation + :math:`\sigma_{xx} = -(\lambda/M)\,F - (2G/M)\,(3K\alpha_s\,\Delta T + b\,p)` + (with :math:`M = K + 4G/3` the oedometric modulus), exhibiting the transient overshoot + driven by the coupled pressure and thermal fields. + + .. plot:: docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/plotTests_tutorial_THM.py ------------------------------------------------------------------ To go further diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/displacementHistory.hdf5 b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/displacementHistory.hdf5 deleted file mode 100644 index 7c6f53df872..00000000000 Binary files a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/displacementHistory.hdf5 and /dev/null differ diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/plotTests_tutorial_THM.py b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/plotTests_tutorial_THM.py new file mode 100755 index 00000000000..0b78369b6c7 --- /dev/null +++ b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/plotTests_tutorial_THM.py @@ -0,0 +1,100 @@ +"""Figure de verification du cas de consolidation thermo-poro-elastique. + +Execute par Sphinx (directive `.. plot::` de Example.rst) a chaque build de la documentation. +Il ne lit que deux fichiers CSV versionnes a cote de lui : + + thermalConsolidation_geos.csv resultats GEOS + thermalConsolidation_analytical.csv solution analytique + +Ces CSV sont produits hors ligne par postprocessTests_tutorial_THM.py, a partir des sorties +TimeHistory (.hdf5) d'un run du cas benchmark. Les .hdf5 ne sont volontairement pas versionnes +(plusieurs Mo) ; ce script n'a donc besoin ni de h5py ni de mpmath, seulement de numpy et +matplotlib, ce qui garde le build de la doc leger et reproductible. + +Pour regenerer les donnees apres une modification du cas : + geosx -i ThermoPoroElastic_consolidation_benchmark_fim.xml + python3 postprocessTests_tutorial_THM.py -o +""" +import os + +import numpy as np +import matplotlib.pyplot as plt + +try: + HERE = os.path.dirname(os.path.abspath(__file__)) +except NameError: + HERE = os.getcwd() + +# sous-echantillonnage des marqueurs GEOS (1 = tous les instants collectes) +N1 = 4 +fsize = 32 +msize = 8 +lw = 4 +mew = 2 +malpha = 0.6 +lalpha = 0.8 + + +def _load(name): + path = os.path.join(HERE, name) + if not os.path.exists(path): + raise FileNotFoundError( + f"{name} introuvable dans {HERE}. Regenerez-le avec postprocessTests_tutorial_THM.py " + f"(voir l'en-tete de ce fichier).") + return np.genfromtxt(path, delimiter=",", names=True) + + +def main(): + geos = _load("thermalConsolidation_geos.csv") + ana = _load("thermalConsolidation_analytical.csv") + + tg = geos["time"][::N1] + ta = ana["time"] + + fig, ax = plt.subplots(2, 2, figsize=(32, 18)) + cmap = plt.get_cmap("tab10") + + # (axe, colonnes analytiques, colonnes GEOS, labels, titre y) + panels = [ + (ax[0, 0], ["p_zh0p0", "p_zh0p6", "p_zh0p8"], ["p_0m", "p_4p2m", "p_5p6m"], + ["z/h = 0.0", "z/h = 0.6", "z/h = 0.8"], ["z = 0.0 m", "z = 4.2 m", "z = 5.6 m"], + "Pore Pressure [Pa]", "lower left"), + (ax[0, 1], ["T_zh0p0", "T_zh0p6", "T_zh0p8"], ["T_0m", "T_4p2m", "T_5p6m"], + ["z/h = 0.0", "z/h = 0.6", "z/h = 0.8"], ["z = 0.0 m", "z = 4.2 m", "z = 5.6 m"], + "Temperature [K]", "upper left"), + (ax[1, 0], ["disp_1p4m", "disp_4p2m", "disp_7m"], ["disp_1p4m", "disp_4p2m", "disp_7m"], + ["z = 1.4 m", "z = 4.2 m", "z = 7.0 m"], ["z = 1.4 m", "z = 4.2 m", "z = 7.0 m"], + "Displacement [m]", "upper left"), + (ax[1, 1], ["sxx_0m", "sxx_4p2m", "sxx_5p6m"], ["sxx_0m", "sxx_4p2m", "sxx_5p6m"], + ["$\\sigma_{xx}$: z = 0.0 m", "$\\sigma_{xx}$: z = 4.2 m", "$\\sigma_{xx}$: z = 5.6 m"], + ["$\\sigma_{xx}$: z = 0.0 m", "$\\sigma_{xx}$: z = 4.2 m", "$\\sigma_{xx}$: z = 5.6 m"], + "Total stress [Pa]", "upper left"), + ] + + for a, acols, gcols, alabels, glabels, ylabel, legloc in panels: + for k in range(3): + a.semilogx(ta, ana[acols[k]], color=cmap(k), lw=lw, alpha=lalpha, + label="Analytical: " + alabels[k]) + a.plot(tg, geos[gcols[k]][::N1], 'o', color=cmap(k), markersize=msize, alpha=malpha, + mec=cmap(k), fillstyle='none', mew=mew, label="GEOS: " + glabels[k]) + a.set_xlabel('Time [s]', size=fsize, weight="bold") + a.set_ylabel(ylabel, size=fsize, weight="bold") + a.grid(True) + a.xaxis.set_tick_params(labelsize=fsize) + a.yaxis.set_tick_params(labelsize=fsize) + + # contrainte verticale totale : constante, egale a la charge appliquee + ax[1, 1].semilogx(ta, ana["syy"], 'k--', lw=lw * 0.6, alpha=0.7, + label='$\\sigma_{yy}$ total = -F (all z)') + ax[1, 1].plot(tg, geos["syy_0m"][::N1], 'k+', markersize=msize, mew=mew, alpha=malpha, + label='GEOS $\\sigma_{yy}$ total') + + for a, *_rest in panels: + a.legend(loc=_rest[-1], fontsize=fsize * 0.5) + + plt.subplots_adjust(left=0.1, bottom=0.1, right=0.9, top=0.9, wspace=0.4, hspace=0.4) + plt.show() + + +if __name__ == "__main__": + main() diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/postprocessTests_tutorial_THM.py b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/postprocessTests_tutorial_THM.py new file mode 100644 index 00000000000..31ed04df454 --- /dev/null +++ b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/postprocessTests_tutorial_THM.py @@ -0,0 +1,235 @@ +"""Post-traitement LOCAL du cas de consolidation thermo-poro-elastique. + +A executer a la main, dans le repertoire ou GEOS a tourne (celui qui contient les .hdf5) : + + python3 postprocessTests_tutorial_THM.py [-o ] + +Il lit les sorties TimeHistory de GEOS et calcule la solution analytique, puis ecrit deux +fichiers CSV legers : + + thermalConsolidation_geos.csv resultats GEOS (une ligne par instant collecte) + thermalConsolidation_analytical.csv solution analytique (100 instants log-espaces) + +Ce sont ces deux CSV qui sont versionnes dans le depot. Les .hdf5 (plusieurs Mo) ne le sont +pas : la figure de la documentation est produite par plotTests_tutorial_THM.py, qui ne lit +que ces CSV et n'a besoin ni de h5py ni de mpmath. + +Entrees attendues (sorties du cas benchmark) : + pressureHistory.hdf5 temperatureHistory.hdf5 displacementHistory.hdf5 stressHistory.hdf5 +La sortie stress suppose que le XML active le Task 'stressCollection' et l'Output +'stressHistoryOutput' (voir ThermoPoroElastic_consolidation_base.xml). +""" +import argparse +import os + +import numpy as np +import h5py +from mpmath import nsum, exp, sin, cos, inf +import math + + +# -------------------------------------------------------------------------------------- +# Solution analytique (identique a la version historique du script de tracé) +# -------------------------------------------------------------------------------------- +def computePP(x, t, h, phi12, phi11, r11, r22): + p = nsum(lambda m: 1 / ((2 * m + 1) * math.pi / 2 / h) * + (phi12 * exp(-r11 * ((2 * m + 1) * math.pi / 2 / h)**2 * t) + + phi11 * exp(-r22 * ((2 * m + 1) * math.pi / 2 / h)**2 * t)) * + sin(((2 * m + 1) * math.pi / 2 / h) * (h - x)), [0, inf]) + return 2 / h * p + + +def computeTemp(x, t, h, phi22, phi21, r11, r22): + p = nsum(lambda m: 1 / ((2 * m + 1) * math.pi / 2 / h) * + (phi22 * exp(-r11 * ((2 * m + 1) * math.pi / 2 / h)**2 * t) + + phi21 * exp(-r22 * ((2 * m + 1) * math.pi / 2 / h)**2 * t)) * + sin(((2 * m + 1) * math.pi / 2 / h) * (h - x)), [0, inf]) + return 2 / h * p + + +def computeDisp(x, t, h, t1, t2, r11, r22, A, M): + p = nsum(lambda m: 1 / ((2 * m + 1) * math.pi / 2 / h)**2 * + (t1 * exp(-r11 * ((2 * m + 1) * math.pi / 2 / h)**2 * t) + + t2 * exp(-r22 * ((2 * m + 1) * math.pi / 2 / h)**2 * t)) * + cos(((2 * m + 1) * math.pi / 2 / h) * (h - x)), [0, inf]) + return -A * x + 2 / h / M * p + + +def analyticalParameters(): + """Parametres du probleme et coefficients de la solution analytique.""" + E = 6.0e3 + nu = 0.4 + Ks = 1.0e27 + a_s = 3.0e-7 + phi = 0.2 + + lamda = E * nu / (1.0 + nu) / (1.0 - 2.0 * nu) + G = E / 2.0 / (1.0 + nu) + K = E / 3.0 / (1.0 - 2.0 * nu) + alpha = 1.0 - K / Ks + M = lamda + 2.0 * G + + a_m = 9.0e-7 + m = 1.672e5 + beta = a_m * (3.0 * lamda + 2.0 * G) / 3.0 + perm = 4.0e-9 + mu = 1.0e-3 + km = perm / mu + h = 7.0 + theta_0 = 0.0 + theta_a = 50.0 + F = 1.0 + Kt = 836 + + nuu = 0.499 + M_biot = 2.0 * G * (nuu - nu) / (alpha**2 * (1.0 - 2.0 * nuu) * (1.0 - 2.0 * nu)) + a_p = 1.0 / M_biot + p0 = alpha * M_biot / (lamda + 2.0 * G + alpha**2 * M_biot) * F + A = (F - alpha * 0 - beta * theta_a) / M + + g11 = 1.0 / km * (a_p + alpha**2 / M) + g12 = 1.0 / km * (alpha * beta / M - a_m) + g21 = alpha * beta * theta_0 / (Kt * M) + g22 = (m + beta**2 * theta_0 / M) / Kt + gs = g11 * g22 - g12 * g21 + r11, r12 = g11 / gs, g12 / gs + r21, r22 = g21 / gs, g22 / gs + r1s = r12 * (theta_0 - theta_a) - r22 * p0 + r2s = r21 * p0 - r11 * (theta_0 - theta_a) + rs = r11 * r22 - r12 * r21 + + phi12 = -(r12 * r2s) / rs + phi11 = -(r11 * r1s) / rs + phi22 = -(r22 * r2s) / rs + phi21 = -(r21 * r1s) / rs + t1 = alpha * phi12 + beta * phi22 + t2 = alpha * phi11 + beta * phi21 + + return dict(h=h, M=M, G=G, lamda=lamda, alpha=alpha, beta=beta, F=F, theta_a=theta_a, A=A, + r11=r11, r22=r22, phi11=phi11, phi12=phi12, phi21=phi21, phi22=phi22, + t1=t1, t2=t2) + + +# -------------------------------------------------------------------------------------- +# Lecture des sorties GEOS +# -------------------------------------------------------------------------------------- +def _valid_length(t): + """Nombre d'echantillons reellement ecrits : le tampon HDF5 peut etre pre-alloue et + complete par des zeros, auquel cas le temps repasse a 0.""" + for j in range(1, len(t)): + if t[j] < 1e-12: + return j + return len(t) + + +def _pick(coords, targets): + """Premier indice dont la cote y atteint chaque cible (meme regle que le script d'origine).""" + out = [] + for c in targets: + idx = next((j for j in range(len(coords)) if coords[j] >= c), len(coords) - 1) + out.append(idx) + return out + + +def readGeos(biot=1.0): + for f in ("pressureHistory.hdf5", "temperatureHistory.hdf5", + "displacementHistory.hdf5", "stressHistory.hdf5"): + if not os.path.exists(f): + raise FileNotFoundError( + f"{f} introuvable. Lancez le cas benchmark (ThermoPoroElastic_consolidation_" + f"benchmark_fim.xml) et executez ce script dans le repertoire de sortie.") + + hp = h5py.File("pressureHistory.hdf5", 'r') + t = np.array(hp.get('pressure Time')).ravel() + n = _valid_length(t) + t = t[:n] + pc = np.array(hp.get('pressure elementCenter'))[0, :, 1] + P = np.array(hp.get('pressure'))[:n] + ip = _pick(pc, (0.0, 4.2, 5.6)) + + ht = h5py.File("temperatureHistory.hdf5", 'r') + tc = np.array(ht.get('temperature elementCenter'))[0, :, 1] + T = np.array(ht.get('temperature'))[:n] + it = _pick(tc, (0.0, 4.2, 5.6)) + + hd = h5py.File("displacementHistory.hdf5", 'r') + dc = np.array(hd.get('totalDisplacement ReferencePosition'))[0, :, 1] + U = np.array(hd.get('totalDisplacement'))[:n] + iu = _pick(dc, (1.4, 4.2, 7.0)) + + hs = h5py.File("stressHistory.hdf5", 'r') + sc = np.array(hs.get('averageStress elementCenter'))[0, :, 1] + S = np.array(hs.get('averageStress'))[:n] + isg = [int(np.argmin(np.abs(sc - d))) for d in (0.0, 4.2, 5.6)] + + # contrainte TOTALE = averageStress - biot * pression, prises dans la MEME cellule + # (maillages identiques -> l'indice de cellule stress est valable pour la pression) + sxx = [0.5 * (S[:, i, 0] + S[:, i, 2]) - biot * P[:, i] for i in isg] + syy = S[:, isg[0], 1] - biot * P[:, isg[0]] + + # positions reellement echantillonnees (utiles pour tracer des labels honnetes) + pos = dict(p=[float(pc[i]) for i in ip], T=[float(tc[i]) for i in it], + u=[float(dc[i]) for i in iu], s=[float(sc[i]) for i in isg]) + + cols = [t] + cols += [P[:, i] for i in ip] + cols += [T[:, i] for i in it] + cols += [-U[:, i, 1] for i in iu] # signe du trace : tassement compte positif + cols += list(sxx) + [syy] + return np.column_stack(cols), pos + + +def computeAnalytical(par): + h, M = par["h"], par["M"] + lamda, G, alpha, beta, F = par["lamda"], par["G"], par["alpha"], par["beta"], par["F"] + time = np.logspace(-2, 5, 100, endpoint=True) + + zp = [0.0, 0.6, 0.8] # z/h pour pression, temperature et sigma_xx + zu = [0.2, 0.6, 1.0] # z/h pour le deplacement (z = 1.4, 4.2, 7.0 m) + + P = np.empty((len(time), 3)); T = np.empty((len(time), 3)) + U = np.empty((len(time), 3)); SXX = np.empty((len(time), 3)) + for k, zc in enumerate(zp): + for i, ti in enumerate(time): + p_an = float(computePP(zc * h, ti, h, par["phi12"], par["phi11"], par["r11"], par["r22"])) + th_an = float(par["theta_a"] + computeTemp(zc * h, ti, h, par["phi22"], par["phi21"], + par["r11"], par["r22"])) + P[i, k] = p_an + T[i, k] = th_an + 273.15 + SXX[i, k] = -(lamda / M) * F - (2.0 * G / M) * (beta * th_an + alpha * p_an) + for k, zc in enumerate(zu): + for i, ti in enumerate(time): + U[i, k] = -float(computeDisp(zc * h, ti, h, par["t1"], par["t2"], + par["r11"], par["r22"], par["A"], M)) + SYY = np.full((len(time), 1), -F) + return np.column_stack([time, P, T, U, SXX, SYY]) + + +GEOS_HEADER = ("time,p_0m,p_4p2m,p_5p6m,T_0m,T_4p2m,T_5p6m," + "disp_1p4m,disp_4p2m,disp_7m,sxx_0m,sxx_4p2m,sxx_5p6m,syy_0m") +ANA_HEADER = ("time,p_zh0p0,p_zh0p6,p_zh0p8,T_zh0p0,T_zh0p6,T_zh0p8," + "disp_1p4m,disp_4p2m,disp_7m,sxx_0m,sxx_4p2m,sxx_5p6m,syy") + + +def main(): + ap = argparse.ArgumentParser(description=__doc__) + ap.add_argument("-o", "--outdir", default=".", help="dossier ou ecrire les CSV") + args = ap.parse_args() + os.makedirs(args.outdir, exist_ok=True) + + geos, pos = readGeos() + fg = os.path.join(args.outdir, "thermalConsolidation_geos.csv") + np.savetxt(fg, geos, delimiter=",", header=GEOS_HEADER, comments="", fmt="%.8e") + print(f"{fg} ({geos.shape[0]} instants)") + print(f" cotes echantillonnees pression/temperature {pos['p']} m, " + f"deplacement {pos['u']} m, contrainte {pos['s']} m") + + par = analyticalParameters() + ana = computeAnalytical(par) + fa = os.path.join(args.outdir, "thermalConsolidation_analytical.csv") + np.savetxt(fa, ana, delimiter=",", header=ANA_HEADER, comments="", fmt="%.8e") + print(f"{fa} ({ana.shape[0]} instants)") + + +if __name__ == "__main__": + main() diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/pressureHistory.hdf5 b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/pressureHistory.hdf5 deleted file mode 100644 index 11f1523fcc6..00000000000 Binary files a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/pressureHistory.hdf5 and /dev/null differ diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/temperatureHistory.hdf5 b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/temperatureHistory.hdf5 deleted file mode 100644 index 207a0237708..00000000000 Binary files a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/temperatureHistory.hdf5 and /dev/null differ diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationDisplacementFigure.py b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationDisplacementFigure.py deleted file mode 100644 index 285a7b15706..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationDisplacementFigure.py +++ /dev/null @@ -1,86 +0,0 @@ -import matplotlib -matplotlib.use('Agg') -import matplotlib.pyplot as plt -import xml.etree.ElementTree as ElementTree -import numpy as np -import h5py -from numpy import genfromtxt - -def main(): - - # File paths - hdf5FilePathDisplacement = "displacementHistory.hdf5" - - - # Read simulation output from HDF5 file - hf = h5py.File(hdf5FilePathDisplacement, 'r') - timeDisplacement = hf.get('totalDisplacement Time') - timeDisplacement = np.asarray(timeDisplacement) - centerDisplacement = hf.get('totalDisplacement ReferencePosition') - centerDisplacement = np.asarray(centerDisplacement) - displacement = hf.get('totalDisplacement') - displacement = np.asarray(displacement) - - time = 1 - posVertex1 = -1 - posVertex2 = -1 - posVertex3 = -1 - last = -1 - - for j in range(0, 284): - if centerDisplacement[0,j,1] >= 1.4 and posVertex1 == -1: - posVertex1 = j - if centerDisplacement[0,j,1] >= 4.2 and posVertex2 == -1: - posVertex2 = j - if centerDisplacement[0,j,1] >= 7 and posVertex3 == -1: - posVertex3 = j - - for j in range(0, timeDisplacement.shape[0]): - if j > 0 and timeDisplacement[j] < 1e-12 and last == -1: - last = j - - displacement_1p4 = genfromtxt('thermoConsolidationDisp_1p4m.csv', delimiter=',') - displacement_4p2 = genfromtxt('thermoConsolidationDisp_4p2m.csv', delimiter=',') - displacement_7 = genfromtxt('thermoConsolidationDisp_7m.csv', delimiter=',') - - plt.plot( timeDisplacement[0:last], - -displacement[0:last,posVertex1,1], - 'r--', - label='GEOSX: z = 1.4 m') - plt.plot( timeDisplacement[0:last], - -displacement[0:last,posVertex2,1], - 'b--', - label='GEOSX: z = 4.2 m') - plt.plot( timeDisplacement[0:last], - -displacement[0:last,posVertex3,1], - 'k--', - label='GEOSX: z = 5.6 m') - - plt.plot( displacement_1p4[:,0], - displacement_1p4[:,1], - 'r-', - label='Analytical: z = 1.4 m') - plt.plot( displacement_4p2[:,0], - displacement_4p2[:,1], - 'b-', - label='Analytical: z = 4.2 m') - plt.plot( displacement_7[:,0], - displacement_7[:,1], - 'k-', - label='Analytical: z = 7 m') - - plt.xscale("log") - - plt.grid() - plt.ylabel('Displacement [m]') - plt.xlabel('Time [s]') - plt.xlim(0.01,100000) - - plt.legend(loc='upper left') - plt.show() - # plt.savefig('displacement.png') - - - -if __name__ == "__main__": - main() diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationPressureFigure.py b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationPressureFigure.py deleted file mode 100644 index d72455e6020..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationPressureFigure.py +++ /dev/null @@ -1,82 +0,0 @@ -import matplotlib -matplotlib.use('Agg') -import matplotlib.pyplot as plt -import xml.etree.ElementTree as ElementTree -import numpy as np -import h5py -from numpy import genfromtxt - -def main(): - - # File paths - hdf5FilePathPressure = "pressureHistory.hdf5" - - - # Read simulation output from HDF5 file - hf = h5py.File(hdf5FilePathPressure, 'r') - timePressure = hf.get('pressure Time') - timePressure = np.asarray(timePressure) - centerPressure = hf.get('pressure elementCenter') - centerPressure = np.asarray(centerPressure) - pressure = hf.get('pressure') - pressure = np.asarray(pressure) - - time = 1 - posElement1 = -1 - posElement2 = -1 - posElement3 = -1 - last = -1 - for j in range(0, centerPressure.shape[1]): - if centerPressure[0,j,1] >= 0 and posElement1 == -1: - posElement1 = j - if centerPressure[0,j,1] >= 4.2 and posElement2 == -1: - posElement2 = j - if centerPressure[0,j,1] >= 5.6 and posElement3 == -1: - posElement3 = j - - for j in range(0, timePressure.shape[0]): - if j > 0 and timePressure[j] < 1e-12 and last == -1: - last = j - - pressure_0 = genfromtxt('thermoConsolidationPressure_0m.csv', delimiter=',') - pressure_4p2 = genfromtxt('thermoConsolidationPressure_4p2m.csv', delimiter=',') - pressure_5p6 = genfromtxt('thermoConsolidationPressure_5p6m.csv', delimiter=',') - - plt.plot( timePressure[0:last], - pressure[0:last,posElement1], - 'r--', - label='GEOSX: z = 0.0 m') - plt.plot( timePressure[0:last], - pressure[0:last,posElement2], - 'b--', - label='GEOSX: z = 4.2 m') - plt.plot( timePressure[0:last], - pressure[0:last,posElement3], - 'k--', - label='GEOSX: z = 5.6 m') - - plt.plot( pressure_0[:,0], - pressure_0[:,1], - 'r-', - label='Analytical: z = 0.0 m') - plt.plot( pressure_4p2[:,0], - pressure_4p2[:,1], - 'b-', - label='Analytical: z = 4.2 m') - plt.plot( pressure_5p6[:,0], - pressure_5p6[:,1], - 'k-', - label='Analytical: z = 5.6 m') - - plt.xscale("log") - - plt.grid() - plt.ylabel('Pressure [Pa]') - plt.xlabel('Time [s]') - plt.xlim(0.01,100000) - - plt.legend(loc='lower left') - plt.show() - -if __name__ == "__main__": - main() diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationTemperatureFigure.py b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationTemperatureFigure.py deleted file mode 100644 index ee60cf94ab5..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidationTemperatureFigure.py +++ /dev/null @@ -1,85 +0,0 @@ -import matplotlib -matplotlib.use('Agg') -import matplotlib.pyplot as plt -import xml.etree.ElementTree as ElementTree -import numpy as np -import h5py -from numpy import genfromtxt - -def main(): - - # File paths - hdf5FilePathTemperature = "temperatureHistory.hdf5" - - - # Read simulation output from HDF5 file - hf = h5py.File(hdf5FilePathTemperature, 'r') - timeTemperature = hf.get('temperature Time') - timeTemperature = np.asarray(timeTemperature) - centerTemperature = hf.get('temperature elementCenter') - centerTemperature = np.asarray(centerTemperature) - temperature = hf.get('temperature') - temperature = np.asarray(temperature) - - time = 1 - posElement1 = -1 - posElement2 = -1 - posElement3 = -1 - last = -1 - - for j in range(0, centerTemperature.shape[1]): - if centerTemperature[0,j,1] >= 0 and posElement1 == -1: - posElement1 = j - if centerTemperature[0,j,1] >= 4.2 and posElement2 == -1: - posElement2 = j - if centerTemperature[0,j,1] >= 5.6 and posElement3 == -1: - posElement3 = j - - for j in range(0, timeTemperature.shape[0]): - if j > 0 and timeTemperature[j] < 1e-12 and last == -1: - last = j - - temperature_0 = genfromtxt('thermoConsolidationTemp_0m.csv', delimiter=',') - temperature_4p2 = genfromtxt('thermoConsolidationTemp_4p2m.csv', delimiter=',') - temperature_5p6 = genfromtxt('thermoConsolidationTemp_5p6m.csv', delimiter=',') - - plt.plot( timeTemperature[0:last], - temperature[0:last,posElement1], - 'r--', - label='GEOSX: z = 0.0 m') - plt.plot( timeTemperature[0:last], - temperature[0:last,posElement2], - 'b--', - label='GEOSX: z = 4.2 m') - plt.plot( timeTemperature[0:last], - temperature[0:last,posElement3], - 'k--', - label='GEOSX: z = 5.6 m') - - plt.plot( temperature_0[:,0], - 273+temperature_0[:,1], - 'r-', - label='Analytical: z = 0.0 m') - plt.plot( temperature_4p2[:,0], - 273+temperature_4p2[:,1], - 'b-', - label='Analytical: z = 4.2 m') - plt.plot( temperature_5p6[:,0], - 273+temperature_5p6[:,1], - 'k-', - label='Analytical: z = 5.6 m') - - plt.xscale("log") - - plt.grid() - plt.ylabel('Temperature [K]') - plt.xlabel('Time [s]') - plt.ylim(273,323) - plt.xlim(0.01,100000) - - plt.legend(loc='upper left') - plt.show() - #plt.savefig(' temperature.png') - -if __name__ == "__main__": - main() diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidation_analytical.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidation_analytical.csv new file mode 100644 index 00000000000..777f6878ad5 --- /dev/null +++ b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidation_analytical.csv @@ -0,0 +1,101 @@ +time,p_zh0p0,p_zh0p6,p_zh0p8,T_zh0p0,T_zh0p6,T_zh0p8,disp_1p4m,disp_4p2m,disp_7m,sxx_0m,sxx_4p2m,sxx_5p6m,syy +1.00000000e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,8.17869011e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +1.17681195e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,8.31934029e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +1.38488637e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,8.47191345e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +1.62975083e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,8.63742920e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +1.91791026e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,8.81698329e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +2.25701972e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,9.01176446e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +2.65608778e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,9.22306539e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +3.12571585e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,9.45228656e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +3.67837977e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,9.70094789e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +4.32876128e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,9.97069823e-06,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +5.09413801e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.02633257e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +5.99484250e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.05807710e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +7.05480231e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.09251387e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +8.30217568e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.12987121e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +9.77009957e-02,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.17039681e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +1.14975700e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.21435938e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +1.35304777e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.26205039e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +1.59228279e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.31378606e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +1.87381742e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.36990941e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +2.20513074e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.43079257e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +2.59502421e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.49683920e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +3.05385551e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.56848723e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +3.59381366e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.64621169e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +4.22924287e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.73052795e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +4.97702356e-01,9.88023952e-01,9.88023952e-01,9.88023952e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.82199503e-05,-9.96007984e-01,-9.96007984e-01,-9.96007984e-01,-1.00000000e+00 +5.85702082e-01,9.88023952e-01,9.88023952e-01,9.88023943e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,1.92121940e-05,-9.96007984e-01,-9.96007984e-01,-9.96007981e-01,-1.00000000e+00 +6.89261210e-01,9.88023952e-01,9.88023952e-01,9.88023833e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,2.02885896e-05,-9.96007984e-01,-9.96007984e-01,-9.96007944e-01,-1.00000000e+00 +8.11130831e-01,9.88023952e-01,9.88023952e-01,9.88022897e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,2.14562739e-05,-9.96007984e-01,-9.96007984e-01,-9.96007632e-01,-1.00000000e+00 +9.54548457e-01,9.88023952e-01,9.88023952e-01,9.88017165e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,2.27229890e-05,-9.96007984e-01,-9.96007984e-01,-9.96005722e-01,-1.00000000e+00 +1.12332403e+00,9.88023952e-01,9.88023952e-01,9.87990595e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,2.40971338e-05,-9.96007984e-01,-9.96007984e-01,-9.95996865e-01,-1.00000000e+00 +1.32194115e+00,9.88023952e-01,9.88023952e-01,9.87893570e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,2.55878194e-05,-9.96007984e-01,-9.96007984e-01,-9.95964523e-01,-1.00000000e+00 +1.55567614e+00,9.88023952e-01,9.88023952e-01,9.87604571e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,2.72049294e-05,-9.96007984e-01,-9.96007984e-01,-9.95868190e-01,-1.00000000e+00 +1.83073828e+00,9.88023952e-01,9.88023952e-01,9.86881173e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217565e-06,2.89591860e-05,-9.96007984e-01,-9.96007984e-01,-9.95627058e-01,-1.00000000e+00 +2.15443469e+00,9.88023952e-01,9.88023950e-01,9.85320154e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217566e-06,3.08622204e-05,-9.96007984e-01,-9.96007983e-01,-9.95106718e-01,-1.00000000e+00 +2.53536449e+00,9.88023952e-01,9.88023918e-01,9.82352065e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217588e-06,3.29266505e-05,-9.96007984e-01,-9.96007973e-01,-9.94117355e-01,-1.00000000e+00 +2.98364724e+00,9.88023952e-01,9.88023593e-01,9.77285710e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91217847e-06,3.51661642e-05,-9.96007984e-01,-9.96007864e-01,-9.92428570e-01,-1.00000000e+00 +3.51119173e+00,9.88023952e-01,9.88021253e-01,9.69397165e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91220039e-06,3.75956102e-05,-9.96007984e-01,-9.96007084e-01,-9.89799055e-01,-1.00000000e+00 +4.13201240e+00,9.88023952e-01,9.88008794e-01,9.58039210e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91233722e-06,4.02310966e-05,-9.96007984e-01,-9.96002931e-01,-9.86013070e-01,-1.00000000e+00 +4.86260158e+00,9.88023952e-01,9.87957599e-01,9.42740431e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91299699e-06,4.30900977e-05,-9.96007984e-01,-9.95985866e-01,-9.80913477e-01,-1.00000000e+00 +5.72236766e+00,9.88023952e-01,9.87788923e-01,9.23270043e-01,2.73150000e+02,2.73150000e+02,2.73150000e+02,1.30405855e-06,3.91554847e-06,4.61915697e-05,-9.96007984e-01,-9.95929641e-01,-9.74423348e-01,-1.00000000e+00 +6.73415066e+00,9.88023952e-01,9.87328743e-01,8.99658433e-01,2.73150000e+02,2.73150000e+02,2.73150003e+02,1.30405855e-06,3.92372163e-06,4.95560764e-05,-9.96007984e-01,-9.95776248e-01,-9.66552821e-01,-1.00000000e+00 +7.92482898e+00,9.88023952e-01,9.86260031e-01,8.72177065e-01,2.73150000e+02,2.73150000e+02,2.73150033e+02,1.30405855e-06,3.94601524e-06,5.32059256e-05,-9.96007984e-01,-9.95420010e-01,-9.57392454e-01,-1.00000000e+00 +9.32603347e+00,9.88023952e-01,9.84096656e-01,8.41289927e-01,2.73150000e+02,2.73150000e+02,2.73150228e+02,1.30405866e-06,3.99903367e-06,5.71653173e-05,-9.96007984e-01,-9.94698885e-01,-9.47097325e-01,-1.00000000e+00 +1.09749877e+01,9.88023952e-01,9.80201408e-01,8.07591207e-01,2.73150000e+02,2.73150000e+02,2.73151190e+02,1.30406017e-06,4.11121207e-06,6.14605035e-05,-9.96007984e-01,-9.93400469e-01,-9.35867304e-01,-1.00000000e+00 +1.29154967e+01,9.88023950e-01,9.73854119e-01,7.71741948e-01,2.73150000e+02,2.73150000e+02,2.73154898e+02,1.30407539e-06,4.32606037e-06,6.61199630e-05,-9.96007983e-01,-9.91284706e-01,-9.23928677e-01,-1.00000000e+00 +1.51991108e+01,9.88023917e-01,9.64354354e-01,7.34414392e-01,2.73150000e+02,2.73150000e+02,2.73166469e+02,1.30418577e-06,4.70407124e-06,7.11745895e-05,-9.96007972e-01,-9.88118118e-01,-9.11520870e-01,-1.00000000e+00 +1.78864953e+01,9.88023546e-01,9.51129513e-01,6.96248569e-01,2.73150000e+02,2.73150000e+02,2.73196600e+02,1.30479056e-06,5.32279538e-06,7.66578970e-05,-9.96007849e-01,-9.83709838e-01,-8.98889324e-01,-1.00000000e+00 +2.10490414e+01,9.88020634e-01,9.33819510e-01,6.57822407e-01,2.73150000e+02,2.73150000e+02,2.73263855e+02,1.30739925e-06,6.27510205e-06,8.26062419e-05,-9.96006878e-01,-9.77939837e-01,-8.86282366e-01,-1.00000000e+00 +2.47707636e+01,9.88003962e-01,9.12320429e-01,6.19634534e-01,2.73150000e+02,2.73150001e+02,2.73395465e+02,1.31657057e-06,7.66608608e-06,8.90590638e-05,-9.96001321e-01,-9.70773479e-01,-8.73947905e-01,-1.00000000e+00 +2.91505306e+01,9.87931038e-01,8.86784226e-01,5.82097685e-01,2.73150000e+02,2.73150011e+02,2.73625695e+02,1.34364188e-06,9.60931207e-06,9.60591471e-05,-9.95977013e-01,-9.62261441e-01,-8.62126313e-01,-1.00000000e+00 +3.43046929e+01,9.87677619e-01,8.57583019e-01,5.45540117e-01,2.73150000e+02,2.73150087e+02,2.73991745e+02,1.41243217e-06,1.22230839e-05,1.03652905e-04,-9.95892540e-01,-9.52527935e-01,-8.51038608e-01,-1.00000000e+00 +4.03701726e+01,9.86954060e-01,8.25252042e-01,5.10212287e-01,2.73150000e+02,2.73150524e+02,2.74528250e+02,1.56609865e-06,1.56272577e-05,1.11890687e-04,-9.95651353e-01,-9.41752254e-01,-8.40872179e-01,-1.00000000e+00 +4.75081016e+01,9.85207153e-01,7.90425256e-01,4.76296369e-01,2.73150000e+02,2.73152429e+02,2.75261937e+02,1.87314178e-06,1.99408581e-05,1.20827112e-04,-9.95069051e-01,-9.30149040e-01,-8.31767934e-01,-1.00000000e+00 +5.59081018e+01,9.81551858e-01,7.53773758e-01,4.43916659e-01,2.73150000e+02,2.73159029e+02,2.76207802e+02,2.42982458e-06,2.52804641e-05,1.30521432e-04,-9.93850619e-01,-9.17951674e-01,-8.23812292e-01,-1.00000000e+00 +6.57933225e+01,9.74781911e-01,7.15950976e-01,4.13149201e-01,2.73150000e+02,2.73177825e+02,2.77367502e+02,3.35659696e-06,3.17590073e-05,1.41037920e-04,-9.91593970e-01,-9.05400468e-01,-8.17035572e-01,-1.00000000e+00 +7.74263683e+01,9.63481754e-01,6.77542053e-01,3.84028354e-01,2.73150000e+02,2.73223106e+02,2.78729926e+02,4.78767417e-06,3.94842353e-05,1.52446268e-04,-9.87827251e-01,-8.92733336e-01,-8.11415896e-01,-1.00000000e+00 +9.11162756e+01,9.46219147e-01,6.39011242e-01,3.56546048e-01,2.73150000e+02,2.73317677e+02,2.80273401e+02,6.85520752e-06,4.85556856e-05,1.64821859e-04,-9.82073049e-01,-8.80173445e-01,-8.06885551e-01,-1.00000000e+00 +1.07226722e+02,9.21764598e-01,6.00645555e-01,3.30638136e-01,2.73150000e+02,2.73492553e+02,2.81968774e+02,9.67147469e-06,5.90588280e-05,1.78245386e-04,-9.73921533e-01,-8.67909509e-01,-8.03335701e-01,-1.00000000e+00 +1.26185688e+02,8.89277215e-01,5.62505636e-01,3.06156723e-01,2.73150000e+02,2.73784054e+02,2.83782693e+02,1.33132510e-05,7.10554910e-05,1.92800372e-04,-9.63092405e-01,-8.56070708e-01,-8.00616988e-01,-1.00000000e+00 +1.48496826e+02,8.48414806e-01,5.24407806e-01,2.82840565e-01,2.73150001e+02,2.74228851e+02,2.85680558e+02,1.78117179e-05,8.45708777e-05,2.08565861e-04,-9.49471605e-01,-8.44705822e-01,-7.98538529e-01,-1.00000000e+00 +1.74752840e+02,7.99355563e-01,4.85962917e-01,2.60309700e-01,2.73150012e+02,2.74858341e+02,2.87628881e+02,2.31493625e-05,9.95787686e-05,2.25600887e-04,-9.33118557e-01,-8.33779330e-01,-7.96873210e-01,-1.00000000e+00 +2.05651231e+02,7.42745888e-01,4.46682751e-01,2.38109259e-01,2.73150105e+02,2.75693908e+02,2.89596955e+02,2.92630638e-05,1.15986907e-04,2.43918737e-04,-9.14248945e-01,-8.23192641e-01,-7.95377286e-01,-1.00000000e+00 +2.42012826e+02,6.79607595e-01,4.06137772e-01,2.15804099e-01,2.73150681e+02,2.76744134e+02,2.91557869e+02,3.60505713e-05,1.33623806e-04,2.63452805e-04,-8.93204574e-01,-8.12828325e-01,-7.93824973e-01,-1.00000000e+00 +2.84803587e+02,6.11245778e-01,3.64124312e-01,1.93094537e-01,2.73153356e+02,2.78004292e+02,2.93488975e+02,4.33764565e-05,1.52227267e-04,2.84020455e-04,-8.70425327e-01,-8.02604313e-01,-7.92048436e-01,-1.00000000e+00 +3.35160265e+02,5.39192187e-01,3.20790211e-01,1.69908215e-01,2.73163152e+02,2.79457800e+02,2.95371950e+02,5.10741701e-05,1.71435084e-04,3.05294010e-04,-8.46436852e-01,-7.92520136e-01,-7.89968589e-01,-1.00000000e+00 +3.94420606e+02,4.65194096e-01,2.76683290e-01,1.46435389e-01,2.73192401e+02,2.81078945e+02,2.97192577e+02,5.89435792e-05,1.90780239e-04,3.26787406e-04,-8.21858569e-01,-7.82681264e-01,-7.87606193e-01,-1.00000000e+00 +4.64158883e+02,3.91222899e-01,2.32717737e-01,1.23103723e-01,2.73265766e+02,2.82836146e+02,2.98940346e+02,6.67473538e-05,2.09695342e-04,3.47864170e-04,-7.97421597e-01,-7.73297682e-01,-7.85072280e-01,-1.00000000e+00 +5.46227722e+02,3.19451682e-01,1.90077892e-01,1.00511579e-01,2.73424356e+02,2.84695163e+02,3.00607988e+02,7.42124025e-05,2.27532504e-04,3.67769524e-04,-7.73973628e-01,-7.64661453e-01,-7.82544491e-01,-1.00000000e+00 +6.42807312e+02,2.52155069e-01,1.50082653e-01,7.93419413e-02,2.73726345e+02,2.86621966e+02,3.02190982e+02,8.10421423e-05,2.43603846e-04,3.85687698e-04,-7.52447392e-01,-7.57110115e-01,-7.80236927e-01,-1.00000000e+00 +7.56463328e+02,1.91514830e-01,1.14024001e-01,6.02687432e-02,2.74242478e+02,2.88585318e+02,3.03687171e+02,8.69417950e-05,2.57243851e-04,4.00823272e-04,-7.33782378e-01,-7.50980621e-01,-7.78367761e-01,-1.00000000e+00 +8.90215085e+02,1.39353643e-01,8.29891068e-02,4.38598184e-02,2.75046990e+02,2.90559400e+02,3.05096656e+02,9.16539453e-05,2.67888292e-04,4.12499894e-04,-7.18808852e-01,-7.46557903e-01,-7.77126575e-01,-1.00000000e+00 +1.04761575e+03,9.68530025e-02,5.76864100e-02,3.04845616e-02,2.76206397e+02,2.92526669e+02,3.06422278e+02,9.49970481e-05,2.75156571e-04,4.20260775e-04,-7.08120193e-01,-7.44025476e-01,-7.76645020e-01,-1.00000000e+00 +1.23284674e+03,6.43340263e-02,3.83134196e-02,2.02442454e-02,2.77769304e+02,2.94480585e+02,3.07670842e+02,9.68965968e-05,2.78917247e-04,4.23946656e-04,-7.01969253e-01,-7.43429562e-01,-7.76977274e-01,-1.00000000e+00 +1.45082878e+03,4.11835933e-02,2.45117858e-02,1.29481760e-02,2.79759702e+02,2.96427114e+02,3.08854677e+02,9.73984330e-05,2.79314203e-04,4.23724605e-04,-7.00223636e-01,-7.44668605e-01,-7.78096758e-01,-1.00000000e+00 +1.70735265e+03,2.59801613e-02,1.54411989e-02,8.15222216e-03,2.82174785e+02,2.98383600e+02,3.09992319e+02,9.66574585e-05,2.76737853e-04,4.20050508e-04,-7.02401075e-01,-7.47514533e-01,-7.79911031e-01,-1.00000000e+00 +2.00923300e+03,1.68093586e-02,9.96630735e-03,5.25688195e-03,2.84986836e+02,3.00374198e+02,3.11106916e+02,9.49028896e-05,2.71743333e-04,4.13569900e-04,-7.07780293e-01,-7.51661362e-01,-7.82289708e-01,-1.00000000e+00 +2.36448941e+03,1.16793252e-02,6.90339295e-03,3.63705245e-03,2.88147662e+02,3.02422311e+02,3.12221913e+02,9.23907029e-05,2.64940330e-04,4.04988156e-04,-7.15552761e-01,-7.56784729e-01,-7.85094756e-01,-1.00000000e+00 +2.78255940e+03,8.89399991e-03,5.24267504e-03,2.75923218e-03,2.91593512e+02,3.04541910e+02,3.13355169e+02,8.93600693e-05,2.56894110e-04,3.94956971e-04,-7.24961870e-01,-7.62589954e-01,-7.88201917e-01,-1.00000000e+00 +3.27454916e+03,7.26074134e-03,4.27264770e-03,2.24724600e-03,2.95248505e+02,3.06730349e+02,3.14513750e+02,8.60086305e-05,2.48072288e-04,3.84016204e-04,-7.35382429e-01,-7.68831928e-01,-7.91506998e-01,-1.00000000e+00 +3.85352859e+03,6.09595145e-03,3.58442748e-03,1.88470954e-03,2.99025667e+02,3.08964804e+02,3.15691343e+02,8.24912495e-05,2.38846676e-04,3.72599742e-04,-7.46325650e-01,-7.75305888e-01,-7.94918933e-01,-1.00000000e+00 +4.53487851e+03,5.09223168e-03,2.99341994e-03,1.57379047e-03,3.02826463e+02,3.11203000e+02,3.16868864e+02,7.89346209e-05,2.29530377e-04,3.61080402e-04,-7.57393468e-01,-7.81823473e-01,-7.98347857e-01,-1.00000000e+00 +5.33669923e+03,4.15364553e-03,2.44149730e-03,1.28358027e-03,3.06541252e+02,3.13387440e+02,3.18017482e+02,7.54541729e-05,2.20417326e-04,3.49815196e-04,-7.68224971e-01,-7.88192819e-01,-8.01696972e-01,-1.00000000e+00 +6.28029144e+03,3.27431143e-03,1.92459736e-03,1.01182179e-03,3.10053276e+02,3.15451917e+02,3.19102874e+02,7.21628955e-05,2.11800459e-04,3.39164015e-04,-7.78467932e-01,-7.94213951e-01,-8.04862563e-01,-1.00000000e+00 +7.39072203e+03,2.47562389e-03,1.45513565e-03,7.65010168e-04,3.13247535e+02,3.17329475e+02,3.20089969e+02,6.91693039e-05,2.03963111e-04,3.29476503e-04,-7.87784479e-01,-7.99690138e-01,-8.07741576e-01,-1.00000000e+00 +8.69749003e+03,1.78153846e-03,1.04716206e-03,5.50525678e-04,3.16023847e+02,3.18961352e+02,3.20947898e+02,6.65673936e-05,1.97151219e-04,3.21056540e-04,-7.95882054e-01,-8.04449778e-01,-8.10243868e-01,-1.00000000e+00 +1.02353102e+04,1.20960217e-03,7.10986318e-04,3.73787628e-04,3.18311591e+02,3.20306055e+02,3.21654849e+02,6.44233602e-05,1.91538067e-04,3.14118302e-04,-8.02554641e-01,-8.08371826e-01,-8.12305811e-01,-1.00000000e+00 +1.20450354e+04,7.66935046e-04,4.50793110e-04,2.36995963e-04,3.20082260e+02,3.21346828e+02,3.22202016e+02,6.27639209e-05,1.87193598e-04,3.08748243e-04,-8.07719091e-01,-8.11407414e-01,-8.13901714e-01,-1.00000000e+00 +1.41747416e+04,4.48627984e-04,2.63696913e-04,1.38633671e-04,3.21355488e+02,3.22095212e+02,3.22595465e+02,6.15706738e-05,1.84069637e-04,3.04886815e-04,-8.11432674e-01,-8.13590203e-01,-8.15049274e-01,-1.00000000e+00 +1.66810054e+04,2.38692924e-04,1.40300180e-04,7.37601698e-05,3.22195228e+02,3.22588799e+02,3.22854959e+02,6.07836841e-05,1.82009271e-04,3.02340063e-04,-8.13881916e-01,-8.15029831e-01,-8.15806131e-01,-1.00000000e+00 +1.96304065e+04,1.13589444e-04,6.67662002e-05,3.51010687e-05,3.22695642e+02,3.22882935e+02,3.23009596e+02,6.03147050e-05,1.80781468e-04,3.00822414e-04,-8.15341456e-01,-8.15887728e-01,-8.16257154e-01,-1.00000000e+00 +2.31012970e+04,4.74039219e-05,2.78633262e-05,1.46486175e-05,3.22960384e+02,3.23038547e+02,3.23091406e+02,6.00665934e-05,1.80131903e-04,3.00019508e-04,-8.16113621e-01,-8.16341595e-01,-8.16495766e-01,-1.00000000e+00 +2.71858824e+04,1.69506457e-05,9.96333956e-06,5.23803759e-06,3.23082197e+02,3.23110147e+02,3.23129048e+02,5.99524323e-05,1.79833025e-04,2.99650075e-04,-8.16468909e-01,-8.16550428e-01,-8.16605556e-01,-1.00000000e+00 +3.19926714e+04,5.05346248e-06,2.97035072e-06,1.56160579e-06,3.23129786e+02,3.23138119e+02,3.23143754e+02,5.99078330e-05,1.79716263e-04,2.99505749e-04,-8.16607710e-01,-8.16632013e-01,-8.16648448e-01,-1.00000000e+00 +3.76493581e+04,1.21635397e-06,7.14954925e-07,3.75874048e-07,3.23145135e+02,3.23147140e+02,3.23148497e+02,5.98934487e-05,1.79678604e-04,2.99459200e-04,-8.16652476e-01,-8.16658326e-01,-8.16662281e-01,-1.00000000e+00 +4.43062146e+04,2.27597547e-07,1.33778481e-07,7.03315098e-08,3.23149090e+02,3.23149465e+02,3.23149719e+02,5.98897421e-05,1.79668900e-04,2.99447205e-04,-8.16664011e-01,-8.16665106e-01,-8.16665846e-01,-1.00000000e+00 +5.21400829e+04,3.16646797e-08,1.86120317e-08,9.78492415e-09,3.23149873e+02,3.23149926e+02,3.23149961e+02,5.98890076e-05,1.79666977e-04,2.99444829e-04,-8.16666297e-01,-8.16666450e-01,-8.16666553e-01,-1.00000000e+00 +6.13590727e+04,3.10832795e-09,1.82702933e-09,9.60526161e-10,3.23149988e+02,3.23149993e+02,3.23149996e+02,5.98889005e-05,1.79666697e-04,2.99444482e-04,-8.16666630e-01,-8.16666645e-01,-8.16666655e-01,-1.00000000e+00 +7.22080902e+04,2.02415416e-10,1.18976797e-10,6.25498036e-11,3.23149999e+02,3.23150000e+02,3.23150000e+02,5.98888896e-05,1.79666669e-04,2.99444447e-04,-8.16666664e-01,-8.16666665e-01,-8.16666666e-01,-1.00000000e+00 +8.49753436e+04,8.13226050e-12,4.78002279e-12,2.51300670e-12,3.23150000e+02,3.23150000e+02,3.23150000e+02,5.98888889e-05,1.79666667e-04,2.99444445e-04,-8.16666667e-01,-8.16666667e-01,-8.16666667e-01,-1.00000000e+00 +1.00000000e+05,1.85073305e-13,1.08783359e-13,5.71907963e-14,3.23150000e+02,3.23150000e+02,3.23150000e+02,5.98888889e-05,1.79666667e-04,2.99444444e-04,-8.16666667e-01,-8.16666667e-01,-8.16666667e-01,-1.00000000e+00 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidation_geos.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidation_geos.csv new file mode 100644 index 00000000000..1457252d061 --- /dev/null +++ b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermalConsolidation_geos.csv @@ -0,0 +1,188 @@ +time,p_0m,p_4p2m,p_5p6m,T_0m,T_4p2m,T_5p6m,disp_1p4m,disp_4p2m,disp_7m,sxx_0m,sxx_4p2m,sxx_5p6m,syy_0m +0.00000000e+00,0.00000000e+00,0.00000000e+00,0.00000000e+00,2.73000000e+02,2.73000000e+02,2.73000000e+02,-0.00000000e+00,-0.00000000e+00,-0.00000000e+00,0.00000000e+00,0.00000000e+00,0.00000000e+00,0.00000000e+00 +1.00000000e-02,1.00000000e+00,1.00000000e+00,1.00000000e+00,2.73000000e+02,2.73000000e+02,2.73000000e+02,-3.97059798e-21,4.37445479e-20,6.93147483e-07,-1.00000000e+00,-1.00000000e+00,-1.00000000e+00,-1.00000000e+00 +1.00000000e+00,1.00000000e+00,9.99994468e-01,9.97394575e-01,2.73000000e+02,2.73000000e+02,2.72999939e+02,3.61756052e-16,7.79466173e-11,1.52594600e-05,-1.00000000e+00,-9.99998812e-01,-9.99440315e-01,-1.00000000e+00 +3.00000000e+00,9.99999999e-01,9.99621614e-01,9.72884388e-01,2.73000000e+02,2.72999991e+02,2.72999564e+02,1.37143557e-12,8.12223584e-09,2.78152035e-05,-1.00000000e+00,-9.99907358e-01,-9.93297274e-01,-1.00000000e+00 +5.00000000e+00,9.99999984e-01,9.98358609e-01,9.35007484e-01,2.73000000e+02,2.72999961e+02,2.72999719e+02,1.25549597e-11,3.97543495e-08,3.66496698e-05,-9.99999995e-01,-9.99584532e-01,-9.83073657e-01,-1.00000000e+00 +7.00000000e+00,9.99999911e-01,9.95725675e-01,8.91763758e-01,2.73000000e+02,2.72999899e+02,2.73001857e+02,6.21646892e-11,1.15085171e-07,4.38232944e-05,-9.99999970e-01,-9.98888742e-01,-9.70754333e-01,-1.00000000e+00 +9.00000000e+00,9.99999638e-01,9.91448024e-01,8.48447768e-01,2.73000000e+02,2.72999799e+02,2.73007979e+02,2.20165638e-10,2.53622396e-07,5.00089386e-05,-9.99999879e-01,-9.97725855e-01,-9.57904833e-01,-1.00000000e+00 +1.10000000e+01,9.99998837e-01,9.85481141e-01,8.07653043e-01,2.73000000e+02,2.72999659e+02,2.73020451e+02,6.24367476e-10,4.70738895e-07,5.55228878e-05,-9.99999612e-01,-9.96063799e-01,-9.45431013e-01,-1.00000000e+00 +1.30000000e+01,9.99996854e-01,9.77956133e-01,7.70325090e-01,2.73000000e+02,2.72999483e+02,2.73041761e+02,1.50692270e-09,7.76284905e-07,6.05432530e-05,-9.99998951e-01,-9.93923276e-01,-9.33757962e-01,-1.00000000e+00 +1.50000000e+01,9.99992549e-01,9.69104287e-01,7.36588371e-01,2.73000000e+02,2.72999277e+02,2.73074256e+02,3.21506383e-09,1.17479197e-06,6.51816908e-05,-9.99997516e-01,-9.91359242e-01,-9.23034122e-01,-1.00000000e+00 +1.70000000e+01,9.99984127e-01,9.59193781e-01,7.06215204e-01,2.73000000e+02,2.72999048e+02,2.73119937e+02,6.22180049e-09,1.66647392e-06,6.95134729e-05,-9.99994707e-01,-9.88443123e-01,-9.13268146e-01,-1.00000000e+00 +1.90000000e+01,9.99968997e-01,9.48488060e-01,6.78857108e-01,2.72999999e+02,2.72998806e+02,2.73180329e+02,1.11247773e-08,2.24843486e-06,7.35921084e-05,-9.99989663e-01,-9.85249804e-01,-9.04405710e-01,-1.00000000e+00 +2.10000000e+01,9.99943678e-01,9.37223504e-01,6.54149027e-01,2.72999999e+02,2.72998558e+02,2.73256427e+02,1.86339358e-08,2.91576718e-06,7.74572119e-05,-9.99981221e-01,-9.81849722e-01,-8.96368869e-01,-1.00000000e+00 +2.30000000e+01,9.99903761e-01,9.25600500e-01,6.31751677e-01,2.72999998e+02,2.72998317e+02,2.73348711e+02,2.95503626e-08,3.66242090e-06,8.11390663e-05,-9.99967912e-01,-9.78304941e-01,-8.89074961e-01,-1.00000000e+00 +2.50000000e+01,9.99843943e-01,9.13782467e-01,6.11365462e-01,2.72999996e+02,2.72998092e+02,2.73457197e+02,4.47395416e-08,4.48183402e-06,8.46614289e-05,-9.99947967e-01,-9.74667828e-01,-8.82445025e-01,-1.00000000e+00 +2.70000000e+01,9.99758122e-01,9.01898926e-01,5.92732032e-01,2.72999994e+02,2.72997898e+02,2.73581516e+02,6.51022662e-08,5.36735994e-06,8.80433405e-05,-9.99919353e-01,-9.70981236e-01,-8.76407079e-01,-1.00000000e+00 +2.90000000e+01,9.99639541e-01,8.90050202e-01,5.75630949e-01,2.72999992e+02,2.72997750e+02,2.73720995e+02,9.15459977e-08,6.31253929e-06,9.13003380e-05,-9.99879815e-01,-9.67279423e-01,-8.70896950e-01,-1.00000000e+00 +3.10000000e+01,9.99480970e-01,8.78312430e-01,5.59874825e-01,2.72999988e+02,2.72997670e+02,2.73874746e+02,1.24958725e-07,7.31126035e-06,9.44452939e-05,-9.99826944e-01,-9.63589255e-01,-8.65858040e-01,-1.00000000e+00 +3.30000000e+01,9.99274905e-01,8.66742191e-01,5.45304347e-01,2.72999983e+02,2.72997678e+02,2.74041728e+02,1.66186625e-07,8.35784304e-06,9.74890144e-05,-9.99758238e-01,-9.59931442e-01,-8.61240643e-01,-1.00000000e+00 +3.50000000e+01,9.99013767e-01,8.55380541e-01,5.31783773e-01,2.72999977e+02,2.72997802e+02,2.74220818e+02,2.16016155e-07,9.44707283e-06,1.00440675e-04,-9.99671169e-01,-9.56321663e-01,-8.57001154e-01,-1.00000000e+00 +3.70000000e+01,9.98690095e-01,8.44256357e-01,5.19197048e-01,2.72999969e+02,2.72998073e+02,2.74410853e+02,2.75160685e-07,1.05742028e-05,1.03308146e-04,-9.99563250e-01,-9.52771559e-01,-8.53101286e-01,-1.00000000e+00 +3.90000000e+01,9.98296713e-01,8.33389055e-01,5.07444545e-01,2.72999960e+02,2.72998524e+02,2.74610673e+02,3.44251423e-07,1.17349369e-05,1.06098237e-04,-9.99432088e-01,-9.49289552e-01,-8.49507359e-01,-1.00000000e+00 +4.10000000e+01,9.97826878e-01,8.22790766e-01,4.96440363e-01,2.72999949e+02,2.72999194e+02,2.74819143e+02,4.23832162e-07,1.29254020e-05,1.08816887e-04,-9.99275435e-01,-9.45881533e-01,-8.46189677e-01,-1.00000000e+00 +4.30000000e+01,9.97274389e-01,8.12468044e-01,4.86110107e-01,2.72999936e+02,2.73000124e+02,2.75035179e+02,5.14357277e-07,1.41421156e-05,1.11469314e-04,-9.99091224e-01,-9.42551404e-01,-8.43121987e-01,-1.00000000e+00 +4.50000000e+01,9.96633682e-01,8.02423215e-01,4.76389057e-01,2.72999921e+02,2.73001356e+02,2.75257756e+02,6.16192362e-07,1.53819511e-05,1.14060128e-04,-9.98877598e-01,-9.39301523e-01,-8.40281020e-01,-1.00000000e+00 +4.70000000e+01,9.95899877e-01,7.92655428e-01,4.67220660e-01,2.72999904e+02,2.73002939e+02,2.75485916e+02,7.29616920e-07,1.66421037e-05,1.16593431e-04,-9.98632932e-01,-9.36133054e-01,-8.37646100e-01,-1.00000000e+00 +4.90000000e+01,9.95068817e-01,7.83161465e-01,4.58555292e-01,2.72999884e+02,2.73004920e+02,2.75718774e+02,8.54828592e-07,1.79200584e-05,1.19072888e-04,-9.98355839e-01,-9.33046244e-01,-8.35198804e-01,-1.00000000e+00 +5.10000000e+01,9.94137079e-01,7.73936375e-01,4.50349229e-01,2.72999862e+02,2.73007348e+02,2.75955519e+02,9.91948457e-07,1.92135590e-05,1.21501791e-04,-9.98045178e-01,-9.30040640e-01,-8.32922680e-01,-1.00000000e+00 +5.30000000e+01,9.93101959e-01,7.64973966e-01,4.42563794e-01,2.72999838e+02,2.73010274e+02,2.76195410e+02,1.14102702e-06,2.05205813e-05,1.23883110e-04,-9.97700047e-01,-9.27115265e-01,-8.30803000e-01,-1.00000000e+00 +5.50000000e+01,9.91961456e-01,7.56267174e-01,4.35164646e-01,2.72999811e+02,2.73013748e+02,2.76437774e+02,1.30205058e-06,2.18393076e-05,1.26219533e-04,-9.97319779e-01,-9.24268749e-01,-8.28826552e-01,-1.00000000e+00 +5.70000000e+01,9.90714233e-01,7.47808359e-01,4.28121183e-01,2.72999782e+02,2.73017820e+02,2.76682005e+02,1.47494775e-06,2.31681045e-05,1.28513508e-04,-9.96903929e-01,-9.21499439e-01,-8.26981454e-01,-1.00000000e+00 +5.90000000e+01,9.89359579e-01,7.39589522e-01,4.21406043e-01,2.72999750e+02,2.73022538e+02,2.76927558e+02,1.65959590e-06,2.45055023e-05,1.30767265e-04,-9.96452258e-01,-9.18805476e-01,-8.25257002e-01,-1.00000000e+00 +6.10000000e+01,9.87897359e-01,7.31602475e-01,4.14994679e-01,2.72999715e+02,2.73027951e+02,2.77173945e+02,1.85582739e-06,2.58501779e-05,1.32982847e-04,-9.95964723e-01,-9.16184865e-01,-8.23643532e-01,-1.00000000e+00 +6.30000000e+01,9.86327962e-01,7.23838969e-01,4.08865000e-01,2.72999678e+02,2.73034104e+02,2.77420728e+02,2.06343558e-06,2.72009380e-05,1.35162129e-04,-9.95441453e-01,-9.13635520e-01,-8.22132305e-01,-1.00000000e+00 +6.50000000e+01,9.84652254e-01,7.16290789e-01,4.02997071e-01,2.72999639e+02,2.73041041e+02,2.77667521e+02,2.28218043e-06,2.85567056e-05,1.37306836e-04,-9.94882736e-01,-9.11155305e-01,-8.20715397e-01,-1.00000000e+00 +6.70000000e+01,9.82871517e-01,7.08949829e-01,3.97372842e-01,2.72999597e+02,2.73048802e+02,2.77913979e+02,2.51179368e-06,2.99165075e-05,1.39418563e-04,-9.94289001e-01,-9.08742064e-01,-8.19385620e-01,-1.00000000e+00 +6.90000000e+01,9.80987407e-01,7.01808138e-01,3.91975935e-01,2.72999553e+02,2.73057428e+02,2.78159799e+02,2.75198367e-06,3.12794629e-05,1.41498783e-04,-9.93660799e-01,-9.06393642e-01,-8.18136432e-01,-1.00000000e+00 +7.10000000e+01,9.79001903e-01,6.94857966e-01,3.86791444e-01,2.72999506e+02,2.73066952e+02,2.78404710e+02,3.00243966e-06,3.26447737e-05,1.43548864e-04,-9.92998790e-01,-9.04107904e-01,-8.16961878e-01,-1.00000000e+00 +7.30000000e+01,9.76917255e-01,6.88091781e-01,3.81805772e-01,2.72999457e+02,2.73077409e+02,2.78648479e+02,3.26283586e-06,3.40117161e-05,1.45570076e-04,-9.92303724e-01,-9.01882746e-01,-8.15856522e-01,-1.00000000e+00 +7.50000000e+01,9.74735951e-01,6.81502293e-01,3.77006482e-01,2.72999406e+02,2.73088827e+02,2.78890899e+02,3.53283487e-06,3.53796325e-05,1.47563603e-04,-9.91576431e-01,-8.99716105e-01,-8.14815396e-01,-1.00000000e+00 +7.70000000e+01,9.72460667e-01,6.75082463e-01,3.72382180e-01,2.72999352e+02,2.73101234e+02,2.79131790e+02,3.81209091e-06,3.67479251e-05,1.49530549e-04,-9.90817803e-01,-8.97605966e-01,-8.13833954e-01,-1.00000000e+00 +7.90000000e+01,9.70094240e-01,6.68825507e-01,3.67922395e-01,2.72999297e+02,2.73114652e+02,2.79370995e+02,4.10025259e-06,3.81160494e-05,1.51471946e-04,-9.90028786e-01,-8.95550369e-01,-8.12908028e-01,-1.00000000e+00 +8.10000000e+01,9.67639627e-01,6.62724899e-01,3.63617493e-01,2.72999239e+02,2.73129102e+02,2.79608380e+02,4.39696535e-06,3.94835092e-05,1.53388760e-04,-9.89210366e-01,-8.93547408e-01,-8.12033790e-01,-1.00000000e+00 +8.30000000e+01,9.65099881e-01,6.56774368e-01,3.59458586e-01,2.72999179e+02,2.73144602e+02,2.79843830e+02,4.70187364e-06,4.08498518e-05,1.55281896e-04,-9.88363561e-01,-8.91595235e-01,-8.11207725e-01,-1.00000000e+00 +8.50000000e+01,9.62478121e-01,6.50967896e-01,3.55437460e-01,2.72999118e+02,2.73161165e+02,2.80077245e+02,5.01462275e-06,4.22146636e-05,1.57152207e-04,-9.87489411e-01,-8.89692064e-01,-8.10426592e-01,-1.00000000e+00 +8.70000000e+01,9.59777509e-01,6.45299715e-01,3.51546513e-01,2.72999054e+02,2.73178803e+02,2.80308542e+02,5.33486047e-06,4.35775663e-05,1.59000493e-04,-9.86588970e-01,-8.87836165e-01,-8.09687406e-01,-1.00000000e+00 +8.90000000e+01,9.57001231e-01,6.39764294e-01,3.47778693e-01,2.72998989e+02,2.73197525e+02,2.80537651e+02,5.66223838e-06,4.49382138e-05,1.60827505e-04,-9.85663300e-01,-8.86025869e-01,-8.08987412e-01,-1.00000000e+00 +9.10000000e+01,9.54152477e-01,6.34356337e-01,3.44127448e-01,2.72998922e+02,2.73217337e+02,2.80764514e+02,5.99641310e-06,4.62962889e-05,1.62633956e-04,-9.84713465e-01,-8.84259566e-01,-8.08324062e-01,-1.00000000e+00 +9.30000000e+01,9.51234425e-01,6.29070776e-01,3.40586682e-01,2.72998853e+02,2.73238241e+02,2.80989084e+02,6.33704722e-06,4.76515008e-05,1.64420516e-04,-9.83740524e-01,-8.82535702e-01,-8.07694999e-01,-1.00000000e+00 +9.50000000e+01,9.48250227e-01,6.23902760e-01,3.37150714e-01,2.72998783e+02,2.73260239e+02,2.81211324e+02,6.68381011e-06,4.90035826e-05,1.66187818e-04,-9.82745530e-01,-8.80852780e-01,-8.07098039e-01,-1.00000000e+00 +9.70000000e+01,9.45203001e-01,6.18847647e-01,3.33814241e-01,2.72998711e+02,2.73283330e+02,2.81431205e+02,7.03637859e-06,5.03522893e-05,1.67936462e-04,-9.81729520e-01,-8.79209358e-01,-8.06531155e-01,-1.00000000e+00 +9.90000000e+01,9.42095814e-01,6.13900999e-01,3.30572308e-01,2.72998638e+02,2.73307510e+02,2.81648705e+02,7.39443746e-06,5.16973957e-05,1.69667015e-04,-9.80693518e-01,-8.77604047e-01,-8.05992468e-01,-1.00000000e+00 +1.01000000e+02,9.38931682e-01,6.09058570e-01,3.27420277e-01,2.72998564e+02,2.73332774e+02,2.81863811e+02,7.75767990e-06,5.30386949e-05,1.71380016e-04,-9.79638530e-01,-8.76035511e-01,-8.05480227e-01,-1.00000000e+00 +1.21000000e+02,9.04771494e-01,5.65555530e-01,3.00006250e-01,2.72997761e+02,2.73642740e+02,2.83883437e+02,1.16140999e-05,6.62070325e-05,1.87645355e-04,-9.68248800e-01,-8.62112381e-01,-8.01503169e-01,-1.00000000e+00 +1.41000000e+02,8.67679000e-01,5.28964650e-01,2.78095786e-01,2.72996890e+02,2.74045668e+02,2.85675952e+02,1.57175300e-05,7.88573494e-05,2.02586864e-04,-9.55881383e-01,-8.50719663e-01,-7.98912672e-01,-1.00000000e+00 +1.61000000e+02,8.29419949e-01,4.97205792e-01,2.59818872e-01,2.72996000e+02,2.74521958e+02,2.87267662e+02,1.98999134e-05,9.09349842e-05,2.16457293e-04,-9.43125034e-01,-8.41140850e-01,-7.97085816e-01,-1.00000000e+00 +1.81000000e+02,7.91138054e-01,4.68940075e-01,2.44042324e-01,2.72995136e+02,2.75052019e+02,2.88686828e+02,2.40543938e-05,1.02422782e-04,2.29423407e-04,-9.30361150e-01,-8.32891007e-01,-7.95679873e-01,-1.00000000e+00 +2.01000000e+02,7.53540092e-01,4.43297873e-01,2.30055207e-01,2.72994359e+02,2.75618878e+02,2.89958936e+02,2.81157848e-05,1.13323580e-04,2.41599717e-04,-9.17825515e-01,-8.25640486e-01,-7.94502436e-01,-1.00000000e+00 +2.21000000e+02,7.17044327e-01,4.19707509e-01,2.17401544e-01,2.72993764e+02,2.76208827e+02,2.91105678e+02,3.20460180e-05,1.23651245e-04,2.53068066e-04,-9.05657843e-01,-8.19163384e-01,-7.93445825e-01,-1.00000000e+00 +2.41000000e+02,6.81883899e-01,3.97787047e-01,2.05783645e-01,2.72993494e+02,2.76811228e+02,2.92145121e+02,3.58242830e-05,1.33425662e-04,2.63889527e-04,-8.93936281e-01,-8.13303039e-01,-7.92451266e-01,-1.00000000e+00 +2.61000000e+02,6.48175072e-01,3.77275776e-01,1.95003554e-01,2.72993753e+02,2.77418032e+02,2.93092222e+02,3.94406011e-05,1.42669792e-04,2.74111793e-04,-8.82700201e-01,-8.07949077e-01,-7.91488365e-01,-1.00000000e+00 +2.81000000e+02,6.15961260e-01,3.57990771e-01,1.84926712e-01,2.72994809e+02,2.78023263e+02,2.93959369e+02,4.28917121e-05,1.51407906e-04,2.83773807e-04,-8.71964875e-01,-8.03022182e-01,-7.90543049e-01,-1.00000000e+00 +3.01000000e+02,5.85241108e-01,3.39799304e-01,1.75459132e-01,2.72996988e+02,2.78622559e+02,2.94756861e+02,4.61784654e-05,1.59664531e-04,2.92908679e-04,-8.61730829e-01,-7.98464027e-01,-7.89610400e-01,-1.00000000e+00 +3.21000000e+02,5.55986292e-01,3.22601316e-01,1.66532970e-01,2.73000663e+02,2.79212792e+02,2.95493304e+02,4.93041746e-05,1.67463826e-04,3.01545513e-04,-8.51989746e-01,-7.94230588e-01,-7.88690412e-01,-1.00000000e+00 +3.41000000e+02,5.28152765e-01,3.06318265e-01,1.58097376e-01,2.73006245e+02,2.79791774e+02,2.96175937e+02,5.22735841e-05,1.74829219e-04,3.09710574e-04,-8.42728169e-01,-7.90287714e-01,-7.87785488e-01,-1.00000000e+00 +3.61000000e+02,5.01687815e-01,2.90886007e-01,1.50112668e-01,2.73014162e+02,2.80358018e+02,2.96810883e+02,5.50922236e-05,1.81783214e-04,3.17428013e-04,-8.33929812e-01,-7.86608162e-01,-7.86898998e-01,-1.00000000e+00 +3.81000000e+02,4.76534484e-01,2.76250250e-01,1.42546625e-01,2.73024845e+02,2.80910572e+02,2.97403354e+02,5.77660062e-05,1.88347301e-04,3.24720333e-04,-8.25576985e-01,-7.83169610e-01,-7.86034481e-01,-1.00000000e+00 +4.01000000e+02,4.52634307e-01,2.62363633e-01,1.35372119e-01,2.73038717e+02,2.81448874e+02,2.97957814e+02,6.03009807e-05,1.94541924e-04,3.31608691e-04,-8.17651466e-01,-7.79953310e-01,-7.85195214e-01,-1.00000000e+00 +4.21000000e+02,4.29928985e-01,2.49183841e-01,1.28565585e-01,2.73056181e+02,2.81972659e+02,2.98478110e+02,6.27031805e-05,2.00386485e-04,3.38113089e-04,-8.10135026e-01,-7.76943167e-01,-7.84384016e-01,-1.00000000e+00 +4.41000000e+02,4.08361395e-01,2.36672375e-01,1.22106045e-01,2.73077606e+02,2.82481876e+02,2.98967571e+02,6.49785323e-05,2.05899374e-04,3.44252506e-04,-8.03009736e-01,-7.74125106e-01,-7.83603173e-01,-1.00000000e+00 +4.61000000e+02,3.87876176e-01,2.24793745e-01,1.15974464e-01,2.73103328e+02,2.82976635e+02,2.99429094e+02,6.71328032e-05,2.11098004e-04,3.50044990e-04,-7.96258147e-01,-7.71486618e-01,-7.82854440e-01,-1.00000000e+00 +4.81000000e+02,3.68420045e-01,2.13514923e-01,1.10153331e-01,2.73133641e+02,2.83457158e+02,2.99865217e+02,6.91715715e-05,2.15998860e-04,3.55507722e-04,-7.89863376e-01,-7.69016452e-01,-7.82139078e-01,-1.00000000e+00 +5.01000000e+02,3.49941956e-01,2.02804973e-01,1.04626376e-01,2.73168794e+02,2.83923752e+02,3.00278167e+02,7.11002110e-05,2.20617543e-04,3.60657069e-04,-7.83809159e-01,-7.66704376e-01,-7.81457904e-01,-1.00000000e+00 +5.21000000e+02,3.32393148e-01,1.92634789e-01,9.93783757e-02,2.73208995e+02,2.84376777e+02,3.00669909e+02,7.29238852e-05,2.24968817e-04,3.65508627e-04,-7.78079860e-01,-7.64541014e-01,-7.80811355e-01,-1.00000000e+00 +5.41000000e+02,3.15727134e-01,1.82976902e-01,9.43950195e-02,2.73254406e+02,2.84816633e+02,3.01042188e+02,7.46475468e-05,2.29066655e-04,3.70077256e-04,-7.72660475e-01,-7.62517714e-01,-7.80199539e-01,-1.00000000e+00 +5.61000000e+02,2.99899657e-01,1.73805336e-01,8.96628112e-02,2.73305150e+02,2.85243741e+02,3.01396553e+02,7.62759404e-05,2.32924287e-04,3.74377113e-04,-7.67536616e-01,-7.60626454e-01,-7.79622295e-01,-1.00000000e+00 +5.81000000e+02,2.84868614e-01,1.65095490e-01,8.51689933e-02,2.73361310e+02,2.85658538e+02,3.01734388e+02,7.78136066e-05,2.36554234e-04,3.78421684e-04,-7.62694499e-01,-7.58859762e-01,-7.79079239e-01,-1.00000000e+00 +6.01000000e+02,2.70593978e-01,1.56824046e-01,8.09014892e-02,2.73422934e+02,2.86061464e+02,3.02056933e+02,7.92648886e-05,2.39968355e-04,3.82223811e-04,-7.58120921e-01,-7.57210660e-01,-7.78569804e-01,-1.00000000e+00 +6.21000000e+02,2.57037708e-01,1.48968888e-01,7.68488552e-02,2.73490036e+02,2.86452957e+02,3.02365303e+02,8.06339377e-05,2.43177882e-04,3.85795723e-04,-7.53803245e-01,-7.55672609e-01,-7.78093279e-01,-1.00000000e+00 +6.41000000e+02,2.44163661e-01,1.41509028e-01,7.30002396e-02,2.73562599e+02,2.86833452e+02,3.02660503e+02,8.19247206e-05,2.46193451e-04,3.89149057e-04,-7.49729376e-01,-7.54239472e-01,-7.77648842e-01,-1.00000000e+00 +6.61000000e+02,2.31937503e-01,1.34424550e-01,6.93453479e-02,2.73640584e+02,2.87203374e+02,3.02943442e+02,8.31410258e-05,2.49025142e-04,3.92294892e-04,-7.45887743e-01,-7.52905480e-01,-7.77235580e-01,-1.00000000e+00 +6.81000000e+02,2.20326621e-01,1.27696544e-01,6.58744104e-02,2.73723923e+02,2.87563139e+02,3.03214944e+02,8.42864706e-05,2.51682505e-04,3.95243766e-04,-7.42267281e-01,-7.51665200e-01,-7.76852517e-01,-1.00000000e+00 +7.01000000e+02,2.09300041e-01,1.21307058e-01,6.25781544e-02,2.73812534e+02,2.87913148e+02,3.03475759e+02,8.53645075e-05,2.54174592e-04,3.98005703e-04,-7.38857409e-01,-7.50513510e-01,-7.76498631e-01,-1.00000000e+00 +7.21000000e+02,1.98828346e-01,1.15239046e-01,5.94477770e-02,2.73906312e+02,2.88253792e+02,3.03726571e+02,8.63784308e-05,2.56509980e-04,4.00590240e-04,-7.35648013e-01,-7.49445578e-01,-7.76172870e-01,-1.00000000e+00 +7.41000000e+02,1.88883597e-01,1.09476323e-01,5.64749206e-02,2.74005142e+02,2.88585447e+02,3.03968004e+02,8.73313834e-05,2.58696802e-04,4.03006443e-04,-7.32629429e-01,-7.48456843e-01,-7.75874163e-01,-1.00000000e+00 +7.61000000e+02,1.79439264e-01,1.04003518e-01,5.36516500e-02,2.74108894e+02,2.88908474e+02,3.04200630e+02,8.82263623e-05,2.60742770e-04,4.05262934e-04,-7.29792427e-01,-7.47542994e-01,-7.75601432e-01,-1.00000000e+00 +7.81000000e+02,1.70470149e-01,9.88060363e-02,5.09704307e-02,2.74217429e+02,2.89223223e+02,3.04424974e+02,8.90662254e-05,2.62655193e-04,4.07367909e-04,-7.27128189e-01,-7.46699959e-01,-7.75353602e-01,-1.00000000e+00 +8.01000000e+02,1.61952325e-01,9.38700167e-02,4.84241082e-02,2.74330601e+02,2.89530027e+02,3.04641521e+02,8.98536967e-05,2.64441004e-04,4.09329158e-04,-7.24628300e-01,-7.45923883e-01,-7.75129605e-01,-1.00000000e+00 +8.21000000e+02,1.53863069e-01,8.91822973e-02,4.60058889e-02,2.74448256e+02,2.89829208e+02,3.04850716e+02,9.05913726e-05,2.66106780e-04,4.11154085e-04,-7.22284725e-01,-7.45211124e-01,-7.74928390e-01,-1.00000000e+00 +8.41000000e+02,1.46180802e-01,8.47303786e-02,4.37093218e-02,2.74570237e+02,2.90121072e+02,3.05052971e+02,9.12817269e-05,2.67658756e-04,4.12849727e-04,-7.20089798e-01,-7.44558233e-01,-7.74748924e-01,-1.00000000e+00 +8.61000000e+02,1.38885031e-01,8.05023906e-02,4.15282809e-02,2.74696383e+02,2.90405913e+02,3.05248667e+02,9.19271165e-05,2.69102848e-04,4.14422771e-04,-7.18036208e-01,-7.43961947e-01,-7.74590196e-01,-1.00000000e+00 +8.81000000e+02,1.31956297e-01,7.64870606e-02,3.94569490e-02,2.74826533e+02,2.90684015e+02,3.05438157e+02,9.25297863e-05,2.70444668e-04,4.15879570e-04,-7.16116978e-01,-7.43419176e-01,-7.74451224e-01,-1.00000000e+00 +9.01000000e+02,1.25376119e-01,7.26736838e-02,3.74898020e-02,2.74960522e+02,2.90955645e+02,3.05621770e+02,9.30918741e-05,2.71689540e-04,4.17226159e-04,-7.14325460e-01,-7.42926996e-01,-7.74331051e-01,-1.00000000e+00 +9.21000000e+02,1.19126947e-01,6.90520941e-02,3.56215941e-02,2.75098188e+02,2.91221063e+02,3.05799813e+02,9.36154158e-05,2.72842513e-04,4.18468274e-04,-7.12655314e-01,-7.42482637e-01,-7.74228752e-01,-1.00000000e+00 +9.41000000e+02,1.13192116e-01,6.56126375e-02,3.38473440e-02,2.75239368e+02,2.91480515e+02,3.05972568e+02,9.41023495e-05,2.73908382e-04,4.19611360e-04,-7.11100497e-01,-7.42083477e-01,-7.74143430e-01,-1.00000000e+00 +9.61000000e+02,1.07555799e-01,6.23461456e-02,3.21623209e-02,2.75383902e+02,2.91734237e+02,3.06140304e+02,9.45545201e-05,2.74891693e-04,4.20660592e-04,-7.09655253e-01,-7.41727032e-01,-7.74074224e-01,-1.00000000e+00 +9.81000000e+02,1.02202966e-01,5.92439119e-02,3.05620328e-02,2.75531632e+02,2.91982454e+02,3.06303267e+02,9.49736837e-05,2.75796761e-04,4.21620883e-04,-7.08314095e-01,-7.41410950e-01,-7.74020300e-01,-1.00000000e+00 +1.00100000e+03,9.71193437e-02,5.62976681e-02,2.90422136e-02,2.75682402e+02,2.92225383e+02,3.06461691e+02,9.53615114e-05,2.76627682e-04,4.22496902e-04,-7.07071800e-01,-7.41133003e-01,-7.73980860e-01,-1.00000000e+00 +1.20100000e+03,6.17118722e-02,3.57756987e-02,1.84559248e-02,2.77357431e+02,2.94327959e+02,3.07778070e+02,9.74542596e-05,2.80750449e-04,4.26509311e-04,-7.00294045e-01,-7.40465143e-01,-7.74390638e-01,-1.00000000e+00 +1.40100000e+03,3.93521047e-02,2.28149956e-02,1.17700312e-02,2.79155190e+02,2.96123987e+02,3.08854214e+02,9.80581618e-05,2.81420073e-04,4.26622771e-04,-6.98233959e-01,-7.41489484e-01,-7.75435915e-01,-1.00000000e+00 +1.60100000e+03,2.52251142e-02,1.46256569e-02,7.54539590e-03,2.81003845e+02,2.97697029e+02,3.09763430e+02,9.77500581e-05,2.80003826e-04,4.24399689e-04,-6.99070970e-01,-7.43490676e-01,-7.76832814e-01,-1.00000000e+00 +1.80100000e+03,1.62931416e-02,9.44737913e-03,4.87400545e-03,2.82855572e+02,2.99103489e+02,3.10553222e+02,9.68954471e-05,2.77359754e-04,4.20818195e-04,-7.01648962e-01,-7.46028846e-01,-7.78406441e-01,-1.00000000e+00 +2.00100000e+03,1.06396278e-02,6.16947687e-03,3.18293763e-03,2.84679484e+02,3.00381726e+02,3.11254989e+02,9.57251980e-05,2.74025125e-04,4.16488460e-04,-7.05236379e-01,-7.48834922e-01,-7.80051049e-01,-1.00000000e+00 +2.20100000e+03,7.05537438e-03,4.09113866e-03,2.11069015e-03,2.86456274e+02,3.01558076e+02,3.11889800e+02,9.43845731e-05,2.70335051e-04,4.11789191e-04,-7.09372211e-01,-7.51745616e-01,-7.81703937e-01,-1.00000000e+00 +2.40100000e+03,4.77745750e-03,2.77016244e-03,1.42915564e-03,2.88174411e+02,3.02650851e+02,3.12471985e+02,9.29644667e-05,2.66497248e-04,4.06953657e-04,-7.13767546e-01,-7.54663035e-01,-7.83328854e-01,-1.00000000e+00 +2.60100000e+03,3.32450908e-03,1.92751619e-03,9.94393423e-04,2.89827527e+02,3.03672992e+02,3.13011436e+02,9.15212797e-05,2.62639191e-04,4.02123860e-04,-7.18242808e-01,-7.57529625e-01,-7.84905619e-01,-1.00000000e+00 +2.80100000e+03,2.39279629e-03,1.38712098e-03,7.15569232e-04,2.91412636e+02,3.04633833e+02,3.13515089e+02,9.00895527e-05,2.58837858e-04,3.97384627e-04,-7.22687796e-01,-7.60312626e-01,-7.86423728e-01,-1.00000000e+00 +3.00100000e+03,1.79065801e-03,1.03785729e-03,5.35357723e-04,2.92928955e+02,3.05540273e+02,3.13987911e+02,8.86899929e-05,2.55138510e-04,3.92785044e-04,-7.27036264e-01,-7.62994457e-01,-7.87878425e-01,-1.00000000e+00 +3.20100000e+03,1.39714072e-03,8.09592888e-04,4.17577192e-04,2.94377102e+02,3.06397562e+02,3.14433544e+02,8.73345720e-05,2.51566537e-04,3.88351934e-04,-7.31249752e-01,-7.65566764e-01,-7.89268333e-01,-1.00000000e+00 +3.40100000e+03,1.13590186e-03,6.58057681e-04,3.39387528e-04,2.95758579e+02,3.07209821e+02,3.14854731e+02,8.60297641e-05,2.48134925e-04,3.84098326e-04,-7.35307314e-01,-7.68026750e-01,-7.90594009e-01,-1.00000000e+00 +3.60100000e+03,9.58745372e-04,5.55300352e-04,2.86367348e-04,2.97075418e+02,3.07980387e+02,3.15253608e+02,8.47786009e-05,2.44848980e-04,3.80028761e-04,-7.39198982e-01,-7.70374923e-01,-7.91857089e-01,-1.00000000e+00 +3.80100000e+03,8.35234252e-04,4.83666959e-04,2.49407783e-04,2.98329958e+02,3.08712049e+02,3.15631884e+02,8.35819758e-05,2.41709290e-04,3.76142630e-04,-7.42921622e-01,-7.72613701e-01,-7.93059777e-01,-1.00000000e+00 +4.00100000e+03,7.46137044e-04,4.32001645e-04,2.22752523e-04,2.99524688e+02,3.09407205e+02,3.15990980e+02,8.24394715e-05,2.38713597e-04,3.72436254e-04,-7.46476299e-01,-7.74746571e-01,-7.94204544e-01,-1.00000000e+00 +4.20100000e+03,6.79293206e-04,3.93249416e-04,2.02761097e-04,3.00662158e+02,3.10067965e+02,3.16332103e+02,8.13498834e-05,2.35857976e-04,3.68904187e-04,-7.49866604e-01,-7.76777574e-01,-7.95293955e-01,-1.00000000e+00 +4.40100000e+03,6.27004508e-04,3.62943553e-04,1.87128512e-04,3.01744912e+02,3.10696223e+02,3.16656310e+02,8.03115510e-05,2.33137564e-04,3.65540027e-04,-7.53097602e-01,-7.78710993e-01,-7.96330570e-01,-1.00000000e+00 +4.60100000e+03,5.84389226e-04,3.38251187e-04,1.74392836e-04,3.02775447e+02,3.11293705e+02,3.16964545e+02,7.93225663e-05,2.30547023e-04,3.62336916e-04,-7.56175161e-01,-7.80551169e-01,-7.97316891e-01,-1.00000000e+00 +4.80100000e+03,5.48343749e-04,3.17371068e-04,1.63624471e-04,3.03756193e+02,3.11862004e+02,3.17257664e+02,7.83809049e-05,2.28080828e-04,3.59287848e-04,-7.59105537e-01,-7.82302391e-01,-7.98255334e-01,-1.00000000e+00 +5.00100000e+03,5.16887409e-04,2.99153530e-04,1.54230066e-04,3.04689497e+02,3.12402602e+02,3.17536456e+02,7.74845069e-05,2.25733433e-04,3.56385855e-04,-7.61895107e-01,-7.83968839e-01,-7.99148215e-01,-1.00000000e+00 +5.20100000e+03,4.88749082e-04,2.82860682e-04,1.45828775e-04,3.05577614e+02,3.12916887e+02,3.17801651e+02,7.66313274e-05,2.23499385e-04,3.53624118e-04,-7.64550217e-01,-7.85554543e-01,-7.99997746e-01,-1.00000000e+00 +5.40100000e+03,4.63106335e-04,2.68015056e-04,1.38174155e-04,3.06422706e+02,3.13406164e+02,3.18053933e+02,7.58193663e-05,2.21373381e-04,3.50996031e-04,-7.67077078e-01,-7.87063376e-01,-8.00806034e-01,-1.00000000e+00 +5.60100000e+03,4.39420794e-04,2.54304081e-04,1.31104873e-04,3.07226841e+02,3.13871666e+02,3.18293944e+02,7.50466856e-05,2.19350302e-04,3.48495234e-04,-7.69481711e-01,-7.88499040e-01,-8.01575087e-01,-1.00000000e+00 +5.80100000e+03,4.17334243e-04,2.41519785e-04,1.24513582e-04,3.07991991e+02,3.14314560e+02,3.18522291e+02,7.43114194e-05,2.17425234e-04,3.46115631e-04,-7.71769919e-01,-7.89865071e-01,-8.02306813e-01,-1.00000000e+00 +6.00100000e+03,3.96603024e-04,2.29520714e-04,1.18327272e-04,3.08720040e+02,3.14735950e+02,3.18739546e+02,7.36117777e-05,2.15593471e-04,3.43851388e-04,-7.73947266e-01,-7.91164837e-01,-8.03003027e-01,-1.00000000e+00 +6.20100000e+03,3.77056606e-04,2.18207885e-04,1.12494859e-04,3.09412780e+02,3.15136887e+02,3.18946252e+02,7.29460487e-05,2.13850518e-04,3.41696941e-04,-7.76019079e-01,-7.92401551e-01,-8.03665455e-01,-1.00000000e+00 +6.40100000e+03,3.58571410e-04,2.07509584e-04,1.06979332e-04,3.10071921e+02,3.15518365e+02,3.19142923e+02,7.23125975e-05,2.12192087e-04,3.39646982e-04,-7.77990443e-01,-7.93578269e-01,-8.04295739e-01,-1.00000000e+00 +6.60100000e+03,3.41054251e-04,1.97371753e-04,1.01752799e-04,3.10699091e+02,3.15881331e+02,3.19330050e+02,7.17098649e-05,2.10614090e-04,3.37696452e-04,-7.79866210e-01,-7.94697900e-01,-8.04895442e-01,-1.00000000e+00 +6.80100000e+03,3.24431859e-04,1.87751904e-04,9.67933375e-05,3.11295838e+02,3.16226686e+02,3.19508096e+02,7.11363649e-05,2.09112633e-04,3.35840537e-04,-7.81651003e-01,-7.95763213e-01,-8.05466047e-01,-1.00000000e+00 +7.00100000e+03,3.08644232e-04,1.78615254e-04,9.20830054e-05,3.11863637e+02,3.16555284e+02,3.19677502e+02,7.05906821e-05,2.07684007e-04,3.34074651e-04,-7.83349226e-01,-7.96776844e-01,-8.06008967e-01,-1.00000000e+00 +7.20100000e+03,2.93640401e-04,1.69932270e-04,8.76065698e-05,3.12403892e+02,3.16867939e+02,3.19838689e+02,7.00714689e-05,2.06324682e-04,3.32394428e-04,-7.84965075e-01,-7.97741301e-01,-8.06525547e-01,-1.00000000e+00 +7.40100000e+03,2.79375738e-04,1.61677098e-04,8.33506974e-05,3.12917939e+02,3.17165426e+02,3.19992055e+02,6.95774422e-05,2.05031299e-04,3.30795715e-04,-7.86502540e-01,-7.98658967e-01,-8.07017065e-01,-1.00000000e+00 +7.60100000e+03,2.65810217e-04,1.53826559e-04,7.93034348e-05,3.13407049e+02,3.17448480e+02,3.20137980e+02,6.91073808e-05,2.03800659e-04,3.29274558e-04,-7.87965423e-01,-7.99532115e-01,-8.07484737e-01,-1.00000000e+00 +7.80100000e+03,2.52907291e-04,1.46359492e-04,7.54538707e-05,3.13872431e+02,3.17717802e+02,3.20276826e+02,6.86601227e-05,2.02629720e-04,3.27827197e-04,-7.89357340e-01,-8.00362903e-01,-8.07929720e-01,-1.00000000e+00 +8.00100000e+03,2.40633159e-04,1.39256327e-04,7.17919149e-05,3.14315236e+02,3.17974058e+02,3.20408936e+02,6.82345619e-05,2.01515586e-04,3.26450050e-04,-7.90681733e-01,-8.01153387e-01,-8.08353116e-01,-1.00000000e+00 +8.20100000e+03,2.28956269e-04,1.32498801e-04,6.83081508e-05,3.14736559e+02,3.18217881e+02,3.20534637e+02,6.78296459e-05,2.00455502e-04,3.25139713e-04,-7.91941876e-01,-8.01905522e-01,-8.08755970e-01,-1.00000000e+00 +8.40100000e+03,2.17846987e-04,1.26069761e-04,6.49937349e-05,3.15137443e+02,3.18449877e+02,3.20654239e+02,6.74443734e-05,1.99446844e-04,3.23892944e-04,-7.93140887e-01,-8.02621169e-01,-8.09139281e-01,-1.00000000e+00 +8.60100000e+03,2.07277361e-04,1.19953028e-04,6.18403261e-05,3.15518879e+02,3.18670617e+02,3.20768039e+02,6.70777915e-05,1.98487120e-04,3.22706659e-04,-7.94281731e-01,-8.03302097e-01,-8.09503996e-01,-1.00000000e+00 +8.80100000e+03,1.97220949e-04,1.14133298e-04,5.88400343e-05,3.15881811e+02,3.18880648e+02,3.20876318e+02,6.67289934e-05,1.97573955e-04,3.21577924e-04,-7.95367230e-01,-8.03949991e-01,-8.09851018e-01,-1.00000000e+00 +9.00100000e+03,1.87652686e-04,1.08596067e-04,5.59853812e-05,3.16227135e+02,3.19080491e+02,3.20979345e+02,6.63971165e-05,1.96705090e-04,3.20503948e-04,-7.96400067e-01,-8.04566454e-01,-8.10181205e-01,-1.00000000e+00 +9.20100000e+03,1.78548787e-04,1.03327567e-04,5.32692700e-05,3.16555707e+02,3.19270638e+02,3.21077373e+02,6.60813400e-05,1.95878376e-04,3.19482073e-04,-7.97382799e-01,-8.05153011e-01,-8.10495373e-01,-1.00000000e+00 +9.40100000e+03,1.69886657e-04,9.83147246e-05,5.06849596e-05,3.16868339e+02,3.19451560e+02,3.21170645e+02,6.57808826e-05,1.95091768e-04,3.18509772e-04,-7.98317856e-01,-8.05711112e-01,-8.10794299e-01,-1.00000000e+00 +9.60100000e+03,1.61644825e-04,9.35451115e-05,4.82260435e-05,3.17165804e+02,3.19623706e+02,3.21259393e+02,6.54950014e-05,1.94343321e-04,3.17584641e-04,-7.99207551e-01,-8.06242137e-01,-8.11078724e-01,-1.00000000e+00 +9.80100000e+03,1.53802874e-04,8.90069129e-05,4.58864303e-05,3.17448839e+02,3.19787500e+02,3.21343835e+02,6.52229891e-05,1.93631184e-04,3.16704391e-04,-8.00054083e-01,-8.06747401e-01,-8.11349351e-01,-1.00000000e+00 +1.00010000e+04,1.46341387e-04,8.46888925e-05,4.36603273e-05,3.17718142e+02,3.19943348e+02,3.21424180e+02,6.49641729e-05,1.92953594e-04,3.15866844e-04,-8.00859548e-01,-8.07228153e-01,-8.11606849e-01,-1.00000000e+00 +1.20010000e+04,9.33947592e-05,5.40482687e-05,2.78639267e-05,3.19629128e+02,3.21049250e+02,3.21994315e+02,6.31276052e-05,1.88145397e-04,3.09923586e-04,-8.06575152e-01,-8.10639585e-01,-8.13434061e-01,-1.00000000e+00 +1.40010000e+04,5.96043655e-05,3.44935067e-05,1.77827073e-05,3.20848715e+02,3.21755035e+02,3.22358173e+02,6.19555109e-05,1.85076814e-04,3.06130609e-04,-8.10222841e-01,-8.12816754e-01,-8.14600184e-01,-1.00000000e+00 +1.60010000e+04,3.80394003e-05,2.20136947e-05,1.13488922e-05,3.21627054e+02,3.22205465e+02,3.22590387e+02,6.12074824e-05,1.83118450e-04,3.03709938e-04,-8.12550789e-01,-8.14206220e-01,-8.15344402e-01,-1.00000000e+00 +1.80010000e+04,2.42766781e-05,1.40491011e-05,7.24284295e-06,3.22123788e+02,3.22492929e+02,3.22738586e+02,6.07300920e-05,1.81868626e-04,3.02165070e-04,-8.14036481e-01,-8.15092975e-01,-8.15819360e-01,-1.00000000e+00 +2.00010000e+04,1.54933331e-05,8.96611153e-06,4.62236959e-06,3.22440803e+02,3.22676388e+02,3.22833166e+02,6.04254222e-05,1.81070990e-04,3.01179138e-04,-8.14984648e-01,-8.15658900e-01,-8.16122478e-01,-1.00000000e+00 +2.20010000e+04,9.88781790e-06,5.72215658e-06,2.94998813e-06,3.22643122e+02,3.22793472e+02,3.22893527e+02,6.02309825e-05,1.80561940e-04,3.00549918e-04,-8.15589766e-01,-8.16020072e-01,-8.16315927e-01,-1.00000000e+00 +2.40010000e+04,6.31038797e-06,3.65187025e-06,1.88267723e-06,3.22772241e+02,3.22868194e+02,3.22932049e+02,6.01068914e-05,1.80237066e-04,3.00148350e-04,-8.15975951e-01,-8.16250572e-01,-8.16439386e-01,-1.00000000e+00 +2.60010000e+04,4.02727849e-06,2.33061716e-06,1.20152129e-06,3.22854644e+02,3.22915882e+02,3.22956634e+02,6.00276967e-05,1.80029731e-04,2.99892071e-04,-8.16222414e-01,-8.16397677e-01,-8.16518177e-01,-1.00000000e+00 +2.80010000e+04,2.57020204e-06,1.48739577e-06,7.66808772e-07,3.22907234e+02,3.22946316e+02,3.22972324e+02,5.99771548e-05,1.79897411e-04,2.99728514e-04,-8.16379706e-01,-8.16491558e-01,-8.16568462e-01,-1.00000000e+00 +3.00010000e+04,1.64029842e-06,9.49253362e-07,4.89376009e-07,3.22940797e+02,3.22965739e+02,3.22982337e+02,5.99448990e-05,1.79812964e-04,2.99624132e-04,-8.16480090e-01,-8.16551474e-01,-8.16600553e-01,-1.00000000e+00 +3.20010000e+04,1.04683556e-06,6.05811825e-07,3.12318907e-07,3.22962217e+02,3.22978135e+02,3.22988728e+02,5.99243135e-05,1.79759071e-04,2.99557516e-04,-8.16544154e-01,-8.16589711e-01,-8.16621034e-01,-1.00000000e+00 +3.40010000e+04,6.68088614e-07,3.86628040e-07,1.99321377e-07,3.22975887e+02,3.22986046e+02,3.22992806e+02,5.99111758e-05,1.79724676e-04,2.99515001e-04,-8.16585040e-01,-8.16614115e-01,-8.16634105e-01,-1.00000000e+00 +3.60010000e+04,4.26372978e-07,2.46745335e-07,1.27206552e-07,3.22984611e+02,3.22991094e+02,3.22995409e+02,5.99027913e-05,1.79702725e-04,2.99487869e-04,-8.16611133e-01,-8.16629689e-01,-8.16642446e-01,-1.00000000e+00 +3.80010000e+04,2.72110484e-07,1.57472438e-07,8.11829976e-08,3.22990179e+02,3.22994316e+02,3.22997070e+02,5.98974404e-05,1.79688716e-04,2.99470553e-04,-8.16627786e-01,-8.16639628e-01,-8.16647770e-01,-1.00000000e+00 +4.00010000e+04,1.73660432e-07,1.00498633e-07,5.18108462e-08,3.22993732e+02,3.22996373e+02,3.22998130e+02,5.98940255e-05,1.79679775e-04,2.99459502e-04,-8.16638414e-01,-8.16645971e-01,-8.16651168e-01,-1.00000000e+00 +4.20010000e+04,1.10829782e-07,6.41380507e-08,3.30655909e-08,3.22996000e+02,3.22997685e+02,3.22998807e+02,5.98918460e-05,1.79674070e-04,2.99452449e-04,-8.16645197e-01,-8.16650020e-01,-8.16653336e-01,-1.00000000e+00 +4.40010000e+04,7.07313719e-08,4.09327911e-08,2.11024019e-08,3.22997447e+02,3.22998523e+02,3.22999238e+02,5.98904551e-05,1.79670428e-04,2.99447948e-04,-8.16649525e-01,-8.16652603e-01,-8.16654720e-01,-1.00000000e+00 +4.60010000e+04,4.51406371e-08,2.61232353e-08,1.34675158e-08,3.22998371e+02,3.22999057e+02,3.22999514e+02,5.98895675e-05,1.79668104e-04,2.99445075e-04,-8.16652288e-01,-8.16654252e-01,-8.16655603e-01,-1.00000000e+00 +4.80010000e+04,2.88086751e-08,1.66718027e-08,8.59494484e-09,3.22998960e+02,3.22999398e+02,3.22999690e+02,5.98890010e-05,1.79666621e-04,2.99443242e-04,-8.16654051e-01,-8.16655304e-01,-8.16656166e-01,-1.00000000e+00 +5.00010000e+04,1.83856458e-08,1.06399152e-08,5.48527869e-09,3.22999336e+02,3.22999616e+02,3.22999802e+02,5.98886394e-05,1.79665675e-04,2.99442072e-04,-8.16655176e-01,-8.16655976e-01,-8.16656526e-01,-1.00000000e+00 +5.20010000e+04,1.17336868e-08,6.79037513e-09,3.50069519e-09,3.22999576e+02,3.22999755e+02,3.22999874e+02,5.98884087e-05,1.79665071e-04,2.99441325e-04,-8.16655894e-01,-8.16656405e-01,-8.16656756e-01,-1.00000000e+00 +5.40010000e+04,7.48841827e-09,4.33360548e-09,2.23413752e-09,3.22999730e+02,3.22999844e+02,3.22999919e+02,5.98882614e-05,1.79664685e-04,2.99440849e-04,-8.16656352e-01,-8.16656678e-01,-8.16656902e-01,-1.00000000e+00 +5.60010000e+04,4.77909541e-09,2.76569942e-09,1.42582265e-09,3.22999828e+02,3.22999900e+02,3.22999949e+02,5.98881674e-05,1.79664439e-04,2.99440545e-04,-8.16656645e-01,-8.16656853e-01,-8.16656996e-01,-1.00000000e+00 +5.80010000e+04,3.05001039e-09,1.76506456e-09,9.09957534e-10,3.22999890e+02,3.22999936e+02,3.22999967e+02,5.98881075e-05,1.79664282e-04,2.99440351e-04,-8.16656831e-01,-8.16656964e-01,-8.16657055e-01,-1.00000000e+00 +6.00010000e+04,1.94651114e-09,1.12646103e-09,5.80733264e-10,3.22999930e+02,3.22999959e+02,3.22999979e+02,5.98880692e-05,1.79664182e-04,2.99440227e-04,-8.16656951e-01,-8.16657035e-01,-8.16657093e-01,-1.00000000e+00 +6.20010000e+04,1.24226002e-09,7.18905470e-10,3.70622960e-10,3.22999955e+02,3.22999974e+02,3.22999987e+02,5.98880447e-05,1.79664118e-04,2.99440148e-04,-8.16657027e-01,-8.16657081e-01,-8.16657118e-01,-1.00000000e+00 +6.40010000e+04,7.92808172e-10,4.58804221e-10,2.36530931e-10,3.22999971e+02,3.22999983e+02,3.22999991e+02,5.98880292e-05,1.79664077e-04,2.99440097e-04,-8.16657075e-01,-8.16657110e-01,-8.16657133e-01,-1.00000000e+00 +6.60010000e+04,5.05968714e-10,2.92808001e-10,1.50953598e-10,3.22999982e+02,3.22999989e+02,3.22999995e+02,5.98880192e-05,1.79664051e-04,2.99440065e-04,-8.16657106e-01,-8.16657128e-01,-8.16657143e-01,-1.00000000e+00 +6.80010000e+04,3.22908257e-10,1.86869493e-10,9.63382893e-11,3.22999988e+02,3.22999993e+02,3.22999997e+02,5.98880129e-05,1.79664034e-04,2.99440045e-04,-8.16657126e-01,-8.16657140e-01,-8.16657149e-01,-1.00000000e+00 +7.00010000e+04,2.06079631e-10,1.19259869e-10,6.14829699e-11,3.22999993e+02,3.22999996e+02,3.22999998e+02,5.98880088e-05,1.79664024e-04,2.99440031e-04,-8.16657138e-01,-8.16657147e-01,-8.16657154e-01,-1.00000000e+00 +7.20010000e+04,1.31519666e-10,7.61114449e-11,3.92383289e-11,3.22999995e+02,3.22999997e+02,3.22999999e+02,5.98880062e-05,1.79664017e-04,2.99440023e-04,-8.16657146e-01,-8.16657152e-01,-8.16657156e-01,-1.00000000e+00 +7.40010000e+04,8.39355610e-11,4.85741606e-11,2.50418163e-11,3.22999997e+02,3.22999998e+02,3.22999999e+02,5.98880046e-05,1.79664013e-04,2.99440018e-04,-8.16657152e-01,-8.16657155e-01,-8.16657158e-01,-1.00000000e+00 +7.60010000e+04,5.35675402e-11,3.09999480e-11,1.59816455e-11,3.22999998e+02,3.22999999e+02,3.22999999e+02,5.98880035e-05,1.79664010e-04,2.99440014e-04,-8.16657155e-01,-8.16657157e-01,-8.16657159e-01,-1.00000000e+00 +7.80010000e+04,3.41867217e-11,1.97841229e-11,1.01994647e-11,3.22999999e+02,3.22999999e+02,3.23000000e+02,5.98880028e-05,1.79664008e-04,2.99440012e-04,-8.16657157e-01,-8.16657158e-01,-8.16657159e-01,-1.00000000e+00 +8.00010000e+04,2.18178806e-11,1.26261747e-11,6.50927039e-12,3.22999999e+02,3.23000000e+02,3.23000000e+02,5.98880024e-05,1.79664007e-04,2.99440011e-04,-8.16657158e-01,-8.16657159e-01,-8.16657160e-01,-1.00000000e+00 +8.20010000e+04,1.39241861e-11,8.05803191e-12,4.15421838e-12,3.22999999e+02,3.23000000e+02,3.23000000e+02,5.98880021e-05,1.79664006e-04,2.99440010e-04,-8.16657159e-01,-8.16657160e-01,-8.16657160e-01,-1.00000000e+00 +8.40010000e+04,8.88640578e-12,5.14264236e-12,2.65122849e-12,3.23000000e+02,3.23000000e+02,3.23000000e+02,5.98880020e-05,1.79664006e-04,2.99440009e-04,-8.16657160e-01,-8.16657160e-01,-8.16657160e-01,-1.00000000e+00 +8.60010000e+04,5.67121843e-12,3.28197650e-12,1.69198394e-12,3.23000000e+02,3.23000000e+02,3.23000000e+02,5.98880019e-05,1.79664005e-04,2.99440009e-04,-8.16657160e-01,-8.16657160e-01,-8.16657160e-01,-1.00000000e+00 +8.80010000e+04,3.61933802e-12,2.09453453e-12,1.07981039e-12,3.23000000e+02,3.23000000e+02,3.23000000e+02,5.98880018e-05,1.79664005e-04,2.99440009e-04,-8.16657160e-01,-8.16657160e-01,-8.16657161e-01,-1.00000000e+00 +9.00010000e+04,2.30991390e-12,1.33676211e-12,6.89149977e-13,3.23000000e+02,3.23000000e+02,3.23000000e+02,5.98880017e-05,1.79664005e-04,2.99440009e-04,-8.16657160e-01,-8.16657160e-01,-8.16657161e-01,-1.00000000e+00 +9.20010000e+04,1.47410565e-12,8.53068065e-13,4.39785839e-13,3.23000000e+02,3.23000000e+02,3.23000000e+02,5.98880017e-05,1.79664005e-04,2.99440008e-04,-8.16657160e-01,-8.16657161e-01,-8.16657161e-01,-1.00000000e+00 +9.40010000e+04,9.40858832e-13,5.44481789e-13,2.80703127e-13,3.23000000e+02,3.23000000e+02,3.23000000e+02,5.98880017e-05,1.79664005e-04,2.99440008e-04,-8.16657161e-01,-8.16657161e-01,-8.16657161e-01,-1.00000000e+00 +9.60010000e+04,6.00438086e-13,3.47477670e-13,1.79137677e-13,3.23000000e+02,3.23000000e+02,3.23000000e+02,5.98880017e-05,1.79664005e-04,2.99440008e-04,-8.16657161e-01,-8.16657161e-01,-8.16657161e-01,-1.00000000e+00 +9.80010000e+04,3.83251259e-13,2.21794598e-13,1.14343826e-13,3.23000000e+02,3.23000000e+02,3.23000000e+02,5.98880017e-05,1.79664005e-04,2.99440008e-04,-8.16657161e-01,-8.16657161e-01,-8.16657161e-01,-1.00000000e+00 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_1p4m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_1p4m.csv deleted file mode 100644 index fe4937dfb4d..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_1p4m.csv +++ /dev/null @@ -1,64 +0,0 @@ -0.010153382770450296, 0.0000013599017727036408 -0.013363773135962386, 0.000001348306210454211 -0.017588889122929694, 0.0000020169836335079535 -0.022454281034521375, 0.00000200667646706404 -0.028665524766994802, 0.000001996369300620018 -0.03715929564047002, 6.24871965666839e-7 -0.049658839683139715, 0.00000197317817612105 -0.06947524531581625, 0.0000019590058222605875 -0.08869330361688064, 0.0000019486986558165656 -0.12990047589490736, 0.0000019325937082477745 -0.1608473432635658, 0.0000019235749376093772 -0.1961499697321569, 0.0000019152003648735755 -0.2786445563969257, 0.0000019003838131104091 -0.37811725071057417, 0.0000018874998550553546 -0.454121484123263, 0.0000018797694802223652 -0.5290052872535868, 0.0000018733275011948922 -0.6651076743801212, 0.0000018636645326536826 -0.8621294473729189, 0.000002532986153610129 -1.0197394947658354, 0.0000025258999766798436 -1.2435509780552996, 0.000002517525403944042 -1.7398272689135508, 0.000001823080064780407 -2.1874042970353043, 0.00000249369008154237 -2.9233188898260445, 0.0000024814503213901277 -3.7320371169840842, 0.0000017908701696429333 -4.912067090735175, 0.0000017792746073935035 -6.270830076016752, 0.0000017689674409494816 -8.128577892877692, 0.0000017580160766027557 -10.377083574864871, 0.0000017477089101587337 -13.868275049482747, 0.0000017354691500064915 -18.819079409319883, 0.0000017225851919515454 -25.930034799879476, 0.0000017090570359937869 -35.186739680618125, 0.0000016961730779388407 -47.747974844883, 0.0000016832891198837862 -65.7858703632878, 0.00000371057991983587 -92.0282264054625, 0.000007097772492491703 -113.93839291593221, 0.00001117039163367299 -147.662236140643, 0.000017962170122358856 -185.6103461800117, 0.00002543550999215352 -233.2963070659185, 0.00003494966881785808 -280.1206955348425, 0.00004310521426666435 -341.5020276966939, 0.000052620661488174316 -403.8330155958889, 0.0000607768511348834 -499.8741418664249, 0.00007165220012909717 -609.4341686323978, 0.00007980710138000058 -720.7280260248605, 0.00008524219908549654 -892.2269839498045, 0.00009271618315319401 -1121.6392214813618, 0.00009678815809647243 -1367.758754488643, 0.00009814032949434319 -1693.6422214882884, 0.0000974510377384014 -2263.6740781353897, 0.00009403743305173287 -3072.3521667520486, 0.00008790209222594845 -3983.205237157777, 0.00008244895697917553 -5406.171038138221, 0.00007631361615339117 -7008.7803886988195, 0.00007154075389192153 -8418.818380186558, 0.00006677111261996563 -11779.82325702495, 0.00006267530235428553 -15038.951357933107, 0.00006130444921723506 -18622.53255082284, 0.00005993488447599005 -25271.073403302787, 0.000059241727532631766 -33262.205588573604, 0.000058549858985079 -45831.57553825898, 0.00005785605784381796 -60324.27921334647, 0.0000571641892962652 -80617.6723282089, 0.000057832222521416235 -98313.606263342, 0.00005714357496337721 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_4p2m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_4p2m.csv deleted file mode 100644 index 1692f9097fc..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_4p2m.csv +++ /dev/null @@ -1,73 +0,0 @@ -0.009998754088411478, 0.000004081637911819468 -0.013362663107590534, 0.000004069398151667226 -0.017058653323225467, 0.000004739363970526485 -0.022797744140105793, 0.000004727124210374242 -0.029551628218007345, 0.0000047161728460275165 -0.038306365970614646, 0.000004705221481680791 -0.04890155591770684, 0.0000053751873005400495 -0.06339007058199968, 0.000004683962950890043 -0.08343332872231671, 0.000004672367388640613 -0.10489686247826191, 0.000005342977405402576 -0.14453607003290095, 0.000004649176264141645 -0.19023681993992214, 0.000004637580701892107 -0.2542387294495813, 0.000004625340941739864 -0.355700314756735, 0.000003930895602576121 -0.4976325182104469, 0.0000045969962340189395 -0.645043899575232, 0.000005266317854975386 -0.8234913350291972, 0.000004575737703228192 -1.1346312476178406, 0.000005242482532573714 -1.6618188048768607, 0.0000045461045997017505 -2.4713393026023933, 0.0000045293554542302555 -3.3535771692660803, 0.000004516471496175201 -4.413850519257689, 0.000005185148919229052 -5.721582692071965, 0.000004493924569579045 -8.127902712897296, 0.0000044791080178157704 -11.029235405122794, 0.000005146497045064105 -14.296968627598796, 0.000004455272695414098 -18.53170873696416, 0.000005804867301673826 -25.532487870755826, 0.000007832158101626018 -35.17431988860397, 0.000013260813828094398 -45.58721884623693, 0.00001869204634617381 -60.905280583255255, 0.00002888390136557067 -72.02171913747563, 0.00003704009101227954 -85.16536268344663, 0.0000458765536442919 -102.25446821781289, 0.00005539264506370463 -117.27276420768364, 0.00006627121504743219 -132.46228782962527, 0.00007647015624375928 -149.6160975249088, 0.00008734937042538955 -174.21514692485636, 0.00010094838815242742 -199.7816911518674, 0.00011522832306267133 -225.64396296383725, 0.0001274680832149082 -251.00420271602334, 0.00013834794159444123 -283.4855907790619, 0.00015194824771728462 -315.30739108928015, 0.00016690974400863728 -350.7522329441862, 0.00017710932940286707 -378.43483589839576, 0.0001886707491635089 -414.5735096208301, 0.0002009117977115513 -461.1772141359031, 0.00021111138310578103 -505.25934791346043, 0.0002206313397126103 -562.0457029828076, 0.00023151119809214338 -625.2143053779638, 0.0002423910564716764 -695.4824591216696, 0.00025327091485120954 -797.6292074144066, 0.0002641494848349371 -887.4041481726439, 0.00027026743231734726 -1049.351459962489, 0.0002791038949493595 -1339.5090086031735, 0.00028181467972412857 -1633.6049792431556, 0.00027976548619548303 -2023.2073126887071, 0.00027295373757181183 -2623.404183329673, 0.0002627386914279161 -3103.9680973001127, 0.0002525275104714368 -3616.858145407688, 0.00024299724669816353 -4279.40626425186, 0.0002327860657416843 -4986.418512926343, 0.0002239360749537143 -5990.713958814154, 0.00021304397681402902 -7307.371084458944, 0.00020487232641765401 -8778.391478736308, 0.00019670132021918182 -11557.383994177571, 0.00018716590286268657 -14532.676355490463, 0.00018307460198232567 -19721.89247021044, 0.00018102089906836088 -26764.03389796784, 0.00017896719615439615 -40413.7429276549, 0.0001789498028110219 -58292.29712146067, 0.0001796146150466592 -82808.22727093354, 0.00017959979849489598 -100982.90057642768, 0.00017959142392216023 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_7m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_7m.csv deleted file mode 100644 index cf467af410f..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationDisp_7m.csv +++ /dev/null @@ -1,70 +0,0 @@ -0.009997508332052525, 0.000008163275823639153 -0.014642076327913747, 0.000008827443861373642 -0.021444828476583546, 0.00000881133891380496 -0.03189056877732696, 0.000009474862753636746 -0.04599952839405873, 0.000009459402003970767 -0.06737104590376612, 0.000009443297056401976 -0.09717115108065555, 0.00001146865526264584 -0.13388534850727565, 0.000012135400091991362 -0.20842514486611707, 0.0000134772643234181 -0.3195496706161579, 0.000014819772752747758 -0.5050870011596059, 0.000017521538756878246 -0.7510360538100744, 0.000021586427523226435 -1.2427248786483838, 0.00002428626093364881 -1.738454866681777, 0.000027673453506304698 -2.469236035679274, 0.000032420547851664226 -3.5611526400792037, 0.00003716699799912116 -4.981302790847094, 0.00004327528251299017 -6.967649163227823, 0.0000500638400121623 -8.893359505790702, 0.000056175989713447803 -12.250740714758136, 0.0000643257373811292 -16.619272688510673, 0.00007383667521732007 -20.88943607464521, 0.0000826705610577213 -26.660064510881778, 0.0000921840756855231 -36.16085673700339, 0.00010713719740414 -43.41772108499518, 0.00011597301583824956 -52.12766613806268, 0.0001268496532282689 -63.54484090256706, 0.00013908619239099193 -75.13525041748713, 0.00015064374696421726 -86.1704985558791, 0.00016152231694794482 -100.34025795932699, 0.00017444106168967947 -118.63712596192376, 0.00018735916223351135 -138.13991587024188, 0.0002016384529458525 -153.65914921413193, 0.00021387885729599204 -170.92188009864876, 0.00022611926164613173 -187.25574825276573, 0.00023631949123826423 -214.7360360117358, 0.00025059942614850814 -238.85544727734836, 0.000263520103483951 -265.678462720635, 0.00027712105380469715 -291.04341093396084, 0.0002900423753380428 -328.6992769184907, 0.0003043229544461894 -360.08844830161087, 0.0003165640029942318 -388.50789829073506, 0.0003281254227548736 -419.17901626332485, 0.00033900656953021215 -473.44285592736765, 0.00035124632968244903 -526.631466256596, 0.0003634867340325886 -576.9102119716844, 0.0003764080555659342 -651.5792437063674, 0.00038932808870347434 -735.9126290692158, 0.00040224812184101445 -856.9784869877932, 0.00041312604762683926 -1013.3732508788464, 0.0004219625102588515 -1254.585997722622, 0.0004273956753706392 -1553.6016044563441, 0.0004246655646587876 -2014.4452252560188, 0.00041513079150019513 -2420.1166613447454, 0.0004049189663458131 -2863.501154352861, 0.00039402751240403053 -3601.867410039547, 0.00037905184375881736 -4197.377438863477, 0.00036680048804433105 -4743.890522157134, 0.00035727151266686327 -6059.026227402935, 0.0003416149268384441 -7505.148700187856, 0.00033004126731765005 -8879.782426667729, 0.000320510359346474 -11167.617583977653, 0.0003109768745836871 -14478.755753288755, 0.0003048434663516108 -18770.45904248312, 0.0003007508770754445 -23243.20866482086, 0.00029938131233419945 -30129.06263159623, 0.0002993703609698527 -42152.968135838084, 0.000298675915630689 -60803.45995643873, 0.0002979801818957197 -85067.13655571481, 0.00029796600954185926 -96115.22470876326, 0.00029796085595863725 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_0m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_0m.csv deleted file mode 100644 index bcd2945855e..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_0m.csv +++ /dev/null @@ -1,83 +0,0 @@ -0.010150923999646412, 0.987603305785124 -0.012900151049529792, 0.9876237877251735 -0.018206356747197474, 0.9876532305139943 -0.023840942487676983, 0.9890536831648694 -0.03169052018799551, 0.9877005950003585 -0.0408812497263315, 0.9877223570616609 -0.05434131097931401, 0.9877466793654696 -0.0690588482277122, 0.9877671613055189 -0.09043149333627863, 0.9864127930197549 -0.11665799236725406, 0.9878119655493769 -0.15740771780960836, 0.986460157506119 -0.19706507612585966, 0.9864793593249154 -0.2619483175442174, 0.9865036816287239 -0.3328931300843456, 0.9865241635687734 -0.46283640251960934, 0.9865523262363411 -0.5708284053867687, 0.9879476584022039 -0.7254288036348263, 0.9865907298739337 -0.9642746363079265, 0.9879924626460619 -1.2817599322880595, 0.9866393744815509 -1.6784448010661022, 0.9880398271324261 -2.2989216146542857, 0.9880667096787409 -3.3432093195047434, 0.988098712710068 -4.511023535244492, 0.9895017256034493 -5.819288499263777, 0.9895234876647517 -8.590429339305476, 0.9881793603490123 -11.41880924135124, 0.988203682652821 -15.64004211629874, 0.9882305651991358 -20.480393182283496, 0.9882536073816913 -28.051454854794255, 0.9882804899280061 -35.648774909408836, 0.9869235613997359 -45.98745810620617, 0.9855679129927187 -55.04353529985595, 0.9828284535111167 -65.88297990528858, 0.9745793521562364 -80.0471225377729, 0.9621992995176504 -92.9825574226782, 0.9470605855786661 -100.21422402961326, 0.9305380605650968 -112.97282209439288, 0.909887144510328 -131.22897594022814, 0.8795969154198287 -141.43523644191845, 0.8589421590013007 -154.73589743096068, 0.8314016303624281 -166.77039565689014, 0.8052372320706218 -179.7408702784364, 0.7832050651837743 -188.00247591190245, 0.7584155171177815 -202.6242637477963, 0.7405155816358926 -211.93768120406273, 0.7143486231015803 -238.92015515371875, 0.6813010128319356 -246.18633359727394, 0.6551327741763702 -269.33785542229856, 0.6220826036642191 -299.11378691303156, 0.5890337132733214 -317.58408897347647, 0.5546035720503447 -337.19493377427756, 0.5229282517640073 -352.6937547165059, 0.4940064722930557 -397.59634115030553, 0.4582040410867718 -428.5191750993205, 0.41964294858008966 -448.21562157565273, 0.3948534005140967 -468.8173951093013, 0.37144126291642354 -512.9052857579319, 0.3535426075557877 -544.577231254386, 0.32186728726945024 -569.6081678166662, 0.29707773920345737 -604.7815211557064, 0.27366688172703724 -632.5796864503114, 0.24749992319272496 -681.7782191808445, 0.2213355249009188 -702.512854131631, 0.1965446967136728 -791.9520453043104, 0.16762931784898671 -853.5456744651754, 0.13595527768390236 -933.8136610032172, 0.10979215951334931 -1037.048950938775, 0.07674326912245133 -1186.7232263265391, 0.050583991315657495 -2328.6293754671296, 0.012074103659098778 -3649.7881600453043, 0.007980275891732669 -4501.380496556548, 0.006620787120956262 -5806.848833698688, 0.005265138713939299 -7490.922707649899, 0.003909490306921892 -9809.247392683612, 0.0025551220211579295 -13435.467666382601, 0.0025820045674727066 -18679.939472771566, 0.0026101672350404836 -25585.42085334695, 0.0026370497813554827 -33005.57936739651, 0.001281401374338298 -45889.15251453425, 0.0026869745102255926 -58317.54813581721, 0.0013300459819556298 -77518.47216056062, 0.0013543682857641848 -98513.19939296534, 0.0013748502258135176 -1399.2998277410964, 0.028559505156328324 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_4p2m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_4p2m.csv deleted file mode 100644 index 751dcd6ea76..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_4p2m.csv +++ /dev/null @@ -1,78 +0,0 @@ -0.010150923999646412, 0.987603305785124 -0.013493093324807878, 0.9890050385572522 -0.01740629420154237, 0.9890268006185546 -0.02313727815405322, 0.9876737124540437 -0.030297907783064737, 0.9890741651049189 -0.03967464149927611, 0.9877197968191548 -0.052737429656322464, 0.9863667086546438 -0.0690588482277122, 0.9877671613055189 -0.09179632160310944, 0.9877914836093276 -0.11665799236725406, 0.9891893760176964 -0.1550673789056952, 0.9892136983215051 -0.2061229711975749, 0.9892380206253137 -0.2823212020209928, 0.9892649031716285 -0.3809386326441193, 0.9906679160650097 -0.5457438465435344, 0.9893212285067642 -0.7363772652851527, 0.989346830931826 -0.9358141142568936, 0.9893673128718753 -1.2254341045244705, 0.9893903550544308 -1.6046870010334717, 0.9880359867686668 -2.133026654795654, 0.9880603090724754 -2.8781127717460526, 0.9894633219658566 -3.7128089493104866, 0.9894850840271592 -4.7183681274551645, 0.9881281554988889 -6.0867661551458605, 0.9881499175601913 -7.852020280454728, 0.9867942691531743 -10.917019579814932, 0.9813127899474638 -13.873732243014787, 0.9716913986092763 -17.63122663137432, 0.9524281339928519 -19.875915308625075, 0.938664270279681 -24.149026509533922, 0.9180197548311777 -28.474818631052823, 0.8959952686718486 -30.68943178797696, 0.8794727436582794 -35.648774909408836, 0.8491825145677799 -40.18733612639283, 0.8257767775763721 -45.3037162999231, 0.7996162196483252 -51.07147943642877, 0.7762104826569174 -60.21986491247263, 0.7349022499411145 -69.95125957922657, 0.7004797894456564 -80.0471225377729, 0.6701882802339039 -90.23818146768905, 0.6426503118375373 -100.21422402961326, 0.6164884737882372 -114.67785311057571, 0.5875743750448044 -123.5968592541141, 0.5655422081579569 -137.26078186851976, 0.5393803701086568 -154.73589743096068, 0.5104649912439708 -169.28736117040538, 0.48430187307341754 -199.61164496439375, 0.45401292410417127 -221.679180387119, 0.41683180230831474 -246.18633359727394, 0.3906699642590147 -294.66656131348304, 0.3521178326011042 -337.19493377427756, 0.3121844501111146 -403.597024155423, 0.273632318453204 -448.21562157565273, 0.23507378618902797 -505.27940685576704, 0.21304545966593957 -604.7815211557064, 0.1689836861347508 -692.0678887519011, 0.13180512458140048 -803.9045023248603, 0.10151489549090109 -919.9297138228432, 0.07811043862074629 -1006.4404983000102, 0.0629666041967496 -1204.6337279055886, 0.031301524880436804 -1337.8088112866583, 0.013404149641053964 -1974.8723687289714, 0.009305201388675632 -2959.30406758484, 0.0065849437258698185 -3993.0165960405707, 0.005233135682612078 -4999.018665795859, 0.003874927033088671 -5983.45009110729, 0.0038902884881255595 -7379.547623360033, 0.0025307997173493746 -8701.438961584417, 0.002544881051133263 -11394.399699165459, 0.0011905127653693004 -16323.950862412092, 0.0025986461437628172 -20745.050584316636, 0.00261912808381215 -29719.9672924797, 0.002649850993886149 -36654.42356166941, 0.002667772691429482 -47284.76005789792, 0.0013121242844122971 -54925.87088594401, 0.0013249254969431856 -69801.72750582406, 0.0013454074369925184 -90045.28283432488, 0.0013671694982950733 -100000.0000000002, 0.0013761303470665176 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_5p6m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_5p6m.csv deleted file mode 100644 index ec3f4cd6998..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationPressure_5p6m.csv +++ /dev/null @@ -1,76 +0,0 @@ -0.010150923999646412, 0.987603305785124 -0.013292477931346833, 0.9890037584359991 -0.01793566452456119, 0.9890293608610609 -0.02313727815405322, 0.9890511229223633 -0.03169052018799551, 0.9890780054686781 -0.04340567041200643, 0.9877274775466733 -0.058567676169505485, 0.987753079971735 -0.07223306769529865, 0.9877710016692782 -0.08776241192181838, 0.9877876432455683 -0.11665799236725406, 0.9864345550810574 -0.18284445081830478, 0.9878503691869694 -0.2504372260790477, 0.9878772517332842 -0.3042785762390852, 0.9878938933095743 -0.37527483474153583, 0.9879118150071176 -0.46283640251960934, 0.9893071471729803 -0.633934824145139, 0.9893340297192951 -0.8056266589439807, 0.9879771011910249 -1.1034456053673303, 0.9880039837373397 -1.5341704235076514, 0.986654735936588 -1.8640008096307419, 0.9852939670445585 -2.4408802220263235, 0.9811847778221554 -3.1019563892948305, 0.9743182074206069 -4.1232686848701805, 0.96056842504122 -5.162086513789205, 0.9385490593669032 -6.0867661551458605, 0.9165245732075743 -7.070373869743989, 0.8931213964586726 -8.090820479054509, 0.8697169395885179 -9.120889089634584, 0.8490660235337492 -10.129224766183718, 0.8270364168894079 -11.766084696229441, 0.7953687773305889 -13.464249666614887, 0.7650772681188363 -14.952754039751238, 0.7402928405378558 -16.358920687741605, 0.7141297223673027 -17.897324155562217, 0.6934762460700279 -19.875915308625075, 0.6686918184890474 -21.745059222245825, 0.6439061107868138 -23.789978636797112, 0.6246300449578586 -26.027203593921122, 0.599844337255625 -29.783634604013276, 0.5764398803854702 -34.596602938686715, 0.5365077780167338 -40.18733612639283, 0.5103497803311932 -48.10122467677227, 0.46628800680000426 -58.44247833740541, 0.4291145657316662 -68.91122382717396, 0.4002030272307393 -78.8569814339672, 0.3767985703605846 -87.57480564424162, 0.3575237846528825 -97.25640576349953, 0.33962640941349975 -111.29314148970882, 0.31622195254334495 -125.46223248835038, 0.2955710364885763 -143.56983359939346, 0.2762988110233803 -174.43582666387255, 0.2528994746382379 -199.61164496439375, 0.23224983870472227 -238.92015515371875, 0.21298145360328535 -281.7176950225821, 0.1909569674439563 -312.8622468107721, 0.1730595922045738 -374.47262439575746, 0.14690415476153917 -422.1479494026821, 0.1317628805800487 -497.76690956741334, 0.11111580488903916 -586.931422054467, 0.08909131872971021 -671.6415364794294, 0.0711965037328337 -768.5776166923116, 0.05743392014091586 -906.2521932534311, 0.03816425491822595 -1068.5882856063554, 0.025781642037133734 -1298.323445043419, 0.007891947525269671 -1945.5099543625452, 0.0010394584575053045 -2959.30406758484, 0.0010753018525915259 -4434.4539440087965, -0.001644955810214288 -5983.45009110729, -0.001619353385152733 -8966.072174173372, -0.0029622005796388073 -13435.467666382601, 0.001204594099153189 -19538.54382325673, -0.00014081333783910743 -26761.428606782032, -0.00011393079152433039 -41321.01717198699, 0.0013006031931346307 -56596.3083539263, 0.0013274857394494077 -78688.41194705585, 0.0027330588753367024 -100000.0000000002, 0.0013761303470665176 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_0m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_0m.csv deleted file mode 100644 index 0de5c5011e3..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_0m.csv +++ /dev/null @@ -1,79 +0,0 @@ -0.009847522613947512, -0.00006538323730609363 -0.012984704504125002, 0.06969853096941847 -0.01738715297707653, 0.002353796543054898 -0.022926761479327634, 0.0035306948145858996 -0.03069940537329468, 0.004772976323422995 -0.041107135481112604, 0.006015257832252985 -0.05420287194642072, 0.07577917203898465 -0.07600262436389589, 0.07721760325974003 -0.10334485319892368, 0.07852526800588322 -0.13627088934753467, 0.07970216627741422 -0.18246952481410988, 0.08094444778625132 -0.2519553695461515, 0.0823174957697006 -0.3271634977522269, 0.08342901080392551 -0.4313989465438657, 0.0846059090754494 -0.6049024877789514, 0.08604434029621189 -0.888199860485518, 0.08767892122888554 -1.1893178892332754, 0.08892120273772264 -1.693504709248258, 0.02183800126058344 -2.2676376520806625, 0.023080282769420535 -3.1796546285384446, 0.024518713990183016 -4.458474019210867, 0.025957145210938393 -6.0624267499036755, 0.02726480995708158 -8.500662906977, 0.028703241177836958 -12.10383360127639, 0.0987940715711062 -17.234641266451327, 0.10029788602916767 -24.166718992123197, 0.03314930131472238 -36.034295475456396, 0.034849265484709235 -47.51595108902488, -0.03256085217896043 -67.65793371178145, -0.03105703772089896 -91.99618979378293, 0.038837642960444896 -119.4568517711176, 0.039949157994669804 -155.1177308394172, -0.02752634290631306 -210.9219156090737, -0.026218678160169873 -300.33158671082236, -0.024714863702101297 -414.6827600431619, 0.11383221615174932 -581.4264838157069, 0.3210316951781138 -743.3282497011799, 0.9393609703918457 -880.0194011992243, 1.625950345354255 -979.7215232722228, 2.3808652033026334 -1074.0431528439112, 3.272888709884114 -1214.0881612533528, 4.43939104668101 -1372.364780608189, 5.674480399413113 -1504.424295895417, 6.703677937864995 -1623.9766943305667, 7.869984124949973 -1726.3726233559503, 8.89905089692725 -1835.1475197500574, 10.065291700774921 -2011.2731934570777, 11.848946414514025 -2171.1949815654075, 12.878078569728608 -2343.5855923579347, 14.250145804619194 -2491.092184401604, 15.622147656272467 -2607.61861700225, 16.85691009281804 -2688.202299554259, 17.81725908238549 -2901.8873934524395, 18.914978253535267 -3037.245573084486, 20.561262785692044 -3328.811058957642, 22.276330483495958 -3704.388051751737, 24.40298566014839 -4059.3980338875435, 26.598162469498714 -4314.171321748316, 28.518860448633614 -4727.719521263503, 30.64545024204874 -5102.4501979037805, 32.429039572550536 -5591.797012988701, 34.418455334095256 -6318.6496999461115, 36.7509369417906 -6716.207497680589, 38.19152580937909 -7249.31455296477, 39.632180060204874 -8067.737044641319, 41.5530741890517 -8978.178124241254, 43.611142349768926 -10147.9915353643, 45.05199275030664 -11295.807152428426, 46.355603735736636 -12768.403810586182, 47.59069308846874 -15115.120130475732, 48.55163052717196 -17621.425778644534, 49.306741534832256 -19924.130714683673, 49.65019968040673 -23956.254936660134, 49.92533234299522 -30632.831165884967, 49.92637847479214 -39776.66444790172, 49.927489989826356 -51648.82340157383, 49.99718852079578 -65036.27075400391, 49.99816926935539 -83163.54530608773, 49.930628385217105 -98477.30173013602, 49.93134760082748 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_4p2m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_4p2m.csv deleted file mode 100644 index 6587dbad78e..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_4p2m.csv +++ /dev/null @@ -1,79 +0,0 @@ -0.010154838320286299, 0.00006538323730609363 -0.013808079447942344, 0.0013730479834492826 -0.018489302502721943, 0.002615329492286378 -0.024008824998343215, -0.06486017140869649 -0.029770349769560688, 0.004642209848803702 -0.03806730645718513, 0.005688341645722517 -0.05019568508168139, 0.006865239917246413 -0.06418515541884809, 0.007911371714165227 -0.07958968428283685, 0.008826737036464749 -0.10334703142811452, 0.009938252070682552 -0.13013472669692205, 0.010919000630295272 -0.16898322037911714, -0.05655650027068759 -0.2160786180238699, -0.05551036847376878 -0.27629361652116846, 0.014122779258343598 -0.3532888317617504, 0.08375592699045598 -0.45874477739519076, 0.08486744202468088 -0.5956790927636327, 0.08597895705889869 -0.7386429247941366, 0.08689432238119821 -0.9591262266046674, 0.08800583741542312 -1.2647072384377274, 0.08918273568695412 -1.6171794552933025, 0.09022886748386583 -2.1990121708282984, 0.022949516294808348 -3.0834287924135517, 0.024387947515563724 -4.257622973771922, 0.025760995499013006 -5.614116728441364, 0.026937893770544008 -7.069303646215659, 0.027918642330149623 -9.61251211588196, 0.029226307076292812 -11.558564609652718, 0.09859792185918081 -15.716801528530866, 0.099905586605324 -22.72565899768972, 0.10147478430069867 -30.430125012396537, 0.10271706580952866 -40.12527058714834, 0.10389396408105966 -54.56163698475552, 0.036614612892002185 -74.18882658893563, 0.10650929357334604 -99.33408413730689, 0.31351262288779225 -141.4178942814683, 0.8637125648274733 -183.5766496504274, 1.825042302954536 -231.06742154147815, 3.12917635428299 -273.51244862121325, 4.364461856727019 -313.956479536012, 5.599616592696428 -360.3885274683878, 6.766184312730637 -420.03105544083894, 8.41292652754857 -482.10000949068575, 9.92242932725879 -520.4440584795021, 10.882974466538158 -579.3224804765924, 12.117998436032963 -654.7640094544381, 13.764609884376277 -763.1724859787392, 15.205591051388609 -889.417536437656, 17.058094314012145 -989.9966098796707, 18.430292315377343 -1085.215965767925, 19.596663885699634 -1189.5936618416808, 20.763035456021917 -1344.6203326069424, 22.13529884062443 -1473.9479289919866, 23.301670410946716 -1717.8790585492452, 24.94841262576465 -2001.715442798373, 27.349612015869813 -2369.063183064449, 29.065006629860253 -2636.744496793452, 30.711552694966265 -2980.2996180407863, 32.152403095503985 -3266.880439608481, 33.38736168176147 -3580.7165156964043, 34.89666833175977 -4047.181345484366, 36.40610574823268 -4504.6637290568715, 37.91547778146829 -5014.070037651899, 39.28767578283349 -5496.214954253563, 40.52263436909098 -6024.595169045325, 41.82617997128367 -6809.570932258092, 43.267030371821384 -7579.945354707946, 44.50205434131618 -8970.990124537786, 46.21744895530662 -10949.807320316113, 47.59003925609567 -12764.636679653275, 48.55091131156158 -14881.82547508411, 49.16884828735147 -18450.373518659027, 49.71845978015539 -22528.21972490492, 49.925070810045995 -29705.1641920149, 49.994834724252726 -38571.27792798286, 50.06453325522215 -50860.22364751764, 50.065710153493676 -62105.12456086096, 50.06656013557867 -81892.03714535893, 50.0677370338502 -99997.89231566602, 50.068587015935194 diff --git a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_5p6m.csv b/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_5p6m.csv deleted file mode 100644 index b0a5c8a2df0..00000000000 --- a/src/docs/sphinx/advancedExamples/validationStudies/thermoPoromechanics/thermalConsolidation/thermoConsolidationTemp_5p6m.csv +++ /dev/null @@ -1,81 +0,0 @@ -0.010312074131115496, 0.00013076647461218727 -0.016101391474977494, 0.0020268803565173243 -0.021560094495253666, 0.0032691618653544197 -0.0284292169405849, 0.004446060136885421 -0.03635240170142669, 0.00549219193379713 -0.046483767464426694, 0.006538323730715945 -0.06129236848963046, 0.07630223793744051 -0.08463298936573758, 0.07767528592088979 -0.10822009037553766, 0.0787214177178086 -0.14052355128648145, 0.07983293275202641 -0.20633573372439612, 0.08146751368470717 -0.28491002042159363, 0.08284056166815645 -0.40569175463007623, 0.01575736019102436 -0.5516408831558136, 0.01706502493716755 -0.761710114134855, 0.018438072920616833 -1.004372484084059, 0.0882019871273414 -1.4747558792731383, 0.08983656806002216 -2.0999039950600644, 0.09134038251809073 -3.0833638034679507, 0.09297496345076439 -4.390399835067726, 0.09447877790883297 -5.969862622996343, 0.09578644265497616 -7.517259374810218, 0.09676719121458177 -9.912075908769697, 0.16653110542132055 -14.114096308388454, 0.09944790394417424 -20.40823250091457, 0.10101710163954891 -27.748415975181043, 0.3080858141913012 -37.723855920645136, 0.9266766223542646 -47.487917451255115, 1.887875594006708 -55.36093450283487, 2.7115736176022125 -64.53513749551234, 3.741032689003326 -75.22331221585263, 5.044839824145235 -83.7404698922964, 6.00551572989923 -90.40080720133093, 6.966060869178605 -100.63643637440106, 7.926736774932586 -110.3227750905378, 8.887347297449267 -120.9337847435934, 10.053718867771558 -132.56817130753677, 11.151503422158633 -145.32490238426186, 12.180700960610523 -161.76567212667092, 13.41572493010532 -174.63178026735855, 14.376270069384681 -188.52120031663478, 15.336815208664056 -209.84881585339235, 16.571839178158854 -230.0226924789056, 17.875384780351546 -260.0154630833952, 19.04188711714845 -289.4191045522492, 20.414085118513647 -327.135916329084, 21.786348503116162 -353.154823930799, 22.74689364239553 -393.1241980306451, 23.844743580019923 -430.94455482528394, 24.942528134407006 -479.7079079876418, 26.1089650879666 -542.2343972238937, 27.412641456633907 -612.9753616051692, 28.3733827456252 -692.8723087084063, 29.67705911429251 -783.2162655995311, 30.84356145108941 -899.0297249815874, 32.078716187058816 -1048.013503299308, 33.10817525845994 -1166.6748064473145, 34.06885116421392 -1298.8262824664105, 34.8923530380975 -1514.0309103321279, 35.99039912543381 -1711.5543563700394, 36.9511404144251 -2026.0797531807355, 37.98066486906352 -2325.8703601749457, 38.941471541292124 -2670.0760150629, 39.83369119758552 -2972.457956261255, 40.7257800874043 -3464.896292992599, 41.89241319067582 -4164.959538113848, 43.059177060421945 -4708.526818887187, 43.88274431754284 -5488.574560273634, 45.04937742081435 -6398.2550652825785, 46.01024947628026 -7809.908666132243, 47.2456657451989 -9245.493608546547, 48.13801616796691 -11287.001535141568, 48.89332332533913 -13780.168450737357, 49.442869434905745 -17618.454765575774, 49.855437662313875 -22528.21972490492, 49.925070810045995 -30165.749763528773, 49.92631309155483 -40391.70829122072, 49.99614238899887 -49322.08108319346, 49.99699237108386 -62107.74260359248, 49.92938610370827 -78206.15645939638, 49.93036685226787 -98477.30173013602, 49.93134760082748