High spatio-temporal resolution measurements of cohesive sediment erosion

In this study, we present a novel approach to measure fundamental processes of cohesive sediment erosion. The experimental setup consists of a laboratory erosion flume (SETEG) and a photogrammetric method to detect sediment erosion (PHOTOSED). Detailed data are presented for three erosion experiments, which were conducted with a natural non-cohesive/cohesive sediment mixture at increasing sediment depths (4, 8, 16 cm). In each experiment, the sediment was exposed to a set of incrementally increasing shear stresses and the erosion was measured dynamically, pixel-based, and approximate to the process scale given the resolution of PHOTOSED. This enables us to distinguish between (i) individual emerging erosion spots caused by surface erosion and (ii) large holes torn open by detached aggregate chunks. Moreover, interrelated processes were observed, such as (iii) propagation of the erosion in the longitudinal and lateral direction leading to merging of disconnected erosion areas and (iv) progressive vertical erosion of already affected areas. By complementing the (bulk) erosion volume profiles with additional quantitative variables, which contain spatial information (erosion area, specific deepening, number of disconnected erosion areas), conclusions on the erosion behaviour (and the dominant processes) can be drawn without requiring qualitative information (such as visual observations). In addition, we provide figures indicating the spatio-temporal erosion variability and the (bulk) erosion rates for selected time periods. We evaluate the variability by statistical quantities and show that significant erosion is mainly confined to only a few events during temporal progression, but then considerably exceeds the time-averaged median of the erosion (factors between 7.0 and 16.0). Further, we point to uncertainties in using (bulk) erosion rates to assess cohesive sediment erosion and particularly the underlying processes. As a whole, the results emphasise the need to measure cohesive sediment erosion with high spatio-temporal resolution to obtain reliable and robust information. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Lattice Boltzmann simulations of the transient shallow water flows

A two-dimensional lattice Boltzmann model (LBM) is presented for transient shallow water flows. The model is based on the shallow water equations coupled with the large eddy simulation model. In order to obtain accurate results efficiently, a multi-block lattice scheme is applied at the area where a local finer grid is needed for strong change in physical variables. The model is verified by applying to five cases with transient processes: (a) a tidal wave over steps; (b) a perturbation over a submerged hump; (c) partial dam break flow; (d) circular dam break flow; (e) interaction between a dam break surge and four square cylinders. The objectives of this study are to validate the two-dimensional LBM in transient flow simulation and provide the detailed transient processes in shallow water flows.     

准地转大洋风生环流的格子Boltzmann数值模拟

建立了求解准地转相当正压涡度方程的格子Boltzmann (LB)模型. 该模型将准地转相当正压涡度方程作为一个平流-扩散-反应方程来加以处理,在整体二阶精度下,通过Chapman_Enskog多尺度分析法,可将格子Boltzmann方程还原到相当正压涡度方程. 在不同Reynolds数、不同边界条件以及不同风应力驱动下的数值解表明,该模型正确反映了风生环流的基本结构和不同边界的耗散特征,并得到风生环流的多平衡态解等非线性特征. 此外,不同Rossby变形半径下的实验证明,小Rossby变形半径更容易激发环流的非线性模态. 通过与同等类型有限差方案的比较,表明本文的LB模型具有稳定性好、精度高等优点.     

Lattice Boltzmann Models for Flow and Transport in Saturated Karst

Flow and transport simulation in karst aquifers remains a significant challenge for the ground water modeling community. Darcy's law–based models cannot simulate the inertial flows characteristic of many karst aquifers. Eddies in these flows can strongly affect solute transport. The simple two-region conduit/matrix paradigm is inadequate for many purposes because it considers only a capacitance rather than a physical domain. Relatively new lattice Boltzmann methods (LBMs) are capable of solving inertial flows and associated solute transport in geometrically complex domains involving karst conduits and heterogeneous matrix rock. LBMs for flow and transport in heterogeneous porous media, which are needed to make the models applicable to large-scale problems, are still under development. Here we explore aspects of these future LBMs, present simple examples illustrating some of the processes that can be simulated, and compare the results with available analytical solutions. Simulations are contrived to mimic simple capacitance-based two-region models involving conduit (mobile) and matrix (immobile) regions and are compared against the analytical solution. There is a high correlation between LBM simulations and the analytical solution for two different mobile region fractions. In more realistic conduit/matrix simulation, the breakthrough curve showed classic features and the two-region model fit slightly better than the advection-dispersion equation (ADE). An LBM-based anisotropic dispersion solver is applied to simulate breakthrough curves from a heterogeneous porous medium, which fit the ADE solution. Finally, breakthrough from a karst-like system consisting of a conduit with inertial regime flow in a heterogeneous aquifer is compared with the advection-dispersion and two-region analytical solutions.     

Using ADV backscatter strength for measuring suspended cohesive sediment concentration

Laboratory experiments were conducted at two institutes to reveal the relationship between acoustic backscatter strength and suspended sediment concentration (SSC). In total, three acoustic Doppler velocimeters (ADVs) with different frequencies (5, 10 and 16 MHz) were tested. Two different commercial clays and one natural sediment from Clay Bank site in the York River were checked for acoustic responses. The SSCs of selected sediments were artificially changed between a selected low and a high value in tap or de-ion water. Each ADV showed quite different backscatter responses depending on the sediment type and SSC. Not all devices had a good linear relationship between backscatter strength and SSC. Within a limited range of SSC, however, the backscatter strength can be well correlated with the SSC. Compared with optical backscattering sensor (OBS), the fluctuation of ADV backscatter signals was too noisy to be directly converted to the instantaneous changes of SSC due to high amplification ratio and small sampling volume. For the more accurate signal conversion for finding the fluctuation of SSC, the ensemble average should be applied to increase the signal-to-noise ratio. There are unexpected responses for the averaged backscatter wave strength: (1) high signals from small particles but low signals from large particles; and (2) two linear segments in calibration slope. These phenomena would be most likely caused by the different gain setting built in ADVs. The different acoustic responses to flocculation might also contribute somewhat if flocs are tightly packed. This study suggests that an ADV could be a useful instrument to estimate suspended cohesive sediment concentration and its fluctuation if the above concerns are clarified.     

The relative contributions of physical and microbiological factors to cohesive sediment stability

The stabilising effects of natural benthic diatom and bacterial assemblages on cohesive sediments were compared with those caused by physico-chemical binding alone. Cohesive sediment beds were reconstructed in 4 annular laboratory miniflumes, using sediment collected at 5–6 m water depth from a local fjord. The sediment was left to stabilise (consolidate) for 1, 2, 5 and 10 days, before being fully resuspended in a series of erosion experiments. The flumes were aerated and subjected to different light/dark conditions; antibiotics were used to isolate diatom from bacteria effects. During consolidation, a constant current velocity was maintained, at a speed well below erosion threshold.     

Significance of particle interaction to the modelling of cohesive sediment transport in rivers

Aggregation processes of fine sediments have rarely been integrated in numerical simulations of cohesive sediment transport in riverine systems. These processes, however, can significantly alter the hydrodynamic characteristics of suspended particulate matter (SPM), modifying the particle settling velocity, which is one of the most important parameters in modelling suspended sediment dynamics. The present paper presents data from field measurements and an approach to integrate particle aggregation in a hydrodynamic sediment transport model. The aggregation term used represents the interaction of multiple sediment classes (fractions) with corresponding multiple deposition behaviour. The k–ε–turbulence model was used to calculate the coefficient of vertical turbulent mixing needed for the two‐dimensional vertical‐plane simulations. The model has been applied to transport and deposition of tracer particles and natural SPM in a lake‐outlet lowland river (Spree River, Germany). The results of simulations were evaluated by comparison with field data obtained for two levels of river discharge. Experimental data for both discharge levels showed that under the prevailing uniform hydraulic conditions along the river reach, the settling velocity distribution did not change significantly downstream, whereas the amount of SPM declined. It was also shown that higher flow velocities (higher fluid shear) resulted in higher proportions of fast settling SPM fractions. We conclude that in accordance with the respective prevailing turbulence structures, typical aggregation mechanisms occur that continuously generate similar distribution patterns, including particles that settle toward the river bed and thus mainly contribute to the observed decline in the total SPM concentration. In order to determine time‐scales of aggregation and related mass fluxes between the settling velocity fractions, results of model simulations were fitted to experimental data for total SPM concentration and of settling velocity frequency distributions. The comparison with simulations for the case of non‐interacting fractions clearly demonstrated the practical significance of particle interaction for a more realistic modelling of cohesive sediment and contaminant transport. Copyright © 2004 John Wiley & Sons, Ltd.     

