PanelFormulas.h

#pragma once
 
#include "PanelLoadTypes.h"
 
#include <cmath>
 
namespace hydrodynamics
{
namespace panel_formulas
{
 
template<class T>
T CylinderDragCoefficient(const T& re)
{
   using std::log10;
 
   const double coefficients[8] = {-78.46675, 254.73873, -327.8864, 223.64577, -87.92234,
         20.00769, -2.44894, 0.12479};
 
   const T x  = log10(re);
   T       cd = T(coefficients[7]);
   for (int i = 6; i >= 0; i--) {
      cd = cd * x + coefficients[i];
   }
   return cd;
}
 
 
template<class T>
T VelocityReductionFactor(const T& solidity, unsigned& flags)
{
   const T ratio = 1.08 - 0.97 * solidity;
   if (ratio > 1.0) {
      flags |= kPanelFlagVelocityRatioClamped;
      return T(1.0);
   }
   return ratio;
}
 
 
template<class T>
T VelocityReductionFactor(const T& solidity)
{
   unsigned ignoredFlags = 0u;
   return VelocityReductionFactor(solidity, ignoredFlags);
}
 
} // namespace panel_formulas
} // namespace hydrodynamics