Basic Usage¶

compflow is available on the Python Package Index, so installation is as simple as,

python3 -m pip install compflow

Important

The setup script uses NumPy distutils, so if you do not already have NumPy installed, you should first run python3 -m pip install numpy;
As the library is based on Fortran subroutines, you will need a working Fortran compiler for the installation to complete successfully.

We can now start doing some calculations. First, an explicit evaluation of stagnation pressure ratio \(p_0/p\) given a Mach number \(\Ma\),

>>> import compflow
>>> ga = 1.4
>>> compflow.Po_P_from_Ma(0.3, ga)
1.0644302861529382

Second, an inversion of flow function \(\dot{m}\sqrt{c_pT_0}/Ap_0\) where iterative solution for \(\Ma\) is required,

>>> compflow.Ma_from_mcpTo_APo(0.8, ga)
0.39659360325173604

The names and symbols of non-dimensional quantities are fairly self-explanatory, but a full list is given in the Nomenclature. All functions and the equations used for the calculations are documented in the API reference.

Numpy arrays are also accepted as inputs,

>>> import numpy
>>> Ma1 = numpy.array([0., 0.5, 1., 2.])
>>> compflow.To_T_from_Ma(Ma1, ga)
array([1.  , 1.05, 1.2 , 1.8 ])
>>> Ma2 = numpy.array([[0.1, 0.2], [0.3, 0.4], [0.5, 0.6]])
>>> compflow.To_T_from_Ma(Ma2, ga)
array([[1.002, 1.008],
       [1.018, 1.032],
       [1.05 , 1.072]])

When solving for Mach number at a given normalised mass flow, it is assumed that we are on the subsonic branch of the curve unless a flag is specified. Where no solution is possible, i.e. if the flow would choke, NaN is returned,

>>> capacity = [0.6, 2.]
>>> compflow.Ma_from_mcpTo_APo(capacity, ga)
array([0.28442265,        nan])
>>> compflow.Ma_from_mcpTo_APo(capacity, ga, sup=True)
array([2.27028708,        nan])