RandomFields (version 3.1.16)

RMkolmogorov: Identical Model

Description

RMkolmogorov corresponds to a vector-valued random fields 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.

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"))