Skip to main content
Log in

A numerical simulation of boundary-layer flows near shelterbelts

  • Published:
Boundary-Layer Meteorology Aims and scope Submit manuscript

Abstract

We have developed a shelterbelt boundary-layer numerical model to study the patterns and dynamic processes relating to flow interaction with shelterbelts. The model simulates characteristics of all three zones of airflow passing over and through shelterbelts: the windward windspeed-reduction zone, the overspeeding zone above the shelterbelt, and the leeward windspeed-reduction zone. Locations of the maximum windspeed reduction and recirculation zone, as well as the leeward windspeed-recovery rate are well simulated by the model. Where comparisons with field measurements and wind-tunnel experiments were possible, the model demonstrated good performance for flows over and through shelters ranging from almost completely open to almost solid.

The dynamic pressure resulting from the convergence and divergence of the flow field alters the perturbation pressure field. The disturbed pressure controls not only the formation of the separated flow but also the location of maximum windspeed reduction, streamline curvature, speed-up over the shelterbelt, and leeward windspeed recovery rate. The interaction of pressure with the flow produces complex flow patterns, the characteristics of which are determined, to a great extent, by shelterbelt structure.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Similar content being viewed by others

Abbreviations

<>:

Spatial averaged value

(−):

Temporal averaged value

(′):

Departure of variable from its averaged value

A :

Leaf surface-area density

ABL:

Atmospheric boundary layer

C :

Drag coefficient for obstacle exerted on air

C d :

Drag coefficient for unit plant area density

c i :

Experimental constants of turbulent closure scheme (Mellor and Yamada, 1982)

E :

Turbulent kinetic energy (TKE)

F i :

Drag force in thei direction exerted on air flow by obstacle elements

f k :

Coriolis parameter

g i :

Acceleration vector due to gravity

H :

Height of shelterbelt

i, j, k, q :

Subscript variables, indicatingx, y andz directions, respectively, and grid numbers in these three directions

K 0 :

Turbulent exchange coefficient for neutral, obstacle-free ABL

K m :

Turbulent exchange coefficient for momentum transport

K E :

Turbulent exchange coefficient for TKE transport

k r :

Resistance coefficient of shelterbelts

l :

Mixing length of turbulence

MKE:

Mean kinetic energy

n :

Timestep of model integration

n,n i :

Vector and its component in thei direction of the interface of the averaging volume

p :

Atmospheric pressure perturbation

S :

Interface surface of the averaging volume

t :

Time

U :

Total mean windspeed

u, w :

Mean windspeed components inx andz directions, respectively

u′, w′ :

Fluctuating windspeed components inx andz directions

u i :

Windspeed in thei direction

u * :

Friction velocity

u *0 :

Friction velocity for obstacle-free ABL

u t ,v t :

u andv at model top

u aux, andw aux :

Intermediate prediction velocities ofu andw without the dynamic pressure perturbation

u n+1 andw n+1 :

Prediction velocities ofu andw at then+1 timestep of model integration

V :

Volume of the spatial averaging process

x :

Horizontal coordinate axis perpendicular to shelterbelt

x i :

i=1, 2, 3-three direction coordinate,x, y, z

z :

Vertical coordinate axis upward

z 0 :

Ground surface roughness length

β:

Weight coefficient for numerical differencing scheme

γ:

Coefficient of air thermal expansion

∈:

Dissipation rate of turbulence

ijk :

Einstein summation symbol

ρ0:

Air density

ϑ:

Potential-temperature departure from its basic state

ν:

Coefficient of air molecular viscosity

t :

Time step of model integration

κ:

von Karman constant

ϕ:

Macry symbol, standing foru, v, w, E andE1

References

  • Baines, W. D. and Peterson, E. G.: 1951, ‘An Investigation of Flow Through Screens’,Trans. Am. Soc. Mech. Eng. 73, 467–480.

    Google Scholar 

  • Bradley, E. F. and Mulhearn, P. J.: 1983, ‘Development of Velocity and Shear Stress Distributions in the Wake of a Porous Shelter Fence’,J. Wind Eng. Ind. Aerodyn. 15, 145–156.

    Google Scholar 

  • Caborn, J. M.: 1957, ‘Shelterbelts and Microclimate’,For. Comm. Bull. No. 29, Edinburgh, 129 pp.

  • Castro, I. P.: 1971, ‘Wake Characteristics of Two-Dimensional Perforated Plates Normal to an Airstream’,J. Fluid Mech. 46, 599–609.

    Google Scholar 

  • Chang, S. C.: 1966, ‘Velocity Distributions in the Separated Flow behind a Wedge-shaped Model Hill’, Thesis, Colorado State University, Ft. Collins, CO, 101 pp.

    Google Scholar 

  • Chorin, A. J.: 1968, ‘Numerical Solution of the Navier-Stokes Equations’,Math. Comp. 23, 341–354.

    Google Scholar 

  • Counihan, J., Hunt, J. C. R., and Jackson, P. S.: 1974, ‘Wakes Behind Two-Dimensional Surface Obstacles in Turbulent Boundary Layers’,J. Fluid Mech. 64, 529–563.

    Google Scholar 

  • Finnigan, J. J.: 1985, ‘Turbulent Transport in Flexible Plant Canopies’, in B. A. Hutchison and B. B. Hicks (eds.),The Forest-Atmosphere Interaction, D. Reidel Publishing Company, Boston, pp. 443–480.

    Google Scholar 

  • Finnigan, J. J. and Bradley, E. F.: 1983, ‘The Turbulent Kinetic Energy Budget Behind a Porous Barrier: an Analysis in Streamline Coordinates’,J. Wind. Eng. Ind. Aerodyn. 15, 145–156.

    Google Scholar 

  • Hagen, L. J. and Skidmore, E. L.: 1971, ‘Turbulent Velocity Fluctuations and Vertical Flow as Affected by Windbreak Porosity’Trans. ASAE 14, 634–637.

    Google Scholar 

  • Hagen, L. J., Skidmore, E. L., Miller, P. L., and Kipp, J. E.: 1981, ‘Simulation of Effect of Wind Barriers on Airflow’,Trans. ASAE 24, 1002–1008.

    Google Scholar 

  • Heisler, G. M. and Dewalle, D. R.: 1988, ‘Effects of Windbreak Structure on Wind Flow’,Agriculture, Ecosystems and Environment 22/23, 41–69.

    Google Scholar 

  • Hoerner, S. F.: 1965, ‘Fluid Dynamic Drag’, Library of Congress Catalog Card Number 64-19666.

  • Jensen, M.: 1974, ‘The Aerodynamics of Shelter’. In:FAO Report on the FAO/DANIDA Inter-regional Training Center on Heathland and Sand Dune Afforestation. FAO/DEN/TF, 123 pp.

  • Kaiser, H.: 1959, ‘Die Stromung an Windschutzstreifen’,Ber. Deut. Wetterdiensters 7, 36 pp.

    Google Scholar 

  • Li, Z., Lin, J. D., and Miller, D. R.: 1989, ‘Air Flow Over and Through a Forest Edge: A Steady-State Numerical Simulation’,Boundary-Layer Meteorol. 46, 333–354.

    Google Scholar 

  • McNaughton, K. G.: 1988, ‘Effects of Windbreaks on Turbulent Transport and Microclimate’,Agriculture, Ecosystems and Environment 22/23, 17–39.

    Google Scholar 

  • Mellor, G. L. and Yamada, T.: 1974, ‘A Hierarchy of Turbulence Closure Models for Planetary Boundary Layers’,J. Atmos. Sci. 31, 1791–1806.

    Google Scholar 

  • Mellor, G. L. and Yamada, T.: 1982, ‘Development of a Turbulent Closure Model for Geophysical Fluid Problems’,Rev. Geophys. Space Sci. 20, 851–875.

    Google Scholar 

  • Meyers, T. and Paw U, K. T.: 1986, ‘Testing of A Higher-Order Closure Model for Modeling Airflow Within and Above Plant Canopies’,Boundary-Layer Meteorol. 37, 297–311.

    Google Scholar 

  • Miller, D. R., Lin, J. D., and Lu, Z. N.: 1991, ‘Air Flow Across an Alpine Forest Clearing: A Model and Field Measurements’,Agricultural and Forest Meteorol. 56, 209–225.

    Google Scholar 

  • Naot, O. and Mahrer, Y.: 1991, ‘Two-Dimensional Microclimate Distribution Within and Above a Crop Canopy in an Arid Environment: Modeling and Observational Studies’,Boundary-Layer Meteorol. 56, 223–244.

    Google Scholar 

  • Paegle, J., Zdunkowski, W. G., and Welch, R. M.: 1976, ‘Implicit Differencing of Predictive Equations of the Boundary Layer’,Mon. Wea. Rev. 104, 1321–1324.

    Google Scholar 

  • Perera, M. D. A. S.: 1981, ‘Shelter Behind Two-Dimensional Solid and Porous Fences’,J. Wind Eng. Ind. Aerodyn. 8, 93–104.

    Google Scholar 

  • Plate, E. J.: 1971, ‘The Aerodynamics of Shelterbelts’,Agricultural Meteorol. 8, 203–222.

    Google Scholar 

  • Raine, J. K. and Stevenson, D. C.: 1977, ‘Wind Protection by Model Fences in a Simulated Atmospheric Boundary Layer’,J. Indust. Aerodyn. 2, 159–180.

    Google Scholar 

  • Raupach, M. R. and Shaw, R. H.: 1982, ‘Averaging Procedures for Flow Within Vegetation Canopies’,Boundary-Layer Meteorol. 22, 79–90.

    Google Scholar 

  • Richards, P. J., Kay, E. F., Russell, D., and Wilson, G. R. C.: 1984, ‘Porous Artificial Windbreaks in Oblique Winds’,Paper 67/84 for IPENZ Conf., Hastings, New Zealand, 10 pp.

  • Rosenberg, N. J.: 1979, ‘Windbreaks for Reducing Moisture Stress’, in B. J. Barfield and J. F. Gerber (eds.),Modification of Aerial Environment of Plants, ASAE, St Joseph, pp. 538.

    Google Scholar 

  • Sturrock, J. W.: 1969, ‘Aerodynamics Studies of Shelterbelts in Nealand-1: Low to Medium Height Shelterbelts in Mid-Canterbury’,N. Z. J. Sci. 12,754–776.

    Google Scholar 

  • Sturrock, J. W.: 1972, ‘Aerodynamics studies of Shelterbelts in Nealand-2: Medium Height to Tall Shelterbelts in Mid-Canterbury’,N. Z. J. Sci. 15, 113–140.

    Google Scholar 

  • Taylor, P. A.: 1988, ‘Turbulent Wakes in the Atmospheric Boundary Layer’, in W. L. Steffen and O. T. Denmead (eds.),Flow and Transport in the Natural Environment: Advances and Applications, Springer-Verlag, Berlin, pp. 270–292.

    Google Scholar 

  • Thom, A.S.: 1975, ‘Momentum, Mass and Heat Exchange of Plant Communities’,Vegetation and the Atmosphere, Vol. 1, Academic Press, pp. 1–278.

    Google Scholar 

  • van Eimern, J., Karschon, R., Razumova, L. A., and Robertson, G. W.: 1964, ‘Windbreaks and Shelterbelts’, World Meteorol. Org. Technical Note No, 59, pp. 188.

    Google Scholar 

  • Wang, H., and Takle, E. S.: 1995, ‘Boundary-Layer Flow and Turbulence Near Porous Obstacles. I. Derivation of a General Equation Set for a Porous Medium’,Boundary-Layer Meteorol. (in press).

  • Whitaker, S.: 1968,Introduction to Fluid Mechanics, Prentice-Hall, Englewood Cliffs.

    Google Scholar 

  • Whitaker, S.: 1969, ‘Advances in Theory of Fluid Motions in Porous Media’,Ind. Eng. Chem. 61, 14–28.

    Google Scholar 

  • Whitaker, S.: 1973, ‘The Transport Equations for a Multi-phase System’,Chem. Eng. Sci. 28, 139–147.

    Google Scholar 

  • Wilson, N. R. and Shaw, R. H.: 1977, ‘A Higher Order Closure Model for Canopy Flow’,J. Applied Meteorol. 16, 1197–1205.

    Google Scholar 

  • Wilson, J. D.: 1985, ‘Numerical Studies of Flow Through a Windbreak’,J. Wind Eng. Ind. Aerodyn. 21, 119–154.

    Google Scholar 

  • Wilson, J. D.: 1987, ‘On the Choice of a Windbreak Porosity Profile’,Boundary-Layer Meteorol. 38, 37–39.

    Google Scholar 

  • Wilson, J. D.: 1988, ‘A Second-Order Closure Model for Flow Through Vegetation’,Boundary-Layer Meteorol. 42, 371–392.

    Google Scholar 

  • Yamada, T.: 1982, ‘A Numerical Model Study of Turbulent Airflow in and above a Forest Canopy’,J. Meteorol. Soc. Japan 60, 438–454.

    Google Scholar 

  • Yamada, T. and Mellor, G. L.: 1975, ‘A Simulation of the Wangara Atmospheric Boundary Layer Data’,J. Atmos. Sci. 32, 2309–2329.

    Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Rights and permissions

Reprints and permissions

About this article

Cite this article

Wang, H., Takle, E.S. A numerical simulation of boundary-layer flows near shelterbelts. Boundary-Layer Meteorol 75, 141–173 (1995). https://doi.org/10.1007/BF00721047

Download citation

  • Revised:

  • Issue Date:

  • DOI: https://doi.org/10.1007/BF00721047

Keywords

Navigation