Determination of critical shear stress of non‐cohesive soils using submerged jet test and turbulent kinetic energy

Laboratory tests using Jet Erosion Testing (JET) apparatus, impinging normally on a horizontal boundary, were conducted to determine the critical shear stress (τ_c) of non‐cohesive soil samples. A three‐dimensional (3D) SonTek/YSI 16 MHz Micro‐Acoustic Doppler Velocimeter (MicroADV) was used to measure turbulent kinetic energy (TKE) at a radial limit of entrainment in the wall jet zone and the measurements were used to calculate τ_c of the samples. The results showed that TKE increases exponentially with increasing particle size. The τ_c from this study were comparable (R² = 0.8) to the theoretical τ_c from Shields diagram after bed roughness scale ratio (D/k_s), due to the non‐uniform bed conditions, was accounted for. This study demonstrated that JET and TKE can be used to determine τ_c of non‐cohesive soils. The use of JET and TKE was found to be faster and easier when compared to the conventional approach of using flumes. A relationship of TKE at the onset of incipient motion (TKE_c) and samples’ D₅₀ developed in this study can be used to predict τ_c of non‐cohesive soils under similar non‐uniform conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Implementation of a second‐order paraxial boundary condition for a water‐saturated layered half‐space in plane strain

A half‐space finite element and a transmitting boundary are developed for a water‐saturated layered half‐space using a paraxial boundary condition. The exact dynamic stiffness of a half‐space in plane strain is derived and a second‐order paraxial approximation of the stiffness is obtained. A half‐space finite element and a transmitting boundary are then formulated. The development is verified by comparison of the dynamic stiffness of impermeable and permeable rigid strip foundations with other published results. The advantage of using the paraxial boundary condition in comparison with the rigid boundary condition is examined. It is shown that the paraxial boundary condition offers significant gain and the resulting half‐space finite element and transmitting boundary can represent the effects of a water‐saturated layered half‐space with good accuracy and efficiency. In addition, the numerical method described herein maintains the strengths and advantages of the finite element method and can be easily applied to demanding problems of soil–structure interaction in a water‐saturated layered half‐space. Copyright © 2010 John Wiley & Sons, Ltd.     

Phase‐shifts in shear stress as an explanation for the maintenance of pool–riffle sequences

The stability of the pool–rif?e sequence is one of the most fundamental features of alluvial streams. For several decades, the process of velocity, or shear stress, reversal has been proposed as an explanation for an increase in the amplitude of pool–rif?e sequence bars during high ?ows, offsetting gradual scour of rif?es and deposition in pools during low ?ows. Despite several attempts, reversal has rarely been recorded in ?eld measurements. We propose that, instead of being reversed, maxima and minima in shear stress are phase‐shifted with respect to the pool–rif?e sequence bedform pro?le, so that maximum shear stress occurs upstream of rif?e crests at high ?ow, and downstream at low ?ow. Such phase‐shifts produce gradients of shear stress that explain rif?e deposition, and pool scour, at high ?ow, in accord with sediment continuity. The proposal is supported by results of a one‐dimensional hydraulic model applied to the surveyed bathymetry of a pool–rif?e sequence in a straight reach of a gravel‐bed river. In the sequence studied, the upstream phase‐shift in shear stress at high ?ow was associated with variations in channel width, with width minima occurring upstream of rif?e crests, approximately coincident with shear stress maxima, and width maxima occurring downstream of rif?e crests. Assuming that the width variation is itself the result of ?ow de?ection by rif?e crests at low ?ow, and associated bank‐toe scour downstream, low and high ?ow can be seen to have complementary roles in maintaining alluvial pool–rif?e sequences. Copyright © 2004 John Wiley & Sons, Ltd.     

Experimental study of extreme shear stress for shallow flow under simulated rainfall

The raindrop impact and overland flow are two major factors causing soil detachment and particle transportation. In this study, the turbulent characteristics of the shallow rain‐impacted water flow were investigated using a 2‐D fibre‐optic laser Doppler velocimetry (FLDV) and an artificial rainfall simulator. The fluctuating turbulent shear stress was computed using digital data processing techniques. The experimental data showed that the Reynolds shear stress follows a probability distribution with heavy tails. The tail probability increases with an increase of rainfall intensity or raindrop diameter, and it decreases with an increase of Reynolds number. A modified empirical equation was derived using both the raindrop diameter and rainfall intensity as independent variables to provide a better prediction of the Darcy‐Weisbach friction coefficient f under rainfall conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

Sensitivity of interfacial hydraulics to the microtopographic roughness of water‐lain gravels

Flow within the interfacial layer of gravel‐bed rivers is poorly understood, but this zone is important because the hydraulics here transport sediment, generate flow structures and interact with benthic organisms. We hypothesized that different gravel‐bed microtopographies generate measurable differences in hydraulic characteristics within the interfacial layer. This was tested using a high density of spatially and vertically distributed, velocity time series measured in the interfacial layers above three surfaces of contrasting microtopography. These surfaces had natural water‐worked textures, captured in the field using a casting procedure. Analysis was repeated for three discharges, with Reynolds numbers between 165000 and 287000, to evaluate whether discharge affected the impact of microtopography on interfacial flows. Relative submergence varied over a small range (3.5 to 8.1) characteristic of upland gravel‐bed rivers. Between‐surface differences in the median and variance of several time‐averaged and turbulent flow parameters were tested using non‐parametric statistics. Across all discharges, microtopographic differences did not affect spatially averaged (median) values of streamwise velocity, but were associated with significant differences in its spatial variance, and did affect spatially averaged (median) turbulent kinetic energy. Sweep and ejection events dominated the interfacial region above all surfaces at all flows, but there was a microtopographic effect, with Q2 and Q4 events less dominant and structures less persistent above the surface with the widest relief distribution, especially at the highest Reynolds number flow. Results are broadly consistent with earlier work, although this analysis is unique because of the focus on interfacial hydraulics, spatially averaged ‘patch scale’ metrics and a statistical approach to data analysis. An important implication is that observable differences in microtopography do not necessarily produce differences in interfacial hydraulics. An important observation is that appropriate roughness parameterizations for gravel‐bed rivers remain elusive, partly because the relative contributions to flow resistance of different aspects of bed microtopography are poorly constrained. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Removal of water‐surface reflection effects with a temporal minimum filter for UAV‐based shallow‐water photogrammetry

The recent development of structure‐from‐motion (SfM) and multi‐view stereo (MVS) photogrammetry techniques has enabled semi‐automatic high‐resolution bathymetry using aerial images taken by consumer‐grade digital cameras mounted on unmanned aerial vehicles (UAVs). However, the applicability of these techniques is sometimes limited by sun and sky reflections at the water surface, which render the point‐cloud density and accuracy insufficient. In this research, we present a new imaging technique to suppress the effect of these water‐surface reflections. In this technique, we order a drone to take a short video instead of a still picture at each waypoint. We then apply a temporal minimum filter to the video. This filter extracts the smallest RGB values in all the video frames for each pixel, and composes an image with greatly reduced reflection effects. To assess the performance of this technique, we applied it at three small shallow‐water sites. Specifically, we evaluated the effect of the technique on the point cloud density and the accuracy and precision of the photogrammetry. The results showed that the proposed technique achieved a far denser point cloud than the case in which a randomly chosen frame was used for each waypoint, and also showed better overall accuracy and precision in estimating water‐bottom elevation. The effectiveness of this new technique should depend on the surface wave state and sky radiance distribution, and this dependence, as well as the applicability to large areas, should be investigated in future research. Copyright © 2018 John Wiley & Sons, Ltd.     

The causes of bedload pulses in a gravel channel: the implications of bedload grain‐size distributions

Unsteady bedload transport was measured in two c. 5 m wide anabranches of a gravel‐bed braided stream draining the Haut Glacier d'Arolla, Switzerland, during the 1998 and 1999 melt seasons. Bedload was directly sampled using 152 mm square Helley–Smith type samplers deployed from a portable measuring bridge, and independent transport rate estimates for the coarser size fractions were obtained from the dispersion of magnetically tagged tracer pebbles. Bedload transport time series show pulsing behaviour under both marginal (1998) and partial (1999) transport regimes. There are generally weak correlations between transport rates and shear stresses determined from velocity data recorded at the measuring bridge. Characteristic parameters of the bedload grain‐size distributions (D₅₀, D₈₄) are weakly correlated with transport rates. Analysis of full bedload grain‐size distributions reveals greater structure, with a tendency for transport to become less size selective at higher transport rates. The bedload time series show autoregressive behaviour but are dif?cult to distinguish by this method. State–space plots, and associated measures of time‐series separation, reveal the structure of the time series more clearly. The measured pulses have distinctly different time‐series characteristics from those modelled using a one‐dimensional sediment routing model in which bed shear stress and grain size are varied randomly. These results suggest a mechanism of pulse generation based on irregular low‐amplitude bedforms, that may be generated in‐channel or may represent the advection of material supplied by bank erosion events. Copyright © 2003 John Wiley & Sons, Ltd.     

