Showing papers in "Environmental Fluid Mechanics in 2009"

Turbulent air–water flows in hydraulic structures: Dynamic similarity and scale effects

Abstract: In hydraulic structures, free-surface aeration is commonly observed: i.e., the white waters. The air bubble entrainment may be localised (hydraulic jumps, plunging jets) or continuous along an interface (water jets, chutes). Despite recent advances, there are some basic concerns about the extrapolation of laboratory results to large size prototype structures. Herein the basic air bubble entrainment processes are reviewed and the relevant dynamic similarities are discussed. Traditionally, physical studies are conducted using a Froude similitude which implies drastically smaller laboratory Reynolds numbers than in the corresponding prototype flows. Basic dimensional analyses are developed for both singular and interfacial aeration processes. The results are discussed in the light of systematic investigations and they show that the notion of scale effects is closely linked with the selection of relevant characteristic air–water flow properties. Recent studies of local air–water flow properties highlight that turbulence levels, entrained bubble sizes and interfacial areas are improperly scaled based upon a Froude similitude even in large-size models operating with the so defined Reynolds numbers up to 5 E+5. In laboratory models, the dimensionless turbulence levels, air–water interfacial areas and mass transfer rates are drastically underestimated.

An equation of motion for particles of finite Reynolds number and size

Abstract: In order to simulate the motion of bubbles, drops, and particles, it is often important to consider finite Reynolds number effects on drag, lift, torque, and history force. Herein, an equation of motion is developed for spherical particles with a no-slip surface based on theoretical analysis, experimental data, and surface-resolved simulations. The equation of motion is then extended to account for finite particle size. This extension is critical for particles which will have a size significantly larger than the grid cell size, particularly important for bubbles, and low-density particles. The extension to finite particle size is accomplished through spatial-averaging (both volume-based and surface-based) of the continuous flow properties. This averaging is consistent with the Faxen limit for solid spheres at small Reynolds numbers and added mass and fluid stress forces at inviscid limits. The finite Re p corrections are shown to have good agreements with experiments and resolved-surface simulations. The finite size corrections are generally fourth-order accurate and an order of magnitude more accurate than point-force expressions (which neglect quadratic and higher spatial gradients) for particles with size on the order of the gradient length-scales. However, further work is needed for more quantitative assessment of the truncation terms and the overall model robustness and accuracy in complex flows.

Dynamics of the buoyant plume off the Pearl River Estuary in summer

Abstract: Field measurements of salinity, wind and river discharge and numerical simulations of hydrodynamics from 1978 to 1984 are used to investigate the dynamics of the buoyant plume off the Pearl River Estuary (PRE), China during summer. The studies have shown that there are four major horizontal buoyant plume types in summer: Offshore Bulge Spreading (Type I), West Alongshore Spreading (Type II), East Offshore Spreading (Type III), and Symmetrical Alongshore Spreading (Type IV). River mouth conditions, winds and ambient coastal currents have inter-influences to the transport processes of the buoyant plume. It is found that all of the four types are surface-advected plumes by analysing the vertical characteristic of the plumes, and the monthly variations of the river discharge affect the plume size dominantly. The correlation coefficient between the PRE plume size and the river discharge reaches 0.85 during the high river discharge season. A wind strength index has been introduced to examine the wind effect. It is confirmed that winds play a significant role in forming the plume morphology. The alongshore wind stress and the coastal currents determine the alongshore plume spreading. The impact of the ambient currents such as Dongsha Current and South China Sea (SCS) Warm Current on the plume off the shelf has also assessed. The present study has demonstrated that both the river discharge and wind conditions affect the plume evolution.

Determination of the Manning roughness coefficient influenced by vegetation in the river Aa and Biebrza river

Abstract: The aim of this study was to investigate the variation of channel bed roughness in two rivers, as important parameter in hydraulic modelling especially with regard to flood control. The universities of Ghent (UG) and Antwerp (UA) are conducting scientific research in the river Aa in Belgium and the Biebrza river in Poland in order to better understand the phenomena involved and to come to a more accurate determination of the different parameters influencing flow. In this paper, the determination of the roughness coefficient ‘n’ from the Manning equation is used. This coefficient is not easy to determine and is varying constantly. It is influenced by the meandering character of the river, the bed material and the average grain size, the channel bed forms, the channel obstructions, the geometry changes between sections and the vegetation in the channel. Furthermore, due to these parameters, the roughness of the channel is not equally distributed over the channel, the banks and the floodplains. So, using literature data does not always lead to satisfactory results, due to the different situation in the field (Werner et al. J Hydrol 314:139–157, 2005). Therefore, measurements are necessary to determine the variation of the Manning coefficient. The Manning coefficient is a function of the discharge, but will also vary over the time due to the mentioned influences. In a multidisciplinary research project on the fundamental exchange processes in river ecosystems, hydraulic measurements were performed on a regular base in the river Aa. During these measurement campaigns, velocity and discharge measurements were carried out in multiple cross-sections. Once a month, the discharge and the water levels were measured at the upstream and the downstream end of the test stretch. On the river Biebrza, similar intensive measurement campaigns took place along a 6 km stretch in the upstream part of the river. An accurate determination of the Manning coefficient according a seasonal variation is an important tool in hydraulic modelling.

The impact of chemical lateral boundary conditions on CMAQ predictions of tropospheric ozone over the continental United States

Abstract: A sensitivity study is performed to examine the impact of lateral boundary conditions (LBCs) on the NOAA-EPA operational Air Quality Forecast Guidance over continental USA. We examined six LBCS: the fixed profile LBC, three global LBCs, and two ozonesonde LBCs for summer 2006. The simulated results from these six runs are compared to IONS ozonesonde and surface ozone measurements from August 1 to 5, 2006. The choice of LBCs can affect the ozone prediction throughout the domain, and mainly influence the predictions in upper altitude or near inflow boundaries, such as the US west coast and the northern border. Statistical results shows that the use of global model predictions for LBCs could improve the correlation coefficients of surface ozone prediction over the US west coast, but could also increase the ozone mean bias in most regions of the domain depending on global models. In this study, the use of the MOZART (Model for Ozone And Related chemical Tracers) prediction for CMAQ (Community Multiscale Air Quality) LBC shows a better surface ozone prediction than that with fixed LBC, especially over the US west coast. The LBCs derived from ozonesonde measurements yielded better O3 correlations in the upper troposphere.

Energy dissipation, flow resistance and gas-liquid interfacial area in skimming flows on moderate-slope stepped spillways

Abstract: With the re-evaluation and revision of a number of design floods, several embankment overtopping protection systems have been developed and a common technique is the construction of a stepped spillway on the downstream slope. For such moderate slope stepped channels, detailed air–water flow measurements were performed in a large facility with a focus on the rate of energy dissipation, flow resistance, air–water interfacial areas and re-aeration rates. Past and present experimental results showed a significant aeration of the flow. The median dimensionless residual head was about 3 × dc for the 21.8° sloping chute and smaller than that for flatter slopes (θ = 3.4° and 15.9°). The flow resistance results yielded an equivalent Darcy friction factor of about 0.25 implying a larger flow resistance for the 21.8° slope angle than for smaller slope angles. The re-aeration rate was deduced from the integration of the mass transfer equation using measured air–water interfacial areas and air–water flow velocities. The results suggested an increasing re-aeration rate with increasing rate of energy dissipation. The stepped invert contributed to intense turbulence production, free-surface aeration and large interfacial areas. The experimental data showed however some distinctive seesaw pattern in the longitudinal distribution of air–water flow properties with a wave length of about two step cavities. While these may be caused by the interactions between successive adjacent step cavities and their interference with the free-surface, the existence of such “instabilities” implies that the traditional concept of normal flow might not exist in skimming flows above moderate-slope stepped spillways.

Two-phase modeling of turbulence in dilute sediment-laden, open-channel flows

Abstract: In this paper, we focus on assessing the performance of diverse turbulence closures in the simulation of dilute sediment-laden, open-channel flows. To that end, we base our analysis on a framework developed in a companion paper of this special issue, which puts forward a standard sediment transport model (SSTM), a partial two-fluid model (PTFM) and a complete two-fluid model (CTFM), in three- and one-dimensional (3D and 1D) versions. First, we propose in this paper extensions of the transport equations for the Rey- nolds stresses, and of the equations of the K-ω model to two-phase flows, starting from the general two-fluid model. We consider the drag force to be the predominant force amongst all the interactions between the two phases (water and sediment). Second, under the framework of models formed by the SSTM, the PTFM and the CTFM, we discuss simulation results obtained by employing the Reynolds stress model (RSM), the algebraic stress model (ASM), and the K-e and the K-ω models (in their standard and extended versions), paired with each member of the framework. To assess the accuracy of the models, we compare numerical results with the experimental datasets of Vanoni, Trans ASCE 111:67-133, 1946; Coleman, Water Resour Res 22(10):1377-1384, 1986; Muste and Patel, J Hydraul Eng 123(9):742- 751, 1997; Nezu and Azuma, J Hydraul Eng 130:988-1001, 2004; Muste et al. Water Resour Res 41:W10402, 2005 . Third, we obtain from those comparisons the values of the Schmidt number that facilitate the agreement of model predictions with data. We conclude that the standard K-e model, the ASM and the K-ω models all provide satisfactory descriptions of flow variables and sediment concentrations in open-channel flows; further, we show that the more complicated RSM does not provide much improvement in dilute sediment transport as compared to those previous models, even when it is paired with the CTFM. We also show that the inclusion of model extensions in the turbulence closures does not improve the predictions for dilute mixtures either. We find that our values for the Schmidt number agree well with

Estimation of urban sensible heat flux using a dense wireless network of observations

Abstract: The determination of the sensible heat flux over urban terrain is challenging due to irregular surface geometry and surface types. To address this, in 2006–07, a major field campaign (LUCE) took place at the Ecole Polytechnique Federale de Lausanne campus, a moderately occupied urban site. A distributed network of 92 wireless weather stations was combined with routine atmospheric profiling, offering high temporal and spatial resolution meteorological measurements. The objective of this study is to estimate the sensible heat flux over the built environment under convective conditions. Calculations were based on Monin–Obukhov similarity for temperature in the surface layer. The results illustrate a good agreement between the sensible heat flux inferred from the thermal roughness length approach and independent calibrated measurements from a scintillometer located inside the urban canopy. It also shows that using only one well-selected station can provide a good estimate of the sensible heat flux over the campus for convective conditions. Overall, this study illustrates how an extensive network of meteorological measurements can be a useful tool to estimate the sensible heat flux in complex urban environments.

Mass exchange in a shallow channel flow with a series of groynes: LES study and comparison with laboratory and field experiments

Abstract: The exchange of dissolved matter between a straight open channel and a series of shallow embayments present at one of its sides is investigated using large eddy simulation (LES). The direct link between the mechanism of mass exchange and the dynamics of coherent structures is demonstrated. It is shown that for the geometrical configuration considered in the present study, the mass exchange process is very non-uniform over the depth in the vicinity of the channel–embayment interface. Most of the contaminant is ejected from the embayments close to the free surface. The amount of contaminant re-entrained into the embayments situated downstream of the one in which contaminant was introduced is quantified. The mass exchange coefficient predicted by LES does not vary significantly with the embayment rank and is in very good agreement with the one predicted by the model proposed by Weitbrecht et al. (J Hydraul Eng 134(2):173–183, 2008) based on the value of a dimensionless morphometric groyne-field parameter. Field experiments were purposely performed in a natural stream with embayments whose length over width ratios were close to the ratio in the LES study. The concentration fields predicted by LES are compared with video-records of colored dye used to visualize the mass exchange in the field experiment. It is shown that, for both LES and the field experiment, the dominant passage frequency of the eddies inside the interfacial mixing layer is well predicted by the analytical model of Sukhodolov and Sukhodolova (in: Cowen et al (eds) Hydraulic measurements & experimental methods. Proceedings of international conference, Lake Placid, USA, pp 172–177, 2007). The model is then used to scale the time in the LES animations and field video-records showing the temporal evolution of the concentration field. The results of the comparison indicate several similarities in the mixing process, despite the differences in the bathymetry and the large difference in the Reynolds number between LES and the field experiment. This proves the usefulness of performing detailed LES and laboratory studies in well-controlled environments to understand mass-exchange processes around river groyne fields.

Experimental investigation of bubbly flow and turbulence in hydraulic jumps

Abstract: Many environmental problems are linked to multiphase flows encompassing ecological issues, chemical processes and mixing or diffusion, with applications in different engineering fields. The transition from a supercritical flow to a subcritical motion constitutes a hydraulic jump. This flow regime is characterised by strong interactions between turbulence, free surface and air–water mixing. Although a hydraulic jump contributes to some dissipation of the flow kinetic energy, it is also associated with increases of turbulent shear stresses and the development of turbulent eddies with implications in terms of scour, erosion and sediment transport. Despite a number of experimental, theoretical and numerical studies, there is a lack of knowledge concerning the physical mechanisms involved in the diffusion and air–water mixing processes within hydraulic jumps, as well as on the interaction between the free-surface and turbulence. New experimental investigations were undertaken in hydraulic jumps with Froude numbers up to Fr = 8.3. Two-phase flow measurements were performed with phase-detection conductivity probes. Basic results related to the distributions of void fraction, bubble frequency and mean bubble chord length are presented. New developments are discussed for the interfacial bubble velocities and their fluctuations, characterizing the turbulence level and integral time scales of turbulence representing a “lifetime” of the longitudinal bubbly flow structures. The analyses show good agreement with previous studies in terms of the vertical profiles of void fraction, bubble frequency and mean bubble chord length. The dimensionless distributions of interfacial velocities compared favourably with wall-jet equations. Measurements showed high turbulence levels. Turbulence time scales were found to be dependent on the distance downstream of the toe as well as on the distance to the bottom showing the importance of the lower (channel bed) and upper (free surface) boundary conditions on the turbulence structure.

Hierarchical modeling of the dilute transport of suspended sediment in open channels

Abstract: We propose, discuss and validate a theoretical and numerical framework for sediment-laden, open-channel flows which is based on the two-fluid-model (TFM) equations of motion. The framework models involve mass and momentum equations for both phases (sediment and water) including the interactive forces of drag, lift, virtual mass and turbulent dispersion. The developed framework is composed by the complete two-fluid model (CTFM), a partial two-fluid model (PTFM), and a standard sediment-transport model (SSTM). Within the umbrella of the Reynolds-Averaged Navier-Stokes (RANS) equations, we apply K–e type closures (standard and extended) to account for the turbulence in the carrier phase (water). We present the results of numerical computations undertaken by integrating the differential equations over control volumes. We address several issues of the theoretical models, especially those related to coupling between the two phases, interaction forces, turbulence closure and turbulent diffusivities. We compare simulation results with various recent experimental datasets for mean flow variables of the carrier as well as, for the first time, mean flow of the disperse phase and turbulence statistics. We show that most models analyzed in this paper predict the velocity of the carrier phase and that of the disperse phase within 10% of error. We also show that the PTFM provides better predictions of the distribution of sediment in the wall-normal direction as opposed to the standard Rousean profile, and that the CTFM is by no means superior to the PTFM for dilute mixtures. We additionally report and discuss the values of the Schmidt number found to improve the agreement between predictions of the distribution of suspended sediment and the experimental data.

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

Abstract: 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 Karman 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 characteristics within sparse and dense canopies

Abstract: Boundary layer interactions with canopies control various environmental processes. In the case of dense and homogeneous canopies, the so-called mixing layer analogy is most generally used. When the canopy becomes sparser, a transition occurs between the mixing layer and the boundary layer perturbed by interactions between element wakes. This transition has still to be fully understood and characterized. The experimental work presented here deals with the effect of the canopy density on the flow turbulence and involves an artificial canopy placed in a fully developed turbulent boundary layer. One and two-component velocity measurements are performed, both within and above the canopy. The influence of the spacing between canopy elements is studied. Longitudinal velocity statistical moments and Reynolds stresses are calculated and compared to literature data. For spacings greater than the canopy height, evidences of this transition are found in the evolution of the skewness factor, shear length scale and mixing length.

Gas-phase and aerosol chemistry interactions in South Europe and the Mediterranean region

Abstract: The atmospheric chemical composition is affected by the interaction mechanisms among gases and particulate matter through a wide range of chemical reactions that can occur with the aid of particulate matter (e.g. particles act as reacting or absorbing surfaces) or be influenced by the presence of particulate matter in the atmosphere (photochemical reactions). Physical and chemical processes are also bonded in an interactive way that often leads to the influence of the radiation budget, cloud physics and the warming or cooling of the lower atmospheric levels. The Euro-Mediterranean region is a key-sensitive area due to the unique climatic and air quality characteristics associated with the regional climatic patterns, geomorphology (land and water contrast) and coexistence of pollutants from different origin. Focusing on this region, the gas-aerosol interactions are studied using state-of-the-art atmospheric and chemical transport modeling tools following the necessary development in the chemical transport model CAMx. Sensitivity and large-scale simulations have shown significant responses of the modeling system to the inclusion of natural species emissions, the direct shading effect of dust particles on photochemical processes and the formation of new types of aerosols through heterogeneous uptake of gases on dust particles. Including such interactions in the chemical transport model often led to the improvement of the model performance compared with available measurements in the region.

Using a sub-grid scale modeling approach to simulate the transport and fate of toxic air pollutants

Abstract: Several air toxics are emitted from mobile sources on roadways and these emissions account for a significant fraction of the health risks to the population. In addition, health effect studies are now becoming more comprehensive and some account for the spatial heterogeneities of air pollutant concentration fields (as is the case near roadways). Standard models can simulate either the near-source concentration fields or the urban background, but no model can handle both spatial scales in the vicinity of roadways in a coherent and scientifically correct manner. Here, we present a model that provides such an integrated treatment by combining a grid-based air quality model of the urban background with a plume-in-grid representation of roadway emissions. The model is applied to simulate near-roadway concentrations due to emissions from a busy interstate highway in New York City. Qualitative comparisons with typical measured concentration profiles show that the model captures the observed features of toxic air pollutant concentrations near roadways.

Probability law of concentration in plumes dispersing in an urban area

Abstract: The relationships between various normalized higher-order concentration moments in plumes dispersing in a built-up (urban) environment have been investigated using a large concentration data set obtained in a boundary-layer water channel. This data set consists of measurements of plume dispersion in a number of idealized obstacle arrays (e.g., cubical and non-cubical obstacles in aligned and staggered arrangements with uniform, random and alternating heights). A remarkably robust feature of all the concentration data was the observed collapse of the third- and fourth-order normalized concentration moments on the second-order normalized concentration moment. The data are shown to collapse to a series of universal curves (independent of the geometry of the obstacle array) and these curves were found to be identical to those observed previously for open-terrain plumes. The results imply that the probability law of concentration in a plume dispersing in either a built-up environment or open terrain has a universal form that can be specified by at most two independent parameters. The universal functions representing the relationships between the normalized concentration moments were found to be well modeled (approximated) using a two-parameter clipped-gamma probability law for the concentration. Finally, the clipped-gamma distribution was found to be in very good conformance with the measured probability distribution of concentration for plumes dispersing in a built-up environment.

Modelling the propagation of an internal solitary wave across double ridges and a shelf-slope

Abstract: Evolution of the internal solitary waves (ISWs) in the northern South China Sea (SCS) has recently attracted the attention of many oceanographers in Taiwan and the United States. These ISWs are believed to have been induced by a branch of the Kuroshio current over Luzon Strait, which propagates westward over two ridges in the Luzon Strait between Taiwan and the Philippines, and further onto the continental margin with a shelf-slope in the SCS. This paper presents some preliminary results for the evolution of a depression ISW across two triangular obstacles using numerical modelling and laboratory experiments. The experimental results confirm that the intervals and relative height between the two obstacles are important factors in the interaction of an ISW with the obstacles. However, in the case of the movement of an ISW of depression-type across the Luzon Strait, the effect of the two ridges on the characteristics of the ISW might be less significant than that from the shelf-slope, due to the variations in relative water depth. Results from numerical experiments also show that the amplitude of an ISW can be augmented once the wave commences its contact with a shelf-slope, where an internal hydraulic jump and wave breaking with vortex motion are evident in the laboratory experiments. Eventually, an ISW of depression-type could become an elevation-type at the edge of the continental shelf landwards beyond the turning point, where the upper layer is larger than the bottom layer in a stratified water column.

Dust emission calculations in open storage piles protected by means of barriers, CFD and experimental tests

Abstract: The importance of open mineral storage piles in bulk solids port terminals have increased considerably in recent years in Europe and USA (in Spain, great extensions of transoceanic ports are being made) to address the increasing demand of raw material importation due to local mining operation closure. These storage piles are affected by the climatologically processes and the dust emission to the atmosphere, causing repercussions to health and environmental, which is intimately related to the air velocity, according to USEPA studies, maximum ratio of dust emission emitted from CFD and software Ansys CFX 10.0 is determined. In addition, an emission variability study based on the pile distance of different height solid barriers is developed, reducing the emission to 66%. The studies have been contrasted by published investigations and industrial measurement.

