feat: Temperature dependent Thermal Expansion Coefficient

inputFiles/thermoPoromechanics/ThermoElasticOedometric_linearTEC_base.xml

@@ -0,0 +1,245 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Solvers>
+    <SinglePhasePoromechanics
+      name="thermoElasticSolver"
+      solidSolverName="solidMechSolver"
+      flowSolverName="flowSolver"
+      logLevel="1"
+      isThermal="1"
+      targetRegions="{ rock }">
+      <NonlinearSolverParameters
+        couplingType="FullyImplicit"
+        maxAllowedResidualNorm="1e19"
+        newtonMaxIter="100"/>
+      <LinearSolverParameters
+        directParallel="0"/>
+    </SinglePhasePoromechanics>
+
+    <SolidMechanicsLagrangianSSLE
+      name="solidMechSolver"
+      timeIntegrationOption="QuasiStatic"
+      logLevel="1"
+      discretization="FE1"
+      targetRegions="{ rock }">
+    </SolidMechanicsLagrangianSSLE>
+
+    <SinglePhaseFVM
+      name="flowSolver"
+      logLevel="1"
+      discretization="tpfaFlow"
+      temperature="0"
+      isThermal="1"
+      targetRegions="{ rock }">
+      <NonlinearSolverParameters
+	normType="L2"/>
+    </SinglePhaseFVM>
+  </Solvers>
+
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace
+        name="FE1"
+        order="1"/>
+    </FiniteElements>
+
+    <FiniteVolume>
+      <TwoPointFluxApproximation
+        name="tpfaFlow"/>
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <Constitutive>
+    <PorousElasticIsotropic
+      name="rock_linearTEC"
+      solidModelName="rockSolid_linearTEC"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm"
+      solidInternalEnergyModelName="rockInternalEnergy"/>
+
+    <!-- SPHINX_PoroElasticProperties -->
+     <BiotPorosity
+      name="rockPorosity"
+      defaultReferencePorosity="1e-6"
+      defaultGrainBulkModulus="5.535714286e9"/>
+    <!-- SPHINX_PoroElasticPropertiesEnd -->
+
+    <!-- SPHINX_HeatCapacityProperties -->
+    <SolidInternalEnergy
+      name="rockInternalEnergy"
+      referenceVolumetricHeatCapacity="4.56e6"
+      referenceTemperature="0"
+      referenceInternalEnergy="0"/>
+    <!-- SPHINX_HeatCapacityPropertiesEnd -->
+
+    <!-- SPHINX_PermeabilityProperties -->
+    <ConstantPermeability
+      name="rockPerm"
+      permeabilityComponents="{ 1.0e-100, 1.0e-100, 1.0e-100 }"/>
+    <!-- SPHINX_PermeabilityPropertiesEnd -->
+
+    <!-- SPHINX_ThermalProperties -->
+    <SinglePhaseThermalConductivity
+      name="rockThermalCond"
+      defaultThermalConductivityComponents="{ 1.66, 1.66, 1.66 }"/>
+    <!-- SPHINX_ThermalPropertiesEnd -->
+
+    <!-- SPHINX_FluidProperties -->
+    <ThermalCompressibleSinglePhaseFluid
+      name="fluid"
+      defaultDensity="1000"
+      defaultViscosity="1e-3"
+      referencePressure="0.0"
+      referenceTemperature="20.0"
+      compressibility="5e-10"
+      thermalExpansionCoeff="1e-10"
+      viscosibility="0.0"
+      specificHeatCapacity="1"
+      referenceInternalEnergy="1"/>
+    <!-- SPHINX_FluidPropertiesEnd -->
+
+    <!-- SPHINX_ThermoElasticProperties -->
+    <ElasticIsotropic
+      name="rockSolid_linearTEC"
+      defaultDensity="2700"
+      defaultBulkModulus="9.12345679e9"
+      defaultShearModulus="6.00813e9"
+      defaultDrainedTEC="0.2e-5"
+      dDrainedTEC_dT="0.0164e-5"
+      referenceTemperature="0"/>
+    <!-- SPHINX_ThermoElasticPropertiesEnd -->
+  </Constitutive>
+
+  <FieldSpecifications>
+
+    <!-- SPHINX_FixedFarFieldBC -->
+    <FieldSpecification
+      name="constraint_x"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="0"
+      scale="0.0"
+      setNames="{ xneg, xpos }"/>
+
+    <FieldSpecification
+      name="constraint_y"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="1"
+      scale="0.0"
+      setNames="{ yneg, ypos }"/>
+    <!-- SPHINX_FixedFarFieldBCEnd -->
+
+    <!-- SPHINX_PlaneStrainAxisymmetryBC -->
+    <FieldSpecification
+      name="zConstraint"
+      objectPath="nodeManager"
+      fieldName="totalDisplacement"
+      component="2"
+      scale="0.0"
+      setNames="{ zneg }"/>
+    <!-- SPHINX_PlaneStrainAxisymmetryBCEnd -->    
+
+    <FieldSpecification
+      name="initialSigma_x_rock"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/rock/rock"
+      fieldName="rockSolid_linearTEC_stress"
+      component="0"
+      scale="0"/>
+    <FieldSpecification
+      name="initialSigma_y_rock"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/rock/rock"
+      fieldName="rockSolid_linearTEC_stress"
+      component="1"
+      scale="0"/>
+    <FieldSpecification
+      name="initialSigma_z_rock"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions/rock/rock"
+      fieldName="rockSolid_linearTEC_stress"
+      component="2"
+      scale="0"/>
+
+    <!-- SPHINX_ZeroPorePressureBC -->
+    <FieldSpecification
+      name="zeroPressure"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="pressure"
+      scale="0e6"/>    
+    <!-- SPHINX_ZeroPorePressureBCEnd -->
+
+    <!-- SPHINX_TemperatureBC -->
+    <FieldSpecification
+      name="initialTemperature"
+      initialCondition="1"
+      setNames="{ all }"
+      objectPath="ElementRegions"
+      fieldName="temperature"
+      scale="0"/>
+
+    <FieldSpecification
+      name="appliedTemperature"
+      setNames="{ all }"     
+      objectPath="ElementRegions"
+      fieldName="temperature"
+      scale="1"
+      functionName="temperatureTable"
+      />
+    <!-- SPHINX_TemperatureBCEnd -->
+  </FieldSpecifications>
+
+  <Functions>
+    <!-- SPHINX_TemperatureTable -->
+    <TableFunction
+      name="temperatureTable"
+      inputVarNames="{ time }"
+      coordinates="{ 0.0, 1.0e5 }"
+      values="{ 0, 90 }"/>
+    <!-- SPHINX_TemperatureTableEnd -->
+  </Functions>  
+
+  <Tasks>
+    <!-- SPHINX_TemperatureCollection -->    
+    <PackCollection
+      name="temperatureCollection_rock"
+      objectPath="ElementRegions/rock/rock"
+      fieldName="temperature"/>
+    <!-- SPHINX_TemperatureCollectionEnd -->
+
+    <!-- SPHINX_SolidStressCollection -->    
+    <PackCollection
+      name="stressCollection_rock"
+      objectPath="ElementRegions/rock/rock"
+      fieldName="rockSolid_linearTEC_stress"/>
+    <!-- SPHINX_SolidStressCollectionEnd --> 
+  </Tasks>
+
+  <Outputs>
+    <!-- SPHINX_TemperatureTimeHistory -->    
+    <TimeHistory
+      name="temperatureHistoryOutput_rock"
+      sources="{ /Tasks/temperatureCollection_rock }"
+      filename="temperatureHistory_rock"/>
+    <!-- SPHINX_TemperatureTimeHistoryEnd -->
+
+    <!-- SPHINX_SolidStressTimeHistory -->    
+    <TimeHistory
+      name="stressHistoryOutput_rock"
+      sources="{ /Tasks/stressCollection_rock }"
+      filename="stressHistory_rock"/>
+    <!-- SPHINX_SolidStressTimeHistoryEnd -->
+
+    <VTK
+      name="vtkOutput"/> 
+
+    <Restart
+      name="restartOutput"/>  
+  </Outputs>
+</Problem>

