On the calibration of Astigmatism particle tracking velocimetryfor suspensions of different volume fractions

In the present study we demonstrate for the first time how Astigmatism Particle Tracking Velocimetry (APTV) can be utilized to measure suspensions dynamics. Measurements were successfully performed in monodisperse, refractive index matched suspensions of up to a volume fraction of $\Phi=19.9\%$. For this, a small percentage of the particles is labeled with fluorescent dye acting as tracers for the particle tracking procedure. Calibration results show, that a slight deviation of the refractive index of liquid and particles leads to a strong shape change of the calibration curve with respect to the unladen case. This effect becomes more severe along the channel height. To compensate the shape change of the calibration curves the interpolation technique developed by Brockmann et al. (Experiments in Fluids, 61(2), 67, \citeyear{brockmann2020utilizing}) is adapted. Using this technique, the interpolation procedure is applied to suspensions with 6 different volume fractions ranging from $\Phi<0.01\%$ to $\Phi=19.9\%$. To determine the effect of volume fraction on the perfomance of the method, the depth reconstruction error $\sigma_z$ and the measurement volume depth $\Delta z$, obtained in different calibration measurements, are estimated. Here, a relative position reconstruction accuracy of $\sigma_z$/$\Delta z$=0.90\% and $\sigma_z$/$\Delta z$=2.53\% is achieved for labeled calibration particles in dilute ($\Phi<0.01\%$) and semi-dilute ($\Phi\approx19.9\%$) suspensions, respectively. The measurement technique is validated for a laminar flow in a straight rectangular channel with a cross-sectional area of 2.55$\times$30\,mm$^2$. Uncertainties of 1.39\% and 3.34\% for the in-plane and 9.04\% and 22.57\% for the out-of-plane velocity with respect to the maximum streamwise velocity are achieved, at solid volume fractions of $\Phi<0.01\%$ and $\Phi=19.9\%$, respectively.


