File Exchange

image thumbnail

Wind field simulation (the fast version)

version 1.0.8 (8.71 MB) by E. Cheynet
A three-variate turbulent wind field (u,v and w components) is simulated in three-dimensions.


Updated 03 Mar 2019

View License

The present submission deals with the simulation of turbulent wind field (u,v,w, components) in 3-D (two dimensions for space and one for the time). The computational efficiency of the simulation relies on Ref. [1], which leads to a significantly shorter simulation time than the function windSim, also available on fileExchange. However, only the case of a regular 2D vertical grid normal to the flow is here considered.

The submission contains:
- An example file Example1 that illustrates simply how the output variables look like.
- An example file Example2, which is more complete, and which simulates a 3-D turbulent wind field on a 7x7 grid.
- A data file exampleData.mat used in Example1.
- The function windSimFast.m, which is used to generate the turbulent wind field. A similar implementation of windSimFast.m was used in ref. [2].
- The function getSamplingpara.m, which computes the time and frequency vectors.
- The function KaimalModel.m, which generates the one-point auto and cross-spectral densities of the velocity fluctuations, following the Kaimal model [3].


[1] Shinozuka, M., & Deodatis, G. (1991). Simulation of stochastic processes by spectral representation. Applied Mechanics Reviews, 44(4), 191-204.

[2] Wang, J., Cheynet, E., Snæbjörnsson, J. Þ., & Jakobsen, J. B. (2018). Coupled aerodynamic and hydrodynamic response of a long span bridge suspended from floating towers. Journal of Wind Engineering and Industrial Aerodynamics, 177, 19-31.

[3] Davenport, A. G. (1961). The spectrum of horizontal gustiness near the ground in high winds. Quarterly Journal of the Royal Meteorological Society, 87(372), 194-211.

Any comment, suggestion or question is welcomed.

Cite As

E. Cheynet (2019). Wind field simulation (the fast version) (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (3)


Thanks a lot, friend. I guess my model is pretty close to the limit if the ABL is over 1000 meters. The height is between 100 to 120 meters. :-)

E. Cheynet

Hei Wei,
Thank for the feedback! Yes, this simulation should be suitable for the wind-induced response of a high-rise building. You may find the appropriate wind parameters in pre-existing standards and codes. However, I am not sure about their applicability at a height greater than 10% of the atmospheric boundary layer depth.


Thanks, Cheynet.
With your simulation code, I had some clue on how to do the wind simulation now. Thank you so much.
By the way, friend. Is this wind simulation suitable to be applied to achieve the response of high-rise buildings?



Added project website


A minor correction for the reconstruction of the fft vector has been applied




Screenshot updated


Summary updated




Updated description



MATLAB Release Compatibility
Created with R2018a
Compatible with any release
Platform Compatibility
Windows macOS Linux