Utilities#

A collection of routines with multiple purposes organized as follows:

unclasified
databases
display
eos
export
extensions
root_finding
thermo
validations

Unclasified utility functions#

A collection of unclasified functions necessary to perform all the calculations in CT.

Database functions#

A collection of functions necessary to obtain generate the databases in CT.

Display functions#

A collection of functions necessary to display the results (command window and plots).

Routines

displaysweepresults(self, mix, xvar, varargin)#

Plot a given variable against the molar fractions of the mixture

Parameters:

self (struct) – Data of the mixture, conditions, and databases
mix (struct) – Properties of the mixture
xvar (float) – Vector with the x data

Optional Args:

config (struct): Struct with the configuration for plots
xscale (string): Scale fot the x-axis (linear or log)
yscale (string): Scale fot the y-axis (linear or log)
xdir (string): Direction of the x-axis (normal or reverse)
ydir (string): Direction of the x-axis (normal or reverse)

Returns:: ax (axes) – Axes object

get_mixtures(PS, pattern)#

Get all non-empty mixture

Parameters:

PS (struct) – Struct with all the data of Problem Solution (PS)
pattern (str) – Pattern/s name of the mixture

Returns:

mixtures (cell) – Cell with all the non-empty mixtures

Examples

mixtures = get_mixtures(self.PS, ‘mix’);
mixtures = get_mixtures(self.PS, ‘strP’);

interpreter_label(property, varargin)#

Interpreter label for properties - returns property name

Note

The ‘interpreter_label.m’ routine considers that the properties are in mass basis. This will be fixed in a future patch.

Parameters:: property (char) – Property name

Optional Args:: type (char): Type of label to return. Can be ‘short’, ‘medium’ or ‘long’ (default: medium)

Returns:: value (char) – Corresponding name of the property

plot_figure(x_field, x_var, y_field, y_var, varargin)#

Plot figure with customizable settings

Note

The ‘interpreter_label.m’ routine considers that the properties are in mass basis. This will be fixed in a future patch.

Parameters:

x_field (char) – Field name for the x-axis data
x_var (cell) – Cell array containing the x-axis data
y_field (char) – Field name for the y-axis data
y_var (cell) – Cell array containing the y-axis data

Optional Name-Value Pair Args:

config (struct): Struct with the configuration for plots
leg or legend (cell): Cell array of strings containing the legend names
legend_location (char): Location of the legend
ax or axes (object): Handle of the axes to plot on
linestyle (char): Line style
linewidth (float): Line width
fontsize (float): Font size
title (char): Title of the figure
labelx, xlabel, label_x, or x_label (char): x-axis label
labely, ylabel, label_y, or ‘y_label (char): y-axis label
label_type (char): Label type
xscale (char): Set x-axis scale (linear or log)
yscale (char): Set y-axis scale (linear or log)
xdir (char): Set x-axis direction (normal or reverse)
ydir (char): Set y-axis direction (normal or reverse)
color (float): Line color [R, G, B]

Returns:

Tuple containing

ax (object): Handle of the axes
dline (object): Handle of the plotted line

plot_figure_set(range_name, range, properties, mix, varargin)#

Plot a set of properties in a tiled layout figure

Parameters:

range_name (char) – Variable name of the x-axis parameter
range (float) – Vector x-axis values
properties (cell) – Cell array of properties to plot
mix (struct) – Mixture

Optional Args:

ax (handle): Handle to the main axis
config (struct): Configuration settings for the figure
basis (cell): Cell array with the basis for each property

Returns:: ** main_ax (obj)* – Handle to the main axis

Examples

plot_figure_set(‘T’, 300:100:1000, {‘cp’, ‘cv’, ‘h’, ‘s’}, mix)
plot_figure_set(‘T’, 300:100:1000, {‘cp’, ‘cv’, ‘h’, ‘s’}, mix, ‘config’, config)
plot_figure_set(‘T’, 300:100:1000, {‘cp’, ‘cv’, ‘h’, ‘s’}, mix, ‘config’, config, ‘basis’, {‘’, ‘’, ‘mi’, ‘mi’});
plot_figure_set(‘T’, 300:100:1000, {‘cp’, ‘cv’, ‘h’, ‘s’}, mix, ‘config’, config, ‘basis’, {‘’, ‘’, ‘mi’, ‘mi’}, ‘ax’, ax);