Introduction
Particle laden flows are widely present in nature and technical settings such as sediment transport in rivers or industrial waste slurry transportation. In such flows, the particle-particle interaction and the hydrodynamic interaction between particles and fluid can give rise to various phenomena such as shear induced migration, inertial migration, shear thinning and shear thickening (Leighton and Acrivos 1987, Ho and Leal 1974, Stickel and Powell 2005. Both, numerical simulations and experiments are performed to reveal the underlying physical mechanisms behind these phenomena (Morita et al. 2017, Shichi et al. 2017. Recent progress in numerical algorithms and computer technology enables researchers to examine the behavior of individual particles in great detail in suspensions even at high solid volume fractions (Kazerooni et al. 2017). Nonetheless, only a few experimental studies address the individual particle dynamics in suspensions beyond the dilute regime by means of optical measurement techniques such as Particle Image Velocimetry (PIV) (Zade et al. 2018;, Zhang et al. 2018 or a combination of PIV and 2D Particle Tracking Velocimetry (PTV) (Baker and Coletti 2019). Optical measurement techniques have been proven to be reliable tools for investigating fluid flow problems with high spatial and temporal resolutions. However, when it comes to suspension flows at high solid volume fractions, they can not be utilized as the flow becomes opaque, given the presence of a large number of particles.
A solution for this issue is refractive index matching (RIM) such that the resulting suspension becomes transparent (Wiederseiner et al. 2011). As the required monodisperse particles as well as the RIM liquids are expensive and in most cases hazardous, it is desirable to perform the experiments in micro-scale setups. Micro-scale experiments in turn usually only provide restricted optical access, such that most three dimensional optical measurement techniques, which require complex calibration procedures, can not be employed. In this regard, Astigmatism Particle Tracking Velocimetry (APTV) appears to be a cost-effective and easy-to-implement method for measuring the threedimensional motion of both liquid and particles.
APTV is a single camera technique that can be applied to micro-scale enviroments with limited optical access (Cierpka et al. 2012). The basic idea of APTV is to introduce a controlled astigmatism to an optical system such that the particle's image deforms based on its depth position. Hence, the particle out-of-plane position can be reconstructed by the shape of the particle's image. While the origin of APTV is found in microfluidic applications (Kao et al. 1994), the technique was improved by various authors (Angarita Jaimes et al. 2006, Chen et al. 2009, Rossi and Kähler 2014 and has been employed to investigate numerous physical phenomena at length scales ranging from 9 µm (Huang et al. 2016) up to 5.1 mm (Buchmann et al. 2014). A detailed overview on APTV, its applications and its background is given by Cierpka and Kähler 2012. Apart from utilizing APTV in combination with small tracer particles to measure the motion of liquids, in recent years the technique is more and more applied to measure the dynamics of the particles themselves in situations where particles undergo their own dynamics and do not follow the flow streamlines. Rossi et al. 2019 used APTV to measure particle velocity and concentration of an electrokinetically induced particle pattern in a suspension with a solid volume fraction of Φ = 0.05% inside a channel of 350×30 µm 2 crossection. They used a mixture of labeled particles with different dyes (d p = 0.245 µm) with the majority of them being invisble during the APTV measurement to allow for higher solid volume fractions. They concluded that APTV is capable of measuring the concentration profile with a better resolution compared to classical segmentation approaches. Using APTV Blahout et al. 2020 investigated the fractionation of d p = 3.55 µm and d p = 9.87 µm particles in a serpentine channel with a crossection of 200×50 µm 2 . They observed, that a transition of particles from four to two equilibrium trajectories takes place at particle size dependant bulk Reynolds numbers.
Most recent works increasingly focus on using machine learning tools such as deep neural networks to further increase the degree of automatization of APTV (Rossi and Barnkob 2019) or to apply it on scenarios with low signal-to-noise ratios (Franchini et al. 2019). König et al. 2020 compared the performance of conventional and neural network supported APTV utilizing a bidisperse suspension in a laminar channel flow. While they achieved good results with both methods for the bidisperse suspension, they concluded that the neural network supported APTV is more robust against optical abberations and can be of great use for investigating suspensions with different particles sizes, shapes and even with particle clusters.
Efforts are also undertaken to further adapt and extend the classical APTV technique. Zhou et al. 2020 introduced a modification utilizing Holographic imaging principles to establish Holographic Astigmatic Particle Tracking Velocimetry (HAPTV). Using a nozzle flow as a test case, they could successfully validate their measurement technique. Brockmann et al. 2020 developed an adapted APTV procedure, referred to as Ball Lens Astigmatism Particle Tracking Velocimetry (BLAPTV), for transparent particles that are large in comparison to the field of view. They extended the 2D Euclidean calibration approach developed by Cierpka et al. 2010 to the 3D space by additionaly considering the particles light intensity. Furthermore, they developed an interpolation method to account for shape changes of the calibration curve that occur in a channel along the channel height. The authors validated their proposed method by measuring the velocity profile of a laminar channel flow.
In the present study we combine conventional APTV with refractive index matching to investigate the dynamics of large particles (d p = 60 µm) for the first time by means of APTV in suspensions of up to Φ = 19.9% volume fraction. Here, a small amount of particles is labeled with a fluorescent dye while the majority of the particles is invisible to the camera. The key to successful measurements is here that the interpolation method developed by Brockmann et al. 2020 is successfully applied to compensate the effect of remaining refractive index mismatches in the dense suspension.
2 Experimental set-up The measurement system consists of a microscope (Nikon Eclipse LV100) illuminated with a 15 W continuous green laser of 532 nm wavelength. The laser is operated at 1.5 W. For image-recording a 12-bit, 1280×800 pixel CMOS high-speed camera (Phantom Miro Lab 110, Vision Research) with 20 µm pixel size is used. A shematic drawing of the experimental setup is shown in Fig. 1a. Measurements are performed with a Nikon Cfi60 objective lens of M =10× magnification. To introduce astigmatism, a cylindrical lens with a focal length of f cyl =200 mm is placed in front of the camera sensor, generating two spatially separated focal planes with a measured distance of approximately 192 µm. The basic configuration has been already used by Brockmann et al. 2020 with a bright filed illumination. For the present experiments a plane channel with a cross sectional area of h×w=2.55×30 mm 2 and a length of 300 mm was realized. Velocity profiles were measured 150 mm downstream of the channel entrance. The flow is generated with a high pressure syringe pump (LA-800, Landgraf HLL GmbH) and a 100 mL syringe (Braun GmbH). For obtaining a density and refractive index matched suspension with d p = 60 µm PMMA particles we have used the receipe proposed by Bailey and Yoda 2003 with a ternary mixture of 24.85wt% water, 36.03wt% glycerin and 39.12 wt% ammonium thiocyanate that has a refractive index of n RIM = 1.4867, a density of ρ RIM = 1.19 g m −3 and a dynamic viscosity of η RIM = 4.99 cP.

Calibration procedure
In APTV, typically a cylindrical lens is implemented in the optical system which alters the light path such that two spatially separated focal planes F xz and F yz are generated. If the particle is located at a z position in the middle of the focal planes the particle's image appears circular and skews into an oblate or prolate ellipsoid when the particle's z position is closer to F xz or F xz , respectively (see Fig. 2 and Fig. 3). This shape change of the particle image can be quantified by the length of the horizontal and the vertical axis of the autocorrelated particle image denoted as a x and a y , respectively. To reconstruct the particle's out of plane position based on the values of a x and a y a calibration function is required (Cierpka et al. 2010). Here, a 2D Euclidean calibration procedure is presented similar to that of Brockmann et al. 2020. The main difference is that we utilize fluorescent particles in the present work, therefore the abberated particle image itself is used here for particle out-of-plane position reconstruction. Furthermore, in contrast to Brockmann et al. 2020 we do not use the 3D Euclidean calibration procedure within this work, as the distribution of the light intensity is not sufficiently homogenous to improve the accuracy through a 3D procedure, here. The major steps of the calibration procedure are summarized in the following, where we : Images of labeled particles for different depth positions z (d p =60 µm, PMMA, z corrected for refractive index of RIM-liquid (n RIM = 1.488)). The particle is located close to F yz at z ≈ −200 µm and located in F xz at z = 0 µm. a) Labeled particle in a suspension with Φ = 0.01%. b) Labeled particle in a suspension with Φ = 19.9%. also outline the differences in the calibration curves for a labeled particle in a dilute (Φ < 0.01%) and in a dense refractive index matched suspension (Φ = 19.9%). To capture the change of a x and a y for different out-of plane positions of particles, labeled and wall attached particles in a suspension of Φ < 0.01% and Φ = 19.9% are scanned in steps of 1 µm over a distance of 1000 µm such that the deformation of the particle image is entirely captured. Fig. 3 shows particle images of such a scan of a 60 µm PMMA particle labeled with Rhodamin B and located at the bottom channel wall. A magnification of M = 10× for Φ < 0.01% (Fig. 3a) and Φ = 19.9% (Fig. 3b) is chosen. The particle shown in Fig. 3a is submerged in a dilute suspension where only labeled tracer particles of d p = 60 µm are present in the RIM-liquid such that Φ < 0.01%. The change of a x and a y is determined as a function of z − z 0 . Results are displayed in Fig. 4a (large colored dots). The coefficient c a at which the isolines of the autocorrelation map are extracted to measure a x and a y will be denoted as autocorrelation coefficient as introduced and described in Brockmann et al. 2020. In Fig. 4b a y is plotted as a function of a x . This presentation poses the base for the Euclidean calibration method as developed by Cierpka et al. 2010 andadapted in Brockmann et al. 2020. The labeled particle displayed in Fig. 3b is submerged in a suspension of RIM-liquid and additional unlabeled particles such that the total volume fraction is Φ = 19.9%. In fact, the suspension was at rest for 24 hours, such that the particle shown in Fig. 3b is covered with appoximately 11 layers of unlabeled, transparent PMMA particles. From Fig. 3a and b it is obvious that the images of the particle in a suspension of Φ = 19.9% (Fig. 3b) are slightly blurred and also exhibit a speckle pattern for relative out-of plane positions ranging between 200 µm < z − z 0 < 400 µm. The blurriness and the speckle patterns are a result of slight deviations of the refractive index of individual transparent particles and the RIM liquid. The refractive index deviation of transparent particles which are located in the optical path between the labeled calibration particle and the objective, induces these distortions of the particle image affecting the values of a x (z − z 0 ) and a y (z − z 0 ). These deviations become evident in Fig. 4a (small colored dots) and Fig. 4b (small colored dots), respectively. Obviously the volume fraction affects the a x and a y values and particles displayed in Fig. 3a and b can not be treaded with the same calibration curve. is given in c. z − z 0 data is corrected for the refractive index of the RIM-liquid (n RIM = 1.488). a) Selecting z − z 0 range of scattered data by light intensity I (light blue dots=a x , dark blue dots=a y , green dots=I, black dots=I). b) Fitting polynominals of degree 14 to a x and a y (black line=polynomials a x , a y ). c) Reconstruction of z − z 0 of scattered a x -a y data (colored dots) by Euclidean distance (black line=polynomials, red dots=outliers).
To obtain a proper calibration function for both the dilute and the dense suspension the previously described calibration scans are repeated for several labeled particles randomly distributed over the field of view. This results in a data set of a x , a y and I as function of z − z 0 for several particles. An exemplary set of such scattered a x , a y and I data for d p = 60 µm particles submerged in a suspension at Φ < 0.01% is given in Fig. 5a. In a first step the median of I as a function of z −z 0 is computed, denoted as I. For the following steps only data points associated with I ≥ c I ·I max are considered, where c I is defined as the intensity coefficient. The value of c I ·I max will be hereafter referred to as intensity threshold I thr . For further information on c I the interested reader is referred to Brockmann et al. 2020. In the following step, a polynomial of 14th order is fitted to the scattered a x and a y data as shown in Fig. 5. When the polynomial fit of a y (denoted as a y ) is plotted as a function of the polynomial fit of a x (denoted as a x ) the 2D calibration curve is obtained. The z-position of a particle can now be reconstructed by assigning its measured a x , a y values to a point on the calibration curve that is given by the minimum Euclidean distance and then reading out the associated z − z 0 value. Pairs of a x , a y where the Euclidean distance exceeds a certain threshold are rejected as an outlier. This threshold is referred to as a D and defined as the mean Euclidean distance of the a x , a y data of all calibration particles multiplied by the factor c D as described in Brockmann et al. 2020. For the case depicted in Fig. 5 c D is set to a value of c D = 2, resulting in an out-of-plane reconstruction accuracy of σ z = 4.67 µm for d p = 60 µm particles. The measurement volume depth for the given example equals ∆z = 514.84 µm for c I = 0.5, such that the relative reconstruction accuracy, hence referred to as relative error, equals σ z /∆z = 0.90%. With c I = 0.4 the measurement volume depth is increased to ∆z = 606.28 µm whereas σ z increases slightly to 4.92 µm resulting in a relative error of σ z /∆z = 0.80%.

Validation measurements
To validate the applicability of APTV to dense suspensions, measurements are carried out in a plane channel flow with 2550 µm channel height. The suspension is a ternary liquid mixture proposed by Bailey and Yoda 2003 and PMMA particles of diameter d p =60 µm (Microbeads). The channel is filled with the RIM-liquid and a small amount of labeled particles such that Φ < 0.01%. Due to a slight density variation of the RIM-liquid and the particles, the particles settle to the channel bottom when rested over night, while a few are stuck to the top channel wall. Settled particles are used to determine the absolute position of the bottom and the top channel wall prior to the experiments. For this, the whole channel is scanned in steps of 1 µm, to record particles that are located at the top and the bottom wall within the field of view. The evolution of a x is used to detect the particle center and thereby the channel walls by considering the particle radius. In fact the particle center is focused in F xz when a x assumes a minimum. The origin of the scanning coordinate, is set to zero at the channel bottom. Hence, a constant suspension flow rate of 20 mL min −1 is induced by a high pressure syringe pumpe (LA-800, Landgraf HLL GmbH). A container is used to collect the liquid driven from the syringe through the channel.
Before the actual flow measurement is started, calibrations measurements are performed with particles located at the bottom channel wall to generate calibration curves for Φ < 0.01% as shown in Fig. 6 (solid line). Hence, flow measurements for suspensions with six different volume fractions ranging from Φ < 0.01% to Φ = 19.9% are performed. After the final measurements at Φ = 19.9% the setup is rested over night such that particles settle to the channel bottom. Then, calibration measurements are performed with settled and labeled particles to generate a calibration curve for Φ = 19.9% as shown in Fig. 6 (dashed line). As already discussed in section 3 both curves differ significantly due to particle image distortions. Test calibrations performed on particles located at the top and bottom of the channel revealed that this effect does not occur for a dilute suspension (Φ < 0.01%) in the present case. While for low volume fractions no effect can be noticed, for volume fractions larger than Φ ≥ 4.73%, we observed that the calibration curve changes with the z-coordinate and hence is an implicit function of the channel height. This is because the number of particles that disturb the light path is higher for labeled particles located closer to the bottom than for labeled particles located closer to the top of the channel. The challenge is to find a calibration function which is valid for labeled particles located at any z-position in between the bottom and top channel wall and for volume fractions in the range from Φ < 0.01% to Φ = 19.9%. To solve this problem we adapt the interpolation method developed by Brockmann et al. 2020. For this, we interpolate and extrapolate the polynomial coefficients of a x , a y based on the calibration curves for Φ < 0.01% and Φ = 19.9%. In this way, we compute 30 intermediate calibration curves. These are presented as colored lines in Fig. 6. The extrapolation allows us to generate calibration curves which are even more skewed than the calibration curve for Φ = 19.9% (see orange to red lines in Fig. 6). As we will show later these extrapolated curves are required at higher volume fractions. The z-range of all interpolated calibration curves is set to ∆z=514 µm. In addition to the calibration curve, which consists of a y as function of a x , the threshold for the maximum allowed Euclidean distance a D is also interpolated linearly. The out-of-plane position reconstruction uncertainties of particles of d p =60 µm diameter for Φ < 0.01% and Φ = 19.9% are σ z = 4.67 µm and σ z = 15.64 µm, respectively. The channel is scanned in steps of 223 µm during the flow measurements and at each measurement plane 2500 images with a resolution of 512×384 pixel, covering a 1.89×1.57 mm 2 field of view, are recorded at 100fps. In the post processing the image size is reduced to a region of interest of 300×300 pixel. By this, marginal areas of insufficient illumination are reduced and the computation time can be reduced significantly. After data acquisition, the particle positions and velocities need to be determined. To determine the out-of-plane positions of particles, based on their a y , a x values, each of the calibration curves displayed in Fig. 6, is compared with the a x , a y scatter data from the corresponding measurement planes. By this we can select an appropiate calibration curve for each measurement plane. To find the best fitting calibration curve, the number of valid a y , a x pairs that fulfill the Euclidean distance criterion is evaluated as described in Sect. 3. In fact, the curve that yields the largest number of valid particles is considered as a match and selected to determine the z-z 0 of the a y , a x pairs in the respective measurement plane.
Figures 7a-f display the best matching calibration curve (blue solid line) and corresponding valid data points (green dots) and outliers (red dots) as defined in section 3, for measurement planes located at z=468.09 µm, z=913.89 µm and z=2028.39 µm for Φ = 9.04% and Φ = 19.9% particle volume fraction, respectively. As can be seen in Fig. 7a it turns out for Φ = 19.9% at a measurement position of z = 468.09 an extrapolated calibration curve (blue line) fits best to the measurement data. The reason for this is, that during the flow measurement the transparent particles are well distributed along the channel height and not settled to a bottom layer as during the calibration. We therefore conclude, that the distortion induced by a static bot- Fig. 7: Best fitting calibration curve (blue line), valid a x -a y data (green dots) and rejected a x -a y data (red dots) for different measurement planes (c a = 0.5547, c I = 0.5, c D = 4). The solid and the dashed line are the calibration curves as presented in Fig. 6 for a static suspension at Φ=0.01% and Φ=19.9%, respectively. a) z = 468.09 µm, Φ = 19.9%, b) z = 913.89 µm, Φ=19.9%, c) z = 2028.39 µm, Φ=19.9%, d) z = 468.09 µm, Φ = 9.04%, e) z = 913.89 µm, Φ = 9.04%, f) z = 2028.39 µm, Φ = 9.04% tom layer of transparent particles on top of a labeled particle is less intense than the distortion created by homogenously distributed transparent particles at the same total volume fraction. Therefore an extrapolation of the calibration curve as displayed in Fig. 6 is needed to capture the calibration curve deformation during the flow measurement. An interpolation of the maximum Euclidean distance threshold a D as defined in section 3) is crucial as the scattering of a x -a y -I data varies along the gap height within the experiments (Brockmann et al. 2020). In a suspension flow the scattering increases the closer labeled particles are located to the channel bottom. This can be seen from Figures 7a-c and d-e where the maximum distance of valid data points (green dots) with respect to the calibration curve (blue line) decreases when the measurement plane is shifted towards the channel top (increasing z values). In the present study we use the light intensity of particle images for outlier detection. This is essential for the algorithm to reliably pick the best matching calibration curve. The importance of additionaly considering the light intensity as an outlier criterion can be better understood from Fig. 8, where we show the best fitting calibration curve and a typical distribution of I among the a x and a y scatter data. As can be seen, the light intensity increases close to the calibration curve. For using the light intensity as an additional outlier criterion we inter-and extrapolate the value of the intensity threshold I thr analogous to a D . Data points with I < I thr are rejected.
The whole procedure is applied to the a x , a y and I data of all measurement planes, so that for every measurement plane the out-of-plane positions of the particles are computed. Hence, the absolute particle positions can be computed with respect to the channel wall by considering the particle out-of-plane positions and the corresponding measurement plane position. In Fig.  9 we show which interpolated calibration curve (see Fig.  6) matches best with respect to the z position of each measurement plane. As can be seen for the lowest vol- Fig. 8: Typical distribution of scattered a x , a y and I data obtained in one measurement plane (z = 1136.79 µm) during the flow measurement for Φ=12.97%. All data gathered in the measurement is displayed. Blue line=best fitting calibration curve.
ume fraction (Φ < 0.01%) the algorithm picks calibration curves 5, 6 and 7 for calculating the out-of-plane positions. For z > 2000 µm and Φ < 0.01% the best fitting calibration curve matches with the calibration curve for Φ < 0.01% (number 5). For lower values of z, curves 6 and 7 are selected for Φ < 0.01%. However, the differences between the best matching calibration curves for Φ < 0.01% are small and it is sufficient to consider curve number 5 only for the whole measurement. On the contrary, the curve number changes significantly with increasing z for all flow measurements at higher particle volume fractions (Φ ≥ 4.71%). This change becomes more pronounced with increasing volume fraction and appears to be nonlinear as can be clearly seen from Fig. 9. Obviously for Φ ≥ 9.04% extrapolated calibration curves are required, as the maximum curve number exceeds 19. For Φ = 19.9% there is a plateau for z ≤ 468 µm, where the curve number is constant. It seems necessary here to extrapolate further to a higher curve number. However, the distortions of the scatter data for z ≤ 245 µm at the highest volume fraction of Φ = 19.9% are too strong to capture them by further extrapolation to higher curve numbers. Hence, for measurement planes z ≤ 245 µm there is a lack of valid data points and the error of the calculated particle position and velocity increases sharply, as will be also discussed hereafter.  Fig. 10a-f is defined as Re = (U Bulk D H ρ RIM )/ν eff with the effective viscosity estimated as ν eff = ν RIM (1 + 2.5Φ + 5.2Φ 2 ) according to Batchelor and Green 1972. The particle Reynolds number is estimated as Re p = Re(d p /D) 2 , where D denotes the gap height, as used for example by Shichi et al. 2017. In general, it can be seen that the measured velocity profile shows a good agreement with the analytical solution for all considered volume fractions. For higher volume fractions the maximum velocity is slightly lower than the velocity obtained from the analytical solution. This effect could be a result of particleparticle interactions. However the effect is small.
The uncertainty of determining the in-plane velocity increases from σ u =1.39% to σ u =3.34% with respect to the maximum streamwise velocity as the volume fraction increases from Φ ≤ 0.01% to Φ = 19.9%. The uncertainty for determining the out-of-plane velocity increases from σ w =9.6% to σ w =22.57% at the same time. We assume that the loss of accuracy with increasing volume fraction is related to two different effects. Firstly, the accuracy of determining the particles out-of-plane position decreases with increasing values of Φ as more layers of transparent particles inducing distortions to the images of the labeled particles. Secondly the accuracy of detecting the in-plane position of the particles centers decreases with increasing volume fractions due to the disturbances induced by the transparent parti- cles. This effect becomes visible in the uncertainty in determining the span-wise velocity component which increases from σ v =0.50% to σ v =1.78% as the volume fraction increases from Φ ≤ 0.01% to Φ = 19.9%.
The accuracies obtained within this work are lower than those achieved in APTV with small particles and dilute suspensions. This accounts especially for the out-of-plane velocity. For instance Cierpka et al. 2010 obtained uncertainties of σ u =0.9% and σ w =3.72% of u max . We assume that the uncertainty is larger due to slight refractive index mismatches between the RIMliquid and the transparent particles.
Regarding the number of detected particles it can be seen from Fig. 10, that there is no obvious lack of data visible among all the planes. Hence, we conclude that the described interpolation procedure is suitable for capturing the shape changes (a x , a y ) and the mean Euclidean scattering distance a D of the calibration curve. In Fig. 10c,d,e a slight decrease in the number of detected particles towards the middle of the channel becomes evident, which may be related to the measurement procedure and not to a physical phenomena. In fact, for all flow measurements the estimated particle Reynolds numbers are small (Re p << 0.1) as can be seen from Fig. 10a-f. Baghat et al. 2009 showed that inertial migration of particles occurs for Re p > 0.1. Hence we conclude that such effects can be neglected within our experiments.
Precisely measuring the particle concentration along the channel height is of great importance for the investigation of physical phenomena, however to developing an APTV based procedure to reliable measure the particle concentration is beyond the scope of this work.
As the particle Reynolds numbers are small and a good agreement could be observed between measured and theoretical velocity profile we conclude that the particle dynamics are negligible here. Thus, particles can approximately be considered as fluidtracers.

