Showing papers in "Environmental Fluid Mechanics in 2007"

Weak-wind mesoscale meandering in the nocturnal boundary layer

Abstract: This study examines the strength and statistical behavior of mesoscale motions on time scales up to 1 h using eight data sets over different surface types. The mesoscale motions include internal gravity waves, microfront-like structures, horizontal modes, and a complex variety of other signatures, perhaps resulting from superposition of different modes. With weak large-scale flow, the mesoscale motions lead to meandering of the wind direction, as found in previous studies. However, the meandering often takes the form of sudden wind shifts rather than oscillation of wind direction.

Numerical simulation of scour around pipelines using an Euler-Euler coupled two-phase model

Abstract: The scour around a long fixed pipeline placed just above a non-cohesive sandy bed is numerically simulated using an Eulerian two-phase model that implements Euler–Euler coupled governing equations for fluid and solid phases and a modified k−ɛ turbulence closure for the fluid phase, the modeling system being a part of the CFD software package FLUENT. Both flow–particle and particle–particle interactions are considered in the model. During the simulations, the interface between sand and water is specified using a threshold volume fraction of sand, and the evolution of the bedforms is studied in detail. The predictions of bedform evolution are in good agreement with previous laboratory measurements. Investigations into the mechanisms of scour reveal that three sediment transport modes (bed-load, suspended-load and laminated-load) are associated with the scour development. While some previously proposed scour development formulae for cylindrical objects are in good general agreement with the simulations, scour prediction based on a commonly used operational mine-burial model (DRAMBUIE) shows disparities with present simulations.

Experimental analysis of Froude number effect on air entrainment in the hydraulic jump

Abstract: A hydraulic jump is characterized by strong energy dissipation and mixing, large-scale turbulence, air entrainment, waves, and spray. Despite recent pertinent studies, the interaction between air bubbles diffusion and momentum transfer is not completely understood. The objective of this paper is to present experimental results from new measurements performed in a rectangular horizontal flume with partially developed inflow conditions. The vertical distributions of the void fraction and the air bubbles count rate were recorded for inflow Froude number Fr1 in the range from 5.2 to 14.3. Rapid detrainment process was observed near the jump toe, whereas the structure of the air diffusion layer was clearly observed over longer distances. These new data were compared with previous data generally collected at lower Froude numbers. The comparison demonstrated that, at a fixed distance from the jump toe, the maximum void fraction Cmax increases with the increasing Fr1. The vertical locations of the maximum void fraction and bubble count rate were consistent with previous studies. Finally, an empirical correlation between the upper boundary of the air diffusion layer and the distance from the impingement point was derived.

Advanced Post-Processing and Correlation Analyses in High-Velocity Air-Water Flows

Abstract: The interest in air-water flows has not diminished in recent years, but it is accompanied by frequent citations of early, sometimes outdated articles. A basic issue is the inadequate, incomplete interpretation of air-water flow instrumentation by hydraulic engineers and researchers. This article comments on high-velocity air-water flow measurements by means of intrusive phase detection probes. This article focus on the bubbly flow structure of high-velocity air-water flow based upon measurements by means of intrusive phase detection probes. It is shown that some advanced post-processing techniques may yields expanded information on the air-water turbulent flow properties and bubbly flow structures. The outcomes demonstrate simple techniques in high-velocity air-water flow analysis.

Surface waves and roughness in self-aerated supercritical flow

Abstract: In high-velocity open channel flows, free-surface aeration is commonly observed. The effects of surface waves on the air-water flow properties are tested herein. The study simulates the air-water flow past a fixed-location phase-detection probe by introducing random fluctuations of the flow depth. The present model yields results that are close to experimental observations in terms of void fraction, bubble count rate and bubble/droplet chord size distributions. The results show that the surface waves have relatively little impact on the void fraction profiles, but that the bubble count rate profiles and the distributions of bubble and chord sizes are affected by the presence of surface waves.

