Changes to HIP example

Author: pabloseleson

examples/thermomechanics/hip_cylinder.cpp

@@ -160,7 +160,7 @@ void HIPCylinderExample( const std::string filename )
     temp = particles.sliceTemperature();
     // This is purposely delayed until after solver init so that ghosted
     // particles are correctly taken into account for lambda capture here.
-    auto temp_func = KOKKOS_LAMBDA( const int pid, const double t )
+    auto force_temp_func = KOKKOS_LAMBDA( const int pid, const double t )
     {
         // ---------------------
         // Isostatic pressure BC
@@ -197,10 +197,9 @@ void HIPCylinderExample( const std::string filename )
         // ---------------------
         //    Temperature BC
         // ---------------------
-        // temp( pid ) = temp0 + 5000.0 * ( x( pid, 1 ) - low_corner_y ) * t;
-        // temp( pid ) = Tmax;
+        temp( pid ) = Tmax;
     };
-    CabanaPD::BodyTerm body_term( temp_func, particles.size(), false );
+    CabanaPD::BodyTerm body_term( force_temp_func, particles.size(), false );
 
     // ====================================================
     //                   Simulation run

examples/thermomechanics/inputs/hip_cylinder.json

@@ -9,8 +9,8 @@
     "thermal_expansion_coeff"  : {"value": 14.56E-6, "unit": "oC^{-1}"},
     "thermal_conductivity"     : {"value": 14.12, "unit": "W/(m.K)"},
     "specific_heat_capacity"   : {"value": 492, "unit": "J/(kg.K)"},
-    "reference_temperature"    : {"value": 20.0, "unit": "oC"},
     "maximum_pressure"         : {"value": 110e6, "unit": "Pa"},
+    "reference_temperature"    : {"value": 1125, "unit": "oC"},
     "maximum_temperature"      : {"value": 1125, "unit": "oC"},
     "cylinder_outer_radius"    : {"value": 0.1187, "unit": "m"},
     "cylinder_inner_radius"    : {"value": 0.0679, "unit": "m"},