LinearConcentrationAndPressureDependentDensity.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 LinearConcentrationAndPressureDependentDensity final
    : public FluidProperty
{
public:
    explicit LinearConcentrationAndPressureDependentDensity(
        const double reference_density, double reference_concentration,
        const double fluid_density_concentration_difference_ratio,
        double reference_pressure,
        const double fluid_density_pressure_difference_ratio)
        : _reference_density(reference_density),
          _reference_concentration(reference_concentration),
          _fluid_density_concentration_difference_ratio(
              fluid_density_concentration_difference_ratio),
          _reference_pressure(reference_pressure),
          _fluid_density_pressure_difference_ratio(
              fluid_density_pressure_difference_ratio)
    {
    }
 
    std::string getName() const override
    {
        return "Linear concentration and pressure dependent density";
    }
 
    double getValue(const ArrayType& var_vals) const override
    {
        const double C = var_vals[static_cast<int>(PropertyVariableType::C)];
        const double p = var_vals[static_cast<int>(PropertyVariableType::p)];
        return _reference_density *
               (1 +
                _fluid_density_concentration_difference_ratio *
                    (C - _reference_concentration) +
                _fluid_density_pressure_difference_ratio *
                    (p - _reference_pressure));
    }
 
    double getdValue(const ArrayType& /*var_vals*/,
                     const PropertyVariableType var) const override
    {
        if (var == PropertyVariableType::C)
        {
            return _reference_density *
                   _fluid_density_concentration_difference_ratio;
        }
        if (var == PropertyVariableType::p)
        {
            return _reference_density *
                   _fluid_density_pressure_difference_ratio;
        }
        return 0;
    }
 
private:
    const double _reference_density;
    const double _reference_concentration;
    const double _fluid_density_concentration_difference_ratio;
    const double _reference_pressure;
    const double _fluid_density_pressure_difference_ratio;
};
 
}  // namespace Fluid
}  // namespace MaterialLib