Numerical simulation of particle-driven gravity currents

Particle-driven gravity currents frequently occur in nature, for instance as turbidity currents in reservoirs. They are produced by the buoyant forces between fluids of different density and can introduce sediments and pollutants into water bodies. In this study, the propagation dynamics of gravity currents is investigated using the FLOW-3D computational fluid dynamics code. The performance of the numerical model using two different turbulence closure schemes namely the renormalization group (RNG) ${k-\epsilon}$ scheme in a Reynold-averaged Navier-Stokes framework (RANS) and the large-eddy simulation (LES) technique using the Smagorinsky scheme, were compared with laboratory experiments. The numerical simulations focus on two different types of density flows from laboratory experiments namely: Intrusive Gravity Currents (IGC) and Particle-Driven Gravity Currents (PDGC). The simulated evolution profiles and propagation speeds are compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting the temporal and spatial evolution of intrusive gravity currents. In particular, the simulated propagation speeds are in excellent agreement with experimental results. The simulation results do not show any considerable discrepancies between RNG ${k-\epsilon}$ and LES closure schemes. The FLOW-3D model coupled with a particle dynamics algorithm successfully captured the decreasing propagation speeds of PDGC due to settling of sediment particles. The simulation results show that the ratio of transported to initial concentration C _o /C _i by the gravity current varies as a function of the particle diameter d _s . We classify the transport pattern by PDGC into three regimes: (1) a suspended regime (d _s is less than about 16 μm) where the effect of particle deposition rate on the propagation dynamics of gravity currents is negligible i.e. such flows behave like homogeneous fluids (IGC); (2) a mixed regime (16 μm < d _s <40 μm) where deposition rates significantly change the flow dynamics; and (3) a deposition regime (d _s ?> 40 μm) where the PDGC rapidly loses its forward momentum due to fast deposition. The present work highlights the potential of the RANS simulation technique using the RNG ${k-\epsilon}$ turbulence closure scheme for field scale investigation of particle-driven gravity currents.     

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.     

Detailed experimental study of hydrodynamic turbulent flows generated in vertical slot fishways

The kinematics of hydrodynamic turbulent flows developed in vertical slot fishways (VSF) was studied in detail in flow patterns not yet published to date for the purposes of modifying existing devices and to allow for the passage of all fishes, particularly the smaller species. A transparent device based on the typical prototype dimensions of VSF in France was constructed for the experiment. The velocity measurements were carried out by Particle Image Velocimetry (PIV). These measurements were used to determine the various kinematics parameters characterizing the flow. From the dimensions and slope of the fishway, two flow topologies highlighting the swirling pattern were proposed. The method of Proper Orthogonal Decomposition (POD) was used to undertake unsteady and energetic analyses to characterize the main phases of flow evolution that fish passing through the passage may encounter.     

Numerical modelling of horizontal sediment-laden jets

Sediment-laden turbulent flows are commonly encountered in natural and engineered environments. It is well known that turbulence generates fluctuations to the particle motion, resulting in modulation of the particle settling velocity. A novel stochastic particle tracking model is developed to predict the particle settling out and deposition from a sediment-laden jet. Particle velocity fluctuations in the jet flow are modelled from a Lagrangian velocity autocorrelation function that incorporates the physical mechanism leading to a reduction of settling velocity. The model is first applied to study the settling velocity modulation in a homogeneous turbulence field. Consistent with basic experiments using grid-generated turbulence and computational fluid dynamics (CFD) calculations, the model predicts that the apparent settling velocity can be reduced by as much as 30 % of the stillwater settling velocity. Using analytical solution for the jet mean flow and semi-empirical RMS turbulent velocity fluctuation and dissipation rate profiles derived from CFD predictions, model predictions of the sediment deposition and cross-sectional concentration profiles of horizontal sediment-laden jets are in excellent agreement with data. Unlike CFD calculations of sediment fall out and deposition from a jet flow, the present method does not require any a priori adjustment of particle settling velocity.     

A quantitative study of intertidepool movement of the wooly sculpin Clinocottus analis

Movements of the wooly sculpin Clinocottus analis between all pools in a designated area of a southern California (USA) intertidal zone were studied. Fish were captured during low tides, finclipped or bead-tagged, and released in the pool of capture. Percent recaptured in pool of first capture declined from 30 after 2 weeks to about 5 after 20 weeks. Percent recaptured in any pool in the study area declined from 50 after 2 weeks to 20 after 8 weeks, but remained at 20 through 16 weeks. A higher fidelity to an area of the intertidal zone than to a specific pool was suggested. Large fish exhibited higher probability of recapture than smaller individuals, indicating they become more sedentary with age. Changes in pool topography influenced percent recurrence in a pool, most likely as a consequence of changes in quality and quantity of cover available. Such responses to changes in specific pools and fidelity to intertidal areas containing several pools appear to be advantageous behavioral patterns in an environment where topographic change is frequent and often massive.     

Efficient numerical computation and experimental study of temporally long equilibrium scour development around abutment

For the abutment bed scour to reach its equilibrium state, a long flow time is needed. Hence, the employment of usual strategy of simulating such scouring event using the 3D numerical model is very time consuming and less practical. In order to develop an applicable model to consider temporally long abutment scouring process, this study modifies the common approach of 2D shallow water equations (SWEs) model to account for the sediment transport and turbulence, and provides a realistic approach to simulate the long scouring process to reach the full scour equilibrium. Due to the high demand of the 2D SWEs numerical scheme performance to simulate the abutment bed scouring, a recently proposed surface gradient upwind method (SGUM) was also used to improve the simulation of the numerical source terms. The abutment scour experiments of this study were conducted using the facility of Hydraulics Laboratory at Nanyang Technological University, Singapore to compare with the presented 2D SGUM–SWEs model. Fifteen experiments were conducted with their scouring flow durations vary from 46 to 546 h. The comparison shows that the 2D SGUM–SWEs model gives good representation to the experimental results with the practical advantage.     

Monitoring and estimating the flow conditions and fish presence probability under various flow conditions at reach scale using genetic algorithms and kriging methods

The combination of current velocity and water depth influences stream flow conditions, and fish activities prefer particular flow conditions. This study develops a novel optimal flow classification method for identifying types of stream flow based on the current velocity and the water depth using a genetic algorithm. It is applied to the Datuan stream in northern Taiwan. Fish were sampled and their habitat investigated at the study site during the spring, summer, fall and winter of 2008-2009. The current velocity, water depth and maps of the presence probability of fish were estimated by ordinary and indicator kriging. The optimal classification results were compared with the classification results obtained using the Froude number and empirical methods. The flow classification results demonstrate that the proposed optimal flow classification method that considers depth-velocity and optimally identified criteria for classifying flow types, yields a current velocity and water depth of 0.32 (m/s) and 0.29 (m), respectively, and classifies the flow conditions in the study area as pool, run, riffle and slack. The variography results of the current velocity and the water depth data reveal that seasonal flows are not spatially stationary among seasons in the study area. Kriging methods and a two-dimensional hydrodynamic model (River 2D) with empirical and optimal flow classification methods are more effective than the Froude number method in classifying flow conditions in the study area. The flow condition classifications and probability maps were generated by River 2D, ordinary kriging and indicator kriging, to quantify the flow conditions preferred by Sicyopterus japonicus in the study area. However, the proposed optimal classification method with kriging and River 2D is an effective alternative method for mapping flow conditions and determining the relationship between flow and the presence probability of target fish in support of stream restoration.     

Study on the concentration distribution in a trapezoidal open-channel flow with a side discharge

The number and distribution of pollutant concentration in a trapezoidal open channel flow with a side discharge is calculated and effects of the bank gradient are investigated in this paper. A sigma-coordinate water quality numerical model is used to simulate the process of both water and pollutant transportation in the trapezoidal channel open flow. The diffusion coefficient used in the prediction is determined by two methods including constant coefficient and the depth-averaged k-epsilon turbulence closure model. The change of the concentration with the bank gradient is acquired based on the simulation of cases with different bank gradients. An analytical formula is derived by using the mirror image method and related diffusion theories, ignoring the discharge momentum and the influence of the opposite bank. The formula can predict the number and distribution of pollutant concentration with some acceptable errors. The results demonstrate that the bank gradient has great influence on the concentration distribution which will decrease with the increase of the bank gradient approximately following a hyperbolic law.     

2D numerical analysis of the influence of near-bank vegetation patches on the bed morphological adjustment

This paper investigates the influence of near-bank vegetation patches on the bed morphological adjustment in open channel flow systems. The 2D depth-averaged hydro-morphological model is adopted for this investigation, which is first validated by laboratory experimental data measured in an open channel with a single near-bank vegetation patch. The validated model is then applied for extensive numerical simulations, with the aim of conducting a systematic analysis of the influence of different geometric controlling parameters on the bed morphological evolution. The controlling parameters taken into account for numerical analysis include the angle of repose value (RAV) of sediment, vegetation density (VD), patch length (PL) and patch width (PW). The numerical results and analysis show that: (1) the RAV of sediment with slope-failure parametrization only influences the shape of the transverse bed topography in the junction region; (2) increase in VD, PL and PW that substantially enhances flow blockage effect encourages the growth of the pool adjacent to the patch in three dimensions; (3) increase in VD, PL and PW produces analogous retrogressive erosion (erosion toward the upstream) in the pool region, presumably due to the increase in flow resistance. Additional numerical experiments suggest that the staggered-order distribution of multiple patches might be an optimal choice for channel restoration and conservation since pools and riffles with larger scales can be produced.

     

Measurements and modeling of open-channel flows with finite semi-rigid vegetation patches

The hydrodynamics of flows through a finite length semi-rigid vegetation patch (VP) were investigated experimentally and numerically. Detailed measurements have been carried out to determine the spatial variation of velocity and turbulence profiles within the VP. The measurement results show that an intrusion region exists in which the peak Reynolds stress remains near the bed. The velocity profile is invariant within the downstream part of the VP while the Reynolds stress profile requires a longer distance to attain the spatially invariant state. Higher vegetation density leads to a shorter adjustment length of the transition region, and a higher turbulence level within the VP. The vegetation density used in the present study permits the passing through of water and causes the peak Reynolds stress and turbulence kinetic energy each the maximum at the downstream end of the patch. A 3D Reynolds-averaged Navier–Stokes model incorporating the Spalart–Allmaras turbulence closure was employed subsequently to replicate the flow development within the VP. The model reproduced transitional flow characteristics well and the results are in good agreement with the experimental data. Additional numerical experiments show that the adjustment length can be scaled by the water depth, mean velocity and maximum shear stress. Empirical equations of the adjustment lengths for mean velocity and Reynolds stress were derived with coefficients quantified from the numerical simulation results.     

Site fidelity and homing behaviour in intertidal fishes

The theory of ecological cognition poses that the brains and behaviour of animals are shaped by the environmental challenges they face in their everyday lives. Site fidelity and homing ability was tested for five species of intertidal rock pool fish by tagging and displacing them to new rock pools at various distances from their ‘home’ rock pools. Three of the species were rock pool specialists whilst the remaining two spend a small proportion of their life in rock pools during early ontogeny. The three specialists showed strong site fidelity with >50 % of individuals found in the same pool 42 days after tagging. In contrast, the non-specialist species showed low fidelity and poor homing abilities. Homing success in the rock pool specialists remained relatively stable as displacement distance increased. The effect of body size on homing ability was species dependent, with only one species showing a significantly greater tendency to home with increasing size.     

Simulation of exchange flow between open water and floating vegetation using a modified RNG <Emphasis Type="Italic">k</Emphasis>-<Emphasis Type="Italic">ε</Emphasis> turbulence model

The paper focuses on the numerical simulation of the exchange flow between open water and floating vegetation, which plays an important role in maintaining the ecological balance by transporting nutrient matter. The simulation was conducted using a new solver developed upon OpenFOAM. A modified RNG k-ε turbulence model, which is expected to model both the high- and low-Reynolds number flows correctly, was used to determine the eddy viscosity. Several particular terms were added into the momentum equations and turbulence model equations to model the effects of vegetation and buoyancy. Among these terms, the term for the effect of vegetation in the ε-equation was re-modelled. The model was validated by properly predicting the profiles of mean velocity and turbulent kinetic energy for flows through suspended canopies. The density flow between open and vegetated water was simulated with the same conditions as those of the experiment conducted by Zhang and Nepf. The predicted results agreed well with the experimental data and provided more detailed information of such exchange flow. The convection between the root layer and the layer beneath the roots, which was not observed in the experiment, was observed in the numerical simulation.     

Elevated CO2 stimulates net accumulations of carbon and nitrogen in land ecosystems: a meta-analysis

The capability of terrestrial ecosystems to sequester carbon (C) plays a critical role in regulating future climatic change yet depends on nitrogen (N) availability. To predict long-term ecosystem C storage, it is essential to examine whether soil N becomes progressively limiting as C and N are sequestered in long-lived plant biomass and soil organic matter. A critical parameter to indicate the long-term progressive N limitation (PNL) is net change in ecosystem N content in association with C accumulation in plant and soil pools under elevated CO2. We compiled data from 104 published papers that study C and N dynamics at ambient and elevated CO2. The compiled database contains C contents, N contents, and C:N ratio in various plant and soil pools, and root:shoot ratio. Averaged C and N pool sizes in plant and soil all significantly increase at elevated CO2 in comparison to those at ambient CO2, ranging from a 5% increase in shoot N content to a 32% increase in root C content. The C and N contents in litter pools are consistently higher in elevated than ambient CO2 among all the surveyed studies whereas C and N contents in the other pools increase in some studies and decrease in other studies. The high variability in CO2-induced changes in C and N pool sizes results from diverse responses of various C and N processes to elevated CO2. Averaged C:N ratios are higher by 3% in litter and soil pools and 11% in root and shoot pools at elevated relative to ambient CO2. Elevated CO2 slightly increases root:shoot ratio. The net N accumulation in plant and soil pools at least helps prevent complete down-regulation of, and likely supports, long-term CO2 stimulation of C sequestration. The concomitant C and N accumulations in response to rising atmospheric CO2 may reflect intrinsic nature of ecosystem development as revealed before by studies of succession over hundreds to millions of years.     

Effect of environmental flows on deep,anoxic pools

A series of field surveys were carried out on two permanent pools of the upper Glenelg River in SW Victoria, Australia. One was representative of the wider and deeper pools while the other was representative of the more-narrow and shallower pools. Both pools showed a typical seasonal cycle of warm, brackish, oxygen-poor, summer conditions and cool, oxygen-rich, low-salinity, winter conditions. The summer salinity increases were larger than expected, suggesting possible saline groundwater inflow from unidentified springs. Both pools contained anoxic water in their deeper sections but this was permanent only in the deeper pool. A simple model of the flushing rate of such anoxic pools subject to flows, such as environmental flow releases, was developed, based on an energy balance between the potential energy required to lift the anoxic layer and the kinetic energy derived from the river flow. The results were tested against and in agreement with the field measurements. The model also suggests that the anoxic layers are resilient to all but the largest environmental flows.     

Modelling nature-like fishway flow around unsubmerged obstacles using a 2D shallow water model

In the scope to create efficient nature like fish ramps using large-scale roughness elements, the present study is an audit of modelling such complex 3D free surface flows using an industrial 2D code solving shallow water equations. Validation procedure is based upon the comparison between numerous experimental measurements and numerical runs around large-scale roughness patterns disposed on the flume bottom in order to determine what 2D reliable numerical results can be expected. In this paper, we focused on cases of unsubmerged obstacles. The results demonstrate that 2D shallow water modelling using an industrial code such as TELEMAC-2D can be a convenient way for the hydraulic engineer to help design a nature-like fishway. This article emphasizes the limitations due to 2D depth integration of velocities and turbulence modelling and gives the domain of validity of the method.     

A comparison of spatial learning and memory capabilities in intertidal gobies

For the majority of animals, the ability to orient in familiar locations is a fundamental part of life, and spatial memory allows individuals to remember key locations such as food patches, shelter, mating sites or areas regularly occupied by predators. This study determined if gobies collected from rocky platforms and sandy beaches differ in their ability to learn and memorise the locations of tide pools in a simulated rocky intertidal zone. Intertidal rock pool gobies show acute homing abilities and, therefore, should be expected to display superior learning and memory capabilities. In contrast, it is unlikely that natural selection would favour similar learning skills in sandy shore fishes because they simply shift back and forth with the tides. The learning abilities of gobies were tested using small replica rock platforms, containing four tide pools that retained varying depths of water at simulated low tide. Gobies were categorised as having learnt the task if they were able to consistently locate the tide pool that retained the most water at simulated low tide as the pool with the most favourable conditions. Rock pool species were able to locate the deepest pool to wait out low tide for ~95 % of the trials, while species from sandy shores were found in the deepest pool ~10 % of trials. Despite repeated stranding, sandy shore fish continued to follow the tide out. We propose that rock pool species memorised the location of rock pools during simulated high tide enabling them to relocate the best refuge for low tide.     

Litter and dead wood dynamics in ponderosa pine forests along a 160-year chronosequence.

Disturbances such as fire play a key role in controlling ecosystem structure. In fire-prone forests, organic detritus comprises a large pool of carbon and can control the frequency and intensity of fire. The ponderosa pine forests of the Colorado Front Range, USA, where fire has been suppressed for a century, provide an ideal system for studying the long-term dynamics of detrital pools. Our objectives were (1) to quantify the long-term temporal dynamics of detrital pools; and (2) to determine to what extent present stand structure, topography, and soils constrain these dynamics. We collected data on downed dead wood, litter, duff (partially decomposed litter on the forest floor), stand structure, topographic position, and soils for 31 sites along a 160-year chronosequence. We developed a compartment model and parameterized it to describe the temporal trends in the detrital pools. We then developed four sets of statistical models, quantifying the hypothesized relationship between pool size and (1) stand structure, (2) topography, (3) soils variables, and (4) time since fire. We contrasted how much support each hypothesis had in the data using Akaike's Information Criterion (AIC). Time since fire explained 39-80% of the variability in dead wood of different size classes. Pool size increased to a peak as material killed by the fire fell, then decomposed rapidly to a minimum (61-85 years after fire for the different pools). It then increased, presumably as new detritus was produced by the regenerating stand. Litter was most strongly related to canopy cover (r2 = 77%), suggesting that litter fall, rather than decomposition, controls its dynamics. The temporal dynamics of duff were the hardest to predict. Detrital pool sizes were more strongly related to time since fire than to environmental variables. Woody debris peak-to-minimum time was 46-67 years, overlapping the range of historical fire return intervals (1 to > 100 years). Fires may therefore have burned under a wide range of fuel conditions, supporting the hypothesis that this region's fire regime was mixed severity.     

RANS simulation of neutral atmospheric boundary layer flows over complex terrain by proper imposition of boundary conditions and modification on the k-ε model

In this study, a modelling methodology is proposed for RANS simulations of neutral Atmospheric Boundary Layer (ABL) flows on the basis of the standard k-ε model, which allows the adoption of an arbitrary shear stress model. This modelling methodology is first examined in the context of an open flat terrain in an empty domain to ascertain there are no substantial changes in the prescribed profiles. The results show that relatively good homogeneity can be achieved with this modelling methodology for various sets of inflow boundary profiles. In addition, to extend the solutions derived from the standard k-ε model to RNG k-ε model, the RNG k-ε model is in detail assembly and tuned. Finally, the topographic effects on surface wind speeds over a complex terrain are assessed with the combined use of the proposed methodology and the modified RNG model. The numerical results are in good agreement with wind tunnel testing results and long-term field observations. A discussion of the effects of horizontal homogeneity and turbulence models on the simulated wind flows over a complex terrain is also given.     

Flow and drag force around a free surface piercing cylinder for environmental applications

This paper investigates flows around a free surface piercing cylinder with Froude number F > 0.5 and Reynolds number around Re = 50,000. The aim of this work is to gain a better understanding of the flow behaviour in environmental systems such as fishways. The advances are based upon experimental and numerical results. Several flow discharges and slopes are tested to obtain both subcritical and supercritical flows. The drag force exerted on the cylinder is measured with the help of a torque gauge while the velocity field is obtained using particle velocimetry. For the numerical part, two URANS turbulence models are tested, the k-\(\omega\) SST and the RNG k-\(\varepsilon\) models using the OpenFOAM software suite for subcritical cases, and then compared with the corresponding experimental results. With fishways applications in mind, the changes in drag coefficient \(C_d\) versus Froude number and water depth are studied and experimental correlations proposed. We conclude that the most suitable URANS turbulence model for reproducing this kind of flow is the k-\(\omega\) SST model.     

The size and function of the internal inorganic carbon pool of the foraminifer Amphistegina lobifera

The existence of an internal inorganic carbon pool in the perforate foraminifer Amphistegina lobifera, as suggested recently (ter Kuile and Erez 1987), has been established by direct measurements using a new ¹⁴C tracer method. The imperforate species Amphisorus hemprichii does not contain such a pool. The size of the pool in A. lobifera is proportional to its calcification rate and approximately equals the amount of carbon incorporated into the skeleton during 24 h. Time course experiments show that inorganic carbon (Ci) is photoassimilated at constant rates by the algal symbionts, that the pool is filled to maximum capacity in ca. 24 h, and that Ci incorporation into the skeleton starts only after the pool is filled up. During the chase phase of pulsechase experiments, all ¹⁴C initially residing in the pool is transferred to the skeleton, indicating that the pool serves for calcification and not for photosynthesis. Uptake of Ci into the pool occurs only in the light, indicating that energy may be required for this process. Furthermore, calculations of the Ci concentration inside the pool suggest that it is higher by 2 to 3 orders of magnitude compared to seawater concentration, suggesting that its accumulation is an energy dependent process.