Horizontal diffusion by submeso motions in the stable boundary layer

Abstract: Four networks of wind data are used to construct the first systematic estimates of the horizontal diffusivity from observations of submeso motions on scales often unresolved in numerical models. Currently, the horizontal diffusivity in numerical models is specified mainly for numerical reasons without observational support. The data analysis in this study emphasizes the stable boundary layer although results are briefly presented for the unstable boundary layer. The horizontal diffusivity is estimated from the horizontal gradient and the observed flux. Horizontal gradients of scalars are generally difficult to directly estimate from observations with sufficient accuracy for much of the data. As an alternative, simulated particles with conservative properties are introduced into the observed wind field in order to estimate the horizontal diffusivity for submeso motions. The sensitivity of the horizontal diffusivity to details of the method is examined. The horizontal diffusivity increases with the range of time and space scales that are included in the evaluation. The horizontal diffusivity is much larger with significant topography and may increase with wind speed, depending on the site location. The coarse station spacing or the small domain size is found to be a major limitation to the analysis.

Verification of the geographical origin of modeled air-mass trajectories by means of the isotope composition of rainwater during the SALLJEX experiment

Abstract: The SALLJEX experiment was held during the summer 2002–2003. It consisted of three-dimensional observation of the atmosphere to study the structure of the low level jet along the eastern slopes of the Andes. Daily precipitation water samples were collected at two stations (Resistencia and Salta) in northern Argentina and isotope content was analyzed. The isotope data were used in conjunction with air parcel trajectories obtained from a 3-D kinematic model (3D-MTC) developed by the University of Sao Paulo, Brazil. Values of deuterium excess were related with air masses of continental origin, whilst low values were associated with air masses with longer oceanic trajectories. Furthermore, although data are scarce, results show that oxygen-18 and deuterium excess in rainwater are related with the occurrence of the low level jet.

Determination of the optimum MM5 configuration for long term CMAQ simulations of aerosol bound pollutants in Europe

Abstract: Realistic meteorological fields are a prerequisite for the determination of pollutant concentrations and depositions by means of a chemistry transport model. Different configurations of the 5th generation NCAR/Penn State University mesoscale meteorological model MM5 were tested to determine the optimum set up for long term hindcasts that cover several months up to years. Four dimensional data assimilation (FDDA) significantly enhances the spatio temporal representation of temperature, humidity and wind. Best agreement with radiosonde observations could be achieved when temperature, humidity and wind were grid nudged every 6 h. The quality of the resulting meteorological fields showed no significant systematic temporal or spatial variation over Europe in a model run of the year 2000. It was found that the hydrological cycle was not correctly reproduced by the model when no nudging was applied. The relevant model run showed too high relative humidity and too high rainfall when compared to observations. This led to considerably lower aerosol concentrations close to ground and a shift in the deposition patterns of particle bound pollutants like the carcinogenic benzo(a)pyrene (B(a)P).

Evaluating the dimensionality and significance of “active periods” in turbulent environmental flows defined using Lipshitz/Hölder regularity

Abstract: Active periods within perturbed boundary-layer flows are considered in terms of the local roughness of measured velocity time series and defined in terms of Holder/Lipshitz exponents. Such events are associated with the passage of energetic, coherent flow structures and are responsible for exerting high turbulent stresses because of the rapid changes in velocity that occur at such times. A method is proposed for assessing the effective dimensionality of such active periods, as well as their significance to the flow field, for a particular choice of flow metric. The method is applied to the turbulent flow through a confluence flow geometry, with velocity samples acquired close to the bed of the channel in a zone of complex mixing. The dimensionality of the active periods is consistent with the observed patterns of sediment entrainment from the bed, with the significance of the active periods decaying away from the erosional zone.

Progress in the observation and modeling of turbulent multi-phase flows

Abstract: Multi-phase flows of environmental import involve the interactions of solid, gas, and/or liquid phases. In the last decades, important advances have taken place in the observation and modeling of multi-phase flows. New accurate instruments to measure the air concentration in air-water flows and to analyse the motion of solid particles in water and gas; new theories to address the interactions of the diverse turbulence scales of the flow; and more powerful computers have all allowed researchers to study multi-phase flows from previously unexplored points of view. In Nature, multi-phase flows are mostly turbulent and, therefore, they are extremely complicated, with a broad range of relevant length and time scales [1]. For example, the time scales range from less than 1ms for the turbulence dissipation in a small stream to about 24 h and 50min for a tidal cycle in coastal zones (Fig. 1), and to more than 50 years for the currents controlling the balances between oxygen and carbon dioxide. Several recent papers have remarked that the difficulties in understanding multi-phase flows stem from two basic facts: the phases do not distribute uniformly, and the small-scale interactions may have profound effects on large-scale behaviours. These difficulties become important for the prediction of the behaviour of flows in many environmentalproblems and applications. The aim of this Special Issue is to present a group of papers that: (a) summarize thestate-of-the-art in the knowledge about multi-phase flows in environmental applications; (b)report recent results of research (of experimental, numerical and/or theoretical nature) on environmental multi-phase flows; and (c) suggest novel pathways of analysis in the area.

Secondary flow within a river bed and contaminant transport

Abstract: We have developed a numerical method to simulate the transport of non-sorbing contaminants within the sediment layer of a stream and the leaching of these contaminants in the steam. Typical stream bottom surfaces are uneven with triangularly shaped undulation forms. The flow of the water above such triangular surfaces causes external pressure changes that result in a “pumping effect” and a secondary flow within the sediment. The latter causes a significant contaminant advection within the sediment layer. The flow field in the porous sediment layer is obtained by solving numerically Darcy’s equations. The unsteady mass transfer equation is solved by using a finite-difference method with an up-wind scheme. The effects of parameters, such as channel slope, hydraulic head and dispersion, are studied by quantitatively comparing the numerical results of the total mass flow rate from the contaminant source, the concentration front propagation, and the contaminant mass flow rate into the water column. The “pumping effect,” increases the flow in the vertical direction and, thus, enhances the vertical advective mass transport of the contaminant. This bedform-shape induced flow is largely responsible for the mass transfer of contaminants into the water column. The numerical results also show that the mechanical dispersion inside the sediment bed will significantly increase the contaminant mass flow rate from the source.

Performance evaluation of NOAA-EPA developmental aerosol forecasts

Abstract: To aid air quality model development and assess air quality forecasts, the Meteorological Development Laboratory (MDL) provided categorical verification metrics for developmental aerosol predictions. The National Air Quality Forecasting Capability (NAQFC) generated 48 h (of) gridded hourly developmental predictions for the lower 48 states (CONUS) domain in 12 km horizontal spacing. The NAQFC uses the North American Mesoscale (NAM) model with EPA’s Community Multiscale Air Quality (CMAQ) model to produce predictions of ground level aerosol concentrations. We used bilinear interpolation to calculate predicted daily maximum values at the location of the observation sites. We compared these interpolated predicted values to the observed daily maximum to produce 2 × 2 contingency tables, with a threshold of 40 μg/m3 during the months of March–August, 2007. The model showed some degree of skill in predicting aerosol exceedances. These results are preliminary as the NAQFC model for aerosol prediction is in the developmental stage. A more comprehensive performance evaluation will be accomplished in 2008, when more data become available. Our verification metrics included categorical analyses for Fraction Correct (FC) or percent correct (FC × 100), Threat Score (TS) or Critical Success Index (CSI), Probability of Detection (POD), and the False Alarm Rate (FAR), Mean Absolute Error (MAE) and mean algebraic error or bias, where bias is forecast minus observation. Graphic products included weekly statistics for the CONUS displayed in the form of bar charts, scatterplots, and graphs. In addition, we split the CONUS into six geographic regions and provided regional statistics on a monthly basis. MDL produced spatial maps of daily 1-h maximum predicted aerosol values overlaid with the corresponding point observations. MDL also provided spatial maps of the daily maximum of the 24-h running average. We derived the 24-h running average from the 1-h average predicted aerosol values and observations.

Impact of consistent boundary layer mixing approaches between NAM and CMAQ

Abstract: 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 O3 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 O3 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 O3 forecast related processes are warranted.