An assessment of a new settling velocity parameterisation for cohesive sediment transport modeling

An important element within the Defra funded Estuary Process Research project “EstProc” was the implementation of the new or refined algorithms, produced under EstProc, into cohesive sediment numerical models. The implementation stage was important as any extension in the understanding of estuarine processes from EstProc was required to be suitable for dissemination into the wider research community with a level of robustness for general applications demonstrated. This report describes work undertaken to implement the new Manning Floc Settling Velocity Model, developed during EstProc. All Manning component algorithms could be combined to provide estimates of mass settling flux. The algorithms are initially assessed in a number of 1-D scenarios, where the Manning model output is compared against both real observations and the output from alternative settling parameterisations. The Manning model is then implemented into a fully 3-D computational model (TELEMAC3D) of estuarine hydraulics and sediment transport of the Lower Thames estuary. The 3-D model results with the Manning algorithm included were compared to runs with a constant settling velocity of 0.5 mm s⁻¹ and settling velocity based on a simple linear multiplier of concentration and with the above mentioned observations of suspended concentration. The findings of the 1-D case studies found the Manning empirical settling model could reproduce 93% of the total mass settling flux observed over a spring tidal cycle. The floc model fit was even better within the turbidity maximum (TM) zone. A constant 0.5 mm s⁻¹ only estimated 15% of the TM mass flux, whereas the fixed 5 mm s⁻¹ settling rate over-predicted the TM mass flux by 47%. Both settling velocity as a simple linear function of concentration, and van Leussen's method, did not fare much better estimating less than half the observed flux during the various tidal and sub-tidal cycle periods. When the Manning-settling model was applied to a layer with suspended concentrations approaching 6 g l⁻¹, it calculated 96% of the observed mass flux. The main conclusions of the implementation exercise were that it was feasible to implement a complex relationship between settling velocity and concentration in a 3-D computational model of estuarine hydraulics, without producing any significant increase in model run times or reducing model stability. The use of the Manning algorithm greatly improved the reproduction of the observed distribution of suspended concentration both in the vertical and horizontal directions compared to the other simulations. During the 1-D assessments, the Manning-settling model demonstrated flexibility in adapting to a wide range of estuarine environmental conditions (i.e. shear stress and concentration), specifically for applied modelling purposes.     

Terminal shape and velocity of a rising bubble by phase-field-based incompressible Lattice Boltzmann model

This article describes the simulation of three-dimensional buoyancy-driven bubble rise using a phase-field-based incompressible Lattice Boltzmann model. The effect of the Cahn–Hilliard mobility parameter, which is the rate of diffusion relaxation from non-equilibrium toward equilibrium state of chemical potential, is evaluated in detail. In contrast with previous work that pursues a high density ratio of binary fluids in the hydrodynamic equation, we apply a large dynamic viscosity ratio, together with a matched density pair and a separate compensating gas phase buoyant force, and the numerical results fit previous experimental results well. Through analysis, it is noted that for cases with moderate Reynolds number, a large value of mobility keeps a relatively sharp interface, while smaller values of mobility would result in diffusive interfacial regions. Moreover, for cases with large Reynolds number, small bubbles at the tail tend to separate more easily when the value of mobility is larger. This article offers some potentially useful details for performing phase-field-based simulations.     

Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system,Taiwan

A laterally averaged two-dimensional numerical model is used to simulate hydrodynamics and cohesive sediment transport in the Tanshui River estuarine system. The model handles tributaries as well as the main stem of the estuarine system. Observed time series of salinity data and tidally averaged salinity distributions have been compared with model results to calibrate the turbulent diffusion coefficients. The overall model verification is achieved with comparisons of residual currents and salinity distribution. The model reproduces the prototype water surface elevation, currents and salinity distributions. Comparisons of the suspended cohesive sediment concentrations calculated by the numerical model and the field data at various stations show good agreement. The validated model is applied to investigate the tidally averaged salinity distributions, residual circulation and suspended sediment concentration under low flow conditions in the Tanshui River estuarine system. The model results show that the limit of salt intrusion in the mainstem estuary is located at Hsin-Hai bridge in Tahan Stream, 26 km from the River mouth under Q75 flow. The null point is located at the head of salt intrusion, using 1 ppt isohaline as an indicator. The tidally averaged sediment concentration distribution exhibits a local maximum around the null point.     

An idealized model study of flocculation on sediment trapping in an estuarine turbidity maximum

Flocculation has an important impact on particle trapping in estuarine turbidity maximum (ETM) through associated increases in particle settling velocity. To quantify the importance of the flocculation processes, a size-resolved flocculation model is implemented into an ocean circulation model to simulate fine-grained particle trapping in an ETM. The model resolves the particle size from robust small flocs, about 30 μm, to very large flocs, over 1000 μm. An idealized two-dimensional model study is performed to simulate along-channel variations of suspended sediment concentrations driven by gravitational circulation and tidal currents. The results indicate that the flocculation processes play a key role in generating strong tidal asymmetrical variations in suspended sediment concentration and particle trapping. Comparison with observations suggests that the flocculation model produces realistic characteristics of an ETM.     

Calibrating Lattice Boltzmann flow simulations and estimating uncertainty in the permeability of complex porous media

A common way to simulate fluid flow in porous media is to use Lattice Boltzmann (LB) methods. Permeability predictions from such flow simulations are controlled by parameters whose settings must be calibrated in order to produce realistic modelling results. Herein we focus on the simplest and most commonly used implementation of the LB method: the single-relaxation-time BGK model. A key parameter in the BGK model is the relaxation time τ which controls flow velocity and has a substantial influence on the permeability calculation. Currently there is no rigorous scheme to calibrate its value for models of real media. We show that the standard method of calibration, by matching the flow profile of the analytic Hagen-Poiseuille pipe-flow model, results in a BGK-LB model that is unable to accurately predict permeability even in simple realistic porous media (herein, Fontainebleau sandstone). In order to reconcile the differences between predicted permeability and experimental data, we propose a method to calibrate τ using an enhanced Transitional Markov Chain Monte Carlo method, which is suitable for parallel computer architectures. We also propose a porosity-dependent τ calibration that provides an excellent fit to experimental data and which creates an empirical model that can be used to choose τ for new samples of known porosity. Our Bayesian framework thus provides robust predictions of permeability of realistic porous media, herein demonstrated on the BGK-LB model, and should therefore replace the standard pipe-flow based methods of calibration for more complex media. The calibration methodology can also be extended to more advanced LB methods.     

格子玻耳兹曼方法计算混合物整体电导率

本文采用格子玻耳兹曼方法计算混合物的整体电导率．整体电导率与各组分的电导率、体积分数和混合物结构有关．数值计算方法的有效性通过两相混合物并联或串联模型的解析解加以验证．对随机分布模型的计算发现格子玻耳兹曼方法得到的电导率落在H-S理论边界内．对Al-Bi合金整体电导率的数值模拟结果和实验结果非常相近．对饱和水岩石的介电常数的计算与实验结果相比误差较小．格子玻耳兹曼方法为混合物的整体电导率的计算提供了一个有效的途径.     

Three-dimensional numerical modeling of cohesive sediment transport and wind wave impact in a shallow oxbow lake

It was observed that in some closed inland lakes sediment transport was dominated by wind-induced currents, and the sediment resuspension was primarily driven by wind-induced waves. This paper presents the development and application of a three-dimensional numerical model for simulating cohesive sediment transport in water bodies where wind-induced currents and waves are important. In the model, the bottom shear stresses induced by currents and waves were calculated, and the processes of resuspension (erosion), deposition, settling, etc. were considered. This model was first verified by a simple test case consisting of the movement of a non-conservative tracer in a prismatic channel with uniform flow, and the model output agreed well with the analytical solution. Then it was applied to Deep Hollow Lake, a small oxbow lake in Mississippi. Simulated sediment concentrations were compared with available field observations, with generally good agreement. The transport and resuspension processes of cohesive sediment due to wind-induced current and wave in Deep Hollow Lake were also discussed.     

致密砂岩水压瞬态致裂液化格子波尔兹曼及有限元数值研究

水压瞬态致裂液化在地应力测量、 地震破坏评估和机理研究、 油气(天然气、 页岩气)及地热资源开发等地学领域都具有重要的理论和应用价值。 自20世纪60年代至今在理论和实际应用方面取得了一些奠基性和开创性成果, 但鉴于问题复杂性, 瞬态致裂液化机理至今尚不清楚。 本研究应用格子波尔兹曼及有限元多孔介质流固耦合物理模型, 对地震波载荷作用下致密砂岩水压瞬态致裂液化过程进行数值模拟研究。 首先, 以鄂尔多斯盆地某油田延长组致密砂岩为例, 利用X射线CT断层成像技术, 应用基于量子力学第一性原理格子波尔兹曼方法, 建立致密砂岩数字岩芯模型。 进而, 推导格子波尔兹曼及有限元多孔介质流固耦合数值模型公式, 建立致密砂岩水压致裂液化物理模型。 最后, 应用地震波载荷边界和初始条件, 模拟瞬态致裂液化流固耦合过程, 讨论了地震波载荷幅值、 频率及作用时间对致密砂岩孔隙结构(孔隙度大小及连通性)、 致密砂岩破裂最大主应力之间关系, 得到了地震波作用下致密砂岩致裂液化准则。     

Composite suspended sediment particles and flocculation in glacial meltwaters: preliminary evidence from Alpine and Himalayan basins

Research over the last decade has shown that the suspended sediment loads of many rivers are dominated by composite particles. These particles are also known as aggregates or flocs, and are commonly made up of constituent mineral particles, which evidence a wide range of grain sizes, and organic matter. The resulting in situ or effective particle size characteristics of fluvial suspended sediment exert a major control on all processes of entrainment, transport and deposition. The significance of composite suspended sediment particles in glacial meltwater streams has, however, not been established. Existing data on the particle size characteristics of suspended sediment in glacial meltwaters relate to the dispersed mineral fraction (absolute particle size), which, for certain size fractions, may bear little relationship to the effective or in situ distribution. Existing understanding of composite particle formation within freshwater environments would suggest that in‐stream flocculation processes do not take place in glacial meltwater systems because of the absence of organic binding agents. However, we report preliminary scanning electron microscopy data for one Alpine and two Himalayan glaciers that show composite particles are present in the suspended sediment load of the meltwater system. The genesis and structure of these composite particles and their constituent grain size characteristics are discussed. We present evidence for the existence of both aggregates, or composite particles whose features are largely inherited from source materials, and flocs, which represent composite particles produced by in‐stream flocculation processes. In the absence of organic materials, the latter may result solely from electrochemical flocculation in the meltwater sediment system. This type of floc formation has not been reported previously in the freshwater fluvial environment. Further work is needed to test the wider significance of these data and to investigate the effective particle size characteristics of suspended sediment associated with high concentration outburst events. Such events make a major contribution to suspended sediment fluxes in meltwater streams and may provide conditions that are conducive to composite particle formation by flocculation. Copyright © 2002 John Wiley & Sons, Ltd.     

A device to examine the in situ response of intertidal cohesive sediment deposits to fluid shear

A self-contained, field-portable recirculating seawater flume was designed and constructed in order to measure in situ the erodibility of cohesive estuarine muds. The flume can be easily deployed by one person and is designed for subaerial use only. Bed shear stresses up to 0.6 Nm⁻² can be generated by the flume. Rates of sediment erosion are assessed in terms of the mass of particulate material eroded with time. Flow calibrations yield a drag coefficient (C_D) of 4.5 × 10⁻³ which enables single mid-depth velocity measurements to be related directly to the imposed bed stress. Water motions within the flume are complex, but secondary (radial) circulations are not considered sufficiently competent to dominate sediment erosion. Stratification effects due to high levels of suspended sediment ( 1.5 gl⁻¹) are negligible. However, the drag-reducing properties of the sediment suspension are uncertain.