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.     

Location and height of intertidal bars on macrotidal ridge and runnel beaches

Previous studies devoted to the morphology and hydrodynamics of ridge and runnel beaches highlight characteristics that deviate from those initially postulated by King and Williams (Geographical Journal, 1949, vol. 113, 70–85) and King (Beaches and Coasts, 1972, Edward Arnold). Disagreements on the morphodynamics of these macrotidal beaches include the position of the ridges relative to the mean neap and spring tide levels, the variation in the height of the ridges across the intertidal profile and, most importantly, whether the ridges are formed by swash or surf zone processes. The morphological characteristics of ridge and runnel beaches from three locations with varying wave, tidal and geomorphic settings were investigated to address these disagreements. Beach profiles from each site were analysed together with water‐level data collected from neighbouring ports. It was found that the ridges occur over the entire intertidal zone. On one site (north Lincolnshire, east England), the ridges are uniformly distributed over the intertidal beach, whereas on the two other sites (Blackpool beach, northwest England, and Leffrinckoucke beach, north France) there is some indication that the ridges appear to occur at preferential locations. Most significantly, the locations of the ridge crests were found to be unrelated to the positions on the intertidal profile where the water level is stationary for the longest time. It was further found that the highest ridges generally occur just above mid‐tide level where tidal non‐stationarity is greatest. These findings argue against the hypothesis that the ridges are formed by swash processes acting at stationary tide levels. It is tentatively suggested that the ridges are the result of a combination of swash and surf zone processes acting across the intertidal zone. Elucidation of the morphodynamic roles of these two types of processes, and other processes such as strong current flows in the runnels, requires further comprehensive field measurements complemented by numerical modelling. Copyright © 2001 John Wiley & Sons, Ltd.     

The formation of large mudwaves by turbidity currents on the levees of the Toyama deep-sea channel, Japan Sea

The development of mudwaves on the levees of the modern Toyama deep‐sea channel has been studied using gravity core samples combined with 3·5‐kHz echosounder data and airgun seismic reflection profiles. The mudwaves have developed on the overbank flanks of a clockwise bend of the channel in the Yamato Basin, Japan Sea, and the mudwave field covers an area of 4000 km². Mudwave lengths range from 0·2 to 3·6 km and heights vary from 2 to 44 m, and the pattern of mudwave aggradation indicates an upslope migration direction. Sediment cores show that the mudwaves consist of an alternation of fine‐grained turbidites and hemipelagites whereas contourites are absent. Core samples demonstrate that the sedimentation rate ranged from 10 to 14 cm ka^?1 on the lee sides to 17–40 cm ka^?1 on the stoss sides. A layer‐by‐layer correlation of the deposits across the mudwaves shows that the individual turbidite beds are up to 20 times thicker on the stoss side than on the lee side, whereas hemipelagite thicknesses are uniform. This differential accretion of turbidites is thought to have resulted in the pattern of upcurrent climbing mudwave crests, which supports the notion that the mudwaves have been formed by spillover turbidity currents. The mudwaves are interpreted to have been instigated by pre‐existing large sand dunes that are up to 30 m thick and were created by high‐velocity (10°ms^?1), thick (c. 500 m) turbidity currents spilling over the channel banks at the time of the maximum uplift of the Northern Japan Alps during the latest Pliocene to Early Pleistocene. Draping of the dunes by the subsequent, lower‐velocity (10^?1ms^?1), mud‐laden turbidity currents is thought to have resulted in the formation of the accretionary mudwaves and the pattern of upflow climbing. The dune stoss slopes are argued to have acted as obstacles to the flow, causing localized loss of flow strength and leading to differential draping by the muddy turbidites, with greater accretion occurring on the stoss side than on the lee slope. The two overbank flanks of the clockwise channel bend show some interesting differences in mudwave development. The mudwaves have a mean height of 9·8 m on the outer‐bank levee and 6·2 m on the inner bank. The turbidites accreted on the stoss sides of the mudwaves are 4–6 times thicker on the outer‐bank levee than their counterparts on the inner‐bank levee. These differences are attributed to the greater flow volume (thickness) and sediment flux of the outer‐bank spillover flow due to the more intense stripping of the turbidity currents at the outer bank of the channel bend. Differential development of mudwave fields may therefore be a useful indicator in the reconstruction of deep‐sea channels and their flow hydraulics.     

ANALYSIS ON MOUTH BAR EVOLUTION OF THE YELLOW RIVER ESTUARY

ＡＮＡＬＹＳＩＳＯＮＭＯＵＴＨＢＡＲＥＶＯＬＵＴＩＯＮＯＦＴＨＥＹＥＬＬＯＷＲＩＶＥＲＥＳＴＵＡＲＹＪＩＺｕｗｅｎ；ＨＵＣｈｕｎｈｏｎｇ（Ｅｎｇｉｎｅｅｒ，ＤｅｐａｒｔｍｅｎｔｏｆＳｅｄｉｍｅｎｔａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇ，ＣｈｉｎａＩｎｓｔｉ...     

变动流域滞时方程

流域滞时是流域汇流计算和预报中的一个重要因素。本文尝试将流域滞时表示成流域重心长度Lca(流域重心到流域出口断面的河道长度)与汇流速度v之比的函数,从而使流域滞时随汇流速度的不同而不同,而不再是一个常数。长度Lca反映了流域尺度和几何形态对滞时的影响,而汇流速度则反映流域动力条件对滞时的作用。5个小流域的实测资料检验表明结果较为满意。     

泥沙起动条件与卵石浅滩航道整治线宽度的确定

在简述泥沙起动条件表达式的基础上,本文分析卵石浅滩整治与泥沙起动条件的关系。根据起动拖曳力和起动流速的概念,分别导出了卵石浅滩航道整治宽度计算公式,并用北江整治河段的实测资料对诸公式作了比较。结果表明,式(16)的计算值与实际采用值有较好的一致性。     

重庆大学危岩锚固工程简介

重庆大学18层高层建筑物拟建在嘉陵江高陡斜坡脚下,斜坡上部陡崖为侏罗系上沙溪庙组巨厚层砂岩,岩石原生及岸剪裂隙发育,危岩直接威胁着建筑物。川地208队采用锚杆加固及涂抹砂浆水泥护坡综合防治工程,效果明显。     

Hydraulic Effects of Shales in Fluvial-Deltaic Deposits: Ground-Penetrating Radar, Outcrop Observations, Geostatistics, and Three-Dimensional Flow Modeling for the Ferron Sandstone, Utah

Ground-penetrating radar (GPR) surveys, outcrop measurements, and cores provide a high-resolution 3D geologic model to investigate the hydraulic effects of shales in marine-influenced lower delta-plain distributary channel deposits within the Cretaceous-age Ferron Sandstone at Corbula Gulch in central Utah, USA. Shale statistics are computed from outcrop observations. Although slight anisotropy was observed in mean length and variogram ranges parallel and perpendicular to pale of low , the anisotropy is not statistically significant and the estimated mean length is 5.4 m. Truncated Gaussian simulation was used to create maps of shales that are placed on variably dipping stratigraphic surfaces interpreted from high-resolution 3D GPR surveys, outcrop interpretations, and boreholes. Sandstone permeability is estimated from radar responses calibrated to permeability measurements from core samples. Experimentally designed flow simulations examine the effects of variogram range, shale coverage fraction, and trends in shale coverage on predicted upscaled permeability, breakthrough time, and sweep efficiency. Approximately 1500 flow simulations examine three different geologic models, flow in the 3 coordinate directions, 16 geostatistical parameter combinations, and 10 realizations for each model. ANOVA and response models computed from the flow simulations demonstrate that shales decrease sweep, recovery, and permeability, especially in the vertical direction. The effect on horizontal flow is smaller. Flow predictions for ideal tracer displacements at Corbula Gulch are sensitive to shale-coverage fraction, but are relatively insensitive to twofold variations in variogram range or to vertical trends in shale coverage. Although the hydraulic effects of shale are statistically significant, the changes in flow responses rarely exceed 20%. As a result, it may be reasonable to use simple models when incorporating analogous shales into models of reservoirs or aquifers.     

现代冲积河流的河型空间转化模式探讨

论述了现代冲积河流系统中常见河型的自然分布及转化特征,归纳了河型从上游向下游发生空间转化的 5种主要模式。这些河型转化模式发生的相对频率是不同的,其中从辫状河流向曲流河流转化的发生频率最大。文中通过具体实例总结了这些转化模式分别发生时的主要影响因素。     

Spatial and temporal variations in bank erosion on sand‐bed streams in the seasonally wet tropics of northern Australia

Bank erosion rates and processes across a range of spatial scales are poorly understood in most environments, especially in the seasonally wet tropics of northern Australia where sediment yields are among global minima. A total of 177 erosion pins was installed at 45 sites on four sand‐bed streams (Tributaries North and Central, East Tributary and Ngarradj) in the Ngarradj catchment in the Alligator Rivers Region. Bank erosion was measured for up to 3·5 years (start of 1998/99 wet season to end of 2001/02 wet season) at three spatial scales, namely a discontinuous gully (0·6 km²) that was initiated by erosion of a grass swale between 1975 and 1981, a small continuous channel (2·5 km²) on an alluvial fan that was formed by incision of a formerly discontinuous channel between 1964 and 1978, and three medium‐sized, continuous channels (8·5–43·6 km²) with riparian vegetation. The bank erosion measurements during a period of average to above‐average rainfall established that substantial bank erosion occurred during the wet season on the two smaller channels by rapid lateral migration (Tributary Central) and by erosion of gully sidewalls due to a combination of within‐gully flows and overland flow plunging over the sidewalls (Tributary North). Minor bank erosion also occurred during the dry season by faunal activity, by desiccation and loss of cohesion of the sandy bank sediments and by dry flow processes. The larger channels with riparian vegetation (East Tributary and Ngarradj) did not generate significant amounts of sediment by bank erosion. Deposition (i.e. negative pin values) was locally significant at all scales. Bank profile form and channel planform exert a strong control on erosion rates during the wet season but not during the dry season. Copyright © 2006 Commonwealth Government of Australia.