[Examples] Update speed of sound examples

Updated sound_speed.py and sound_speed_units.py to check against the new built-in function
Cantera · Jun 14, 2023 · 931089b · 931089b
1 parent ae38de4
commit 931089b
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 10 deletions.
diff --git a/samples/python/thermo/sound_speed.py b/samples/python/thermo/sound_speed.py
@@ -1,7 +1,7 @@
 """
 Compute the "equilibrium" and "frozen" sound speeds for a gas
 
-Requires: cantera >= 2.5.0
+Requires: cantera >= 3.0.0
 Keywords: thermodynamics, equilibrium
 """
 
@@ -42,10 +42,8 @@ def equilSoundSpeeds(gas, rtol=1.0e-6, max_iter=5000):
     # equilibrium sound speed
     aequil = math.sqrt((p1 - p0)/(gas.density - r0))
 
-    # compute the frozen sound speed using the ideal gas expression as a check
-    gamma = gas.cp/gas.cv
-    afrozen2 = math.sqrt(gamma * ct.gas_constant * gas.T /
-                         gas.mean_molecular_weight)
+    # check against the built-in sound speed function
+    afrozen2 = gas.sound_speed
 
     return aequil, afrozen, afrozen2
 

diff --git a/samples/python/thermo/sound_speed_units.py b/samples/python/thermo/sound_speed_units.py
@@ -46,11 +46,8 @@ def equilibrium_sound_speeds(gas, rtol=1.0e-6, max_iter=5000):
     # equilibrium sound speed
     aequil = np.sqrt((p1 - p0)/(gas.density - r0)).to("ft/s")
 
-    # compute the frozen sound speed using the ideal gas expression as a check
-    gamma = gas.cp/gas.cv
-    gamma * ctu.units.molar_gas_constant
-    afrozen2 = np.sqrt(gamma * ctu.units.molar_gas_constant * gas.T /
-                         gas.mean_molecular_weight).to("ft/s")
+    # check against the built-in sound speed function
+    afrozen2 = gas.sound_speed.to("ft/s")
 
     return aequil, afrozen, afrozen2