version 3.10.0
Loading...
Searching...
No Matches
diffusivityconstanttortuosity.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// SPDX-FileCopyrightText: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
5// SPDX-License-Identifier: GPL-3.0-or-later
6//
12
13#ifndef DUMUX_MATERIAL_DIFFUSIVITY_CONSTANT_TORTUOSITY_HH
14#define DUMUX_MATERIAL_DIFFUSIVITY_CONSTANT_TORTUOSITY_HH
15
17
18namespace Dumux {
19
24
38template<class Scalar>
40{
41public:
54 template<class VolumeVariables>
55 static Scalar effectiveDiffusionCoefficient(const VolumeVariables& volVars,
56 const int phaseIdx,
57 const int compIdxI,
58 const int compIdxJ)
59 {
60 static const Scalar tau = getParam<Scalar>("SpatialParams.Tortuosity", 0.5);
61 const Scalar diffCoeff = volVars.diffusionCoefficient(phaseIdx, compIdxI, compIdxJ);
62 return volVars.porosity() * volVars.saturation(phaseIdx) * tau * diffCoeff;
63 }
64
65};
66} // end namespace Dumux
67
68#endif
Relation for the effective diffusion coefficient with constant tortuosity.
Definition diffusivityconstanttortuosity.hh:40
static Scalar effectiveDiffusionCoefficient(const VolumeVariables &volVars, const int phaseIdx, const int compIdxI, const int compIdxJ)
Returns the effective diffusion coefficient ( )
Definition diffusivityconstanttortuosity.hh:55
T getParam(Args &&... args)
A free function to get a parameter from the parameter tree singleton.
Definition parameters.hh:139
Definition adapt.hh:17
The infrastructure to retrieve run-time parameters from Dune::ParameterTrees.