A Multidisciplinary Study of Spatial and Temporal Scales Containing Information in Turbulent Chemical Plume Tracking

This report describes the results of a multidisciplinary study of turbulent chemical plume tracking of blue crabs and autonomous agents. The study consists of a coordinated investigation of animal behavior, fluid mechanics, strategy simulations, and chemical sensing. The objective is to provide a comprehensive understanding of chemical plume tracking in a single biological system and to prescribe strategies that are effective for autonomous agents. The consensus of the study is that spatial variation in the plume, measured by sampling at multiple locations simultaneously, yields information that is useful for plume tracking. Behavioral investigations reveal that blue crabs demonstrate the ability to detect the chemical plume and use lateral movements to avoid losing contact with the odor. Blue crabs move rapidly towards the source, strongly suggesting that temporal comparisons of odor properties are not employed during navigation. Analysis of the concentration fields reveals that a spatial correlation between spanwise-separated sensors indicates the relative direction of the plume centerline over short time periods provided the sensor spacing is scaled appropriately relative to the plume. Similarly, simulations of tracking strategies reveal an optimal separation for the sensors at a distance roughly equal to the plume width; both smaller and larger sensor spans degrade tracking performance. The simulations further reveal an optimal sensor size above which the fine details of the concentration distribution are obscured and below which there is insufficient contact with the odor to enable effective navigation. Finally, analysis of the chemical signal shows that the frequency dependent correlation function between two (or more) sensors indicates the relative position of the source.     

Measurement of Odor-Plume Structure in a Wind Tunnel Using a Photoionization Detector and a Tracer Gas

The patterns of stimulus available to moths flying along pheromone plumes in a 3-m-long wind tunnel were characterized using a high frequency photoionization detector in conjunction with an inert tracer gas. Four contrasting flow regimes and source conditions were produced: odor released in pulses from a vertical and horizontal array of four sources, odor released continuously from a point source, and odor released continuously from a point source into an oscillatory wake. Although the four flow regimes produced plumes of intermittent and fluctuating concentration, there were considerable differences in the structure of the signal presented to the sensor. Pulses of tracer gas released at 10 Hz retained most of their longitudinal and lateral separation. The plume growing in the disturbed flow (`oscillatory'), was broader in its lateral extent than the plume growing in an undisturbed flow (`continuous'), and the concentrations in the former were lower at each downstream position. The signal recorded in the disturbed flow had higher intermittency, but the ratio between the peak concentration and the signal mean was lower than in the continuous plume. Time scales were typically longer in the tunnel than in a field setting, but length scales and the main features of intermittency and fluctuation were similar. Moths flying along plumes of pheromone in this and similar wind tunnels typically slow their velocity and narrow the lateral excursions of their flight track as they approach a pheromone source. Which features of the plumes measured in this study account for these behavioral reactions remains to be determined.     

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

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.     

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.     

The Information Content of a Scalar Plume – A Plume Tracing Perspective

The ability of many animals and insects to track a plume to its source is a particularly impressive feat when the fluid dynamics is considered. Inspired by this observation this research seeks to identify the information in a passive scalar plume suitable for developing robust and efficient plume tracing algorithms. The subject of this study is a scalar plume emanating from a point source in a turbulent boundary layer which has been modeled in a laboratory facility built specifically for this purpose. A coupled PIV-LIF technique is used to measure the velocity and scalar field in a time resolved fashion. This data set is analyzed and the convergence rates of five single-point statistics, suitable as kernels of plume tracing algorithms, are investigated. The experimental data shows that the scalar fluctuations over long downstream distances from the source are characterized by filamentary structures that lead to relatively slow convergence rates for any statistic that is based on mean concentrations. The scalar intermittency, however, converges rapidly toward its true value, in fact converging to a testable hypothesis for source location direction faster than the time scale of the larger scale plume meander.     

Embedded large eddy simulation approach for pollutant dispersion around a model building in atmospheric boundary layer

In the present article, the potential of embedded large eddy simulation (ELES) approach to reliably predict pollutant dispersion around a model building in atmospheric boundary layer is assessed. The performance of ELES in comparison with large eddy simulation (LES) is evaluated in several ways. These include a number of qualitative and quantitative comparisons of time-averaged and instantaneous results with wind tunnel measurements supplemented by statistical data analyses using scatter plots and standard evaluation metrics. Results obtained by both LES and ELES approaches show very good agreement with the experiment. However, addition of turbulence to mean flow at Reynolds averaged Navier–Stokes (RANS)–LES interface in ELES approach not only increases the turbulence intensity, it also results in larger values of turbulent kinetic energy (TKE) as well as a shorter reattachment length in the wake region. Accordingly, higher levels of TKE predicted by ELES increase the local intensity of concentration leading to shorter plume shapes as compared with LES. In general, ELES shows better agreement with experiment on the surfaces of model building and also in the downstream wake region. In terms of computational costs, the CPU time required to obtain statistical values in ELES is about 49 % lower than that of LES and the number of iterations per time step is also reduced by 55 % as compared with LES.     

Numerical Simulation of Interaction of the Heavy Gas Cloud with the Atmospheric Surface Layer

The numerical time-dependent three-dimensional model [Kovalets, I.V. and Maderich, V.S.: 2001, Int. J. Fluid Mech. Res. 30, 410–429] of the heavy gas dispersion in the atmospheric boundary layer has been improved by parameterizing momentum and heat fluxes on the surface of Earth using Monin–Obukhov similarity theory. Three parameterizations of heat exchange with the surface of Earth were considered: (A) formula of Yaglom A.M. and Kader B.A. [1974, J. Fluid Mech. 62, 601–623] for forced convection, (B) interpolation formula for mixed convection and (C) similarity relationship for mixed convection [Kader, B.A. and Yaglom, A.M.: 1990, J. Fluid Mech. 212, 637–662]. Two case studies were considered. In the first study based on experiment of Zhu et al., J. Hazard Mater 62, 161–186], the interaction of an isothermal heavy gas plume with an atmospheric surface layer was simulated. It was found that stable stratification in the cloud essentially suppresses the turbulence in the plume, reducing the turbulent momentum flux by a factor of down to 1/5 in comparison with the undisturbed value. This reduction essentially influences velocities in the atmospheric boundary layer above the cloud, increasing the mean velocity by a factor of up to 1.3 in comparison with the undisturbed value. A simulation of cold heavy gas dispersion was carried out in the second case based on field experiment BURRO 8. It was shown that both forced and free convections under moderate wind speeds significantly influence the plume. The relative rms and bias errors of prediction the plume’s height were σ_H ≈ 30% and ɛ_H = − 10%, respectively, for parameterization B, while for A and C the errors were σ_H ≈ 80% and ɛ_H ≈ − 65%. It is therefore advised to use the simple parameterization B in dense gas dispersion models.     

