Scanning stereo-PLIF method for free surface measurements in large 3D domains

In this work, we extend a planar laser-induced fluorescence method for free surface measurements to a three-dimensional domain using a stereo-camera system, a scanning light sheet, and a modified self-calibration procedure. The stereo-camera set-up enables a versatile measurement domain with self-calibration, improved accuracy, and redundancy (e.g., possibility to overcome occlusions). Fluid properties are not significantly altered by the fluorescent dye, which results in a non-intrusive measurement technique. The technique is validated by determining the free surface of a hydraulic flow over an obstacle and circular waves generated after droplet impact. Free surface waves can be accurately determined over a height of L=100\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L=100$$\end{document} mm in a large two-dimensional domain (y(x,z)=120×62\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y(x,z) = 120\times 62$$\end{document} mm2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^2$$\end{document}), with sufficient accuracy to determine small amplitude variations (η≈0.2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \approx 0.2$$\end{document} mm). The temporal resolution (Δt=19\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varDelta }t = 19$$\end{document} ms) is only limited by the available scanning equipment (f=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f = 1$$\end{document} kHz rate). For other applications, this domain can be scaled as needed.


Introduction
Small-scale free surface dynamics play a significant role in many applications. For example, they strongly influence the response of structures to wave impacts (Lafeber et al. 2012), and the transfer of heat, momentum, mass, and energy between the ocean and atmosphere (Buckley and Veron 2016;Jähne and Haußecker 1998). Therefore, measurements of the small-scale free surface dynamics are required. These measurements are typically performed with intrusive point measurements (e.g., resistive wave probes).
However, instantaneous wave height measurements can be obtained with optical techniques that are non-intrusive, and therefore preferential in many applications. These optical techniques can be categorized as refraction, stereo-correlation or projection based (Fig. 1).
Techniques based on refraction relate the slope change of a liquid free surface to the refraction of light. A one-dimensional collimated light source (e.g., laser beam) is typically used to determine the deflection on position-sensing photodiodes (Liu et al. 1993) or a camera. A collection of points can be obtained by rapidly scanning the light source over a line (Savelsberg et al. 2006). The technique can be extended to a two-dimensional domain by imaging a pattern through a refractive medium. The patterns are either created by a random-dot pattern (Moisy et al. 2009) or a distribution of particles in a flow (Gomit et al. 2013). Advantages of refraction-based techniques are the ease-of-implementation, the relatively large domain size, and the small error (Table 1). However, ray crossing, which occurs when either the surface curvature or the surface-to-pattern distance becomes too large, limits the measurement of steep and large amplitude waves (Moisy et al. 2009).
Stereo-correlation techniques either use stereo-vision to triangulate characteristic features, such as ocean whitecaps (Benetazzo et al. 2012), or use particle image velocimetry (PIV)/particle tracking velocimetry (PTV) approaches to Fig. 1 General measurement principle of free surface measurement techniques. a Refraction-based techniques relate the slope change to a change in height. b Stereo-correlationbased techniques determine the three-dimensional world position of particles or naturally present features. c Projectionbased techniques determine the free surface height with a projected pattern Table 1 Overview of selected free surface measurement techniques Several parameters are estimated from data presented in each paper, which is denoted by an asterisk ( * ) A parameter that was not defined and could not be estimated is denoted by (-)

Refraction-based techniques
Laser beam refraction Savelsberg et al. (2006) 0 O(1-10) * Duncan et al. (1999) 0.12 500 O(50) * O(1) * determine the free surface height. The particles used in these approaches are often buoyant (Turney et al. 2009) or neutrally buoyant (Douxchamps et al. 2005). The effect of particles on surface tension can be minimized by ensuring that the particles are hydrophobic. However, particle clustering can occur in regions of convergence of interfacial velocity, which can result in local variance of fluid properties (Turney et al. 2009). The stereo-correlation techniques can obtain free surface heights over large spatial domains (Table 1). However, considerable effort is required to obtain appropriate illumination. Furthermore, the sensitivity to small amplitude free surface disturbances is limited (Zavadsky et al. 2017). Lastly, projection-based techniques project an image or pattern on a liquid free surface. No surface seeding is required when there are naturally present contaminants (Gomit et al. 2015) or when the reflection of colored light is used (Dabiri and Gharib 2001). Otherwise, a dye is required to make the liquid fully opaque (Cobelli et al. 2009;Tsubaki and Fujita 2005) or to make the liquid fluorescent (André and Bardet 2014;Buckley and Veron 2017;Duncan et al. 1999). For fully opaque fluids, fringe projection techniques can be applied to obtain accurate, two-dimensional free surface height measurements in large three-dimensional domains (Cobelli et al. 2009). For fluorescent fluids, accurate line measurements can easily be obtained with planar laser-induced fluorescence (PLIF) techniques (Duncan et al. 1999). Furthermore, the fluorescent dye does not change the fluid properties as the concentration of dye is typically low.
So far, PLIF-based free surface measurement techniques have been limited to line measurements. The technique offers accurate free surface height measurements, without altering the fluid properties or the need for fixed patterns. Furthermore, per-pixel wave height measurements can be obtained along a line (Buckley and Veron 2017). In contrast, stereo-correlation-based approaches use particles or features and are limited by the distribution of these particles or features over the surface. The accuracy of PLIF-based techniques can even be improved with detection methods relying on light sheet properties (e.g., André and Bardet 2014). However, the PLIF-based approaches are currently limited to planar domains.
This study presents a two-dimensional PLIF-based approach, which applies a stereo-camera set-up with a scanned light sheet (Brücker 1996). Conventional free surface measurements approaches are summarized in Table 1. The current technique can be used in experiments with limited optical access, due to the small camera separation angle. A multi-step calibration procedure is used to reduce the calibration requirements (Hori and Sakakibara 2004;Wieneke 2008). These optical limitations are imposed by future work in the newly build Multiphase Wave Lab at MARIN in The Netherlands. The technique is also evaluated for typical optical conditions. Measurements can be obtained in relatively large three-dimensional domains ( (x, y, z) = (120, 100, 62) mm) with sufficient accuracy to resolve small amplitude waves ( ≈ 0.2 mm).
The paper is organized as follows. The experimental setup and equipment are introduced in Sect. 2. This section also describes the calibration procedure. Section 3 shows and discusses the measurements of a hydraulic flow over an obstacle and the waves generated by a droplet impact are shown. Finally, the last section presents some concluding remarks.

Method
The method allows one to extend one-dimensional height measurements (y(x, t)) typically obtained with PLIF (Table 1) to a two-dimensional domain (y(x, z, t)). In essence, this could be achieved with a single camera and a multi-plane calibration procedure. However, a stereo-camera set-up allows one to apply a broad calibration that can be refined with an adapted self-calibration procedure as the z-direction is constrained by the scanned light sheet. Furthermore, the accuracy of the method can be improved with a stereo-camera set-up. This results in a robust method that can, with some adjustments, be used in applications with occlusions.
The experimental set-up and equipment are introduced in this section. Furthermore, the method to obtain stereo-PLIF measurements is detailed. First, the relevant equipment is introduced. Then, the edge detection procedure is introduced, which is required for the adapted self-calibration procedure and measurements. Finally, the stereo-PLIF calibration is detailed, which uses an adapted self-calibration procedure. Figure 2 shows the experimental set-up used in this study. The experiments are performed in the water tunnel of the Laboratory for Aero-and Hydrodynamics at the Delft University of Technology. The test section has an area of 0.6 × 0.6 m 2 and a length of 5 m. A false bottom, 190 mm above the channel bottom, is used to generate a defined boundary layer, and to allow the water tunnel to be operated at a reduced water depth.

Flow facility
Free surface waves are generated behind an obstacle, which is mounted on the false bottom at a distance of 0.85 m from the leading edge. The shape of the obstacle is defined by a fourth-order polynomial with H b = 0.117 m the obstacle height, and L = 0.295 m the obstacle half-length (Gui et al. 2014). (1) In the current work, a number of free surface wave patterns is obtained by varying the Froude number, while keeping the initial water depth constant. The free surface wave pattern is defined by the water depth above the obstacle ( H 0 − H b ) and the upstream Froude number Fr = U 0 ∕ √ gH 0 , where U 0 is the upstream bulk velocity, g is the gravitational acceleration, and H 0 is the initial water depth. The upstream Froude number is always sub-critical ( Fr < 1 ). However, transition to critical conditions ( Fr = 1 ) occurs at or near the obstacle (Gui et al. 2014).
The upstream flow conditions, at constant water depth ( H 0 ), are determined prior to the free surface wave height measurements ( Table 2). The liquid velocity is measured with a disk-type programmable electromagnetic liquid velocity meter (P-EMS E30, Deltares) where one of the axes is aligned with the flow. The probe is calibrated for a velocity range of ± 1 m s −1 with an accuracy of 1 % . The initial water depth ( H 0 ) is measured with a ruler.

Reference measurement
A commonly applied PLIF system is used as a reference for the stereo-PLIF free surface measurements (Buckley and Veron 2016;Duncan et al. 1999). The light sheet from a Nd:YLF laser (LDY 304PIV laser, Litron) illuminates the liquid containing a fluorescent dye (Rhodamine WT at 120 mg m −3 ). The concentration of fluorescent dye is low enough such that it does not influence the static surface tension of the fluid (Appendix 1). The dynamic surface tension is more appropriate for steep waves where compression can locally alter the surfactant concentration (Duncan et al. 1999). However, for the current application, the static surface tension is sufficient, but the effect of compression on the surfactant concentration at a time scale similar to the wave action needs to be investigated in future work. Images are acquired with a high-speed CMOS camera (Imager HS 4M, LaVision) equipped with a 105 mm Micro-Nikkor objective and a high-pass filter ( OG570 , Schott). The magnification ( M 0 ) at the center plane is about M 0 = 0.1 , with an object distance of Z = 1 m . A large depth-of-field is obtained, which with an aperture of f / 11, and a wavelength of = 527 nm , results in z ≅ 4(1 + 1∕M 0 ) 2 f # 2 ≈ 30 mm (Adrian and Westerweel 2011). The camera (Fig. 2) is placed at an angle ( ) of 15 • with respect to the light sheet (xy-plane) to avoid interference from the liquid meniscus (Belden and Techet 2011).
An inverse, third-order polynomial is used to determine the mapping from image to world coordinates (Soloff et al. 1997). The camera is calibrated using a two-plane dot-pattern target (Type 22, LaVision). A resolution of 10.1 pixels per millimeter is obtained over a field-of-view of approximately 180 × 180 mm 2 . The wave pattern is determined by the incoming flow ( U 0 ), the initial water depth ( H 0 ), and the obstacle height ( H b ). The cameras for the PLIF and stereo-PLIF are aligned on a yz-plane. b Front view of the set-up, orientation of the PLIF (reference) and stereo-PLIF system (camera 1 and 2), with respect to the light sheet location. Two-dimensional free surface measurements are obtained with the stereo-PLIF system. Reference measurements are obtained at the central light sheet location. c The light sheet is scanned with an oscillating mirror driven by a galvanometer over an angle interval of ( ) spanning a separation angle ( ) Table 2 Flow conditions for the free surface waves generated behind the obstacle Fr × 10 2 1 5.23 ± 0.17 0.144 0.189 4.40 ± 0.14 2 6.40 ± 0.17 0.144 0.189 5.39 ± 0.14 3 11.2 ± 0.19 0.144 0.189 9.43 ± 0.16

Stereo-PLIF measurement
The stereo-PLIF system extends a PLIF system to a threedimensional domain with a scanning light sheet and a stereocamera set-up. The method is comparable to conventional techniques such as scanning-PIV (Brücker 1996). Images are acquired with two high-speed CMOS cameras (Imager HS 4M, LaVision) equipped with a 105 mm Micro-Nikkor objective and a high-pass filter ( OG570, Schott). The magnification at the center plane is approximately M 0 = 0.1 , which with an aperture of f/11 results in a depth-of-field of z ≈ 30 mm. The cameras are placed in a stereo configuration (Fig. 2), with a full separation angle ( 2 ) of either 30 • or 50 • . These two full separation angles are imposed by optical limitations in the experimental facilities, respectively, the Multiphase Wave Lab 1 at MARIN in The Netherlands and the water tunnel of the Laboratory for Aero-and Hydrodynamics at the Delft University of Technology. Optimal accuracy, for conventional stereo-PIV applications, is obtained for a full separation angle of 60 • ≤ 2 ≤ 90 • (Lawson and Wu 1997). The two angles are used to determine the influence of the separation angle on the measurement accuracy at comparable focal points ( L f ≈ 1 m ). Consequently, the nominal magnification remains the same for both camera separation angles. Therefore, the error is expected to scale as e r ∝ z ∕ x ∝ tan( ) −1 (Lawson and Wu 1997).
Two-dimensional free surface height measurements are obtained by scanning the light sheet in a sawtooth profile over the three-dimensional measurement domain with an oscillating mirror driven by a galvanometer ( CT-6210H , Cambridge Technology) located at y c = 1745 mm above the false bottom (Fig. 2a). The large distance of the scanning mirror ensures almost vertical measurement planes (Fig. 5a), with a maximum angle deviation, with respect to the y-axis, of approximately 2.3 • at the outer edges of the measurement domain ( z d ≈ 62 mm ) (Fig. 5b).

Edge detection procedure
An accurate edge detection method is required for the calibration procedure and the free surface reconstruction. Large variations in image intensity (I) are observed as a result of the liquid properties (Fig. 3a). The variations are a result of refraction, light focusing, and reflection at the air-water interface. Traditional edge detection can result in detection errors due to the semi-reflective properties of the air-water interface (André and Bardet 2014).
The edge detection error depends also on the camera separation angle. The specular bias is reduced at small camera separation angles, but the quantization errors increase. Hence, there is trade-off between detection errors and quantization errors (Benetazzo 2006). Furthermore, the variation in intensity along the laser light sheet can introduce other angle dependent errors (André and Bardet 2014). These errors can be reduced with other detection methods, such as hyperbolic tangent fit methods (Hwung et al. 2009). However, the detection accuracy is not necessarily improved, and often these methods require an increased computational effort. Instead, a multi-step intensity-based detection procedure is used in the current work.
The edge information is obtained from the raw image ( 2016 × 2016 pixels ) with a multi-step intensity-based detection procedure (Fig. 3a). First, the background intensity variation, based on the windowed mean and standard deviation, Image processing steps for edge detection of the free surface waves. a The raw image ( 2016 × 2016 pixels ) is cropped to show the surface waves. b The threshold distribution is determined over the top part of the image ( 500 pixels ), which does not contain free surface information. c The image is thresholded and smoothed with a ( 3 × 3 ) Gaussian filter. The edge is detected based on the maximum value per column is determined over the first part of the image ( 500 pixels ). Then, a discontinuous threshold distribution is determined per window, which is defined as the mean plus two times the standard deviation (Fig. 3b). Finally, a continuous threshold distribution is obtained by fitting a second-order polynomial to the discontinuous data.
The raw image is binarized with the continuous threshold distribution. Then, the image is smoothed with a twodimensional 3 × 3 Gaussian filter to suppress the small-scale image noise (Fig. 3c). Finally, the edge is defined as the pixel with a gray level above the threshold with the lowest y-value (i.e. topmost "white" pixel in Fig. 3c).
Detection errors can occur with a per-column threshold operation. However, pixel data cannot be easily filtered, due to the variation in magnification over the domain. This would result in a variable filter length in world coordinates. Therefore, pixel coordinate data are processed without further filtering.

Stereo-PLIF calibration
The scanning stereo-PLIF system is calibrated with a selfcalibration procedure after an initial calibration at the outer edges of the measurement domain (i.e., front and aft plane). An updated mapping function is determined at each light sheet location for the two-dimensional free surface measurements.
The self-calibration procedure is similar to typical stereo-PIV applications (e.g., Hori and Sakakibara 2004;Wieneke 2008). There is, however, a distinct difference. A stereo-PLIF measurement contains only curvilinear lines (i.e., the intersection of the light sheet with the surface) at a specific free surface height, whereas a stereo-PIV measurement contains information over the entire field-of-view. Images of the liquid free surface are obtained at several unique still water heights, hereby sampling the entire field-of-view. The whole process of the stereo-PLIF calibration procedure is detailed in the following sections.
First, a broad calibration is performed at two planes enclosing the measurement domain (i.e., positive and negative z), shown in Fig. 4, with a two-plane dot pattern target (Type 22, LaVision). An inverse third-order polynomial mapping function = F −1 0 i,j ( ) is determined for each camera (j) and plane (i), which maps pixel coordinates = (X, Y) to world coordinates = (x, y, z) (Soloff et al. 1997).
Then, images of the free surface are acquired at several still water heights ( N = 14 ) for each light sheet location (k). The edge detection procedure is applied to obtain curvilinear lines ( L k,n ) for each light sheet location (k) and still water height (n). These lines span the entire field-of-view.
Compared to stereo-PIV, the lines do not contain identifiable points (i.e., matching world and pixel coordinates).
Furthermore, identifiable points can not be obtained with the approach of Hori and Sakakibara (2004), as the scanning mirror is not calibrated. Identifiable points are obtained by distance minimization between camera projection lines ( R j ( ) ), which are the rays formed by back-projection of single pixels to the enclosing calibration planes (Fig. 4). The projection lines of both cameras are matched based on the minimal distance between skew lines with a maximum distance threshold (Gellert et al. 1989). The point of minimal distance defines the world coordinate ( ) with corresponding pixel coordinates for both cameras (Fig. 4).
Next, an updated inverse polynomial mapping is determined for each camera and light sheet location (k). The matched features are all constrained on a light sheet (i.e., a single plane), which enables the use of back-projection methods. The inverse mapping defines world coordinates (x, y, 0) = F −1 j,k ( ) on a light sheet plane (Adrian and Westerweel 2011). The mapping function is completed by fitting a plane to the z-coordinates ( z = f (x, y)).
Finally, the calibration domain is aligned with the liquid free surface, as the rough calibration is not necessarily aligned with the still water level. The two-dimensional free surface height y = y(x, z) is determined with the calibrated stereo-PLIF system. A plane is fitted to the reconstructed L k,n ∼ still water level P la n e 2 P la n e 1 L ig h t sh ee t ), at the outer edges of the domain (i.e., Plane 1 and 2). The projection lines ( R j ) formed by backprojection are used to triangulate the still water level (i.e., curvilinear line L k,n ). An updated mapping F −1 j,k is determined at each light sheet location (k) with the triangulated world ( ) and pixel ( ) coordinates. The inset shows a side view of the intersecting projection lines and the calibration planes ( P 1 , P 2 ) at the outer edges of the domain two-dimensional free surface, which is subtracted as a correction for the misalignment.

