Dominant features in three-dimensional turbulence structure: comparison of non-uniform accelerating and decelerating flows

The results are presented from an experimental study to investigate three-dimensional turbulence structure profiles, including turbulence intensity and Reynolds stress, of different non-uniform open channel flows over smooth bed in subcritical flow regime. In the analysis, the uniform flow profiles have been used to compare with those of the non-uniform flows to investigate their time-averaged spatial flow turbulence structure characteristics. The measured non-uniform velocity profiles are used to verify the von Karman constant κ and to estimate sets of log-law integration constant B_r and wake parameter П, where their findings are also compared with values from previous studies. From κ, B_r and П findings, it has been found that the log-wake law can sufficiently represent the non-uniform flow in its non-modified form, and all κ, B_r and П follow universal rules for different bed roughness conditions. The non-uniform flow experiments also show that both the turbulence intensity and Reynolds stress are governed well by exponential pressure gradient parameter β equations. Their exponential constants are described by quadratic functions in the investigated β range. Through this experimental study, it has been observed that the decelerating flow shows higher empirical constants, in both the turbulence intensity and Reynolds stress compared to the accelerating flow. The decelerating flow also has stronger dominance to determine the flow non-uniformity, because it presents higher Reynolds stress profile than uniform flow, whereas the accelerating flow does not.     

Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

Numerical simulation and optimization study of the wind flow through a porous fence

Three turbulence closure models (RNG k-ε, SST k-ω and RSM) were used to investigate the flow characteristics around a two-dimensional isolated porous fence. The comparison between the numerical results and the experimental measurements indicated that RSM model shows a better performance than the other two models. The aim of this paper is to accurately and efficiently determine the optimum porosity that attain the best shelter effect of the wind fence in the near wake region (0–4h_b) and in the far wake region (4h_b–10h_b) respectively, where h_b is the height of the fence. The gradient algorithm was adopted as the optimization algorithm and the RSM model was used to model turbulent features of the flow. The shelter effect was parameterized by the peak velocity ratio involving velocity and turbulence. The objective was to reduce the peak velocity ratio in the near or far wake region by changing the design variable porosity (?) of the fence, which ranged between 2 and 60%. The results revealed that a porosity of 10.2% was found as the optimum value giving rise to the best shelter effect in the near wake region, and ? = 22.1% was determined in the case of the far wake region. In addition, based on the proposed optimization method, it is found that the recirculating bubble behind the fence can only be detected when ? < 29.9%.     

Characteristics of flow structure of free-surface flow in a partly obstructed open channel with vegetation patch

Free-surface flows over patchy vegetation are common in aquatic environments. In this study, the hydrodynamics of free-surface flow in a rectangular channel with a bed of rigid vegetation-like cylinders occupying half of the channel bed was investigated and interpreted by means of laboratory experiments and numerical simulations. The channel configurations have low width-to-depth aspect ratio (1.235 and 2.153). Experimental results show that the adjustment length for the flow to be fully developed through the vegetation patch in the present study is shorter than observed for large-aspect-ratio channels in other studies. Outside the lateral edge of the vegetation patch, negative velocity gradient (\(\partial \overline{u}/\partial z < 0\)) and a local velocity maximum are observed in the vertical profile of the longitudinal velocity in the near-bed region, corresponding to the negative Reynolds stress (\(- \overline{{u^{\prime}w^{\prime}}} < 0\)) at the same location. Assuming coherent vortices to be the dominant factor influencing the mean flow field, an improved Spalart–Allmaras turbulence model is developed. The model improvement is based on an enhanced turbulence length scale accounting for coherent vortices due to the effect of the porous vegetation canopy and channel bed. This particular flow characteristic is more profound in the case of high vegetation density due to the stronger momentum exchange of horizontal coherent vortices. Numerical simulations confirmed the local maximum velocity and negative gradient in the velocity profile due to the presence of vegetation and bed friction. This in turn supports the physical interpretation of the flow processes in the partly obstructed channel with vegetation patch. In addition, the vertical profile of the longitudinal velocity can also be explained by the vertical behavior of the horizontal coherent vortices based on a theoretical argument.     

Mixing at a sediment concentration interface in turbulent open channel flow

In this work we address the role of turbulence on mixing of clear layer of fluid with sediment-laden layer of fluid at a sediment concentration interface. This process can be conceived as the entrainment of sediment-free fluid into the sediment-laden layer, or alternatively, as the transport of sediment into the top sediment-free flow. This process is governed by four parameters—Reynolds number of the flow \(Re_\tau\), non-dimensional settling velocity of the sediment (proxy for sediment size) \(\tilde{V}\), Richardson number \(Ri_\tau\) and Schmidt number Sc. For this work we have performed direct numerical simulations for fixed Reynolds and Schmidt numbers while varying the values of Richardson number and particle settling velocity. In the simple model considered here, the flow’s momentum and turbulence pre-exists over the entire layer of fluid, while the sediment is initially confined to a layer close to the bed. Mixing of sediment-free fluid with the sediment-laden layer is associated primarily with upward transport of sediment and buoyancy. There is no simultaneous upward transport of fluid momentum and turbulence into the sediment-free fluid layer, which is already in motion and turbulent. The analysis performed shows that the ability of the flow to transport a given sediment size decreases with the distance from the bottom, and thus only fine enough sediment particles are transported across the sediment concentration interface. For these cases, the concentration profiles evolve to a final steady state in good agreement with the well-known Rouse profile. The approach towards the Rouse profile happens through a transient self-similar state. This behavior of the flow is not seen for larger particles. Detailed analysis of the three dimensional structure of the sediment concentration interface shows the mechanisms by which sediment particles are lifted up by tongues of sediment-laden fluid with positive correlation between vertical velocity and sediment concentration. Finally, the mixing ability of the flow is addressed by monitoring the time evolution of the center of mass of the sediment-laden layer and the vertical location of the sediment-free/sediment-laden interface.     

Model of turbulence penetration into a suspension layer on a sediment bed and effect on vertical solute transfer

The vertical diffusional mass (solute) transfer through a suspended sediment layer, e.g. at the bottom of a lake, reservoir or estuary, by the propagation of velocity fluctuations from above was investigated. The attenuation of the velocity fluctuations in the suspension layer and the associated effect on solute transfer through the suspension layer was simulated. To represent large eddies traveling downstream in water over a high-concentration suspended sediment layer, a streamwise velocity fluctuation moving in downstream direction was imposed along the upper boundary of the suspension layer. Velocity fluctuations and downstream velocity were normalized by the shearvelocity (U_*) at the top of the suspension layer. Streamwise and vertical velocity components inside the suspension layer, were obtained from the 2-D continuity and the Navier–Stokes equations. The persistence of turbulence with depth—as it penetrates from the overlying water into the suspension layer—was found to depend on its amplitude, its period, and on the apparent viscosity of the suspension. The turbulence was found to propagate efficiently into the suspension layer when its frequency is low, and the apparent viscosity of the suspension is high. Effects on vertical mass transfer were parameterized by penetration depth and effective diffusion coefficient, and related to apparent viscosity of the suspension, Schmidt number and shear velocity on top of the suspension layer. The enhancement of turbulence penetration by viscosity is similar to the flow near an oscillating flat plate (Stokes’ second problem), but is opposite to turbulence penetration into a stationary porous and permeable sediment bed. The information is applicable to water quality modeling mear the sediment/water interface of lakes, river impoundments and estuaries.     

Turbulent velocity profile in fully-developed open channel flows

The determination of velocity profile in turbulent narrow open channels is a difficult task due to the significant effects of the anisotropic turbulence that involve the Prandtl’s second type of secondary flow occurring in the cross section. With these currents the maximum velocity appears below the free surface that is called dip phenomenon. The well-known logarithmic law describes the velocity distribution in the inner region of the turbulent boundary layer but it is not adapted to define the velocity profile in the outer region of narrow channels. This paper relies on an analysis of the Navier–Stokes equations and yields a new formulation of the vertical velocity profile in the center region of steady, fully developed turbulent flows in open channels. This formulation is able to predict time averaged primary velocity in the outer region of the turbulent boundary layer for both narrow and wide open channels. The proposed law is based on the knowledge of the aspect ratio and involves a parameter C_Ar depending on the position of the maximum velocity (ξ_dip). ξ_dip may be derived, either from measurements or from an empirical equation given in this paper. A wide range of longitudinal velocity profile data for narrow open channels has been used for validating the model. The agreement between the measured and the computed velocities is rather good, despite the simplification used.     

Two-phase flow insights into open-channel flows with suspended particles of different densities

The effect of particle density on the turbulent open-channel flow carrying dilute particle suspensions is investigated using two specific gravities and three concentrations of solid particles. The particles, identical in size and similar in shape, were natural sand and a neutrally buoyant plastic. The particles were fully suspended, and formed no particle streaks on the channel’s bed. Accordingly, the changes in the flow are attributed to the interactions between suspended particles and flow turbulence structures. Measurements were obtained by means of image velocimetry enabling simultaneous, but distinct, measurement of liquid and particle velocities. The experimental results show that, irrespective of particle specific gravity, particle suspension influences bulk velocity of flow and the Kármán coefficient, while friction velocity essentially remains constant. The results also show that particles in suspension modify local water turbulence over the flow depth, but in ways not accurately predicted using the customary parameters for characterizing turbulence modification.     

Turbulence structure and fluid–particle interaction in sediment-laden flows over developing sand dunes

In this study, we conducted the fluid–particle simultaneous measurements in order to reveal the fluid–particle interaction mechanism in the developing stages of sand dunes. We measured the instantaneous velocities over the growing sand dunes by the use of both the discriminator particle-image velocimetory (D-PIV) and the discriminator particle-tracking velocimetory (D-PTV). In this D-PIV&PTV system, fluid tracers and sediment particles are discriminated accurately by their occupied particle sizes, and thus, the particle velocity U _p and fluid velocity U _f can be measured simultaneously. It was found from the present measurements that turbulence intensities in the trough area became larger after the generation of the reverse flow, i.e., U _f < 0. There are two kinds of coherent eddies behind dunes, i.e., one is convected along the shear layer behind dune and the other is lifted up from the reattachment point due to the kolk-boil vortex.     

Laboratory measurements and multi-block numerical simulations of the mean flow and turbulence in the non-aerated skimming flow region of steep stepped spillways

We present and discuss the results of a comprehensive study addressing the non-aerated region of the skimming flow in steep stepped spillways. Although flows in stepped spillways are usually characterized by high air concentrations concomitant with high rates of energy dissipation, the non-aerated region becomes important in small dams and/or spillways with high specific discharges. A relatively large physical model of such spillway was used to acquire data on flow velocities and water levels and, then, well-resolved numerical simulations were performed with a commercial code to reproduce those experimental conditions. The numerical runs benefited from the ability of using multi-block grids in a Cartesian coordinate system, from capturing the free surface with the TruVOF method embedded in the code, and from the use of two turbulence models: the k-e{k{-}\varepsilon} and the RNGk-e{k{-}\varepsilon} models. Numerical results are in good agreement with the experimental data corresponding to three volumetric flow rates in terms of the time-averaged velocities measured at diverse steps in the spillway, and they are in very satisfactory agreement for water levels along the spillway. In addition, the numerical results provide information on the turbulence statistics of the flow. This work also discusses important aspects of the flow, such as the values of the exponents of the power-law velocity profiles, and the characteristics of the development of the boundary layer in the spillway.     

Impact of consistent boundary layer mixing approaches between NAM and CMAQ

Discrepancies in grid structure, dynamics and physics packages in the offline coupled NWS/NCEP NAM meteorological model with the U.S. Environmental Protection Agency Community Multiscale Air Quality (CMAQ) model can give rise to inconsistencies. This study investigates the use of three vertical mixing schemes to drive chemistry tracers in the National Air Quality Forecast Capability (NAQFC). The three schemes evaluated in this study represent various degrees of coupling to improve the commonality in turbulence parameterization between the meteorological and chemistry models. The methods tested include: (1) using NAM predicted TKE-based planetary boundary height, h, as the prime parameter to derive CMAQ vertical diffusivity; (2) using the NAM mixed layer depth to determine h and then proceeding as in (1); and (3) using NAM predicted vertical diffusivity directly to parameterize turbulence mixing within CMAQ. A two week period with elevated surface O₃ concentrations during the summer 2006 has been selected to test these schemes in a sensitivity study. The study results are verified and evaluated using the EPA AIRNow monitoring network and other ozonesonde data. The third method is preferred a priori as it represents the tightest coupling option studied in this work for turbulent mixing processes between the meteorological and air quality models. It was found to accurately reproduce the upper bounds of turbulent mixing and provide the best agreement between predicted h and ozonesonde observed relative humidity profile inferred h for sites investigated in this study. However, this did not translate into the best agreement in surface O₃ concentrations. Overall verification results during the test period of two weeks in August 2006, did not show superiority of this method over the other 2 methods in all regions of the continental U.S. Further efforts in model improvement for the parameterizations of turbulent mixing and other surface O₃ forecast related processes are warranted.     

Compound channel flow with a longitudinal transition in hydraulic roughness over the floodplains

Flows in a compound open-channel (two-stage geometry with a main channel and adjacent floodplains) with a longitudinal transition in roughness over the floodplains are experimentally investigated in an 18 m long and 3 m wide flume. Transitions from submerged dense vegetation (meadow) to emergent rigid vegetation (wood) and vice versa are modelled using plastic grass and vertical wooden cylinders. For a given roughness transition, the upstream discharge distribution between main channel and floodplain (called subsections) is also varied, keeping the total flow rate constant. The flows with a roughness transition are compared to flows with a uniformly distributed roughness over the whole length of the flume. Besides the influence of the downstream boundary condition, the longitudinal profiles of water depth are controlled by the upstream discharge distribution. The latter also strongly influences the magnitude of the lateral net mass exchanges between subsections, especially upstream from the roughness transition. Irrespective of flow conditions, the inflection point in the mean velocity profile across the mixing layer is always observed at the interface between subsections. The longitudinal velocity at the main channel/floodplain interface, denoted \(U_{int}\), appeared to be a key parameter for characterising the flows. First, the mean velocity profiles across the mixing layer, normalised using \(U_{int}\), are superimposed irrespective of downstream position, flow depth, floodplain roughness type and lateral mass transfers. However, the profiles of turbulence quantities do not coincide, indicating that the flows are not fully self-similar and that the eddy viscosity assumption is not valid in this case. Second, the depth-averaged turbulent intensities and Reynolds stresses, when scaled by the depth-averaged velocity \(U_{d,int}\) exhibit two plateau values, each related to a roughness type, meadow or wood. Lastly, the same results hold when scaling by \(U_{d,int}\) the depth-averaged lateral flux of momentum due to secondary currents. Turbulence production and magnitude of secondary currents are increased by the presence of emergent rigid elements over the floodplains. The autocorrelation functions show that the length of the coherent structures scales with the mixing layer width for all flow cases. It is suggested that coherent structures tend to a state where the magnitude of velocity fluctuations (of both horizontal vortices and secondary currents) and the spatial extension of the structures are in equilibrium.     

Turbulence and aeration in hydraulic jumps: free-surface fluctuation and integral turbulent scale measurements

In an open channel, a change from a supercritical to subcritical flow is a strong dissipative process called a hydraulic jump. Herein some new measurements of free-surface fluctuations of the impingement perimeter and integral turbulent time and length scales in the roller are presented with a focus on turbulence in hydraulic jumps with a marked roller. The observations highlighted the fluctuating nature of the impingement perimeter in terms of both longitudinal and transverse locations. The results showed further the close link between the production and detachment of large eddies in jump shear layer, and the longitudinal fluctuations of the jump toe. They highlighted the importance of the impingement perimeter as the origin of the developing shear layer and a source of vorticity. The air–water flow measurements emphasised the intense flow aeration. The turbulent velocity distributions presented a shape similar to a wall jet solution with a marked shear layer downstream of the impingement point. The integral turbulent length scale distributions exhibited a monotonic increase with increasing vertical elevation within 0.2 < L_z/d₁ < 0.8 in the shear layer, where L_z is the integral turbulent length scale and d₁ the inflow depth, while the integral turbulent time scales were about two orders of magnitude smaller than the period of impingement position longitudinal oscillations.     

Observation of bioturbation and hyporheic flux in streambeds

In the Elkhorn River, burrows, tubes, and sediment mounds created by invertebrate bioturbation were observed in the exposed streambed and commonly concentrated on the fine-sediment patches, which consist of silt, clay, and organic matter. These invertebrate activities could loosen the thin layer of clogging sediments and result in an increase of pore size in the sediments, leading to greater vertical hydraulic conductivity of the streambed (K _v ). The measurements of the vertical hydraulic gradient across the submerged streambed show that vertical flux in the hyporheic zone can alter directions (upward versus downward) for two locations only a few meters apart. In situ permeameter tests show that streambed K _v in the upper sediment layer is much higher than that in the lower sediment layer, and the calculated K _v in the submerged streambed is consistently greater than that in the clogged sediments around the shorelines of the sand bars. Moreover, a phenomenon of gas bubble release at the water-sediment interface from the subsurface sediments was observed in the groundwater seepage zone where flow velocity is extremely small. The bursting of gas bubbles can potentially break the thin clogging layer of sediments and enhance the vertical hydraulic conductivity of the streambed.     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

Variations of bed elevations due to turbulence around submerged cylinder in sand beds

This paper presents the spatio-temporal variations in bed elevations and the near-bed turbulence statistics over the deformed bed generated around the submerged cylindrical piers embedded vertically on loose sediment bed at a constant flow discharge. Experiments were carried out in a laboratory flume for three blockage ratios in the range of 0.04–0.06 using three different sizes of submerged cylinders individually placed vertically at the centerline of the flume. Clear-water experimental conditions were maintained over the smooth sediment bed surface with a constant flow discharge (\(Q = 0.015\,{\rm m}^3/{\rm sec}\)), thereby giving three different cylinder Reynolds numbers \(Re_{D_c} = \frac{U_mD_c}{\nu }\) (=10200, 12750, 15300) away from the cylinder locations, where \(U_m\) is the maximum mean velocity, \(D_c\) is the cylinder diameter and \(\nu\) is the kinematic viscosity of fluid. Instantaneous sand bed elevations around the cylinders were recorded using a SeaTek 5MHz ultrasonic ranging system of net 24 transducers to estimate bed form migration, and the near-bed velocity data at transducer locations over the stable deformed bed around the pier-like structures were collected using down-looking three-dimensional (3D) Micro-acoustic Doppler velocimeter to estimate the bottom Reynolds shear stresses and the contributions of bursting events to the dominant shear stress component. The flow perturbation generated due to relatively lower flow blockage ratio favored to achieve the stable bed condition more rapidly than the others, and larger upstream scour-depth and deformed areas were noticed for greater flow blockage ratio due to larger cylinder diameter. For larger blockage ratio in the upstream of scour-hole near the bed, occurrences of probabilities of both boundary-ward interactions (Q1 and Q3) were the dominant; whereas in the downstream of the scoured region, occurrences of probabilities of second and third quadrant events (Q2 and Q4) were dominant. On the other hand, for the lower blockage ratio, quadrant (Q2) was dominant over Q4 in the downstream of scour-hole, and in the upstream of scour-hole, quadrant Q4 was the dominant.     

Turbulent flow characteristics and drag over 2-D forward-facing dune shaped structures with two different stoss-side slopes

The present paper explores the characteristics of turbulent flow and drag over two artificial 2-D forward-facing waveform structures with two different stoss side slopes of $50^{\circ }$ and $90^{\circ },$ respectively. Both structures possessed a common slanted lee side slope of $6^{\circ }.$ Flume experiments were conducted at the Fluvial Mechanics Laboratory of Indian Statistical Institute, Kolkata. The velocity data were analyzed to identify the spatial changes in turbulent flow addressing the flow separation region with recirculating eddy, the Reynolds stresses, the turbulent events associated with burst-sweep cycles and the drag over two upstream-facing bedforms for Reynolds number $Re_h=1.44\times 10^5.$ The divergence at the stoss side slope between the two structures revealed significant changes in the mean flow and turbulence. Comparison showed that during the flood-tide condition there was no flow separation region on the gentle lee side of the structure with smaller slope at the stoss side, while for the other structure with vertical stoss side slope a thick flow separation region with recirculating eddy was observed at the gentle lee side just downstream of the crest. The recirculating eddy induced on the lee-side had a strong influence on the resistance that the structure exerts to the flow due to loss of energy through turbulence. In contrast, a great amount of reduction in drag was observed in the case of smaller stoss side sloped structure as there was no flow separation. The quadrant analysis was also used to highlight the turbulent event evolution along the bed form structures under flood-tide conditions.     

Analysis of the velocity field in a large rectangular channel with lateral shockwave

In this work the authors describe the main characteristics of the velocity field of hydraulic jumps in a very large channel where lateral shockwaves occur. Experiments were carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). Extensive flow velocity measurements were investigated in order to have a clearer understanding of both hydraulic jump development and lateral shockwave formation in a very large channel. Eight experiments were performed in a 4m wide rectangular channel; the experiments differed in the inlet Froude number F ₀ and the jump type. Seven tests were carried out with undular jumps and one with a roller jump. The flow velocity and the flow free surface measurements were taken using a two-dimensional Acoustic Doppler Velocimeter (ADV) and an ultrasonic profiler, respectively. The experimental results can be summarized as follow: (i) the formation of well developed lateral shockwaves similar to those of oblique jumps were observed; (ii) the comparison of the experimental and theoretical data shows that the classic shockwave theory is sufficiently confirmed in the analyzed range of Reynolds number, taking into account the experimental errors and the difference between the theoretical and experimental assumptions; (iii) the transversal flow velocity profiles in the recirculating zone show a good agreement with the numerical simulations presented in literature in the case of a separated turbulent boundary layer over a flat plate. This conclusion enables us to confirm the hypothesis that the lateral shockwaves in the channel are the result of a boundary layer which, as observed, forms on the channel sidewalls.     

Transport and Diffusion of Ozone in the Nocturnal and Morning Planetary Boundary Layer of the Phoenix Valley

The evolution of ozone (O ₃) in the nocturnal and morning-transitional planetary boundary layer (PBL) of the Phoenix valley was measured as a part of the `Phoenix Sunrise Experiment 2001' of the U.S. Department of Energy conducted in June 2001. The goal of the field program was to study the transport, distribution and storage of ozone and its precursors in the urban boundary layer over a diurnal cycle. The ground level O <sub>3</sub> as well as mean meteorological variables and turbulence were measured over the entire period, and vertical profiling (using a tethered balloon) was made during the morning transition period. Approximately half of the observational days showed the usual diurnal cycle of high O <sub>3</sub> during the day and low O <sub>3</sub> at night, with nitrogen oxides (NO <sub>x</sub> = NO <sub>2</sub> + NO) showing an out of phase relationship with O <sub>3</sub>. The rest of the days were signified by an anomalous increase of O <sub>3</sub> in the late evening ( 2200 LST), concomitant with a sudden drop of temperature, an enhancement of wind speed and Reynolds stresses, a positive heat flux and a change of wind direction. NO <sub>x</sub> measurements indicated the simultaneous arrival of an `aged' air mass, which was corroborated by the wind predictions of a mesoscale numerical model. In all, the results indicate that the recirculation of O ₃ rich air masses is responsible for the said high-O ₃ events. Such air masses are produced during the transport of O ₃ precursors by the upslope flow toward mountainous suburbs during the day, and they return back to the city at night via downslope winds (i.e. mountain breeze). The corresponding flow patterns, and hence the high-O ₃ events, are determined by background meteorological conditions. The vertical profiling of O ₃ and flow variables during the morning transition points to a myriad of transport, mixing and chemical processes that determine the fate of tropospheric O ₃. How well such processes are incorporated and resolved in predictive O ₃ models should determine the accuracy of their predictions.