Desert ants adjust their approach to a foraging site according to experience

In the habitat of desert ants, Cataglyphis fortis, a constant wind is usually blowing during the daytime. When visiting a familiar food source, the ants steer some distance downwind of the feeder, rather than attempting a direct approach that might miss small food sources, in particular. In the downwind area, the ants pick up the odor plume emanating from the food and follow it upwind to the prey. This strategy saves considerable walking distance and time. The additional path necessitated by the downwind strategy is only about 0.75 to 2 m, depending on nest–feeder distance, while missing the food on the upwind side results in much longer search trajectories. During the initial three to five visits to a feeding site, downwind distance and length of the approach path are shortened notably, and the approach trajectory is straightened. Desert ants further exhibit considerable short-term flexibility in their approach. Experienced individuals are evidently able to decide upon leaving the nest which direction to choose toward the feeder, depending on current wind direction (that fluctuates slightly during the day). Notable changes in wind direction occur primarily overnight. For larger nest–feeder distances, the animals adjust their approach en route to the altered wind direction during their first foraging trip in the morning.     

Evaluation of the QUIC-URB fast response urban wind model for a cubical building array and wide building street canyon

This paper describes the QUIC-URB fast response urban wind modeling tool and evaluates it against wind tunnel data for a 7 × 11 cubical building array and wide building street canyon. QUIC-URB is based on the Röckle diagnostic wind modeling strategy that rapidly produces spatially resolved wind fields in urban areas and can be used to drive urban dispersion models. Röckle-type models do not solve transport equations for momentum or energy; rather, they rely heavily on empirical parameterizations and mass conservation. In the model-experiment comparisons, we test two empirical building flow parameterizations within the QUIC-URB model: our implementation of the standard Röckle (SR) algorithms and a set of modified Röckle (MR) algorithms. The MR model attempts to build on the strengths of the SR model and introduces additional physically based, but simple parameterizations that significantly improve the results in most regions of the flow for both test cases. The MR model produces vortices in front of buildings, on rooftops and within street canyons that have velocities that compare much more favorably to the experimental results. We expect that these improvements in the wind field will result in improved dispersion calculations in built environments.     

