OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

fixed some minor bugs in tutorial

Browse Source
This commit is contained in:
Benjamin Faigle 2010-11-19 14:43:58 +00:00 committed by Andreas Lauser
parent 5d3a12e7d9
commit be1a31aa25
3 changed files with 8 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
examples/tutorial_decoupled.cc
View File

@ -94,7 +94,7 @@ int main(int argc, char** argv)
Problem problem(gridPtr->leafView()); /*@\label{tutorial-decoupled:instantiate-problem}@*/
// load restart file if necessary
if (restart) /*@\label{tutorial-decoupled:restart}@*/
if (restart) /*@\label{tutorial-decoupled:mainRestart}@*/
problem.deserialize(restartTime);
// define simulation parameters

18
examples/tutorialproblem_decoupled.hh
View File

@ -97,11 +97,6 @@ public:
// Set the spatial parameters
SET_PROP(TutorialProblemDecoupled, SpatialParameters) /*@\label{tutorial-decoupled:set-spatialparameters}@*/
{
private:
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Grid)) Grid;
typedef typename GET_PROP_TYPE(TypeTag, PTAG(Scalar)) Scalar;
public:
typedef Dumux::TutorialSpatialParametersDecoupled<TypeTag> type;
};
@ -203,7 +198,8 @@ public:
//! Returns a constant pressure to enter material laws
/* For incrompressible simulations, a constant pressure is necessary
* to enter the material laws to gain a constant density etc.
* to enter the material laws to gain a constant density etc. In the compressible
* case, the pressure is used for the initialization of material laws.
*/
Scalar referencePressure(const GlobalPosition& globalPos, const Element& element) const /*@\label{tutorial-decoupled:refPressure}@*/
{
@ -225,7 +221,7 @@ public:
*/
typename BoundaryConditions::Flags bctypePress(const GlobalPosition& globalPos, const Intersection& intersection) const /*@\label{tutorial-decoupled:bctypePress}@*/
{
if ((globalPos[0] < this->bboxMin()[0] + eps_))
if (globalPos[0] < this->bboxMin()[0] + eps_)
return BoundaryConditions::dirichlet;
// all other boundaries
return BoundaryConditions::neumann;
@ -248,21 +244,15 @@ public:
*/
Scalar dirichletPress(const GlobalPosition& globalPos, const Intersection& intersection) const /*@\label{tutorial-decoupled:dirichletPress}@*/
{
if (globalPos[0] < this->bboxMin()[0] + eps_)
return 2e5;
// all other boundaries
return 0;
}
//! Value for transport dirichlet boundary condition (dimensionless).
/*! In case of a dirichlet BC for the transport equation, a saturation
* have to be defined on boundaries.
* has to be defined on boundaries.
*/
Scalar dirichletSat(const GlobalPosition& globalPos, const Intersection& intersection) const /*@\label{tutorial-decoupled:dirichletSat}@*/
{
if (globalPos[0] < this->bboxMin()[0] + eps_)
return 1;
// all other boundaries
return 0;
}
//! Value for pressure neumann boundary condition \f$ [\frac{kg}{m^3 \cdot s}] \f$.
/*! In case of a neumann boundary condition, the flux of matter

2
examples/tutorialspatialparameters_decoupled.hh
View File

@ -92,7 +92,7 @@ public:
// parameters for the Brooks-Corey Law
// entry pressures
materialLawParams_.setPe(10000);
materialLawParams_.setPe(1000);
// Brooks-Corey shape parameters
materialLawParams_.setAlpha(2);