LineKernel.h

#pragma once
 
#include <marenv/wake/Kernel.h>
#include <marenv/wake/WakeTypes.h>
 
namespace marenv::wake
{
struct LineCoordinates
{
    double x; 
    double s; 
    double y; 
};
 
class LineKernelProfile
{
public:
    LineKernelProfile(const LineSource& source, const KernelParams& params);
 
    double XStart() const { return m_xStart; }
 
    const double* Position() const { return m_pos; }
 
    LineCoordinates Coordinates(const double pos[3]) const;
 
    KernelCentreline At(double x) const;
 
    double Amplitude(double sigma, double envelope) const;
 
    double Ramp(double x) const;
 
    double PeakAmplitude(double xFrom, double xTo) const;
 
    double HalfSpan() const { return m_halfSpan; }
 
private:
    double LinearSigma(double x) const;
 
    double Sigma(double x) const;
 
    double m_pos[3]; 
    double m_flow[3]; 
    double m_span[3]; 
    double m_cross[3]; 
    double m_halfSpan; 
    double m_width; 
    double m_expansion; 
    double m_xStart; 
    double m_rampLength; 
    double m_momentumPerLength; 
    double m_nearWakeDeficit; 
    double m_epsilon; 
    double m_farWakeStart; 
    double m_farWakeSigma; 
};
 
double LineDeficit(const LineKernelProfile& profile, const double pos[3]);
 
double LineDeficit(const LineSource& source, const KernelParams& params, const double pos[3]);
} // namespace marenv::wake