turbine.hpp

Go to the documentation of this file.

#pragma once
 
#include <span>
#include <vector>
 
#include "interfaces/components/beam.hpp"
#include "interfaces/components/turbine_input.hpp"
#include "interfaces/constraint_data.hpp"
#include "interfaces/host_state.hpp"
#include "interfaces/node_data.hpp"
 
namespace kynema {
class Model;
}
 
namespace kynema::interfaces::components {
class Turbine {
public:
    //--------------------------------------------------------------------------
    // Types and Constants
    //--------------------------------------------------------------------------
 
    using DeviceType =
        Kokkos::Device<Kokkos::DefaultExecutionSpace, Kokkos::DefaultExecutionSpace::memory_space>;
 
    static constexpr size_t invalid_id{9999999};
 
    static constexpr double kMinHubDiameter{1e-8};
 
    static constexpr double kZeroTolerance{1e-12};
 
    //--------------------------------------------------------------------------
    // Elements
    //--------------------------------------------------------------------------
 
    std::vector<Beam> blades;                        //< Blades in the turbine
    Beam tower;                                      //< Tower in the turbine
    size_t hub_mass_element_id{invalid_id};          //< Hub mass element ID
    size_t nacelle_cm_mass_element_id{invalid_id};   //< Nacelle center of mass mass element ID
    size_t yaw_bearing_mass_element_id{invalid_id};  //< Yaw bearing mass element ID
 
    //--------------------------------------------------------------------------
    // Nodes
    //--------------------------------------------------------------------------
 
    std::vector<NodeData> apex_nodes;  //< Blade root nodes
    NodeData hub_node;                 //< Hub node
    NodeData azimuth_node;             //< Azimuth node
    NodeData shaft_base_node;          //< Shaft base node
    NodeData nacelle_cm_node;          //< Nacelle center of mass node
    NodeData yaw_bearing_node;         //< Yaw bearing node
 
    //--------------------------------------------------------------------------
    // Constraints
    //--------------------------------------------------------------------------
 
    ConstraintData tower_base;                 //< Tower base constraint
    ConstraintData tower_top_to_yaw_bearing;   //< Tower top to yaw bearing constraint
    ConstraintData yaw_bearing_to_shaft_base;  //< Yaw bearing to shaft base constraint
    ConstraintData yaw_bearing_to_nacelle_cm;  //< Yaw bearing to nacelle center of mass constraint
    ConstraintData shaft_base_to_azimuth;      //< Nacelle mass to generator constraint
    ConstraintData azimuth_to_hub;             //< Azimuth to hub constraint
    std::vector<ConstraintData> blade_pitch;   //< Blade root to apex constraints
    std::vector<ConstraintData> apex_to_hub;   //< Apex to hub constraints
 
    //--------------------------------------------------------------------------
    // Control inputs
    //--------------------------------------------------------------------------
 
    std::vector<double> blade_pitch_control;  //< Blade pitch angles
    double torque_control{0.};                //< Torque control value
    double yaw_control{0.};                   //< Yaw control value
 
    Turbine(const TurbineInput& input, Model& model);
 
    Turbine(Turbine& other) = delete;
 
    Turbine(Turbine&& other) = delete;
 
    Turbine& operator=(const Turbine&) = delete;
 
    Turbine& operator=(Turbine&&) = delete;
 
    ~Turbine() = default;
 
    void GetMotion(const HostState<DeviceType>& host_state);
 
    void GetLoads(const HostConstraints<DeviceType>& host_constraints);
 
    void SetLoads(HostState<DeviceType>& host_state) const;
 
    [[nodiscard]] const TurbineInput& GetTurbineInput() const;
 
private:
    //--------------------------------------------------------------------------
    // Node ID collections of turbine components
    //--------------------------------------------------------------------------
 
    std::vector<size_t> tower_node_ids;       
    std::vector<size_t> drivetrain_node_ids;  
    std::vector<size_t> blade_node_ids;       
    std::vector<size_t> rotor_node_ids;       
    std::vector<size_t> nacelle_node_ids;     
    std::vector<size_t>
        all_turbine_node_ids;  
 
    //--------------------------------------------------------------------------
    // Turbine inputs
    //--------------------------------------------------------------------------
 
    TurbineInput turbine_input;  
 
    [[nodiscard]] static std::vector<Beam> CreateBlades(
        std::span<const BeamInput> blade_inputs, Model& model
    );
 
    static void ValidateInput(const TurbineInput& input);
 
    void PositionNodes(const TurbineInput& input, Model& model);
 
    void CreateIntermediateNodes(const TurbineInput& input, Model& model);
 
    void AddMassElements(const TurbineInput& input, Model& model);
    void AddConstraints(const TurbineInput& input, Model& model);
 
    void SetInitialConditions(const TurbineInput& input, Model& model);
    void SetInitialDisplacements(const TurbineInput& input, Model& model);
 
    void SetInitialRotorVelocity(const TurbineInput& input, Model& model);
};
 
}  // namespace kynema::interfaces::components