Skip to content
Snippets Groups Projects
Commit a664f7b6 authored by Felix Weinhardt's avatar Felix Weinhardt Committed by Katharina Heck
Browse files

added plots of the results and added to the README.md

parent d1e2361e
No related branches found
No related tags found
1 merge request!1547[examples] Set up problem for documentation (2p, adaptive, point source)
examples/1ptracer/Plots/Tracer_2500.png

14.2 KiB

examples/1ptracer/Plots/Tracer_5000.png

15.7 KiB

examples/1ptracer/Plots/VelocityProfile.png

24.2 KiB

This tutorial was copied from dumux/test/porousmediumflow/tracer/1ptracer.
# One-phase flow with with random permeability distribution (and a tracer model)
# One-phase flow with with random permeability distribution and a tracer model
## Problem set-up
A contaminant tracer is diluted by diffusion and a base groundwater flow from the bottom to the top. The permeability within the domain is randomly distributed.
### main file
This example combines a stationary One-phase flow problem with a tracer model. In the first step, the groundwater-velocity is evaluated under stationary conditions. Based on the volume fluxes, the tracer model is solved instationary. Therefore both, the problem_1p.hh and the problem_tracer.hh have to be included in the main file.
## Tracer model
This example combines a stationary One-phase flow problem with a tracer model. In the first step, the groundwater-velocity is evaluated under stationary conditions. Based on the volume fluxes, the tracer model is solved instationary.
### main.cc
Both, the problem_1p.hh and the problem_tracer.hh have to be included in the main file.
```C++
#include "problem_1p.hh"
......@@ -15,11 +18,14 @@ This example combines a stationary One-phase flow problem with a tracer model. I
In lines 83-192, the stationary 1p problem is setup and soved and the volume fluxes are calculated for the tracer problem.
In lines 193-292, the tracer problem is set up on the same grid and solved instationary.
## 1p problem
### problem_1p.hh
??
### problem_tracer.hh
The molar mass of the component and the binary diffusion coefficient can modified in lines 117-128.
### spatialparams_tracer.hh
The density, the molar mass of the fluid as well as the prosity and dispersivity can be adapted in the lines 63-93.
## Random permeability generation
### spatialparams_1p.hh
#### Random permeability generation
The follwing code can be found in lines 64-72.
```C++
......@@ -36,10 +42,12 @@ The follwing code can be found in lines 64-72.
```
Two lognormal distributions are defined: K and KLens with different means and standard deviations. In this case the mean is defined as the log of the permeability given from the input file and the standard deviation is 10% of the mean.
Within a loop through all elements the randomly generated permeabilities are assigned according to their position in the domain (inside or outside the lense).
### params.input
The two mean values of the permeability are specified in the params.input file.
```C++
Permeability = 1e-10 # [m^2]
PermeabilityLens = 1e-11 # [m^2]
```
## Results
## Tracer model
### problem_tracer.hh
The molar mass of the component and the binary diffusion coefficient can modified in lines 117-128.
### spatialparams_tracer.hh
The density, the molar mass of the fluid as well as the prosity and dispersivity can be adapted in the lines 63-93.
\ No newline at end of file
![](/Plots/VelocityProfile.png) ![](/Plots/Tracer_2500.png) ![](/Plots/Tracer_5000.png)
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment