fixup: add horizon tolerance for all thermal problems

ORNL · Jan 27, 2025 · aa77cf6 · aa77cf6
1 parent 9df6df1
commit aa77cf6
Show file tree

Hide file tree

Showing 4 changed files with 4 additions and 0 deletions.
diff --git a/examples/thermomechanics/thermal_crack.cpp b/examples/thermomechanics/thermal_crack.cpp
@@ -43,6 +43,7 @@ void thermalCrackExample( const std::string filename )
     double K = E / ( 3 * ( 1 - 2 * nu ) );
     double G0 = inputs["fracture_energy"];
     double delta = inputs["horizon"];
+    delta += 1e-10;
     double alpha = inputs["thermal_expansion_coeff"];
 
     // Problem parameters

diff --git a/examples/thermomechanics/thermal_deformation.cpp b/examples/thermomechanics/thermal_deformation.cpp
@@ -41,6 +41,7 @@ void thermalDeformationExample( const std::string filename )
     double nu = 0.25;
     double K = E / ( 3 * ( 1 - 2 * nu ) );
     double delta = inputs["horizon"];
+    delta += 1e-10;
     double alpha = inputs["thermal_expansion_coeff"];
 
     // Problem parameters

diff --git a/examples/thermomechanics/thermal_deformation_heat_transfer.cpp b/examples/thermomechanics/thermal_deformation_heat_transfer.cpp
@@ -41,6 +41,7 @@ void thermalDeformationHeatTransferExample( const std::string filename )
     double nu = 0.25;
     double K = E / ( 3 * ( 1 - 2 * nu ) );
     double delta = inputs["horizon"];
+    delta += 1e-10;
     double alpha = inputs["thermal_expansion_coeff"];
     double kappa = inputs["thermal_conductivity"];
     double cp = inputs["specific_heat_capacity"];

diff --git a/examples/thermomechanics/thermal_deformation_heat_transfer_prenotched.cpp b/examples/thermomechanics/thermal_deformation_heat_transfer_prenotched.cpp
@@ -43,6 +43,7 @@ void thermalDeformationHeatTransferPrenotchedExample(
     double K = E / ( 3 * ( 1 - 2 * nu ) );
     double G0 = inputs["fracture_energy"];
     double delta = inputs["horizon"];
+    delta += 1e-10;
     double alpha = inputs["thermal_expansion_coeff"];
     double kappa = inputs["thermal_conductivity"];
     double cp = inputs["specific_heat_capacity"];