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

adapted tutorial decoupled to unique neumann-BC definition

Browse Source
This commit is contained in:
Benjamin Faigle 2011-01-03 10:56:26 +00:00 committed by Andreas Lauser
parent f4208fd27b
commit 0d7f8adc41
2 changed files with 3 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
doc/handbook/tutorial-decoupled.tex
View File

@ -166,7 +166,7 @@ Hence, they all feature a common argument list:
\end{itemize} \end{itemize}
In the following, there are the methods for general parameters, source- or In the following, there are the methods for general parameters, source- or
sinkterms, boundary conditions (lines \ref{tutorial-decoupled:bctypePress} to sinkterms, boundary conditions (lines \ref{tutorial-decoupled:bctypePress} to
\ref{tutorial-decoupled:neumannSat}) and initial values for the transported \ref{tutorial-decoupled:neumann}) and initial values for the transported
quantity in line \label{tutorial-decoupled:initSat}. For more information quantity in line \label{tutorial-decoupled:initSat}. For more information
on the functions, it is referred to the documentation in the code. on the functions, it is referred to the documentation in the code.

12
examples/tutorialproblem_decoupled.hh
View File

@ -253,11 +253,11 @@ public:
{ {
return 1; return 1;
} }
//! Value for pressure neumann boundary condition \f$ [\frac{kg}{m^3 \cdot s}] \f$. //! Value for neumann boundary condition \f$ [\frac{kg}{m^3 \cdot s}] \f$.
/*! In case of a neumann boundary condition, the flux of matter /*! In case of a neumann boundary condition, the flux of matter
* is returned as a vector. * is returned as a vector.
*/ */
std::vector<Scalar> neumannPress(const GlobalPosition& globalPos, const Intersection& intersection) const /*@\label{tutorial-decoupled:neumannPress}@*/ std::vector<Scalar> neumann(const GlobalPosition& globalPos, const Intersection& intersection) const /*@\label{tutorial-decoupled:neumann}@*/
{ {
std::vector<Scalar> neumannFlux(2,0.0); std::vector<Scalar> neumannFlux(2,0.0);
if (globalPos[0] > this->bboxMax()[0] - eps_) if (globalPos[0] > this->bboxMax()[0] - eps_)
@ -266,14 +266,6 @@ public:
} }
return neumannFlux; return neumannFlux;
} }
//! Value for transport neumann boundary condition \f$ [\frac{kg}{m^3 \cdot s}] \f$.
/*! In case of a neumann boundary condition for the transport equation
* the flux of matter for the primary variable is returned as a scalar.
*/
Scalar neumannSat(const GlobalPosition& globalPos, const Intersection& intersection, Scalar factor) const /*@\label{tutorial-decoupled:neumannSat}@*/
{
return 0;
}
//! Saturation initial condition (dimensionless) //! Saturation initial condition (dimensionless)
/*! The problem is initialized with the following saturation. /*! The problem is initialized with the following saturation.
*/ */