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.
We introduce Neural Flow Maps, a novel simulation method bridging the emerging paradigm of implicit neural representations with fluid simulation based on the theory of flow maps, to achieve state-of-the-art simulation of in-viscid fluid phenomena. We devise a novel hybrid neural field representation, Spatially Sparse Neural Fields (SSNF), which fuses small neural networks with a pyramid of overlapping, multi-resolution, and spatially sparse grids, to compactly represent long-term spatiotemporal velocity fields at high accuracy. With this neural velocity buffer in hand, we compute long-term, bidirectional flow maps and their Jacobians in a mechanistically symmetric manner, to facilitate drastic accuracy improvement over existing solutions. These long-range, bidirectional flow maps enable high advection accuracy with low dissipation, which in turn facilitates high-fidelity incompressible flow simulations that manifest intricate vortical structures. We demonstrate the efficacy of our neural fluid simulation in a variety of challenging simulation scenarios, including leapfrogging vortices, colliding vortices, vortex reconnections, as well as vortex generation from moving obstacles and density differences. Our examples show increased performance over existing methods in terms of energy conservation, visual complexity, adherence to experimental observations, and preservation of detailed vortical structures.
more » « less- PAR ID:
- 10518997
- Publisher / Repository:
- ACM Transactions on Graphics (TOG)
- Date Published:
- Journal Name:
- ACM Transactions on Graphics
- Volume:
- 42
- Issue:
- 6
- ISSN:
- 0730-0301
- Page Range / eLocation ID:
- 1 to 21
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We propose a novel Clebsch method to simulate the free-surface vortical flow. At the center of our approach lies a level-set method enhanced by a wave-function correction scheme and a wave-function extrapolation algorithm to tackle the Clebsch method's numerical instabilities near a dynamic interface. By combining the Clebsch wave function's expressiveness in representing vortical structures and the level-set function's ability on tracking interfacial dynamics, we can model complex vortex-interface interaction problems that exhibit rich free-surface flow details on a Cartesian grid. We showcase the efficacy of our approach by simulating a wide range of new free-surface flow phenomena that were impractical for previous methods, including horseshoe vortex, sink vortex, bubble rings, and free-surface wake vortices.more » « less
-
Abstract The identification of vortices in a fluid flow is a dynamically interesting problem that has practical applications in oceanography due to the outsized role eddies play in water mass, heat, and tracer transport. Here a new Eulerian scheme is developed to detect both vortices and strongly strained fronts, which are both ubiquitous in the world ocean. The new scheme is conceptually linked to the well-known Okubo-Weiss parameter, but is extended to quasigeostrophic flows by recognizing the strong role played by vertical shear in ocean dynamics. Adapted from the λ 2 -criterion for vortex identification (Jeong and Hussain 1995), the scheme considers the curvature of the pressure field as the differentiator between vortical and strained flow structures, and it is shown that its underlying geometry also exhibits characteristics of quasigeostrophic flow. The uses and skill of the scheme are demonstrated using a high-resolution regional ocean simulation, and prospects for its use with observational products are discussed.more » « less
-
Flow-induced vibrations of flexible surfaces driven by coherent vortical structures are ubiquitous in engineering and biological flows; from the extraction of fluidic energy via oscillating electro-active polymers to vocal fold dynamics during voiced speech production. These scenarios may involve either discrete or periodic loading conditions due to the advection of vortices past the structure. This work considers, as a function of the vortex production frequency, the fluid-structure interaction that occurs as vortices are propagated tangentially over flexible plates with variable structural properties. Velocity fields are acquired with particle image velocimetry and used to compute the vorticity and pressure fields, while the plate energy is estimated from its kinematics. Primary and secondary peaks in plate deflection amplitude and the plate energy as a function of vortex production frequency are observed at integer fractions of the fundamental plate frequency. At resonance conditions, plate energy relative to discrete vortex loading is increased by approximately three orders of magnitude, while the sub-harmonics increase the plate energy by about two orders of magnitude. Additional physical influences on the energy exchange process, including vortex-to-plate spacing and Strouhal number, are also investigated, detailing the importance of spatial and temporal interactions. The magnitude of the initial plate deflection as the vortex ring approaches the plate, due to persistent vibrations from previous interactions, is shown to retard the time at which the maximum load is applied as the increased relative vortex-to-plate spacing weakens cross-sign vorticity interactions. Finally, plate properties are scaled to model the structural properties of the vocal folds and the effect of intra-glottal vortices on vocal fold dynamics is quantified, where a negligible influence is observed.more » « less
-
Modal analysis on micro-vortex generator (MVG)-controlled supersonic flow at different Mach numbers is performed in this paper. The purpose of this investigation is to clarify the different properties of streamwise and ring-like vortical modes, and the effects of different Mach numbers on these modes, to further understand the vortical structures as they travel from MVG down to the shock wave/boundary-layer interaction (SWBLI) region. To this end, a high order and high resolution large eddy simulation (LES) was carried out, which identified the vortical structures behind the MVG and in the shock wave/boundary-layer interaction (SWBLI) region in the supersonic ramp flow with flow speeds of three different Mach numbers 1.5, 2.0, and 2.5. The proper orthogonal decomposition (POD) then was adopted to investigate the modes of the fluctuation flow field. It emerged that the streamwise and ring-like vortical modes were disparate in energy distribution, structural order, frequency and amplitude. Furthermore, it showed that as the Mach number increased, the energy of the streamwise modes increased while the opposite was true for ring-like modes; and the streamwise modal structures were altered more significantly than the ring-like modes, and the frequency of each mode scarcely varied. It was also found that the streamwise vortices absorbed energy from the ring-like vortices while they traveled from the MVG down to the SWBLI region, but the dominant frequency of each mode rarely changed during this process.more » « less