Fixed heat transer solver and add temporary example

examples/thermomechanics/inputs/thermal_deformation_cube_temp.json

@@ -1,15 +1,15 @@
 {
-    "num_cells"               : {"value": [51, 51, 51]},
-    "system_size"             : {"value": [0.5, 0.5, 0.5], "unit": "m"},
-    "density"                 : {"value": 3980,  "unit": "kg/m^3"},
-    "elastic_modulus"         : {"value": 370e+9, "unit": "Pa"},
-    "thermal_expansion_coeff" : {"value": 7.5E-6, "unit": "oC^{-1}"},
-    "thermal_conductivity"    : {"value": 31, "unit": "W/(m.K)"},
-    "specific_heat_capacity"  : {"value": 880, "unit": "J/(kg.K)"},
+    "num_cells"               : {"value": [101, 101, 101]},
+    "system_size"             : {"value": [0.1, 0.1, 0.1], "unit": "m"},
+    "density"                 : {"value": 8915,  "unit": "kg/m^3"},
+    "elastic_modulus"         : {"value": 115e+9, "unit": "Pa"},
+    "thermal_expansion_coeff" : {"value": 17e-6, "unit": "oC^{-1}"},
+    "thermal_conductivity"    : {"value": 387, "unit": "W/(m.K)"},
+    "specific_heat_capacity"  : {"value": 385, "unit": "J/(kg.K)"},
     "reference_temperature"   : {"value": 100.0, "unit": "oC"},
-    "horizon"                 : {"value": 0.03, "unit": "m"},  
-    "final_time"              : {"value": 0.01, "unit": "s"},
-    "timestep"                : {"value": 1.0E-7,  "unit": "s"},
+    "horizon"                 : {"value": 0.003, "unit": "m"},  
+    "final_time"              : {"value": 8, "unit": "s"},
+    "timestep"                : {"value": 0.004,  "unit": "s"},
     "output_frequency"        : {"value": 100},
     "output_reference"        : {"value": true}
 }

src/CabanaPD_HeatTransfer.hpp

@@ -59,8 +59,7 @@ class HeatTransfer : public Force<ExecutionSpace, BaseForceModel>
             double xi, r, s;
             getDistance( x, u, i, j, xi, r, s );
 
-            const double coeff =
-                model.thermal_coeff * model.conductivity_function( xi );
+            const double coeff = model.conductivity_function( xi );
             conduction( i ) +=
                 coeff * ( temp( j ) - temp( i ) ) / xi / xi * vol( j );
         };