plot_hugoniot(self, mix1, mix2, varargin)#

Plot the Hugoniot curve for a a given pre-shock state (mix1) and post-shock state (mix2)

Parameters:

self (struct) – Data of the mixture, conditions, and databases
mix1 (struct) – Pre-shock mixture
mix2 (struct) – Post-shock mixture

Optional Args:: ax (object): Axis handle to plot on

Returns:: ax (object) – Axis handle to plot on

Examples

ax = plot_hugoniot(self, mix1, mix2)
ax = plot_hugoniot(self, mix1, mix2, ax)

plot_molar_fractions(self, x_var, x_field, y_field, varargin)#

Plot molar fractions againts any variable

Parameters:

self (struct) – Data of the mixture, conditions, and databases
x_var (cell) – Properties of the mixture for all the cases
x_field (char) – Field name for the x-axis data
y_field (char) – Field name for the y-axis data

Optional Name-Value Pair Args:

validation (struct): Struct that contains validations with (x_field, y_field)
nfrec (float): Frequency points to plot validations
mintol (float): Minimum limit i-axis with the composition of the mixture
display_species (cell): List of species to plot
y_var (cell): Get y-axis data from a different mixture
config (struct): Struct with the configuration for plots
axis_x (char): Set x-axis limits
axis_y (char): Set y-axis limits
xscale (char): Set x-axis scale (linear or log)
yscale (char): Set y-axis scale (linear or log)
xdir (char): Set x-axis direction (normal or reverse)
ydir (char): Set y-axis direction (normal or reverse)
ax (object): Handle of the axes to plot on

Returns:

Tuple containing

ax (object): Handle of the axes
fig (object): Handle of the figure

Examples

[ax, fig] = plot_molar_fractions(self, mix1, ‘phi’, ‘Xi’)
[ax, fig] = plot_molar_fractions(self, mix1, ‘phi’, ‘Xi’, ‘y_var’, mix2)
[ax, fig] = plot_molar_fractions(self, mix1, ‘phi’, ‘Xi’, ‘y_var’, mix2, ‘validation’, results_CEA)
[ax, fig] = plot_molar_fractions(self, mix1, ‘phi’, ‘Xi’, ‘y_var’, mix2, ‘validation’, results_CEA, ‘display_species’, display_species)

plot_shock_polar(varargin)#

Routine to obtain shock polar plots

Plot (pressure, deflection)
Plot (wave angle, deflection)
Plot velocity components

polynomial_regression(x, y, n)#

Obtain polynomial regression for the given dataset (x, y) and polynomial order

Parameters:

x (float) – x values
y (float) – y values
n (float) – polynomial order

Returns:

y_poly (float) – y values of the polynomial regression

print_mixture(self, varargin)#

Print properties and composition of the given mixtures in the command window

Parameters:: self (struct) – Data of the mixture, conditions, and databases

Optional Args:

mix1 (struct): Struct with the properties of the mixture
mix2 (struct): Struct with the properties of the mixture
mixi (struct): Struct with the properties of the mixture
mixN (struct): Struct with the properties of the mixture

Examples

print_mixture(self, mix1)
print_mixture(self, mix1, mix2)
print_mixture(self, mix1, mix2, mix3)
print_mixture(self, mix1, mix2, mix3, mix4)

print_stoichiometric_matrix(self, varargin)#

Print stoichiometric matrix

Parameters:: self (struct) – Data of the mixture, conditions, and databases

Optional args:

type (char): ‘transpose’

Optional returns:

A0 (table): Stoichiometric matrix. In case type == ‘transpose’: it returns the transpose of stoichiometric matrix

results(self, i)#

Display results in the command window

Parameters:

self (struct) – Data of the mixture, conditions, and databases
i (float) – index of the current iteration

set_figure(varargin)#

Initialize figure with a standard composition

Optional Args:

ax (axis): Figure axis
config (struct): Struct with default plot parameters

Returns:

Tuple containing

ax (axis): Axis of the standard figure
config (struct): Struct with default plot parameters
fig (figure): Standard figure

set_legends(ax, legend_name, varargin)#

Set legend to the given axes

Parameters:

ax (object) – Handle to the axes
legend_name (cell) – Cell array of char containing the legend names

Optional Name-Value Pairs Args:

config (struct): Struct containing the configuration parameters for the plots
obj (object): Handle to the plotted objects (e.g. lines, patches, etc.)

set_title(ax, varargin)#: Set legend to the given axes

species2latex(species, varargin)#

Convert species name into LateX format

Parameters:: species (char) – Species name

Optional Args:: FLAG_BURCAT (bool): If true, do not remove Burcat database prefix (default: true)

Returns:: speciesLatex (char) – Species name in LateX format

Examples

species2latex(‘H2ObLb’)
species2latex(‘Si2H6_M’)
species2latex(‘Si2H6_M’, false)

Equation of State functions#

A collection of Equation of States (EoS) implemented in CT.

Routines

eos_PengRobinson(self, T, p, species, Xi)#

Compute molar volume of the mixture considering Peng-Robinson Equation of State (EoS)

Parameters:

self (struct) – Data of the mixture, conditions, and databases
species (cell) – List of the species in the mixture
Xi (float) – Molar fractions of the mixture
T (float) – Temperature [K]
p (float) – Pressure [bar]

Returns:

Tuple containing

V (float): Molar volume of the mixture [m3/mol]
Vi (float): Molar volume of the components [m3/mol]

eos_VanderWaals(self, species, Xi, T, p)#

Compute molar volume of the mixture considering Van der Waals Equation of State (EoS)

Parameters:

self (struct) – Data of the mixture, conditions, and databases
species (cell) – List of the species in the mixture
Xi (float) – Molar fractions of the mixture
T (float) – Temperature [K]
p (float) – Pressure [bar]

Returns:

Tuple containing

V (float): Molar volume of the mixture [m3/mol]
Vi (float): Molar volume of the components [m3/mol]

eos_ideal(self, T, p, varargin)#

Compute pressure considering ideal Equation of State (EoS)

Parameters:

self (struct) – Data of the mixture, conditions, and databases
T (float) – Temperature of the mixture [K]
p (float) – Pressure of the mixture [Pa]

Returns:

V (float) – Molar volume of the mixture [m3/mol]

eos_ideal_p(self, n, T, v, varargin)#

Compute pressure considering ideal Equation of State (EoS)

Parameters:

self (struct) – Data of the mixture, conditions, and databases
n (float) – Number of moles of the mixture in gaseous phase [mol]
T (float) – Temperature of the mixture [K]
v (float) – Volume of the mixture [m3]

Returns:

p (float) – Pressure of the mixture [Pa]

mu_ex_eos(self, Xi, T, p, V, a_mix, b_mix, a, b)#

Compute non ideal contribution (excess) of the chemical potential assuming cubic Equation of State (EoS) [J/mol]

Parameters:

self (struct) – Data of the mixture, conditions, and databases
T (float) – Temperature of the mixture [K]
p (float) – Pressure of the mixture [bar]
V (float) – Molar volume of the mixture [m3/mol]
a_mix (float) – Atraction factor mixture of the cubic EoS
b_mix (float) – Repulsion factor mixture of the cubic EoS
a (float) – Atraction factor components of the cubic EoS
b (float) – Repulsion factor components of the cubic EoS

Returns:

chemical_potential_ex (float) – Chemical potential excess [J/mol]

mu_ex_ideal(self, moles, temperature, volume)#

Compute non ideal contribution (excess) of the chemical potential assuming ideal Equation of State [J/mol]

Parameters:

self (struct) – Data of the mixture, conditions, and databases
moles (float) – Number of moles of the mixture in gaseous phase [mol]
temperature (float) – Temperature of the mixture [K]
volume (float) – Volume of the mixture [m3]

Returns:

pressure (float) – Pressure of the mixture [Pa]

mu_ex_vanderwaals(self, moles, temperature, volume)#

Compute non ideal contribution (excess) of the chemical potential assuming Van der Waal’s Equation of State [J/mol]

Parameters:

self (struct) – Data of the mixture, conditions, and databases
moles (float) – Number of moles of the mixture in gaseous phase [mol]
temperature (float) – Temperature of the mixture [K]
volume (float) – Volume of the mixture [m3]

Returns:

pressure (float) – Pressure of the mixture [Pa]

mu_ex_virial(self, moles, temperature, volume)#

Compute non ideal contribution (excess) of the chemical potential assuming Virial Equation of State [J/mol]

Parameters:

self (struct) – Data of the mixture, conditions, and databases
moles (float) – Number of moles of the mixture in gaseous phase [mol]
temperature (float) – Temperature of the mixture [K]
volume (float) – Volume of the mixture [m3]

Returns:

pressure (float) – Pressure of the mixture [Pa]

Export functions#

A collection of functions to export results.

Extensions functions#

A collection of external functions from other repositories.

