The Urban Canopy and the Plant Canopy
|
|
- Felicity Fisher
- 6 years ago
- Views:
Transcription
1 The Urban Canopy and the Plant Canopy John Finnigan: CSIRO Australia Margi Bohm: University of Canberra Roger Shaw: U C Davis Ned Patton: NCAR Ian Harman: CSIRO Australia
2 Properties of turbulent shear flows over deep rough surfaces 1. Turbulence structure in plant canopies and the Roughness Sub Layer Current understanding of ISL eddy structure New evidence of RSL eddy structure How does it differ from the Inertial Sub Layer (Logarithmic layer) 3. A phenomenological model for Canopy-RSL turbulence 3. Turbulence structure in urban canopies and the Roughness Sub Layer How similar is it to that in a vegetation canopy? 4. Extending the model for large eddies in vegetation canopies to urban canopies When should it work and when should it fail? A new hypothesis 5. Conclusions
3 Urban canopies and plant canopies absorb momentum over a finite depth, not at a surface plane
4 Boundary layer Velocity Profiles Smooth wall Rough wall U U
5 Monin-Obukhov Similarity Theory (MOST) fails close to rough surfaces like plant or urban canopies Wind Potential temperature observations MOST Flow over tall (plant) canopies does not conform to MOST to within 2-3 canopy heights this is the roughness sublayer. NB these are not errors in choice of d Chen and Schwerdtfeger (1989)
6 Flow statistics in the RSL are distinctly different from the ISL above: eddy coherence EOF spectra converge more rapidly than in the ISL, suggesting turbulence is dominated by coherent eddies with a distinct spatial scale (typically >80% of tke in first 5 eigenmodes) EOF analysis of wind tunnel canopy model: Finnigan and Shaw (2000)
7 Flow statistics in the RSL are distinctly different from the ISL above: correlations, higher moments, momentum and scalar exchange The correlation coefficient is ~0.5 in the neutrally stratified RSL/canopy as compared to ~0.35 in the ISL. Together with the smaller gradients for the same flux, this shows that RSL/canopy turbulence is in some sense more efficient at transport than that in the ISL. Skewnesses are O[1] and transport terms in covariance budgets are large-the turbulence is not in local equilibrium The turbulent Schmidt and Prandtl numbers are ~0.5 at the canopy- RSL interface compared with ~1.0 in the ISL. (Raupach, Finnigan and Brunet (1996))
8 Flow statistics in the RSL are distinctly different from the ISL above: ejection: sweep ratios Transport of momentum and scalars in the RSL/ canopy layer is dominated by the sweep quadrant whereas through the rest of the ISL, ejections dominate. Circles: almond orchard diamonds: cork oak plantation triangles: wind tunnel model canopy squares: DNS of urban blocks Ratios Q2/Q4 of the contributions to momentum flux from ejections and sweeps from a range of canopies
9 What is the origin of these differences? It is reasonable to assume that turbulence in the RSL must have a structure that is distinctly different from that in the ISL We will look first at the current thinking on eddy structure in the ISL over smooth walls This is obtained mainly from analysis of DNS data in boundary layers and channel flows. We will compare this with earlier and very recent analysis of eddy structure in vegetation canopies and RSLs
10 With a series of collaborators, Adrian has shown that the dominant, ejection-producing structures in the ISL are packets of head-up hairpin vortices The initial hairpin is presumed to formed by the upward deflection of a spanwise vortex line of the mean shear by a vertical gust. The initial HU hairpin generates an ejection that deflects another vortex line, eventually producing a packet aligned around an elongated low speed region (Gerz et al, 1994; Zhou et al, 1999, Tomkins and Adrian, 2003 )
11 What is the symmetry breaking mechanism that allows HU hairpins and hence ejections to dominate? In a homogeneous shear flow, head-up and head-down deflections are equally probable. Near a wall in contrast, inward gusts are blocked so that outward deflections dominate. The larger the initial deflection, the more the vortex loop is stretched and rotated by the mean shear and its vorticity amplified. U(z) Questions remain about the mechanism of initial scale selection. Over a smooth wall, a single dominant scale does not seem to be strongly selected at any level. Schematics from Gerz et al. 1994
12 Structure of Canopy Turbulence Time-height traces from single towers in tall canopies give information about the x-z plane
13 Scalar ramps correlated through the depth of the canopy show wholesale flushing of the canopy airspace by large scale gusts Gao et al (1989), Camp Borden, Canada
14 Compositing shows that these ramps are signals of a scalar microfront compressed between downwind ejections and upwind sweeps Gao et al (1989), Camp Borden, Canada
15 Structure of Canopy Turbulence 3d structure from detailed wind tunnel simulations 1. Tombstones Raupach, Coppin and Le
16 hm, Finnigan and Raupach (2000) 2. The Black Forest canopy
17 unet, Finnigan and Raupach (1994) 3. The waving wheet canopy
18 4. Water flume experiment over multiple hills covered with a deep canopy Poggi, D., Katul, G.G. and Finnigan, J.J., 2006
19 Finnigan and Hughes (2008) 5. Stably Stratified canopy model on a 2D Hill
20 3D Structure. The eddy structure of the RSL has been reconstructed from multi-point analysis of a WT model canopy and LES simulation of the same canopy using conditional sampling (compositing) and EOF analysis These are iso-surfaces of λ2, a measure constructed from the invariants of the du i /dx j tensor that captures aspects of the local vorticity and global rotation associated with the vortices (Jeong and Hussain, 1995). Data from canopy LES using NCAR code of Patton and Sullivan. Sampling triggered by p +
21 The dominant eddy structure of the Canopy-RSL has been reconstructed from multi-point analysis of a WT model canopy (wheet) and LES simulation of the same canopy using conditional sampling (compositing) and EOF analysis These are iso-surfaces of λ2, a measure constructed from the invariants of the dui/dxj tensor that captures aspects of the local vorticity and global rotation associated with the vortices (Jeong and Hussain, 1995). Coincidence of λ2 and vortex core is seen in this plot of velocity vectors projected on the y-z plane
22 y-z slice through the composite eddy
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40 u -w vectors from the same composite eddy on the x-z plane of symmetry
41 Convergence between the underlying ejection and overlying sweep produces a scalar microfront. Shear stress <u w > is concentrated between the hairpin legs The scalar is released from the canopy at a uniform rate (independent of local wind velocity). Within a structure, a sloping microfront is formed with high concentration below and in advance of the front, while low concentration follows and is above the microfront. Blue- λ2 isosurface Green-scalar microfront Red- u w sweep Orange-u w ejection
42 We have evidence that the Head-up and Head-Down hairpins are formed simultaneously as the linear instability theory suggests
43 Compositing using the high pressure trigger above the RSL produces a less coherent structure with the HU hairpin and ejection much stronger than the HD and sweep This composite is formed at z=3hc, about twice the height above the canopy where ejections first equal sweep contributions to Blue- λ2 isosurface Green-scalar microfront Red- u w sweep Orange-u w ejection
44 So what happens at the top of a canopy that is different to the ISL above in order to produce this distinct eddy structure? Unlike the boundary layer profile, the inflected velocity profile at canopy top is inviscidly unstable, leading to rapid growth and strong selection for a single scale, proportional to the vorticity thickness δ ω. Spanwise Stuart vortices develop which can be deflected into HU and HD hairpins. This is the mixing layer analogy (Raupach et al, 1996).
45 What happens at the top of a canopy that is different to the ISL above? Trains of transverse Stuart vortices are unstable and subject to a helical pairing instability which was described by Pierrehumbert and Widnall (1982). Non-linear development of the linear eigenmodes results in trains of Head-up and Head-down vortices spaced at twice the wavelength of the original Kelvin-Helmholtz inflection point instability
46 We have modelled a sequence of Stuart vortices in an inflected (tanh-1) velocity profile with an initial (isotropic) background turbulence intensity of 0.05% U
47 The top view shows the spanwise scale selection. Note, only the central part of the calculation domain is shown
48 Two opposing symmetry breaking mechanisms explain the shift from sweep dominance to ejection dominance as we move above the canopy The presence of the porous canopy allows HD hairpins to be deflected downwards -as long as their spanwise scale is <h c In the canopy-top shear flow, HD hairpins are stretched and rotated faster than HU hairpins so HD s dominate Further from the canopy top, large scale upward deflections become dominant again as downward deflections are blocked by the solid wall so that HU hairpins begin to dominate.
49 Vegetation canopy/rsl Eddy structure: Summary 1 The inflected velocity profile at the canopy top provides strong scale selection by an inviscid instability mechanism (unlike in a BL profile) The secondary instability produces trains of coherent spanwise vortices that are unstable to small perturbations The tertiary (helical pairing) instability preferentially produces HUs and HDs in pairs The convergence between the sweeps and ejections produces intense scalar microfronts
50 Vegetation canopy/rsl Eddy structure: Summary 2 Two symmetry-breaking mechanisms act in opposition to determine whether HUs or HDs dominate Near solid walls, large downward deflections are blocked so energetic HD hairpins cannot be formed. A porous canopy layer allows HD deflections of order h c near the canopy top. In a boundary-layer shear flow HD hairpins are stretched and rotated faster than HUs and so sweeps dominate ejections. Further from the canopy top the larger amplitude deflection of HUs overrides the greater local stretching of HDs and so ejections dominate
51 RSL model obtained by using length in MOST theory as an extra scaling LAI ~ 2 β = 0.31 n=17 Tumbarumba LAI ~ 3 β = 0.39 n=42 Duke LAI ~3.8 β = 0.28 n=27 Harman and Finnigan (2007, 2008) increasingly dense
52 Urban canopy statistics Q2/Q4 Castro, Cheng and Reynolds (2006) Various roughness types Bohm, Finnigan and Raupach (2009) Black Forest vs LES of Wheet WT canopy Finnigan, Shaw and Patton (2009), Coceal et al. (2006) Blocks vs vegetation
53 But, not all urban canopies have Q4>Q2 Diagonal Array of blocks Square Array of blocks Kanda, 2006 LES data
54 And not all urban canopies have r uw >0.35 Diagonal or Staggered Array of blocks Random array Square Array of blocks Cheng and Castro (2001) WT data
55 And not all urban canopies have r uw >0.35 The Black Forest WT model canopy of smooth bluff objects shows characteristics of canopy (K type roughness) and boundary layer (D type roughness) statistics Bohm, Finnigan and Raupach (2009)
56 A Hypothesis: Urban canopies that behave like plant canopies have the same statistical eigenmodes as attractors for their large eddies In vegetation canopies, the inflection point in the mean velocity is inviscidly unstable and generates energetic coherent eddies The wakes shed by plant elements are much smaller in scale than the eigenmodes generated by the inflected mean velocity profile Because of this, the inflected velocity profile is dynamically significant even though instantaneously it is distorted by plant wakes Under what circumstances could these dynamics be relevant to urban canopies, where the size of the element wakes is of the same order as the large eddies? We will examine the equation governing the relationship of a local velocity perturbation to the spatially averaged, inflected mean velocity profile
57 A note on the use of linear theory in fully turbulent flows Perhaps surprisingly, there is a long history of successful application of linear theory to explain the eddy structure of fully developed turbulent shear flows It has been applied to understand large-scale structures in free turbulent shear flows like mixing layers where the coherent eddies correspond to unstable eigenmodes (eg. Liu, 1988) And to boundary layer and channel flows where the coherent eddies correspond to singular neutral modes. We call this Rapid distortion Theory (Townsend, 1976; Hunt and Carruthers, 1990). In each case the central assumption is that the dominant eddy structures are determined by the interaction of the turbulence with the mean flow and that non-linear, turbulent-turbulent interactions act primarily to destroy the large eddies.
58 Averaging the velocity and Reynolds stresses in time and space x x
59 Equation for the local (in space and time) velocity perturbation Terms in black are considered in linearized approaches to analysing canopy large eddy structure such as hydrodynamic stability and rapid distortion theory. Terms in red represent the distortion of the local velocity perturbation by local variations in the time-averaged strain rates around buildings or blocks The non linear term in blue contains both the spatially averaged Reynolds stress and the local variations around this.
60 Mean Velocity and its spatial gradient in a canopy of blocks Coceal, Thomas and Belcher (2007)
61 Reynolds Shear Stress and its spatial gradient in a canopy of blocks Coceal, Thomas and Belcher (2007)
62 Hypothesis If the mean velocity gradients and Reynolds stress gradients vary on a scale that is similar or larger than the spatial scale of the inflexionpoint eigenmode, then we will not observe coherent eddies similar to those in vegetation canopies. Instead the turbulence will be dominated entirely by the eddies shed from individual obstacles and will have no universal features If the opposite is true, double hairpin type eddies will emerge and the characteristic features of vegetation canopy turbulence will be seen. The eigenmode size is linked to the vorticity thickness Urban canopies with large relative to element size as a result of varying building height or complex element geometry, may, therefore, have dynamically significant inflection point profiles Note that this hypothesis breaks down if the canopy elements are too sparse
63 Example from Black Forest
64 Conclusions The turbulence in and just above vegetation canopies is significantly more coherent than in the ISL or Log layer above and more efficient at transporting momentum and scalars. We can extract the space-time structure of canopy RSL eddies by compositing and we find that they have a double hairpin structure that is quite different to the Head-up hairpins or attached eddies generally assumed in smooth wall shear layers. A sequence of instability processes triggered by the inflected shear layer at the canopy top together with two opposing symmetry-breaking mechanisms successfully explains the eddy structure as well as a very wide range of observations in vegetation canopies
65 Conclusions In urban canopies, the crucial inflexion point is also present in the spatially averaged time-mean velocity profile but may not be dynamically significant Furthermore, not all urban canopy data show the universal features seen in vegetation canopies such as Q2<Q4 or r uw ~0.5 However, by comparing the magnitudes of terms in the equation for a local velocity perturbation, we have suggested a simple criterion to distinguish urban canopies that behave in a similar way to natural canopies from those that do not This criterion remains to be tested against the data.
Roughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman
Roughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman 1. Characteristics of the Roughness Sub layer With well understood caveats, the time averaged statistics of flow in the atmospheric
More informationJ2.1 EDDY STRUCTURE NEAR THE PLANT CANOPY INTERFACE
J. EDDY STRUTURE NEAR THE PLANT ANOPY INTERFAE Roger H. Shaw * University of alifornia, Davis John J. Finnigan SIRO, Marine and Atmospheric Research, anberra, Australia Edward G. Patton National enter
More informationThe plant canopy layer and roughness sublayer. 16 Oct., 2012
The plant canopy layer and roughness sublayer 16 Oct., 2012 the surface layer over a rough boundary must be considered in two parts: an inertial sublayer in which height above the effective surface provides
More informationTURBULENT STATISTICS OF NEUTRALLY STRATIFIED FLOW WITHIN AND ABOVE A SPARSE FOREST FROM LARGE-EDDY SIMULATION AND FIELD OBSERVATIONS
TURBULENT STATISTICS OF NEUTRALLY STRATIFIED FLOW WITHIN AND ABOVE A SPARSE FOREST FROM LARGE-EDDY SIMULATION AND FIELD OBSERVATIONS HONG-BING SU 1,, ROGER H. SHAW 1,KYAWTHAPAWU 1, CHIN-HOH MOENG 2 and
More informationA NOTE ON THE CONTRIBUTION OF DISPERSIVE FLUXES TO MOMENTUM TRANSFER WITHIN CANOPIES. Research Note
A NOTE ON THE CONTRIBUTION OF DISPERSIVE FLUXES TO MOMENTUM TRANSFER WITHIN CANOPIES Research Note D. POGGI Dipartimento di Idraulica, Trasporti ed Infrastrutture Civili, Politecnico di Torino, Torino,
More informationCanopy flow. Acknowledged contribution from Corey Markfort, Iowa State U
Canopy flow Examples of canopy flows canopy characterization Failure of log-law and K-theory in and over canopies Roughness sublayer and canopy flow dynamics Mixing layer analogy for canopies Potential
More informationThe structure of canopy turbulence and its implication to scalar dispersion
The structure of canopy turbulence and its implication to scalar dispersion Gabriel Katul 1,2,3 & Davide Poggi 3 1 Nicholas School of the Environment and Earth Sciences, Duke University, USA 2 Department
More informationTurbulent boundary layer
Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic
More informationLES, lab and empirical representations of a neutral flow over canopy
LES, lab and empirical representations of a neutral flow over canopy 23-06-2016, Huug Ouwersloot, Arnold Moene, Jisk Attema & Jordi Vilà-Guerau de Arellano Based on Ouwersloot et al. (Boundary-Layer Meteorology,
More informationTemperature fronts and vortical structures in turbulent stably stratified atmospheric boundary layers
VIIIth International Symposium on Stratified Flows August 29 - September 1 2016, San Diego, CA Temperature fronts and vortical structures in turbulent stably stratified atmospheric boundary layers Peter
More informationESPM 129 Biometeorology Wind and Turbulence, Part 2, Canopy Air Space: Observations and Principles
Lecture 21 Wind and Turbulence, Part 2, Canopy Air Space: Observations and Principles Instructor: Dennis Baldocchi Professor of Biometeorology Ecosystem Science Division Department of Environmental Science,
More informationStructure of turbulent flow over regular arrays of cubical roughness
Structure of turbulent flow over regular arrays of cubical roughness Journal: Manuscript ID: mss type: Date Submitted by the Author: Complete List of Authors: Keyword: Journal of Fluid Mechanics JFM-6-S-35.R3
More informationHORSEHOE VORTICES IN UNIFORMLY SHEARED TURBULENCE
HORSEHOE VORTICES IN UNIFORMLY SHEARED TURBULENCE Christina Vanderwel Department of Mechanical Engineering University of Ottawa Ottawa, Ontario, Canada cvand072@uottawa.ca Stavros Tavoularis Department
More informationarxiv: v1 [physics.flu-dyn] 7 Dec 2018
This draft was prepared using the LaTeX style file belonging to the Journal of Fluid Mechanics arxiv:8.8v [physics.flu-dyn] 7 Dec 8 Contrasts Between Momentum and Scalar Transport Over Very Rough Surfaces
More informationTWO-POINT CORRELATION ANALYSIS OF NEUTRALLY STRATIFIED FLOW WITHIN AND ABOVE A FOREST FROM LARGE-EDDY SIMULATION
TWO-POINT CORRELATION ANALYSIS OF NEUTRALLY STRATIFIED FLOW WITHIN AND ABOVE A FOREST FROM LARGE-EDDY SIMULATION HONG-BING SU 1,2, ROGER H. SHAW 1 and KYAW THA PAW U 1 1 Atmospheric Science Program, University
More informationKeywords: Large-eddy simulation, Turbulent coherent structure, Four quadrant analysis, Integral scale
The Eighth Asia-Pacific Conference on Wind Engineering, December 4, 3, Chennai, India NUMERICAL ANALYSIS OF THE MOMENTUM TRANSPORT AND TEMPORAL AND SPATIAL SCALES OF TURBULENT COHERENT STRUCTURES IN THE
More informationLarge-Eddy Simulation for Turbulent Nature of Flow and Pressure Fields over Urban Building Arrays C. Hirose*, A. Hagishima, N. Ikegaya, and J. Tanimot
Large-Eddy Simulation for Turbulent Nature of Flow and Pressure Fields over Urban Building Arrays C. Hirose*, A. Hagishima, N. Ikegaya, and J. Tanimoto Interdisciplinary Graduate School of Engineering
More informationRegularity diagnostics applied to a turbulent boundary layer
Center for Turbulence Research Proceedings of the Summer Program 208 247 Regularity diagnostics applied to a turbulent boundary layer By H. J. Bae, J. D. Gibbon, R. M. Kerr AND A. Lozano-Durán Regularity
More informationFirst order turbulence closure for modelling complex canopy flows
First order turbulence closure for modelling complex canopy flows Article Accepted Version Finnigan, J., Harman, I., Ross, A. and Belcher, S. (015) Firstorder turbulence closure for modelling complex canopy
More informationLarge Eddy Simulation of an Inhomogeneous Atmospheric Boundary Layer under Neutral Conditions
2479 Large Eddy Simulation of an Inhomogeneous Atmospheric Boundary Layer under Neutral Conditions CHING-LONG LIN Department of Mechanical and Industrial Engineering and IIHR Hydroscience and Engineering,
More informationDual Vortex Structure Shedding from Low Aspect Ratio, Surface-mounted Pyramids
Dual Vortex Structure Shedding from Low Aspect Ratio, Surface-mounted Pyramids Robert J. Martinuzzi Department of Mechanical and Manufacturing Engineering Schulich School of Engineering University of Calgary
More informationChapter 5 Phenomena of laminar-turbulent boundary layer transition (including free shear layers)
Chapter 5 Phenomena of laminar-turbulent boundary layer transition (including free shear layers) T-S Leu May. 3, 2018 Chapter 5: Phenomena of laminar-turbulent boundary layer transition (including free
More informationLogarithmic velocity profile in the atmospheric (rough wall) boundary layer
Logarithmic velocity profile in the atmospheric (rough wall) boundary layer P =< u w > U z = u 2 U z ~ ε = u 3 /kz Mean velocity profile in the Atmospheric Boundary layer Experimentally it was found that
More informationAtm S 547 Boundary Layer Meteorology
Lecture 5. The logarithmic sublayer and surface roughness In this lecture Similarity theory for the logarithmic sublayer. Characterization of different land and water surfaces for surface flux parameterization
More informationLecture 2. Turbulent Flow
Lecture 2. Turbulent Flow Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of this turbulent water jet. If L is the size of the largest eddies, only very small
More informationTurbulent eddies in the RANS/LES transition region
Turbulent eddies in the RANS/LES transition region Ugo Piomelli Senthil Radhakrishnan Giuseppe De Prisco University of Maryland College Park, MD, USA Research sponsored by the ONR and AFOSR Outline Motivation
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More informationThe mean velocity profile in the smooth wall turbulent boundary layer : 1) viscous sublayer
The mean velocity profile in the smooth wall turbulent boundary layer : 1) viscous sublayer u = τ 0 y μ τ = μ du dy the velocity varies linearly, as a Couette flow (moving upper wall). Thus, the shear
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationLOW SPEED STREAKS INSTABILITY OF TURBULENT BOUNDARY LAYER FLOWS WITH ADVERSE PRESSURE GRADIENT
LOW SPEED STREAKS INSTABILITY OF TURBULENT BOUNDARY LAYER FLOWS WITH ADVERSE PRESSURE GRADIENT J.-P. Laval (1) CNRS, UMR 8107, F-59650 Villeneuve d Ascq, France (2) Univ Lille Nord de France, F-59000 Lille,
More informationTransition to turbulence in plane Poiseuille flow
Proceedings of the 55th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, February 25-26, 2015 ThL2T5.1 Transition to turbulence in plane Poiseuille flow F. Roizner, M. Karp and
More informationThis is the first of several lectures on flux measurements. We will start with the simplest and earliest method, flux gradient or K theory techniques
This is the first of several lectures on flux measurements. We will start with the simplest and earliest method, flux gradient or K theory techniques 1 Fluxes, or technically flux densities, are the number
More informationHomogeneous Turbulence Dynamics
Homogeneous Turbulence Dynamics PIERRE SAGAUT Universite Pierre et Marie Curie CLAUDE CAMBON Ecole Centrale de Lyon «Hf CAMBRIDGE Щ0 UNIVERSITY PRESS Abbreviations Used in This Book page xvi 1 Introduction
More informationThe relative importance of ejections and sweeps to momentum transfer in the atmospheric boundary layer
DOI.7/s546-6-964-6 ORIGINAL ARTICLE The relative importance of ejections and sweeps to momentum transfer in the atmospheric boundary layer Gabriel Katul Davide Poggi Daniela Cava John Finnigan Received:
More informationNowadays, the rapid development of computer resources has enabled the numerical simulation based on the computational fluid dynamics (CFD) techniques
Large-Eddy Simulation On The Gust Probability In Urban Pedestrian Spaces Y. Ikeda 1,*, A. Hagishima 1, N. Ikegaya 1, and J. Tanimoto 1 1 Interdisciplinary Graduate School of Engineering Science, Kyushu
More informationTurbulent flow over anisotropic porous media
Turbulent flow over anisotropic porous media Alfredo Pinelli School of Mathematics, Computer Science and Engineering City, University of London U.K. Work done in collaboration with: M. Omidyeganeh, A.
More informationThe role of coherent structures in subfilter-scale dissipation of turbulence measured in the atmospheric surface layer
JOT J OURNAL OF TURBULENCE http://jot.iop.org/ The role of coherent structures in subfilter-scale dissipation of turbulence measured in the atmospheric surface layer Matthew A Carper 1 and Fernando Porté-Agel
More informationSimulation of Flow around a Surface-mounted Square-section Cylinder of Aspect Ratio Four
Simulation of Flow around a Surface-mounted Square-section Cylinder of Aspect Ratio Four You Qin Wang 1, Peter L. Jackson 2 and Jueyi Sui 2 1 High Performance Computing Laboratory, College of Science and
More informationFirst-order turbulence closure for modelling complex canopy flows
QuarterlyJournalof theroyalmeteorologicalsociety Q. J. R. Meteorol. Soc. 141: 2907 2916, October 2015 A DOI:10.1002/qj.2577 First-order turbulence closure for modelling complex canopy flows John Finnigan,
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More information6 VORTICITY DYNAMICS 41
6 VORTICITY DYNAMICS 41 6 VORTICITY DYNAMICS As mentioned in the introduction, turbulence is rotational and characterized by large uctuations in vorticity. In this section we would like to identify some
More informationτ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization
τ xz = τ measured close to the the surface (often at z=5m) these three scales represent inner unit or near wall normalization Note that w *3 /z i is used to normalized the TKE equation in case of free
More informationTurbulent Pressure and Velocity Perturbations Induced by Gentle Hills Covered with Sparse and Dense Canopies
Boundary-Layer Meteorol DOI.7/s546-9-947-x ARTICLE Turbulent Pressure and Velocity Perturbations Induced by Gentle Hills Covered with Sparse and Dense Canopies Edward G. Patton Gabriel G. Katul Received:
More informationWall turbulence with arbitrary mean velocity profiles
Center for Turbulence Research Annual Research Briefs 7 Wall turbulence with arbitrary mean velocity profiles By J. Jiménez. Motivation The original motivation for this work was an attempt to shorten the
More informationNUMERICAL INVESTIGATION OF THE FLOW OVER A GOLF BALL IN THE SUBCRITICAL AND SUPERCRITICAL REGIMES
NUMERICAL INVESTIGATION OF THE FLOW OVER A GOLF BALL IN THE SUBCRITICAL AND SUPERCRITICAL REGIMES Clinton Smith 1, Nikolaos Beratlis 2, Elias Balaras 2, Kyle Squires 1, and Masaya Tsunoda 3 ABSTRACT Direct
More informationConvective Fluxes: Sensible and Latent Heat Convective Fluxes Convective fluxes require Vertical gradient of temperature / water AND Turbulence ( mixing ) Vertical gradient, but no turbulence: only very
More informationSurface Roughness Effects in Near-Bed Turbulence: Implications to Sediment Entrainment
Surface Roughness Effects in Near-Bed Turbulence: Implications to Sediment Entrainment by A.N. Papanicolaou, P. Diplas, C. L. Dancey, and M. Balakrishnan Journal of Engineering Mechanics, Vol. 127, No.
More informationChallenges of modelling wind engineering problems
Challenges of modelling wind engineering problems Zheng-Tong Xie With thanks to: Vladimir Fuka, Paul Hayden, Ian Castro, Alan Robins, Janet Barlow, Yusik Kim, Bob Plant, Omduth Coceal, Denise Hertwig,
More informationPolymer maximum drag reduction: A unique transitional state
Center for Turbulence Research Annual Research Briefs 47 Polymer maximum drag reduction: A unique transitional state By Y. Dubief, C. M. White, E. S. G. Shaqfeh AND V. E. Terrapon. Motivation and objectives
More information2013 Annual Report for Project on Isopycnal Transport and Mixing of Tracers by Submesoscale Flows Formed at Wind-Driven Ocean Fronts
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. 2013 Annual Report for Project on Isopycnal Transport and Mixing of Tracers by Submesoscale Flows Formed at Wind-Driven
More informationProbability density functions of turbulent velocity and temperature in the atmospheric surface layer
WATER RESOURCES RESEARCH, VOL. 32, NO. 6, PAGES 68 688, JUNE 996 Probability density functions of turbulent velocity and temperature in the atmospheric surface layer Chia R. Chu,,2 Marc B. Parlange, Gabriel
More information$%&i&fr 4e. ?s~ /zy REPORT DOCUMENTATION PAGE. tw'*/<?/ //s. VS^/I'M /V'9%s>lJ /I F5^c#>c7^ Crf /f^z-e^st^sy
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 9^^^ «or rlw.no Instructions, searchmg.,,,, data ^ - collection of information, including suggestions for reducing this burden to Washmotor«2,?nl?
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
More informationBoundary Layer Meteorology The wind that shakes the buildings
Boundary Layer Meteorology The wind that shakes the buildings Prof. Janet F. Barlow Department of Meteorology Director Technologies for Sustainable Built Environments (TSBE) Centre University of Reading
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationLecture 3: The Navier-Stokes Equations: Topological aspects
Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without
More informationAnalysis of multi-plane PIV measurements in a turbulent boundary layer: large scale structures, coupled and decoupled motions
Perry Fest 2004 Queen s University, Kingston Analysis of multi-plane PIV measurements in a turbulent boundary layer: large scale structures, coupled and decoupled motions Ivan Marusic, Nick Hutchins, Will
More informationThe interaction of vorticity and rate-of-strain in homogeneous sheared turbulence
PHYSICS OF FLUIDS VOLUME 12, NUMBER 4 APRIL 2000 The interaction of vorticity and rate-of-strain in homogeneous sheared turbulence K. K. Nomura a) and P. J. Diamessis Department of Mechanical and Aerospace
More informationA NUMERICAL STUDY OF NOCTURNAL WAVELIKE MOTION IN FORESTS. 1. Introduction
A NUMERICAL STUDY OF NOCTURNAL WAVELIKE MOTION IN FORESTS XINZHANG HU 1, XUHUI LEE 1,, DAVID E. STEVENS 2 and RONALD B. SMITH 3 1 School of Forestry and Environmental Studies, Yale University, Connecticut,
More informationProceedings of the 4th Joint US-European Fluids Engineering Division Summer Meeting ASME-FEDSM2014 August 3-7, 2014, Chicago, Illinois, USA
Proceedings of the 4th Joint US-European Fluids Engineering Division Summer Meeting ASME-FEDSM4 August 3-7, 4, Chicago, Illinois, USA FEDSM4-38 SUPPRESSION OF UNSTEADY VORTEX SHEDDING FROM A CIRCULAR CYLINDER
More informationFor example, for values of A x = 0 m /s, f 0 s, and L = 0 km, then E h = 0. and the motion may be influenced by horizontal friction if Corioli
Lecture. Equations of Motion Scaling, Non-dimensional Numbers, Stability and Mixing We have learned how to express the forces per unit mass that cause acceleration in the ocean, except for the tidal forces
More information(Wind profile) Chapter five. 5.1 The Nature of Airflow over the surface:
Chapter five (Wind profile) 5.1 The Nature of Airflow over the surface: The fluid moving over a level surface exerts a horizontal force on the surface in the direction of motion of the fluid, such a drag
More informationA canopy model of mean winds through urban areas
Q. J. R. Meteorol. Soc. (24), 3, pp. 349 372 doi:.256/qj.3.4 A canopy model of mean winds through urban areas By O. COCEAL and S. E. BELCHER University of Reading, UK (Received 3 March 23; revised 24 September
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationA Note on Spatial Averaging and Shear Stresses Within Urban Canopies
Boundary-Layer Meteorol (2018) 167:171 179 https://doi.org/10.1007/s10546-017-0321-7 NOTES AND COMMENTS A Note on Spatial Averaging and Shear Stresses Within Urban Canopies Zheng-Tong Xie 1 Vladimir Fuka
More informationDirect Numerical Simulations of Transitional Flow in Turbomachinery
Direct Numerical Simulations of Transitional Flow in Turbomachinery J.G. Wissink and W. Rodi Institute for Hydromechanics University of Karlsruhe Unsteady transitional flow over turbine blades Periodic
More informationPECASE: MULTI-SCALE EXPERIMENTS AND MODELING
PECASE: MULTI-SCALE EXPERIMENTS AND MODELING IN WALL TURBULENCE AFOSR # FA9550-09-1-0701 (P.M. J. Schmisseur) Beverley J. McKeon Graduate Aerospace Laboratories California Institute of Technology http://mckeon.caltech.edu
More informationFootprints: outline Üllar Rannik University of Helsinki
Footprints: outline Üllar Rannik University of Helsinki -Concept of footprint and definitions -Analytical footprint models -Model by Korman and Meixner -Footprints for fluxes vs. concentrations -Footprints
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationTurbulent kinetic energy budgets in a model canopy: comparisons between LES and wind-tunnel experiments
DOI.7/s65-7-949- ORIGINAL ARTICLE Turbulent kinetic energy budgets in a model canopy: comparisons between LES and wind-tunnel experiments Wusi Yue Charles Meneveau Marc B. Parlange Weihong Zhu Hyung Suk
More informationLARGE-EDDY SIMULATIONS OF A WIND TURBINE WAKE ABOVE A FOREST
LARGE-EDDY SIMULATIONS OF A WIND TURBINE WAKE ABOVE A FOREST Josef Schröttle a,c, Zbigniew Piotrowski b, Thomas Gerz a, Antonia Englberger a, Andreas Dörnbrack a a Institute for Atmospheric Physics, German
More informationThe Atmospheric Boundary Layer. The Surface Energy Balance (9.2)
The Atmospheric Boundary Layer Turbulence (9.1) The Surface Energy Balance (9.2) Vertical Structure (9.3) Evolution (9.4) Special Effects (9.5) The Boundary Layer in Context (9.6) Fair Weather over Land
More informationMULTIDIMENSIONAL TURBULENCE SPECTRA - STATISTICAL ANALYSIS OF TURBULENT VORTICES
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 MULTIDIMENSIONAL TURBULENCE SPECTRA - STATISTICAL ANALYSIS OF TURBULENT VORTICES
More informationPrototype Instabilities
Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies
More informationSurface waves and spatially coherent structures in the near-surface layer of the ocean
Surface waves and spatially coherent structures in the near-surface layer of the ocean Alexander Soloviev 1,2 and Cayla Dean 1 1 Oceanographic Center, Nova Southeastern University, USA 2 Rosenstiel School
More informationChapter (3) TURBULENCE KINETIC ENERGY
Chapter (3) TURBULENCE KINETIC ENERGY 3.1 The TKE budget Derivation : The definition of TKE presented is TKE/m= e = 0.5 ( u 2 + v 2 + w 2 ). we recognize immediately that TKE/m is nothing more than the
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationWALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION
WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION Joongnyon Kim, Kyoungyoun Kim, Hyung Jin Sung Department of Mechanical Engineering, Korea Advanced Institute of Science
More informationWind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the
1 2 Wind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the canopy. 3 Next we discuss turbulence in the canopy.
More information6.4 EXPERIMENTAL DETERMINATION OF THE TURBULENT KINETIC ENERGY BUDGET WITHIN AND ABOVE AN URBAN CANOPY
6.4 EXPERIMENTAL DETERMINATION OF THE TURBULENT KINETIC ENERGY BUDGET WITHIN AND ABOVE AN URBAN CANOPY Andreas Christen () *, Mathias W. Rotach (), Roland Vogt () () University of Basel, Institute of Meteorology,
More informationThree-dimensional Floquet stability analysis of the wake in cylinder arrays
J. Fluid Mech. (7), vol. 59, pp. 79 88. c 7 Cambridge University Press doi:.7/s78798 Printed in the United Kingdom 79 Three-dimensional Floquet stability analysis of the wake in cylinder arrays N. K.-R.
More informationEvapotranspiration: Theory and Applications
Evapotranspiration: Theory and Applications Lu Zhang ( 张橹 ) CSIRO Land and Water Evaporation: part of our everyday life Evapotranspiration Global Land: P = 800 mm Q = 315 mm E = 485 mm Evapotranspiration
More informationMultichoice (22 x 1 2 % = 11%)
EAS572, The Atmospheric Boundary Layer Final Exam 7 Dec. 2010 Professor: J.D. Wilson Time available: 150 mins Value: 35% Notes: Indices j = (1, 2, 3) are to be interpreted as denoting respectively the
More informationValidation of FIRETEC wind-flows over a canopy and a fuel-break
CSIRO PUBLISHING www.publish.csiro.au/journals/ijwf International Journal of Wildland Fire 29, 18, 775 79 Validation of FIRETEC wind-flows over a canopy and a fuel-break François Pimont A,D, Jean-Luc Dupuy
More information7C.6 MOMENTUM FLUX STRUCTURES AND STATISTICS IN LOW-WIND MARINE SURFACE LAYERS: OBSERVATIONS AND LARGE-EDDY SIMULATIONS
7C.6 MOMENTUM FLUX STRUCTURES AND STATISTICS IN LOW-WIND MARINE SURFACE LAYERS: OBSERVATIONS AND LARGE-EDDY SIMULATIONS Peter P. Sullivan 1, James B. Edson 2, Tihomir Hristov 3, and James C. McWilliams
More informationSymmetry of Turbulent Characteristics Inside Urban Intersection
Colloquium FLUID DYNAMICS 2007 Institute of Thermomechanics AS CR, v. v. i., Prague, October 24-26, 2007 p.1 Symmetry of Turbulent Characteristics Inside Urban Intersection Radka Kellnerová 1,2 Zbyněk
More informationPeriodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.
Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.
More informationPart I: Overview of modeling concepts and techniques Part II: Modeling neutrally stratified boundary layer flows
Physical modeling of atmospheric boundary layer flows Part I: Overview of modeling concepts and techniques Part II: Modeling neutrally stratified boundary layer flows Outline Evgeni Fedorovich School of
More informationShear instabilities. Chapter Energetics of shear instabilities
Chapter 7 Shear instabilities In this final Chapter, we continue our study of the stability of fluid flows by looking at another very common source of instability, shear. By definition, shear occurs whenever
More informationLIDAR OBSERVATIONS OF FINE-SCALE ATMOSPHERIC GRAVITY WAVES IN THE NOCTURNAL BOUNDARY LAYER ABOVE AN ORCHARD CANOPY
LIDAR OBSERVATIONS OF FINE-SCALE ATMOSPHERIC GRAVITY WAVES IN THE NOCTURNAL BOUNDARY LAYER ABOVE AN ORCHARD CANOPY Tyson N. Randall, Elizabeth R. Jachens, Shane D. Mayor California State University, Chico
More informationApril 1990 T. Watanabe and J. Kondo 227. The Influence of Canopy Structure and Density. upon the Mixing Length within and above Vegetation
April 1990 T. Watanabe and J. Kondo 227 The Influence of Canopy Structure and Density upon the Mixing Length within and above Vegetation By Tsutomu Watanabe and Junsei Kondo Geophysical Institute, Tohoku
More informationEffects of Free-Stream Vorticity on the Blasius Boundary Layer
17 th Australasian Fluid Mechanics Conference Auckland, New Zealand 5-9 December 2010 Effects of Free-Stream Vorticity on the Boundary Layer D.A. Pook, J.H. Watmuff School of Aerospace, Mechanical & Manufacturing
More informationLand/Atmosphere Interface: Importance to Global Change
Land/Atmosphere Interface: Importance to Global Change Chuixiang Yi School of Earth and Environmental Sciences Queens College, City University of New York Outline Land/atmosphere interface Fundamental
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More informationTurbulence Laboratory
Objective: CE 319F Elementary Mechanics of Fluids Department of Civil, Architectural and Environmental Engineering The University of Texas at Austin Turbulence Laboratory The objective of this laboratory
More informationPIV STUDY OF LONGITUDINAL VORTICES IN A TURBULENT BOUNDARY LAYER FLOW
ICAS CONGRESS PIV STUDY OF LONGITUDINAL VORTICES IN A TURBULENT BOUNDARY LAYER FLOW G. M. Di Cicca Department of Aerospace Engineering, Politecnico di Torino C.so Duca degli Abruzzi, 4 - I 19 Torino, ITALY
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationTilting Shear Layers in Coastal Flows
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Tilting Shear Layers in Coastal Flows Karl R. Helfrich Department of Physical Oceanography, MS-21 Woods Hole Oceanographic
More informationChuichi Arakawa Graduate School of Interdisciplinary Information Studies, the University of Tokyo. Chuichi Arakawa
Direct Numerical Simulations of Fundamental Turbulent Flows with the Largest Grid Numbers in the World and its Application of Modeling for Engineering Turbulent Flows Project Representative Chuichi Arakawa
More information