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

APPENDIX B. PROGRAM LISTINGS 190
/*
Acoustics.h
By Paul Dickens, 2005
Listing B.1: Acoustics.h
Acoustics library.
*/
#ifndef ACOUSTICS_H_PROTECTOR
#define ACOUSTICS_H_PROTECTOR
#include "Complex.h"
#include "TransferMatrix.h"
double saturationVapourPressureWater(double T);
/*
Calculates the saturation vapour pressure of water.
Parameters:
T: the temperature in Kelvin.
Returns:
The vapour pressure in Pa.
*/
double speedSound(double t, double p, double h, double x_c);
/*
Calculates the speed of sound.
Based on Owen Cramer (1993) J. Acoust. Soc. Am. 93(5) p2510-2616;
formula at p. 2514
Parameters:
t: the temperature in Celsius (between 0 and 30 degrees C).
h: the relative humidity expressed as a fraction (between 0 and
1).
x_c: the mole fraction of carbon dioxide (between 0 and 0.01).
(note that the parameter ranges are not enforced)
(note that Cramer gives x_c = 0.000314 for air)
Returns:
The speed of sound in m/s.
*/
double speedSoundCramer(double t, double p, double x_w, double x_c);
/*
Calculates the speed of sound.
Based on Owen Cramer (1993) J. Acoust. Soc. Am. 93(5) p2510-2616;
formula at p2514.
Parameters:
t: the temperature in Celsius (between 0 and 30 degrees C).
p: the pressure in Pa (between 75 and 102 kPa).
x_w: the mole fraction of water vapour (between 0 and 0.06).
x_c: the mole fraction of carbon dioxide (between 0 and 0.01).
(note that the parameter ranges are not enforced)
Returns:
The speed of sound in m/s.
*/
double densityAirIdeal(double t, double p, double h, double x_c);