SpecificHeatCapacityWithLatentHeat.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
 
namespace MaterialPropertyLib
{
class SpecificHeatCapacityWithLatentHeat final : public Property
{
    struct PhaseProperties
    {
        Property const* liquid = nullptr;
        Property const* frozen = nullptr;
        Property const* porous = nullptr;
    };
 
public:
    SpecificHeatCapacityWithLatentHeat(std::string name, double const l);
 
    void checkScale() const override;
 
    // initialize container with pointers to the phases properties
    void setProperties(
        std::vector<std::unique_ptr<Phase>> const& phases) override;
 
    double effectiveVolumetricHeatCapacity(
        VariableArray const& variable_array,
        ParameterLib::SpatialPosition const& pos,
        double const t,
        double const dt) const;
 
    PropertyDataType value(VariableArray const& variable_array,
                           ParameterLib::SpatialPosition const& pos,
                           double const t,
                           double const dt) const override;
 
    PropertyDataType dValue(VariableArray const& variable_array,
                            Variable const variable,
                            ParameterLib::SpatialPosition const& pos,
                            double const t,
                            double const dt) const override;
 
private:
    double const l_;  //< specific latent heat of melting

    PhaseProperties densities_;
    PhaseProperties spec_heat_capacities_;
};
}  // namespace MaterialPropertyLib