Combustion Toolbox’s color palette is obtained from the following repository: Stephen (2021). ColorBrewer: Attractive and Distinctive Colormaps (https://github.com/DrosteEffect/BrewerMap), GitHub. Retrieved December 3, 2021.
For validations, Combustion Toolbox uses CPU Info from the following repository: Ben Tordoff (2022). CPU Info (https://github.com/BJTor/CPUInfo/releases/tag/v1.3), GitHub. Retrieved March 22, 2022.
Combustion Toolbox’s splash screen is based on a routine from the following repository: Ben Tordoff (2022). SplashScreen (https://www.mathworks.com/matlabcentral/fileexchange/30508-splashscreen), MATLAB Central File Exchange. Retrieved October 15, 2022.

Root finding algorithms#

Roots algorithm used to obtain the temperature at equilibrium for a given thermochemical transformation. The methods implemented are:

Newton-Raphson method.
Steffensen-Aitken method.

Routines

newton_2(f, fprime, x0)#

Find the temperature [K] (root) for the set chemical transformation at equilibrium using the Newton-Raphson method

Parameters:

self (struct) – Data of the mixture, conditions, and databases
mix1 (struct) – Properties of the initial mixture
pP (float) – Pressure [bar]
field (str) – Fieldname in Problem Description (PD)
x0 (float) – Guess temperature [K]
guess_moles (float) – Guess moles final mixture

Returns:

Tuple containing

x (float): Temperature at equilibrium [K]
STOP (float): Relative error [-]
guess_moles (struct): Guess moles final mixture

print_error_root(it, itMax, T, STOP)#

Print error of the method if the number of iterations is greater than maximum iterations allowed

Parameters:

it (float) – Number of iterations executed in the method
itMax (float) – Maximum nNumber of iterations allowed in the method
T (float) – Temperature [K]
STOP (float) – Relative error [-]

regula_guess(self, mix1, pP, field)#

Find a estimate of the temperature for the set chemical equilibrium transformation using the regula falsi method

Parameters:

self (struct) – Data of the mixture, conditions, and databases
mix1 (struct) – Properties of the initial mixture
pP (float) – Pressure [bar]
field (str) – Fieldname in Problem Description (PD)

Returns:

x0 (float) – Guess temperature [K]

newton(self, mix1, pP, field, x0, guess_moles)#

Find the temperature [K] (root) for the set chemical transformation at equilibrium using the second-order Newton-Raphson method

Parameters:

self (struct) – Data of the mixture, conditions, and databases
mix1 (struct) – Properties of the initial mixture
pP (float) – Pressure [bar]
field (str) – Fieldname in Problem Description (PD)
x0 (float) – Guess temperature [K]
guess_moles (float) – Guess moles final mixture

Returns:

Tuple containing

x (float): Temperature at equilibrium [K]
STOP (float): Relative error [-]
guess_moles (struct): Guess moles final mixture

get_gpoint(self, mix1, pP, field, x0, guess_moles)#

Get fixed point of a function based on the chemical transformation

Parameters:

self (struct) – Data of the mixture, conditions, and databases
mix1 (struct) – Properties of the initial mixture
pP (float) – Pressure [bar]
field (str) – Fieldname in Problem Description (PD)
x0 (float) – Guess temperature [K]

Returns:

Tuple containing

gpoint (float): Fixed point of the function [kJ] (HP, EV) or [kJ/K] (SP, SV)
gpoint_relative (float): Fixed relative point of the function [kJ] (HP, EV) or [kJ/K] (SP, SV)

get_point(x_vector, f_vector)#

Get point using the regula falsi method

Parameters:

x_vector (float) – Guess temperature [K]
f_vector (struct) – evaluated functions [kJ] (HP, EV) or [kJ/K] (SP, SV)

Returns:

point (float) – Point of the function [K]

get_point_aitken(x0, g_vector)#

Get fixed point of a function based on the chemical transformation using the Aitken acceleration method

Parameters:

x0 (float) – Guess temperature [K]
g_vector (struct) – Fixed points of the function [kJ] (HP, EV) or [kJ/K] (SP, SV)

Returns:

point (float) – Point of the function [K]

steff(self, mix1, pP, field, x0, guess_moles)#

Find the temperature [K] (root) for the set chemical transformation at equilibrium using the Steffenson-Aitken method

Parameters:

self (struct) – Data of the mixture, conditions, and databases
mix1 (struct) – Properties of the initial mixture
pP (float) – Pressure [bar]
field (str) – Fieldname in Problem Description (PD)
x0 (float) – Guess temperature [K]
guess_moles (float) – Guess moles final mixture

Returns:

Tuple containing

x (float): Temperature at equilibrium [K]
STOP (float): Relative error [-]
guess_moles (float): Guess moles final mixture

Thermodynamic properties#

Functions to obtain thermodynamic properties from a given mixture.

Routines

MolecularWeight(mix)#

Get the molecular weight [g/mol] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Molecular weight [g/mol] of the mixture

adiabaticIndex(mix)#

Get the adiabatic index [-] of the mixture from the ratio of the specific heat capacities

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Adiabatic index of the mixture [-]

adiabaticIndex_sound(mix)#

Get the adiabatic index [-] of the mixture from definition of sound velocity

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Adiabatic index [-] of the mixture

compressibility_factor(mix)#

Compute compressibility factor of the mixture [-]

Parameters:: mix (struct) – Properties of the mixture
Returns:: Z (float) – Compressibility factor [-]

compute_heatrelease(mix1, mix2)#

Compute heat release [J/kg] of the chemical transformation of the mixture 1 to mixture 2

Parameters:

mix1 (struct) – Properties of the initial mixture
mix2 (struct) – Properties of the final mixture

Returns:

q (float) – heat release [J/kg] == [m^2/s^2] Q (float): dimensionless heat release

compute_sound(T, p, species, composition, varargin)#

Routine to compute sound velocity [m/s] for a given temperature-pressure profile

Parameters:

T (float) – Temperature [K]
p (float) – Pressure [bar]
species (cell) – List of species
composition (float) – composition of species (mol)

Optional Name-Value Pairs Args:: self (struct): Data of the mixtures, conditions, databases

Returns:: sound (float) – Sound velocity [m/s]

Examples

sound = compute_sound(300, 1, {‘H2’, ‘O2’}, [1, 1]) sound = compute_sound(300, 1, {‘H2’, ‘O2’}, [1, 1], ‘self’, self)

cp_mass(mix)#

Get the mass-basis specific heat at constant pressure [kJ/kg-K] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mass-basis specific heat at constant pressure [kJ/kg-K] of the mixture

cp_mole(mix)#

Get the mole-basis specific heat at constant pressure [kJ/mol-K] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mole-basis specific heat at constant pressure [kJ/mol-K] of the mixture

cv_mass(mix)#

Get the mass-basis specific heat at constant volume [kJ/kg-K] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mass-basis specific heat at constant volume [kJ/kg-K] of the mixture

cv_mole(mix)#

Get the mole-basis specific heat at constant volume [kJ/mol-K] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mole-basis specific heat at constant volume [kJ/mol-K] of the mixture

density(mix)#

Get the density [kg/m3] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – density [kg/m3] of the mixture

enthalpy_formation_mass(mix)#

Get the mass specific enthalpy formation [kJ/kg] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mass-basis specific enthalpy formation [kJ/kg] of the mixture

enthalpy_formation_mole(mix)#

Get the mole specific enthalpy formation [kJ/mol] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mole-basis specific enthalpy formation [kJ/mol] of the mixture

enthalpy_mass(mix)#

Get the mass specific enthalpy [kJ/kg] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mass-basis specific enthalpy [kJ/kg] of the mixture

enthalpy_mole(mix)#

Get the mole specific enthalpy [kJ/mol] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mole-basis specific enthalpy [kJ/mol] of the mixture

entropy_mass(mix)#

Get the mass specific entropy [kJ/kg-K] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mass-basis specific entropy [kJ/kg-K] of the mixture

entropy_mole(mix)#

Get the mole specific entropy [kJ/mol-K] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mole-basis specific entropy [kJ/mol-K] of the mixture

equivalenceRatio(mix)#

Get the equivalence ratio of the initial mixture [-]

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Equivalence ratio of the initial mixture [-]

gibbs_mass(mix)#

Get the mass specific gibbs free energy [kJ/kg] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mass-basis specific gibbs free energy [kJ/kg] of the mixture

gibbs_mole(mix)#

Get the mole specific gibbs free energy [kJ/mol] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mole-basis specific gibbs free energy [kJ/mol] of the mixture

humidity_specific(T, p, humidity_relative)#

Get the specific humidity of air [kg_w/kg_da] at a given temperature, pressure, and relative humidity

Parameters:

T (float) – Temperature [K]
p (float) – Pressure [bar]
humidity_relative (float) – Relative humidity [%]

Returns:

value (float) – Specific humidity of air [kg_w/kg_da]

intEnergy_mass(mix)#

Get the mass specific internal energy [kJ/kg] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mass-basis specific internal energy [kJ/kg] of the mixture

intEnergy_mole(mix)#

Get the mole specific internal energy [kJ/mol] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mole-basis specific internal energy [kJ/mol] of the mixture

mass(mix)#

Get the mass [kg] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – mass [kg] of the mixture

massFractions(mix)#

Get the mass fractions of all the species in the mixture [-]

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mass fractions of all the species in the mixture [-]

meanMolecularWeight(mix)#

Get the mean molecular weight [g/mol] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mean molecular weight [g/mol] of the mixture

moleFractions(mix)#

Get the mole fractions of all the species in the mixture [-]

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Mole fractions of all the species in the mixture [-]

moles(mix)#

Get the moles [mol] of all the species in the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Moles [mol] of all the species in the mixture

molesGas(mix)#: Get the moles of the gases in the mixture [mol]

pressure(mix)#

Get the pressure [bar] in the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Pressure [bar] in the mixture

soundspeed(mix)#

Get the speed of sound [m/s] in the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Speed of sound [m/s] in the mixture

temperature(mix)#

Get the temperature [K] in the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Temperature [K] in the mixture

velocity_relative(mix)#

Get the velocity of the gases relative to the shock front [m/s] in the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – Velocity of the gases relative to the shock front [m/s] in the mixture

volume(mix)#

Get the volume [m3] of the mixture

Parameters:: mix (struct) – Properties of the mixture
Returns:: value (float) – volume [m3] of the mixture

Validations functions#

A collection of functions to generate the validations automatically.