A Computational Fluid Dynamics Approach for Urban Area Transport and Dispersion Modeling

A simulation tool has been developed to model the wind fields, turbulence fields, and the dispersion of Chemical, Biological, Radiological and Nuclear (CBRN) substances in urban areas on the building to city blocks scale. A Computational Fluid Dynamics (CFD) approach has been taken that naturally accounts for critical flow and dispersion processes in urban areas, such as channeling, lofting, vertical mixing and turbulence, by solving the steady-state, Reynolds-Averaged Navier–Stokes (RANS) equations. Rapid generation of high quality cityscape volume meshes is attained by a unique voxel-based model generator that directly interfaces with common Geographic Information Systems (GIS) file formats. The flow and turbulence fields are obtained by solving the steady-state RANS equations using a collocated, pressure-based approach formulated for unstructured and polyhedral mesh elements. Turbulence modeling is based upon the Renormalization Group variant of the k–ε model (k–ε RNG). Neutrally buoyant simulations are made by prescribing velocity boundary condition profiles found by a power–law relationship, while turbulence quantities boundary conditions are defined by a prescribed mixing length in conjunction with the assumption of turbulence equilibrium. Dispersion fields are computed by solving an unsteady transport equation of a dilute gas, formulated in a Eulerian framework, using the velocity and turbulence fields found from the steady-state RANS solution. In this paper the model is explained and detailed comparisons of predicted to experimentally obtained velocity, turbulence and dispersion fields are made to neutrally stable wind tunnel and hydraulic flume experiments.     

Aphid population response to agricultural landscape change: A spatially explicit,individual-based model

A spatially explicit individual-based simulation model has been developed to represent aphid population dynamics in agricultural landscapes. The application of the model to Rhopalosiphum padi (L.) population dynamics is detailed, including an outline of the construction of the model, its parameterisation and validation. Over time, the aphids interact with the landscape and with one another. The landscape is modified by varying a simple pesticide regime, and the multi-scale spatial and temporal implications for a population of aphids is analysed. The results show that a spatial modelling approach that considers the effects on the individual of landscape properties and factors such as wind speed and wind direction provides novel insight into aphid population dynamics both spatially and temporally. This forms the basis for the development of further simulation models that can be used to analyse how changes in landscape structure impact upon important species distributions and population dynamics.     

Effects of the Resolution and Kinematics of Olfactory Appendages on the Interception of Chemical Signals in a Turbulent Odor Plume

A variety of animals use olfactory appendages bearing arrays of chemosensory neurons to detect chemical signatures in the water or air around them. This study investigates how particular aspects of the design and behavior of such olfactory appendages on benthic aquatic animals affect the patterns of intercepted chemical signals in a turbulent odor plume. We use virtual olfactory `sensors' and `antennules' (arrays of sensors on olfactory appendages) to interrogate the concentration field from an experimental dataset of a scalar plume developing in a turbulent boundary layer. The aspects of the sensors that we vary are: (1) The spatial and temporal scales over which chemical signals arriving at the receptors of a sensor are averaged (e.g., by subsequent neural processing), and (2) the shape and orientation of a sensor with respect to ambient water flow. Our results indicate that changes in the spatial and temporal resolution of a sensor can dramatically alter its interception of the intermittency and variability of the scalar field in a plume. By comparing stationary antennules with those sweeping through the flow (as during antennule flicking by the spiny lobster, Panulirus argus), we show that flicking alters the frequency content of the scalar signal, and increases the likelihood that the antennule encounters peak events. Flicking also enables a long, slender (i.e., one-dimensional) antennule to intercept two-dimensional scalar patterns.     

Concentration fluctuations in a downtown urban area. Part II: analysis of Joint Urban 2003 wind-tunnel measurements

