Settling of cohesive sediment: particle-resolved simulations

We develop a physical and computational model for performing fully coupled, particle-resolved Direct Numerical Simulations of cohesive sediment, based on the Immersed Boundary Method. The model distributes the cohesive forces over a thin shell surrounding each particle, thereby allowing for the spatial and temporal resolution of the cohesive forces during particle-particle interactions. The influence of the cohesive forces is captured by a single dimensionless parameter in the form of a cohesion number, which represents the ratio of cohesive and gravitational forces acting on a particle. We test and validate the cohesive force model for binary particle interactions in the Drafting-Kissing-Tumbling (DKT) configuration. The DKT simulations demonstrate that cohesive particle pairs settle in a preferred orientation, with particles of very different sizes preferentially aligning themselves in the vertical direction, so that the smaller particle is drafted in the wake of the larger one. To test this mechanism in a system of higher complexity, we perform large simulations of 1,261 polydisperse settling particles starting from rest. These simulations reproduce several earlier experimental observations by other authors, such as the accelerated settling of sand and silt particles due to particle bonding. The simulations demonstrate that cohesive forces accelerate the overall settling process primarily because smaller grains attach to larger ones and settle in their wakes. For the present cohesion number values, we observe that settling can be accelerated by up to 29%. We propose physically based parametrization of classical hindered settling functions proposed by earlier authors, in order to account for cohesive forces. An investigation of the energy budget shows that the work of the collision forces can substantially modify the relevant energy conversion processes.Comment: 39 page

Rheology of mobile sediment beds sheared by viscous, pressure-driven flows

We present a detailed comparison of the rheological behaviour of sheared sediment beds in a pressure-driven, straight channel configuration based on data that was generated by means of fully coupled, grain-resolved direct numerical simulations and experimental measurements reviously published by Aussillous {\it et al.} (J. Fluid Mech., vol. 736, 2013, pp. 594-615). The highly-resolved simulation data allows to compute the stress balance of the suspension in the streamwise and vertical directions and the stress exchange between the fluid and particle phase, which is information needed to infer the rheology, but has so far been unreachable in experiments. Applying this knowledge to the experimental and numerical data, we obtain the statistically-stationary, depth-resolved profiles of the relevant rheological quantities. The scaling behavior of rheological quantities such as the shear and normal viscosities and the effective friction coefficient are examined and compared to data coming from rheometry experiments and from widely-used rheological correlations. We show that rheological properties that have previously been inferred for annular Couette-type shear flows with neutrally buoyant particles still hold for our setup of sediment transport in a Poiseuille flow and in the dense regime we found good agreement with empirical relationships derived therefrom. Subdividing the total stress into parts from particle contact and hydrodynamics suggests a critical particle volume fraction of 0.3 to separate the dense from the dilute regime. In the dilute regime, i.e., the sediment transport layer, long-range hydrodynamic interactions are screened by the porous media and the effective viscosity obeys the Einstein relation

Confronting Grand Challenges in Environmental Fluid Dynamics

Environmental fluid dynamics underlies a wealth of natural, industrial and, by extension, societal challenges. In the coming decades, as we strive towards a more sustainable planet, there are a wide range of grand challenge problems that need to be tackled, ranging from fundamental advances in understanding and modeling of stratified turbulence and consequent mixing, to applied studies of pollution transport in the ocean, atmosphere and urban environments. A workshop was organized in the Les Houches School of Physics in France in January 2019 with the objective of gathering leading figures in the field to produce a road map for the scientific community. Five subject areas were addressed: multiphase flow, stratified flow, ocean transport, atmospheric and urban transport, and weather and climate prediction. This article summarizes the discussions and outcomes of the meeting, with the intent of providing a resource for the community going forward

Confronting Grand Challenges in environmental fluid mechanics

Environmental fluid mechanics underlies a wealth of natural, industrial and, by extension, societal challenges. In the coming decades, as we strive towards a more sustainable planet, there are a wide range of grand challenge problems that need to be tackled, ranging from fundamental advances in understanding and modeling of stratified turbulence and consequent mixing, to applied studies of pollution transport in the ocean, atmosphere and urban environments. A workshop was organized in the Les Houches School of Physics in France in January 2019 with the objective of gathering leading figures in the field to produce a road map for the scientific community. Five subject areas were addressed: multiphase flow, stratified flow, ocean transport, atmospheric and urban transport, and weather and climate prediction. This article summarizes the discussions and outcomes of the meeting, with the intent of providing a resource for the community going forward

Parameter subset reduction for imaging-based digital twin generation of patients with left ventricular mechanical discoordination

Background: Integration of a patient’s non-invasive imaging data in a digital twin (DT) of the heart can provide valuable insight into the myocardial disease substrates underlying left ventricular (LV) mechanical discoordination. However, when generating a DT, model parameters should be identifiable to obtain robust parameter estimations. In this study, we used the CircAdapt model of the human heart and circulation to find a subset of parameters which were identifiable from LV cavity volume and regional strain measurements of patients with different substrates of left bundle branch block (LBBB) and myocardial infarction (MI). To this end, we included seven patients with heart failure with reduced ejection fraction (HFrEF) and LBBB (study ID: 2018-0863, registration date: 2019–10–07), of which four were non-ischemic (LBBB-only) and three had previous MI (LBBB-MI), and six narrow QRS patients with MI (MI-only) (study ID: NL45241.041.13, registration date: 2013–11–12). Morris screening method (MSM) was applied first to find parameters which were important for LV volume, regional strain, and strain rate indices. Second, this parameter subset was iteratively reduced based on parameter identifiability and reproducibility. Parameter identifiability was based on the diaphony calculated from quasi-Monte Carlo simulations and reproducibility was based on the intraclass correlation coefficient (ICC) obtained from repeated parameter estimation using dynamic multi-swarm particle swarm optimization. Goodness-of-fit was defined as the mean squared error (χ2) of LV myocardial strain, strain rate, and cavity volume. Results: A subset of 270 parameters remained after MSM which produced high-quality DTs of all patients (χ2 < 1.6), but minimum parameter reproducibility was poor (ICCmin = 0.01). Iterative reduction yielded a reproducible (ICCmin = 0.83) subset of 75 parameters, including cardiac output, global LV activation duration, regional mechanical activation delay, and regional LV myocardial constitutive properties. This reduced subset produced patient-resembling DTs (χ2 < 2.2), while septal-to-lateral wall workload imbalance was higher for the LBBB-only DTs than for the MI-only DTs (p < 0.05). Conclusions: By applying sensitivity and identifiability analysis, we successfully determined a parameter subset of the CircAdapt model which can be used to generate imaging-based DTs of patients with LV mechanical discoordination. Parameters were reproducibly estimated using particle swarm optimization, and derived LV myocardial work distribution was representative for the patient’s underlying disease substrate. This DT technology enables patient-specific substrate characterization and can potentially be used to support clinical decision making

A settling-driven instability in two-component, stably stratified fluids

We analyse the linear stability of stably stratified fluids whose density depends on two scalar fields where one of the scalar fields is unstably stratified and involves a settling velocity. Such conditions may be found, for example, in flows involving the transport of sediment in addition to heat or salt. A linear stability analysis for constant-gradient base states demonstrates that the settling velocity generates a phase shift between the perturbation fields of the two scalars, which gives rise to a novel, settling-driven instability mode. This instability mechanism favours the growth of waves that are inclined with respect to the horizontal. It is active for all density and diffusivity ratios, including for cases in which the two scalars diffuse at identical rates. If the scalars have unequal diffusivities, it competes with the elevator mode waves of the classical double-diffusive instability. We present detailed linear stability results as a function of the governing dimensionless parameters, including for lateral gradients of the base state density fields that result in predominantly horizontal intrusion instabilities. Highly resolved direct numerical simulation results serve to illustrate the nonlinear competition of the various instabilities for such flows in different parameter regimes