adVeloc: Advection velocities in porous medium

In dkneis/diatools: Utilities for Diagenetic Modeling

Description

Yields advection velocities of a solid or liquid component in a porous medium to fulfil the constraint of mass conservation (continuity eqn.).

Usage

adVeloc(uRef, pRef, p, solid)

Arguments

uRef	Reference advection velocity (e.g. at a certain depth).
pRef	Porosity corresponding to uRef.
p	Vector of porosities for which advection velocities are to be returned.
solid	Logical. Should velocities be returned for solids or liquid?

Value

Vector of the same length as p.

Note

The underlying principle of continuity is also known as steady-state compaction; see ReacTran::setup.compaction.1D.

Author(s)

David Kneis [e-mail]

Examples

z= 0:10
p= seq(0.99, 0.95, length.out=length(z))
plot(z, adVeloc(uRef=1, pRef=p[length(p)], p=p, solid=TRUE), type="l")
lines(z, adVeloc(uRef=1, pRef=p[length(p)], p=p, solid=FALSE), lty=2)
legend("top", bty="n", lty=1:2, legend=c("solid","liquid"))

dkneis/diatools documentation built on May 15, 2019, 9:12 a.m.