Observations of the cross-wind velocity variance in the stable boundary layer

Abstract: Nine tower datasets over grassland, brush rangeland, snow covered plain, the ocean, three different pine forests, an aspen forest and an urban site, are used to document the scale-dependence of the cross-wind velocity variance in the stable boundary layer. The turbulence velocity variance scales with the surface momentum flux, as reported in previous studies. Such scaling removes the stability dependence of the variance at a given site, and also removes most of the differences between sites. The scaling is more effective with use of a record-dependent averaging time for defining the turbulent fluctuations. The variable averaging time is the timescale associated with the gap region in the heat flux multiresolution cospectra. On scales larger than turbulence and less than a few hours (mesoscale), variations in the cross-wind velocity variance at a given site are not related to local variables such as the friction velocity. Possible exceptions include suppression of turbulence and mesoscale motions in well-defined drainage flows and enhancement of turbulence and mesoscale motions in stronger winds downstream of a ridge. Larger mesoscale variance is associated with complex terrain and forested sites compared to the more homogeneous sites in flat terrain with short or no vegetation. These differences between sites are related to the absence of a gap region in the velocity spectra at the complex terrain and forested sites. The observed probability distribution functions of the total variance and the mesoscale variance are documented for different averaging times, stability classes and site characteristics.

Estimation of total sediment load concentration obtained by experimental study using artificial neural networks

Abstract: Estimation of sediment concentration in rivers is very important for water resources projects planning and managements. The sediment concentration is generally determined from the direct measurement of sediment concentration of river or from sediment transport equations. Direct measurement is very expensive and cannot be conducted for all river gauge stations. However, sediment transport equations do not agree with each other and require many detailed data on the flow and sediment characteristics. The main purpose of the study is to establish an effective model which includes nonlinear relations between dependent (total sediment load concentration) and independent (bed slope, flow discharge, and sediment particle size) variables. In the present study, by performing 60 experiments for various independent data, dependent variables were obtained, because of the complexity of the phenomena, as a soft computing method artificial neural networks (ANNs) which is the powerful tool for input–output mapping is used. However, ANN model was compared with total sediment transport equations. The results show that ANN model is found to be significantly superior to total sediment transport equations.

Development of an empirical equation for the transverse dispersion coefficient in natural streams

Abstract: In this study, a new empirical equation for the transverse dispersion coefficient has been developed based on the hydraulic and geometric parameters in natural streams using a regression technique. First, a total of 32 data sets in 16 streams were collected. Among those sets, 16 sets were used for deriving the new equation, and the other 16 sets were used for verifying the equation. Then, through dimensional analysis, it was found that the normalized transverse dispersion coefficient is associated with several parameters such as sinuosity, aspect ratio, and a friction term. The robust least square method was applied to estimate regression coefficients. The newly proposed equation was proven to be superior in explaining the dispersion characteristics of natural streams more precisely compared to the existing equations.

Large-scale turbulent structure of uniform shallow free-surface flows

Abstract: The hydrodynamics of super- and sub-critical shallow uniform free-surface flows are assessed using laboratory experiments aimed at identifying and quantifying flow structure at scales larger than the flow depth. In particular, we provide information on probability distributions of horizontal velocity components, their correlation functions, velocity spectra, and structure functions for the near-water-surface flow region. The data suggest that for the high Froude number flows the structure of the near-surface layer resembles that of two-dimensional turbulence with an inverse energy cascade. In contrast, although large-scale velocity fluctuations were also present in low Froude number flow its behaviour was different, with a direct energy cascade. Based on our results and some published data we suggest a physical explanation for the observed behaviours. The experiments support Jirka’s [Jirka GH (2001) J Hydraul Res 39(6):567–573] hypothesis that secondary instabilities of the base flow may generate large-scale two-dimensional eddies, even in the absence of transverse gradients in the time-averaged flow properties.

Comparison of dynamic subgrid-scale models for simulations of neutrally buoyant shear-driven atmospheric boundary layer flows

Abstract: Several non-dynamic, scale-invariant, and scale-dependent dynamic subgrid-scale (SGS) models are utilized in large-eddy simulations of shear-driven neutral atmospheric boundary layer (ABL) flows. The popular Smagorinsky closure and an alternative closure based on Kolmogorov’s scaling hypothesis are used as SGS base models. Our results show that, in the context of neutral ABL regime, the dynamic modeling approach is extremely useful, and reproduces several establised results (e.g., the surface layer similarity theory) with fidelity. The scale-dependence framework, in general, improves the near-surface statistics from the Smagorinsky model-based simulations. We also note that the local averaging-based dynamic SGS models perform significantly better than their planar averaging-based counterparts. Lastly, we find more or less consistent superiority of the Smagorinsky-based SGS models (over the corresponding Kolmogorov’s scaling hypothesis-based SGS models) for predicting the inertial range scaling of spectra.

Landslides into reservoirs and their impacts on banks

Abstract: Mass wasting processes, like slope failures, on the margins of dam reservoirs, lakes, bays and oceans may generate large water waves that can produce disasters due to flooding over the banks, run up along the shoreline and overtopping dam crests. Therefore, the study of slope failures, the subsequent generation of impulse waves and their consequences are of paramount importance for safety. In this paper the generation and propagation of water waves in reservoirs induced by landslides and their impact on banks were investigated by means of a laboratory study carried out at University of Coimbra wave channel, in a flume measuring 12.0 m × 1.5 m × 1.0 m (L × H × W), where two banks with variable slope were placed. The study considered the sliding of calcareous blocks over a sliding slope bank into the reservoir, the generation of impulse waves, their propagation in the reservoir and their impact on the downstream bank. A number of waves were generated by different fallings of calcareous blocks, considering different volumes, sliding slopes, initial positions and reservoir depths. All fallings were recorded by video-camera and the results were processed afterwards to obtain the time history of the falling. The water surface variations due to transient waves were measured at five gauges placed between the banks. The waves overtopping and breaking on the downstream bank were also filmed using a video camera, and the hydrodynamic forces on this bank were also measured using four pressure transducers.

A statistical model to forecast the profile of the index structure constant {C_N^2}

Abstract: A statistical study was performed over 145 profiles of meteorological balloons, equipped with microthermal sensors, from the ground to the midstratosphere. This study put into evidence the lognormal distribution of the fluctuations of the structure constant of temperature, the fluctuations of the buoyancy force and the vertical shear of the horizontal wind speed. We show that these quantities, computed over a large scale (100 m), are correlated up to the midstratosphere. A model is adjusted to estimate the optical turbulence strength $$C_N^2$$ from the macroscopic meteorological parameters. The model performances, to estimate and forecast the strength and the altitude of the optical turbulent layers, are quantified and compared with other already defined models. This model was shown to have the best performances and put into evidence a new relation to describe the turbulence on a large scale.

Energy dissipation estimates in oscillating grid setup: LDV and PIV measurements

Abstract: The dissipation of turbulent kinetic energy has been increasingly used as a scaling parameter to integrate microbiological accrual and metabolic rates with fluid-flow motion in natural and engineered aquatic ecosystems. The estimation of turbulent kinetic energy under field conditions and the generation of energy dissipation rates under controlled laboratory conditions with microbiological organisms are necessities required to integrate environmental/ecological laboratory protocols with a moving fluid in the environment. Turbulent fluid-flow conditions were generated in an oscillating grid setup, and turbulence variables were quantified using laser-Doppler velocimetry (LDV) and particle image velocimetry (PIV) measuring techniques. The rate of dissipation of the turbulent kinetic energy in the setup ranged from 10−9 to 10−4 m2/s3 and was similar to the levels of energy dissipation commonly reported in engineered and natural aquatic ecosystems. Time-averaged velocities were close to zero with the root-mean-square velocity ratios about 1, indicating nearly isotropic fluid-flow conditions in the setup. The velocity spectra, obtained by stationary LDV measurements for the vertical and horizontal velocity components across the setup revealed the existence of inertial subrange with the frequency power scaling law of “ω −5/3.” The estimated Eulerian frequency spectrum followed the theoretical functional relation and confirmed the applicability of inertial dissipation method for the estimation of turbulent kinetic energy dissipation rates. PIV was used for a direct estimation of dissipation by evaluating spatially distributed velocity gradients. The direct dissipation estimate in conjunction with the estimated Eulerian frequency spectrum provided evaluation of a “universal” constant, α, commonly used for the estimation of an energy dissipation rate over the inertial subrange of the Eulerian spectrum. The results demonstrated a range of values, rather than a universal constant, of α with a lognormal probability distribution for vertical and horizontal velocity components. In order to encompass a 0.955 probability range under the lognormal distribution $${({\frac{{\bar {\alpha}}}{{\sigma }^{2}} > {\alpha } > {\bar{\alpha}\sigma }^{2})}}$$ the universal constant, α, should be in the range 2.91 ≥ α u ≥ 0.43 and 4.44 ≥ α w ≥ 0.42 for horizontal and vertical velocity components, respectively.

The impact on tropospheric ozone formation on the implementation of a program for mobile emissions control : a case study in São Paulo, Brazil

Abstract: The main sources of reactive hydrocarbons (RHC) and nitrogen oxides (NOx), ozone precursors, in the Metropolitan Area of Sao Paulo (MASP) in the southeast of Brazil are emissions from vehicles fleets. Ambient surface ozone and particulate matter concentrations are air quality problem in the MASP. This study examined the impact that implementing a control program for mobile emissions has on ozone concentrations, An episode of high surface ozone concentrations occurring in the MASP during the March 13–15, 2000 period was used as a case study that was modeled for photochemical oxidants using the California Institute of Technology/Carnegie Mellon University three-dimensional photochemical model. Different scenarios were analyzed in relationship to the implementation of the Programa Nacional de Controle de Poluicao por Veiculos Automotores (PROCONVE, National Program to Control Motor Vehicle Pollution). Scenario 1 assumed that all vehicles were operating within PROCONVE guidelines. Scenarios 2 and 3 considered hypothetical situations in which the PROCONVE was not implemented. Scenario 2 set the premise that vehicles were using pre-1989 technology, whereas scenario 3 allowed for technological advances. A base case scenario, in which the official emission inventory for the year 2000 was employed, was also analyzed. The CIT model results show agreement with most measurements collected during 13–15 March 2000 modeling episode. Mean normalized bias for ozone, CO, RHC and NO x are approximately 9.0, 6.0, −8.3, 13.0%, respectively. Tropospheric ozone concentrations predicted for scenario 2 were higher than those predicted for scenarios 1, 3 and base case. This study makes a significant contribution to the evaluation of air quality improvement and provides data for use in evaluating the economic costs of implementing a program of motor vehicle pollution control aimed at protecting human health.

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

Abstract: 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 0 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.

Clustering of aerosols in atmospheric turbulent flow

Abstract: A mechanism of formation of small-scale inhomogeneities in spatial distributions of aerosols and droplets associated with clustering instability in the atmospheric turbulent flow is discussed. The particle clustering is a consequence of a spontaneous breakdown of their homogeneous space distribution due to the clustering instability, and is caused by a combined effect of the particle inertia and a finite correlation time of the turbulent velocity field. In this paper a theoretical approach proposed in Elperin et al. (2002) Phys Rev E 66:036302 is further developed and applied to investigate the mechanisms of formation of small-scale aerosol inhomogeneities in the atmospheric turbulent flow. The theory of the particle clustering instability is extended to the case when the particle Stokes time is larger than the Kolmogorov time scale, but is much smaller than the correlation time at the integral scale of turbulence. We determined the criterion of the clustering instability for the Stokes number larger than 1. We discussed applications of the analyzed effects to the dynamics of aerosols and droplets in the atmospheric turbulent flow.

An analytical solution of the advection-diffusion equation considering non-local turbulence closure

Abstract: Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection-diffusion equation. Moreover, large eddies are able to mix scalar quantities in a manner that is counter to the local gradient. We present a general solution of a two-dimension steady state advection-diffusion equation, considering non-local turbulence closure using the General Integral Laplace Transform Technique. We show some examples of applications of the new solution with different vertical diffusion parameterisations.

Quality control and flux gap filling strategy for Bowen ratio method: revisiting the Priestley–Taylor evaporation model

Abstract: Micrometeorological measurements were made over an irrigated rice paddy during the Huaihe River Basin Experiment (HUBEX) in 1999. This study addresses the quality control and gap filling strategy for the heat fluxes with the Bowen ratio energy balance (BREB) method. It also endeavors to benefit future studies through comparing five methodologies to estimate the net water exchange. First, a three-step quality control strategy is constructed. Its first two steps guarantee the correct flux directions and reject suspicious data, respectively. The third step forbids supersaturation by considering the Bowen ratio ranges for different flux combinations (termed as “scenarios”). The quality-controlled latent heat (LE) and sensible heat (H) fluxes fall in three scenarios, namely I (LE > 0, H > 0), II (LE > 0, H < 0), and IV (LE < 0, H < 0). Second, the Priestley–Taylor evaporation model (PTEM) is applied to fill the gaps of LE, while the energy balance relationship is used to gap-fill H, namely H = R N − G − LE (R N : net radiation; G: soil heat flux). Central to the success of this strategy is the idea to calibrate the Priestley–Taylor parameter (α) in a scenario-specific manner. On average, α values are calibrated as 1.20 and 1.35 for scenario I and II, respectively. For scenario IV, most α values lie in a narrow range, namely from 0.9 to 1.0. Then, α is calibrated as 0.97 to extend the applicability of the PTEM to condensation (negative LE). The scenario-specific treatment explicitly explains the diurnal variation of α derived without distinguishing the scenarios. Third, five methodologies are compared in the calculation of net water exchange, including PTEM-gap-filled BREB method (M1), energy-balance-based eddy covariance method (M2), and three simplified models to estimate LE, respectively as R N − H, R N − G, and R N . A major finding is that G and H have a similar effect of about 3% in the net water exchange. Thus, if either is neglected, the net evaporation should be slightly lowered to avoid overestimate, by 3% as a rule of thumb from this work.

River water quality management model using genetic algorithm

Abstract: Conventional mathematical programming methods, such as linear programming, non linear programming, dynamic programming and integer programming have been used to solve the cost optimization problem for regional wastewater treatment systems. In this study, a river water quality management model was developed through the integration of a genetic algorithm (GA). This model was applied to a river system contaminated by three determined discharge sources to achieve the water quality goals and wastewater treatment cost optimization in the river basin. The genetic algorithm solution, described the treatment plant efficiency, such that the cost of wastewater treatment for the entire river basin is minimized while the water quality constraints in each reach are satisfied. This study showed that genetic algorithm can be applied for river water quality modeling studies as an alternative to the present methods.

Renormalized numerical simulation of flow over planar and non-planar fractal trees

Abstract: We review the fundamentals of a new numerical modeling technique called Renormalized Numerical Simulation (RNS). The goal of RNS is to model the drag force produced by high Reynolds-number turbulent flow over objects that display scale-invariant properties, objects such as tree-like fractals. The hallmark of RNS in this application is that the drag of the unresolved tree branches is modeled using drag coefficients measured from the resolved branches and unresolved branches (as modeled in previous iterations of the procedure). In the present paper, RNS is used to study the effects of branch orientation on the drag force generated by highly idealized trees in which trunk and branches have square cross-section, and the branches all lie in a plane perpendicular to the incoming flow. Then, the procedure is generalized to the more general case of non-planar branch arrangements. Results illustrate that RNS may enable numerical modeling of environmental flow processes associated with fractal geometries using affordable computational resolution.

Experimental investigations on the stability of riprap layers on overtoppable earthdams

Abstract: In Germany small dams (<10 m) on flood retention reservoirs are often provided with an overtoppable dam section for flood relief. In case of overtopping a protection of the dam body is absolutely mandatory as the mainly cohesive dam material is not capable to withstand the affecting erosive forces of the flow. There is a big number of possibilities how to protect the dam whilst overtopping. Some of the most common construction types of slope protection layers for a sufficient dam protection such as riprap and placed stones, have been tested in large scale physical models at the Hydraulic Laboratory at Universitat Stuttgart. Thereby different experiments for the failure scenarios “erosion of single stones”, “sliding of the protection layer” as well as “disruption of the protection layer” have been conducted and analysed aiming to find some design criteria for the construction of such dams.

OHP02 gravity wave campaign in relation to optical turbulence

Abstract: Herein we present a campaign dedicated to the detection and the characterization of Gravity Waves (GW) in the Earth’s atmosphere in relation to the generation of Optical Turbulence (OT). The observations took place in France from 17 to 24 July 2002 at the Haute Provence Observatory (OHP) and simultaneously at the Sirene Observatory, some 20 km apart. From both sites, several balloons were launched that measured the classical PTU-Wind profiles and additionally the structure constant of the temperature field \(C_T^2\) vertical profiles. A Generalized Scidar (GS) technique was implemented at the 1.93 m-diameter OHP telescope, providing \(C_N^2(h)\) profiles every minute. From our observations, a significant amount of GW activity was observed at both sites, but without clear evidence of correlation between the two sites. It seems from our observations that a wide spectrum of GW is present at a given altitude and that this could result in a lack of correlation between observations made from two sites 20 km apart. Most GW are non-stationary with long horizontal wavelengths (λ ∼ 100–200 km), kilometric vertical wavelengths (λ ~ 0.5–2 km) and long intrinsic period (T ~ 2–15 h). They belong in the category of “hydrostatic rotating or non-rotating waves”. Layers of optical turbulence detected by balloons and the Scidar technique correlate well with regions of GW activity.

A performance comparison between nonlinear similarity functions in bulk parameterization for very stable conditions

Abstract: Three sets of nonlinear similarity functions for strong stability are selected to compare their performance in bulk parameterization. To uncover their advantages and disadvantages, theoretical and measurement analyses are made with four profile metrics and the Deacon number technique. Main disadvantages include the negligence of the different transfer efficiency between momentum and heat, the flux cutoff due to the upper limit in gradient Richardson number (Ri) and the ignorance of limited stability range where the dimensionless gradient functions (\({\varphi_{{\rm M}}}\) and \({\varphi_{{\rm H}}}\)) approach constants. Accordingly, three suggestions are made for future improvement. First, the functions for wind velocity and potential temperature should have the same function form, but with different coefficients. Second, \({\varphi_{\rm M}}\) and \({\varphi _{\rm H}}\) should approach constants only within a certain stability range. Third, the limit value in Ri should be avoided to widen their applicability in flux modeling. Furthermore, quantitative comparisons in transfer coefficients for moment and sensible heat (CD and CH) are made among the similarity functions in the bulk Richardson number (RiB) range 0 < RiB < 1. Generally, significant discrepancy is found, which may approach a factor of two and three at large RiB in CD and CH, respectively. Finally, a new recommendation is made to one of the three sets, mainly because of its ability to predict CD and CH that decrease rather slowly in very stable conditions.

