setup.compaction.1D {ReacTran} | R Documentation |
This function calculates the advective velocities of the pore water and the solid phase in a sediment based on the assumption of steady state compaction.
The velocities of the pore water (u
) and the solid phase (v
)
are calculated in the middle (mid
) of the grid cells and the
interfaces (int
).
One needs to specify the porosity at the interface (por.0
), the
porosity at infinite depth (por.inf
), the porosity profile
(por.grid
) encoded as a 1D grid property
(see setup.prop.1D
, as well as the advective
velocity of the solid phase at one particular depth (either at the sediment
water interface (v.0
) or at infinite depth (v.inf
)).
setup.compaction.1D(v.0 = NULL, v.inf = NULL, por.0, por.inf, por.grid)
v.0 |
advective velocity of the solid phase at the sediment-water
interface (also referred to as the sedimentation velocity); if NULL
then v.inf must not be NULL [L/T]
|
v.inf |
advective velocity of the solid phase at infinite depth
(also referred to as the burial velocity); if NULL then
v.0 must not be NULL [L/T]
|
por.0 |
porosity at the sediment-water interface |
por.inf |
porosity at infinite depth |
por.grid |
porosity profile specified as a 1D grid property
(see setup.prop.1D for details on the structure of this list)
|
A list containing:
u |
list with pore water advective velocities at the middle of the
grid cells (mid ) and at the grid cell interfaces (int ).
|
v |
list with solid phase advective velocities at the middle of the
grid cells (mid ) and at the grid cell interfaces (int ).
|
Filip Meysman <f.meysman@nioo.knaw.nl>, Karline Soetaert <k.soetaert@nioo.knaw.nl>
Meysman, F. J. R., Boudreau, B. P., Middelburg, J. J. (2005) Modeling Reactive Transport in Sediments Subject to Bioturbation and Compaction. Geochimica Et Cosmochimica Acta 69, 3601-3617
# setup of the 1D grid L <-10 grid <- setup.grid.1D(x.up=0,L=L,N=20) # attaching an exponential porosity profile to the 1D grid # this uses the "p.exp" profile function por.grid <- setup.prop.1D(func=p.exp,grid=grid,y.0=0.9,y.inf=0.5,x.att=3) # calculate the advective velocities dummy <- setup.compaction.1D(v.0=1, por.0=0.9, por.inf=0.5, por.grid=por.grid) u.grid <-dummy$u v.grid <-dummy$v # plotting the results par(mfrow=c(2,1),cex.main=1.2,cex.lab=1.2) matplot(por.grid$int,grid$x.int,pch=19,ylim=c(L,0), xlim=c(0,1), xlab="",ylab="depth [cm]",main=expression("porosity"),axes=FALSE) abline(h = 0) axis(pos=0, side=2) axis(pos=0, side=3) matplot(u.grid$int,grid$x.int,type="l",lwd=2,col="blue",ylim=c(L,0), xlim=c(0,max(u.grid$int,v.grid$int)), xlab="",ylab="depth [cm]",main=expression("advective velocity [cm yr-1]"), axes=FALSE) abline(h = 0) axis(pos=0, side=2) axis(pos=0, side=3) lines(v.grid$int,grid$x.int,lwd="2",col="red") legend(x="bottomright", legend=c("pore water","solid phase"), col=c("blue","red"),lwd=c(2,2))