In this paper, a three‐dimensional numerical solver is developed for suspensions of rigid and soft particles and droplets in viscoelastic and elastoviscoplastic (EVP) fluids. The presented algorithm is designed to allow for the first time three‐dimensional simulations of inertial and turbulent EVP fluids with a large number particles and droplets. This is achieved by combining fast and highly scalable methods such as an FFT‐based pressure solver, with the evolution equation for non‐Newtonian (including EVP) stresses. In this flexible computational framework, the fluid can be modeled by either Oldroyd‐B, neo‐Hookean, FENE‐P, or Saramito EVP models, and the additional equations for the non‐Newtonian stresses are fully coupled with the flow. The rigid particles are discretized on a moving Lagrangian grid, whereas the flow equations are solved on a fixed Eulerian grid. The solid particles are represented by an immersed boundary method with a computationally efficient direct forcing method, allowing simulations of a large numbers of particles. The immersed boundary force is computed at the particle surface and then included in the momentum equations as a body force. The droplets and soft particles on the other hand are simulated in a fully Eulerian framework, the former with a level‐set method to capture the moving interface and the latter with an indicator function. The solver is first validated for various benchmark single‐phase and two‐phase EVP flow problems through comparison with data from the literature. Finally, we present new results on the dynamics of a buoyancy‐driven drop in an EVP fluid.
This content will become publicly available on July 19, 2025
We propose a novel Particle Flow Map (PFM) method to enable accurate long-range advection for incompressible fluid simulation. The foundation of our method is the observation that a particle trajectory generated in a forward simulation naturally embodies a perfect flow map. Centered on this concept, we have developed an Eulerian-Lagrangian framework comprising four essential components: Lagrangian particles for a natural and precise representation of bidirectional flow maps; a dual-scale map representation to accommodate the mapping of various flow quantities; a particle-to-grid interpolation scheme for accurate quantity transfer from particles to grid nodes; and a hybrid impulse-based solver to enforce incompressibility on the grid. The efficacy of PFM has been demonstrated through various simulation scenarios, highlighting the evolution of complex vortical structures and the details of turbulent flows. Notably, compared to NFM, PFM reduces computing time by up to 49 times and memory consumption by up to 41%, while enhancing vorticity preservation as evidenced in various tests like leapfrog, vortex tube, and turbulent flow.
more » « less- Award ID(s):
- 1919647
- PAR ID:
- 10546830
- Publisher / Repository:
- ACM Transactions on Graphics (SIGGRAPH 2024)
- Date Published:
- Journal Name:
- ACM Transactions on Graphics
- Volume:
- 43
- Issue:
- 4
- ISSN:
- 0730-0301
- Page Range / eLocation ID:
- 1 to 20
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Summary -
more » « less
The correlation between helicity and turbulent transport in turbulent flows is probed with the use of direct numerical simulation and Lagrangian scalar tracking. Channel flow and plane Couette flow at friction Reynolds number 300 and Lagrangian data along the trajectories of fluid particles and passive particles with Schmidt numbers 0.7 and 6 are used. The goal is to identify characteristics of the flow that enhance turbulent transport from the wall, and how flow regions that exhibit these characteristics are related to helicity. The relationship between vorticity and relative helicity along particle trajectories is probed, and the relationship between the distribution of helicity conditioned on Reynolds stress quadrants is also evaluated. More importantly, the correlation between relative helicity density and the alignment of vorticity with velocity vectors and eigenvectors of the rate of strain tensor is presented. Separate computations for particles that disperse the farthest into the flow field and those that disperse the least are conducted to determine the flow structures that contribute to turbulent dispersion. The joint distribution of helicity and vertical velocity, and helicity and vertical vorticity depends on the location of particle release and the Schmidt number. The trajectories of particles that disperse the least are characterized by a correlation between the absolute value of the relative helicity density and the absolute value of the cosine between the vorticity vector and the eigenvectors of the rate of strain tensor, while the value of this correlation approaches zero for the particles that disperse the most.
-
more » « less
Abstract In shallow coastal oceans, turbulent flows driven by surface winds and waves and constrained by a solid bottom disperse particles. This work examines the mechanisms driving horizontal and vertical dispersion of buoyant and sinking particles for times much greater than turbulent integral time scales. Turbulent fields are modeled using a wind‐stress driven large eddy simulation (LES), incorporating wave‐driven Langmuir turbulence, surface breaking wave turbulent kinetic energy inputs, and a solid bottom boundary. A Lagrangian stochastic model is paired to the LES to incorporate Lagrangian particle tracking. Within a subset of intermediate buoyant rise velocities, particles experience synergistic vertical mixing in which breaking waves (BW) inject particles into Langmuir downwelling velocities sufficient to drive deep mixing. Along‐wind dispersion is controlled by vertical shear in mean along‐wind velocities. Wind and bottom friction‐driven vertical shear enhances dispersion of buoyant and sinking particles, while energetic turbulent mixing, such as from BW, dampens shear dispersion. Strongly rising and sinking particles trapped at the ocean surface and bottom, respectively, experience no vertical shear, resulting in low rates of along‐wind dispersion. Crosswind dispersion is shaped by particle advection in wind‐aligned fields of counter‐rotating Langmuir and Couette roll cells. Langmuir cells enhance crosswind dispersion in neutrally to intermediately buoyant particles through enhanced cell hopping. Surface trapping restricts particles to Langmuir convergence regions, strongly inhibiting crosswind dispersion. In shallow coastal systems, particle dispersion depends heavily on particle buoyancy and wave‐dependent turbulent effects.
-
We introduce Aquila-LCS, GPU and CPU optimized object-oriented, in-house codes for volumetric particle advection and 3D Finite-Time Lyapunov Exponent (FTLE) and Finite-Size Lyapunov Exponent (FSLE) computations. The purpose is to analyze 3D Lagrangian Coherent Structures (LCS) in large Direct Numerical Simulation (DNS) data. Our technique uses advanced search strategies for quick cell identification and efficient storage techniques. This solver scales effectively on both GPUs (up to 62 Nvidia V100 GPUs) and multi-core CPUs (up to 32,768 CPU cores), tracking up to 8-billion particles. We apply our approach to four turbulent boundary layers at different flow regimes and Reynolds numbers.more » « less
-
more » « less
In this work, we introduce a scalable and efficient GPU-accelerated methodology for volumetric particle advection and finite-time Lyapunov exponent (FTLE) calculation, focusing on the analysis of Lagrangian coherent structures (LCS) in large-scale direct numerical simulation (DNS) datasets across incompressible, supersonic, and hypersonic flow regimes. LCS play a significant role in turbulent boundary layer analysis, and our proposed methodology offers valuable insights into their behavior in various flow conditions. Our novel owning-cell locator method enables efficient constant-time cell search, and the algorithm draws inspiration from classical search algorithms and modern multi-level approaches in numerical linear algebra. The proposed method is implemented for both multi-core CPUs and Nvidia GPUs, demonstrating strong scaling up to 32,768 CPU cores and up to 62 Nvidia V100 GPUs. By decoupling particle advection from other problems, we achieve modularity and extensibility, resulting in consistent parallel efficiency across different architectures. Our methodology was applied to calculate and visualize the FTLE on four turbulent boundary layers at different Reynolds and Mach numbers, revealing that coherent structures grow more isotropic proportional to the Mach number, and their inclination angle varies along the streamwise direction. We also observed increased anisotropy and FTLE organization at lower Reynolds numbers, with structures retaining coherency along both spanwise and streamwise directions. Additionally, we demonstrated the impact of lower temporal frequency sampling by upscaling with an efficient linear upsampler, preserving general trends with only 10% of the required storage. In summary, we present a particle search scheme for particle advection workloads in the context of visualizing LCS via FTLE that exhibits strong scaling performance and efficiency at scale. Our proposed algorithm is applicable across various domains, requiring efficient search algorithms in large, structured domains. While this article focuses on the methodology and its application to LCS, an in-depth study of the physics and compressibility effects in LCS candidates will be explored in a future publication.
