RMkolmogorov: Identical Model
In RandomFields: Simulation and Analysis of Random Fields

Description Usage Arguments Value References See Also Examples

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)

1	RMkolmogorov(var, scale, Aniso, proj)

var,scale,Aniso,proj

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

RMkolmogorov returns an object of class RMmodel.

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

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

RMmodel, RMcurlfree, RMdivfree, RMvector.

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