RMkolmogorov {RandomFields}

R Documentation

Identical Model

Description

RMkolmogorov corresponds to a vector-valued random field with covariance function

\gamma_{ij}(h) = \|h\|^{2/3}\left(\frac43 \delta_{ij} - \frac13\frac{h_ih_j}{\|h\|^2}\right)

Usage

 RMkolmogorov(var, scale, Aniso, proj)

Arguments

var, scale, Aniso, proj

optional arguments; same meaning for any RMmodel. If not passed, the above covariance function remains unmodified.

Value

RMkolmogorov returns an object of class RMmodel.

Author(s)

Martin Schlather, schlather@math.uni-mannheim.de, https://www.wim.uni-mannheim.de/schlather/

References

The above formula is from eq. (6.32) of section 6.2 in

Pope, S.B. (2000) Turbulent Flows. Cambridge: Cambridge University Pess.

See Also

RMmodel, RMcurlfree, RMdivfree, RMvector.

Examples

RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set
##                   RFoptions(seed=NA) to make them all random again

x <- y <- seq(-2, 2, len=20)
model <- RMkolmogorov()
plot(model, dim=3, MARGIN=1:2, fixed.MARGIN=1)

simu <- RFsimulate(model, x, y, z=0)
plot(simu, select.variables=list(c(1,2)), col=c("red"))

---