pipeline.hpp

#pragma once
 
#include <armadillo>
 
#include <vector>
 
#include "composition.hpp"
#include "heattransfer/pipewall.hpp"
#include "heattransfer/burialmedium.hpp"
#include "heattransfer/ambientfluid.hpp"
#include "heattransfer/heattransferbase.hpp"
#include "advection/batchtrackingstate.hpp"
 
class BoundaryConditions;
 
class Pipeline
{
public:
    class State
    {
    public:
        friend class Pipeline;
 
        State(
                const arma::vec& gridPoints,
                const arma::vec& pressure,
                const arma::vec& temperature,
                const arma::vec& flow,
                const std::vector<Composition>& composition);
 
        State(
                const arma::vec& gridPoints,
                const arma::vec& pressure,
                const arma::vec& temperature,
                const arma::vec& flow,
                const Composition& composition = Composition::defaultComposition);
 
        // getters
        // main/governing properties
        const arma::vec& flow() const                           { return m_flow; }
        const arma::vec& pressure() const                       { return m_pressure; }
        const arma::vec& temperature() const                    { return m_temperature; }
        const std::vector<Composition>& composition() const     { return m_composition; }
 
        // derived properties
        const arma::vec& heatCapacityConstantVolume() const     { return m_heatCapacityConstantVolume; }
        const arma::vec& heatCapacityConstantPressure() const   { return m_heatCapacityConstantPressure; }
        const arma::vec& density() const                        { return m_density; }
        const arma::vec& viscosity() const                      { return m_viscosity; }
        const arma::vec& specificGasConstant() const            { return m_specificGasConstant; }
        const arma::vec& molarMass() const                      { return m_molarMass; }
 
        const arma::vec& compressibilityFactor() const          { return m_compressibilityFactor; }
        const arma::vec& dZdtAtConstantPressure() const         { return m_temperatureDerivativeConstantPressure; }
        const arma::vec& dZdpAtConstantTemperature() const      { return m_pressureDerivativeConstantTemperature; }
        const arma::vec& dZdtAtConstantDensity() const          { return m_temperatureDerivativeConstantDensity; }
 
        const arma::vec& velocity() const                       { return m_velocity; }
        const arma::vec& frictionFactor() const                 { return m_frictionFactor; }
        const arma::vec& reynoldsNumber() const                 { return m_reynoldsNumber; }
 
        const arma::vec& ambientTemperature() const             { return m_ambientTemperature; }
        const arma::vec& heatFlow() const                       { return m_heatFlow; }
 
//        const std::vector<HeatTransferState>& heatTransferState() const { return m_heatTransferState; }
//        const BatchTrackingState& batchTrackingState() const  { return m_batchTrackingState; }
//        bool batchTrackingIsInitialized() const                 { return m_batchTrackingIsInitialized; }
 
//        // setters
//        // main/governing properties
//        arma::vec& flow()                           { return m_flow; }
//        arma::vec& pressure()                       { return m_pressure; }
//        arma::vec& temperature()                    { return m_temperature; }
//        std::vector<Composition>& composition();
 
//        // derived properties
//        arma::vec& heatCapacityConstantVolume()     { return m_heatCapacityConstantVolume; }
//        arma::vec& heatCapacityConstantPressure()   { return m_heatCapacityConstantPressure; }
//        arma::vec& density()                        { return m_density; }
//        arma::vec& viscosity()                      { return m_viscosity; }
//        arma::vec& specificGasConstant()            { return m_specificGasConstant; }
//        arma::vec& molarMass()                      { return m_molarMass; }
 
//        arma::vec& compressibilityFactor()          { return m_compressibilityFactor; }
//        arma::vec& dZdtAtConstantPressure()         { return m_temperatureDerivativeConstantPressure; }
//        arma::vec& dZdpAtConstantTemperature()      { return m_pressureDerivativeConstantTemperature; }
//        arma::vec& dZdtAtConstantDensity()          { return m_temperatureDerivativeConstantDensity; }
 
//        arma::vec& velocity()                       { return m_velocity; }
//        arma::vec& frictionFactor()                 { return m_frictionFactor; }
//        arma::vec& reynoldsNumber()                 { return m_reynoldsNumber; }
 
//        arma::vec& ambientTemperature()             { return m_ambientTemperature; }
//        arma::vec& heatFlow()                       { return m_heatFlow; }
 
//        std::vector<HeatTransferState>& heatTransferState() { return m_heatTransferState; }
//        BatchTrackingState& batchTrackingState()  { return m_batchTrackingState; }
//        bool& batchTrackingIsInitialized()          { return m_batchTrackingIsInitialized; }
 
