adjusted time input for heat() to accept relative time

kaelyndunnell · Nov 1, 2024 · 76d0aca · 76d0aca
1 parent 8de5ae1
commit 76d0aca
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/mb_scenario.py b/mb_scenario.py
@@ -270,7 +270,9 @@ def T_surface(t: dolfinx.fem.Constant) -> float:
             monoblock surface temperature in K
         """
         pulse_type = my_scenario.get_pulse_type(float(t))
-        return 1.1e-4 * heat(pulse_type, t=t) + COOLANT_TEMP
+        pulse_row = my_scenario.get_row(float(t))
+        time_elapsed = my_scenario.get_time_till_row(pulse_row)
+        return 1.1e-4 * heat(pulse_type, t=t-time_elapsed) + COOLANT_TEMP
 
     def T_rear(t: dolfinx.fem.Constant):
         """Monoblock surface temperature
@@ -282,7 +284,9 @@ def T_rear(t: dolfinx.fem.Constant):
             monoblock surface temperature in K
         """
         pulse_type = my_scenario.get_pulse_type(float(t))
-        return 2.2e-5 * heat(pulse_type, t=t) + COOLANT_TEMP
+        pulse_row = my_scenario.get_row(float(t))
+        time_elapsed = my_scenario.get_time_till_row(pulse_row)
+        return 2.2e-5 * heat(pulse_type, t=t-time_elapsed) + COOLANT_TEMP
 
     def T_function(x, t: Constant):
         """Monoblock temperature function