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

handbook: a few more changes to the coupled tutorial

Browse Source
This commit is contained in:
Andreas Lauser 2012-02-13 19:16:29 +00:00 committed by Andreas Lauser
parent 25aba6a3f2
commit 0c3e573194
1 changed files with 11 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -313,7 +313,7 @@ access / set functions as well as the rest of the implementation of
the material description can be found in the material description can be found in
\verb+dumux/material/fluidmatrixinteractions/2p+. \verb+dumux/material/fluidmatrixinteractions/2p+.
\subsection{The definition of run-time parameters}\label{tutorial-coupled:runtime-parameters} \subsection{Definition of Run-Time Parameters}\label{tutorial-coupled:runtime-parameters}
Some parameters need to be specified at runtime. These can either be Some parameters need to be specified at runtime. These can either be
specified directly via command line optioms or the can be collected in specified directly via command line optioms or the can be collected in
@ -339,7 +339,7 @@ to make only some small changes in the tutorial files.
\begin{enumerate} \begin{enumerate}
\item \textbf{Runing the Program} \\ \item \textbf{Runing the Program} \\
To get an impression what the results should look like you can first run the original version of the coupled tutorial model by typing \texttt{./tutorial\_coupled}. Note, that the time-step size is automatically adapted during the simulation. For the visualization of the results using paraview please refer to \ref{quick-start-guide}.\\ To get an impression what the results should look like you can first run the original version of the coupled tutorial model by typing \texttt{./tutorial\_coupled}. Note, that the time-step size is automatically adapted during the simulation. For the visualization of the results using paraview please refer to section \ref{quick-start-guide}.\\
\item \textbf{Changing the Model Domain and the Boundary Conditions} \\ \item \textbf{Changing the Model Domain and the Boundary Conditions} \\
Change the size of the model domain so that you get a rectangle with Change the size of the model domain so that you get a rectangle with
@ -364,13 +364,13 @@ Now you can change the fluids. Use DNAPL instead of Oil and Brine instead of Wat
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
\texttt{/dumux/material/components} and read chapter~\ref{sec:fluidframework}. \texttt{/dumux/material/components} and read chapter~\ref{sec:fluidframework}.
\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/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.\\
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} \\
Use an unregularized linear law with an entry pressure of $p_e = 0.0$ and maximal capillary pressure of e.g. $p_{c_{max}} = 2000.0$ instead of using a Use an unregularized linear law with an entry pressure of $p_e = 0.0\;\text{Pa}$ and maximal capillary pressure of e.g. $p_{c_{max}} = 2000.0\;\text{Pa}$ instead of using a
regularized Brooks-Corey law for the regularized Brooks-Corey law for the
relative permeability and for the capillary pressure saturation relationship. To do that you have relative permeability and for the capillary pressure saturation relationship. To do that you have
to change the file \texttt{tutorialspatialparameters\_coupled.hh}. to change the file \texttt{tutorialspatialparameters\_coupled.hh}.
@ -385,20 +385,20 @@ of the linear law and the respective \texttt{set}-functions can be found
\ref{tutorial-coupled:exercise1_d}. Adjust the boundary conditions \ref{tutorial-coupled:exercise1_d}. Adjust the boundary conditions
so that water is again flowing from the left to the right of the so that water is again flowing from the left to the right of the
\begin{figure}[ht] \begin{figure}[ht]
\psfrag{K1 =}{K $= 10^{-8}\text{ m}^2$} \psfrag{K1 =}{$\mathbf{K} = 10^{-8}\;\text{m}^2$}
\psfrag{phi1 =}{$\phi = 0.15$} \psfrag{phi1 =}{$\phi = 0.15$}
\psfrag{K2 =}{\textcolor{white}{K $= 10^{-9}\text{ m}^2$}} \psfrag{K2 =}{\textcolor{white}{$\mathbf{K} = 10^{-9}\;\text{m}^2$}}
\psfrag{phi2 =}{\textcolor{white}{$\phi = 0.3$}} \psfrag{phi2 =}{\textcolor{white}{$\phi = 0.3$}}
\psfrag{600 m}{600 m} \psfrag{600 m}{$600 \;\text{m}$}
\psfrag{300 m}{300 m} \psfrag{300 m}{$300 \;\text{m}$}
\centering \centering
\includegraphics[width=0.5\linewidth,keepaspectratio]{EPS/exercise1_c.eps} \includegraphics[width=0.5\linewidth,keepaspectratio]{EPS/exercise1_c.eps}
\caption{Exercise 1f: Set-up of a model domain with a heterogeneity. $\Delta \text{x} = 20$ m $\Delta \text{y} = 20$ m.}\label{tutorial-coupled:exercise1_d} \caption{Exercise 1f: Set-up of a model domain with a heterogeneity. $\Delta x = 20 \;\text{m}$ $\Delta y = 20\;\text{m}$.}\label{tutorial-coupled:exercise1_d}
\end{figure} \end{figure}
domain. You can use the fluids of exercise 1c).\\ domain. You can use the fluids of exercise 1c).\\
\textbf{Hint:} The current position of the control volume can be obtained via \texttt{element.geometry().corner(scvIdx);}.\\ \textbf{Hint:} The current position of the control volume can be obtained using \texttt{element\allowbreak.geometry()\allowbreak.corner(scvIdx)}.\\
When does the front cross the material border? In paraview, the When does the front cross the material border? In paraview, the
animation inspector (\textit{View} $\rightarrow$ \textit{Animation animation view (\textit{View} $\rightarrow$ \textit{Animation
View}) is a convenient way to get a rough feeling of the time-step View}) is a convenient way to get a rough feeling of the time-step
sizes. sizes.
\end{enumerate} \end{enumerate}