Discussion and Conclusion
In the present study we show that APTV can be applied to measure the particle dynamics of suspensions of up to 19.9% volume fraction. Measurements have been performed at six different volume fractions ranging from Φ ≤ 0.01% to Φ = 19.9%.
To make the suspension optically accessible we use a refractive index matched liquid (RIM-liquid) and transparent particles of which just a small portion is labeled with a fluorescent dye. Firstly, we study the effect of remaining image aberrations due to slight refractive index mismatches. We observe that slight deviations of the refractive index of individual transparent particles induce optical distortions that result in a shape change of the calibration curve of a labeled calibration particle. Thus, at high volume fractions of transparent particles, images of labeled particles get increasingly distorted the closer they are located to the channel bottom as more transparent particles disturb the optical path. To overcome this effect we adapt the interpolation method developed by Brockmann et al. 2020. Inter-and extrapolated calibration curves for labeled particles are generated from calibration measurements at Φ ≤ 0.01% and Φ = 19.9% and are related to the measurement data in a best fit procedure. In the present study, we used the light intensity as a simple outlier criterion. We found this additional outlier criterion to be crucial for a stable fitting procedure. Depth reconstruction accuracies of σ z /∆z = 0.90% and σ z /∆z = 2.53% were achieved for labeled static calibration particles of d p = 60 µm and a magnification of M = 10× in a suspension of Φ ≤ 0.01% and Φ = 19.9% volume fraction, respectively. Ultimately, we validated the interpolation technique successfully by measuring the laminar velocity profile in a rectangular duct with a 2.550×30 mm 2 cross section for six individual volume fractions ranging from Φ ≤ 0.01% to Φ = 19.9%. The uncertainty of the measured in-plane velocity was found to be σ u = 1.39% and σ u = 3.34% while the uncertainty for the out-ofplane velocity was σ w = 9.06% and σ w = 22.57% for Φ ≤ 0.01% and Φ = 19.9%, respectively. These uncertainties are higher compared to those in APTV with small particles and dilute suspensions. However, we are convinced that by using a improved illumination technique and the use of a 3D calibration the uncertainties can be further reduced. Furthermore we only used two calibration curves generated at Φ ≤ 0.01% and Φ = 19.9% for our interpolation technique. In future works the method could be improved by considering a higher number of volume fractions for generating interand extrapolated curves. In fact, preliminary tests (not shown here) indicate that it is benefitial to perform the calibration at solid volume fractions which are actually higher than those used during flow measurements. Finally, to date APTV has been scarcely applied to large fluorescent particles -following studies can help to gain further understanding of how particle size and properties of the optical system affect the measurement accuracy.