From 030587c3502033230f408d410c55c807a29115cc Mon Sep 17 00:00:00 2001
From: Klaus Mosthaf <klmos@env.dtu.dk>
Date: Wed, 11 Apr 2012 06:22:08 +0000
Subject: [PATCH] reimplemented commit #7961: removed calculation of unused
velocityDivergence (commented) and massFractionGradient (moleFractionGradient
is used for the diffusion); unfortunately, deprecation was not possible in
this case
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@8067 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
dumux/freeflow/stokes/stokesfluxvariables.hh | 20 +++++++++----------
.../stokes2c/stokes2cfluxvariables.hh | 13 +-----------
2 files changed, 11 insertions(+), 22 deletions(-)
diff --git a/dumux/freeflow/stokes/stokesfluxvariables.hh b/dumux/freeflow/stokes/stokesfluxvariables.hh
index e2520ca192..cd0937d634 100644
--- a/dumux/freeflow/stokes/stokesfluxvariables.hh
+++ b/dumux/freeflow/stokes/stokesfluxvariables.hh
@@ -98,7 +98,7 @@ protected:
velocityAtIP_ = Scalar(0);
pressureGradAtIP_ = Scalar(0);
velocityGradAtIP_ = Scalar(0);
- velocityDivAtIP_ = Scalar(0);
+// velocityDivAtIP_ = Scalar(0);
for (int idx = 0;
idx < fvGeom_.numVertices;
@@ -136,7 +136,7 @@ protected:
tmp *= elemVolVars[idx].velocity()[dimIdx];
velocityGradAtIP_[dimIdx] += tmp;
- velocityDivAtIP_ += face().grad[idx][dimIdx]*elemVolVars[idx].velocity()[dimIdx];
+// velocityDivAtIP_ += face().grad[idx][dimIdx]*elemVolVars[idx].velocity()[dimIdx];
}
}
@@ -148,7 +148,7 @@ protected:
Valgrind::CheckDefined(velocityAtIP_);
Valgrind::CheckDefined(pressureGradAtIP_);
Valgrind::CheckDefined(velocityGradAtIP_);
- Valgrind::CheckDefined(velocityDivAtIP_);
+// Valgrind::CheckDefined(velocityDivAtIP_);
};
void determineUpwindDirection_(const ElementVolumeVariables &elemVolVars)
@@ -238,12 +238,12 @@ public:
const VectorGradient &velocityGradAtIP() const
{ return velocityGradAtIP_; }
- /*!
- * \brief Return the divergence of the normal velocity at the
- * integration point.
- */
- Scalar velocityDivAtIP() const
- { return velocityDivAtIP_; }
+// /*!
+// * \brief Return the divergence of the normal velocity at the
+// * integration point.
+// */
+// Scalar velocityDivAtIP() const
+// { return velocityDivAtIP_; }
/*!
* \brief Return the local index of the upstream sub-control volume.
@@ -274,7 +274,7 @@ protected:
Scalar viscosityAtIP_;
Scalar pressureAtIP_;
Scalar normalVelocityAtIP_;
- Scalar velocityDivAtIP_;
+// Scalar velocityDivAtIP_;
VelocityVector velocityAtIP_;
// gradients at the IPs
diff --git a/dumux/freeflow/stokes2c/stokes2cfluxvariables.hh b/dumux/freeflow/stokes2c/stokes2cfluxvariables.hh
index 631b6b2e5d..a457d1a351 100644
--- a/dumux/freeflow/stokes2c/stokes2cfluxvariables.hh
+++ b/dumux/freeflow/stokes2c/stokes2cfluxvariables.hh
@@ -97,12 +97,6 @@ public:
Scalar diffusionCoeffAtIP() const
{ return diffusionCoeffAtIP_; }
- /*!
- * \brief Return the gradient of the mass fraction at the integration point.
- */
- const ScalarGradient &massFractionGradAtIP() const
- { return massFractionGradAtIP_; }
-
/*!
* \brief Return the gradient of the mole fraction at the integration point.
*/
@@ -117,7 +111,6 @@ protected:
{
massFractionAtIP_ = Scalar(0);
diffusionCoeffAtIP_ = Scalar(0);
- massFractionGradAtIP_ = Scalar(0);
moleFractionGradAtIP_ = Scalar(0);
// calculate gradients and secondary variables at IPs
@@ -133,9 +126,6 @@ protected:
// the gradient of the mass fraction at the IP
for (int dimIdx=0; dimIdx<dim; ++dimIdx)
{
- massFractionGradAtIP_ +=
- this->face().grad[idx][dimIdx]*
- elemVolVars[idx].fluidState().massFraction(phaseIdx, lCompIdx);
moleFractionGradAtIP_ +=
this->face().grad[idx][dimIdx] *
elemVolVars[idx].fluidState().moleFraction(phaseIdx, lCompIdx);
@@ -144,12 +134,11 @@ protected:
Valgrind::CheckDefined(massFractionAtIP_);
Valgrind::CheckDefined(diffusionCoeffAtIP_);
- Valgrind::CheckDefined(massFractionGradAtIP_);
+ Valgrind::CheckDefined(moleFractionGradAtIP_);
}
Scalar massFractionAtIP_;
Scalar diffusionCoeffAtIP_;
- ScalarGradient massFractionGradAtIP_;
ScalarGradient moleFractionGradAtIP_;
};
--
GitLab