WakeTypes.h

#pragma once
 
#include <limits>
#include <vector>
 
namespace marenv::wake
{
struct WakeSample
{
    double deficit; 
    double sourceVelocity[3]; 
};
 
inline constexpr int kUntaggedSource = -1;
 
struct ElementarySource
{
    double pos[3]; 
    double flowDir[3]; 
    double area; 
    double momentumDeficit; 
    double solidity = 0.0; 
    int tag         = kUntaggedSource; 
};
 
struct LineSource
{
    double pos[3]; 
    double flowDir[3]; 
    double axis[3]; 
    double length; 
    double width; 
    double momentumDeficit; 
};
 
struct WakeFrame
{
    double origin[3]; 
    double axes[3][3]; 
    double velocity[3]; 
    double t0; 
};
 
struct GridSpec
{
    double lengthDownstream; 
    double halfWidth; 
    double halfHeight; 
    double upstreamMargin; 
    int n[3]; 
};
 
enum class Superposition
{
    Product, 
    MomentumFlux 
};
 
struct KernelParams
{
    double expansion            = 0.05; 
    double cutoff               = 1e-3; 
    double xStartFactor         = 0.5; 
    double xStartMin            = 0.02; 
    double rampFactor           = 0.0; 
    double farWakeStartFactor   = 10.0; 
    bool surfaceImage           = true; 
    double seabedDepth          = std::numeric_limits<double>::infinity(); 
    Superposition superposition = Superposition::Product; 
};
 
struct SourceSet
{
    std::vector<ElementarySource> points; 
    std::vector<LineSource> lines; 
};
 
struct GriddedSnapshot;
 
class SourceSetWake;
 
class GriddedWake;
struct WakeComposition
{
    double deficit; 
    double sourceVelocity[3]; 
    int numSources; 
    double jetVelocity[3] = {}; 
    int numJets           = 0; 
 
    void Apply(double vel[3]) const;
};
} // namespace marenv::wake