3.5-git
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
simplefluidlumping.hh
Go to the documentation of this file.
1// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2// vi: set et ts=4 sw=4 sts=4:
3/*****************************************************************************
4 * See the file COPYING for full copying permissions. *
5 * *
6 * This program is free software: you can redistribute it and/or modify *
7 * it under the terms of the GNU General Public License as published by *
8 * the Free Software Foundation, either version 3 of the License, or *
9 * (at your option) any later version. *
10 * *
11 * This program is distributed in the hope that it will be useful, *
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
14 * GNU General Public License for more details. *
15 * *
16 * You should have received a copy of the GNU General Public License *
17 * along with this program. If not, see <http://www.gnu.org/licenses/>. *
18 *****************************************************************************/
24#ifndef DUMUX_MATERIAL_FLUIDMATRIX_THERMALCONDUCTIVITY_SIMPLE_FLUID_LUMPING_HH
25#define DUMUX_MATERIAL_FLUIDMATRIX_THERMALCONDUCTIVITY_SIMPLE_FLUID_LUMPING_HH
26
27#include <assert.h>
28#include <algorithm>
29
30namespace Dumux {
31
36template<class Scalar>
38{
39public:
46 template<class VolumeVariables>
47 static Scalar effectiveThermalConductivity(const VolumeVariables& volVars)
48 {
49 using FluidSystem = typename VolumeVariables::FluidSystem;
50 const Scalar sw = volVars.saturation(FluidSystem::phase0Idx);
51 const Scalar lambdaW = volVars.fluidThermalConductivity(FluidSystem::phase0Idx);
52 const Scalar lambdaN = volVars.fluidThermalConductivity(FluidSystem::phase1Idx);
53 const Scalar lambdaSolid = volVars.solidThermalConductivity();
54 const Scalar porosity = volVars.porosity();
55
56 return effectiveThermalConductivity_(sw, lambdaW, lambdaN, lambdaSolid, porosity);
57 }
58
59private:
71 static Scalar effectiveThermalConductivity_(const Scalar sw,
72 const Scalar lambdaW,
73 const Scalar lambdaN,
74 const Scalar lambdaSolid,
75 const Scalar porosity)
76 {
77 // Franz Lindner / Shi & Wang 2011
78 using std::max;
79 const Scalar satW = max<Scalar>(0.0, sw);
80 return porosity * ( (1. - satW) * lambdaN + satW * lambdaW ) + (1.0 - porosity) * lambdaSolid ; ; // arithmetic
81 }
82};
83} // end namespace Dumux
84#endif
Definition: adapt.hh:29
std::string porosity() noexcept
I/O name of porosity.
Definition: name.hh:139
Relation for the saturation-dependent effective thermal conductivity.
Definition: simplefluidlumping.hh:38
static Scalar effectiveThermalConductivity(const VolumeVariables &volVars)
Effective thermal conductivity .
Definition: simplefluidlumping.hh:47