From 7bcdcdb9455c3fb5433a9520a646017d1396bfe2 Mon Sep 17 00:00:00 2001 From: Leon Keim <leon.keim@googlemail.com> Date: Fri, 31 Mar 2023 18:33:55 +0200 Subject: [PATCH] specify the execise README for possible non multiphase experts --- exercises/exercise-fluidsystem/README.md | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/exercises/exercise-fluidsystem/README.md b/exercises/exercise-fluidsystem/README.md index cf206580..3546dd6f 100644 --- a/exercises/exercise-fluidsystem/README.md +++ b/exercises/exercise-fluidsystem/README.md @@ -1,6 +1,6 @@ # Exercise Fluidsystem (DuMuX Course) -The aim of this exercise is to get familiar with the _DuMu<sup>x</sup>_ way of implementing new components (fluids) and fluid systems (mixtures). In the scope of this exercise, a new fictitious component is implemented (Section 2) as well as its mixture with water (Section 3). +The aim of this exercise is to get familiar with the _DuMu<sup>x</sup>_ way of implementing new components and using them in fluid systems (immiscible/mixture). In the scope of this exercise, different two-phase (liquid-liquid) fluid systems are implemented. A new fictitious component is implemented and used to describe a immiscible liquid phase. The second phase will only consist of water (Section 2). Furthermore, the mixture of the two components is implemented as two miscible phases (Section 3). ## Problem set-up @@ -63,14 +63,14 @@ while the `BoxModel` _TypeTag_ can be found in the `discretization/box.hh` heade For a cell-centered scheme, you could derive from `CCTpfaModel` or `CCMpfaModel` instead (and, of course, include the right headers). -As one of the two phases we want to use water and we want to precompute tables on which the properties are then interpolated in order to save computational time. Since we need this part in both problem file and properties file, we have to include the following two headers in our problem file, i.e. `2pproblem.hh` file, and the properties file has access to them through problem file. +One of the two phases should only contain the component water. We want to precompute tables on which the properties are then interpolated in order to save computational time. Since we need this part in both problem file and properties file, we have to include the following two headers in our problem file, i.e. `2pproblem.hh` file, and the properties file has access to them through problem file. ```c++ // The water component #include <dumux/material/components/tabulatedcomponent.hh> #include <dumux/material/components/h2o.hh> ``` -The other phase will be created from our new component, where we want to implement an incompressible and a compressible variant. +The other phase will be created only contain our new component, where we want to implement an incompressible and a compressible variant. The respective headers are prepared, but still incomplete. The compressible variant is still commented so that compilation does not fail when finishing the incompressible variant. ```c++ @@ -78,7 +78,7 @@ The respective headers are prepared, but still incomplete. The compressible vari #include "components/myincompressiblecomponent.hh" // #include "components/mycompressiblecomponent.hh" ``` -As mentioned above, we want to simulate two non-mixing components. The respective fluid system is found in: +As mentioned above, we want to simulate two non-mixing phases (2pimmiscible). The respective fluid system is found in: ```c++ // The two-phase immiscible fluid system @@ -93,7 +93,7 @@ This fluid system expects __phases__ as input and so far we have only included t ``` which creates a _liquid phase_ from a given component. Finally, using all of the included classes we set the fluid system property by choosing that the water phase is liquid (`OnePLiquid`) and consists of the tabulated water component, and -the other phase is liquid as well and consists of the incompressible fictitious component. Both will make up the immiscible fluid system (`TwoPImmiscible`): +the other phase is liquid as well and consists of the incompressible fictitious component. Both will make up the immiscible fluid system (`TwoPImmiscible`), consiting of two liquids each consisting of one component: ```c++ @@ -176,7 +176,7 @@ Compare the gnuplot output to the following plot of the density function from ab ### 3. Implement a new fluid system The problem file and properties file for this part of the exercise are `2p2cproblem.hh` and `2p2cproperties.hh`, respectively. -We now want to implement a new fluid system consisting of two liquid phases,which are water and the previously implemented compressible component. We will consider compositional effects, which is why we now have to derive our _TypeTag_ (`ExerciseFluidsystemTwoPTwoC`) from a _TypeTag_ (`TwoPTwoC`) that holds the miscible two-phase +We now want to implement a new fluid system, which still consists of two liquid phases. However, one phase consists of mainly of water and the other consists mainly of the previously implemented compressible component. We will now consider compositional effects, which is why we have to derive our _TypeTag_ (`ExerciseFluidsystemTwoPTwoC`) from a _TypeTag_ (`TwoPTwoC`) that holds the miscible two-phase two-component model properties: ```c++ @@ -222,7 +222,7 @@ In the `fluidsystems/h2omycompressiblecomponent.hh` file, your implemented compr __Task__: -Under the assumption that one molecule of `MyCompressibleComponent` displaces exactly one molecule of water, the water phase density can be expressed as follows: +Under the assumption that one molecule of `MyCompressibleComponent` displaces exactly one molecule of water, the density of the mostly-water-phase can be expressed as follows: $` \rho_{w} = \frac{ \rho_{w, pure} }{ M_{H_2O} }*(M_{H_2O}*x_{H_2O} + M_{MyComponent}*x_{MyComponent}) `$ @@ -257,7 +257,7 @@ int wettingPhaseAtPos(const GlobalPosition& globalPos) const } ``` -Change this function such that the phase of our new component is the wetting phase __only__ within the lenses. Execute the program of task 3 again: +Change this function such that the phase of our mostly-new-component-phase is the wetting phase __only__ within the lenses. Execute the program of task 3 again: ```bash cd build-cmake/exercises/exercise-fluidsystem -- GitLab