Reactive transport modelling in 1D, 2D and 3D

Description

R-package ReacTran contains routines that enable the development of reactive transport models in aquatic systems (rivers, lakes, oceans), porous media (floc aggregates, sediments,...) and even idealized organisms (spherical cells, cylindrical worms,...).

The geometry of the model domain is either one-dimensional, two-dimensional or three-dimensional.

The package contains:

  • Functions to setup a finite-difference grid (1D or 2D)

  • Functions to attach parameters and properties to this grid (1D or 2D)

  • Functions to calculate the advective-diffusive transport term over the grid (1D, 2D, 3D)

  • Utility functions

Details

Package: ReacTran
Type: Package
Version: 1.3.2
Date: 2011-07-01
License: GNU Public License 2 or above

Author(s)

Karline Soetaert (Maintainer)

Filip Meysman

See Also

Functions ode.1D, ode.2D, ode.3D from package deSolve to integrate the reactive-transport model

Functions steady.1D, steady.2D, steady.3D from package rootSolve to find the steady-state solution of the reactive-transport model

tran.1D, tran.2D, tran.3D for a discretisation of the general transport equations

tran.volume.1D for discretisation of the 1-D transport equations using finite volumes

tran.cylindrical, tran.spherical for a discretisation of 3-D transport equations in cylindrical and spherical coordinates

tran.polar, for a discretisation of 2-D transport equations in polar coordinates

setup.grid.1D, setup.grid.2D for the creation of grids in 1-D and in 2-D

setup.prop.1D, setup.prop.2D for defining properties on these grids.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
## Not run: 

## show examples (see respective help pages for details)
## 1-dimensional transport
example(tran.1D)
example(tran.volume.1D)

## 2-dimensional transport
example(tran.2D)
example(tran.polar)

## 3-dimensional transport
example(tran.3D)
example(tran.cylindrical)
example(tran.spherical)

## open the directory with documents
browseURL(paste(system.file(package="ReacTran"), "/doc", sep=""))

## open the directory with fortran codes of the transport functions
browseURL(paste(system.file(package="ReacTran"), "/doc/fortran", sep=""))

## show package vignette with how to use ReacTran and how to solve PDEs
## + source code of the vignettes
vignette("ReacTran")
vignette("PDE")
edit(vignette("ReacTran"))

## a directory with fortran implementations of the transport
browseURL(paste(system.file(package="ReacTran"), "/doc/fortran", sep=""))

## End(Not run)