inputFiles/thermoPoromechanics/ThermoElasticOedometric_linearTEC_benchmark.xml

@@ -0,0 +1,68 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Included>
+    <File name="./ThermoElasticOedometric_linearTEC_base.xml" />
+  </Included>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ 0, 0.1 }"
+      yCoords="{ 0, 0.1 }"
+      zCoords="{ 0, 0.1 }"
+      nx="{ 2 }"
+      ny="{ 2 }"
+      nz="{ 2 }"
+      cellBlockNames="{ rock }"/>
+  </Mesh>
+
+  <ElementRegions>
+  <CellElementRegion
+      name="rock"
+      cellBlocks="{ rock }"
+      materialList="{ fluid, rock_linearTEC, rockThermalCond }"/>
+  </ElementRegions>
+
+  <Events
+    maxTime="1e5">
+
+
+    <PeriodicEvent
+      name="stressHistoryCollection_rock"
+      endTime="1e5"
+      forceDt="1e4"
+      target="/Tasks/stressCollection_rock"/>
+
+    <PeriodicEvent
+      name="stressTimeHistoryOutput_rock"
+      endTime="1e5"
+      forceDt="1e4"
+      target="/Outputs/stressHistoryOutput_rock"/>
+
+
+    <PeriodicEvent
+      name="temperatureHistoryCollection_rock"
+      endTime="1e5"
+      forceDt="1e4"
+      target="/Tasks/temperatureCollection_rock"/>
+
+    <PeriodicEvent
+      name="temperatureTimeHistoryOutput_rock"
+      endTime="1e5"
+      forceDt="1e4"
+      target="/Outputs/temperatureHistoryOutput_rock"/>
+
+    <PeriodicEvent
+      name="outputs"
+      timeFrequency="1e4"
+      target="/Outputs/vtkOutput"/>
+
+    <PeriodicEvent
+      name="solverApplications"
+      maxEventDt="1e3"
+      target="/Solvers/thermoElasticSolver"/>
+  </Events>
+
+</Problem>

