Fixed solver by changing order of temp b.c. and temp gathering and mi…

…nor change to thermal_deformation_heat_transfer.json file
ORNL · Nov 6, 2024 · ab896df · ab896df
1 parent c061c7e
commit ab896df
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 4 deletions.
diff --git a/examples/thermomechanics/inputs/thermal_deformation_heat_transfer.json b/examples/thermomechanics/inputs/thermal_deformation_heat_transfer.json
@@ -10,7 +10,7 @@
     "horizon"                 : {"value": 0.00615, "unit": "m"},  
     "final_time"              : {"value": 8, "unit": "s"},
     "timestep"                : {"value": 4.2e-07, "unit": "s"},
-    "thermal_subcycle_steps"  : {"value": 1e5},
+    "thermal_subcycle_steps"  : {"value": 5e+4},
     "output_frequency"        : {"value": 10},
     "output_reference"        : {"value": true}
 }
diff --git a/src/CabanaPD_Solver.hpp b/src/CabanaPD_Solver.hpp
@@ -281,14 +281,15 @@ class SolverElastic
                     computeHeatTransfer( *heat_transfer, *particles, *neighbors,
                                          neigh_iter_tag{}, dt );
             }
-            if constexpr ( is_temperature_dependent<
-                               typename force_model_type::thermal_type>::value )
-                comm->gatherTemperature();
 
             // Add non-force boundary condition.
             if ( !boundary_condition.forceUpdate() )
                 boundary_condition.apply( exec_space(), *particles, step * dt );
 
+            if constexpr ( is_temperature_dependent<
+                               typename force_model_type::thermal_type>::value )
+                comm->gatherTemperature();
+
             // Compute internal forces.
             updateForce();