3.6-git
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
freeflow/compositional/staggered/localresidual.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_FREEFLOW_NC_STAGGERED_LOCAL_RESIDUAL_HH
25#define DUMUX_FREEFLOW_NC_STAGGERED_LOCAL_RESIDUAL_HH
26
30
31namespace Dumux {
32
38// forward declaration
39template<class TypeTag, class DiscretizationMethod>
40class FreeflowNCResidualImpl;
41
42template<class TypeTag>
43class FreeflowNCResidualImpl<TypeTag, DiscretizationMethods::Staggered>
44: public NavierStokesResidual<TypeTag>
45{
47
51 using GridView = typename GridGeometry::GridView;
52 using Element = typename GridView::template Codim<0>::Entity;
53 using FVElementGeometry = typename GetPropType<TypeTag, Properties::GridGeometry>::LocalView;
54 using SubControlVolume = typename FVElementGeometry::SubControlVolume;
55 using CellCenterPrimaryVariables = GetPropType<TypeTag, Properties::CellCenterPrimaryVariables>;
57
58 static constexpr int numComponents = ModelTraits::numFluidComponents();
59 static constexpr bool useMoles = getPropValue<TypeTag, Properties::UseMoles>();
60
61 using EnergyLocalResidual = typename ParentType::EnergyLocalResidual;
62
63public:
64 using ParentType::ParentType;
65
67 template<class VolumeVariables>
68 CellCenterPrimaryVariables computeStorageForCellCenter(const Problem& problem,
69 const SubControlVolume& scv,
70 const VolumeVariables& volVars) const
71 {
72 CellCenterPrimaryVariables storage(0.0);
73
74 const Scalar density = useMoles ? volVars.molarDensity() : volVars.density();
75
76 // compute storage term of all components
77 for (int compIdx = 0; compIdx < numComponents; ++compIdx)
78 {
79 const int eqIdx = compIdx;
80
81 const Scalar massOrMoleFraction = useMoles ? volVars.moleFraction(compIdx) : volVars.massFraction(compIdx);
82 const Scalar s = density * massOrMoleFraction;
83
84 if (eqIdx != ModelTraits::replaceCompEqIdx())
85 storage[eqIdx] += s;
86 }
87
88 // in case one balance is substituted by the total mass balance
89 if(ModelTraits::replaceCompEqIdx() < numComponents)
90 storage[ModelTraits::replaceCompEqIdx()] = density;
91
92 EnergyLocalResidual::fluidPhaseStorage(storage, volVars);
93
94 return storage;
95 }
96};
97} // end namespace Dumux
98
99#endif
The available discretization methods in Dumux.
VolumeVariables::PrimaryVariables::value_type massOrMoleFraction(const VolumeVariables &volVars, ReferenceSystemFormulation referenceSys, const int phaseIdx, const int compIdx)
returns the mass or mole fraction to be used in Fick's law based on the reference system
Definition: referencesystemformulation.hh:66
Adaption of the non-isothermal two-phase two-component flow model to problems with CO2.
Definition: adapt.hh:29
EnergyLocalResidualImplementation< TypeTag, GetPropType< TypeTag, Properties::ModelTraits >::enableEnergyBalance()> EnergyLocalResidual
Definition: porousmediumflow/nonisothermal/localresidual.hh:39
typename GetProp< TypeTag, Property >::type GetPropType
get the type alias defined in the property
Definition: propertysystem.hh:180
std::string density(int phaseIdx) noexcept
I/O name of density for multiphase systems.
Definition: name.hh:65
Element-wise calculation of the multi-component free-flow residual for models using the staggered dis...
Definition: freeflow/compositional/localresidual.hh:36
CellCenterPrimaryVariables computeStorageForCellCenter(const Problem &problem, const SubControlVolume &scv, const VolumeVariables &volVars) const
Evaluate fluxes entering or leaving the cell center control volume.
Definition: freeflow/compositional/staggered/localresidual.hh:68
Element-wise calculation of the Navier-Stokes residual for models using the staggered discretization.
Definition: freeflow/navierstokes/localresidual.hh:35
Declares all properties used in Dumux.