Skip to content

Latest commit

 

History

History
160 lines (119 loc) · 6.08 KB

README.md

File metadata and controls

160 lines (119 loc) · 6.08 KB
Raw

Speed of Sound

CircleCI Coverage Status

Speed of Sound is a library to efficiently compute the speed of sound in air. It is intended for use in embedded systems that require a high-precision computation of the speed of sound in air based a variety of ambient conditions.

Contents

Quick facts

  • Computes the speed of sound in air based on temperature, relative humidity, pressure, and carbon dioxide mole fraction.
  • Runtime of subsequent speed of sound computations reduced by 67.7% through the use of calculus linear approximations (reference).
  • Approximation maintains precision of at least 0.05% with environmental factors varying up to 20% (reference).

Environmental parameters

The speed of sound computation takes on the following parameters and is only valid in within the following range.

Variable (units) Minimum Default Maximum
temperature (°C) 0.0 20.0 30.0
relative humidity (dimensionless) 0.0 0.5 1.0
pressure (Pa) 75000.0 101325.0 102000.0
Carbon dioxide mole fraction (dimensionless) 0.0 0.000314 0.01

Note: This library is intended for use with embedded systems and therefore does not throw an exception if passed an invalid quantity. Instead, methods are provided to validate that values are within the given constraints (see Usage).

Usage

#include "speed_of_sound.h"

Create Environment object to hold ambient conditions (initializes values to default).

speedofsound::Environment ambient_conitions;

Access Environment object members.

ambient_conitions.temperature_;
ambient_conitions.humidity_;
ambient_conitions.pressure_;
ambient_conitions.co2_mole_fraction_;

Adjust environment variables.

ambient_conitions.temperature_ = 25.0;
...

Check if environmental parameters are within specified constraints (see Constraints for details).

bool is_valid;

// Validate individual parameters
is_valid = ambient_conitions.ValidateTemperature();
is_valid = ambient_conitions.ValidateHumidity();
is_valid = ambient_conitions.ValidatePressure();
is_valid = ambient_conitions.ValidateCO2MoleFraction();

// Validate all parameters
is_valid = ambient_conitions.ValidateEnvironment();

Create SpeedOfSound object and set initial conditions.

speedofsound::SpeedOfSound speed_of_sound(ambient_conitions);

speedofsound::SpeedOfSound speed_of_sound(); // Default Environment values

Compute speed of sound.

double sound_speed;

// Full computation, reset initial conditions
sound_speed = speed_of_sound.Compute(ambient_conitions);

// Full computation, do not reset initial conditions
sound_speed = speed_of_sound.QuickCompute(ambient_conitions);

// Linear approximation, do not reset initial conditions
sound_speed = speed_of_sound.Approximate(ambient_conitions);

Example

#include <iostream>
#include "speed_of_sound.h"
#include "some_sensor.h"

int main() {
  speedofsound::Environment ambient_conitions;
  ambient_conitions.temperature_ = sensor::MeasureTemperature();

  // Initial "full" computation
  speedofsound::SpeedOfSound speed_of_sound(ambient_conitions);

  double sound_speed;
  const double distance = 100.0; // meters
  while (true) {
    ambient_conitions.temperature_ = sensor::MeasureTemperature();
    if (!ambient_conitions.ValidateTemperature()) break;

    sound_speed = speed_of_sound.Approximate(ambient_conitions);

    std::cout << "It would take a sound wave approximately " << \
                 distance/sound_speed << " seconds to travel " << \
                 distance << " meters." << std::endl;
  }

  return 0;
}

Notes on notation

The following abbreviations are used in theory-related computations.

Lower-case variables (t, h, etc.) represent measured environmental parameters. Variables starting with an uppercase character (C, Psv, etc.) represent computed values.

Variable Meaning Units
C speed of sound meters per second (m/s)
t temperature Celsius (°C)
T thermodynamic temperature Kelvin (K)
h relative humidity fraction (dimensionless)
xc Carbon dioxide mole fraction fraction (dimensionless)
p pressure Pascals (Pa)
Xw mole fraction of water in air fraction (dimensionless)
Psv saturation vapor pressure Pascals (Pa)
F Enhancement factor for moist air dimensionless
k00, k01, ... experimental constants various

dy_dx denotes the partial derivative of y with respect to x.

Testing

Tests are written using GoogleTest.

Run all tests:

$ test/test

Attributions

The equation for computing the speed of sound in air uses Owen Cramer's research.

J. Acoust. Soc. Am. 93, 2510 (1993); http://dx.doi.org/10.1121/1.405827