Efficient Computation of Large-Scale Statistical Solutions to Incompressible Fluid Flows
Article No.: 8, Pages 1 - 11
Abstract
This work presents the development, performance analysis and subsequent optimization of a GPU-based spectral hyperviscosity solver for turbulent flows described by the three dimensional incompressible Navier-Stokes equations. The method solves for the fluid velocity fields directly in Fourier space, eliminating the need to solve a large-scale linear system of equations in order to find the pressure field. Special focus is put on the communication intensive transpose operation required by the fast Fourier transform when using distributed memory parallelism. After multiple iterations of benchmarking and improving the code, the simulation achieves close to optimal performance on the Piz Daint supercomputer cluster, even outperforming the Cray MPI implementation on Piz Daint in its communication routines. This optimal performance enables the computation of large-scale statistical solutions of incompressible fluid flows in three space dimensions.
References
[1]
Utkarsh Ayachit, Andrew Bauer, Berk Geveci, Patrick O'Leary, Kenneth Moreland, Nathan Fabian, and Jeffrey Mauldin. 2015. Paraview catalyst: enabling in situ data analysis and visualization. In Proceedings of the First Workshop on In Situ Infrastructures for Enabling Extreme-Scale Analysis and Visualization (ISAV2015). Association for Computing Machinery, Austin, TX, USA, 25--29. isbn: 9781450340038.
[2]
David L. Brown. 1995. Performance of under-resolved two-dimensional incompressible flow simulations. Journal of Computational Physics, 122, 1, 165--183.
[3]
CSCS. 2023. Piz daint. Retrieved October 04, 2023 from https://www.cscs.ch/computers/piz-daint/.
[4]
Ronald J. DiPerna and Andrew J. Majda. 1987. Oscillations and concentrations in weak solutions of the incompressible fluid equations. Communications in Mathematical Physics, 108, 4, 667--689.
[5]
Ulrich Drepper. 2007. What every programmer should know about memory. Retrieved February 20, 2024 from https://www.akkadia.org/drepper/cpumemory.pdf.
[6]
Truong Vinh Truong Duy and Taisuke Ozaki. 2014. A decomposition method with minimum communication amount for parallelization of multi-dimensional ffts. Computer Physics Communications, 185, 1, 153--164.
[7]
Niklas Fehn, Martin Kronbichler, Peter Munch, and Wolfgang A. Wall. 2022. Numerical evidence of anomalous energy dissipation in incompressible euler flows: towards grid-converged results for the inviscid taylor-green problem. Journal of Fluid Mechanics, 932, A40.
[8]
Ulrik S. Fjordholm, Samuel Lanthaler, and Siddhartha Mishra. 2017. Statistical solutions of hyperbolic conservation laws: foundations. Archive for Rational Mechanics and Analysis, 226, 2, (Nov. 2017), 809--849.
[9]
Ulrik S. Fjordholm, Kjetil O. Lye, Siddhartha Mishra, and Franziska Weber. 2020. Statistical solutions of hyperbolic systems of conservation law: Numerical approximation. Mathematical Models and Methods in Applied Sciences, 30, 3, 539--609.
[10]
Ulrik S. Fjordholm, Siddhartha Mishra, and Franziska Weber. 2022. On the vanishing viscosity limit of statistical solutions of the incompressible navier-stokes equations. arXiv: 2110.04674.
[11]
Ciprian Foias, Oscar Manley, Ricardo Rosa, and Roger Temam. 2001. Navier-Stokes Equations and Turbulence. Encyclopedia of Mathematics and its Applications. Cambridge University Press.
[12]
Matteo Frigo and Steven G. Johnson. 2005. The design and implementation of fftw3. Proceedings of the IEEE, 93, 2, 216--231.
[13]
Uriel Frisch. 1995. Turbulence: The Legacy of A.N. Kolmogorov. Cambridge University Press. isbn: 9780521451031.
[14]
Salman Habib et al. 2016. Hacc: simulating sky surveys on state-of-the-art supercomputing architectures. New Astronomy, 42, 49--65.
[15]
Samuel Lanthaler and Siddhartha Mishra. 2015. Computation of measure-valued solutions for the incompressible Euler equations. Mathematical Models and Methods in Applied Sciences, 25, 11, 2043--2088.
[16]
Samuel Lanthaler and Siddhartha Mishra. 2020. On the convergence of the spectral viscosity method for the two-dimensional incompressible euler equations with rough initial data. Foundations of Computational Mathematics, 20, 5, (Oct. 2020), 1309--1362.
[17]
Samuel Lanthaler, Siddhartha Mishra, and Carlos Parés-Pulido. 2021. Statistical solutions of the incompressible euler equations. Mathematical Models and Methods in Applied Sciences, 31, 02, (Feb. 2021), 223--292.
[18]
Kjetil O. Lye. 2020. Computation of statistical solutions of hyperbolic systems of conservation laws. Ph.D. Dissertation.
[19]
Siddhartha Mishra, Carlos Parés-Pulido, and Kyle G. Pressel. 2020. Arbitrarily high-order (weighted) essentially non-oscillatory finite difference schemes for anelastic flows on staggered meshes. (2020). arXiv: 1905.13665.
[20]
NVIDIA. 2023. Cufft library user's guide. NVIDIA. Retrieved October 04, 2023 from https://docs.nvidia.com/cuda/pdf/CUFFT_Library.pdf.
[21]
Steven A. Orszag. 1971. On the elimination of aliasing in finite-difference schemes by filtering high-wavenumber components. Journal of Atmospheric Sciences, 28, 6, 1074--1074. <1074:OTEOAI>2.0.CO;2.
[22]
Carlos Pares-Pulido. 2021. Statistical solutions of the incompressible Euler equations with finite volume methods. Ph.D. Dissertation.
[23]
Eitan Tadmor. 2004. Burgers' Equation with Vanishing Hyper-Viscosity. Communications in Mathematical Sciences, 2, 2, 317--324.
[24]
Eitan Tadmor. 1989. Convergence of spectral methods for nonlinear conservation laws. SIAM Journal on Numerical Analysis, 26, 1, 30--44.
[25]
2002. From semidiscrete to fully discrete: stability of runge-kutta schemes by the energy method. ii. D.J. Estep and S. Tavener. Collected Lectures on the Preservation of Stability Under Discretization. Chap. 3.
[26]
Geoffrey I. Taylor and Albert E. Green. 1937. Mechanism of the production of small eddies from large ones. Proceedings of the Royal Society of London. Series A - Mathematical and Physical Sciences, 158, 895, 499--521.
Index Terms
- Efficient Computation of Large-Scale Statistical Solutions to Incompressible Fluid Flows
Recommendations
A fast pressure-correction method for incompressible two-fluid flows
We have developed a new pressure-correction method for simulating incompressible two-fluid flows with large density and viscosity ratios. The method's main advantage is that the variable coefficient Poisson equation that arises in solving the ...
Efficient computation of compressible and incompressible flows
The combination of explicit Runge-Kutta time integration with the solution of an implicit system of equations, which in earlier work demonstrated increased efficiency in computing compressible flow on highly stretched meshes, is extended toward ...
Method SMIF for Incompressible Fluid Flows Modeling
NAA 2012: Revised Selected Papers of the 5th International Conference on Numerical Analysis and Its Applications - Volume 8236For solving of the Navier-Stokes equations describing 3D incompressible viscous fluid flows the Splitting on physical factors Method for Incompressible Fluid flows SMIF with hybrid explicit finite difference scheme second-order accuracy in space, ...
Comments
Please enable JavaScript to view thecomments powered by Disqus.Information & Contributors
Information
Published In
June 2024
296 pages
ISBN:9798400706394
DOI:10.1145/3659914
- Proceedings Chairs:
- Timothy Robinson,
- Michèle Woodtli,
- Program Chairs:
- Axel Huebl,
- Cristina Silvano
Copyright © 2024 Copyright is held by the owner/author(s). Publication rights licensed to ACM.
Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than the author(s) must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from [email protected].
Sponsors
Publisher
Association for Computing Machinery
New York, NY, United States
Publication History
Published: 03 June 2024
Check for updates
Author Tags
Qualifiers
- Research-article
Conference
PASC '24
Sponsor:
Acceptance Rates
PASC '24 Paper Acceptance Rate 26 of 36 submissions, 72%;
Overall Acceptance Rate 109 of 221 submissions, 49%
Contributors
Other Metrics
Bibliometrics & Citations
Bibliometrics
Article Metrics
- 0Total Citations
- 55Total Downloads
- Downloads (Last 12 months)55
- Downloads (Last 6 weeks)6
Reflects downloads up to 23 Feb 2025
Other Metrics
Citations
View Options
Login options
Check if you have access through your login credentials or your institution to get full access on this article.
Sign in