fullflow 2.0.3__tar.gz → 3.0.0__tar.gz

{fullflow-2.0.3 → fullflow-3.0.0}/CHANGELOG.md

@@ -1,5 +1,24 @@
 # Changelog
 
+## 3.0.0
+
+FullFlow 3.0.0 focuses on solver performance, transient robustness, consistent
+examples, and concise user-facing documentation.
+
+### Solver performance
+
+* Added a safeguarded modified-Newton transient solver with numerical Jacobian
+  construction, Broyden updates, backtracking, and automatic SciPy fallback.
+* Added guarded Jacobian reuse between compatible transient timesteps.
+* Added accepted-step prediction and validation to improve transient initial
+  guesses without reusing stale state or control information.
+* Reduced repeated state discovery, allocation, conversion, and component
+  evaluation work in steady-state and transient residual calculations.
+* Added compact state snapshots and cached runtime metadata with network-version
+  invalidation.
+* Improved invalid-trial handling so temporary property-evaluation failures can
+  be rejected without immediately terminating a nonlinear solve.
+
 ## 2.0.3
 
 ### Sequence inputs

{fullflow-2.0.3 → fullflow-3.0.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: fullflow
-Version: 2.0.3
+Version: 3.0.0
 Summary: Python-native fluid, thermal, propulsion, control, and transient network simulation.
 Keywords: fluid-dynamics,fluid-network,network-analysis,thermodynamics,thermal-analysis,heat-transfer,propulsion,rocket-propulsion,rocket-engine,compressible-flow,turbomachinery,control-systems,transient-simulation,simulation,engineering,scientific-computing
 Author: Saaketh Ramoju

{fullflow-2.0.3 → fullflow-3.0.0}/examples/steady_state/0function_lookup.py

@@ -39,6 +39,13 @@ import math
 from fullflow import *
 
 
+# -----------------------------------------------------------------------------
+# Run options
+# -----------------------------------------------------------------------------
+
+filename = '0function_lookup'
+
+
 # -----------------------------------------------------------------------------
 # User-defined callable lookup
 # -----------------------------------------------------------------------------
@@ -141,9 +148,8 @@ ColumnDesign.track("Applied Stress [Pa]", AppliedStress)
 # Solve
 # -----------------------------------------------------------------------------
 
-filename = "0function_lookup"
 
 SteadyState(ColumnDesign).solve(
     verbose=True,
     filename=filename,
-)
+)

{fullflow-2.0.3 → fullflow-3.0.0}/examples/steady_state/2compressible_flow.py

@@ -74,6 +74,7 @@ Demonstrates
 """
 
 import numpy as np
+from scipy.constants import *
 
 from fullflow import *
 from thermoprop import *
@@ -125,11 +126,10 @@ NodeFluid = Lookup(
 # These values are only the initial geometry guesses / fixed downstream geometry.
 # -----------------------------------------------------------------------------
 
-IN_TO_M = 0.0254
 
-D_diff_in = 0.75 * IN_TO_M
-D_diff_out = 1.00 * IN_TO_M
-D_nozzle = 0.75 * IN_TO_M
+D_diff_in = 0.75 * inch
+D_diff_out = 1.00 * inch
+D_nozzle = 0.75 * inch
 
 A_diff_in = np.pi * D_diff_in**2 / 4.0
 A_diff_out = np.pi * D_diff_out**2 / 4.0
@@ -226,4 +226,4 @@ InletBalance = Balance(
 # Solve network
 # -----------------------------------------------------------------------------
 
-SteadyState(CompressibleNetwork).solve(verbose=True)
+SteadyState(CompressibleNetwork).solve(verbose=True)

{fullflow-2.0.3 → fullflow-3.0.0}/examples/steady_state/3custom_component.py

@@ -71,6 +71,7 @@ Demonstrates
 """
 
 import numpy as np
+from scipy.constants import *
 
 from fullflow import *
 from thermoprop import *
@@ -235,7 +236,7 @@ class CustomPump(Component):
         density: State,
         upstream_pressure: State,
         discharge_pressure: State,
-        gravitational_acceleration: float = 9.80665
+        gravitational_acceleration: float = g
     ):
         # The balances property below reads pressure_error.  Initialize it once
         # so the component is safe to inspect before the first evaluate_states()
@@ -531,4 +532,4 @@ LOXPumpSystem.track("Pressure Rise [psid]", dP/6894.76)
 # verbose=True prints the solver summary and final network solution.
 # -----------------------------------------------------------------------------
 
-SteadyState(LOXPumpSystem).solve(verbose=True)
+SteadyState(LOXPumpSystem).solve(verbose=True)

{fullflow-2.0.3 → fullflow-3.0.0}/examples/steady_state/4cavitating_venturi.py

@@ -96,13 +96,19 @@ from fullplot import Axis, generate_map
 from thermoprop import *
 
 
+# -----------------------------------------------------------------------------
+# Run options
+# -----------------------------------------------------------------------------
+
+filename = '4cavitating_venturi'
+
+
 # -----------------------------------------------------------------------------
 # Map file settings
 # -----------------------------------------------------------------------------
 # The chamber gas property map is expected to live in RP1_LOX.h5. The map is
 # generated separately from the solve because equilibrium calculations are much
 # more expensive than interpolation during each network evaluation.
-filename = "4cavitating_venturi"
 
 
 # -----------------------------------------------------------------------------

{fullflow-2.0.3 → fullflow-3.0.0}/examples/steady_state/5heat_exchanger.py

@@ -104,6 +104,7 @@ Approximate cross section:
 """
 
 import numpy as np
+from scipy.constants import *
 
 from fullflow import *
 from thermoprop import *
@@ -233,22 +234,21 @@ gasoline_specific_heat_cp = (
 # The current network uses three hydraulic tube elements per side, so Nseg is 3.
 # -----------------------------------------------------------------------------
 
-IN_TO_M = 0.0254
 
 # Total heat-exchanger length and length of each model segment.
-L_total = 40.0 * IN_TO_M
+L_total = 40.0 * inch
 Nseg = 3
 L = L_total / Nseg
 
 # Inner-tube inside diameter. This is the hot-fluid flow diameter.
-D_inner = 0.50 * IN_TO_M
+D_inner = 0.50 * inch
 
 # Tube wall thickness and resulting tube outside diameter.
-wall_thickness = 0.035 * IN_TO_M
+wall_thickness = 0.035 * inch
 D_tube_outer = D_inner + 2.0 * wall_thickness
 
 # Outer diameter of the annulus. Coolant flows between D_tube_outer and this ID.
-D_annulus_outer = 0.875 * IN_TO_M
+D_annulus_outer = 0.875 * inch
 
 # Basic geometry check to avoid a negative annular area.
 if D_tube_outer >= D_annulus_outer:
@@ -952,15 +952,14 @@ CoolantNode2.heat_rate = -Coolant2Solid1Convection.heat_rate
 # exchanger.
 # -----------------------------------------------------------------------------
 
-PA_TO_PSIA = 1.0 / 6894.757293
 
-HeatExchanger.track("Hot Fluid Source Pressure [psia]", hot_fluid_source_pressure * PA_TO_PSIA)
-HeatExchanger.track("Hot Fluid Node 1 Pressure [psia]", hot_fluid_node1_pressure * PA_TO_PSIA)
-HeatExchanger.track("Hot Fluid Node 2 Pressure [psia]", hot_fluid_node2_pressure * PA_TO_PSIA)
+HeatExchanger.track("Hot Fluid Source Pressure [psia]", hot_fluid_source_pressure * (1.0 / psi))
+HeatExchanger.track("Hot Fluid Node 1 Pressure [psia]", hot_fluid_node1_pressure * (1.0 / psi))
+HeatExchanger.track("Hot Fluid Node 2 Pressure [psia]", hot_fluid_node2_pressure * (1.0 / psi))
 
-HeatExchanger.track("Coolant Source Pressure [psia]", CoolantSource.pressure * PA_TO_PSIA)
-HeatExchanger.track("Coolant Node 1 Pressure [psia]", CoolantNode1Fluid.pressure * PA_TO_PSIA)
-HeatExchanger.track("Coolant Node 2 Pressure [psia]", CoolantNode2Fluid.pressure * PA_TO_PSIA)
+HeatExchanger.track("Coolant Source Pressure [psia]", CoolantSource.pressure * (1.0 / psi))
+HeatExchanger.track("Coolant Node 1 Pressure [psia]", CoolantNode1Fluid.pressure * (1.0 / psi))
+HeatExchanger.track("Coolant Node 2 Pressure [psia]", CoolantNode2Fluid.pressure * (1.0 / psi))
 
 HeatExchanger.track("Hot Fluid Source Temperature [K]", hot_fluid_source_temperature)
 HeatExchanger.track("Hot Fluid Node 1 Temperature [K]", hot_fluid_node1_temperature)

{fullflow-2.0.3 → fullflow-3.0.0}/examples/steady_state/8turbopump.py

@@ -115,6 +115,7 @@ user-generated performance maps rather than detailed hardware-calibrated maps.
 """
 
 import numpy as np
