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

corrections in the coupled tutorial, spell check & beautification

Browse Source
This commit is contained in:
Klaus Mosthaf 2012-02-23 14:48:40 +00:00 committed by Andreas Lauser
parent ef35f75d0e
commit e8323070d0
3 changed files with 75 additions and 67 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

128
doc/handbook/tutorial-coupled.tex
View File

@ -9,7 +9,7 @@ The process of setting up a problem using \Dumux can be roughly divided into fou
\end{enumerate} \end{enumerate}
The problem being solved in this tutorial is illustrated in Figure \ref{tutorial-coupled:problemfigure}. The problem being solved in this tutorial is illustrated in Figure \ref{tutorial-coupled:problemfigure}.
A rectangular domain with no flow boundaries on the top and on the bottom, which is initially saturated with oil, is considered. A rectangular domain with no-flow boundaries on the top and on the bottom, which is initially saturated with oil, is considered.
Water infiltrates from the left side into the domain and replaces the oil. Gravity effects are neglected here. Water infiltrates from the left side into the domain and replaces the oil. Gravity effects are neglected here.
\begin{figure}[ht] \begin{figure}[ht]
@ -62,23 +62,26 @@ above.
\end{lst} \end{lst}
From line \ref{tutorial-coupled:include-begin} to line From line \ref{tutorial-coupled:include-begin} to line
\ref{tutorial-coupled:include-end} the headers required are included. \ref{tutorial-coupled:include-end} the required headers are included.
At line \ref{tutorial-coupled:set-type-tag} the type tag of the At line \ref{tutorial-coupled:set-type-tag} the type tag of the
problem, which is going to be simulated, is specified. All other data problem, which is going to be simulated, is specified. All other data
types can be retrieved via the \Dumux property system and only depend types can be retrieved via the \Dumux property system and only depend
on this single type tag. For a more thourough introduction to the on this single type tag. For a more thorough introduction to the
\Dumux property system, see chapter~\ref{sec:propertysystem}. \Dumux property system, see chapter~\ref{sec:propertysystem}.
After this \Dumux' default startup routine \texttt{Dumux::start()} is After this, the default startup routine \texttt{Dumux::start()} is
called on line \ref{tutorial-coupled:call-start}. This function deals called on line \ref{tutorial-coupled:call-start}. This function deals
with parsing the command line arguments, reading the parameter file, with parsing the command line arguments, reading the parameter file,
setting up the infrastructure necessary for \Dune, loads the grid, and setting up the infrastructure necessary for \Dune, loading the grid, and
starts the simulation. When it comes to parameters, all parameters can starting the simulation.
be either specified by command line arguments of the form Required parameters for the start of the simulation,
such as the initial time-step size, the simulation time or details of the grid,
can be either specified by command line arguments of the form
(\texttt{-ParameterName ParameterValue}), in the file specified by the (\texttt{-ParameterName ParameterValue}), in the file specified by the
\texttt{-parameterFile} argument, or if the latter is not specified, \texttt{-parameterFile} argument, or if the latter is not specified,
in the file \texttt{tutorial\_coupled.input}. If a parameter gets in the file \mbox{\texttt{tutorial\_coupled.input}}.
If a parameter is
specified on the command line as well as in the parameter file, the specified on the command line as well as in the parameter file, the
values provided in the command line have values provided in the command line have
precedence. Listing~\ref{tutorial-coupled:parameter-file} shows the precedence. Listing~\ref{tutorial-coupled:parameter-file} shows the
@ -91,10 +94,10 @@ default parameter file for the tutorial problem.
To provide an error message, the usage message which is displayed to To provide an error message, the usage message which is displayed to
the user if the simulation is called incorrectly, is printed via the the user if the simulation is called incorrectly, is printed via the
custom function which is defined on custom function which is defined on
line~\ref{tutorial-coupled:usage-function}. In this function the usage line~\ref{tutorial-coupled:usage-function}.
message is customized to the problem at hand. This means that at least In this function the usage message is customized to the problem at hand.
the necessary parameters are listed here. For more information about This means that at least the necessary parameters are listed here.
the input file please refer to section \ref{sec:inputFiles}. For more information about the input file please refer to section \ref{sec:inputFiles}.
\subsection{The Problem Class} \subsection{The Problem Class}
@ -127,13 +130,12 @@ this grid creator the physical domain of the grid is specified via the
run-time parameters \texttt{Grid.upperRightX}, run-time parameters \texttt{Grid.upperRightX},
\texttt{Grid.upperRightY}, \texttt{Grid.numberOfCellsX} and \texttt{Grid.upperRightY}, \texttt{Grid.numberOfCellsX} and
\texttt{Grid.numberOfCellsY}. These parameters can be specified via \texttt{Grid.numberOfCellsY}. These parameters can be specified via
the command-line or in a parameter file, the command-line or in a parameter file.
see~\ref{tutorial-coupled:runtime-parameters}.
Next, the appropriate fluid system, which specifies the thermodynamic Next, the appropriate fluid system, which specifies the thermodynamic
relations of the fluid phases, has to be chosen. By default, the relations of the fluid phases, has to be chosen. By default, the
two-phase model uses the \texttt{TwoPImmiscibleFluidSystem}, which two-phase model uses the \texttt{TwoPImmiscibleFluidSystem}, which
assumes immiscibility of the phases, but requires that the components assumes immiscibility of the phases, but requires the components
used for the wetting and non-wetting phases to be explicitly set. In used for the wetting and non-wetting phases to be explicitly set. In
this case, liquid water which uses the relations from this case, liquid water which uses the relations from
IAPWS'97~\cite{IAPWS1997} is chosen as the wetting phase on line IAPWS'97~\cite{IAPWS1997} is chosen as the wetting phase on line
@ -181,18 +183,18 @@ available:
specified for a vertex, which represents a control volume in the box specified for a vertex, which represents a control volume in the box
discretization. This avoids the specification of two different discretization. This avoids the specification of two different
boundary condition types for one equation at different control boundary condition types for one equation at different control
volume. Be aware that the second parameter is a Dune grid entity volumes. Be aware that the second parameter is a Dune grid entity
with the codimension \texttt{dim}. with the codimension \texttt{dim}.
\end{description} \end{description}
To ensure that no boundaries are undefined, a small safeguard value To ensure that no boundaries are undefined, a small safeguard value
\texttt{eps\_} is usually added when comparing spatial \texttt{eps\_} is usually added when comparing spatial
coordinates. The left boundary is hence detected by comparing the coordinates. The left boundary is hence not detected by checking, if the
first coordinate is equal with zero, but by testing whether it is first coordinate of the global position is equal to zero, but by testing whether it is
smaller than a very small value \texttt{eps\_}. smaller than a very small value \texttt{eps\_}.
Methods which make statements about boundary segments of the grid Methods which make statements about boundary segments of the grid
(i.e. \texttt{neumann()}) get called with six arguments: (such as \texttt{neumann()}) are called with six arguments:
\begin{description} \begin{description}
\item[values:] A vector \texttt{neumann()}, in which the mass fluxes per area unit \item[values:] A vector \texttt{neumann()}, in which the mass fluxes per area unit
over the boundary segment are specified. over the boundary segment are specified.
@ -202,7 +204,7 @@ Methods which make statements about boundary segments of the grid
finite element by the box scheme. finite element by the box scheme.
\item[isIt:] The \texttt{Intersection} of the boundary segment as given by the grid. \item[isIt:] The \texttt{Intersection} of the boundary segment as given by the grid.
\item[scvIdx:] The index of the sub-control volume in \item[scvIdx:] The index of the sub-control volume in
\texttt{fvElementGeometry} adjacent to the boundary segment. \texttt{fvElementGeometry} which is assigned to the boundary segment.
\item[boundaryFaceIdx:] The index of the boundary face in \item[boundaryFaceIdx:] The index of the boundary face in
\texttt{fvElementGeometry} which represents the boundary segment. \texttt{fvElementGeometry} which represents the boundary segment.
\end{description} \end{description}
@ -219,9 +221,9 @@ in \textsc{Kelvin} of the fluids and the rock matrix in the
domain. This temperature is then used by the model to calculate fluid domain. This temperature is then used by the model to calculate fluid
properties which possibly depend on it, e.g. density. The properties which possibly depend on it, e.g. density. The
\texttt{bboxMax()} (``\textbf{max}imum coordinated of the grid's \texttt{bboxMax()} (``\textbf{max}imum coordinated of the grid's
\textbf{b}ounding \textbf{b}ox'') method that is used here to \textbf{b}ounding \textbf{b}ox'') method is used here to
determine the extend of the phsical domain, returns a vector with the determine the extend of the physical domain. It returns a vector with the
maximum values of each coordinate of the grid's physical. This method maximum values of each global coordinate of the grid. This method
and the analogous \texttt{bboxMin()} method are provided by the base and the analogous \texttt{bboxMin()} method are provided by the base
class \texttt{Dumux::BoxProblem<TypeTag>}. class \texttt{Dumux::BoxProblem<TypeTag>}.
@ -243,20 +245,20 @@ interactions are defined by {\em fluid systems}, which are located in
% In this example, a class for the definition of a two-phase system is used. This allows for the choice % In this example, a class for the definition of a two-phase system is used. This allows for the choice
% of the two components oil and water and for access of the parameters that are relevant for the two-phase model. % of the two components oil and water and for access of the parameters that are relevant for the two-phase model.
\subsection{Definiting Spatially Dependent Parameters}\label{tutorial-coupled:description-spatialParameters} \subsection{Defining Spatially Dependent Parameters}\label{tutorial-coupled:description-spatialParameters}
In \Dumux, many properties of the porous medium can depend on the In \Dumux, many properties of the porous medium can depend on the
spatial location. Such properties are the \textit{intrinsic spatial location. Such properties are the \textit{intrinsic
permeability}, the parameters of the \textit{capillary pressure} and permeability}, the parameters of the \textit{capillary pressure} and
the \textit{relative permeability}, the \textit{porosity}, the the \textit{relative permeability}, the \textit{porosity}, the
\textit{heat capacity} as well as the \textit{heat conductivity}. Such \textit{heat capacity} as well as the \textit{heat conductivity}. Such
parameters are define using a so-called \textit{spatial parameters} parameters are defined using a so-called \textit{spatial parameters}
class. class.
The the box discretization is to be used, the spatial paramters class If the box discretization is to used, the spatial parameters class
should be derived from the base class should be derived from the base class
\texttt{Dumux::BoxSpatialParameters<TypeTag>}. Listing \texttt{Dumux::BoxSpatialParameters<TypeTag>}. Listing
\ref{tutorial-coupled:spatialparametersfile} shows the file \ref{tutorial-coupled:spatialparametersfile} shows the file \\
\verb+tutorialspatialparameters_coupled.hh+: \verb+tutorialspatialparameters_coupled.hh+:
\begin{lst}[File tutorial/tutorialspatialparameters\_coupled.hh]\label{tutorial-coupled:spatialparametersfile} \mbox{} \begin{lst}[File tutorial/tutorialspatialparameters\_coupled.hh]\label{tutorial-coupled:spatialparametersfile} \mbox{}
@ -266,11 +268,11 @@ numberstyle=\tiny, numbersep=5pt, firstline=28]{../../tutorial/tutorialspatialpa
First, the spatial parameters type tag is created on line First, the spatial parameters type tag is created on line
\ref{tutorial-coupled:define-spatialparameters-typetag}. The type tag \ref{tutorial-coupled:define-spatialparameters-typetag}. The type tag
for the problem then derives from it. The \Dumux properties defined on for the problem is then derived from it. The \Dumux properties defined on
the type tag for the spatial parameters are for example, the spatial the type tag for the spatial parameters are, for example, the spatial
parameters class itself (line parameters class itself (line
\ref{tutorial-coupled:set-spatialparameters}) or the capillary \ref{tutorial-coupled:set-spatialparameters}) or the capillary
pressure/relative permability relations\footnote{Taken together, the pressure/relative permeability relations\footnote{Taken together, the
capillary pressure and the relative permeability relations are capillary pressure and the relative permeability relations are
called \textit{material law}.} which ought to be used by the called \textit{material law}.} which ought to be used by the
simulation (line simulation (line
@ -279,17 +281,17 @@ simulation (line
\verb+dumux/material/fluidmatrixinteractions+. The selected one -- \verb+dumux/material/fluidmatrixinteractions+. The selected one --
here it is a relation according to a regularized version of here it is a relation according to a regularized version of
\textsc{Brooks} \& \textsc{Corey} -- is included in line \textsc{Brooks} \& \textsc{Corey} -- is included in line
\ref{tutorial-coupled:rawLawInclude}. After the selection, an adapter \ref{tutorial-coupled:rawLawInclude}.
class is specified ob line \ref{tutorial-coupled:eff2abs} to After the selection, an adapter class is specified in line \ref{tutorial-coupled:eff2abs} to
translates between effective and absolute saturations. This way, translate between effective and absolute saturations. Like this,
residual saturations can be specified in a generic way. As only used residual saturations can be specified in a generic way. As only the employed
material law knows which the names of the parameters which it material law knows the names of the parameters which it
requires, it provides a parameter class requires, it provides a parameter class
\texttt{RegularizedBrooksCoreyParams} which is exported by the type \texttt{RegularizedBrooksCoreyParams} which has the type
\texttt{Params} and defined in line \texttt{Params} and which is defined in line
\ref{tutorial-coupled:matLawObjectType}. In this case, the spatial \ref{tutorial-coupled:matLawObjectType}. In this case, the spatial
parameters only require a single set of parameters which means that it parameters only require a single set of parameters which means that it
only requires a single material parameter object as can be seen on only requires a single material parameter object as can be seen in
line~\ref{tutorial-coupled:matParamsObject}. line~\ref{tutorial-coupled:matParamsObject}.
In line \ref{tutorial-coupled:permeability}, a method returning the In line \ref{tutorial-coupled:permeability}, a method returning the
@ -298,18 +300,18 @@ to be called with three arguments:
\begin{description} \begin{description}
\item[\texttt{element}:] Just like for the problem itself, this \item[\texttt{element}:] Just like for the problem itself, this
parameter describes the considered element by means of a \Dune parameter describes the considered element by means of a \Dune
entity. Elements provide information about its geometry and entity. Elements provide information about their geometry and
position and can be mapped to a global index. position and can be mapped to a global index.
\item[\texttt{fvElemGeom}:] Holds information about the finite-volume \item[\texttt{fvElemGeom}:] It holds information about the finite-volume
geometry induced by the box-method on the element. geometry of the element induced by the box method.
\item[\texttt{scvIdx}:] Is the index of the sub-control volume of the \item[\texttt{scvIdx}:] This is the index of the sub-control volume of the
element which is considered. This is equivalent to the local index element which is considered. It is equivalent to the local index
of the vertex which corrosponts to the considered control volume in of the vertex which corresponds to the considered control volume in
the element. the element.
\end{description} \end{description}
The intrinsic permeability is a tensor and is thus the method returnes The intrinsic permeability is usually a tensor. Thus the method returns
a $\texttt{dim} \times \texttt{dim}$-matrix where \texttt{dim} is the a $\texttt{dim} \times \texttt{dim}$-matrix, where \texttt{dim} is the
dimension of the grid. dimension of the grid.
The method \texttt{porosity()} defined in line The method \texttt{porosity()} defined in line
@ -317,9 +319,9 @@ The method \texttt{porosity()} defined in line
\texttt{intrinsicPermeability()} and returns a scalar value for \texttt{intrinsicPermeability()} and returns a scalar value for
porosity dependent on the position in the domain. porosity dependent on the position in the domain.
Next, the method \texttt{materialLawParams()}, defined on line Next, the method \texttt{materialLawParams()}, defined in line
\ref{tutorial-coupled:matLawParams}, specifies \ref{tutorial-coupled:matLawParams}, returns the
\verb+materialLawParams+ object applies at the specified \verb+materialLawParams+ object that is applied at the specified
position. Although in this case only one object is returned, in position. Although in this case only one object is returned, in
general, the problem may be heterogeneous, which necessitates general, the problem may be heterogeneous, which necessitates
returning different objects at different positions in space. While returning different objects at different positions in space. While
@ -363,13 +365,13 @@ to make only some small changes in the tutorial files.
Compile the main file by typing \texttt{make tutorial\_coupled} and Compile the main file by typing \texttt{make tutorial\_coupled} and
run the model as explained above. run the model as explained above.
\item \textbf{Changing the shape of the discrete elements} \\ \item \textbf{Changing the Shape of the Discrete Elements} \\
Change the types of elements used for discretizing the domain. In line \ref{tutorial-coupled:set-gridcreator} of the problem file the type of gridcreator is chosen. By choosing a different grid creator you can discretize the domain with different elements. Hint: You can find gridcreators in \texttt{dumux/common/}. The shape of the employed elements can be visualized in paraview by choosing \texttt{Surface with Edges}. Change the types of elements used for discretizing the domain. In line \ref{tutorial-coupled:set-gridcreator} of the problem file the type of gridcreator is chosen. By choosing a different grid creator you can discretize the domain with different elements. Hint: You can find gridcreators in \texttt{dumux/common/}. The shape of the employed elements can be visualized in paraview by choosing \texttt{Surface with Edges}.
\item \textbf{Changing Fluids} \\ \item \textbf{Changing Fluids} \\
Now you can change the fluids. Use DNAPL instead of Oil and Brine instead of Water. To do that, you have to select different components via the property system in the problem file: Now you can change the fluids. Use DNAPL instead of Oil and Brine instead of Water. To do that, you have to select different components via the property system in the problem file:
\begin{enumerate} \begin{enumerate}
\item Brine: Brine is thermodynamically very similar to pure water but also considers a fixed amount of salt in the liquid phase. Hence, the class \texttt{Dumux::Brine} uses a pure water class, such as \texttt{Dumux::H2O}, as a second template argument after the data type \texttt{<Scalar>} as a template argument. \item Brine: Brine is thermodynamically very similar to pure water but also considers a fixed amount of salt in the liquid phase. Hence, the class \texttt{Dumux::Brine} uses a pure water class, such as \texttt{Dumux::H2O}, as a second template argument after the data type \texttt{<Scalar>}.
\item DNAPL: A standard set of chemical substances, such as Oil and Brine, is already included in the problem (via a list of \texttt{\#include ..} statements) and hence easily accessible by default. However, this is not the case for the class \texttt{Dumux::SimpleDNAPL}, which describes a simple \textbf{d}ense \textbf{n}on-\textbf{a}queous \textbf{p}hase \textbf{l}iquid and is located in the folder \texttt{dumux/material/components/}. Try to include the file and select the component as the non-wetting phase via the property system. \item DNAPL: A standard set of chemical substances, such as Oil and Brine, is already included in the problem (via a list of \texttt{\#include ..} statements) and hence easily accessible by default. However, this is not the case for the class \texttt{Dumux::SimpleDNAPL}, which describes a simple \textbf{d}ense \textbf{n}on-\textbf{a}queous \textbf{p}hase \textbf{l}iquid and is located in the folder \texttt{dumux/material/components/}. Try to include the file and select the component as the non-wetting phase via the property system.
\end{enumerate} \end{enumerate}
If you want to take a closer look on how the fluid classes are defined and which substances are already available please browse through the files in the directory If you want to take a closer look on how the fluid classes are defined and which substances are already available please browse through the files in the directory
@ -377,7 +379,7 @@ If you want to take a closer look on how the fluid classes are defined and which
\item \textbf{Use a Full-Fledged Fluid System} \\ \item \textbf{Use a Full-Fledged Fluid System} \\
\Dumux usually describes fluid mixtures via \textit{fluid systems}, see also chapter \ref{sec:fluidframework}. In order to include a fluid system, you first have to comment out lines \ref{tutorial-coupled:2p-system-start} to \ref{tutorial-coupled:2p-system-end} in the problem file. If you use eclipse, this can easily be done by pressing \textit{Ctrl + Shift + 7} -- the same as to cancel the comment later on.\\ \Dumux usually describes fluid mixtures via \textit{fluid systems}, see also chapter \ref{sec:fluidframework}. In order to include a fluid system, you first have to comment out lines \ref{tutorial-coupled:2p-system-start} to \ref{tutorial-coupled:2p-system-end} in the problem file. If you use eclipse, this can easily be done by pressing \textit{Ctrl + Shift + 7} -- the same as to cancel the comment later on.\\
Now include the file \texttt{fluidsystems/h2oairsystem.hh} in the material folder, and set a property \texttt{FluidSystem} with the appropriate type, i.e. \texttt{Dumux::H2OAirFluidSystem<TypeTag>}. However, this is a rather complicated fluid system considers mixtures of components and also uses tabulated components that need to be initialized -- i.e. the tables need to be filled with values. Initializating the fluid system is normally done in the constructor of the problem by calling \texttt{GET\_PROP\_TYPE(TypeTag, FluidSystem)::init();}. As water flow replacing a gas is much faster, test your simulation only until $2000$ seconds and start with a time step of $1$ second.\\ Now include the file \texttt{fluidsystems/h2oairfluidsystem.hh} in the material folder, and set a property \texttt{FluidSystem} with the appropriate type, i.e. \texttt{Dumux::H2OAirFluidSystem<TypeTag>}. However, this is a rather complicated fluid system which considers mixtures of components and also uses tabulated components that need to be initialized -- i.e. the tables need to be filled with values. The initialization of the fluid system is normally done in the constructor of the problem by calling \texttt{GET\_PROP\_TYPE(TypeTag, FluidSystem)::init();}. As water flow replacing a gas is much faster, test your simulation only until $2000$ seconds and start with a time step of $1$ second.\\
Please reverse the changes made in this part of the exercise, as we will continue to use immiscible phases from here on and hence do not need a complex fluid system. Please reverse the changes made in this part of the exercise, as we will continue to use immiscible phases from here on and hence do not need a complex fluid system.
\item \textbf{Changing Constitutive Relations} \\ \item \textbf{Changing Constitutive Relations} \\
@ -420,21 +422,23 @@ the file \texttt{tutorialproblem\_coupled.hh} (e.g. with the name
\texttt{ex2\_tutorialproblem\_coupled.hh} and new spatial parameters \texttt{ex2\_tutorialproblem\_coupled.hh} and new spatial parameters
just like \texttt{tutorialspatialparameters\_coupled.hh}. The new just like \texttt{tutorialspatialparameters\_coupled.hh}. The new
problem file needs to problem file needs to
be included in the file \texttt{tutorial\_coupled.cc}.\\ be included in the file \texttt{tutorial\_coupled.cc}.
The new files should contain the definition of new classes with names The new files should contain the definition of new classes with names
that relate to the file name, such as that relate to the file name, such as
\texttt{Ex2TutorialProblemCoupled}. Make sure that you also adjust the \texttt{Ex2TutorialProblemCoupled}. Make sure that you also adjust the
guardian macros in lines \ref{tutorial-coupled:guardian1} and guardian macros in lines \ref{tutorial-coupled:guardian1} and
\ref{tutorial-coupled:guardian1} \ref{tutorial-coupled:guardian1}
in the header files (e.g. change \\ in the header files (e.g. change
\texttt{DUMUX\_TUTORIALPROBLEM\_COUPLED\_HH} to \mbox{\texttt{DUMUX\_TUTORIALPROBLEM\_COUPLED\_HH}} to\\
\texttt{DUMUX\_EX2\_TUTORIALPROBLEM\_COUPLED\_HH}). Besides also \mbox{\texttt{DUMUX\_EX2\_TUTORIALPROBLEM\_COUPLED\_HH}}).
adjusting the guardian macros, the new problem file should define and Besides adjusting the guardian macros, the new problem file should define and
use a new type tag for the problem as well as a new problem class use a new type tag for the problem as well as a new problem class
e.g. \texttt{Ex2TutorialProblemCoupled}. Make sure to assign your e.g. \mbox{\texttt{Ex2TutorialProblemCoupled}}. Make sure to assign your
newly defined spatial parameter class to the newly defined spatial parameter class to the
\texttt{SpatialParameters} property for the new \texttt{SpatialParameters} property for the new
type tag. \\ type tag.
After this, change the run-time parameters so that they match the After this, change the run-time parameters so that they match the
domain described by figure \ref{tutorial-coupled:ex2_Domain}. Adapt domain described by figure \ref{tutorial-coupled:ex2_Domain}. Adapt
the problem class so that the boundary conditions are consistent with the problem class so that the boundary conditions are consistent with
@ -453,7 +457,7 @@ compile the program.
\begin{figure}[ht] \begin{figure}[ht]
\psfrag{K1}{K $= 10^{-7}\;\text{m}^2$} \psfrag{K1}{K $= 10^{-7}\;\text{m}^2$}
\psfrag{phi1}{$\phi = 0.2$} \psfrag{phi1}{$\phi = 0.2$}
\psfrag{Lin}{Brooks-Corey Law} \psfrag{Lin}{\textsc{Brooks}-\textsc{Corey} Law}
\psfrag{Lin2}{$\lambda = 1.8$, $p_e = 1000\;\text{Pa}$} \psfrag{Lin2}{$\lambda = 1.8$, $p_e = 1000\;\text{Pa}$}
\psfrag{K2}{K $= 10^{-9}\;\text{m}^2$} \psfrag{K2}{K $= 10^{-9}\;\text{m}^2$}
\psfrag{phi2}{$\phi = 0.15$} \psfrag{phi2}{$\phi = 0.15$}
@ -490,7 +494,7 @@ compile the program.
As you have experienced, compilation takes quite some time. Therefore, As you have experienced, compilation takes quite some time. Therefore,
\Dumux provides a simple method to read in parameters at run-time \Dumux provides a simple method to read in parameters at run-time
via \textit{paramter input files}.\\ via \textit{parameter input files}.\\
In the code, parameters can be read via the macro In the code, parameters can be read via the macro
\texttt{GET\_RUNTIME\_PARAM(TypeTag, Scalar, \texttt{GET\_RUNTIME\_PARAM(TypeTag, Scalar,
@ -505,7 +509,7 @@ happens if they are modified. For more information about the input file please r
Create a new file for the benzene component called \texttt{benzene.hh} Create a new file for the benzene component called \texttt{benzene.hh}
and implement a new component. (You may get a hint by looking at and implement a new component. (You may get a hint by looking at
existing fluid systems in the directory \verb+/dumux/material/components+). \\ existing components in the directory \verb+/dumux/material/components+). \\
Use benzene as a new fluid and run the model of Exercise 2 with water Use benzene as a new fluid and run the model of Exercise 2 with water
and benzene. Benzene has a density of $889.51 \, \text{kg} / and benzene. Benzene has a density of $889.51 \, \text{kg} /
\text{m}^3$ and a viscosity of $0.00112 \, \text{Pa} \; \text{s}$. \text{m}^3$ and a viscosity of $0.00112 \, \text{Pa} \; \text{s}$.

1
examples/tutorialproblem_coupled.hh
View File

@ -34,6 +34,7 @@
// The DUNE grid used // The DUNE grid used
#include <dune/grid/yaspgrid.hh> #include <dune/grid/yaspgrid.hh>
// Spatially dependent parameters
#include "tutorialspatialparameters_coupled.hh" #include "tutorialspatialparameters_coupled.hh"
// The components that are used // The components that are used

13
examples/tutorialspatialparameters_coupled.hh
View File

@ -46,7 +46,8 @@ namespace Properties
NEW_TYPE_TAG(TutorialSpatialParametersCoupled);/*@\label{tutorial-coupled:define-spatialparameters-typetag}@*/ NEW_TYPE_TAG(TutorialSpatialParametersCoupled);/*@\label{tutorial-coupled:define-spatialparameters-typetag}@*/
// Set the spatial parameters // Set the spatial parameters
SET_TYPE_PROP(TutorialSpatialParametersCoupled, SpatialParameters, Dumux::TutorialSpatialParametersCoupled<TypeTag>); /*@\label{tutorial-coupled:set-spatialparameters}@*/ SET_TYPE_PROP(TutorialSpatialParametersCoupled, SpatialParameters,
Dumux::TutorialSpatialParametersCoupled<TypeTag>); /*@\label{tutorial-coupled:set-spatialparameters}@*/
// Set the material law // Set the material law
SET_PROP(TutorialSpatialParametersCoupled, MaterialLaw) SET_PROP(TutorialSpatialParametersCoupled, MaterialLaw)
@ -114,8 +115,9 @@ public:
* \param fvElemGeom The finite-volume geometry in the box scheme * \param fvElemGeom The finite-volume geometry in the box scheme
* \param scvIdx The local vertex index * \param scvIdx The local vertex index
* *
* Alternatively, the function porosityAtPos(const GlobalPosition& globalPos) could be defined, * Alternatively, the function porosityAtPos(const GlobalPosition& globalPos)
* where globalPos is the vector including the global coordinates of the finite volume. * could be defined, where globalPos is the vector including the global coordinates
* of the finite volume.
*/ */
Scalar porosity(const Element &element, /*@\label{tutorial-coupled:porosity}@*/ Scalar porosity(const Element &element, /*@\label{tutorial-coupled:porosity}@*/
const FVElementGeometry &fvElemGeom, const FVElementGeometry &fvElemGeom,
@ -129,8 +131,9 @@ public:
* \param fvElemGeom The finite-volume geometry in the box scheme * \param fvElemGeom The finite-volume geometry in the box scheme
* \param scvIdx The local vertex index * \param scvIdx The local vertex index
* *
* Alternatively, the function materialLawParamsAtPos(const GlobalPosition& globalPos) could be defined, * Alternatively, the function materialLawParamsAtPos(const GlobalPosition& globalPos)
* where globalPos is the vector including the global coordinates of the finite volume. * could be defined, where globalPos is the vector including the global coordinates
* of the finite volume.
*/ */
const MaterialLawParams& materialLawParams(const Element &element, /*@\label{tutorial-coupled:matLawParams}@*/ const MaterialLawParams& materialLawParams(const Element &element, /*@\label{tutorial-coupled:matLawParams}@*/
const FVElementGeometry &fvElemGeom, const FVElementGeometry &fvElemGeom,