LinearTemperatureDependentDensity.h

Go to the documentation of this file.

 
#pragma once
 
#include <string>
 
#include "BaseLib/ConfigTree-fwd.h"
#include "MaterialLib/Fluid/FluidProperty.h"
 
namespace MaterialLib
{
namespace Fluid
{
class LinearTemperatureDependentDensity final : public FluidProperty
{
public:
    explicit LinearTemperatureDependentDensity(const double rho0, double T0,
                                               const double beta)
        : _rho0(rho0), _temperature0(T0), _beta(beta)
    {
    }
 
    std::string getName() const override
    {
        return "Linear temperature dependent density";
    }
 
    double getValue(const ArrayType& var_vals) const override
    {
        const double T = var_vals[static_cast<int>(PropertyVariableType::T)];
        return _rho0 * (1 - _beta * (T - _temperature0));
    }
 
    double getdValue(const ArrayType& var_vals,
                     const PropertyVariableType var) const override
    {
        (void)var_vals;
        if (var != PropertyVariableType::T)
        {
            return 0.0;
        }
        return -_rho0 * _beta;
    }
 
private:
    const double _rho0;          
    const double _temperature0;  
    const double _beta;          
};
 
}  // namespace Fluid
}  // namespace MaterialLib