Skip to content
Snippets Groups Projects
Commit 736402b0 authored by Simon Emmert's avatar Simon Emmert
Browse files

[2p1c] cleanup fluidsystem and problem

parent b7120f26
No related branches found
No related tags found
1 merge request!1140Feature/2p1c wettability
......@@ -259,7 +259,7 @@ public:
if (Component::isTabulated)
{
Component::init(tempMin, tempMax, nTemp,
pressMin, pressMax, nPress);
pressMin, pressMax, nPress);
}
}
......@@ -287,8 +287,8 @@ public:
{
return Component::gasDensity(temperature, pressure);
}
else DUNE_THROW(Dune::NotImplemented,
"wrong index");
else
DUNE_THROW(Dune::NotImplemented, "wrong index");
}
using Base::molarDensity;
......@@ -310,7 +310,6 @@ public:
return Component::liquidMolarDensity(temperature, pressure);
return Component::gasMolarDensity(temperature, pressure);
}
/*!
......@@ -330,15 +329,16 @@ public:
Scalar pressure = fluidState.pressure(phaseIdx);
// liquid phase
if (phaseIdx == liquidPhaseIdx) {
if (phaseIdx == liquidPhaseIdx)
{
return Component::liquidViscosity(temperature, pressure);
}
else if (phaseIdx == gasPhaseIdx) // gas phase
{
return Component::gasViscosity(temperature, pressure) ;
return Component::gasViscosity(temperature, pressure);
}
else DUNE_THROW(Dune::NotImplemented,
"wrong index");
else
DUNE_THROW(Dune::NotImplemented, "wrong index");
}
/*!
......@@ -352,9 +352,9 @@ public:
const unsigned int phaseIdx)
{
assert(0 <= phaseIdx && phaseIdx < numPhases);
Scalar pressure = fluidState.pressure(gasPhaseIdx) ;
Scalar pressure = fluidState.pressure(gasPhaseIdx);
return Component::vaporTemperature(pressure) ;
return Component::vaporTemperature(pressure);
}
/*!
......@@ -444,9 +444,7 @@ public:
int compIIdx,
int compJIdx)
{
DUNE_THROW(Dune::NotImplemented, "Binary Diffusion coefficients");
}
{ DUNE_THROW(Dune::NotImplemented, "Binary Diffusion coefficients"); }
/*!
* \brief Calculate specific enthalpy [J/kg].
......@@ -462,7 +460,8 @@ public:
assert(0 <= phaseIdx && phaseIdx < numPhases);
// liquid phase
if (phaseIdx == liquidPhaseIdx) {
if (phaseIdx == liquidPhaseIdx)
{
return Component::liquidEnthalpy(fluidState.temperature(phaseIdx),
fluidState.pressure(phaseIdx));
}
......@@ -471,8 +470,8 @@ public:
return Component::gasEnthalpy(fluidState.temperature(phaseIdx),
fluidState.pressure(phaseIdx));
}
else DUNE_THROW(Dune::NotImplemented,
"wrong index");
else
DUNE_THROW(Dune::NotImplemented, "wrong index");
}
/*!
......@@ -499,8 +498,8 @@ public:
return Component::gasThermalConductivity(fluidState.temperature(phaseIdx),
fluidState.pressure(phaseIdx)); //0.0248;
}
else DUNE_THROW(Dune::NotImplemented,
"wrong index");
else
DUNE_THROW(Dune::NotImplemented, "wrong index");
}
/*!
......
......@@ -118,9 +118,7 @@ public:
*/
InjectionProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry)
: ParentType(fvGridGeometry)
{
FluidSystem::init();
}
{ FluidSystem::init(); }
/*!
* \name Problem parameters
......@@ -130,12 +128,10 @@ public:
//! \copydoc Dumux::FVProblem::source()
NumEqVector source(const Element &element,
const FVElementGeometry& fvGeometry,
const ElementVolumeVariables& elemVolVars,
const SubControlVolume &scv) const
{
return NumEqVector(0.0);
}
const FVElementGeometry& fvGeometry,
const ElementVolumeVariables& elemVolVars,
const SubControlVolume &scv) const
{ return NumEqVector(0.0); }
/*!
* \name Boundary conditions
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment