Upstream stagnation points in stratified flow past obstacles

An experimental study has been made of stagnation points and flow splitting on the upstream side of obstacles in uniformly stratified flow. A range from small to large values of Nh/U (where N is the buoyancy frequency, h_m is the maximum obstacle height and U is the undisturbed fluid velocity) has been covered, for three obstacle shapes which are, respectively, axisymmetric, and elongated in the across-stream and in the downstream directions. Upstream stagnation for the first two of these models does not occur until Nh_m/U > 1.05, where it occurs at z ≈ h_m/2. On the central line below this point the flow descends and diverges, and we term this ‘flow splitting’. For the third model (elongated in the downstream direction), stagnation upstream first occurs at Nh_m/U ≈ 1.43, at z ≈ 0. Results for this obstacle are not consistent with the ‘Sheppard criterion’, and this upstream flow stagnation is not apparently related to lee wave overturning, in contrast to flow over two-dimensional obstacles.     

Stratified flow over three-dimensional topography

In order to investigate flows over topography in an atmospheric context, we have studied experimentally the wake structure of axi-symmetric Gaussian obstacles towed through a linearly stratified fluid. Three dimensionless parameters govern the flow dynamics: F, the Froude number based on the topography height h; Re, the Reynolds number and the aspect ratio r = h/L, where L is the topography horizontal scale. Two-dimensional (2-D), saturated lee wave (SLW) and three-dimensional (3-D) regimes, as defined in Chomaz et al. (1993), are found to be functions of F and r only (Fig. 1) as soon as Re is larger than Re_c ≈ 2000. For F < 0.7 the flow goes around the obstacle and the motion in the wake is quasi-two-dimensional. This 2-D layer is topped by a region affected by lee wave motions with amplitude increasing with r and F. For 0.7 < F < 1/r, the flow is entirely dominated by a lee wave of saturated amplitude which suppresses the separation of the boundary layer from the obstacle. Above the critical value 1/r, the lee wave amplitude decreases with F and a recirculating zone appears behind the obstacle. Simultaneously, coherent large-scale vortices start to be shed periodically from the wake at a Strouhal number which decreases as 1/F until it reaches its neutral asymptotic value.     

Mesoscale vortex shedding from large islands: A comparison with laboratory experiments of rotating stratified flows

Summary Mesoscale vortex shedding in the wake of large islands has been analysed with respect to the influence of stratification and background rotation. Comparison with results from laboratory experiments on rotating-stratified flows past isolated obstacles based on Rossby and Froude numbers have shown that forFr<0.4 vortex shedding can be expected for moderate rotation rates. Rossby numbers for cases of atmospheric vortex streets were found to be much larger than in laboratory experiments. Depending on the Froude number different flow pattern in the wake of three-dimensional obstacles have been identified in both, laboratory experiments and atmospheric observations.With 5 Figures     

Quasi-geostrophic stratified flow over isolated finite amplitude topography

The quasi-geostrophic response of a stratified stream incident upon isolated finite amplitude topography on a f-plane is examined in the limit of a Boussinesq, incompressible, inviscid fluid. Compact solutions are derived subject to the following stipulations: uniform upstream velocity and stratification, a circular obstacle and an entirely isentropic/isopycnic lower surface.It is shown that for a semi-infinite flow domain the criterion for Taylor cap formation (i.e., a region of closed streamlines) is . However, for the isentropic lower boundary condition the solutions exist (i.e., have physical validity) only if R₀F⁻¹ < 0.5. (Here R₀ and F refer to the Rossby and Froude numbers defined respectively in terms of the mountain half-width and height.) Also considered are the modifications both to the flow response and to the foregoing existence criterion that are induced by the introduction of an upstream profile comprising two layers of uniform but different stratification. In addition, the relationship of the derived solutions to the results obtained in previous studies is explored, and in particular an outline is given of the impact of adopting the ‘traditional’ simplified lower boundary condition.     

Effects of shear and sharp gradients in static stability on two-dimensional flow over an isolated mountain ridge

Summary ?We have investigated the effects of shear and sharp gradients in static stability and demonstrated how a mountain wave and its associated surface winds can be strongly influenced. Linear theory for two-dimensional, nonrotating stratified flow over an isolated mountain ridge with positive shear and constant static stability shows that the horizontal wind speeds on both the lee and upslope surfaces are suppressed by positive shear. The critical F(=U/Nh where U is the basic wind speed, N the Brunt-Vaisala frequency, and h the mountain height) for the occurrence of wave breaking decreases when the strength of the positive shear increases, while the location for the wave-induced critical level is higher in cases with larger positive shear. The linear theory is then verified by a series of systematic nonlinear numerical experiments. Four different flow regimes are found for positive shear flow over a two-dimensional mountain. The values of critical F which separate the flow regimes are lower when the strength of the positive shear is larger. The location of stagnation aloft from numerical simulations is found to be quite consistent with those predicted by linear theory. We calculate the strongest horizontal wind speed on the lee surface (U _max), the smallest horizontal wind speed on the upslope surface (U _min), the reflection (Ref), and the transmission (Tran) coefficients for different combinations of the stability ratio between the upper and lower layers (i.e. and z ₁ (interface height) in a two-layer atmosphere from linear analytical solutions. Both Ref and Tran are found to be functions of log() but not the interface height (z ₁). Ref is larger when is much different from 1, no matter whether it is larger or smaller than 1. However, Tran decreases when log() increases and approaches 0 when log() is large. The magnitude of the largest U _max (smallest U _min) increases (decreases) as the absolute value of log() increases. It is found that the largest U _max occurs when the nondimensional z ₁ is near for cases with a less stable upper layer or when z ₁ is near for cases with a more stable upper layer. These results are confirmed by nonlinear numerical simulations. We find that linear theory is very useful in qualitative analysis of the possibility of high-drag state for different stability profiles. The location of stagnation aloft in a two-layer atmosphere from numerical simulations agrees very well with those predicted by linear theory. The above findings are applied to investigate the Boulder severe downslope windstorm of 11 January 1972. We find that the windstorm cannot develop if the near mountain-top inversion is located at a higher altitude (e.g.,  km). However, if there exists a less stable layer right below the tropopause, the windstorm can develop in the absence of a low-level inversion. These results indicate the importance of partial reflection due to the structured atmosphere in influencing the possibility of severe downslope windstorms, although partial reflection may not be the responsible mechanism for the generation of windstorms. Received September 25, 1999/Revised February 9, 2000     

Wind profile constants in a neutral atmospheric boundary layer over complex terrain

The roughness height z ₀ and the zero-plane displacement height d ₀ were determined for a region of complex terrain in the Pre-Alps of Switzerland. This region is characterized by hills of the order of 100 m above the valley elevations, and by distances between ridges of the order of 1 km; it lies about 20 to 30 km north from the Alps. The experimental data were obtained from radiosonde observations under near neutral conditions. The analysis was based on the assumption of a logarithmic profile for the mean horizontal wind existing over one half of the boundary layer. The resulting (z ₀/h) and (d ₀/h) (where h is the mean height of the obstacles) were found to be in reasonable agreement with available relationships in terms of placement density and shape factor of the obstacles, which were obtained in previous experiments with h-scales 2 to 4 orders of magnitude smaller than the present ones.     

On atmospheric vortex streets in the wake of large islands

Summary Results from laboratory experiments on stably stratified flows over three-dimensional obstacles are related to atmospheric vortex streets formed in the lee of large islands. A quasi horizontal flow around the island can be expected if stable stratification favours the formation of a so-called dividing streamline below the islands top. The subsequent shedding of vortices with vertical axis from islands may then be due to viscous boundary layer separation, but also other possible mechanisms of vortex shedding in stably stratified flows are discussed.With 6 Figures     

Stratified airflow over one or two hills

Methods for calculating, interpolating and idealising air flow in complex terrain are reviewed. Then the general structure of stratified airflow over a single hill of height H and length L₁ is studied in detail and shown to be determined by the upwind velocity profile, the magnitude of a characteristic Froude number and the dimensions of the hill. Let N(L₁) be the buoyancy frequency upwind at a height L₁, and u_* and U₀ be the upwind friction and mean velocity respectively; then the flow is effectively neutral if u^*/NL₁>1. But if u_*/NL₁>1 and u₀/NL₁>1, the stratification is weak enough to affect the upwind turbulence and velocity profile but not the dynamics of the flow over the hill. If U₀/NL₁ <1 but U₀/NH>1 the buoyancy forces are strong enough to affect the mean flow over the hill but not strong enough to prevent it passing over the top. In this regime the flow is very sensitive to the form of the upwind temperature profile. If U₀/NH>1, much of the flow passes round the hill. A similar classification, with different flow patterns, is appropriate for unstably stratified flows. When the wind is weak enough, local slope winds can dominate. Results from the analysis of these different regimes are described and compared with laboratory of field measurements where possible. It is shown how some of these results can be extended to groups of hills.     

Interaction between atmospheric flows and obstacles: Experiments in a rotating channel

We performed an experimental study using scale models in a hydrodynamic rotating channel, concerning the interactions between fluid flows and obstacles of different shapes. The study was meant to analyze the characteristics of the wakes observed on the lee side of quasi-bidimensional obstacles, in a neutral atmosphere.The obstacles were half-cylinders (with aspect ratio 0.87), placed transversally on the channel bottom and totally submerged in the fluid. We call them quasi-bidimensional since their width was a little smaller than the channel width, thus allowing the flow to partially go round their edges.The simulations were performed in the rotating hydraulic channel of ICG-CNR in Turin, and included various conditions of rotation period and flow speed. An interesting behaviour of the wakes was found on the lee side of subsynoptic-scale obstacles, modelled in conditions of Reynolds-Rossby similitude. More precisely, if a given threshold of flow velocity is exceeded, wake size is constant and is fully determined by the height of the obstacle.     

Wind-tunnel Modelling of Dispersion from a Scalar Area Source in Urban-Like Roughness

A wind-tunnel study was conducted to investigate ventilation of scalars from urban-like geometries at neighbourhood scale by exploring two different geometries a uniform height roughness and a non-uniform height roughness, both with an equal plan and frontal density of λ _p = λ _f = 25%. In both configurations a sub-unit of the idealized urban surface was coated with a thin layer of naphthalene to represent area sources. The naphthalene sublimation method was used to measure directly total area-averaged transport of scalars out of the complex geometries. At the same time, naphthalene vapour concentrations controlled by the turbulent fluxes were detected using a fast Flame Ionisation Detection (FID) technique. This paper describes the novel use of a naphthalene coated surface as an area source in dispersion studies. Particular emphasis was also given to testing whether the concentration measurements were independent of Reynolds number. For low wind speeds, transfer from the naphthalene surface is determined by a combination of forced and natural convection. Compared with a propane point source release, a 25% higher free stream velocity was needed for the naphthalene area source to yield Reynolds-number-independent concentration fields. Ventilation transfer coefficients w _T /U derived from the naphthalene sublimation method showed that, whilst there was enhanced vertical momentum exchange due to obstacle height variability, advection was reduced and dispersion from the source area was not enhanced. Thus, the height variability of a canopy is an important parameter when generalising urban dispersion. Fine resolution concentration measurements in the canopy showed the effect of height variability on dispersion at street scale. Rapid vertical transport in the wake of individual high-rise obstacles was found to generate elevated point-like sources. A Gaussian plume model was used to analyse differences in the downstream plumes. Intensified lateral and vertical plume spread and plume dilution with height was found for the non-uniform height roughness.     

Momentum and Heat Transfer over Urban-like Surfaces

Momentum and heat transfer was examined for the urban-like surfaces used within the Comprehensive Outdoor Scale MOdel (COSMO) experiments. Simultaneous and comparative meteorological measurements were made over a pair of scale models with different block geometries. These data were analyzed to investigate the influence of height variations, obstacle elongation, and packing density, λ _p , of blocks on the aerodynamic properties. In addition, the robustness of theoretical expressions of bulk transfer coefficients for momentum and heat with respect to geometric parameters was examined. Our analyses showed: (1) the theoretical framework for the bulk transfer coefficient for momentum, C _m , and that for heat, C _h , was applicable for homogeneous building arrays, (2) the sensitivity of C _h to the surface geometry was smaller than that of C _m , (3) the transfer coefficients were increased by variations of block heights, but not by elongation of blocks, (4) first-order approximations of C _m and C _h for an array of blocks with two different heights can be made by applying simple theoretical assumptions to include the effects of height variation, and (5) variations of block heights increased the momentum flux significantly, but caused little change in the sensible heat flux. This can be explained by the feedback mechanism of aerodynamic– thermal interaction; aerodynamic mixing decreased both the advective velocity and the vertical temperature gradient.     

The Entrainment Rate for Buoyant Plumes in a Crossflow

We consider large-eddy simulations (LES) of buoyant plumes from a circular source with initial buoyancy flux F ₀ released into a stratified environment with constant buoyancy frequency N and a uniform crossflow with velocity U. We make a systematic comparison of the LES results with the mathematical theory of plumes in a crossflow. We pay particular attention to the limits [(U)\tilde] << 1{\tilde{U}\ll1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}, where [(U)\tilde]=U/(F₀ N)^1/4{\tilde{U}=U/(F_0 N)^{1/4}}, for which analytical results are possible. For [(U)\tilde] >> 1{\tilde{U}\gg 1}, the LES results show good agreement with the well-known two-thirds law for the rise in height of the plume. Sufficiently far above the source, the centreline vertical velocity of the LES plumes is consistent with the analytical z ^−1/3 and z ^−1/2 scalings for respectively [(U)\tilde] << 1{\tilde{U}\ll 1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}. In the general case, where the entrainment is assumed to be the sum of the contributions from the horizontal and vertical velocity components, we find that the discrepancy between the LES data and numerical solutions of the plume equations is largest for [(U)\tilde]=O(1){\tilde{U}=O(1)}. We propose a modified additive entrainment assumption in which the contributions from the horizontal and vertical velocity components are not equally weighted. We test this against observations of the plume generated by the Buncefield fire in the U.K. in December 2005 and find that the results compare favourably. We also show that the oscillations of the plume as it settles down to its final rise height may be attenuated by the radiation of gravity waves. For [(U)\tilde] << 1{\tilde{U}\ll 1} the oscillations decay rapidly due to the transport of energy away from the plume by gravity waves. For ${\tilde{U}>rsim 1}${\tilde{U}>rsim 1} the gravity waves travel in the same direction and at the same speed as the flow. In this case, the oscillations of the plume do not decay greatly by radiation of gravity waves.     

Three-dimensional simulations and wind-tunnel experiments on airflow over isolated forest stands

Prediction of windthrow risk to individual or groups of retained trees in harvested stands requires an improved understanding of canopy airflow dynamics. Large-eddy simulations were used to simulate wind-tunnel experiments in two and three dimensions to compare with observations for model validation and to address parameter space considerations for the design of subsequent retention pattern experiments. The three-dimensional simulations were similar to the observed wind-tunnel data for the statistical profiles for but there were greater differences in skewness and kurtosis. These results were obtained using a common leaf-area drag formulation without either skin friction or speed dependent drag that enables scaling with U ₀ (ambient wind speed) and h (height of the canopy). This scaling results in a single non-dimensional parameter h/h _c where h _c (x, y, z) is the momentum range resulting from the canopy drag. The validity of the model scaling was tested using two-dimensional simulations. The irrotational component of the flow (potential flow) was found to be important when defining vertical domain limitations and has significant implications for time dependent flow (i.e. turbulent conditions) when considering retention pattern design. The sudden onset of drag associated with the isolated stand presents some unexpected challenges. The horizontal scales of the shearing instabilities were simulated in two dimensions and found to range between 2h for early times to 7h for later times. The early-time horizontal scales are in the range of logical retention pattern scales and as such need to be taken into account as part of the parameter space, i.e. a range of retention pattern lengths need consideration.     

Effects of Topographical Slope Angle and Atmospheric Stratification on Surface-Layer Turbulence

The effect of topographical slope angle and atmospheric stratification on turbulence intensities in the unstably stratified surface layer have been parameterized using observations obtained from a three-dimensional sonic anemometer installed at 8 m height above the ground at the Seoul National University (SNU) campus site in Korea for the years 1999–2001. Winds obtained from the sonic anemometer are analyzed according to the mean wind direction, since the topographical slope angle changes significantly along the azimuthal direction. The effects of the topographical slope angle and atmospheric stratification on surface-layer turbulence intensity are examined with these data. It is found that both the friction velocity and the variance for each component of wind normalized by the mean wind speed decrease with increase of the topographical slope angle, having a maximum decreasing rate at very unstable stratification. The decreasing rate of the normalized friction velocity (u _* /U) is found to be much larger than that of the turbulence intensity of each wind component due to the reduction of wind shear with increase in slope angle under unstable stratification. The decreasing rate of the w component of turbulence intensity (σ _w /U) is the smallest over the downslope surface whereas that of the u component (σ _u /U) has a minimum over the upslope surface. Consequently, σ _w /u _* has a maximum increasing rate with increase in slope angle for the downslope wind, whereas σ _u /u _* has its maximum for the upslope wind. The sloping terrain is found to reduce both the friction velocity and turbulence intensity compared with those on a flat surface. However, the reduction of the friction velocity over the sloping terrain is larger than that of the turbulence intensity, thereby enhancing the turbulence intensity normalized by the friction velocity over sloping terrain compared with that over a flat surface.     

Effect of atmospheric stability on the growth of surface gravity waves

In this paper we study the effect of atmospheric stability on the growth of surface gravity waves. To that end we numerically solved the Taylor-Goldstein equation for wind profiles which deviate from a logarithmic form because stratification affects the turbulent momentum transport. Using Charnock's relation for the roughness height z ₀ of the wind profile, it is argued that the growth rate of the wave depends on the dimensionless phase velocity c/u _* (where u _* is the friction velocity) and a measure of the effect of atmospheric stability, namely the dimensionless Obukhov length gL/u _* ², whereas it only depends weakly on gz _t /u _* ² (where z _t is the roughness height of the temperature profile). Remarkably for a given value of u _* /c, the growth rate is larger for a stable stratification (L > 0) than for an unstable one (L < 0). We explain why this is the case. If, on the other hand, one considers the growth rate as a function of c/U ₁₀ (where U ₁₀ is the windspeed at 10 m), the situation reverses for c/U ₁₀ < 1. For practical application in wave prediction models, we propose a new parameterization of the growth rate of the waves which is an improvement of the Snyder et al. (1981) proposal because the effect of stability is taken into account.     

A simple model of drainage flow on a slope

The concept of a cold air ‘Parcel’ is introduced for describing the bulk properties of drainage flow. By means of a model based on the momentum and sensible heat transports under calm conditions, the thickness h and velocity u of the Parcel are derived in simple forms. It is shown that h and u correspond to the inversion height and maximum velocity of actual drainage flow. The governing parameters for h and u are the length and vertical drop of the slope, potential temperature difference between the ambient atmosphere and the Parcel, aerodynamic condition of the slope surface expressed by the mean bulk coefficients, and ambient stability. The mean bulk coefficients depend on the roughness lengths for the velocity and potential temperature profiles and are decreasing functions of the slope length. The Parcel Model agrees qualitatively with Manins and Sawford's (1979) model under neutral ambient stratification. But agreement is not so good under stable conditions. The thickness and velocity of drainage flow predicted by the Parcel Model agree with observations on slopes several tens of meters to several hundred kilometers long.     

Two-dimensional asymmetric vortex merger: merger dynamics and critical merger distance

Two-dimensional asymmetric merger of two like-signed vorticity monopoles with different sizes and vorticities is examined by combining simplified analytical models and contour dynamics experiments. The model results can capture the key dynamics and hence allow the prediction of the critical merger distance in a number of the situations. The models ignore deformation of one of the two vortices, replacing it with a point vortex, and employ a corotating frame of reference with a rotation rate estimated by point vortices. Thus, the two vortex problem becomes two separate problems of a single vortex in a background shear flow. Vortex merger is found to happen when the vortex cannot resist the background shear flow. Vortex merger and merging processes depend on the centroid distance d, the circulation ratio, (q_i and r_i are the vorticity and radius, respectively) and initial conditions. In the lowest order, the background flow is approximated by a uniform shear field, and the behavior of an elliptical vortex can be described by the Kida (1981) equation supplemented with one describing the time evolution of the centroid distance. This model reveals that merger takes place because the natural rotation of an elliptical vortex is overcome by the background uniform shear flow; the ellipse inversely rotates and is drawn out by the background straining field. The vortex deformation in a background flow field induces an inward flow at the position of the other vortex; as a result, the centroid distance decreases and two vortices merge. The critical merger distance from this model agrees quite well with the results from contour dynamics experiments for two vortices. Inclusion of higher order non-uniform shear in the background flow extends the critical merger distance, which gives almost perfect estimates for the experiment. In the non-uniform shear flow, partial merger occurs, where the vortex sheds off a filament, but the remaining part of the vortex resumes its natural rotation.     

The Influence of Hilly Terrain on Canopy-Atmosphere Carbon Dioxide Exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO₂ (N_ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO₂ sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO₂ exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO₂. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N_ee = F_c, the vertical turbulent flux of CO₂. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F_c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F_c-at canopy top differs by more than 20% compared to the flat-terrain case.     

A variable K planetary boundary-layer model

The steady-state, homogeneous and barotropic equations of motion within the planetary boundary layer are solved with the assumption that the coefficient of eddy viscosity varies as K(Z) = K _O(1–Z/h) ^p , where h is the height of the bounday layer and p is a parameter which depends on atmospheric stability. The solutions compare favourably with observed velocity profiles based on the Wangara data.