propeller.cc

#include <vtestbed/base/blitz.hh>
#include <vtestbed/base/constants.hh>
#include <vtestbed/config/real_type.hh>
#include <vtestbed/core/propeller.hh>

template class vtb::core::Propeller<VTB_REAL_TYPE>;

template <class T>
T vtb::core::Propeller<T>::thrust() const {
    using blitz::pow2;
    using blitz::pow4;
    using vtb::base::constants;
    return constants<T>::water_density() * pow2(revolutions()) *
        pow4(diameter()) * thrust_coefficient();
}

template <class T>
T vtb::core::Propeller<T>::torque() const {
    using blitz::pow2;
    using blitz::pow5;
    using vtb::base::constants;
    return constants<T>::water_density() * pow2(revolutions()) *
        pow5(diameter()) * torque_coefficient();
}

template <class T>
T vtb::core::Propeller<T>::advance_coefficient(T velocity) const {
    return velocity / (revolutions() * diameter());
}

template <class T>
T vtb::core::Propeller<T>::efficiency(T velocity) const {
    using vtb::base::constants;
    return advance_coefficient(velocity) * thrust_coefficient() /
        (constants<T>::pi(2)*torque_coefficient());
}