        friend std::ostream& operator << (std::ostream &out, const State& c);
 
    protected:
        State(const arma::vec& gridPoints);
 
        // input/governing properties
        arma::vec m_flow; 
        arma::vec m_pressure; 
        arma::vec m_temperature; 
        std::vector<Composition> m_composition; 
 
        // output/derived properties
        arma::vec m_heatCapacityConstantVolume; 
        arma::vec m_heatCapacityConstantPressure; 
        arma::vec m_density; 
        arma::vec m_viscosity; 
        arma::vec m_specificGasConstant; 
        arma::vec m_molarMass; 
 
        arma::vec m_compressibilityFactor; 
        arma::vec m_temperatureDerivativeConstantPressure;
        arma::vec m_pressureDerivativeConstantTemperature;
        arma::vec m_temperatureDerivativeConstantDensity;
 
        arma::vec m_velocity; 
        arma::vec m_frictionFactor; 
        arma::vec m_reynoldsNumber; 
 
        arma::vec m_ambientTemperature; 
        arma::vec m_heatFlow; 
 
        std::vector<HeatTransferState> m_heatTransferState; 
        bool m_heatTransferIsInitialized; 
 
        BatchTrackingState m_batchTrackingState; 
        bool m_batchTrackingIsInitialized; 
    };
 
    Pipeline(const arma::uword nGridPoints = 100, const double length = 100e3);
 
    arma::uword size() const                                { return m_gridPoints.n_elem; }
 
    BoundaryConditions getBoundaryConditions() const;
 
    void enableBatchTracking();
 
    void initializeBatchTracking();
 
    void setLength(const double length);
 
    arma::uword timestamp() const { return m_timestamp; }
    arma::uword& timestamp() { return m_timestamp; }
 
    const Pipeline::State& state() const { return m_state; }
 
    // getters
    // constants
    double length() const                                   { return m_length; }
    const arma::vec& gridPoints() const                     { return m_gridPoints; }
    const arma::vec& diameter() const                       { return m_diameter; }
    const arma::vec& height() const                         { return m_height; }
    const arma::vec& elevation() const                      { return m_height; }
    const arma::vec& roughness() const                      { return m_roughness; }
 
    const arma::vec& burialDepth() const                    { return m_burialDepth; }
    const std::vector<PipeWall>& pipeWall() const           { return m_pipeWall; }
    const std::vector<BurialMedium>& burialMedium() const   { return m_burialMedium; }
    const std::vector<AmbientFluid>& ambientFluid() const   { return m_ambientFluid; }
 
    // main/governing properties
    const arma::vec& flow() const                           { return prop().m_flow; }
    const arma::vec& pressure() const                       { return prop().m_pressure; }
    const arma::vec& temperature() const                    { return prop().m_temperature; }
    const std::vector<Composition>& composition() const     { return prop().m_composition; }
 
    // derived properties
    const arma::vec& heatCapacityConstantVolume() const     { return prop().m_heatCapacityConstantVolume; }
    const arma::vec& heatCapacityConstantPressure() const   { return prop().m_heatCapacityConstantPressure; }
    const arma::vec& density() const                        { return prop().m_density; }
    const arma::vec& viscosity() const                      { return prop().m_viscosity; }
    const arma::vec& specificGasConstant() const            { return prop().m_specificGasConstant; }
    const arma::vec& molarMass() const                      { return prop().m_molarMass; }
 
    const arma::vec& compressibilityFactor() const          { return prop().m_compressibilityFactor; }
    const arma::vec& dZdtAtConstantPressure() const         { return prop().m_temperatureDerivativeConstantPressure; }
    const arma::vec& dZdpAtConstantTemperature() const      { return prop().m_pressureDerivativeConstantTemperature; }
    const arma::vec& dZdtAtConstantDensity() const          { return prop().m_temperatureDerivativeConstantDensity; }
 
    const arma::vec& velocity() const                       { return prop().m_velocity; }
    const arma::vec& frictionFactor() const                 { return prop().m_frictionFactor; }
    const arma::vec& reynoldsNumber() const                 { return prop().m_reynoldsNumber; }
 
