Commit 512465bd authored by Andreas Lauser's avatar Andreas Lauser
Browse files

2pni box model: adapt to new structure

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@3959 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent b0b79fbe
// $Id: 2pnielementdata.hh 3357 2010-03-25 13:02:05Z lauser $
/*****************************************************************************
* Copyright (C) 2009 by Andreas Lauser *
* Institute of Hydraulic Engineering *
* University of Stuttgart, Germany *
* email: <givenname>.<name>@iws.uni-stuttgart.de *
* *
* This program 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, as long as this copyright notice *
* is included in its original form. *
* *
* This program is distributed WITHOUT ANY WARRANTY. *
*****************************************************************************/
/*!
* \file
*
* \brief This template class contains the quantities which are
* constant within a finite element in the non-isothermal
* two-phase model.
*/
#ifndef DUMUX_2PNI_ELEMENT_DATA_HH
#define DUMUX_2PNI_ELEMENT_DATA_HH
#include <dumux/boxmodels/2p/2pelementdata.hh>
namespace Dumux
{
/*!
* \ingroup TwoPNIBoxModel
* \brief This template class contains the quantities which are
* constant within a finite element in the non-isothermal two-phase model.
*
* For the plain non-isothermal two-phase model everything is given on the finite
* volumes, so this class is empty.
*/
template <class TypeTag>
class TwoPNIElementData : public TwoPElementData<TypeTag>
{
};
} // end namepace
#endif
// $Id: 2pnifluxdata.hh 3793 2010-06-25 15:38:42Z melanie $
// $Id: 2pnifluxvars.hh 3793 2010-06-25 15:38:42Z melanie $
/*****************************************************************************
* Copyright (C) 2008-2009 by Melanie Darcis *
* Copyright (C) 2008-2009 by Andreas Lauser *
......@@ -49,18 +49,18 @@ namespace Dumux
* the integration point, etc.
*/
template <class TypeTag>
class TwoPNIFluxData : public TwoPFluxData<TypeTag>
class TwoPNIFluxVars : public TwoPFluxVars<TypeTag>
{
typedef TwoPFluxData<TypeTag> ParentType;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef TwoPFluxVars<TypeTag> ParentType;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(GridView)) GridView;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Problem)) Problem;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(VertexData)) VertexData;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Problem)) Problem;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(SecondaryVars)) SecondaryVars;
typedef typename GridView::ctype CoordScalar;
typedef typename GridView::ctype CoordScalar;
typedef typename GridView::template Codim<0>::Entity Element;
typedef std::vector<VertexData> VertexDataArray;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(ElementSecondaryVars)) ElementSecondaryVars;
enum {
dim = GridView::dimension,
......@@ -70,17 +70,17 @@ class TwoPNIFluxData : public TwoPFluxData<TypeTag>
};
typedef typename GET_PROP_TYPE(TypeTag, PTAG(FVElementGeometry)) FVElementGeometry;
typedef typename FVElementGeometry::SubControlVolume SCV;
typedef typename FVElementGeometry::SubControlVolumeFace SCVFace;
typedef typename FVElementGeometry::SubControlVolume SCV;
typedef typename FVElementGeometry::SubControlVolumeFace SCVFace;
typedef Dune::FieldVector<CoordScalar, dimWorld> Vector;
typedef Dune::FieldVector<CoordScalar, dimWorld> Vector;
public:
TwoPNIFluxData(const Problem &problem,
TwoPNIFluxVars(const Problem &problem,
const Element &element,
const FVElementGeometry &elemGeom,
int scvfIdx,
const VertexDataArray &elemDat)
const ElementSecondaryVars &elemDat)
: ParentType(problem, element, elemGeom, scvfIdx, elemDat)
{
// calculate temperature gradient using finite element
......
// $Id: 2pniboxjacobian.hh 3840 2010-07-15 10:14:15Z bernd $
// $Id: 2pnilocalresidual.hh 3840 2010-07-15 10:14:15Z bernd $
/*****************************************************************************
* Copyright (C) 2008-2009 by Melanie Darcis *
* Copyright (C) 2008-2009 by Andreas Lauser *
......@@ -27,11 +27,11 @@
#include "2pniproperties.hh"
#include <dumux/boxmodels/2p/2pboxjacobian.hh>
#include <dumux/boxmodels/2p/2plocalresidual.hh>
#include <dumux/boxmodels/2pni/2pnielementdata.hh>
#include <dumux/boxmodels/2pni/2pnivertexdata.hh>
#include <dumux/boxmodels/2pni/2pnifluxdata.hh>
#include <dumux/boxmodels/2pni/2pnisecondaryvars.hh>
#include <dumux/boxmodels/2pni/2pnifluxvars.hh>
namespace Dumux
......@@ -43,17 +43,17 @@ namespace Dumux
*/
template<class TypeTag>
class TwoPNIBoxJacobian : public TwoPBoxJacobian<TypeTag>
class TwoPNILocalResidual : public TwoPLocalResidual<TypeTag>
{
typedef TwoPNIBoxJacobian<TypeTag> ThisType;
typedef TwoPBoxJacobian<TypeTag> ParentType;
typedef TwoPNILocalResidual<TypeTag> ThisType;
typedef TwoPLocalResidual<TypeTag> ParentType;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Problem)) Problem;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(GridView)) GridView;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Problem)) Problem;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(GridView)) GridView;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef typename GET_PROP(TypeTag, PTAG(SolutionTypes)) SolutionTypes;
typedef typename SolutionTypes::PrimaryVarVector PrimaryVarVector;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(PrimaryVarVector)) PrimaryVarVector;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(TwoPIndices)) Indices;
......@@ -70,21 +70,16 @@ class TwoPNIBoxJacobian : public TwoPBoxJacobian<TypeTag>
};
typedef typename GET_PROP_TYPE(TypeTag, PTAG(VertexData)) VertexData;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(FluxData)) FluxData;
typedef std::vector<VertexData> VertexDataArray;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(SecondaryVars)) SecondaryVars;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(FluxVars)) FluxVars;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(ElementSecondaryVars)) ElementSecondaryVars;
typedef Dune::FieldVector<Scalar, dim> LocalPosition;
typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
typedef Dune::FieldVector<Scalar, dim> LocalPosition;
typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
static const Scalar mobilityUpwindAlpha = GET_PROP_VALUE(TypeTag, PTAG(MobilityUpwindAlpha));
public:
TwoPNIBoxJacobian(Problem &problem)
: ParentType(problem)
{
};
/*!
* \brief Evaluate the amount all conservation quantites
* (e.g. phase mass) within a sub-control volume.
......@@ -102,8 +97,8 @@ public:
// used. The secondary variables are used accordingly. This
// is required to compute the derivative of the storage term
// using the implicit euler method.
const VertexDataArray &vertDatArray = usePrevSol ? this->prevElemDat_ : this->curElemDat_;
const VertexData &vertDat = vertDatArray[scvIdx];
const ElementSecondaryVars &vertDatArray = usePrevSol ? this->prevSecVars_() : this->curSecVars_();
const SecondaryVars &vertDat = vertDatArray[scvIdx];
// compute the energy storage
result[temperatureIdx] =
......@@ -125,7 +120,7 @@ public:
* This method is called by compute flux (base class)
*/
void computeAdvectiveFlux(PrimaryVarVector &flux,
const FluxData &fluxData) const
const FluxVars &fluxData) const
{
// advective mass flux
ParentType::computeAdvectiveFlux(flux, fluxData);
......@@ -134,8 +129,8 @@ public:
flux[energyEqIdx] = 0;
for (int phase = 0; phase < numPhases; ++phase) {
// vertex data of the upstream and the downstream vertices
const VertexData &up = this->curElemDat_[fluxData.upstreamIdx(phase)];
const VertexData &dn = this->curElemDat_[fluxData.downstreamIdx(phase)];
const SecondaryVars &up = this->curSecVars_(fluxData.upstreamIdx(phase));
const SecondaryVars &dn = this->curSecVars_(fluxData.downstreamIdx(phase));
flux[energyEqIdx] +=
fluxData.KmvpNormal(phase) * (
......@@ -156,7 +151,7 @@ public:
* the face of a sub-control volume.
*/
void computeDiffusiveFlux(PrimaryVarVector &flux,
const FluxData &fluxData) const
const FluxVars &fluxData) const
{
// diffusive mass flux
ParentType::computeDiffusiveFlux(flux, fluxData);
......
......@@ -24,12 +24,12 @@
* \brief Adaption of the BOX scheme to the non-isothermal twophase flow model.
*/
#ifndef DUMUX_NEW_2PNI_BOX_MODEL_HH
#define DUMUX_NEW_2PNI_BOX_MODEL_HH
#ifndef DUMUX_2PNI_MODEL_HH
#define DUMUX_2PNI_MODEL_HH
#include <dumux/boxmodels/2p/2pboxmodel.hh>
#include "2pniboxjacobian.hh"
#include "2pniboxproblem.hh"
#include <dumux/boxmodels/2p/2pmodel.hh>
#include "2pnilocalresidual.hh"
#include "2pniproblem.hh"
namespace Dumux {
......@@ -81,39 +81,8 @@ namespace Dumux {
*/
template<class TypeTag>
class TwoPNIBoxModel: public BoxScheme<TypeTag, TwoPNIBoxModel<TypeTag> > {
typedef TwoPNIBoxModel<TypeTag> ThisType;
typedef BoxScheme<TypeTag, ThisType> ParentType;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Problem)) Problem;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(LocalJacobian)) LocalJacobian;
public:
TwoPNIBoxModel(Problem &prob)
: ParentType(prob)
{
}
/*!
* \brief Append all quantities of interest which can be derived
* from the solution of the current time step to the VTK
* writer.
*/
template <class MultiWriter>
void addOutputVtkFields(MultiWriter &writer)
{
this->localJacobian().addOutputVtkFields(writer, this->curSol());
}
/*!
* \brief Calculate the phase masses in the system for
* the current timestep.
*/
void calculateMass(Dune::FieldVector<Scalar, 2> &mass)
{
this->localJacobian().calculateMass(this->curSol(), mass);
}
class TwoPNIModel: public TwoPModel<TypeTag>
{
};
}
......
......@@ -19,34 +19,34 @@
* \brief Base class for all problems which use the non-isothermal
* two-phase box model
*/
#ifndef DUMUX_2PNI_BOX_PROBLEM_HH
#define DUMUX_2PNI_BOX_PROBLEM_HH
#ifndef DUMUX_2PNI_PROBLEM_HH
#define DUMUX_2PNI_PROBLEM_HH
#include <dumux/boxmodels/2p/2pboxproblem.hh>
#include <dumux/boxmodels/2p/2pproblem.hh>
namespace Dumux
{
/*!
* \ingroup TwoPNIProblems
* \brief Base class for all problems which use the non-isothermal
* \brief Base class for all problems which use the non-isothermal
* two-phase box model.
*
* \todo Please doc me more!
*/
template<class TypeTag, class Implementation>
class TwoPNIBoxProblem : public TwoPBoxProblem<TypeTag, Implementation>
template<class TypeTag>
class TwoPNIProblem : public TwoPProblem<TypeTag>
{
typedef TwoPBoxProblem<TypeTag, Implementation> ParentType;
typedef TwoPProblem<TypeTag> ParentType;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(GridView)) GridView;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(TimeManager)) TimeManager;
public:
TwoPNIBoxProblem(const GridView &gridView)
: ParentType(gridView)
{
}
TwoPNIProblem(TimeManager &timeManager, const GridView &gridView)
: ParentType(timeManager, gridView)
{}
/*!
* \name Problem parameters
......
......@@ -25,9 +25,9 @@
#define DUMUX_2PNI_PROPERTIES_HH
#include <dumux/boxmodels/2p/2pproperties.hh>
#include "2pnivertexdata.hh"
#include "2pnielementdata.hh"
#include "2pnifluxdata.hh"
#include "2pnisecondaryvars.hh"
#include "2pnifluxvars.hh"
namespace Dumux
{
......@@ -40,19 +40,16 @@ namespace Dumux
// forward declarations
////////////////////////////////
template<class TypeTag>
class TwoPNIBoxModel;
class TwoPNIModel;
template<class TypeTag>
class TwoPNIBoxJacobian;
template <class TypeTag>
class TwoPNIVertexData;
class TwoPNILocalResidual;
template <class TypeTag>
class TwoPNIElementData;
class TwoPNISecondaryVars;
template <class TypeTag>
class TwoPNIFluxData;
class TwoPNIFluxVars;
/*!
* \brief Enumerations for the non-isothermal two-phase model
......@@ -92,20 +89,17 @@ SET_INT_PROP(BoxTwoPNI, NumEq, 3); //!< set the number of equations to 3
//! Use the 2pni local jacobian operator for the 2pni model
SET_TYPE_PROP(BoxTwoPNI,
LocalJacobian,
TwoPNIBoxJacobian<TypeTag>);
LocalResidual,
TwoPNILocalResidual<TypeTag>);
//! the Model property
SET_TYPE_PROP(BoxTwoPNI, Model, TwoPNIBoxModel<TypeTag>);
//! the VertexData property
SET_TYPE_PROP(BoxTwoPNI, VertexData, TwoPNIVertexData<TypeTag>);
SET_TYPE_PROP(BoxTwoPNI, Model, TwoPNIModel<TypeTag>);
//! the ElementData property
SET_TYPE_PROP(BoxTwoPNI, ElementData, TwoPNIElementData<TypeTag>);
//! the SecondaryVars property
SET_TYPE_PROP(BoxTwoPNI, SecondaryVars, TwoPNISecondaryVars<TypeTag>);
//! the FluxData property
SET_TYPE_PROP(BoxTwoPNI, FluxData, TwoPNIFluxData<TypeTag>);
//! the FluxVars property
SET_TYPE_PROP(BoxTwoPNI, FluxVars, TwoPNIFluxVars<TypeTag>);
//! The indices required by the non-isothermal two-phase model
SET_TYPE_PROP(BoxTwoPNI, TwoPIndices, TwoPNIIndices<0>);
......
// $Id: 2pnivertexdata.hh 3736 2010-06-15 09:52:10Z lauser $
// $Id: 2pnisecondaryvars.hh 3736 2010-06-15 09:52:10Z lauser $
/*****************************************************************************
* Copyright (C) 2008-2009 by Melanie Darcis *
* Copyright (C) 2009 by Andreas Lauser *
......@@ -19,10 +19,10 @@
* \brief Contains the quantities which are are constant within a
* finite volume in the non-isothermal two-phase model.
*/
#ifndef DUMUX_2PNI_VERTEX_DATA_HH
#define DUMUX_2PNI_VERTEX_DATA_HH
#ifndef DUMUX_2PNI_SECONDARY_VARS_HH
#define DUMUX_2PNI_SECONDARY_VARS_HH
#include <dumux/boxmodels/2p/2pvertexdata.hh>
#include <dumux/boxmodels/2p/2psecondaryvars.hh>
namespace Dumux
{
......@@ -33,14 +33,14 @@ namespace Dumux
* finite volume in the non-isothermal two-phase model.
*/
template <class TypeTag>
class TwoPNIVertexData : public TwoPVertexData<TypeTag>
class TwoPNISecondaryVars : public TwoPSecondaryVars<TypeTag>
{
typedef TwoPVertexData<TypeTag> ParentType;
typedef TwoPSecondaryVars<TypeTag> ParentType;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(GridView)) GridView;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(FVElementGeometry)) FVElementGeometry;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Problem)) Problem;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Problem)) Problem;
typedef typename GridView::template Codim<0>::Entity Element;
......@@ -54,33 +54,33 @@ class TwoPNIVertexData : public TwoPVertexData<TypeTag>
typedef typename GET_PROP_TYPE(TypeTag, PTAG(TwoPIndices)) Indices;
enum { temperatureIdx = Indices::temperatureIdx };
typedef typename GET_PROP_TYPE(TypeTag, PTAG(FluidSystem)) FluidSystem;
typedef typename GET_PROP(TypeTag, PTAG(SolutionTypes)) SolutionTypes;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(FluidSystem)) FluidSystem;
typedef typename GET_PROP(TypeTag, PTAG(ReferenceElements)) RefElemProp;
typedef typename RefElemProp::Container ReferenceElements;
typedef typename RefElemProp::Container ReferenceElements;
typedef typename SolutionTypes::PrimaryVarVector PrimaryVarVector;
typedef Dune::FieldVector<Scalar, numPhases> PhasesVector;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(PrimaryVarVector)) PrimaryVarVector;
typedef Dune::FieldVector<Scalar, numPhases> PhasesVector;
public:
/*!
* \brief Update all quantities for a given control volume.
*/
void update(const PrimaryVarVector &sol,
void update(const PrimaryVarVector &priVars,
const Problem &problem,
const Element &element,
const FVElementGeometry &elemGeom,
int vertIdx,
const Problem &problem,
int scvIdx,
bool isOldSol)
{
typedef Indices I;
// vertex update data for the mass balance
ParentType::update(sol,
ParentType::update(priVars,
problem,
element,
elemGeom,
vertIdx,
problem,
scvIdx,
isOldSol);
// the internal energies and the enthalpies
......@@ -101,14 +101,14 @@ public:
}
// this method gets called by the parent class
void updateTemperature_(const PrimaryVarVector &sol,
void updateTemperature_(const PrimaryVarVector &priVars,
const Element &element,
const FVElementGeometry &elemGeom,
int scvIdx,
int scvIdx,
const Problem &problem)
{
// retrieve temperature from solution vector
this->temperature_ = sol[temperatureIdx];
// retrieve temperature from primary variables
this->temperature_ = priVars[temperatureIdx];
heatCapacity_ =
problem.spatialParameters().heatCapacity(element, elemGeom, scvIdx);
......
......@@ -18,7 +18,7 @@
/**
* @file
* \ingroup TwoPTwoCBoxProblems
* @brief Definition of a problem, where air is injected under a low permeable layer
* @brief Definition of a problem, where air is injected under a low permeable layer
* @author Klaus Mosthaf, Andreas Lauser, Bernd Flemisch
*/
#ifndef DUMUX_INJECTIONPROBLEM_HH
......@@ -62,8 +62,8 @@ SET_PROP(InjectionProblem, LocalFEMSpace)
enum{dim = GridView::dimension};
public:
typedef Dune::PDELab::Q1LocalFiniteElementMap<Scalar,Scalar,dim> type; // for cubes
// typedef Dune::PDELab::P1LocalFiniteElementMap<Scalar,Scalar,dim> type; // for simplices
typedef Dune::PDELab::Q1LocalFiniteElementMap<Scalar,Scalar,dim> type; // for cubes
// typedef Dune::PDELab::P1LocalFiniteElementMap<Scalar,Scalar,dim> type; // for simplices
};
#endif
......@@ -108,13 +108,13 @@ SET_INT_PROP(InjectionProblem, NewtonLinearSolverVerbosity, 0);
* This problem uses the \ref TwoPTwoCBoxModel.
*/
template <class TypeTag = TTAG(InjectionProblem) >
class InjectionProblem : public TwoPTwoCBoxProblem<TypeTag, InjectionProblem<TypeTag> >
class InjectionProblem : public TwoPTwoCBoxProblem<TypeTag>
{
typedef InjectionProblem<TypeTag> ThisType;
typedef TwoPTwoCBoxProblem<TypeTag, ThisType> ParentType;
typedef InjectionProblem<TypeTag> ThisType;
typedef TwoPTwoCBoxProblem<TypeTag> ParentType;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(GridView)) GridView;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(GridView)) GridView;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
enum {
// Grid and world dimension
......@@ -135,18 +135,18 @@ class InjectionProblem : public TwoPTwoCBoxProblem<TypeTag, InjectionProblem<Typ
contiGEqIdx = Indices::contiGEqIdx,
};
typedef typename GET_PROP(TypeTag, PTAG(SolutionTypes)) SolutionTypes;
typedef typename SolutionTypes::PrimaryVarVector PrimaryVarVector;
typedef typename SolutionTypes::BoundaryTypeVector BoundaryTypeVector;
typedef typename GridView::template Codim<0>::Entity Element;
typedef typename GridView::template Codim<dim>::Entity Vertex;
typedef typename GridView::Intersection Intersection;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(PrimaryVarVector)) PrimaryVarVector;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(BoundaryTypes)) BoundaryTypes;
typedef typename GridView::template Codim<0>::Entity Element;
typedef typename GridView::template Codim<dim>::Entity Vertex;
typedef typename GridView::Intersection Intersection;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(FVElementGeometry)) FVElementGeometry;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(FluidSystem)) FluidSystem;
typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
public:
InjectionProblem(const GridView &gridView)
......@@ -177,7 +177,7 @@ public:
*/
Scalar temperature(const Element &element,
const FVElementGeometry &fvElemGeom,
int scvIdx) const
int scvIdx) const
{
return temperature_;
};
......@@ -193,12 +193,12 @@ public:
* \brief Specifies which kind of boundary condition should be
* used for which equation on a given boundary segment.
*/
void boundaryTypes(BoundaryTypeVector &values,
void boundaryTypes(BoundaryTypes &values,
const Element &element,
const FVElementGeometry &fvElemGeom,
const Intersection &is,
int scvIdx,
int boundaryFaceIdx) const
int scvIdx,
int boundaryFaceIdx) const
{
const GlobalPosition &globalPos = element.geometry().corner(scvIdx);
......@@ -218,8 +218,8 @@ public:
const Element &element,
const FVElementGeometry &fvElemGeom,
const Intersection &is,
int scvIdx,
int boundaryFaceIdx) const
int scvIdx,
int boundaryFaceIdx) const
{
const GlobalPosition &globalPos = element.geometry().corner(scvIdx);
......@@ -238,8 +238,8 @@ public:
const Element &element,
const FVElementGeometry &fvElemGeom,
const Intersection &is,
int scvIdx,
int boundaryFaceIdx) const
int scvIdx,
int boundaryFaceIdx) const
{
const GlobalPosition &globalPos = element.geometry().corner(scvIdx);
......@@ -268,7 +268,7 @@ public:
void source(PrimaryVarVector &values,
const Element &element,
const FVElementGeometry &fvElemGeom,
int scvIdx) const
int scvIdx) const
{
values = Scalar(0.0);
}
......@@ -282,7 +282,7 @@ public:
void initial(PrimaryVarVector &values,
const Element &element,