24#ifndef DUMUX_MATERIAL_THERMALCONDUCTIVITY_JOHANSEN_HH
25#define DUMUX_MATERIAL_THERMALCONDUCTIVITY_JOHANSEN_HH
64template<
class Scalar,
class Indices = JohansenIndices>
71 template<
class VolumeVariables,
class SpatialParams,
class Element,
class FVGeometry>
72 [[deprecated(
"Signature deprecated. Use signature with volume variables only!")]]
74 const SpatialParams& spatialParams,
75 const Element& element,
76 const FVGeometry& fvGeometry,
77 const typename FVGeometry::SubControlVolume& scv)
97 template<
class VolumeVariables>
100 using FluidSystem =
typename VolumeVariables::FluidSystem;
101 static_assert(FluidSystem::numPhases == 2,
"ThermalConductivitySomerton only works for two-phase fluid systems!");
104 const Scalar sw = volVars.saturation(Indices::wPhaseIdx);
105 const Scalar lambdaW = volVars.fluidThermalConductivity(Indices::wPhaseIdx);
106 const Scalar lambdaN = volVars.fluidThermalConductivity(Indices::nPhaseIdx);
107 const Scalar lambdaSolid = volVars.solidThermalConductivity();
108 const Scalar
porosity = volVars.porosity();
109 const Scalar rhoSolid = volVars.solidDensity();
127 const Scalar lambdaW,
128 const Scalar lambdaN,
129 const Scalar lambdaSolid,
131 const Scalar rhoSolid)
134 const Scalar satW = max<Scalar>(0.0, Sw);
136 const Scalar kappa = 15.6;
137 const Scalar rhoBulk = rhoSolid*
porosity;
141 const Scalar lSat = lambdaSolid * pow(lambdaW / lambdaSolid,
porosity);
142 const Scalar lDry = (0.135*rhoBulk + 64.7)/(rhoSolid - 0.947*rhoBulk);
143 const Scalar Ke = (kappa*satW)/(1+(kappa-1)*satW);
145 return lDry + Ke * (lSat - lDry);
make the local view function available whenever we use the grid geometry
Definition: adapt.hh:29
std::string porosity() noexcept
I/O name of porosity.
Definition: name.hh:139
Definition: thermalconductivityjohansen.hh:33
static const int wPhaseIdx
Definition: thermalconductivityjohansen.hh:34
static const int nPhaseIdx
Definition: thermalconductivityjohansen.hh:35
Relation for the saturation-dependent effective thermal conductivity.
Definition: thermalconductivityjohansen.hh:66
static Scalar effectiveThermalConductivity(const VolumeVariables &volVars)
Returns the effective thermal conductivity after Johansen (1975) .
Definition: thermalconductivityjohansen.hh:98
static Scalar effectiveThermalConductivity(const Scalar Sw, const Scalar lambdaW, const Scalar lambdaN, const Scalar lambdaSolid, const Scalar porosity, const Scalar rhoSolid)
Returns the effective thermal conductivity after Johansen (1975) .
Definition: thermalconductivityjohansen.hh:126
static Scalar effectiveThermalConductivity(const VolumeVariables &volVars, const SpatialParams &spatialParams, const Element &element, const FVGeometry &fvGeometry, const typename FVGeometry::SubControlVolume &scv)
effective thermal conductivity after Johansen (1975)
Definition: thermalconductivityjohansen.hh:73