Flute acoustics: measurement, modelling and design - School of ...

APPENDIX B. PROGRAM LISTINGS 191
/*
Calculates the density of air using the ideal gas law.
For a more accurate treatment, see Giacomo (1982), Metrologia 18,
33--40.
Parameters:
t: the temperature in Celsius.
p: the pressure in Pa.
h: the relative humidity (between 0 and 1).
x_c: the mole fraction of carbon dioxide.
Returns:
The density of air in kg/m3.
*/
double densityAirGiacomo(double t, double p, double h, double x_c);
/*
Calculates the density of air.
Based on Giacomo (1982), Metrologia 18, 33--40.
Parameters:
t: the temperature in Celsius.
p: the pressure in Pa.
h: the relative humidity (between 0 and 1).
x_c: the mole fraction of carbon dioxide.
Returns:
The density of air in kg/m3.
*/
double compressibilityAir(double t, double p, double x_w);
/*
Calculates the compressibility of air.
Parameters:
t: the temperature in Celsius.
p: the pressure in Pa.
h: the relative humidity (between 0 and 1).
x_w: the mole fraction of water vapour.
Returns:
The compressibility of air.
*/
complex waveNum(double f, double c, double a, double alphacorrection);
/*
Calculates k, the complex wave number for a given segment which
includes wall losses. Refer to Fletcher and Rossing (1998).
Parameters:
f: the frequency.
c: the speed of sound
a: the radius of the tube in metres.
alphacorrection: the multiplicative attenuation coefficient
factor.
Returns:
The complex wave number k, including wall losses.
*/
double phaseVel(double f, double c, double a);
/*
Calculates the phase velocity of a given pipe, part of the wave
number.
Refer to Fletcher and Rossing (1998).
Parameters: