Bulletin of the American Physical Society
76th Annual Meeting of the Division of Fluid Dynamics
Sunday–Tuesday, November 19–21, 2023; Washington, DC
Session A41: Particle-Laden Flows: Particle-Resolved Simulations |
Hide Abstracts |
Chair: M. Houssem Kasbaoui, Arizona State University Room: 206 |
Sunday, November 19, 2023 8:00AM - 8:13AM |
A41.00001: Particle-Resolved DNS of heavy and neutrally-buoyant particles with the volume-filtering method M. Houssem H Kasbaoui, Marcus Herrmann We present a novel computational method for direct numerical simulations of particle-laden with fully-resolved particles (PR-DNS). The method is based on the recently developed Volume-Filtering Immersed Boundary (VFIB) method derived by volume-filtering the transport equations. This approach is mathematically and physically rigorous, in contrast to other PR-DNS methods which rely on ad-hoc numerical schemes to impose no-slip boundary conditions on the surface of particles. With the present PR-DNS based on VFIB, we show that the ratio of filter size to particle diameter acts as a slider that controls the level of fidelity. In the limit where this ratio is very small, a well-resolved PR-DNS is obtained. Conversely, when the ratio of filter size to particle diameter is large, a classic point-particle method is obtained. The discretization of the filtered equations is discussed and compared to other PR-DNS strategies based on direct-forcing immersed boundary methods. Numerical examples with sedimenting resolved spheres show that the method is very robust, stable, and accurate, even with neutrally buoyant particles. The volume-filtering formulations enables multi-resolution simulations where, depending on the local resolution, the particles are represented as resolved spheres or Lagrangian particles. |
Sunday, November 19, 2023 8:13AM - 8:26AM |
A41.00002: Particle-resolved simulations of the wall effect on the sedimentation of a single sphere in yield-stress fluids Mohammad Sarabian, Marco E. Rosti, Luca Brandt We perform high-fidelity numerical simulations to investigate the confinement effect on the sedimentation of a single spherical particle in an otherwise quiescent yield stress fluid. The simulations are performed in the presence of finite elasticity and weak inertia. The carrier fluid is modeled utilizing the elastoviscoplastic constitutive laws proposed by Saramito (2009). The additional elastic stress tensor is fully coupled with the flow equation, while the rigid particle is represented by an immersed boundary method. The particle-resolved simulations demonstrate the faster relaxation of the fluid velocity and the progressive translation of the location of the negative wake downstream of the sphere as the bounding walls are brought closer to the particle. Furthermore, the drag force exerted from the surrounding fluid on the particle decreases by escalating the particle-wall distance. We show that the confinement ratio (ratio of the gap between rigid confining walls and the sphere radius) reaches a critical value beyond which the wall-effect on the particle and flow dynamics becomes negligible. The key finding here is that the critical confinement ratio and the maximum variation of the Stokes drag with confinement ratio are weakly dependent on the level of material elasticity and plasticity for a certain range of material parameters. Finally, we propose an expression for the Stokes drag coefficient, as a function of material plasticity and confinement ratio. |
Sunday, November 19, 2023 8:26AM - 8:39AM |
A41.00003: Discrete Transport from a Tumbling Prolate Spheroid in a Simple Shear Flow Yanxing Wang, Hui Wan, Tie Wei, Fangjun Shu Through high-fidelity numerical simulation based on a multi-grid strategy with the lattice Boltzmann framework, an in-depth study on the heat and mass transport from a prolate spheroid suspended in a simple shear flow has been carried out. In the simulation, the temperature and mass concentration are modeled as a passive scalar released at the spheroid's surface. The interaction between the carrier fluid and the suspended spheroid, as well as the resultant scalar transport process, have been investigated. The interaction creates several particle-scale flow modes in the fluid layer surrounding the spheroid. When the advective scalar transport is weak, the passive scalar released at the spheroid surface is transported through the fluid layer to the outer passing fluid via diffusion, and then is transported downstream by the passing fluid via flow advection. When the advective scalar transport is strong, most of the passive scalar is gained by the fluid layer surrounding the spheroid. The multiple flow modes of the fluid layer and the outer passing fluid lead to multiple advective scalar transport modes, constituting a spatiotemporally discrete set of scalar transport mechanisms for the prolate spheroid in a simple shear flow. Factors such as spheroid-to-fluid density ratio, spheroid aspect ratio, and Reynolds number, affect the behaviors of the fluid and spheroid, thus have influences on the scalar transport. |
Sunday, November 19, 2023 8:39AM - 8:52AM |
A41.00004: Interface-resolved simulation of particles with different shapes in sediment transport Yuanqing Liu, Lian Shen Sediment transport is a dynamic process with significant implications for landscapes and human activities. The movement and collision of sediment particles have a profound impact on the surrounding flow field. Among various flow configurations and particle properties, the shape of particles plays a crucial role in the fluid-particle interaction process. In this work, we employ numerical simulations to investigate sedimentation transport for different particle shapes, including spherical particles, oblate particles, and their mixtures. |
Sunday, November 19, 2023 8:52AM - 9:05AM |
A41.00005: Resolved CFD-DEM simulation of irregularly shaped cuttings removal in deep hard-rock drilling Pasha Piroozmand, Dany Suter, Davood Farshi Deep geothermal energy is a promising renewable energy source that harnesses the Earth's natural heat from deep beneath the surface. However, drilling in deep geothermal reservoirs, especially those exceeding 5 km in depth, poses significant challenges to conventional rotary mechanical drilling due to the presence of hard crystalline rocks. To address this, Plasma-Pulse Geo-Drilling (PPGD) has been proposed, which uses high-voltage pulses to remove rocks. Nevertheless, the significant size of the cuttings generated during the drilling, reaching up to 2 cm in diameter, poses a difficult hurdle for effective removal through injection of pure water. In this study, we examine the cutting removal efficiency of PPGD by employing a Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) approach. A highly turbulent flow regime necessitates highly refined grid cells, which are significantly smaller than particle sizes. Therefore, a resolved CFD-DEM approach based on a fictitious domain method is used in the open-source C++ package CFDEM, which integrates the URANS solver in OpenFOAM with the DEM solver in LIGGGHTS. The arbitrary shape of the cuttings is represented using a multi-sphere method. The simulation enables the exploration of critical design considerations, including feasibility assessments, determination of the optimal power required for efficient particle removal under diverse operating conditions, and geometrical modification for improving the removal process. |
Sunday, November 19, 2023 9:05AM - 9:18AM |
A41.00006: Towards a more realistic configuration to simulate the settling of finite-size particles in the dilute regime Manuel Moriche Guerrero, Manuel García-Villalba, Markus Uhlmann Simulating the settling of particles in the dilute regime is a very complex problem which, thanks to the relatively recent maturity reached by immersed boundary methods, we can analyze by means of particle-resolved direct numerical simulations. The cost of these simulations is still very high and one of the assumptions typically made by the authors to reduce the complexity of the problem is the use of a triply periodic configuration. This configuration has shown to be very successful to understand the onset of clustering, but the analysis of cluster dynamics is not possible because of the strong correlation of the solution in the vertical direction. In this talk we present preliminary results of the settling of particles in a configuration where the periodicity in the vertical direction is removed. The advantages and challenges of using this configuration will be also briefly discussed. In order to exploit our previous knowledge on the settling of particles we set the Galileo number to 121 and 178 and the particle-to-fluid density ratio to 1.5 based on the work of Uhlmann and Doychev (2014). We observe a similar clustering behavior as in the triply periodic counterpart, but at a higher particle concentration. |
Sunday, November 19, 2023 9:18AM - 9:31AM |
A41.00007: High Resolution Simulations of Porous Fractal Aggregates Settling through Density Gradients Zachary Maches, Eckart Meiburg We perform Particle-Resolved Direct Numerical Simulations (PR-DNS) on rigidly-bonded fractal aggregates of spherical particles settling across varying density gradients. These simulations are performed to study the impact of the complex, fractal geometry on the settling of the aggregate, particularly as it passes into a region of high density fluid while retaining less dense fluid in the pore space between individual particles. The study of these aggregates and their behaviors has applications in a broad range of problems in environmental fluid dynamics. We investigate over many simulations the velocity and rotation of the aggregates, varying the number of particles N in an aggregate, the aggregate’s fractal dimension D_{f}, and the associated fractal prefactor k_{f}, which together all characterize the overall structure of the cluster. We compare these results with the settling of spheres of comparable volume to the aggregate, to discuss the effect that the shape of the aggregate has on its motion. From these results, we intend to develop general relationships between the size of the aggregate and its fractal parameters to its overall settling velocity during sedimentation. |
Sunday, November 19, 2023 9:31AM - 9:44AM |
A41.00008: Flocculation of Finite Sized Particles in Turbulent Channel Flow Alexandre D Leonelli, Eckart Meiburg We present a computational study of the flocculation dynamics and turbulence modulation in a channel flow laden with small, finite size, cohesive particles. We conduct a series of four way coupled, grain-resolved direct numerical simulations, across which the strength of the short-range inter-particle cohesive force is varied, in order to determine its effect on the flocculation process. The strength of cohesion is shown to govern the size (O(1-100) particle diameters) and number of aggregates. Additionally, we analyze specific floc properties such as size, shape, and lifespan as a function of the wall normal distance. The higher mean velocity gradient and turbulent activity near the wall is found to both facilitate aggregation by bringing particles together, while also prompting breakup. Finally, enabled by the fully coupled nature of the simulations, we quantify the effect of aggregation on the turbulence by reporting classical first and second order turbulence statistics for both the fluid and particles. |
Sunday, November 19, 2023 9:44AM - 9:57AM |
A41.00009: Particle-Resolved DNS of flow modulation induced by finite-size particles dispersed in a turbulent channel Jonathan Van Doren, Houssem Kasbaoui The modification of turbulent structures by finite size particles is of interest to pipeline design for oil and slurry transport. The interaction of finite size particles with turbulent flow structures is investigated in Particle-Resolved Direct Numerical Simulations of turbulent channel flow at friction Reynolds number 180 laden with spherical particles with diameter equal to 12 wall units and dispersed at the volume fraction 1%. To understand the role of particle inertia in modulating turbulence, particle-to-fluid density ratio is varied from 1, i.e., neutrally buoyant particles, up to 4. About 2320 resolved particles are tracked throughout the computational domain and represented using the recently developed volume-filtering immersed boundary method to account for stresses at the fluid solid interface. Additional turbulence statistics are collected for a particle free channel for comparison with the particle laden channels. |
Follow Us |
Engage
Become an APS Member |
My APS
Renew Membership |
Information for |
About APSThe American Physical Society (APS) is a non-profit membership organization working to advance the knowledge of physics. |
© 2024 American Physical Society
| All rights reserved | Terms of Use
| Contact Us
Headquarters
1 Physics Ellipse, College Park, MD 20740-3844
(301) 209-3200
Editorial Office
100 Motor Pkwy, Suite 110, Hauppauge, NY 11788
(631) 591-4000
Office of Public Affairs
529 14th St NW, Suite 1050, Washington, D.C. 20045-2001
(202) 662-8700