+from scipy.constants import *
 
 from fullflow import *
 from fullplot import Axis, generate_map
@@ -132,8 +133,6 @@ Turbopump = Network("Turbopump")
 # Units
 # -----------------------------------------------------------------------------
 
-bar = 100000.0
-
 
 # -----------------------------------------------------------------------------
 # Design point
@@ -211,8 +210,6 @@ Preburner = Equilibrium(
 
 LOX = Propellant("lox", temperature=90)
 
-g = 9.80665
-
 LOXPumpDensity = State(LOX.density)
 LOXPumpDesignVolumetricFlow = LOXPumpDesignMassFlow / LOX.density
 LOXPumpDesignHead = (LOXPumpDesignDischargePressure - LOXPumpDesignInletPressure) / (LOX.density * g)
@@ -255,8 +252,9 @@ OxRichTurbineDesignMassFlow = LOXPumpDesignShaftPower / (OxRichTurbineDesignEffi
 
 turbopump_map_filename = "8turbopump"
 
+# Set this to True to generate both maps in 8turbopump.h5 before loading them.
+generate_maps = False
 
-'''
 # -----------------------------------------------------------------------------
 # Turbine map generation
 # -----------------------------------------------------------------------------
@@ -327,7 +325,7 @@ def turbine_map_point(
     }
 
 
-generate_map(
+generate_maps and generate_map(
     turbopump_map_filename,
     group="turbine",
     axes=[
@@ -421,7 +419,7 @@ def lox_pump_map_point(rotor_speed, volumetric_flow):
     }
 
 
-generate_map(
+generate_maps and generate_map(
     turbopump_map_filename,
     group="lox_pump",
     axes=[
@@ -452,7 +450,6 @@ generate_map(
     overwrite=True,
     raise_errors=True,
 )
-'''
 
 
 # -----------------------------------------------------------------------------

{fullflow-2.0.3 → fullflow-3.0.0}/examples/steady_state/9equilibrium_nozzle.py

@@ -97,6 +97,7 @@ from concurrent.futures import ProcessPoolExecutor
 import os
 
 import numpy as np
+from scipy.constants import *
 
 from fullflow import *
 from fullplot import Axis, generate_map
@@ -107,7 +108,6 @@ from thermoprop import *
 # Constants and script options
 # ---------------------------------------------------------------------------
 
-psia_to_pa = 6894.76
 
 # All generated map groups are stored in this one HDF5 file.
 filename = "9equilibrium_nozzle"
@@ -275,7 +275,6 @@ def rp1_lox_sp_map(pressure, entropy, mixture_ratio, guess_temperature):
         "gas_constant": eq.gas_constant,
         "density": eq.density,
         "enthalpy": eq.enthalpy,
-        "entropy": eq.entropy,
         "speed_of_sound": eq.speed_of_sound,
         "temperature": eq.temperature,
     }
@@ -312,7 +311,6 @@ def rp1_lox_tp_map(pressure, temperature, mixture_ratio):
         "enthalpy": eq.enthalpy,
         "entropy": eq.entropy,
         "speed_of_sound": eq.speed_of_sound,
-        "temperature": eq.temperature,
     }
 
 