    const arma::vec& ambientTemperature() const             { return prop().m_ambientTemperature; }
    const arma::vec& heatFlow() const                       { return prop().m_heatFlow; }
 
    const std::vector<HeatTransferState>& heatTransferState() const { return prop().m_heatTransferState; }
    bool heatTransferIsInitialized() const                  { return prop().m_heatTransferIsInitialized; }
 
    const BatchTrackingState& batchTrackingState() const    { return prop().m_batchTrackingState; }
    bool batchTrackingIsInitialized() const                 { return prop().m_batchTrackingIsInitialized; }
 
    // settings
    bool constantComposition() const                        { return m_constantComposition; }
 
    // helpers
    const arma::vec& inletComposition() const;
    const arma::vec& outletComposition() const;
 
    // setters
    // constants
    arma::vec& gridPoints()                     { return m_gridPoints; }
    arma::vec& diameter()                       { return m_diameter; }
    arma::vec& height()                         { return m_height; }
    arma::vec& elevation()                      { return m_height; }
    arma::vec& roughness()                      { return m_roughness; }
 
    arma::vec& burialDepth()                    { return m_burialDepth; }
    std::vector<PipeWall>& pipeWall()           { return m_pipeWall; }
    std::vector<BurialMedium>& burialMedium()   { return m_burialMedium; }
    std::vector<AmbientFluid>& ambientFluid()   { return m_ambientFluid; }
 
    // main/governing properties
    arma::vec& flow()                           { return prop().m_flow; }
    arma::vec& pressure()                       { return prop().m_pressure; }
    arma::vec& temperature()                    { return prop().m_temperature; }
 
    // derived properties
    arma::vec& heatCapacityConstantVolume()     { return prop().m_heatCapacityConstantVolume; }
    arma::vec& heatCapacityConstantPressure()   { return prop().m_heatCapacityConstantPressure; }
    arma::vec& density()                        { return prop().m_density; }
    arma::vec& viscosity()                      { return prop().m_viscosity; }
    arma::vec& specificGasConstant()            { return prop().m_specificGasConstant; }
    arma::vec& molarMass()                      { return prop().m_molarMass; }
 
    arma::vec& compressibilityFactor()          { return prop().m_compressibilityFactor; }
    arma::vec& dZdtAtConstantPressure()         { return prop().m_temperatureDerivativeConstantPressure; }
    arma::vec& dZdpAtConstantTemperature()      { return prop().m_pressureDerivativeConstantTemperature; }
    arma::vec& dZdtAtConstantDensity()          { return prop().m_temperatureDerivativeConstantDensity; }
 
    arma::vec& velocity()                       { return prop().m_velocity; }
    arma::vec& frictionFactor()                 { return prop().m_frictionFactor; }
    arma::vec& reynoldsNumber()                 { return prop().m_reynoldsNumber; }
 
    arma::vec& ambientTemperature()             { return prop().m_ambientTemperature; }
    arma::vec& heatFlow()                       { return prop().m_heatFlow; }
 
    std::vector<HeatTransferState>& heatTransferState() { return prop().m_heatTransferState; }
    bool& heatTransferIsInitialized()           { return prop().m_heatTransferIsInitialized; }
 
    BatchTrackingState& batchTrackingState()    { return prop().m_batchTrackingState; }
    bool& batchTrackingIsInitialized()          { return prop().m_batchTrackingIsInitialized; }
 
    // settings
    bool& constantComposition()                 { return m_constantComposition; }
 
    void updateComposition(const std::vector<Composition>& composition);
 
    void updateComposition(const Composition& composition);
 
    void setCompositionUnsafe(const std::vector<Composition>& composition);
 
private:
    // constants
    double m_length; 
    arma::vec m_gridPoints; 
    arma::vec m_diameter; 
    arma::vec m_height; 
    arma::vec m_roughness; 
 
    arma::vec m_burialDepth; 
    std::vector<PipeWall> m_pipeWall; 
    std::vector<BurialMedium> m_burialMedium; 
    std::vector<AmbientFluid> m_ambientFluid; 
 
    // settings
    bool m_constantComposition; 
 
    // other
    State m_state; 
    arma::uword m_timestamp = 0; 
 
    // shorthand for internal use
    const State& prop() const { return m_state; }
    State& prop() { return m_state; }
};