From c8652107f1e60d82c12fbd89870609610c7fe8a1 Mon Sep 17 00:00:00 2001
From: Bernd Flemisch <bernd@iws.uni-stuttgart.de>
Date: Thu, 2 Dec 2010 11:14:57 +0000
Subject: [PATCH] adapted tutorial such that we have advection-dominated flow
again
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@4806 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
tutorial/tutorialproblem_coupled.hh | 10 +++++-----
tutorial/tutorialproblem_decoupled.hh | 4 ++--
tutorial/tutorialspatialparameters_coupled.hh | 2 +-
tutorial/tutorialspatialparameters_decoupled.hh | 2 +-
4 files changed, 9 insertions(+), 9 deletions(-)
diff --git a/tutorial/tutorialproblem_coupled.hh b/tutorial/tutorialproblem_coupled.hh
index bf94e53b54..c7651b62a3 100644
--- a/tutorial/tutorialproblem_coupled.hh
+++ b/tutorial/tutorialproblem_coupled.hh
@@ -57,8 +57,8 @@ SET_PROP(TutorialProblemCoupled, Grid) /*@\label{tutorial-coupled:set-grid}@*/
Dune::FieldVector<int, 2> cellRes;
Dune::FieldVector<ctype, 2> lowerLeft(0.0);
Dune::FieldVector<ctype, 2> upperRight;
- cellRes[0] = 30;
- cellRes[1] = 10;
+ cellRes[0] = 100;
+ cellRes[1] = 1;
upperRight[0] = 300;
upperRight[1] = 60;
return new Dune::SGrid<2,2>(cellRes,
@@ -192,10 +192,10 @@ public:
fvElemGeom.boundaryFace[boundaryFaceIdx].ipGlobal;
Scalar right = this->bboxMax()[0];
// extraction of oil on the right boundary for approx. 1.e6 seconds
- if (pos[0] > right - eps_ && this->timeManager().time() <= 1.e6) {
- // oil outflux of 0.3 g/(m * s) on the right boundary.
+ if (pos[0] > right - eps_) {
+ // oil outflux of 30 g/(m * s) on the right boundary.
values[Indices::contiWEqIdx] = 0;
- values[Indices::contiNEqIdx] = 3e-4;
+ values[Indices::contiNEqIdx] = 3e-2;
} else {
// no-flow on the remaining Neumann-boundaries.
values[Indices::contiWEqIdx] = 0;
diff --git a/tutorial/tutorialproblem_decoupled.hh b/tutorial/tutorialproblem_decoupled.hh
index a913708264..2398b274b1 100644
--- a/tutorial/tutorialproblem_decoupled.hh
+++ b/tutorial/tutorialproblem_decoupled.hh
@@ -184,7 +184,7 @@ public:
bool shouldWriteOutput() const /*@\label{tutorial-decoupled:output}@*/
{
return this->timeManager().timeStepIndex() > 0 &&
- (this->timeManager().timeStepIndex() % 1 == 0);
+ (this->timeManager().timeStepIndex() % 20 == 0);
}
//! Returns the temperature within the domain.
@@ -262,7 +262,7 @@ public:
std::vector<Scalar> neumannFlux(2,0.0);
if (globalPos[0] > this->bboxMax()[0] - eps_)
{
- neumannFlux[nPhaseIdx] = 3e-4;
+ neumannFlux[nPhaseIdx] = 3e-2;
}
return neumannFlux;
}
diff --git a/tutorial/tutorialspatialparameters_coupled.hh b/tutorial/tutorialspatialparameters_coupled.hh
index 88c6161670..b6f62f29aa 100644
--- a/tutorial/tutorialspatialparameters_coupled.hh
+++ b/tutorial/tutorialspatialparameters_coupled.hh
@@ -104,7 +104,7 @@ public:
materialParams_.setSnr(0.0);
//linear material law
- materialParams_.setPe(1000.0);
+ materialParams_.setPe(500.0);
materialParams_.setAlpha(2);
}
diff --git a/tutorial/tutorialspatialparameters_decoupled.hh b/tutorial/tutorialspatialparameters_decoupled.hh
index 997d7dc247..22479f852b 100644
--- a/tutorial/tutorialspatialparameters_decoupled.hh
+++ b/tutorial/tutorialspatialparameters_decoupled.hh
@@ -92,7 +92,7 @@ public:
// parameters for the Brooks-Corey Law
// entry pressures
- materialLawParams_.setPe(1000);
+ materialLawParams_.setPe(500);
// Brooks-Corey shape parameters
materialLawParams_.setAlpha(2);
--
GitLab