@@ -327,7 +325,7 @@ if __name__ == "__main__":
     # -----------------------------------------------------------------------
 
     mixture_ratio = State(2.3)
-    injector_face_pressure = State(400 * psia_to_pa)
+    injector_face_pressure = State(400 * psi)
     ambient_pressure = State(101325.0)
 
     fuel = Propellant("rp-1", temperature=298.15)
@@ -373,8 +371,8 @@ if __name__ == "__main__":
             axes=[
                 Axis.linear(
                     "pressure",
-                    start=385 * psia_to_pa,
-                    stop=405 * psia_to_pa,
+                    start=385 * psi,
+                    stop=405 * psi,
                     count=41,
                     units="Pa",
                 ),
@@ -408,8 +406,8 @@ if __name__ == "__main__":
             axes=[
                 Axis.log(
                     "pressure",
-                    start=1 * psia_to_pa,
-                    stop=500 * psia_to_pa,
+                    start=1 * psi,
+                    stop=500 * psi,
                     count=18,
                     units="Pa",
                 ),
@@ -451,12 +449,12 @@ if __name__ == "__main__":
 
     # Chamber pressure and temperature are solved by the combustion chamber
     # Volume. These are only initial guesses.
-    chamber_pressure = State(393 * psia_to_pa)
+    chamber_pressure = State(393 * psi)
     chamber_temperature = State(3458.0)
 
     # Throat and exit pressures are also solved by the network.
-    throat_pressure = State(300 * psia_to_pa)
-    exit_pressure = State(10 * psia_to_pa)
+    throat_pressure = State(300 * psi)
+    exit_pressure = State(10 * psi)
 
 
     # -----------------------------------------------------------------------
@@ -689,11 +687,11 @@ if __name__ == "__main__":
     # Tracked outputs
     # -----------------------------------------------------------------------
 
-    EquilibriumNozzle.track("Injector Face Pressure [psia]", InjectorFace.pressure / psia_to_pa)
-    EquilibriumNozzle.track("Chamber Pressure [psia]", chamber_pressure / psia_to_pa)
-    EquilibriumNozzle.track("Throat Pressure [psia]", throat_pressure / psia_to_pa)
-    EquilibriumNozzle.track("Exit Pressure [psia]", exit_pressure / psia_to_pa)
-    EquilibriumNozzle.track("Ambient Pressure [psia]", ambient_pressure / psia_to_pa)
+    EquilibriumNozzle.track("Injector Face Pressure [psia]", InjectorFace.pressure / psi)
+    EquilibriumNozzle.track("Chamber Pressure [psia]", chamber_pressure / psi)
+    EquilibriumNozzle.track("Throat Pressure [psia]", throat_pressure / psi)
+    EquilibriumNozzle.track("Exit Pressure [psia]", exit_pressure / psi)
+    EquilibriumNozzle.track("Ambient Pressure [psia]", ambient_pressure / psi)
 
     EquilibriumNozzle.track("Injector Face Temperature [K]", InjectorFace.temperature)
     EquilibriumNozzle.track("Chamber Temperature [K]", chamber_temperature)
@@ -736,4 +734,4 @@ if __name__ == "__main__":
 
     print(f"Predicted Line Loss: {ratio_predicted.value:.3f}")
     print(f"Rayleigh Line Loss: {ratio_rayleigh.value:.3f}")
-    print(f"Empirical Loss: {ratio_empirical:.3f}")
+    print(f"Empirical Loss: {ratio_empirical:.3f}")

{fullflow-2.0.3 → fullflow-3.0.0}/examples/testing/1test_data_anchoring.py

@@ -3,6 +3,7 @@ import math
 
 import h5py
 import numpy as np
+from scipy.constants import *
 
 from fullflow import *
 import fullplot as fplt
@@ -51,6 +52,14 @@ loss, heat input, etc.
 """
 
 
+# -----------------------------------------------------------------------------
+# Run options
+# -----------------------------------------------------------------------------
+
+dt = 0.02
+t_final = 6.5
+
+
 # ---------------------------------------------------------------------------
 # User options
 # ---------------------------------------------------------------------------
@@ -63,7 +72,6 @@ HERE = Path(__file__).resolve().parent
 TEST_DATA_FILE = HERE / "test_data_anchoring_hotfire.h5"
 RESULT_FILE = HERE / "test_data_anchoring_result.h5"
 
-psia_to_pa = 6894.76
 
 
 # ---------------------------------------------------------------------------
@@ -112,7 +120,7 @@ def generate_test_data(filename):
     density = 810.0
     area = (math.pi / 4.0) * (0.25 / 39.37) ** 2
     true_cd = 0.18 + 0.54 * run_level
-    pressure_drop = np.maximum((ftpt - fipt) * psia_to_pa, 0.0)
+    pressure_drop = np.maximum((ftpt - fipt) * psi, 0.0)
     fuel_mdot = true_cd * area * np.sqrt(2.0 * density * pressure_drop)
     fuel_mdot += rng.normal(0.0, 0.006, len(time))
     fuel_mdot = np.maximum(fuel_mdot, 0.0)
@@ -155,8 +163,8 @@ fflow = data.trace(y="FUEL_MDOT", x=time, name="Fuel Flow Test Data", role="data
 
 # A user can process a test-data trace before anchoring to it. The filtered trace
 # keeps the same role as the source trace, so it is still valid Sensor.data.
-ftpt = ftpt_raw.filter("moving_average", window=0.04) * psia_to_pa
-fipt = fipt_raw.filter("moving_average", window=0.04) * psia_to_pa
+ftpt = ftpt_raw.filter("moving_average", window=0.04) * psi
+fipt = fipt_raw.filter("moving_average", window=0.04) * psi
 
 # All three traces share this TimeAxis, so one window selects the useful part of
 # the test for every attached trace.
@@ -181,7 +189,7 @@ fflow_redline = fplt.Trace.constant(
 # Plot the imported data before building the model. This is the normal first
 # sanity check when working with test data.
 data.plot(
-    y=[ftpt / psia_to_pa, fipt / psia_to_pa],
+    y=[ftpt / psi, fipt / psi],
     y2=[fflow, fflow_warning, fflow_redline],
     xlabel="Time [s]",
     ylabel="Pressure [psia]",
@@ -255,8 +263,8 @@ Test.track("Estimated Fuel Cd", Line.discharge_coefficient)
 # ---------------------------------------------------------------------------
 
 Transient(Test).solve(
-    dt=0.02,
-    t_final=6.5,
+    dt=dt,
+    t_final=t_final,
     verbose=True,
     filename=str(RESULT_FILE),
 )

{fullflow-2.0.3 → fullflow-3.0.0}/examples/testing/2sensor_events.py

@@ -32,6 +32,14 @@ Sequence commands and aborts can respond to those events.
 """
 
 
+# -----------------------------------------------------------------------------
+# Run options
+# -----------------------------------------------------------------------------
+
+dt = 0.02
+t_final = 5.2
+
+
 HERE = Path(__file__).resolve().parent
 RESULT_FILE = HERE / "sensor_events_result.h5"
 
@@ -109,8 +117,8 @@ Test.track("CHPT", chamber_pressure)
 # ---------------------------------------------------------------------------
 
 Transient(Test).solve(
-    dt=0.02,
-    t_final=5.2,
+    dt=dt,
+    t_final=t_final,
     verbose=True,
     filename=str(RESULT_FILE),
 )

{fullflow-2.0.3 → fullflow-3.0.0}/examples/testing/3sequences_and_commands.py

@@ -35,6 +35,14 @@ combustion, nozzles, turbomachinery, etc. The Sequence API is the same.
 """
 
 
+# -----------------------------------------------------------------------------
+# Run options
+# -----------------------------------------------------------------------------
+
+dt = 0.01
+t_final = 5.0
+
+
 HERE = Path(__file__).resolve().parent
 RESULT_FILE = HERE / "sequence_commands_abort_result.h5"
 
@@ -182,8 +190,8 @@ Test.track("Chamber Pressure", chamber_pressure)
 # ---------------------------------------------------------------------------
 
 Transient(Test).solve(
-    dt=0.01,
-    t_final=5.0,
+    dt=dt,
+    t_final=t_final,
     verbose=True,
     filename=str(RESULT_FILE),
 )

{fullflow-2.0.3 → fullflow-3.0.0}/examples/transient/10custom_component.py

@@ -48,6 +48,15 @@ monotonically drains.
 from fullflow import *
 
 
+# -----------------------------------------------------------------------------
+# Run options
+# -----------------------------------------------------------------------------
+
+dt = 0.05
+t_final = 20.0
+filename = '10custom_component'
+
+
 # ---------------------------------------------------------------------------
 # Custom component 1: tank dynamics
 # ---------------------------------------------------------------------------
@@ -159,8 +168,8 @@ TankNetwork.track("Outlet Flow", outlet_flow)
 # Transient solve
 # ---------------------------------------------------------------------------
 Transient(TankNetwork).solve(
-    dt=0.05,
-    t_final=20.0,
-    filename="10custom_component",
+    dt=dt,
+    t_final=t_final,
+    filename=filename,
     verbose=True,
 )

{fullflow-2.0.3 → fullflow-3.0.0}/examples/transient/11heat_transfer.py

@@ -59,12 +59,20 @@ from fullflow import *
 from thermoprop import *
 
 
+# -----------------------------------------------------------------------------
+# Run options
+# -----------------------------------------------------------------------------
+
+dt = 1.0
+t_final = 900.0
+filename = '11heat_transfer'
+
+
 # ---------------------------------------------------------------------------
 # Network
 # ---------------------------------------------------------------------------
 HeatTransferNetwork = Network("Solid Heat Transfer Transient")
 
-filename = "11heat_transfer"
 
 
 # ---------------------------------------------------------------------------
@@ -203,9 +211,9 @@ HeatTransferNetwork.track("Biot Number", biot_number)
 # Transient solve
 # ---------------------------------------------------------------------------
 Transient(HeatTransferNetwork).solve(
-    dt=1.0,
-    t_final=900.0,
+    dt=dt,
+    t_final=t_final,
     filename=filename,
     save_dt=5.0,
     verbose=True,
-)
+)

{fullflow-2.0.3 → fullflow-3.0.0}/examples/transient/12heat_exchanger.py

@@ -137,16 +137,25 @@ import re
 import h5py
 import matplotlib.pyplot as plt
 import numpy as np
+from scipy.constants import *
 
 from fullflow import *
 from thermoprop import *
 
 
+# -----------------------------------------------------------------------------
+# Run options
+# -----------------------------------------------------------------------------
+
+dt = 0.1
+t_final = 4.0
+filename = '12heat_exchanger'
+
+
 # -----------------------------------------------------------------------------
 # Output file
 # -----------------------------------------------------------------------------
 
-filename = "12heat_exchanger"
 
 
 def harmonic_mean(a, b):
@@ -273,22 +282,21 @@ gasoline_specific_heat_cp = (
 # The current network uses three hydraulic tube elements per side, so Nseg is 3.
 # -----------------------------------------------------------------------------
 
-IN_TO_M = 0.0254
 
 # Total heat-exchanger length and length of each model segment.
-L_total = 40.0 * IN_TO_M
+L_total = 40.0 * inch
 Nseg = 3
 L = L_total / Nseg
 
 # Inner-tube inside diameter. This is the hot-fluid flow diameter.
-D_inner = 0.50 * IN_TO_M
+D_inner = 0.50 * inch
 
 # Tube wall thickness and resulting tube outside diameter.
-wall_thickness = 0.035 * IN_TO_M
+wall_thickness = 0.035 * inch
 D_tube_outer = D_inner + 2.0 * wall_thickness
 
 # Outer diameter of the annulus. Coolant flows between D_tube_outer and this ID.
-D_annulus_outer = 0.875 * IN_TO_M
+D_annulus_outer = 0.875 * inch
 
 # Basic geometry check to avoid a negative annular area.
 if D_tube_outer >= D_annulus_outer:
@@ -1040,15 +1048,14 @@ CoolantNode2.heat_rate = -Coolant2Solid1Convection.heat_rate
 # exchanger.
 # -----------------------------------------------------------------------------
 
-PA_TO_PSIA = 1.0 / 6894.757293
 
-HeatExchanger.track("Hot Fluid Source Pressure [psia]", hot_fluid_source_pressure * PA_TO_PSIA)
-HeatExchanger.track("Hot Fluid Node 1 Pressure [psia]", hot_fluid_node1_pressure * PA_TO_PSIA)
-HeatExchanger.track("Hot Fluid Node 2 Pressure [psia]", hot_fluid_node2_pressure * PA_TO_PSIA)
+HeatExchanger.track("Hot Fluid Source Pressure [psia]", hot_fluid_source_pressure * (1.0 / psi))
+HeatExchanger.track("Hot Fluid Node 1 Pressure [psia]", hot_fluid_node1_pressure * (1.0 / psi))
+HeatExchanger.track("Hot Fluid Node 2 Pressure [psia]", hot_fluid_node2_pressure * (1.0 / psi))
 
-HeatExchanger.track("Coolant Source Pressure [psia]", CoolantSource.pressure * PA_TO_PSIA)
-HeatExchanger.track("Coolant Node 1 Pressure [psia]", CoolantNode1Fluid.pressure * PA_TO_PSIA)
-HeatExchanger.track("Coolant Node 2 Pressure [psia]", CoolantNode2Fluid.pressure * PA_TO_PSIA)
+HeatExchanger.track("Coolant Source Pressure [psia]", CoolantSource.pressure * (1.0 / psi))
+HeatExchanger.track("Coolant Node 1 Pressure [psia]", CoolantNode1Fluid.pressure * (1.0 / psi))
+HeatExchanger.track("Coolant Node 2 Pressure [psia]", CoolantNode2Fluid.pressure * (1.0 / psi))
 
 HeatExchanger.track("Hot Fluid Source Temperature [K]", hot_fluid_source_temperature)
 HeatExchanger.track("Hot Fluid Node 1 Temperature [K]", hot_fluid_node1_temperature)
@@ -1104,10 +1111,9 @@ SteadyState(HeatExchanger).solve(
 # A small timestep is used because the fluid residence time in this short,
 # small-diameter heat exchanger is much faster than the wall thermal response.
 Transient(HeatExchanger).solve(
-    dt=0.01,
-    t_final=4.0,
-    save_dt=0.02,
+    dt=dt,
+    t_final=t_final,
     filename=filename,
     verbose=True,
     statistics=True,
-)
+)