Learn R Programming

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

Run this code
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"))

Run the code above in your browser using DataLab