unsteady.hpp

#pragma once
 
#include <armadillo>
 
#include "heattransfer/radial.hpp"
 
class UnsteadyHeatTransfer : public RadialHeatTransfer
{
public:
    UnsteadyHeatTransfer(
            const double diameter,
            const PipeWall& pipeWall,
            const double burialDepth,
            const BurialMedium& burialMedium,
            const AmbientFluid& ambientFluid);
 
    virtual HeatTransferState evaluate(
            const HeatTransferState& current,
            const double timeStep,
            const double ambientTemperature,
            const double gasPressure,
            const double gasTemperature,
            const double gasReynoldsNumber,
            const double gasHeatCapacity,
            const double gasViscosity) const override;
 
    HeatTransferState evaluateInternal(
            const arma::vec& shellTemperature,
            const double timeStep,
            const double ambientTemperature,
            const double gasPressure,
            const double gasTemperature,
            const double gasReynoldsNumber,
            const double gasHeatCapacity,
            const double gasViscosity) const;
 
    HeatTransferState thermalizeToSteadyState(
            const double ambientTemperature,
            const double gasPressure,
            const double gasTemperature,
            const double gasReynoldsNumber,
            const double gasHeatCapacity,
            const double gasViscosity) const;
 
private:
    arma::vec m_heatTransferCoefficient;
 
    arma::vec solveEquations(
            const arma::vec& shellTemperature,
            const double timeStep,
            const double gasPressure,
            const double gasTemperature,
            const double ambientTemperature,
            const double gasReynoldsNumber,
            const double gasHeatCapacity,
            const double gasViscosity) const;
};