CNetStructureWakeEffect.h

#pragma once
 
//#define CNetStructureWakeEffect_DO_PROFILE
 
// Includes
#include <string>
 
#include <CEnvironment.h>
 
#include <SimObject.h>
#include <ICurrentVel.h>
#include "sfh/timers/ProfilerWallclock.h"
#include <IBasicBottom.h>
 
#include <environment/DeepSeaGravityWaves.h>
#include "CNetElement3NWakeEffect.h"
#include "CNetCableElementSpec.h"
#include "CNetCableElement.h"
#include "CNetElementDef.h"
#include "CNetBottomContactFilter.h"
#include "CNetStructureSpec.h"
 
 
#ifdef FH_VISUALIZATION
    #include "sfh/ogre/NMeshGenerators.h"
    #include "sfh/ogre/C3DLine.h"
    #include <CTextRenderer.h>
#endif
 
//Class definition
class CNetStructureWakeEffect : public SimObject
{
public:
    CNetStructureWakeEffect(std::string sSimObjectName, ISimObjectCreator *pCreator);
    virtual ~CNetStructureWakeEffect();
    virtual void OdeFcn(const double dT, const double *const adX, double *const adXDot, const bool bIsMajorTimeStep);
    virtual void InitialConditionSetup(const double dT, const double *const adCurrentIC, double *const adUpdatedIC, ISimObjectCreator *const pCreator);
    virtual void FinalSetup(const double dT, const double *const adX, ISimObjectCreator *const pCreator);
    const double * outExternalNodeForce(const double dT, const double *const adX, int iExtNode);
 
#ifdef FH_VISUALIZATION
    virtual void RenderInit(Ogre::Root *const pOgreRoot, ISimObjectCreator *const pCreator);
    virtual void RenderUpdate(const double dT, const double *const adX);
 
    void CreateCapsule( Ogre::ManualObject *capsule_manual);
    C3DLine *myline;
#endif
    virtual const double * NodePos(const double dT, const double *const adX, int iExtNode);
    virtual const double * NodeVel(const double dT, const double *const adX, int iExtNode);
    virtual const double * BottomContact(const double dT, const double *const adX);
    void SetAddedDragPerMPS(double AddedDragPerMPS);
    const double* GetNetTotalForceForCheck(const double dT, const double *const adX); 
protected:
    virtual void CalcStateDerivatives(const double *const adX, double *const adXDot);
    virtual void AddExternalForces(const double dT, const double *const adX);
    virtual void AddNetForces(const double dT, const double *const adX);
    virtual void AddCableForces(const double dT, const double *const adX);
    virtual double AddNodeBottomForce(const double *const adX, int iNode, double dDiam, double adForceInOut[3]);
    virtual void AddBottomForces(const double *const adX);
    virtual void AddInputPorts(ISimObjectCreator *pCreator);
    virtual void CreateNetPanels(std::vector<const CNetElementDef*> VNetElementDefs, ISimObjectCreator *pCreator);
    virtual void CreateNetCables(std::vector<const CNetCableElementSpec*> VCableElements, ISimObjectCreator *pCreator);
    virtual void CreateDataStructures();
    virtual void CreateOutputPorts(ISimObjectCreator *pCreator);
    VSPairStrInt CreateExternalNodeMap( std::string sTag, int* piUsed, CNetStructureSpec & NetSpec, ISimObjectCreator *pCreator);
    void PRINTALLNODEFORCES(std::string spec, double dT);
 
 
    ISignalPort **m_apInExternalForces;
    int *m_aIStatePos;
    int *m_aIStateVel;
    double*         m_NodeExternalForces;
 
    int m_iNumNodesInputForce;
    int m_iNumNodesOutputPosAndVel;
    int m_iNumPanelsCalc;
    int m_iNumCablesCalc;
    int m_iNumNodesOriginal;
    int m_iNumNodesCalc;
    int m_iNumNodeOutputs;
 
 
    CNetBottomContactFilter m_BottomContactFilter;
    double                  m_adBottomContact[6];
 
    IBasicBottom**          m_ppBottom;
    ICurrentVel**           m_ppCurrentVel;
    CEnvironment*           m_pEnvironment;
 
 
    std::vector<CNetCableElement*>  m_VpCableElements;
    std::vector<CNetElement3NWakeEffect*>   m_VpNetElements;
    float*                          m_afInertia;
    int*                            m_aiNetConn;
    int**                           m_aaiCableConn;
    SAdaptationSpec_t               m_SAdaptations;
    bool                            m_bFoldingPointRefinementUsed;
 
 
 
    VSPairStrInt            m_VsNodesOutputPosAndVel;
    VSPairStrInt            m_VsNodesInputForce;
    bool                    m_bUseTagForExternalNodes;
 
    double*                 m_adNodeForce;
    sfh::timers::StopWatch m_OdeTimer;
 
    // Screen model
    bool m_UseScreenModel;
    bool m_UseMoeFore2021; // Use drag and lift coefficients from model tests published in Moe-F�re et al. 2021.
    double m_SnImg; // Use solidity derived from image analysis. Only to be used with drag formulation from Moe-F�re et al. 2021.
    double m_RhoWater;
    double m_NetForcesTotalForOutputCheck[3];
 
    // Wake effect inclusion
    bool                    m_ExternalElementWakeEffectsActive;                     
    bool                    m_InternalElementWakeEffectsActive;                     
    bool                    m_UseActualElementSizeInWakeEffect;                     
    int                     m_NElementBarsInWakeEffectCalc;                         
    bool                    m_useReynoldsDependentDragCoefficients;
 
#ifdef CNetStructureWakeEffect_DO_PROFILE
    sfh::timers::ProfilerWallclock m_Profiler;
#endif
 
 
};