Cyclic tests of steel plate shear walls using box‐shape vertical boundary elements with or without infill concrete

The steel plate shear wall (SPSW) is an effective lateral force resisting system in which unstiffened steel infill plates are connected to the horizontal and vertical boundary elements (VBEs) on all sides of the plates. The boundary elements must be designed to resist the tension field force of the infill panels. When the VBEs are made from a steel box section, the flange of each box VBE connected with the infill panels can be pulled out‐of‐plane by significant panel forces, called pull‐out action. This study investigates capacity design methods for box VBEs in SPSWs. Simplified fixed beam and portal frame models aim to estimate the pull‐out responses of the flange of the box sections with and without infill concrete, respectively. In this study, cyclic tests of three full‐scale two‐story SPSWs using box VBEs with or without the infill concrete are conducted. Inelastic pushover analyses of the finite element models are conducted. The tests and analytical results confirm that the proposed design methods, which aimed to prevent the full yield of the flange under the pull‐out action, are applicable. Furthermore, the test and analytical results suggest that the initial yielding of the flange of box VBEs under the collective effects of the pull‐out action on the flange, the gravity load, and the sway action on the SPSW represents a local yielding. A strict prohibition of the initial yielding on the flange under the aforementioned collective effects is not recommended for pursuing a cost‐effective design. Copyright © 2017 John Wiley & Sons, Ltd.     

Coupled shear‐bending formulation for seismic analysis of stacked wood shear wall systems

A summary of the development of a new coupled shear‐bending model for analysis of stacked wood shear walls and multi‐story wood‐frame buildings is presented in this paper. The model focuses on dynamic response of light‐frame wood structures under seismic excitation. The formulation is intended to provide a more versatile option than present pure shear models in that the new model is capable of accurately capturing the overall lateral response of each story diaphragm and separates the inter‐story shear deformation and the deformation associated with rotation of the diaphragm due to rod elongation, which is an analogue to the bending deformation in an Euler–Bernoulli beam model. Modeling the coupling of bending and shear deformation is shown to provide more accurate representation of stacked shear wall system behavior than a pure shear model, particularly for the upper stories in the assembly. The formulation is coupled with the newly developed evolutionary parameter hysteretic model for wood shear walls. Existing data from a shake table test of an isolated three‐story wood shear wall were used to verify the accuracy of the model prediction. The numerical results agreed very well with shake table test measurements. The influence of a continuous rod hold‐down system on the dynamic behavior of the three‐story stacked wood shear wall was also successfully simulated. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical investigation of the resistance of precast RC pinned beam‐to‐column connections under shear loading

In precast technology, the effective design and construction are related to the behaviour of the connections between the structural members in order to cater for all service, environmental and earthquake load conditions. Therefore, the design and detailing of the connections should be undertaken consistently and with awareness of the desired structural response. In the research presented herein, an analytical expression is proposed for the prediction of the resistance of precast pinned connections under shear monotonic and cyclic loading. The proposed formula addresses the case where the failure of the connection occurs with simultaneous flexural failure of the dowel and compression failure of the concrete around the dowel, expected to occur either when (i) adequate concrete cover of the dowels is provided (d > 6 D) or (ii) adequate confining reinforcement (as defined in the article) is foreseen around the dowels in the case of small concrete covers (d < 6 D). The expression is calibrated against available experimental data and numerical results derived from a nonlinear numerical investigation. Emphasis is given to identifying the effect of several parameters on the horizontal shear resistance of the connection such as: the number and diameter of the dowels; the strength of materials (concrete, grout, steel); the concrete cover of the dowels; the thickness of the elastomeric pad; the type of shear loading (monotonic or cyclic); the pre‐existing axial stress in the dowels; and the rotation of the joint. In addition, recommendations for the design of precast pinned beam‐to‐column connections are given, especially when the connections are utilised in earthquake resistant structures. Copyright © 2017 John Wiley & Sons, Ltd.     

Use of Large-Eddy Simulation for the bed shear stress estimation over a dune

Environmental flows are generally characterized by complex bed morphology and high current speeds. Such configurations favor the formation of vortex structures that strongly affect hydrody-namics and sediment transport. Large-Eddy Simulation (LES) enables investigation of the dynam-ics of the largest turbulence scales and, thanks to enhanced calculation resources, has now become applicable for simulating environmental flows. In this paper, a LES approach is developed in a CFD code (TELEMAC-3D), which was originally developed to simulate free surface flows using RANS methods. The present developments involve implementing subgrid models, boundary con-ditions and numerical schemes suitable for LES. The LES version of TELEMAC-3D was validated by comparing results on the model with experimental data for flow past a cylinder. Then, the model was applied to a test case representing flow over dunes. After validating the hydrodynamics, the model was used to assess the bottom shear stress, using both a RANS and a LES approach. Com-parison highlighted the potential contribution of LES to investigating the hydrodynamic forces acting on the bottom.     

Direct finite element method for nonlinear analysis of semi‐unbounded dam–water–foundation rock systems

A direct finite element method is presented for nonlinear earthquake analysis of interacting dam–water–foundation rock systems. The analysis procedure applies viscous damper absorbing boundaries to truncate the semi‐unbounded fluid and foundation‐rock domains and specifies at these boundaries effective earthquake forces determined from the design ground motion defined at a control point on the free surface. The analysis procedure is validated numerically by computing the frequency response functions and transient response of an idealized dam–water–foundation rock system and comparing with results from the substructure method. Because the analysis procedure is applicable to nonlinear systems, it allows for modeling of concrete cracking, as well as sliding and separation at construction joints, lift joints, and at concrete–rock interfaces. Implementation of the procedure is facilitated by commercial finite element software with nonlinear material models that permit modeling of viscous damper boundaries and specification of effective earthquake forces at these boundaries. Copyright © 2017 John Wiley & Sons, Ltd.     

Time‐harmonic response of non‐homogeneous elastic unbounded domains using the scaled boundary finite‐element method

The scaled boundary finite‐element method is extended to simulate time‐harmonic responses of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. The unbounded domains and the elasticity matrices are transformed to the scaled boundary coordinates. The scaled boundary finite‐element equation in displacement amplitudes are derived directly from the governing equations of elastodynamics. To enforce the radiation condition at infinity, an asymptotic expansion of the dynamic‐stiffness matrix for high frequency is developed. The dynamic‐stiffness matrix at lower frequency is obtained by numerical integration of ordinary differential equations. Only the boundary is discretized yielding a reduction of the spatial dimension by one. No fundamental solution is required. Material anisotropy is modelled without additional efforts. Examples of two‐ and three‐dimensional non‐homogeneous isotropic and transversely isotropic unbounded domains are presented. The results demonstrate the accuracy and simplicity of the scaled boundary finite‐element method. Copyright © 2005 John Wiley & Sons, Ltd.     

Surfactant effects on the water‐stable aggregation of wettable soils from the continental USA

Surfactants may affect soil structure differently depending upon the soil or the quality of rainfall or irrigation water. This study examined whether the water‐stable aggregation of 11 wettable soils was affected by surfactants and the water in which the soils were sieved. The study also examined whether the wettable soils' water drop penetration time (WDPT) was affected by surfactants, water drop quality, and elapsed time since the surfactants were applied. Two nonionic surfactants and a surfactant‐free water control were sprayed (by misting) upon air‐dry soil, then WDPT was measured 1 and 72 h thereafter. Subsequently, this treated soil was slowly wetted with an aerosol to its water content at a matric potential of ?3 kPa, then immediately sieved for 600 s in water that contained either appreciable or few electrolytes. Water‐stable aggregation, quantified as mean weight diameter (MWD), varied widely among soils, ranging from 0.10 to 1.36 mm. The MWDs were affected (at p = 0.06) by surfactant treatments, depending upon the soil but not sieving water quality. Surfactants affected the MWD of an Adkins loamy sand and Feltham sand, two of the three coarsest‐textured soils. Although WDPTs never exceeded 5 s, depending upon the soil WDPTs were affected by surfactant treatments but not by water drop quality. After surfactant application, WDPTs generally decreased with time for three soils but increased with time for one soil. Findings suggested that surfactants interacted (1) with clay mineralogy to affect MWD and (2) with soluble calcium to affect WDPT for certain soils. Surfactant treatments but not water quality affected both MWD and WDPT for some but not all of 11 wettable, US soils. Published 2012. This article is a US Government work and is in the public domain in the USA.     

根据环境应力场方法确定地震危险区的3年实践检验

地震是在剪应力作用下发生的岩石大规模破裂过程。因此，环境剪应力值τ０是控制地震发生的主要因素，我们根据环境应力场确定大震危险区，经过３年的时间检验，结果表明，大约７０％的地震都发生在我们事先确定的危险区内。     

Comparison of web plate numerical models for self‐centering steel plate shear walls

Previous research has shown that self‐centering steel plate shear walls (SC‐SPSWs) are capable of achieving enhanced seismic performance at multiple hazard levels, including recentering following design‐level earthquakes. When modeling SC‐SPSWs numerically, these studies considered an idealized tension‐only steel plate shear wall (SPSW) web plate behavior. Research has shown that web plate behavior is more complex than predicted by the idealized model, and web plates can provide more strength, stiffness, and energy dissipation than predicted by the idealized model. The idealized model of web plate behavior is used widely in SPSW numerical models where the moment‐resisting boundary frame provides supplemental hysteretic damping and stiffness; however, in SC‐SPSWs, where the post‐tensioned boundary frame is designed to remain elastic during an earthquake, accounting for the more complex web plate behavior can have a significant impact on seismic performance estimates from numerical simulation. This paper presents different methods for modeling SC‐SPSWs. Responses from these models are compared with experimental results. A simple modification of the tension‐only model, referred to as the tension‐compression strip model, is shown to provide a reasonable approximation of SC‐SPSW behavior. Results from nonlinear response history analyses of SC‐SPSWs with the tension‐only and tension‐compression web plate models are compared to assess how the approximation of web plate behavior affects SC‐SPSW seismic performance. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluating the dual‐boundary forcing concept in subsurface–land surface interactions of the hydrological cycle

Subsurface and land surface processes (e.g. groundwater flow, evapotranspiration) of the hydrological cycle are connected via complex feedback mechanisms, which are difficult to analyze and quantify. In this study, the dual‐boundary forcing concept that reveals space–time coherence between groundwater dynamics and land surface processes is evaluated. The underlying hypothesis is that a simplified representation of groundwater dynamics may alter the variability of land surface processes, which may eventually affect the prognostic capability of a numerical model. A coupled subsurface–land surface model ParFlow.CLM is applied over the Rur catchment, Germany, and the mass and energy fluxes of the coupled water and energy cycles are simulated over three consecutive years considering three different lower boundary conditions (dynamic, constant, and free‐drainage) based on groundwater dynamics to substantiate the aforementioned hypothesis. Continuous wavelet transform technique is applied to analyze scale‐dependent variability of the simulated mass and energy fluxes. The results show clear differences in temporal variability of latent heat flux simulated by the model configurations with different lower boundary conditions at monthly to multi‐month time scales (~32–91 days) especially under soil moisture limited conditions. The results also suggest that temporal variability of latent heat flux is affected at even smaller time scales (~1–3 days) if a simple gravity drainage lower boundary condition is considered in the coupled model. This study demonstrates the importance of a physically consistent representation of groundwater dynamics in a numerical model, which may be important to consider in local weather prediction models and water resources assessments, e.g. drought prediction. Copyright © 2015 John Wiley & Sons, Ltd.     

Computational fluid dynamics modelling of boundary roughness in gravel‐bed rivers: an investigation of the effects of random variability in bed elevation

Results from a series of numerical simulations of two‐dimensional open‐channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open‐channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra‐grid‐scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h = 0·1–0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra‐grid‐scale roughness also allows a reduction in the height of the near‐bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near‐bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel‐bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh‐independent solutions. Copyright © 2001 John Wiley & Sons, Ltd.