7b8951ed83
thermal laws are the heat conduction laws plus "solid energy laws" which can be used to specify the relations which govern the volumetric internal energy of the solid matrix of the porous medium.
65 lines
1.9 KiB
C++
65 lines
1.9 KiB
C++
// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
|
|
// vi: set et ts=4 sw=4 sts=4:
|
|
/*
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 2 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
Consult the COPYING file in the top-level source directory of this
|
|
module for the precise wording of the license and the list of
|
|
copyright holders.
|
|
*/
|
|
/*!
|
|
* \file
|
|
* \copydoc Opm::ConstantSolidHeatCapLaw
|
|
*/
|
|
#ifndef OPM_CONSTANT_SOLID_HEAT_CAP_LAW_HPP
|
|
#define OPM_CONSTANT_SOLID_HEAT_CAP_LAW_HPP
|
|
|
|
#include "ConstantSolidHeatCapLawParams.hpp"
|
|
|
|
#include <opm/material/densead/Math.hpp>
|
|
|
|
namespace Opm
|
|
{
|
|
/*!
|
|
* \ingroup material
|
|
*
|
|
* \brief Implements a solid energy storage law which assumes constant heat capacity.
|
|
*/
|
|
template <class ScalarT,
|
|
class ParamsT = ConstantSolidHeatCapLawParams<ScalarT> >
|
|
class ConstantSolidHeatCapLaw
|
|
{
|
|
public:
|
|
typedef ParamsT Params;
|
|
typedef typename Params::Scalar Scalar;
|
|
|
|
/*!
|
|
* \brief Given a fluid state, compute the volumetric internal energy of the solid matrix [W/m^3].
|
|
*/
|
|
template <class FluidState, class Evaluation = typename FluidState::Scalar>
|
|
static Evaluation solidInternalEnergy(const Params& params, const FluidState& fluidState)
|
|
{
|
|
const Evaluation& T = fluidState.temperature(/*phaseIdx=*/0);
|
|
|
|
Scalar cp = params.solidHeatCapacity();
|
|
|
|
return T*cp;
|
|
}
|
|
};
|
|
} // namespace Opm
|
|
|
|
#endif
|