An analysis of concentration time series measured in a boundary-layer wind tunnel at the University of Hamburg is presented. The measurements were conducted with a detailed aerodynamic model of the Oklahoma City (OKC) central business district (CBD) at the scale of 1:300 and were part of the Joint Urban 2003 (JU2003) project. Concentration statistics, as well as concentration probability density (PDF) and exceedance probability (EDF) functions were computed for street- and roof-level sites for three different wind directions. Taking into account the different length scales and wind speeds in the wind-tunnel (WT) and full-scale experiments, dimensionless concentrations and a dimensionless time scale are computed for the comparison with data from the JU2003 full-scale tracer experiments, conducted in OKC in 2003. Using such dimensionless time, the WT time series cover a ~20 times longer time span than the JU2003 full-scale time series, which are analysed in detail in an accompanying, first part of this paper. The WT time series are thus divided into 20 consecutive blocks of equal length and the statistical significance of parameters based on relatively short records is assessed by studying the variability of the concentration statistics and probability functions for the different blocks. In particular at sites closer to the plume edge, the results for the individual blocks vary significantly and at such sites statistics from short records are not very representative. While the location of three sampling sites in the WT closely matched the sites during the full-scale experiments, the prevailing wind directions during the JU2003 releases were not exactly matched. The comparison between full-scale and WT concentration parameters should thus primarily be interpreted in a qualitative rather than direct quantitative sense. Given the differences in mean wind directions and concerns about the representativeness of full-scale concentration statistics, the WT and full-scale results compared well. The 98 percentile concentrations for almost all full-scale releases analyzed are within the scatter of the percentiles observed in the block analysis of the WT time series. Furthermore, the concentration percentiles appear linearly correlated with the fluctuation intensities and the linear relationships determined in the wind tunnel agree well with full-scale results.     

Statistical properties of concentration fluctuations in two merging plumes

The statistics of the fluctuating concentration field within a plume is important in the analysis of atmospheric dispersion of toxic, inflammable and odorous gases. Previous work has tended to focus on concentration fluctuations in single plumes released in the surface layer or at ground level and there is a general lack of information about the mixing of two adjacent plumes and how the statistical properties of the concentration fluctuations are modified in these circumstances. In this work, data from wind tunnel experiments are used to analyse the variance, skewness, kurtosis, intermittency, probability density function and power spectrum of the concentration field during the mixing of two identical plumes and results are compared with those obtained for an equivalent single plume. The normalised variance, skewness and kurtosis on the centre-lines of the combined plume increase with distance downwind of the stack and, in the two-source configuration, takes lower values than those found in the single plumes. The results reflect the merging process at short range, which is least protracted for cases in which the sources are in-line or up to 30 \(^{\circ }\) off-line. At angles of 45 \(^{\circ }\) and more, the plumes are effectively side-by-side during the merging process and the interaction between the vortex pairs in each plume is strong. Vertical asymmetry is observed between the upper and the lower parts of the plumes, with the upper part having greater intermittency (i.e. the probability that no plume material is present) and a more pronounced tail to the concentration probability distribution. This asymmetry tends to diminish at greater distances from the source but occurs in both buoyant and neutral plumes and is believed to be associated with the ‘bending-over’ of the emission in the cross-flow and the vortex pair that this generates. The results allowed us to identify three phases in plume development. The first, very near the stack, is dominated by turbulence generated within the plume and characterised by concentration spectra with distinct peaks corresponding to scales comparable with those of the counter-rotating vortex pair. A second phase follows at somewhat greater distances downwind, in which there are significant contributions to the concentration fluctuations from both the turbulence internal to the plume and the external turbulence. The third phase is one in which the concentration fluctuations appear to be controlled by the external turbulence present in the ambient flow.     

Plume rise and spread in buoyant releases from elevated sources in the lower atmosphere

This study focuses on the influence of emission conditions—velocity and temperature—on the dynamics of a buoyant gas release in the atmosphere. The investigations are performed by means of wind tunnel experiments and numerical simulations. The aim is to evaluate the reliability of a Lagrangian code to simulate the dispersion of a plume produced by pollutant emissions influenced by thermal and inertial phenomena. This numerical code implements the coupling between a Lagrangian stochastic model and an integral plume rise model being able to estimate the centroid trajectory. We verified the accuracy of the plume rise model and we investigated the ability of two Lagrangian models to evaluate the plume spread by means of comparisons between experiments and numerical solutions. A quantitative study of the performances of the models through some suitable statistical indices is presented and critically discussed. This analysis shows that an additional spread has to be introduced in the Lagrangian trajectory equation in order to account the dynamical and thermal effects induced by the source conditions.     

Sensitivity of atmospheric dispersion simulations by HYSPLIT to the meteorological predictions from a meso-scale model