Post-processing
Post-processing is applied to the world coordinates obtained with both stereo cameras. The post-processing is applied to obtain uniformly distributed coordinates. Furthermore, the obtained data are filtered and averaged to obtain a combined stereo-PLIF measurement.
First, the pixel data obtained with the edge processing procedure are mapped to the respective measurement plane (i.e., light sheet location k) using the updated mapping function. The stereo camera data cannot be directly combined as the world points are non-uniformly distributed. Therefore, the world data are interpolated, with a linear interpolation method, to a uniformly distributed grid over the x-direction with a spacing of x = 0.1 mm . Higher-order interpolation methods are in the current application not required, as a grid point is displaced by only || − i || ≈ 0.04 mm on average.
Then, the now uniformly distributed stereo data is filtered. First, a Hampel filter, with a filter length of L h,f = 2.7 mm , is applied, which removes values that deviate more than three standard deviations from the median over the filter length (Liu et al. 2004). Then, to smooth the data, a second-order Savitzky-Golay finite impulse response filter is applied, with a filter length of L h,s = 2.7 mm (Orfanidis 1995). Finally, the stereo-PLIF data are obtained by averaging the data from both stereo cameras.

Calibration accuracy
The resolution of the initial polynomial mapping is not constant over the field-of-view (Table 3). The large angle of the stereo cameras, with respect to the z-axis, results in a magnification change over the field-of-view. Hence, there is a significant difference between the vertical and horizontal resolutions, where the horizontal resolution changes approximately 20% over the field-of-view. Furthermore, the accuracy changes as function of the camera separation angle (Lawson and Wu 1997).
A measure of accuracy for the self-calibration procedure is the standard deviation per still water height. Therefore, a selection of still water levels used in the self-calibration procedure for a calibration domain of z d ≈ 62 mm , measurement plane spacing of z ≈ 3 mm , an initial domain size of z d = 80 mm , and camera separation angle of 2 = 50 • is shown in Fig. 5a. The standard deviation varies from up = 0.52 mm to down = 0.30 mm over the height −110 ≤ y ≤ 50 mm with an average of = 0.38 mm ; see Table 3. The higher standard deviation of up is either caused by the increased magnification at the top of the domain or Table 3 Resolution change over the field-of-view for the initial polynomial mapping. The initial mapping is defined at the outer edges of the domain (i.e., plane 1 and 2). The resolution change of the stereo-PLIF system is presented per separation angle ( 2 ) and averaged over the two initial calibration planes the limited depth-of-field (i.e., the edge of the focal domain), which is evident from the larger spread over the x-direction (Fig. 5b). The measurement planes are positioned at a small angle, which remains minimal ( 2.3 • ) even at the outer edges of the domain (Fig. 5b). Furthermore, the angle could be removed easily by interpolating the data to a vertical plane. The virtual location of the scanning mirror, that is the intersection point of the measurement planes, is located at y ≈ 1743 mm above the false bottom, which corresponds with the measured value of y c ≈ 1745 mm.
The accuracy of the self-calibration procedure is further evaluated by comparing the still water levels reconstructed with the reference and stereo measurement system (Fig. 6). Measurements are performed with the stereo cameras at two different separation angles ( 2 ), while maintaining a constant focal point. The zero level correction is performed at equal still water level for both the stereo and reference measurement. However, the zero level is not equal for the two separation angles.
The free surface height is reconstructed accurately over the domain size of − 100 ≤ y ≤ 100 mm without significant influence of the domain size ( z d ) and inter-plane spacing ( z ) for both camera separation angles (Fig. 6). The difference between y stereo − y ref over the domain size, for either separation angle 2 = 30 • (indicated by open squares) or 2 = 50 • (indicated by open circles), is relatively small, as shown in the inset of Fig. 6. For both separation angles, the error remains within 1% of the target domain size (L) of 100 mm . A systematic error is observed, which clearly displays a parabolic behavior.

Hydraulic flow
The free surface waves behind an obstacle are measured at two separation angles ( 2 ) for three different Froude numbers ( Table 2). The free surface wave patterns vary over a significant height, which is ideal for a proof of principle measurement. Furthermore, the wave patterns display a variety of length scales from capillary waves to large free surface undulations as shown in Fig. 7. These free surface undulations are clearly observed to travel downstream as shown in the animation (Online Resource 2).
An instantaneous snapshot of the free surface height (y) is shown in Fig. 8a, d, g for the three different Froude numbers. The measurements are obtained for a domain size of z d ≈ 62 mm , plane spacing of z = 3 mm , and a separation angle of 2 = 50 • . The free surface height is reconstructed over N p = 19 measurement planes with an acquisition frequency of f a = f ∕N p ≈ 53 Hz per scan, where the cameras and scanning mirror operate at f = 1 kHz.
The differences in amplitude and wavelength of the free surface waves are apparent. The temporal variance of ( (x, z, t) = y − y ), with respect to the temporal mean y , is clearly observed in a space-time diagram. These space-time diagrams are obtained at two locations, spanning the z-direction, as shown in Fig. 8b-f, h, i).
Limited free surface disturbances, without a repetitive pattern, are expected at x = − 40 mm (Fig. 8b, e, h), as the upstream flow is relatively uniform. There is indeed no apparent pattern visible over the z-direction. However, the temporal fluctuations ( y −ȳ ) for the first two cases are still significant ( 1 = 0.34 and 2 = 0.48 mm ). These periodic wave shape variations over time (Fig. 10b, e) are introduced by fluctuations in the upstream flow. For Fr = 0.094 (Fig. 10h) the velocity fluctuations have a negligible influence on the drowned hydraulic jump. Therefore, the standard deviation for the measurement at Fr = 0.094 is considered as a representable error measure, which has a magnitude of 3 ≈ 0.15 mm . The bias error, between the stereo and reference measurement, defined as the difference between the temporal average of the stereo and reference measurement, is 1.0 mm for Fr = 0.094 , or 1% of the target domain size.
The second location is defined for each Froude number, to capture the changes in the free surface behavior. For Fr = 0.044 at x = 20 mm (Fig. 8c), small variations in spanwise wave behavior are observed. For the Fr = 0.054 at x = 10 mm (Fig. 8f), the spanwise variations are more pronounced with an increased disturbance amplitude. Furthermore, large, spanwise moving, surface disturbances are observed. For the final measurement at Fr = 0.094 , a The shape of the crescent-shaped disturbance, with preceding capillary waves, observed in Fig. 8d resembles the instabilities observed on spilling breakers (Duncan et al. 1999;Su et al. 1982). The surface disturbance shown by Gui et al. (2014) is, however, not observed with the current free surface measurements.
It is found that bubbles on the surface seriously deteriorate the signal quality, so that the measurement technique is not suited for flows with a lot of air bubbles. An example of such a detection error is shown in Fig. 9. The viewing lines of stereo-camera 1 cross the interface of an air bubble (Fig. 9a). This is registered as an additional interface above the true interface, due to the reflection of light in the air bubble. Camera 2 registers a shadow below the true interface, as the fluorescent light is reflected (Fig. 9b). However, the error introduced by the air bubble is limited (Fig. 9c). In the current work, a quality indicator (e.g., the vector norm of the distance between triangulated and mapped data) is not included, but more elaborate validation and/or processing can easily be implemented for particular situations (e.g., do not average if one of the cameras has "step-like" results).
A comparison between the temporal mean of the stereo and reference measurements is shown in Fig. 10a. The average difference is negligible ( 0.2 mm or 0.2% of the target domain size) with a maximum deviation of 0.45 mm . The error for Fr = 0.044 seems to be considerably larger. However, the standard deviation of the difference is, for both cases (i.e., 1 = 0.45 and 2 = 0.28 mm ), comparable to the previous error definition. The amplitude of the free surface disturbance is estimated as A 1 ≈ 3.2 mm and A 2 ≈ 8.6 mm , which results in a relative error ( i ∕A i ), as defined in Table 1, of 14% and 10% respectively. However, the relative error compared to the target domain size ( i ∕L ) is negligible ( ∼ 0.4% ). For the smaller camera separation angle ( 2 = 30 • ), the averaged difference remains the same; whereas the maximum increases to 0.81 mm . The reference camera system also has an inherent error, which for a typical PLIF system is approximately 0.5 mm (Duncan et al. 1999). In this case, the typical error of the reference system is estimated as the error induced by a single pixel shift, which corresponds to 0.1 mm or approximately 1-3% of the amplitude. The errors of the stereo-camera system are comparable to typical PLIF systems (e.g., Buckley and Veron 2017;Duncan et al. 1999).
The space-time diagrams of the disturbances ( (x, 0, t) = y − y ) are, for Fr = 0.044 and Fr = 0.054 , shown in Fig. 10b, c. These results are obtained at the center line of the domain. The growth and decay of the free surface disturbances are visible. First of all, the amplitude of the free surface disturbances is observed to be larger for Fr = 0.044 , whereas the occurrence frequency is higher for case one. Second, for Fr = 0.054 a crescent type of disturbance is observed, where the small amplitude waves are observed to travel upstream (towards negative x).
The local phase speed of the disturbances is, for both cases, estimated with a line as c p,1 ≈ 24.7 and c p,2 ≈ 29.7 cm s −1 , respectively, (Fig. 10b, c). The estimate is checked by comparing the wavelength of the disturbances on the wave crest for Fr = 0.054 (Fig. 7a) Fig. 7 a Instantaneous free surface height for Fr = 0.054 . The obstacle center is located at x ≈ −80 mm . The free surface displays typical crescent shaped disturbances, which can be clearly observed in the animation (Online Resource 2). The thick black line represents the data shown in b. b The back-projected reference image is shown, with the stereo and reference data at the corresponding location. The inset shows the difference between the stereo and reference measurement, with an indication of the local peaks (open squares) and troughs (open circles). b Will be discussed later the difference between crests and troughs, is approximately 1 = 4.8 mm ( k = 1.3 × 10 3 m −1 ). The dispersion relation for gravity-capillary waves on still water is, with the surface tension, k = 2 ∕ the wavenumber, and H the water depth. For deep-water surface waves, the depth does not influence the dispersion relation (i.e., tanh(kH) ≅ 1 ) (Lamb 1993;Whitham 1999). The local phase velocity, (2) 2 = gk + k 3 tanh(kH), on still water, is defined as c = ∕k . Based on the wavelength of the first disturbance ( k = 1.3 × 10 3 m −1 ) with surface tension and density values of = 0.072 N m −1 , = 998 kg m −3 , respectively, the phase velocity is estimated as c p, 1 ≈ 32 cm s −1 . The wavelength of the following crest is estimated as 2 = 7.3 mm , which results in a phase velocity of c p, 2 ≈ 28 cm s −1 . The determined phase velocities are comparable to the estimate of c p,2 ≈ 29.7 cm s −1 (Fig. 10c). A small difference is expected as the space-time diagram is determined in a fixed frame of reference, whereas  Table 2). The ability of the method to measure capillary waves is indicated by the estimated wavelength and corresponding phase velocity.

Wave swell by water drop impact
As a second test case, a droplet impact is generated by releasing a water drop of approximately 5 mm in diameter from a height of 757 mm to impact on a deep-liquid pool ( H 0 = 141 mm ). The resulting free surface disturbances are recorded with the stereo PLIF system in exactly the same configuration as in the first test case, with a domain size of z d ≈ 62 mm , a plane spacing of z = 3 mm , and a camera separation angle of 2 = 50 • . The measurements are timed to obtain two impact events in one recording. However, the droplet release time is not synchronized with the recording time.
The free surface disturbances resulting from the droplet impact are clearly observed in the animation (Online Resource 4). Three instantaneous snapshots of the ring waves, generated by the second impact, are shown in Fig. 11a-c. The two markers indicate the locations where time series are extracted, as shown in Fig. 11d.
The first time instance (t1) is approximately 228 ms after the droplet impact at p ≈ (− 60.0, 0.0) mm (Fig. 11a). The maximum free surface disturbance (x, z, t) = y −ȳ is limited to ± 1.5 mm to remove the unphysical amplitude introduced by interference of either the incoming water drop or the rebounding liquid jet with the light sheet.
At the second time instance (t2), the marker ( • ) is located in the trough of the first wave train (Fig. 11b). The small amplitude and short wavelength waves preceding the second wave train are also observed. The marker ( △ ) observes the small wavelength waves preceding the first wave train. Furthermore, the second wave train generated by the droplet rebound is observed in the instantaneous snapshot of the free surface height.
Finally, the third time instance (t3) shows clearly the amplitude decay (Fig. 11c), as energy is conserved, while the front and back group velocities diverge. The amplitude decay is proportional to time as ∝ t −1∕2 (Whitham 1999).
A space-time diagram displays the free surface variation along the center line ( z = 0 mm ) of the measurement domain, which is used to estimate the local phase velocity ( c p ) (Fig. 12). The minimum group velocity ( c g,min ) can, however, not be estimated according to Moisy et al. (2009), as the wave amplitude is below the measurement accuracy of the current stereo-PLIF configuration.
The local phase velocities of the wave crests, resulting from the droplet impact and rebound, are estimated as c p ≈ 27.5 cm s −1 using the space-time diagram in Fig. 12. The local phase velocity is larger than the minimum phase velocity for surface waves c p,min = (4g ∕ ) 1∕4 = 23.1 cm s −1 with a wavelength at minimum phase speed of p,min = 2 √ ∕( g) = 17.3 mm . The wavelength of the surface waves, determined as the distance between the first peak and trough at t 3 (Fig. 11c), is approximately ≈ 44 mm . The observed waves are deep-water gravity waves ( < H 0 ∕2 ) for which the local phase velocity c p = √ (g + −1 k 2 )∕k ≈ 28.1 cm s −1 is comparable to the local phase velocity estimated from the space-time diagram shown in Fig. 12.
After, approximately, t r ≈ 2.6 s another free surface disturbance occurs (Figs. 11d,12). This free surface disturbance is induced by the reflection of the initial wave train  Fig. 9 Image of an edge detection error due to the presence of an air bubble on the liquid surface at the recirculation area of the drowned hydraulic jump ( Fr = 0.094 ). a Stereo camera 1 looks through the interface of the air bubble, which results in a detection error close to the liquid free surface. b Stereo camera 2 does not look through the bubble, but the bubble appears as a shadow below the interface. c The local free surface height per camera (cam 1 / 2) and filtered free surface height (comb.) are shown at the detection error location. The inset shows a zoom at the detection error location, which shows an error of approximately 0.4 mm from the side walls of the channel. The reflected wave traveled a distance of z ≈ 600 mm during t r ≈ 2.6 s , which corresponds to a local phase velocity of c p ≈ 23 cm s −1 . This local phase velocity is close to the minimum phase velocity ( c p,min ) of surface gravity waves.
A skew is introduced in the two-dimensional free surface measurements ( (x, z, t) ), as is clearly observed in Fig. 13a. The scanning system acquires line measurements over the z-coordinate with a temporal spacing of t = 1 ms . In the current example, there are N p = 19 light sheet locations. The difference between the first and last measurement plane is significant ( 18 ms ). Based on the local phase velocity, this introduces a relative displacement of x = c p ( tN p ) ≈ 5.2 mm . The measurements are, therefore, non-instantaneous, as the measurement time is too large which introduces a significant displacement.
Temporal interpolation can be applied to reduce the skew (Wellander et al. 2014) (Fig. 13d). The instantaneous free surface height (Fig. 13a) is interpolated to the center plane time (t), where a linear change in amplitude is assumed (Fig. 13b). The absolute difference between original ( ) and interpolated ( i ) data shows an amplitude deviation of | i − | ≈ 0.3 mm with a horizontal shift of x ≈ 2.5 mm (Fig. 13c). The horizontal shift corresponds to the displacement based on the

Conclusion
This article presents a scanning stereo-PLIF measurement system. The measurement system uses a stereo-camera set-up with a self-calibration procedure adapted for free surface flows. Two-dimensional free surface height measurements can be acquired with a temporal resolution of 19 ms ( f = 53 Hz ), limited only by the available scanning equipment ( 1 kHz rate).
To the best of our knowledge, this system is unique as it allows us to measure time-dependent two-dimensional free surface height (y(x, z, t)) with disturbance amplitudes of ≈ 0.2 mm , over a domain height L of 100 mm . The accuracy with respect to the domain height (L) is high (0.2%). On the other hand, accuracy with respect to the wave height is 12%, which is comparable with the accuracy of other projection-based approaches ( Table 1). The temporal resolution is lower than other projection-based approaches, but the temporal resolution is reduced in favor of the spanwise resolution with currently 19 measurement planes at a spanwise spatial resolution of z ≈ 3 mm . The streamwise spatial resolution is high ( x = 0.1 mm ) compared to correlation-based approaches.
The temporal resolution is currently limited as is evident from the introduced skew, but results can be skewcorrected with temporal interpolation. Furthermore, the technique is versatile in its application. Therefore, the stereo-PLIF technique shows promise in the reconstruction of small amplitude variations over large two-dimensional surfaces with minimal influence on fluid properties.  Fig. 11a-c. The local phase velocity ( c p = 27.5 cm s −1 ) of the two wave crests, generated by the droplet impact and rebound, are indicated with two dashed lines, and the interference resulting from the reflected wave is shown by the dash-dotted line at t r Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.

Appendix 1
The change in fluid properties due to the addition of the particular fluorescent dye used for the stereo-PLIF is investigated. The static surface tension is measured with an adapted pendant drop method (Saad and Neumann 2016). The method minimizes the difference between the parameterized Young-Laplace equation and an experimentally obtained profile to estimate the static surface tension value. The experimental profile is determined from the shadowgraph of a pendant drop.
The pendant drop is suspended from a blunt tip needle (Fig. 14a). The droplet volume is controlled with a syringe pump (101 syringe infusion pump, KD Scientific), allowing volume adjustments of 600 L s −1 . The droplet is imaged with a CCD camera (VC-Imager Pro X 4M, LaVision) equipped with a long-distance microscope (QM 1 Long Distance Microscope, Questar). The back light illumination (OSL2 High-Intensity Fiber-Coupled Illuminator, Thorlabs) is uniform over the entire field-of-view of 5 × 5 mm . The Fig. 13 Temporal interpolation to reduce the skew introduced by scanning the three-dimensional domain. a Instantaneous free surface height ( (x, z, t) = y −ȳ ), with ȳ the still-water level, for a droplet impact. b Skew-corrected instantaneous free surface height ( i (x, z, t) = y i −ȳ i ). c Absolute difference between original ( ) and interpolated data ( i ). d Method used to interpolate data to common time (t) edge is easily obtained with a gradient-based edge detection procedure (Fig. 14c). The static surface tension is determined by minimizing the difference between the theoretical and experimental profile (Fig. 14b). The static surface tension of water-air at T = 21.3 • C is measured as = 72.58 mN m −1 . The influence of Rhodamine 6G with increasing molar concentration ( c i ) has previously been investigated, where it was shown that there is a negligible effect for practical concentrations (André and Bardet 2015). The influence of Rhodamine WT is investigated over a similar concentration span (Fig. 15). The static surface tension value is slightly overestimated, but the trend is similar to the previous work. The effect on the static surface tension is limited for the current application of the stereo-PLIF system, as indicated by the filled marker ( ▪ ) (Fig. 15). However, the dynamic surface tension needs to investigated for cases where compression can locally alter the surfactant concentration.
Fig. 14 Image processing to obtain the static surface tension from pendant drop images. a The experimental equipment used to obtain the shadowgraphy images. b The static surface tension of a water droplet is determined by minimizing the difference between the theoretical and experimental profile. c Image used to obtain the experimental profile for the minimization