Characterising strongly advected discharges in the initial dilution zone

Abstract: Mean concentration fields of strongly advected non-buoyant discharges are characterised with a double-Gaussian assumption. Comparisons with experimental data show that the approximation provides a reasonable representation of the cross-sectional profiles. The self-similarity of these profiles enables their form to be represented by two additional parameters, one describing the relative separation of the peaks and the other the ratio of the cross-sectional spreads. Values for these additional parameters are determined from experimental data. This systematic approach to characterising the strongly advected flows provides a consistent framework for determining spreading rates and concentration ratios, such as the peak to centreline maximum and the peak to top hat. The double-Gaussian framework also provides a basis for comparisons with the CorJet and VisJet numerical models. In addition the double-Gaussian assumption is employed to interpret data obtained using the Light Attenuation technique. This is a relatively simple measuring system, which provides depth integrated concentration information. The data obtained using this technique is shown to be generally consistent with that from previous studies.

Large-eddy simulation of low-level jet-like flow in a canopy

Abstract: The computational method of Large-Eddy Simulations has been used to study the weak, neutrally stable drainage flow within tree canopies. The computational results show that a secondary velocity maximum that resembles a jet is formed within the canopy under the nocturnal flow conditions. This jet-like flow is important in the analysis and measurements of the net ecosystem-atmosphere exchange (NEE) for carbon dioxide (CO2). A uniformly distributed, plane source was placed within the canopy in order to simulate the nocturnal production of CO2. The NEE is calculated as the sum of the integration of the rate of change of the concentration of CO2 over the computational domain, the vertical turbulent flux measured directly by eddy-covariance (EC) method, and the advection terms, which are not taken into account in the EC method. Numerical results of the velocity and concentration fields, within and above the canopy, are presented and their impact on the CO2 transport is investigated in detail. The computational results show that 15–20% of NEE is drained out by the advection process under the canopy. The results also show that the turbulent fluctuations in the lateral direction are also significant and may result in 2–5% CO2 transport.

The inertial range “outer scale” and optical turbulence

Abstract: The crucial parameter used to calculate turbulence effects upon light waves propagating through the atmosphere is known as the structure constant, $$c_n^2$$ . As Tatarski has shown, this parameter depends upon the “outer scale” of the inertial sub-range of the turbulence. Recently there have been successful predictions of astronomical “seeing” conditions at Mauna Kea Astronomical Observatory which have increased interest in this subject and in the use of the so-called “Dewan Optical Turbulence Model”. In the case of the Air Force, there has been a longstanding need for such optical turbulence prediction, especially in the stratosphere. In the past researchers have used a relation due to Tatarski, (which plays a prominent role in this model) in order to deduce values of the “outer scale” from $$c_n^2$$ measurements. When doing this, they have been surprised to find values very much smaller than expected. The goal of the paper is to explain this unexpected result. As we will show, this result can be explained by two factors: (a) the average turbulent layer thicknesses are smaller than originally believed, and, more importantly, (b) only a minor fraction of the stratosphere is turbulent. In order to arrive at this conclusion, we used the high-resolution (10 m) wind profiles that were originally used to formulate the previously mentioned optical turbulence model.

Automatic sensitivity analysis of a finite volume model for two-dimensional shallow water flows

Abstract: Given a numerical model for solving two-dimensional shallow water equations, we are interested in the robustness of the simulation by identifying the rate of change of the water depths and discharges with respect to a change in the bottom friction coefficients. Such a sensitivity analysis can be carried out by computing the corresponding derivatives. Automatic differentiation (AD) is an efficient numerical method, free of approximation errors, to evaluate derivatives of the objective function specified by the computer program, Rubar20 for example. In this paper AD software tool Tapenade is used to compute forward derivatives. Numerical tests were done to show the robustness of the model and to demonstrate the efficiency of these AD-derivatives.

Estimation of far-field horizontal and vertical turbulent diffusion coefficients from the concentration field of a wastewater plume near the Akashi Strait

Abstract: We carried out a field study of the plume discharged by a near-shore wastewater outfall near the Akashi Strait, Japan. Using an Acoustic Doppler Current Profiler and a tow-body CTD, we measured the near-surface salinity and temperature fields in the region throughout an M2 tidal cycle. We filtered the data in T–S space to remove water masses other than the wastewater, and then used the adiabatic mixing assumption to calculate the concentration of wastewater in the far field of this plume. Averaging the T–S fields of repeated surveys over a time period during which the tidal regime did not change substantially, allowed comparison of the time-averaged plume with the analytical solution for a plume diffusing in both the horizontal and vertical dimensions. The resulting vertical turbulent diffusion coefficients agreed well with those resulting from Thorpe scales determined via a vertically-profiling CTD, as well as with the canonical value for open channel flow of Dz = 0.067hu*. The corresponding horizontal turbulent diffusion coefficients, however, were two orders of magnitude larger than those typically observed in straight channels, and an order of magnitude larger than those observed in meandering rivers. This is likely a result of enhanced horizontal mixing due to barotropic eddies generated by the interaction of strong tidal flow with headlands and levees, as well as due to the time-varying nature of tidal flow, and baroclinic spreading of the buoyant wastewater plume.

Analysis of a two-layer cloud system with RAMS model and comparison to airborne observations

Abstract: A three-dimensional numerical model (Regional Atmospheric Modeling System—-RAMS) was used to study the formation and evolution of water forms in a two-layer cloud structure observed during a field campaign over Brest (France). The model performance in regular operations, using conventional meteorological data as initial and lateral boundary conditions, was also examined. Remote sensing observations of the cloud system and in-situ aircraft data, selected during the campaign, were used to validate the model outputs. The model simulations showed that the lower cloud formation was characterized by high number concentration of pristine ice and snow, while the concentration of aggregates, graupel and hail were considerably lower. Hydrometeors in liquid phase appeared demonstrating high number concentration and water content on the top of this layer. The upper cloud layer consisted only of frozen water substances in lower amounts. The qualitative and quantitative comparison of the model-calculated meteorological and microphysical fields to the available observational data revealed that the model reflected fairly well the cloud structure (e.g., the spatio-temporal variability of the cloud parameters, the geometry of the cloud system). However, there were deviations as far as the model underestimating the ice water content (IWC) and number concentration (Nt) fields is concerned, especially at the atmospheric layer between 2.5 and 4 km of altitude. These deviations of the model simulated quantities from the measured ones may be attributed either to the performance of the model’s microphysics scheme, to instrument inaccuracies and to the local disturbance caused by the aircraft.

OHP-APT 2002 Gravity Wave campaign : Waves, turbulence and forecasts

Abstract: In July 2002, a multi-sensor campaign was conducted in southern France to investigate the hypothesized connection between gravity waves and optical turbulence. A generalized scidar (GS) was mounted on the 1.93 m diameter telescope at l’Observatoire de Haute-Provence (OHP). The GS provides continuous profiles of optical turbulence with 300 m vertical resolution from telescope altitude up to 25 km. Thermosondes, which provide in situ measurement of optical turbulence by measuring temperature variance, were launched at the OHP site and at a site approximately 20 km west-northwest of OHP. Gravity wave activity was deduced from temperature and wind velocity measured by radiosondes, which are part of the thermosonde system. In this paper, gravity waves were analyzed using techniques for simple two-dimensional mountain waves, with only fair results. Mesoscale models were run at moderately high resolution for the period. The forecasts were analyzed for wave activity, and a post processor model was used to diagnose the optical turbulence. Mountain waves were evident in the forecast, but quantitative comparison showed the forecast to be inadequate in predicting wave strength. The forecast optical turbulence was in fair agreement with measurements with notable exceptions.