Mesoscale transport and dispersion of air pollutants from a few major point sources in the Mississippi Gulf coastal region is calculated using a coupled modeling system consisting of the atmospheric dynamical model WRF and the lagrangian particle model HYSPLIT. The sensitivity of the dispersion model results to the meteorological fields is studied by conducting an ensemble of simulations using the WRF model for the same dispersion case. Several parameterization schemes for the physical processes of boundary layer turbulence and land surface temperature/moisture prediction in WRF are used in various combinations to produce different meteorological members which are then used for dispersion simulation. The uncertainty in the simulated concentration probabilities to the meteorological model configurations and the ensemble mean are presented. The parameters used for determining the uncertainties include the wind fields, temperature, area of concentration and the levels of concentration. The results indicate that dispersion model results are influenced by the choices made in respect of the planetary boundary layer and land surface schemes in the mesoscale model to produce the meteorological forecast thereby leading to certain amount of uncertainty in the resultant concentrations. Results show that the specific choices made about the atmospheric model configuration can significantly after the simulated concentrations.     

Air-mediated pollen flow from genetically modified to conventional crops.

Tools for estimating pollen dispersal and the resulting gene flow are necessary to assess the risk of gene flow from genetically modified (GM) to conventional fields, and to quantify the effectiveness of measures that may prevent such gene flow. A mechanistic simulation model is presented and used to simulate pollen dispersal by wind in different agricultural scenarios over realistic pollination periods. The relative importance of landscape-related variables such as isolation distance, topography, spatial configuration of the fields, GM field size and barrier, and environmental variation are examined in order to find ways to minimize gene flow and to detect possible risk factors. The simulations demonstrated a large variation in pollen dispersal and in the predicted amount of contamination between different pollination periods. This was largely due to variation in vertical wind. As this variation in wind conditions is difficult to control through management measures, it should be carefully considered when estimating the risk of gene flow from GM crops. On average, the predicted level of gene flow decreased with increasing isolation distance and with increasing depth of the conventional field, and increased with increasing GM field size. Therefore, at a national scale and over the long term these landscape properties should be accounted for when setting regulations for controlling gene flow. However, at the level of an individual field the level of gene flow may be dominated by uncontrollable variation. Due to the sensitivity of pollen dispersal to the wind, we conclude that gene flow cannot be summarized only by the mean contamination; information about the frequency of extreme events should also be considered. The modeling approach described in this paper offers a way to predict and compare pollen dispersal and gene flow in varying environmental conditions, and to assess the effectiveness of different management measures.     

轨迹烟流模式计算方法

集合了高斯烟流模式和烟团模式的优点,提出一个新的模型———轨迹烟流模式,既保留了烟流模式计算简单,计算结果连续的特点,又可以根据当时风场,确定不同的烟流轨迹,模拟烟团的弯曲、折回和重叠运动,分段进行计算,提高了计算的准确性扩大了模式应用范围。     

A simplified computational analysis of turbulent plumes and jets

The application of computational fluid dynamics (CFD), particularly Large Eddy Simulation, for the modelling of buoyant turbulent plumes, has been demonstrated to be very accurate, but computationally expensive. Here a more basic, and therefore more generally practicable, approach is presented. Commercial CFD software is used to model such plumes using Reynolds-Averaged Navier-Stokes (RANS) turbulence models. A careful comparison is made between the numerical predictions and well-established results regarding the bulk properties of plumes. During this process, we are able to observe the well-known approximate Gaussian nature of the plume and achieve quantitative agreement with empirical plume spread coefficients. The use of numerical modelling allows for the investigation of the flow field and turbulence in those regions of the plume of most interest—the plume edge and near source regions. A comprehensive sensitivity study is conducted to identify the limits of applicability of this modelling approach. It is shown that the standard modelling approach of Morton, Taylor and Turner, which introduced the well-known entrainment assumption, pertains in a region well above the source region. At the plume edge, the levels of turbulence are contrasted with the value of the entrainment parameter. Finally, the effects of forcing the plumes with additional momentum at the source are considered, including the case of a pure jet. We show how these forced plumes eventually lose their momentum excess and tend to the behaviour of a pure, buoyant plume.     

重气连续泄漏扩散的盐水模拟实验

根据相似理论,分析了用盐水模拟方法研究意外泄漏的重气在大气中扩散过程的可行性,通过假设重气为不可压缩气体以及在均匀温度场中扩散等条件,推导出模拟实验的准则数.采用缩比模型的盐水模拟实验,对重气在大气中的扩散速度及浓度变化进行了分析,证实了重气扩散过程中的重力沉降、密度分层以及近源区分叉等现象.实验表明,在满足主导准则数相等的条件下,盐水在清水中的扩散可以较好地再现重气在大气中的扩散特征.