More Like this

1. Video: High-Resolution 4D Lagrangian Coherent Structures

   https://doi.org/10.1103/APS.DFD.2022.GFM.V0025  

   Lagares, Christian; Araya, Guillermo (November 2022, 75th APS-DFD November 2022)

   High-speed, spatially-evolving turbulent boundary layers are of great importance across civilian and military applications. Furthermore, compressible boundary layers present additional challenges for energy and active scalar transport. Understanding transport phenomena is critical to efficient high-speed vehicle designs. Although at any instantaneous point in time a flow field may seem random, regions within the flow can exhibit coherency across space and time. These coherent structures play a key role in momentum and energy transport within the boundary layer. The two main categories for coherent structure identification are Eulerian and Lagrangian approaches. In this video, we focus on 4D (3D+Time) Lagrangian Coherent Structure (LCS), and the effect of wall curvature/temperature on these structures. We present the finite-time Lyapunov exponent (FTLE) for three wall thermal conditions (cooling, quasi-adiabatic and heating) for a concave wall curvature that builds on the experimental study by Donovan et al. (J. Fluid Mech., 259, 1-24, 1994). The flow is subject to a strong concave curvature (δ/R ~ -0.083, R is the curvature radius) followed by a very strong convex curvature (δ/R = 0.17). A GPU-accelerated particle simulation forms the basis for the 3-D FTLE where particles are advected over flow fields obtained via Direct Numerical Simulation (DNS) with high spatial/temporal resolution. We also show the cross-correlation between Q2 events (ejections) and the FTLE. The video is available at: https://gfm.aps.org/meetings/dfd-2022/63122e0e199e4c2da9a946a0 

   more » « less
2. A GPU-Accelerated Particle Advection Methodology for 3D Lagrangian Coherent Structures in High-Speed Turbulent Boundary Layers

   https://doi.org/10.3390/en16124800  

   Lagares, Christian; Araya, Guillermo (June 2023, Energies)

   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. 

   more » « less
3. A GPU-Accelerated Particle Advection Methodology for 3D Lagrangian Coherent Structures in High-Speed Turbulent Boundary Layers

   Lagares C. and Araya G. (May 2023, Preprints)

   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 manuscript 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. 

   more » « less
4. Aquila-LCS: GPU/CPU-Accelerated Particle Advection Schemes for Large-Scale Simulations

   Christian Lagares and Guillermo Araya (March 2024, Software X)

   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
5. Characterizing the onset of transitional and turbulent flow regimes in pipe flows using instantaneous time-frequency-based analysis

   https://doi.org/10.1063/5.0226070  

   Shirdade, Nikhil; Kolliyil, Jibin Joy; El-Khader, Baha_Al-Deen T; Brindise, Melissa C (October 2024, Physics of Fluids)

   Accurately identifying the onset of transitional and turbulent flow within any pipe flow environment is of great interest. Most often, the critical Reynolds number (Re) is used to pinpoint the onset of turbulence. However, the critical Re is known to be highly variable, depending on the specifics of the flow system. Thus, for flows (e.g., blood flows), where only one realization (i.e., one mean Re) exists, the presence of transitional and turbulent flow behaviors cannot be accurately determined. In this work, we aim to address this by evaluating the extent to which instantaneous time-frequency (TF)-based analysis of the fluctuating velocity field can be used to evaluate the onset of transitional and turbulent flow regimes. Because current TF analysis methods are not suitable for this, we propose a novel “wavelet-Hilbert time-frequency” (WHTF) method, which we validate herein. Using the WHTF method, we analyzed the instantaneous dominant frequency of three planar particle image velocimetry-captured pipe flows, which included one steady and two pulsatile with Womersley numbers of 4 and 12. For each case, data were captured at Re's spanning 800–4500. The instantaneous dominant frequency analysis of these flows revealed that the magnitude, size, and coherence of two-dimensional spatial frequency structures were uniquely different across flow regimes. Specifically, the transitional regime maintained the most coherent, but lowest magnitude frequency structures, while the laminar regime had the highest magnitude, lowest coherence, and smallest frequency structures. Overall, this study demonstrates the efficacy of TF-based metrics for characterizing the progression of transition and turbulent flow development. 

   more » « less