Commit b78a1ff1 authored by Gabi Seitz's avatar Gabi Seitz
Browse files

Merge branch 'feature/biomin-readme' into 'master'

Feature/biomin readme

See merge request !19
parents b429c869 8ea575b9
......@@ -9,4 +9,4 @@ Click on the exercise to go the description
* [Exercise on properties](./exercise-properties/README.md)
* [Exercise on fluid systems](./exercise-fluidsystem/README.md)
* [Exercise on how to create a new Dune module](./exercise-dunemodule/README.md)
* [Exercise on biomineralization (SFB PA C)](./exercise-biomineralization/README.md)
# Exercise #4 (DuMuX course)
# Exercise Biomineralization
The aim of this exercise is to a first glimpse with the _DuMuX_ way of implementing mineralization and reaction processes. In the scope of this exercise, the setting of boundary conditions is revisited and a new reaction term is implemented.
The aim of this exercise is to get a first glimpse with the _DuMuX_ way of implementing mineralization and reaction processes. In the scope of this exercise, the setting of boundary conditions is revisited and a new reaction term is implemented.
## Problem set-up
The domain has a size of 20 x 15 m and contains a sealing aquitard in the middle. The aquitard is interrupted by a "fault zone" and thereby connects the upper (drinking water) aquifer and the lower CO2-storage aquifer. Initially, the domain is fully water saturated and biofilm is present in the lower CO2-storage aquifer. Calcium and urea are injected in the upper drinking water aquifer by means of a Neumann boundary condition. The remaining parts of the upper and the entire lower boundary are Neumann no-flow while on the right side a Dirichlet boundary conditions is applied (according to the initial values).
The domain has a size of 20 x 15 m and contains a sealing aquitard in the middle. The aquitard is interrupted by a "fault zone" and thereby connects the upper drinking water aquifer and the lower CO2-storage aquifer. Initially, the domain is fully water saturated and biofilm is present in the lower CO2-storage aquifer. Calcium and urea are injected in the upper drinking water aquifer by means of a Neumann boundary condition. The remaining parts of the upper and the entire lower boundary are Neumann no-flow while on the right side a Dirichlet boundary conditions is applied according to the initial values.
Disclaimer: Please note, that this is not a realistic scenario. One does not want to store gaseous CO2 in such a subcritical setting.
......@@ -33,11 +33,11 @@ Furthermore you will find the following folders:
To see more chemistry, components, fluidsystems and solidsystems implementations, have a look at the folder `dumux/material`.
__Special note on solidsystems:__
There are two types of solid components. Reactive and inert. For each reactive component one mass balances is solved. The inert components compose the "unchanging" (inert) rock matrix.
There are two types of solid components. Reactive and inert. For each reactive component one mass balance is solved. The inert components compose the "unchanging" (inert) rock matrix.
### 2. Implement a chemical equation
In the following the basic steps required to set the new chemical equation are outlined. Here, this is done in the __chemistry__ folder in the prepared file: 'simplebiominreactions.hh' starting in line 96.
In the following the basic steps required to set the new chemical equation are outlined. Here, this is done in the __chemistry__ folder in the prepared file: `simplebiominreactions.hh` starting in line 96.
Please be aware, that the chemistry file already provides some convenience functions (e.g. ``moleFracToMolality()``).
__Task__
......@@ -58,11 +58,11 @@ $`\displaystyle Z_{urease,biofilm} = k_{urease,biofilm} * mass_{biofilm}`$
The last step is defining the source term for each component according to the chemical equation:
$`\displaystyle Ca^{2+} + CO(NH)_{2} + 2 H_{2}O --> 2 NH_{4}^{+} + CaCO_{3}`$
$`\displaystyle Ca^{2+} + CO(NH)_{2} + 2 H_{2}O -> 2 NH_{4}^{+} + CaCO_{3}`$
which is:
Calcium ion + Urea + 2 Water --> 2 Ammonium ion + Calcite
Calcium ion + Urea + 2 Water → 2 Ammonium ions + Calcite
### 3. Make use of your newly created chemical equation
......@@ -119,7 +119,7 @@ Now the sealed aquitard is tested with a CO2-Injection into the lower CO2-storag
__Task:__
Implement a new boundary condition from the left boundary, injecting CO2 from 2 m to 3 m. Make sure, that the injection time for the calcium and urea is finished. You can use the predefined value `gasFlux` directly and divide it by the molar mass of CO2.
Implement a new boundary condition on the left boundary, injecting CO2 from 2 m to 3 m from the bottom. Make sure, that the injection time for the calcium and urea is finished. You can use the predefined value `gasFlux` directly and divide it by the molar mass of CO2.
Run two simulations and compare them side by side by creating two input files, or overwriting the input file in the command line:
```bash
./exercisebiomin -Problem.Name biominNoUrea -Injection.ConcUrea 0
......
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