Wind turbulence in the neutral and unstable surface-layer

Version 1.4 (4.45 MB) by E. Cheynet
A short set of functions to study surface layer fluxes and wind velocity spectra (simulated data) for wind engineering applications.
311 Downloads
Updated 31 Oct 2019

View License

A meteorological tower with 3 sonic anemometers at 3 different heights is "simulated" and 4 samples corresponding to both neutral and unstable atmospheric stratifications are provided for each height.
The goal is to use the eddy-covariance method to retrieve the different turbulence statistics and the associated flux of momentum and heat in the surface layer.
The submission contains

• 4 .mat files with the different samples (simulated time series)

• One example file Example.mlx

• The functions “verticalSpectrum_Hojstrup” and “longitudinalSpectrum_Hojstrup” to compute the spectral model from Højstrup [1].

• The function “getStability” to estimate the Obukhov length, the heat flux, the scaling temperature and friction velocity.

• The function “obukhovLength” to compute the Obukhov length.

• The function “similarityFun” to get the non-dimensional vertical profile of heat and momentum.

• The functions “plotSimilarityFun” and “plotVelocitySpectra” to plot more clearly the velocity spectra and the non-dimensional profiles

Any comments or suggestion to improve this submission is warmly welcomed.

[1] Højstrup, J. (1981). A simple model for the adjustment of velocity spectra in unstable conditions downstream of an abrupt change in roughness and heat flux. Boundary-Layer Meteorology, 21(3), 341-356.

Cite As

E. Cheynet (2024). Wind turbulence in the neutral and unstable surface-layer (https://www.mathworks.com/matlabcentral/fileexchange/63355-wind-turbulence-in-the-neutral-and-unstable-surface-layer), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2018b
Compatible with any release
Platform Compatibility
Windows macOS Linux
Categories
Find more on Weather and Atmospheric Science 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
1.4

Example updated

1.3

The example file has been converted into an interactive example

1.2

Added project website

1.1

Correct version loaded
The possibility to include the across-wind component "v" is added
Some typo are corrected + the description is improved

1.0.0.0

Title modified
Description added
typo