A fast, matrix-based method to perform omnidirectional pressure integration
Abstract Experimentally-measured pressure fields play an important role in understanding many fluid dynamics problems. Unfortunately, pressure fields are difficult to measure directly with non-invasive, spatially resolved diagnostics, and calculations of pressure from velocity have proven sensitive to error in the data. Omnidirectional line integration methods are usually more accurate and robust to these effects as compared to implicit Poisson equations, but have seen slower uptake due to the higher computational and memory costs, particularly in 3D domains. This paper demonstrates how omnidirectional line integration approaches can be converted to a matrix inversion problem. This novel formulation uses an iterative approach so that the boundary conditions are updated each step, preserving the convergence behavior of omnidirectional schemes while also keeping the computational efficiency of Poisson solvers. This method is implemented in Matlab and also as a GPU-accelerated code in CUDA-C++. The behavior of the new method is demonstrated on 2D and 3D synthetic and experimental data. Three-dimensional grid sizes of up to 125 million grid points are tractable with this method, opening exciting opportunities to perform volumetric pressure field estimation from 3D PIV measurements.
- Sponsoring Organization:
- USDOE
- Grant/Contract Number:
- 89233218CNA000001
- OSTI ID:
- 2319074
- Alternate ID(s):
- OSTI ID: 2315071; OSTI ID: 2350641
- Journal Information:
- Measurement Science and Technology, Journal Name: Measurement Science and Technology Journal Issue: 6 Vol. 35; ISSN 0957-0233
- Publisher:
- IOP PublishingCopyright Statement
- Country of Publication:
- United Kingdom
- Language:
- English
Similar Records
Singular value decomposition utilizing parallel algorithms on graphical processors
One-shot omnidirectional pressure integration through matrix inversion