Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA coupled finite volume and material point method for two-phase simulation of liquid-sediment and gas-sediment flows
84
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Mixtures of fluids and granular sediments play an important role in many industrial, geotechnical, and aerospace engineering problems, from waste management and transportation (liquid--sediment mixtures) to dust kick-up below helicopter rotors (gas--sediment mixtures). These mixed flows often involve bulk motion of hundreds of billions of individual sediment particles and can contain both highly turbulent regions and static, non-flowing regions. This breadth of phenomena necessitates the use of continuum simulation methods, such as the material point method (MPM), which can accurately capture these large deformations while also tracking the Lagrangian features of the flow (e.g. the granular surface, elastic stress, etc.). Recent works using two-phase MPM frameworks to simulate these mixtures have shown substantial promise; however, these approaches are hindered by the numerical limitations of MPM when simulating pure fluids. In addition to the well-known particle ringing instability and difficulty defining inflow/outflow boundary conditions, MPM has a tendency to accumulate quadrature errors as materials deform, increasing the rate of overall error growth as simulations progress. In this work, we present an improved, two-phase continuum simulation framework that uses the finite volume method (FVM) to solve the fluid phase equations of motion and MPM to solve the solid phase equations of motion, substantially reducing the effect of these errors and providing better accuracy and stability for long-duration simulations of these mixtures.
rate research
Read More
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.
Capturing the interaction between objects that have an extreme difference in Young s modulus or geometrical scale is a highly challenging topic for numerical simulation. One of the fundamental questions is how to build an accurate multi-scale method with optimal computational efficiency. In this work, we develop a material-point-spheropolygon discrete element method (MPM-SDEM). Our approach fully couples the material point method (MPM) and the spheropolygon discrete element method (SDEM) through the exchange of contact force information. It combines the advantage of MPM for accurately simulating elastoplastic continuum materials and the high efficiency of DEM for calculating the Newtonian dynamics of discrete near-rigid objects. The MPM-SDEM framework is demonstrated with an explicit time integration scheme. Its accuracy and efficiency are further analysed against the analytical and experimental data. Results demonstrate this method could accurately capture the contact force and momentum exchange between materials while maintaining favourable computational stability and efficiency. Our framework exhibits great potential in the analysis of multi-scale, multi-physics phenomena
We present a thermodynamically consistent constitutive model for fluid-saturated sediments, spanning dense to dilute regimes, developed from the basic balance laws for two phase-mixtures. The model can represent various limiting cases, such as pure fluid and dry grains. It is formulated to capture a number of key behaviors such as: (i) viscous inertial rheology of submerged wet grains under steady shearing flows, (ii) the critical state behavior of grains, which causes granular Reynolds dilation/contraction due to shear, (iii) the viscous thickening of the fluid response due to the presence of suspended grains, and (iv) the Darcy-like drag interaction observed in both dense and dilute mixtures, which gives rise to complex fluid-grain interactions under dilation and flow. The full constitutive model is combined with the basic equations of motion for each mixture phase and implemented in the material point method (MPM) to accurately model the coupled dynamics of the mixed system. Qualitative results show the breadth of problems which this model can address. Quantitative results demonstrate the accuracy of this model as compared with analytical limits and experimental observations of fluid and grain behaviors in inhomogeneous geometries.
A hybrid parallel approach for fully resolved simulations of particle-laden flows in sediment transport is proposed. To overcome the challenges of load imbalance in the traditional domain decomposition method when encountering highly uneven distributions of particles in space, we develop a hybrid parallel approach adopting the domain decomposition method for the carrier phase and a mirror domain technique for the disperse phase. We modify the mirror domain technique originally developed for point particles to fully resolved particle simulations, which are more challenging since a finite-sized particle may be split into different subdomains; thus, more complex treatments of particle-fluid interactions are needed. By utilizing the mirror domain technique, in which each processor handles nearly the same number of particles regardless of the particle spatial distribution, excellent load balance is achieved. The present hybrid parallel approach also shows strong scalability and high parallel efficiency in a test of a fully resolved simulation case of sediment transport. Furthermore, a novel memory optimization method is proposed for spherical particles of equal size, which can substantially reduce the memory cost and enable the simulation of millions of fully resolved particles on a common highly parallel computing platform. Our code is validated by several benchmark cases, and the results show good agreement with experimental and computational data in the literature.
We investigate a mathematical theory for the erosion of sediment which begins with the study of a non-linear, parabolic, weighted 4-Laplace equation on a rectangular domain corresponding to a base segment of an extended landscape. Imposing natural boundary conditions, we show that the equation admits entropy solutions and prove regularity and uniqueness of weak solutions when they exist. We then investigate a particular class of weak solutions studied in previous work of the first author and produce numerical simulations of these solutions. After introducing an optimal transportation problem for the sediment flow, we show that this class of weak solutions implements the optimal transportation of the sediment.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