inputFiles/thermoPoromechanics/ThermoElasticOedometric_linearTEC_smoke.xml

@@ -0,0 +1,47 @@
+<?xml version="1.0" ?>
+
+<Problem>
+  <Included>
+    <File name="./ThermoElasticOedometric_linearTEC_base.xml" />
+  </Included>
+
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ 0, 0.1 }"
+      yCoords="{ 0, 0.1 }"
+      zCoords="{ 0, 0.1 }"
+      nx="{ 2 }"
+      ny="{ 2 }"
+      nz="{ 2 }"
+      cellBlockNames="{ rock }"/>
+  </Mesh>
+
+  <ElementRegions>
+    <CellElementRegion
+      name="rock"
+      cellBlocks="{ rock }"
+      materialList="{ fluid, rock_linearTEC, rockThermalCond }"/>
+  </ElementRegions>
+
+  <Events
+    maxTime="1e4">
+
+    <PeriodicEvent
+      name="outputs"
+      timeFrequency="1e4"
+      target="/Outputs/vtkOutput"/>
+
+    <PeriodicEvent
+      name="solverApplications"
+      maxEventDt="1e3"
+      target="/Solvers/thermoElasticSolver"/>
+
+    <PeriodicEvent
+      name="restarts"
+      timeFrequency="5e3"
+      target="/Outputs/restartOutput"/>
+  </Events>
+
+</Problem>