Wind field simulation (text-based input)

Simulation of spatially correlated wind velocity time-histories
Updated 27 Mar 2022

Wind turbulence generation using text-based input files

For a more robust and time-efficient Matlab implementation, see

View Wind field simulation (the user-friendly version) on File Exchange Donation


A method to simulate spatially correlated turbulent wind histories is implemented following [1,2]. Two possible vertical wind profiles and two possible wind spectra are implemented. The user is free to implement new ones. The wind co-coherence is a simple exponential decay as done by Davenport [3]. If the wind field is simulated in a grid, the function windSim.m should be used (cf. Examples 1 and 2). For a more complex geometry, such as a radial grid, the function windSim.m has an optional parameter to include two inputs (cf. Example3.mlx): The first one contains the wind properties, and the second one contains the coordinates of the nodes where wind histories are simulated (cf. Example 3).


The folder contains:

  • 1 input file INPUT.txt for Example1.m
  • 1 input file INPUT_MAST.txt for Example2.m
  • 2 input files windData.txt and circle.txt for Example3.m
  • The function windSim.m
  • 4 examples files Example1.m, Example2.m, Example3.m and Example4.m
  • The function coherence.m that computes the co-coherence. Notes:
  • Simulating the wind field in a high number of points with a high sampling frequency may take a lot of time.
  • This code aims to be highly customizable
  • A faster version of the present submission has been used to simulate the turbulent wind load on a floating suspension bridge [4].


[1] Shinozuka, M., Monte Carlo solution of structural dynamics, Computers and Structures, Vol. 2, 1972, pp. 855 – 874

[2] Deodatis, G., Simulation of ergodic multivariate stochastic processes, Journal of Engineering Mechanics, ASCE, Vol. 122 No. 8, 1996, pp. 778 – 787.

[3] Davenport, A. G. (1961), The spectrum of horizontal gustiness near the ground in high winds. Q.J.R. Meteorol. Soc., 87: 194–211

[4] 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.

Cite As

Cheynet, E. Wind Field Simulation (Text-Based Input). Zenodo, 2020, doi:10.5281/ZENODO.3823864.

View more styles
MATLAB Release Compatibility
Created with R2020b
Compatible with any release
Platform Compatibility
Windows macOS Linux
Find more on Numerical Integration and Differential Equations in Help Center and MATLAB Answers

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
Version Published Release Notes

See release notes for this release on GitHub:


See release notes for this release on GitHub:


See release notes for this release on GitHub:


See release notes for this release on GitHub:


See release notes for this release on GitHub:


See release notes for this release on GitHub:


See release notes for this release on GitHub:


See release notes for this release on GitHub:


See release notes for this release on GitHub:


description in the txt files






performance improved using the IFFT (same method as in windSim - the fast version)


Added the project website

The example files have been re-written as live-script

Updated the title

Major update on windSim2.m: The integral length scale in the von Karman spectral model is now a function of the height. its value is defined using the asymptotic behaviour of the spectrum at high frequencies. + New Example file.

Major update: The standard deviation of the wind fluctuations is now asked in the input files instead of the turbulence intensity. This offers a more realistic simulation, where the turbulence intensity decreases for increasing mean wind velocities.

- new Example, and updated description

summary updated
new title
updated picture
- new example included

- Updated example (need 90 sec to run instead of 2-3 sec ) but it is more complete

-description updated
- The cross-wind component v has been added !
- Example updated
- removed "2D" from title because it is not actual any more
-description updated

- update following @Paolo's feedback

- I have replaced the name of the function 'cohere' with 'coherence' since 'cohere' is a Matlab function that already exist.

-correct update now
- New definition of the vector frequency to optimize the computation time
- introduction of the co-coherence function
- The co-coherence is used instead of the magnitude squared coherence to be more rigorous
- title of example file is changed

HTML example introduced + update: one single function to run the whole simulation
-dexcription updated

-picture updated
- erratum at the end of Run_Me is corrected

- included verification.m for a preview of the output ( spectra and time series)


- A description of the output of the simulation is included at the end of the .m file "Run_Me"

- description updated

-description updated

To view or report issues in this GitHub add-on, visit the GitHub Repository.
To view or report issues in this GitHub add-on, visit the GitHub Repository.