Skip to content
Snippets Groups Projects
Commit 79b0318c authored by Natalie Schröder's avatar Natalie Schröder
Browse files

[darcyfluxvaribles] change flux and gradients from DimVector GlobalPosition...

[darcyfluxvaribles] change flux and gradients from DimVector GlobalPosition (from dim to dimWorld) (reviewed by gruenich)

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@12479 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 75a34a3f
No related branches found
No related tags found
No related merge requests found
...@@ -186,13 +186,13 @@ protected: ...@@ -186,13 +186,13 @@ protected:
idx++) // loop over adjacent vertices idx++) // loop over adjacent vertices
{ {
// FE gradient at vertex idx // FE gradient at vertex idx
const DimVector &feGrad = face().grad[idx]; const GlobalPosition &feGrad = face().grad[idx];
// index for the element volume variables // index for the element volume variables
int volVarsIdx = face().fapIndices[idx]; int volVarsIdx = face().fapIndices[idx];
// the pressure gradient // the pressure gradient
DimVector tmp(feGrad); GlobalPosition tmp(feGrad);
tmp *= elemVolVars[volVarsIdx].fluidState().pressure(phaseIdx); tmp *= elemVolVars[volVarsIdx].fluidState().pressure(phaseIdx);
potentialGrad_[phaseIdx] += tmp; potentialGrad_[phaseIdx] += tmp;
} }
...@@ -201,7 +201,7 @@ protected: ...@@ -201,7 +201,7 @@ protected:
if (GET_PARAM_FROM_GROUP(TypeTag, bool, Problem, EnableGravity)) if (GET_PARAM_FROM_GROUP(TypeTag, bool, Problem, EnableGravity))
{ {
// ask for the gravitational acceleration at the given SCV face // ask for the gravitational acceleration at the given SCV face
DimVector g(problem.gravityAtPos(face().ipGlobal)); GlobalPosition g(problem.gravityAtPos(face().ipGlobal));
// calculate the phase density at the integration point. we // calculate the phase density at the integration point. we
// only do this if the wetting phase is present in both cells // only do this if the wetting phase is present in both cells
...@@ -218,7 +218,7 @@ protected: ...@@ -218,7 +218,7 @@ protected:
const Scalar density = (fI*rhoI + fJ*rhoJ)/(fI + fJ); const Scalar density = (fI*rhoI + fJ*rhoJ)/(fI + fJ);
// make gravity acceleration a force // make gravity acceleration a force
DimVector f(g); GlobalPosition f(g);
f *= density; f *= density;
// calculate the final potential gradient // calculate the final potential gradient
......
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