Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hybrid parallel approach for fully resolved simulations of particle-laden flows in sediment transport
73
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
A general, two-way coupled, point-particle formulation that accounts for the disturbance created by the dispersed particles in obtaining the undisturbed fluid flow field needed for accurate computation of the force closure models is presented. Specifically, equations for the disturbance field created by the presence of particles are first derived based on the inter-phase momentum coupling force in a finite-volume formulation. Solution to the disturbance field is obtained using two approaches: (i) direct computation of the disturbance velocity and pressure using the reaction force due to particles at computational control volumes, and (ii) a linearized, approximate computation of the disturbance velocity field, specifically applicable for low Reynolds number flows. In both approaches, the computed disturbance field is used to obtain the undisturbed fluid velocity necessary to model the aerodynamic forces on the particle. The two approaches are thoroughly evaluated for a single particle in an unbounded and wall-bounded flow on uniform, anisotropic, as well as unstructured grids to show accurate computation of the particle motion and inter-phase coupling. The approach is straightforward and can be applied to any numerical formulation for particle-laden flows including Euler-Lagrange as well as Euler-Euler formulations.
We present a comparison of different particles velocity and acceleration statistics in two paradigmatic turbulent swirling flows: the von Karman flow in a laboratory experiment, and the Taylor-Green flow in direct numerical simulations. Tracers, as well as inertial particles, are considered. Results indicate that, in spite of the differences in boundary conditions and forcing mechanisms, scaling properties and statistical quantities reveal similarities between both flows, pointing to new methods to calibrate and compare models for particles dynamics in numerical simulations, as well as to characterize the dynamics of particles in simulations and experiments.
The hydrodynamics of a liquid-vapour interface in contact with an heterogeneous surface is largely impacted by the presence of defects at the smaller scales. Such defects introduce morphological disturbances on the contact line and ultimately determine the force exerted on the wedge of liquid in contact with the surface. From the mathematical point of view, defects introduce perturbation modes, whose space-time evolution is governed by the interfacial hydrodynamic equations of the contact line. In this paper we derive the response function of the contact line to such generic perturbations. The contact line response may be used to design simplified 1+1 dimensional models accounting for the complexity of interfacial flows coupled to nanoscale defects, yet offering a more tractable mathematical framework to include thermal fluctuations and explore thermally activated contact line motion through a disordered energy landscape.
In this paper a fully Eulerian solver for the study of multiphase flows for simulating the propagation of surface gravity waves over submerged bodies is presented. We solve the incompressible Navier-Stokes equations coupled with the volume of fluid technique for the modeling of the liquid phases with the interface, an immersed body method for the solid bodies and an iterative strong-coupling procedure for the fluid-structure interaction. The flow incompressibility is enforced via the solution of a Poisson equation which, owing to the density jump across the interfaces of the liquid phases, has to resort to the splitting procedure of Dodd & Ferrante [12]. The solver is validated through comparisons against classical test cases for fluid-structure interaction like migration of particles in pressure-driven channel, multiphase flows, water exit of a cylinder and a good agreement is found for all tests. Furthermore, we show the application of the solver to the case of a surface gravity wave propagating over a submerged reversed pendulum and verify that the solver can reproduce the energy exchange between the wave and the pendulum. Finally the three-dimensional spilling breaking of a wave induced by a submerged sphere is considered.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
