Commit 6c33aa50 authored by Andreas Lauser's avatar Andreas Lauser
Browse files

fix some comments

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@4860 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 6b6813ab
......@@ -52,13 +52,13 @@ class TwoPProblem : public BoxProblem<TypeTag>
typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
public:
/*!
* \brief The constructor
*
* \param timeManager The time manager
* \param gridView The grid view
* \param verbose Turn verbosity on or off
*/
/*!
* \brief The constructor
*
* \param timeManager The time manager
* \param gridView The grid view
* \param verbose Turn verbosity on or off
*/
TwoPProblem(TimeManager &timeManager,
const GridView &gridView,
bool verbose = true)
......
......@@ -29,7 +29,12 @@ namespace Dumux
{
/*!
* \ingroup BoxModel
* \brief Base class for all problems which use the box scheme
* \brief Base class for all problems which use the box scheme.
*
* \note All quantities are specified assuming a threedimensional
* world. Problems discretized using 2D grids are assumed to be
* extruded by \f$1 m\f$ and 1D grids are assumed to have a
* cross section of \f$1m \times 1m\f$.
*/
template<class TypeTag>
class BoxProblem
......
......@@ -54,15 +54,15 @@ inline Scalar henryIAPWS(Scalar E,
Scalar Tr = temperature/H2O::criticalTemperature();
Scalar tau = 1 - Tr;
const Scalar c[6] = {
static const Scalar c[6] = {
1.99274064, 1.09965342, -0.510839303,
-1.75493479,-45.5170352, -6.7469445e5
};
const Scalar d[6] = {
static const Scalar d[6] = {
1/3.0, 2/3.0, 5/3.0,
16/3.0, 43/3.0, 110/3.0
};
const Scalar q = -0.023767;
static const Scalar q = -0.023767;
Scalar f = 0;
for (int i = 0; i < 6; ++i) {
......
......@@ -136,14 +136,17 @@ SET_BOOL_PROP(LensProblem, EnableGravity, true);
*
* The domain is sized 6m times 4m and features a rectangular lens
* with low permeablility which spans from (1 m , 2 m) to (4 m, 3 m)
* and is surrounded by a medium with higher permability.
* and is surrounded by a medium with higher permability. Note that
* this problem is discretized using only two dimensions, so from the
* point of view of the two-phase model, the depth of the domain
* implicitly is 1 m everywhere.
*
* On the top and the bottom of the domain neumann boundary conditions
* are used, while dirichlet conditions apply on the left and right
* boundaries.
*
* DNAPL is injected at the top boundary from 3m to 4m at a rate of
* 0.04 kg/(s m), the remaining neumann boundaries are no-flow
* 0.04 kg/(s m^2), the remaining neumann boundaries are no-flow
* boundaries.
*
* The dirichlet boundaries on the left boundary is the hydrostatic
......@@ -153,12 +156,12 @@ SET_BOOL_PROP(LensProblem, EnableGravity, true);
*
* This problem uses the \ref TwoPBoxModel.
*
* This problem should typically simulated until \f$t_{\text{end}} =
* 50\,000\;s\f$ is reached. A good choice for the initial time step size
* is \f$t_{\text{inital}} = 1\,000\;s\f$.
* This problem should typically be simulated until \f$t_{\text{end}}
* \approx 50\,000\;s\f$ is reached. A good choice for the initial time step
* size is \f$t_{\text{inital}} = 250\;s\f$.
*
* To run the simulation execute the following line in shell:
* <tt>./test_2p 50000 1000</tt>
* <tt>./test_2p 50000 250</tt>
*/
template <class TypeTag >
class LensProblem : public TwoPProblem<TypeTag>
......@@ -348,7 +351,7 @@ public:
* \brief Evaluate the boundary conditions for a neumann
* boundary segment.
*
* \param values The neumann values for the conservation equations
* \param values The neumann values for the conservation equations [kg / (m^2 *s )]
* \param element The finite element
* \param fvElemGeom The finite-volume geometry in the box scheme
* \param is The intersection between element and boundary
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment