Pressure spectra from single-snapshot tomographic PIV

The power spectral density and coherence of temporal pressure fluctuations are obtained from low-repetition-rate tomographic PIV measurements. This is achieved by extension of recent single-snapshot pressure evaluation techniques based upon the Taylor’s hypothesis (TH) of frozen turbulence and vortex-in-cell (VIC) simulation. Finite time marching of the measured instantaneous velocity fields is performed using TH and VIC. Pressure is calculated from the resulting velocity time series. Because of the theoretical limitations, the finite time marching can be performed until the measured flow structures are convected out of the measurement volume. This provides a lower limit of resolvable frequency range. An upper limit is given by the spatial resolution of the measurements. Finite time-marching approaches are applied to low-repetition-rate tomographic PIV data of the flow past a straight trailing edge at 10 m/s. Reference results of the power spectral density and coherence are obtained from surface pressure transducers. In addition, the results are compared to state-of-the-art experimental data obtained from time-resolved tomographic PIV performed at 10 kHz. The time-resolved approach suffers from low spatial resolution and limited maximum acquisition frequency because of hardware limitations. Additionally, these approaches strongly depend upon the time kernel length chosen for pressure evaluation. On the other hand, the finite time-marching approaches make use of low-repetition-rate tomographic PIV measurements that offer higher spatial resolution. Consequently, increased accuracy of the power spectral density and coherence of pressure fluctuations are obtained in the high-frequency range, in comparison to the time-resolved measurements. The approaches based on TH and VIC are found to perform similarly in the high-frequency range. At lower frequencies, TH is found to underestimate coherence and intensity of the pressure fluctuations in comparison to time-resolved PIV and the microphone reference data. The VIC-based approach, on the other hand, returns results on the order of the reference.


Introduction
Measurements of the unsteady pressure fluctuations, their time spectrum and their coherence are relevant in aerodynamics, aero-elasticity and aero-acoustics, since they provide information for estimation of the unsteady loads and noise sources. In the past two decades, pressure from particle image velocimetry (PIV) has received significant attention from the scientific community. It is an alternative to wind tunnel model instrumentation with pressure transducers, or pressure-sensitive paint (PSP;Bitter et al. 2012), which has limited sensitivity in the low-speed regime. Many research groups have focused their efforts on all parts of the pressurefrom-PIV measurement chain, as summarized in Table 1. Starting from the measurement techniques (Table 1a), researchers have considered time-resolved high-repetitionrate PIV measurements, but also multi-pulse and dual-pulse (resp. single-snapshot) approaches. The measurement data is subsequently processed using either a cross-correlation or Lagrangian particle tracking (LPT) approach (Table 1b). All existing approaches obtain the pressure gradient from the velocity measurements by invoking the momentum equation. Depending on the measurement setup and data analysis procedure, many techniques to obtain the pressure gradient have been proposed (Table 1c). Finally, also the subsequent integration of the pressure gradient remains a topic of current research (Table 1d). In Table 1e salient comparative assessments that compare part of the above procedures are listed. Note that the majority of these works use volumetric measurement data, following recent research that shows this 1 3 57 Page 2 of 15 is a requirement to accurately obtain the pressure gradient (Charonko et al. 2010;de Kat and van Oudheusden 2012;Ghaemi et al. 2012). It should be remarked that some works have considered pressure evaluation in compressible flows (e.g., van Oudheusden et al. 2007;Ragni et al. 2009), but the present discussion and work is limited to incompressible, isothermal flows and Newtonian fluids.
Despite the above developments, recent works have not been able to measure pressure spectra and coherence at high frequencies (de Kat and van Oudheusden 2012; Ghaemi et al. 2012;Pröbsting et al. 2015). Even though these studies used time-resolved tomographic PIV at frequencies up to 10 kHz, the pressure spectra were found to be dominated by measurement errors well before the Nyquist frequency. This is a direct consequence of the limited spatial resolution of available time-resolved pressurefrom-PIV techniques (see e.g. Charonko et al. 2010;de Kat and van Oudheusden 2012;McClure and Yarusevych 2017). The largest resolvable frequency of pressure fluctuations that can be estimated depends on the smallest resolvable scale. The latter is approximately twice the interrogation volume (IV) size used for PIV cross-correlation analysis (Schrijer and Scarano 2008). Accordingly, the largest frequency of pressure fluctuations that can be measured is given by: where u c,p is the convective speed of the pressure fluctuations. For example, in case the turbulent boundary layer measurements by Pröbsting et al. (2015), a window size of 1 mm was achieved. This resulted in a maximum resolved frequency of approximately 3 kHz, despite the use of a timeresolved measurement system operating at 10 kHz. The comparative assessment in van Gent et al. (2017) shows that a better result (i.e., higher spatial resolution) could be obtained using Lagrangian particle tracking (Table 1b). However, due to limited measurement rates, Lagrangian particle tracking (LPT) remains limited to relatively low-speed flows to allow accurate particle pairing and tracking. Multi-pulse approaches (Table 1a) overcome these measurement rate limitations. However, these techniques return uncorrelated velocity and pressure fields that do not directly allow for evaluation of the spectrum and moreover require relatively expensive setups with an extended number of imagers and lasers.
Avoiding the need for additional hardware, a similar result can be obtained using single-snapshot approaches, that have been introduced to approximate instantaneous pressure from standard low-repetition-rate tomographic PIV (de Kat and Ganapathisubramani 2012;Laskari et al. 2016;. These techniques invoke modeling approaches to approximate the velocity material derivative that otherwise would be measured by time-resolved or multi-pulse approaches. The comparative assessment by van Gent et al. (2017) shows that single-snapshot approaches  Gurka et al. (1999); Hosokawa et al. (2003); de Kat and Ganapathisubramani (2012) (2017) are less accurate than the time-resolved ones. This conclusion, however, does not take into account that singlesnapshot approaches can make use of a low-repetition-rate tomographic PIV system, which can benefit from increased measurement dynamic range due to the availability of larger format CCD cameras (see e.g. Hain et al. 2008;Schanz et al. 2012) and more powerful lasers with respect to the ones used for time-resolved measurements. This allows to employ a smaller interrogation volume size for cross-correlation analysis. However, the established single-snapshot approaches (Table 1a) return only uncorrelated pressure fields, and therefore, require an additional modeling assumption to obtain pressure spectra and coherence. Due to the increased spatial resolution with respect to a time-resolved approach, the result could, however, have higher accuracy.
To extract temporal pressure spectra and coherence from the velocity snapshots, a logical first step is to invoke Taylor's hypothesis of frozen turbulence. This was already done for single-snapshot pressure evaluation (de Kat and Ganapathisubramani 2012; Laskari et al. 2016). In addition, Scarano and Moore (2011) proposed a local advection-only approximation of the Navier-Stokes equations to approximate the temporal velocity fluctuations between correlated PIV measurement time instants, yielding an increased temporal resolution. This was called as time-supersampling. Schneiders et al. (2014) later generalized this technique by avoiding the linearization and use the incompressible vorticity transport equation as model for time-supersampling.
It should be remarked that the problem of supersampling measurement data has also been approached with more complex methods that take further into account the non-linear dynamics in the temporal evolution of the velocity field. In this respect, several methods have been proposed that use a variational approach to assimilate data from measurements and numerical simulations (Suzuki et al. 2009;Cuzol and Mémin 2009;Lemke and Sesterhenn 2016;Yegavian et al. 2015;Gesemann et al. 2016). Such variational approaches have proven to be successful in increasing the spatial resolution of three-dimensional time-resolved measurements Gesemann et al. 2016;Schneiders et al. 2017). However, the literature shows no successful application of variational data assimilation to increase the temporal resolution of volumetric flow measurements.
Single-snapshot pressure evaluation techniques (de Kat and Ganapathisubramani 2012; Laskari et al. 2016;) so far have not included a time-marching scheme for velocity and pressure, and therefore, do not allow for spectral information to be obtained when applied to lowrepetition-rate PIV data. Combination of the single-snapshot pressure evaluation methods and the abovementioned timesupersampling techniques could allow for evaluation of short time segments of the temporal pressure fluctuations from uncorrelated tomographic PIV measurement data. This is different from, for example, , which only evaluates the instantaneous velocity material derivative from the approximated vorticity temporal derivative. The instantaneous velocity material derivative allows calculation of instantaneous pressure, but the scheme used by  does not include time marching and no pressure time series is obtained when the technique is applied to instantaneous data.
The current work, instead, investigates the feasibility to obtain finite time segments of the pressure from low-repetition-rate volumetric velocity measurements using Taylor's hypothesis and a non-linear technique based on solution of the vorticity transport equation. The approaches are assessed for low-repetition-rate (less than 1 Hz) tomographic PIV measurements in the turbulent boundary layer near the trailing edge of a flat plate. The proposed methods are used to obtain the pressure spectrum in the kilo-Hertz range. For validation, the results are compared to a reference spectrum obtained by surface pressure transducer. In addition, the methods are compared to state-of-the-art time-resolved tomographic PIV results.

Finite time marching
To obtain a finite time segment of the velocity and pressure temporal evolution from a single velocity field, information from the spatial domain can be used to increase temporal resolution. The measured velocity field (Fig. 2a, blue volume) is marched forward in time over a finite time segment. At time t + ΔT the volume has marched downstream over a distance of approximately u c,u ΔT ( Fig. 2b, orange volume), where u c,u is the convective speed of the velocity fluctuations. Accordingly, erosion of the measurement region starting from the inflow boundary is expected.
Finite time marching of the measured velocity field can be performed using the Taylor's hypothesis of frozen turbulence as follows  where u c,u is the local convective speed of velocity. This speed is not generally known and needs to be selected. Laskari et al. (2016) compared several choices for determining the appropriate convective velocity for single-snapshot pressure evaluation in case of a turbulent boundary layer. It was concluded that the use of a streamwise mean of the velocity field yields most accurate results.
Alternatively, the work of Schneiders et al. (2014) provides a generalization of the above procedure to flows where the convective speed cannot be estimated, by performing a time-marching simulation of the incompressible vorticity transport equation using the vortex-in-cell technique (VIC, Christiansen 1973): where in the present study the viscous term is included. Within the VIC framework, velocity is related to vorticity through: Boundary conditions for the time-marching simulation are obtained by finite-difference solution of the vorticity transport equation. Boundary values of velocity are set following the padding procedure used in . The VIC technique in Schneiders et al. (2014) requires imposing boundary values of vorticity, which are obtained by finitedifference solution of the vorticity transport equation. The measurement volume is padded by a buffer region (gray area in Fig. 1), where the streamwise average velocity is initially prescribed. During the time-marching simulation, the information in the measurement volume region convects into this buffer region as illustrated in Fig. 1. The extent L p of domain padding is determined by the integration time, T, and a typical velocity value representative of flow convection: Subsequently, the pressure temporal evolution in each time segment is obtained using an established pressurefrom-PIV approach that solves the Poisson equation for pressure (Table 1d). Boundary conditions on the pressure gradient (von Neumann type) are obtained from the momentum equation and are prescribed on all boundaries except one (i.e., top boundary in Fig. 1), where Dirichlet-type boundary conditions are assigned. The latter can be assigned as the free-stream mean pressure level, if part of the measurement volume extends up to the free stream, or it can be approximated from a model such as the Bernoulli equation or isentropic relations (Kurtulus et al. 2007;Ragni et al. 2009; de Kat and van Oudheusden 2012; among others).
The velocity material derivative, Du/Dt, is calculated in Eulerian formulation (Table 1c).
where the temporal derivative is obtained using a secondorder central difference scheme from the velocity time series. The Eulerian formulation can be used because the methods allow for selection of a sufficiently small time step Δt for time-marching to allow Eulerian derivative evaluation.
The statistical convergence of the pressure spectrum estimation is achieved by averaging the spectra obtained from time-segments from multiple uncorrelated tomographic PIV snapshots. This is analogous to Welch's power spectral density estimation procedure (Welch 1967). The short time segments require a windowing function. To avoid a discontinuity at the edges, a Hann window is used.

Dynamic frequency range
The theoretical maximum duration of a time segment, T max , corresponds to the time when all the measured fluid parcels are convected out of the measurement volume. This can be estimated as the ratio of the streamwise measurement domain length, L x , and a typical convective velocity in the flow, u c , Accordingly, the lowest resolvable frequency is It should be remarked that in a shear flow, different parts of the volume are convected out of the measurement in different moments. In such a case, the maximum velocity is to be taken as characteristic velocity. Additionally, numerical instabilities or end effects may reduce the maximum time-segment time. As a consequence, the above equations serve only as an indication of the achievable f min .
The highest resolvable frequency is limited by the available spatial resolution. An estimate can be obtained considering the ratio of the smallest resolved length scale, λ min , and a typical convective velocity, Considering spatial cross-correlation for the analysis of the particle images, the smallest resolved length scale of velocity fluctuations is a function of the size of the interrogation volume, IV. According to Schrijer and Scarano (2008), a spatial modulation up to 50% of the velocity fluctuation is achieved when λ = 2IV. Equation (3) can, therefore, be rewritten as: In conclusion, the theoretical frequency dynamic range (DFR), given as the ratio of the highest to lowest resolvable frequency, reads as Interestingly, the above expression corresponds to the dynamic spatial range, DSR, introduced by Adrian (1997) multiplied by a coefficient. The latter is a direct consequence of the fact that the proposed technique evaluates the pressure spectrum from measurement data distributed in the spatial domain only.

Modeling limitations
In addition to limitations on the achievable frequency range (Sect. 3.1), the theoretical limitations of single-snapshot pressure from PIV modeling approaches also hold for the present work . It was shown that only pressure fluctuations caused by measured flow structures are reconstructed. External effects caused by, for example, pulsating jets or a flapping wing are not predicted by the modeling approach, and therefore, the corresponding pressure fluctuations are not measured.

Description of experiments
The approaches discussed in Sect. 2 are applied to a tomographic PIV measurement at the trailing edge of a flat plate (Fig. 3) installed in a low-speed wind tunnel of 40 × 40 cm 2 cross-section operating at free-stream velocity of 10 m/s. The flat plate is 60 cm long, 40 cm wide and is 1 cm thick. The flat plate thickness decreases to 0.2 mm at the trailing edge with an angle of 2.4° on both sides, which creates a mild pressure gradient. The boundary layer is tripped 10 cm downstream an elliptical leading edge using a 1 cm strip of 3D carborundum roughness elements of 0.8 mm diameter. The boundary layer parameters were measured by planar PIV. The boundary layer thickness, δ 99 , equals approximately 11.9 mm at x = − 15 mm (Fig. 2). Other relevant boundarylayer parameters are summarized in Table 2, where the wall shear velocity was estimated by fitting with the law of the wall with κ = 0.41 and A 2 = 5.0 following Pröbsting et al. (2016). The values are in good agreement with those reported by Pröbsting et al. (2015) in a similar experiment where measurements were performed using high-repetitionrate tomographic PIV. The results from the latter study are used for comparison to those obtained by the present lowrepetition-rate experiment.  Fig. 2, captures 30 mm in streamwise direction, 10 mm in spanwise direction and 5.6 mm in wall-normal direction, which corresponds to 50 percent of the boundary layer thickness. The average magnification equals M = 0.38 and maximum particle displacement equals 15 voxels. A summary of the relevant experiment parameters is given in Table 3. A series of 50 particle image pairs is reconstructed using the MART technique (Elsinga et al. 2006). Crosscorrelation is done using a multi-pass volume deformation (Scarano and Riethmuller 2000) algorithm with final interrogation volume size equal to 40 × 40 × 40 voxels (0.78 × 0.78 × 0.78 mm 3 ) and with 75% overlap. The resulting vector spacing equals 0.19 mm. Universal outlier detection (Westerweel and Scarano 2005) is applied to the vector fields to remove spurious vectors. The origin of the axis system is positioned at the plate-trailing edge (Fig. 2).

Microphone surface pressure measurements
A Sonion 8010T microphone is mounted below a pinhole positioned 30 mm upstream of the trailing edge (red dot in Fig. 2), to measure surface pressure fluctuations to be compared with those extracted from tomographic PIV. Figure 3 shows a schematic of the pinhole-microphone arrangement. The microphone has a flat frequency response in the range of 0.25-11.5 kHz and the pinhole arrangement is designed such that the resonance frequency of the cavity is beyond this range (i.e., 12 kHz, Alster 1972). The microphone signal is recorded over 90 s at 51.2 kHz using an NI 9234 signal acquisition module that has a built-in anti-aliasing module. The microphone measurements are carried out after the PIV measurements. A frequency-dependent calibration (Zuckerwar et al. 2006) was carried out using a LinearX M53 reference microphone, which was calibrated by means of a GRAS 42AA piston phone. No further signal processing is applied. Spectral analysis is performed using a periodogram method with Hamming windows of 80 ms, resulting in a frequency sampling of 12.5 Hz.

Results and discussion
Two consecutive low-repetition-rate tomographic PIV velocity measurements are plotted in Fig. 4. The figure shows isosurfaces of high-and low-speed streaks (blue, u −ū = u � = − 1.25 m/s; green, u ′ = + 1.25 m/s) and vortices identified by the Q-Criterion (orange, Q = 10 7 /s 2 ). The free-stream flow is directed along the x-axis. A typical turbulent boundary layer flow organization can be seen, with vortices positioned on top of the low-speed streaks. Two vortical structures, S1 and S2, are indicated in the figure by black arrows. At the measurement rate of 0.4 Hz, there is 2.5 s in between the two consecutive velocity measurements. Taking the free-stream velocity as convective speed, the vortices are expected to travel approximately 25 m in this time corresponding to 400 measurement volume lengths. The actual convective speed in the turbulent boundary layer is expected to be lower, however. Taking the friction velocity as convective speed, structures are expected to travel over approximately 1.1 m, corresponding to 30 measurement volume lengths. Therefore, at t = 2500 ms, the structures S1 and S2 are not present in the measurement volume anymore.

Finite time marching
The measurement at time t = 0 ms ( Fig. 4-left) is animated using both Taylor's hypothesis and the vorticity equation.
For the VIC-based procedure, a time step of 10 µs is required (Schneiders et al. 2014). This corresponds to a measurement frequency of 100 kHz. However, based on the interrogation volume size and free-stream velocity, the expected maximum resolved frequency (Eq. 9) of velocity fluctuations equals 6.4 kHz. Despite the small integration time step, it is, therefore, not expected that fluctuations beyond 6.4 kHz are resolved. To reduce memory requirements, the results are stored at a reduced sampling rate of 33 kHz, which is still well beyond the expected maximum resolved frequency. Time marching is performed over 1.5 ms, which equals half T max (Eq. 6). The minimum resolved frequency equals accordingly approximately 0.667 kHz (Eq. 7). Note that, the VIC-based procedure inherently applies a divergencefree filter to the measurement data. The recent literature has shown that this allows for important regularization of the measurement data (de Silva et al. 2013;Schiavazzi et al. 2014;Wang et al. 2016a, b). Therefore, the same is applied before using the Taylor's hypothesis to allow direct comparison of the time-marching techniques. The computational time of the VIC-based approach is on the order of the computational time required for processing of a similar time-resolved measurement dataset. This is approximately an order of magnitude larger than the computational cost of application of Taylor hypothesis. The computational cost of the divergence-free filter is negligible (order of seconds per snapshot) in comparison to the other processing steps in the measurement chain. The structures S1 and S2 are indicated in the perspective view, and have marched downstream for both the Taylor's hypothesis and VIC results. The side-view visualizations of both results show an empty region without vortical structures in the upstream flow region for t > 0, as indicated by the dashed line. This eroded flow region is growing with time marching because the finite time marching techniques only leverage the measured flow information and do not provide a boundary condition for new structures entering into the domain, as discussed in Sect. 2.
A closer inspection to the perspective view at t = 0.75 ms shows less vortical structures in the near-wall region for the VIC result in comparison to the result from Taylor's hypothesis. The side views indicate that increased coherence of the vortical structures is obtained by the VIC-based result. In addition, the VIC-based approach employs side-padding of the measurement volume (Fig. 1), which results in marginally increased coherence of the vortical structures on the domain boundaries as visible in Fig. 5. The video-animations provided as supplementary material further illustrate the development of the near-wall vortical structures.
The coherence of velocity fluctuations is calculated and compared to the time-resolved tomographic PIV results obtained by Pröbsting et al. (2015). The results in the domain center along the full spanwise extent of the measurement volume are taken to increase convergence of the coherence estimation. The results are evaluated at the same wall-normal location (y = 3 mm) as was used in Pröbsting et al. (2015). Note that, by Pröbsting et al. (2015), the highrepetition-rate results were evaluated 30 mm downstream at the location of the trailing edge (Fig. 2). In Fig. 6, the spanwise coherence of wall-normal velocity fluctuations at 2 kHz (left) and 5 kHz (right) is plotted. Pröbsting et al. (2015) found the spatial resolution of time-resolved tomographic PIV (window size indicated in Fig. 6) to be sufficient to capture fluctuations at 2 kHz, providing a reference for the current results.
The results from the low-repetition-rate system (blue and red lines) are found to be in good correspondence with the TR-PIV results. At the higher frequency of 5 kHz (right figure), on the Nyquist limit of the time-resolved system, Pröbsting et al. (2015) reported that the coherence is overestimated due to spatial filtering. Because the low-repetition-rate system uses a smaller correlation volume, a sharper coherence function than obtained by the time-resolved system is expected. This is confirmed by the results in Fig. 6-right. The result from Taylor's hypothesis and VIC-based finite time marching yields similar coherence, with the VIC-based result being marginally higher. Both techniques yield a coherence of 0.2 at approximately 25% smaller distance than the high-repetition-rate result.

Pressure time evolution from a snapshot
Pressure is calculated from the velocity temporal evolution in the time segments. Figure 7 shows the resulting isosurfaces of pressure (colored by z, p = − 0.5 Pa) in a visualization similar to Fig. 6. In both the TH (left) and VIC (right) results, low-pressure regions can be seen to coincide with the vortical structures S1 and S2 highlighted earlier. First comparison of the TH and VIC results shows that the latter yields a smoother pressure field. In addition, the VIC-based result shows a large low-pressure region near the domain The results in Fig. 7 show instantaneous pressure fields at the measurement time-instant (t 0 = 0) and the results obtained from the time-marching approaches at two later time instants (t = 0.75 ms and t = 1.5 ms). The latter two correspond to time-marching over ΔT = t 0 − t = 0.75 and 1.5 ms. The root mean squared (rms) surface pressure fluctuations can be evaluated at each ΔT by considering the data ensemble of all low-repetition-rate measurements. Inspection of the rms surface pressure fluctuations calculated at each time instant at the wall (y = 0) shows that the VIC-based result settles to a constant rms level of 0.6 Pa after 0.1 ms time marching (Fig. 8). This level compares well to the rms pressure fluctuation level obtained from the microphone measurements, that were band-pass filtered in the region where they return a flat response . This compares well to the rms pressure level of 0.57 Pa that was measured in a similar turbulent boundary layer by Ghaemi et al. (2012). In addition, the time-resolved result by Pröbsting et al. (2015) using N = 5 returns this rms value. The value N indicates the number of consecutive snapshots used for evaluation of the velocity material derivative. Increasing this number corresponds to stronger filtering and thereby returns a lower rms level. The result obtained by Taylor's hypothesis yields the lowest rms pressure level (initially 0.46 Pa), but does not suffer from the initial high rms peak value returned by the VIC analysis. In addition, at approximately 1.5 ms, the VICbased result returns again an increased rms pressure value.

Pressure coherence
The magnitude of the spanwise coherence of surface (y = 0) pressure fluctuations is plotted in Fig. 9. This represents the cross-power spectral density of pressure fluctuations in two points at a spanwise separation Δz . Again, for calculation of the coherence, the results in the domain center along the full spanwise extent of the measurement volume are taken to increase convergence of the coherence estimation. In addition, again the high-repetition-rate results (Pröbsting et al. 2015) were evaluated 30 mm downstream at the location of the trailing edge (Fig. 2). The time-resolved PIV result (left figure) shows data up to a frequency of 5000 Hz, which equals half the sampling frequency. The coherence is found to be highest in the low-frequency range and decreasing for increased frequencies. The low-repetition-rate results (middle and right figures) provide data up to higher frequencies, but lack the low-frequency range (dashed in the figures) due to the limited time segments that can be obtained from instantaneous velocity snapshots (Sect. 3.1). The results are similar to those obtained by the time-resolved approach, but a significantly smaller coherence is found in the high frequency range.
The lower coherence at high frequency is attributed to the smaller interrogation volume size that is achieved using the low-repetition-rate measurement system. To investigate the  Pröbsting et al. (2015, green lines) are indicated in the figure   Fig. 9 Magnitude of the spanwise coherence of surface (y = 0) pressure fluctuations as obtained from time-resolved tomographic PIV by Pröbsting et al. (2015, left figure) and low-repetition-rate PIV using Taylor's hypothesis (middle) and VIC (right) effect of the interrogation volume size on the results, the data is also processed using a 50 and 100% larger interrogation volume size. In Fig. 10 the resulting coherence obtained using the VIC-based approach is given. Taking a threshold of 0.2 as the coherence length, a minimum coherence length is found at approximately l min = 3Δz , corresponding to 75% of the interrogation volume size. This is indicated by the dashed red line in the figure. The vertical dashed black line indicates the corresponding expected maximum resolved frequency calculated from Eq. (8). As expected, the minimum coherence length and the maximum resolved frequency increases when interrogation volume size is reduced. When the interrogation volume size is reduced to 0.8 mm, noise is found to start dominating the high-frequency result, indicating that the interrogation volume size cannot be reduced further. The latter interrogation volume size is used in the remainder of this work. Pröbsting et al. (2015) computed the correlation length from the results by means of an exponential fit. These results are reproduced in Fig. 11. The results from the models by Corcos (1964) and Efimtsov (1982) that were used for comparison by Pröbsting et al. (2015) are also reproduced in the figure. The result from high-repetition-rate PIV is found to correspond well to the models around 1 kHz, but flattens for higher frequencies.
In the high-frequency range, the low-repetition-rate results using both Taylor's hypothesis and VIC show good correspondence with the Corcos and Efimtsov models, as expected from the fact that these methods benefit from increased spatial resolution. The approach using Taylor's hypothesis shows a lower coherence length than the VIC result for frequencies below 2 kHz. The VIC result follows the Efimtsov model and returns the same coherence length as high-repetition-rate PIV in the range of 0.7-1 kHz. Reduced accuracy of the result of Taylor's hypothesis at low frequencies could be expected, as the assumption of frozen turbulence is not valid for large structures corresponding to low frequencies. The non-linear VIC-based model appears to improve upon this, as also observed by Schneiders et al. (2014).

Pressure spectrum
The pressure spectrum calculated from the microphone measurements is represented with a gray line in Fig. 12. The resonance frequency of the pinhole microphone system was estimated at 12 kHz (Sect. 3.2), which is close to the peak region in the present experiments. The high-repetitionrate tomographic PIV experiments performed by Pröbsting et al. (2015) at a measurement rate of 10 kHz yielded pressure fluctuations up to 5 kHz. These are shown in Fig. 8 by the blue markers. Note that, the high-repetition-rate results were evaluated at the trailing edge location, whereas the microphone location is 30 mm upstream of the trailing edge (Fig. 2).
The high-repetition-rate results showed a significant dependence upon the length of the temporal kernel used  to estimate the material derivative (N = 3, 5, 7 and 9 snapshots). The pressure spectrum estimated with a relatively short kernel overestimates the pressure fluctuation amplitude in the high-frequency range, which is ascribed to the effect of measurement noise. The latter is plausible causing the noise floors evident at higher frequencies. A longer kernel does filter out the noise but at the cost of damping part of the fluctuations due to physical phenomena. This is particularly evident for N = 9. With a kernel of N = 5 or 7 time-resolved tomographic PIV results follow the reference microphone data reasonably well up to a frequency of 3.5 kHz. Reliable results are, therefore, obtained up to a frequency 30% below the Nyquist frequency (5 kHz), which is ascribed to the limited spatial resolution of the high-repetition-rate PIV measurement system. The work of Pröbsting et al. (2015) concluded that the results using kernel sizes N = 5 and 7 are most accurate. Accordingly, these should be considered most accurate for comparison with the current low-repetition-rate approaches.
The pressure spectrum is also evaluated at the microphone location (x = − 30 mm, Fig. 2) for the approaches using low-repetition-rate PIV. As expected from the discussion on the coherence length, the Taylor's hypothesis (blue line) approach shows acceptable results in the high-frequency range. In the low-frequency range, an underestimation of the pressure is found. The VIC-based result follows the microphone measurements with similar accuracy as the high-repetition-rate tomographic PIV data with kernel sizes N = 5 and 7. Beyond 3.5 kHz, the VIC result maintains a good agreement with the microphone data. Based on Eq. (9) the maximum resolvable frequency was estimated around 6.4 kHz. A departure of the microphone results is found at approximately 5 kHz. This is attributed to the convective velocity of the pressure fluctuations being typically lower than the free-stream velocity (Panton and Linebarger 1974).
A more quantitative inspection of the results is possible by considering the difference between the reference spectrum obtained from the transducer and the PIV results. In Fig. 13 this difference is plotted for the relevant range of frequencies. As expected from the discussion above, the highrepetition-rate tomographic PIV result shows a noise floor starting from approximately 3 kHz. The smallest resolved power equals 10 − 6 Pa 2 /Hz, as indicated by the dotted line in Fig. 13. The Taylor's hypothesis and VIC-based results that use low-repetition rate measurement results show a more than an order of magnitude smaller noise floor at 5 × 10 −8 Pa 2 /Hz. In the lower frequency range (1-2 kHz) the highrepetition rate and low-repetition rate VIC results show similar error levels on the order of 20% of the reference. As expected from the discussion above, the Taylor's hypothesis results shows increased errors of around 50% of the reference.

Conclusions
The work shows that the power spectral density and coherence of temporal pressure fluctuations can be obtained from low-repetition-rate tomographic PIV measurements. Theoretical analysis shows that spatial resolution (i.e., interrogation volume size) is inversely proportional to the maximum resolved frequency. Higher spatial resolution of low-repetition-rate measurements in comparison to high-repetition-rate measurements is, therefore, expected to allow for a linear increase in maximum resolved frequency. Spectral information is obtained from the low-repetition-rate measurements Fig. 12 Pressure power spectral density at x = − 30 mm. The dashed lines indicate the expected minimum and maximum resolvable frequencies with finite time-marching analysis. The time-resolved PIV data are taken at x = 0 mm and reproduced from a previous study in a similar flow configuration (Pröbsting et al. 2015) Page 13 of 15 57 by finite time marching using Taylor's hypothesis or a VICbased approach. The experimental assessment shows that both the approach based on Taylor's hypothesis and VIC allow for acceptable results in the high-frequency range. The noise floor of the high-repetition-rate result is found at 10 −6 Pa 2 /Hz. Instead, the low-repetition-rate techniques show a lower noise floor at 5 × 10 −8 Pa 2 /Hz. For lower frequencies (below 2 kHz) corresponding to larger structures, the assumption of frozen turbulence does not hold and an underestimation of the pressure fluctuations is found by the results based on Taylor's hypothesis (on the order of 50% error). The VIC-based approach is able to obtain spectral information in this range with similar accuracy as a highrepetition-rate measurement system (20% error). Fig. 13 Difference between the reference pressure spectrum obtained from the microphone measurements and the PIV results. The time-resolved PIV data is reproduced from a previous study in a similar flow configuration (Pröbsting et al. 2015)