LinearConcentrationDependentDensity.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 LinearConcentrationDependentDensity final : public FluidProperty
{
public:
    explicit LinearConcentrationDependentDensity(
        const double reference_density, double reference_concentration,
        const double fluid_density_difference_ratio)
        : _reference_density(reference_density),
          _reference_concentration(reference_concentration),
          _fluid_density_difference_ratio(fluid_density_difference_ratio)
    {
    }
 
    std::string getName() const override
    {
        return "Linear concentration dependent density";
    }
 
    double getValue(const ArrayType& var_vals) const override
    {
        const double C = var_vals[static_cast<int>(PropertyVariableType::C)];
        return _reference_density * (1 +
                                     _fluid_density_difference_ratio *
                                         (C - _reference_concentration));
    }
 
    double getdValue(const ArrayType& var_vals,
                     const PropertyVariableType var) const override
    {
        (void)var_vals;
        if (var != PropertyVariableType::C)
        {
            return 0.0;
        }
        return _reference_density * _fluid_density_difference_ratio;
    }
 
private:
    const double _reference_density;
    const double _reference_concentration;
    const double _fluid_density_difference_ratio;
};
 
}  // namespace Fluid
}  // namespace MaterialLib