setup.grid.2D: Creates a Finite Difference Grid over a Two-Dimensional... In ReacTran: Reactive Transport Modelling in 1d, 2d and 3d

Description

Creates a finite difference grid over a rectangular two-dimensional model domain starting from two separate one-dimensional grids (as created by `setup.grid.1D`).

Usage

 `1` ```setup.grid.2D(x.grid = NULL, y.grid = NULL) ```

Arguments

 `x.grid ` list containing the one-dimensional grid in the vertical direction - see `setup.grid.1D` for the structure of the list `y.grid ` list containing the one-dimensional grid in the horizontal direction - see `setup.grid.1D` for the structure of the list

Value

a list of type `grid.2D` containing:

 `x.up ` position of the upstream interface in x-direction (i.e. if x is vertical, the upper boundary); one value `x.down ` position of the downstream interface in x-direction (i.e. if x is vertical, the lower boundary); one value `x.mid ` position of the middle of the grid cells in x-direction; vector of length `x.N` `x.int ` position of the interfaces of the grid cells in x-direction; vector of length `x.N+1` `dx ` distance between adjacent cell interfaces in x-direction (thickness of grid cells); vector of length `x.N` `dx.aux ` auxiliary vector containing the distance between adjacent cell centers; at the upstream and downstream boundary calculated as (`x[1]-x.up`) and (`x.down-x[x.N]`) respectively; vector of length `x.N+1` `x.N ` total number of grid cells in the x direction; one value `y.left ` position of the upstream interface in y-direction (i.e. if y us the horizontal, the left boundary); one value `y.right ` position of the downstream interface in y-direction (i.e. if y us the horizontal, the right boundary); one value `y.mid ` position of the middle of the grid cells in y-direction; vector of length `y.N` `y.int ` position of the interfaces of the grid cells in y-direction; vector of length `y.N+1` `dy ` distance between adjacent cell interfaces in y-direction (thickness of grid cells); vector of length `y.N` `dy.aux ` auxiliary vector containing the distance between adjacent cell centers; at the upstream and downstream boundary calculated as (`y[1]-y.up`) and (`y.down-y[y.N]`) respectively; vector of length `y.N+1` `y.N ` total number of grid cells in the y direction; one value

Author(s)

Filip Meysman <[email protected]>, Karline Soetaert <[email protected]>

`tran.2D`, for a discretisation of the general transport equation in 2-D

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

`setup.prop.2D` for defining properties on the 2-D grid.

Examples

 ```1 2 3 4``` ```# test of the setup.grid.2D functionality x.grid <- setup.grid.1D(x.up = 0, L = 10, N = 5) y.grid <- setup.grid.1D(x.up = 0, L = 20, N = 10) (grid2D <- setup.grid.2D(x.grid, y.grid)) ```

ReacTran documentation built on Aug. 15, 2017, 3:01 p.m.