DISCONTINUOUS FLOW OF TURBID DENSITY CURRENTS Ⅱ. INTERNAL HYDRAULIC JUMP

Traveling and stationary internal hydraulic jumps in density currents with positive or negative entrainment coefficients were analyzed based on simple assumptions. An expression of internal hydraulic jumps with entrainment coefficients was derived. Experimental data, published in literature, of stationary internal hydraulic jumps in turbid, thermal and saline density currents including measured values of water entrainment were used to compare with theory. Comparison was also made of traveling internal hydraulic jumps between measured data and theory.     

DISCONTINUOUS FLOW OF TURBID DENSITY CURRENTS II.INTERNAL HYDRAULIC JUMP

1 INTERNAL HYDRAULIC JUMPS AND LOCAL WATER ENTRAINMENT Traveling and stationary internal hydraulic jumps in two-layer flow can occur in the following situations (Fig. 1). When a stratified fluid flows over an obstacle created by a sudden blockage, the layer of increased depth within the heavier fluid will propagate upstream as a traveling hydraulic jump. When the heavier fluid overtops the obstacle it will flow along the downstream face of the obstacle as a supercritical flow…     

MODELING UNSTEADY TURBIDITY CURRENTS IN A STEEP NARROW RESERVOIR

1.IN~DUCTIONTurbiditycurrentisoneclassofflowsnameddensitycurrentorgravitycurrent(therHunterRouse(Yih(1980)),whichdePictstheintmsionofheaVyfluidintoalighterone.Usually,thedensitydifferencebetWeentWonuidisrelativelysmallandmixingacrosstheimerfaceoccurs.ThedrivingforceofdensitycurrentsisnotdensitydifferenceitselfbutthedifferenceinspeCmcweights.Turbiditycurrentisnamedwhenthedensitydifferenceisespeciallycausedbysuspendedfinesedimentparticles.Sincesediment-ladenflowcaninteraCtwiththelowerbou…     

Impacts of climate change and land use change on runoff of forest catchment in northeast China

Hydrological processes change from the impacts of climate variability and human activities. Runoff in the upper reaches of the Hun‐Taizi River basin, which is mainly covered by forests in northeast China, decreased from 1960 to 2006. The data used in this study were based on runoff records from six hydrological stations in the upper reaches of the Hun‐Taizi River basin. Nonparametric Mann–Kendall statistic was used to identify change trends and abrupt change points and consequently analyze the change characteristics in hydrological processes. The abrupt change in the annual runoff in most subcatchments appeared after 1975. Finally, the effects of climate change and land cover change on water resources were identified using regression analysis and a hydrology model. Results of the regression analysis suggest that the correlation coefficients between precipitation and runoff prior to the abrupt change were higher compared with those after the abrupt change. Moreover, using hydrology model analysis, the water yield was found to increase because of the decrease in forest land. Copyright © 2012 John Wiley & Sons, Ltd.     

Interactions among waves,current, and mud: Numerical and laboratory studies

Interactions between waves, current, mud and turbulence are very complicated in the coastal and estuarine turbid waters. It is still necessary to improve our understanding of the fundamental physical processes governing the cohesive sediment transport in the coastal and estuarine waters. A numerical model is developed to study the interactions among waves, current, and mud. An eddy viscosity model for wave and current is proposed in order to close the equations of wave motion or of current motion in a combined flow, respectively. The equations of mud transport are derived based on the visco-elastic properties of mud. Coupling the equations of wave motion or of current motion for water layer with those of mud layer can give (1) wave height; (2) distributions of current velocities in the water layer; (3) distributions of transport velocities at the water–mud interface; and (4) distributions of mass transport velocities within the mud layer. These modeled results are in a reasonable agreement with experimental results. Results suggest that (1) the rate of wave attenuation increases in the opposing currents (currents against in the direction in which the waves propagate) and decreases in the following currents (currents in the same direction as that in which the waves propagate); (2) the opposing currents would have more significant effects on the rate of wave height attenuation than the following currents; (3) the effect of current on the rate of wave attenuation on the muddy bottom is larger than that on the rigid bottom; (4) mud transport rate increased in the following currents but decreased in the opposing currents; and (5) the rate of wave height attenuation on the mud bottom is one order of magnitude larger than that on the rigid bottom.     

Dispersal and air entrainment in unconfined dilute pyroclastic density currents

Unconfined scaled laboratory experiments show that 3D structures control the behavior of dilute pyroclastic density currents (PDCs) during and after liftoff. Experiments comprise heated and ambient temperature 20 μm talc powder turbulently suspended in air to form density currents within an unobstructed 8.5?×?6?×?2.6-m chamber. Comparisons of Richardson, thermal Richardson, Froude, Stokes, and settling numbers and buoyant thermal to kinetic energy densities show good agreement between experimental currents and dilute PDCs. The experimental Reynolds numbers are lower than those of PDCs, but the experiments are fully turbulent; thus, the large-scale dynamics are similar between the two systems. High-frequency, simultaneous observation in three orthogonal planes shows that the currents behave very differently than previous 2D (i.e., confined) currents. Specifically, whereas ambient temperature currents show radial dispersal patterns, buoyancy reversal, and liftoff of heated currents focuses dispersal along narrow axes beneath the rising plumes. The aspect ratios, defined as the current length divided by a characteristic width, are typically 2.5–3.5 in heated currents and 1.5–2.5 in ambient temperature currents, reflecting differences in dispersal between the two types of currents. Mechanisms of air entrainment differ greatly between the two currents: entrainment occurs primarily behind the heads and through the upper margins of ambient temperature currents, but heated currents entrain air through their lateral margins. That lateral entrainment is much more efficient than the vertical entrainment, >0.5 compared to ～0.1, where entrainment is defined as the ratio of cross-stream to streamwise velocity. These experiments suggest that generation of coignimbrite plumes should focus PDCs along narrow transport axes, resulting in elongate rather than radial deposits.     

Strategies for managing reservoir sedimentation

Sediment deposition in reservoirs has caused the loss of 66% of the reservoir capacity in China. The main sedimentation control strategies are： 1） storing the clear water and releasing the turbid water; 2） releasing turbidity currents; 3） Draw-down flushing and empty flushing; and 4） dredging, The paper summarizes these strategies with examples. Sediment transport in many Chinese rivers occurs mostly during the 2-4 month flood season, that is, 80-90% of the annual sediment load is transported with 50-60% of the annual runoff. By storing the clear water after the flood season and releasing the turbid water during the flood season, less sediment deposits in the reservoir while the reservoir is still able to store enough water for power generation in the low flow season. The Three Gorges and Sanmenxia reservoirs apply this strategy and control sedimentation effectively. Turbidity currents have become the main sedimentation control strategy for the Xiaolangdi Reservoir. Empty flushing involves reservoir draw-down to temporarily establish riverine flow along the impound reach, flushing the eroded sediment through the outlets. Case studies with the Hengshan Reservoir and Zhuwo Reservoir are presented. Jet dredgers have been used to agitate the reservoir deposit so that the deposit is released from the reservoir with currents. The sediment releasing efficiency is 30-100% for storing the clear and releasing the turbid; 6-40 % for turbidity current; and 2,400-5,500% for empty flushing. Empty flushing causes high ecological stress on the ecosystem to the downstream reaches. Storing the clear and releasing the turbid is the best strategy to control reservoir sedimentation while achieving hydro-power benefit and maintain ecological stability.     

Influence of shallowness,bank inclination and bank roughness on the variability of flow patterns and boundary shear stress due to secondary currents in straight open-channels

Boundary shear stress and flow variability due to its interaction with main flow and secondary currents were investigated under conditions that extend previous research on trapezoidal channels. Secondary currents that scale with the flow depth were found over the entire width in all experiments. These findings contradict the widespread perception that secondary currents die out at a distance of 2.5 times the flow depth from the bank, a perception which is largely based on experiments with smooth boundaries. The reported results indicate that a stable pattern of secondary currents over the entire channel width can only be sustained over a fixed horizontal bed if the bed's roughness is sufficient to provide the required transverse oscillations in the turbulent shear stresses. Contrary to laboratory flumes, alluvial river bed always provide sufficient roughness. The required external forcing of this hydrodynamic instability mechanism is provided by the turbulence-generated near-bank secondary currents. The pattern of near-bank secondary currents depends on the inclination and the roughness of the bank. In all configurations, secondary currents result in a reduction of the bed shear stress in the vicinity of the bank and a heterogeneous bank shear stress that reaches a maximum close to the toe of the bank. Moreover, these currents cause transverse variability of 10–15% for the streamwise velocities and 0.2u_*²–0.3u_*?² for the bed shear stress. These variations are insufficient to provide the flow variability required in river restoration projects, but nevertheless must be accounted for in the design of stable channels.     

Numerical modeling of lock-exchange gravity/turbidity currents by a high-order upwinding combined compact difference scheme

This study presents two-dimensional direct numerical simulations for sediment-laden current with higher density propagating forward through a lighter ambient water.The incompressible NavierStokes equations including the buoyancy force for the density difference between the light and heavy fluids are solved by a finite difference scheme based on a structured mesh.The concentration transport equations are used to explore such rich transport phenomena as gravity and turbidity currents.Within the framework of an Upwinding Combined Compact finite Difference(UCCD)scheme,rigorous determination of weighting coefficients underlies the modified equation analysis and the minimization of the numerical modified wavenumber.This sixth-order UCCD scheme is implemented in a four-point grid stencil to approximate advection and diffusion terms in the concentration transport equations and the first-order derivative terms in the Navier-Stokes equations,which can greatly enhance convective stability and increase dispersive accuracy at the same time.The initial discontinuous concentration field is smoothed by solving a newly proposed Heaviside function to prevent numerical instabilities and unreasonable concentration values.A two-step projection method is then applied to obtain the velocity field.The numerical algorithm shows a satisfying ability to capture the generation,development,and dissipation of the Kelvin-Helmholz instabilities and turbulent billows at the interface between the current and the ambient fluid.The simulation results also are compared with the data in published literatures and good agreements are found to prove that the present numerical model can well reproduce the propagation,particle deposition,and mixing processes of lock-exchange gravity and turbidity currents.     

Gravity currents over a rigid and emergent vegetated slope

Lock-exchange experiments are conducted to investigate the effects of emergent vegetation on gravity currents flowing down a slope. Rigid and emergent cylinders are used to represent vegetation such as reeds in aquatic environments. The results show that the head of gravity currents without cylinders forms a semi-elliptic shape, similar to the flat bed case. For the slope-induced gravity, however, the head of gravity currents gradually grows and accelerates. A steeper slope without cylinders causes the evident entrainment and subsequent energy dissipation of gravity currents with ambient fluid. When the cylinder population becomes denser, developments of the head become slower, and the semi-elliptic head can be only seen at the very front end. As the density reaches 6.9%, the interface between the saltwater and freshwater performs an inclined straight line, corresponding to the pervious theoretical derivation. The gravity current event could go through two different routes under the similar toe speeds. The current head without cylinders grows along the downslope path, and the enrollment and mixing with ambient fluid could occur at the edge of the head; while the currents within dense cylinders (∼6.9%) has a nearly constant thickness and a smaller and more streamlined head at the very front end. As the cylinder density increases, the front speed of gravity currents transforms from acceleration to deceleration phases. The experimental results reveal that the head would accelerate over the downslope course if the cylinder density is less than 2%.     

长江中下游浅水湖泊悬浮颗粒物吸收系数测定方法探讨

对两种水体悬浮颗粒物吸收系数测定方法及相关计算进行对比研究.通过长江中下游湖泊典型藻类的实验室培养,利用T方法和T-R方法分别对藻类颗粒物、藻类泥沙混合悬浊液进行吸收系数测定.通过颗粒物光谱吸收系数与叶绿素a之间的相关性关系,对比了两种方法的测量稳定性.通过对不同比例的藻类和无机悬浮颗粒物(ISS)的混合悬浊液进行分析,获得了不同浊度水体悬浮物吸收光谱的变化情况.结果表明,在纯藻或者泥沙含量较少的水体进行颗粒物吸收系数光谱测定时,T方法和T-R方法均可以采用,并且均具有较高的测定精度.然而,在泥沙含量相对较高的浑浊水体,应尽量选取T-R方法进行颗粒物吸收光谱的测定,以提高测定精度.长江中下游浅水湖泊由于底泥易受风浪影响发生再悬浮,因此在颗粒物吸收系数光谱测定中,当水体中ISS含量超过30 mg/L时,应选择T-R方法.     

Sill dynamics and fjord deep water renewal: Idealized modeling study

Fjord exchange circulation and its response to abrupt changes in forcing is examined by means of an idealized modeling experiment. Puget Sound, a fjord-type estuary in western North America (State of Washington), is the main context for this study. Parameters of the idealized model are representative of the entrance sill at Admiralty Inlet and the Main Basin of Puget Sound. Sensitivity to some of the model parameters relevant to a 3D realistic model is discussed. An idealized tidal forcing with fortnightly modulation drives a qualitatively realistic cycle of exchange circulation while the other boundary conditions are kept fixed in time. The cycle is characterized by fortnightly pulses of deep water intrusions with a sharp front at the leading edge and reversed circulation cells below the sill depth developing between the intrusions. This basic state is then perturbed and response of the circulation to abrupt changes in oceanic salinity and river discharge is examined.     

A MATHEMATICAL MODEL OF RESERVOIR SEDIMENTATION

1 INTRODUCTIONBed aggradation and degradation haPpen to be the most imPOrtant aspects of the alluvial processes instreams if the equilibrium conditions among water discharge, sediment flow, and channel shaPe areclistuIbed by natural or man-made factOfs, e.g., the constrUchon of a dam, change in the sediment suPplyrate, or base level changing. Reliable and quanhtative estimation of the bed aggradation or degradation isimPortant in rivertalning engineering and water management projects. …     

Numerical modeling of the propagation and morphological changes of turbidity currents using a cost-saving strategy of solution updating

Existing layer-averaged numerical models for turbidity currents have mostly adopted the global minimum time step (GMiTS) for solution updating, which confines their computational efficiency and limits their attractiveness for field applications. This paper presents a highly efficient layer-averaged numerical model for turbidity currents by implementing the combined approach of the local graded-time-step (LGTS) and the global maximum-time-step (GMaTS). The governing equations are solved for unstructured triangular meshes by the shock-capturing finite volume method along with a set of well-balanced evaluations of the numerical flux and geometrical slope source terms. The quantitative accuracy of the model, given reasonably estimated empirical and model parameters (e.g., bed friction, water entrainment, sediment deposition and erosion coefficients), is demonstrated by comparing the numerical solutions against laboratory data of the current front positions and deposition profiles, as well as field data of the current front positions. The improved computational efficiency is demonstrated by comparing the computational cost of the present model against that of a traditional model that uses a GMiTS. For the present simulated cases, the maximum reduction of the computational cost is approximately 80% (e.g., a simulation that cost 1 h before will only require 12 min with the new model).     

Patterns of bedload entrainment and transport in forested headwater streams of the Columbia Mountains,Canada

We monitor bedload transport and water discharge at six stations in two forested headwater streams of the Columbia Mountains, Canada. The nested monitoring network is designed to examine the effects of channel bed texture, and the influence of alluvial (i.e. step pools and riffle pools) and semialluvial morphologies (i.e. boulder cascades and forced step pools) on bedload entrainment and transport. Results indicate that dynamics of bedload entrainment are influenced by differences in flow resistance attributable to morphology. Scaled fractional analysis shows that in reaches with high form resistance most bedload transport occurs in partial mobility fashion relative to the available bed material, while calibers finer than 16 mm attain full mobility during bankfull flows. Equal mobility transport for a wider range of grain sizes is achieved in reaches exhibiting reduced form resistance. Our findings confirm that the Shields value for mobilization of the median surface grain size depends on channel gradient and relative submergence; however, we also find that these relations vary considerably for cobble and gravel bed channels due to proportionality between dimensionless shear stress and grain size. Exponents of bedload rating curves across sites correlate most with the D_90s of the mobile bed, however, where grain effects are controlled (i.e. along individual streams), differences in form resistance across morphologies exert a primary control on bedload transport dynamics. Application of empirical formulae developed for use in steep alpine channels present variable success in predicting transport rates in forested snowmelt streams. Formulae that explicitly account for reductions in mobile bed area and high morphological resistance associated with woody debris provide the best approximation to observed empirical data. Copyright © 2014 John Wiley & Sons, Ltd.     

Formation of high-grade ignimbrites Part II. A pyroclastic suspension current model with implications also for low-grade ignimbrites

 Analogue experiments in part I led to the conclusion that pyroclastic flows depositing very high-grade ignimbrite move as dilute suspension currents. In the thermo–fluid–dynamical model developed, the degree of cooling of expanded turbulent pyroclastic flows dynamically evolves in response to entrainment of air and mass loss to sedimentation. Initial conditions of the currents are derived from column-collapse modeling for magmas with an initial H₂O content of 1–3 wt.% erupting through circular vents and caldera ring-fissures. The flows spread either longitudinally or radially from source up to a runout distance that increases with higher mass flux but decreases with higher gas content, temperature, bottom slope and coarser initial grain size. Progressive dilution by entrainment and sedimentation causes pyroclastic currents to transform into buoyant ash plumes at the runout distance. The ash plumes reach stratospheric heights and distribute 30–80% of the erupted material as widespread co-ignimbrite ash. Pyroclastic suspension currents with initial mass fluxes of 10⁷-10¹² kg/s can spread for tens of kilometers with only limited cooling, although they move as supercritical, strongly entraining currents for the eruption conditions considered here. With increasing eruption mass flux, cooling during passage through the fountain diminishes while cooling during flow transport increases. The net effect is that eruption temperature exerts the prime control on emplacement temperature. Pyroclastic suspension currents can form welded ignimbrite across their entire extent if eruption temperature is T_o>1.3^.T_mw, the minimum welding temperature. High eruption rates, a large fraction of fine ash, and a ring-fissure vent favor the formation of extensive high-grade ignimbrite. For very hot eruptions producing sticky, partially molten pyroclasts, analysis of particle aggregation systematics shows that factors favoring longer runout also favor more efficient aggregation, which reduces runout. As a result, very high-grade ignimbrites cannot spread more than a few tens of kilometers from their source. In cooler pyroclastic currents, particles do not aggregate, and the sedimentation process may involve re-entrainment of particles, which potentially leads to more extensive cooling and longer runout; such effects, however, are only significant when net erosion of substrate occurs. Model results can be employed to estimate mass flux and duration of ignimbrite eruptions from measured ignimbrite masses and aspect ratios. The model also provides an alternative explanation of the observed decrease in H/Lratios with ignimbrite mass. Received: 10 May 1998 / Accepted: 21 October 1998     

Assessing entrainment of bed material in a debris‐flow event: a theoretical approach incorporating Monte Carlo method

Entrainment of underlying bed sediment by a debris flow can significantly increase the debris‐flow magnitude. To study this phenomenon, a theoretical approach to assessing bed‐sediment entrainment is presented. The approach is based on a static approximation that bed‐sediment entrainment occurs when the shearing stress of the flow is sufficiently high to overcome the basal resistance of the bed sediment. In order to delineate erodible zones in a channel, we analyze the critical condition of this static equilibrium model, and subsequently propose a new concept of a critical line to detect the entrainment reaches in a channel. Considering the spatial and temporal uncertainties of the input parameter, the approach is further incorporated within a Monte Carlo method, and the distribution of entrainment zones and post‐entrainment volumes can be analyzed. This approach is illustrated by back‐analysis of the 2010 Yohutagawa debris‐flow event, Japan. Results from 10 000 trials of Monte Carlo simulation are compared with the in situ surveys. It is shown that the present approach can be satisfactorily used to delineate erodible zones and estimate possible entrainment volume of the event. Discussion regarding the sensitivities and limitations of the approach concludes the paper. Copyright © 2015 John Wiley & Sons, Ltd.     

Modeling Monthly Mean Flow in a Poorly Gauged Basin by Fuzzy Logic

The estimation of the monthly mean flow is a critical issue in many water resource development projects. However, in practice the mean flow is not easily determined in ungauged and poorly gauged basins. Therefore, in the literature, various flow estimation methods have been developed recently for mountainous regions which are generally ungauged or poorly gauged basins. In this study a fuzzy logic model based on the Mamdani approach was developed to estimate the flow for poorly gauged mountainous basins. This model was applied to the Solakli Basin which is located in the Eastern Black Sea Region of Turkey. Limited rainfall and flow data are available for this basin. In addition to these variables, the stream and time coefficients were introduced and used as variables for modeling. The data was divided into training and testing phases. The model results were compared with the measured data. The comparison depends on seven statistical characteristics, four different error modes and the contour map method. It was observed that the fuzzy model developed in this study yielded reliable results.     

Two-dimensional density currents in a confined basin

We present new experimental results on the mechanisms through which steady two-dimensional density currents lead to the formation of a stratification in a closed basin. A motivation for this work is to test the underlying assumptions in a diffusive “filling box” model that describes the oceanic thermohaline circulation (Hughes, G.O. and Griffiths, R.W., A simple convective model of the global overturning circulation, including effects of entrainment into sinking regions, Ocean Modeling, 2005, submitted.). In particular, they hypothesized that a non-uniform upwelling velocity is due to weak along-slope entrainment in density currents associated with a large horizontal entrainment ratio of E _eq ?～?0.1. We experimentally measure the relationship between the along-slope entrainment ratio, E, of a density current to the horizontal entrainment ratio, E _eq, of an equivalent vertical plume. The along-slope entrainment ratios show the same quantitative decrease with slope as observed by Ellison and Turner (, 6, 423–448.), whereas the horizontal entrainment ratio E _eq appears to asymptote to a value of E eq?=?0.08 at low slopes. Using the measured values of E _eq we show that two-dimensional density currents drive circulations that are in good agreement with the two-dimensional filling box model of Baines and Turner (Baines, W.D. and Turner, J.S., Turbulent buoyant convection from a source in a confined region, J. Fluid. Mech., 1969, 37, 51–80.). We find that the vertical velocities of density fronts collapse onto their theoretical prediction that U =-2^?2/3 B ^1/3 E _eq ^2/3 (H/R) ζ, where U is the velocity, H the depth, B the buoyancy flux, R the basin width, E _eq the horizontal entrainment ratio and?ζ?= z/H the dimensionless depth. The density profiles are well fitted with?Δ?= 2^?1/3 B ^2/3 E _eq ^?2/3 H _-1 [ln(ζ )?+?τ ], where?τ?is the dimensionless time. Finally, we provide a simple example of a diffusive filling box model, where we show how the density stratification of the deep Caribbean waters (below 1850?m depth) can be described by a balance between a steady two-dimensional entraining density current and vertical diffusion in a triangular basin.     

Flow characteristics of a gravity current induced by differential cooling in a small lake

Extensive field data of a density current driven by the differences in the cooling rates between the two lobes of a small lake are presented. These data illustrate the fact that this gravitational motion affects the whole system and dominates its lakewide circulation. Moreover, data are used to evaluate the entrainment into the density current and also to discuss the regime of the flow from a scaling analysis of the horizontal momentum equations. This analysis is specifically applied to the central area of the lake, where the density current flows from one lobe to the other. Results of this study show that in the longitudinal direction there is a balance between the pressure gradient and the turbulent viscous term. Further, it is found that geostrophy correctly describes the transversal pattern of the passage of the gravity current from the northern lobe to the southern, where it flows confined to the western shoreline.