Evaluation of wind environment around a residential complex using a PIV velocity field measurement technique

Abstract: Wind-tunnel simulations were employed to evaluate the wind environment around a tested residential area located near industrial complexes. The scaled-down geomorphological model of the test area was placed in the test section of a boundary layer wind tunnel. Particle image velocimetry (PIV) measurements were made in five vertical planes and one horizontal plane around the test area for two prevailing wind directions. The results showed that the wind speed decreased in the near surface layer and the velocity fluctuations increased in the upper region due to the presence of hills and high-rise buildings around the test area. Regions of flow separation and low-speed flow were found inside the test area for both the wind directions. The result suggests that the high-rise buildings should be well arranged with respect to the main wind directions to increase the natural ventilation inside the residential complex at the initial design stage.

Sensitivity of air quality model prediction to parameterization of vertical eddy diffusivity

Abstract: This paper investigates the effects of vertical eddy diffusivities derived from the 3 different planetary boundary layer (PBL) schemes on predictions of chemical components in the troposphere of East Asia. Three PBL schemes were incorporated into a regional air quality model (RAQM) to represent vertical mixing process and sensitivity simulations were conducted with the three schemes while other options are identical. At altitudes <2km, all schemes exhibit similar skill for predicting SO2 and O3, but more difference in the predicted NOx concentration. The Gayno–Seaman scheme produces the smallest vertical eddy diffusivity (Kz) among all schemes, leading to higher SO2 and NOx concentrations near the surface than that from the other 2 schemes. However, the effect of vertical mixing on O3 concentration is complex and varies spatially due to chemistry. The Gayno–Seaman scheme predicts lower O3 concentrations than the other two schemes in the parts of northern China (generally VOC-limited) and higher ones in most parts of southern China (NOx-limited). The Byun and Dennis scheme produces the largest mixing depth in the daytime, which bring more NOx into upper levels, and the mixing depth predicted by the Gayno–Seaman scheme is the smallest, leading to higher NOx and lower O3 concentrations near the surface over intensive emission regions.

A numerical (LES) investigation of a shallow-water, mid-latitude, tidally-driven boundary layer

Abstract: We investigate turbulent mixing in a tidally driven, mid-latitude, shallow-water basin The study is carried out numerically at a laboratory-scale, using large-eddy simulation We compared the results of the simulation with those of a correspondent purely oscillatory flow (Stokes boundary layer) The effect of rotation on the flow dynamics is twofold First, rotation gives rise to a mean spanwise flow that concurs to redistribute the turbulent energy among the Reynolds stresses, in particular between the horizontal directions, thus increasing the mixing across the water column and thickening the layer where developed turbulence is observable Second, the presence of the horizontal component of the background vorticity (latitude effect) breaks the symmetry between the two semi-cycles of the oscillation, since turbulence results suppressed/enhanced during the first/second semi-cycle These two effects significantly modify the turbulent characteristics with respect to the purely oscillating flow, although the mechanisms that generates turbulence present similar features The qualitative agreement between our results and some measurements carried out in two sites with characteristics similar to the case analyzed suggests that the outcomes here provided may be of general use for the analysis of mid-latitude, neutrally stratified, shallow-water basins mainly driven by semi-diurnal tidal currents

Energy conservation and second-order statistics in stably stratified turbulent boundary layers

Abstract: We address the dynamical and statistical description of stably stratified turbulent boundary layers with the important example of the atmospheric boundary layer with a stable temperature stratification in mind. Traditional approaches to this problem, based on the profiles of mean quantities, velocity second-order correlations, and dimensional estimates of the turbulent thermal flux run into a well-known difficulty, predicting the suppression of turbulence at a small critical value of the Richardson number, in contradiction with observations. Phenomenological attempts to overcome this problem suffer from various theoretical inconsistencies. Here we present a closure approach taking into full account all the second-order statistics, which allows us to respect the conservation of total mechanical energy. The analysis culminates in an analytic solution of the profiles of all mean quantities and all second-order correlations removing the unphysical predictions of previous theories. We propose that the approach taken here is sufficient to describe the lower parts of the atmospheric boundary layer, as long as the Richardson number does not exceed an order of unity. For much higher Richardson numbers the physics may change qualitatively, requiring careful consideration of the potential Kelvin-Helmoholtz waves and their interaction with the vortical turbulence.