JAcousticsSupportkit.hh

Go to the documentation of this file.

#ifndef __JACOUSTICS__JACOUSTICSSUPPORTKIT__
#define __JACOUSTICS__JACOUSTICSSUPPORTKIT__

#include <cmath>

/**
 * \file
 *
 * Acoustics toolkit.
 * \author mdejong
 */
namespace JACOUSTICS {}
namespace JPP { using namespace JACOUSTICS; }

namespace JACOUSTICS {

  static const double  TOAMIN_S                  =    0.0000000;   //!< Minimal allowed time-of-arrival [s]
  static const double  TOAMAX_S                  =    0.6310912;   //!< Maximal allowed time-of-arrival [s]


  /**
   * Auxiliary data structure for calculation of attenuation length.
   *
   * Attenuation accourding to Ainslie and McColm, J. Acoust. Soc. AM. 103 (3) 1671 1998.
   * See also http://resource.npl.co.uk/acoustics/techguides/seaabsorption/physics.html for numerical values.
   */
  struct JAttenuationLength {
    /**
     * Constructor.
     *
     * \param  T                    temperature [C]
     * \param  S                    salinity [pm]
     * \param  pH                   pH
     */
    JAttenuationLength(const double T,
                       const double S,
                       const double pH) :
      T (T),
      S (S),
      pH(pH)
    {}


    /**
     * Get attentuation length.
     *
     * \param  D_m        depth [m]
     * \param  f_kHz      frequency [kHz]
     * \return            attenuation length [m]
     */
    double operator()(const double D_m,
                      const double f_kHz) const
    {
      const double z  = D_m * 1.0e-3;
      const double f1 = 0.78    * sqrt(S/35.0) * exp(T/26.0);
      const double f2 = 42.0    * exp(T/17.0);
      const double a1 = 0.106   * f1*(f_kHz*f_kHz) / (f_kHz*f_kHz + f1*f1) * exp((pH-8.0)/0.56);
      const double a2 = 0.52    * (1.0 + T/43.0) * (S/35.0) * f2 * f_kHz*f_kHz / (f_kHz*f_kHz + f2*f2) * exp(-z/6.0);
      const double a0 = 0.00049 * exp(-T/27.0 + z/17.0) * f_kHz*f_kHz;
      const double A  = a0 + a1 + a2;
      
      return 10.0 / (A*log(10.0));
    }

    double T;
    double S;
    double pH;
  };


  /**
   * Function object to calculate attenutation length.
   */
  static const JAttenuationLength getAttenuationLength(13.2,      // temperature [deg]
                                                       38.0,      // salinity [pm]
                                                       8.0);      // pH
}

#endif