RMkolmogorov: Identical Model

Description Usage Arguments Value References See Also Examples

View source: R/RMmodels.R

Description

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

γ_{ij}(h) = |h|^{2/3}(4/3 * δ_{ij} - 1/3 * h_i h_j / \|h\|^2)

Usage

1

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.

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

1
2
3
4
5
6
7
8
9
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"))

RandomFields documentation built on Jan. 19, 2022, 1:06 a.m.