Calibration Routine for Quantitative Three-Dimensional Flow Field Measurements in Drying Polymer Solutions Subject to Marangoni Convection

Surface-tension induced flows may have a significant impact on the surface topography of thin films or small printed structures derived from polymer solution processing. Despite a century of research on Marangoni convection, the community lacks quantitative experimental flow field data, especially from within drying solutions. We utilize multifocal micro particle tracking velocimetry (μPTV) to obtain these data and show a calibration routine based on point spread function (PSF) simulations as well as experimental data. The results account for a varying sample refractive index, beneficial cover-glass correction collar settings as well as a multifocal lens system. Finally, the calibration procedure is utilized exemplarily to reconstruct a three-dimensional, transient flow field within a poly(vinyl acetate)-methanol solution dried with inhomogeneous boundary conditions.


Introduction
Creating homogeneous thin polymer films from solution is a key processing step in the production of coatings, adhesive tapes, displays, and printed organic electronic devices, such as OLEDs, solar cells, or biosensors.Besides well-established slot die coating processes for large area deposition, in the recent decade, inkjet printing for the selective application of small structures evolved from pure graphics application towards the deposition of functional materials [1].During the subsequent drying of small sessile structures, a hydrodynamic effect occurs, which transports the solute preferably towards the contact line, resulting in elevated edges of the deposit.Deegan et al. (1997) were the first to describe this so-called coffee-ring-effect (CRE) mathematically [2].Large area coatings may also suffer from deformations induced by inhomogeneities in the drying boundary conditions.Either lateral variations of the heat conductivity of the substrate [3] or lateral variations of the solvent mass transport in the gas phase above the drying film cause the liquid-gas interface to deform [4].
It is well established that in thin liquid films, buoyancy (i.e., Rayleigh-Bénard convection) is negligible and surface-tension driven convection is the dominant convective mass transport mechanism.A surface-tension gradient along the liquid gas interface results in a lateral flow from regions with low surface tension towards regions with high surface tension.This phenomenon is known as Marangoni convection and it may result in deformation of the liquid gas interface as well as in the formation of convection cells (Bénard-Marangoni convection [5,6]).The surface tension of a polymer solution depends on its temperature and composition at the free surface.While drying progresses, the temperature of the solution decreases due to the latent heat of evaporation, as does the solvent concentration within the film.This results in an increase in viscosity and ultimately in solidification of the film, including its potentially deformed free surface.Therefore, Marangoni-convection is undesirable while drying large area polymer films, but desirable in small printed structures to counteract the CRE.

Marangoni Convection in Thin Films
Although Marangoni convection has been investigated for more than a century, in combination with polymer film drying it remains a challenging task.In pure liquid films without evaporation, different geometries of convection cells [7], surface deformation due to (locally) heated substrates [8,9], or substrates with non-uniform thermal conductivity [10] have been investigated.The reviews by Oron et al. (1997), Craster and Matar (2009), and Gambaryan-Roisman (2015) give an overview of the experimental, theoretical as well as numerical work on the dynamics of liquid films [11][12][13].Their focus, however, lies on pure liquids, mainly without evaporation.Marangoni convection in pure liquid films with evaporation has been studied experimentally [14] as well as theoretically [15][16][17][18][19] and numerically [15].
Only a few authors have dealt with the complexity of Marangoni flows in drying polymer solutions.They experimentally investigated the influences of the solvent mixture composition [20], polymer molecular weight [21], initial polymer concentration [22,23], and film thickness [23] as well as the impact of the substrate inclination [24], evaporation rate [25], and local heating [26] on the shape and size of convection cell patterns.They all used either shadowgraphy or schlieren photography, which reveals lateral convection cell distributions; however, it does not provide quantitative data of the fluid flow.In a recent comprehensive review on multiphase Marangoni convection by Wang et al. (2016), the authors acknowledge the lack of work on polymer film drying and emphasize the need for new experiments providing quantitative data [27].

Mitigating the Coffee Ring Effect in Sessile Droplets by Means of Marangoni Convection
Following Deegan's work, research has derived numerous strategies controlling the CRE [28][29][30].Among these are methods inducing Marangoni convection within a sessile droplet.A shape change contributed to surface-tension induced convection cells by means of the latent heat of evaporation [31], utilizing different solvent-polymer combinations [32,33], in the presence of surfactants [34,35] or by using binary solvent mixtures with different boiling points [36][37][38] have been reported.These strategies appear to work in suspensions [31,[35][36][37] as well as in polymer solutions [32][33][34]38].Despite all experimental efforts, the transport mechanism of the CRE is still an ongoing matter of investigation.Only recently, an alternative theoretical model for the solute transport, which deviates from past assumptions, was proposed [39].

Measurement Techniques for Surface-Tension Induced Flows
In need of a clearer picture of Marangoni convection in drying polymer solutions, it is essential to acquire quantitative measurements of the flow field within the drying film or droplet.Therefore, several authors have utilized particle image velocimetry (PIV) and were able to reconstruct streamlines in films [23,40] and droplets [31,38,41,42], using sheet illumination to obtain vertical cross sections of the flow field.In the case of droplets, the obtained images had to be corrected due to the refraction of light at the curved droplet surface, resulting in a cropped field-of-view.Flow velocities have been calculated either by the cross-correlation of images [40,41] or by tracking a few selected particles to obtain representative velocities [23,31,38,42,43].This measurement technique, however, is limited to two dimensions.In principle, scanning the third dimension with consecutive measurements in different cross-sections could overcome this shortcoming.Due to the transient nature of the flow field, however, this is not feasible for solution drying.Furthermore, for PIV measurements, the sample must be optically accessible from two almost perpendicular directions.An illumination from above renders it virtually impossible to add defined boundary conditions for controlled drying.A brief discussion of further microfluidic measurement techniques can be found in our recent work [44] or in several comprehensive reviews [45][46][47].
In our group, we have established a measurement setup avoiding the shortcomings of previous PIV measurements.It is based on microscopic particle tracking velocimetry (µPTV) with only one optical access that was implemented using an inverse microscope and a transparent substrate, thus observing the flow field from below [4].The third velocity component of the tracer particles can be obtained from the convoluted optical response.This so-called point spread function (PSF) is shift variant towards the distance between the tracer particle and the focal plane of the microscope.A sharp image is obtained if a particle is in focus.With an increasing distance to the focal plane, a concentric ring structure with an increasing diameter can be observed.Correlating the ring diameter with the distance to the focal plane enables 3D measurements with only one camera [48,49].With an increasing distance between an observed particle and the focal plane, the signal-to-noise ratio of the observed diffraction ring deteriorates.Therefore, detection along the optical axis is limited.Recently, we have combined this approach with a multifocal microscopic setup.Equipped with beam splitters, multiple cameras, and individually motorized tube lens systems in front of each camera, we are able to adapt the observed focal plane for each camera individually, which allows us to expand the vertical field-of-view [44].Despite the observation from below, we utilize a fully temperature controlled substrate and no devices occupying the space above the sample, providing an opportunity for controlled convective drying.Therefore, our setup is ideal for the investigation of flow fields within drying polymer solutions.
The scope of this article is to demonstrate the calibration of such an optical setup for quantitative, transient, three-dimensional flow field measurements.We consider the influence of multifocal lens movements, refractive index of the sample solution as well as the choice of objective lenses and cover-glass correction.The evaluation routine is demonstrated on a simple drying experiment with inhomogeneous boundary conditions subject to Marangoni convection.

Materials and Methods
The measurement setup is the same as in our previous work [44].It consists of an inverse Nikon Ti fluorescence microscope, a CoolLED light source, and a MultiCam beam-splitter unit by Cairn Research Ltd (Faversham, United Kingdom).One camera (Andor Zyla 5.5 sCMOS, 2560 × 2160 px 2 , 6.5 µm/px) is connected to the backport of the microscope with fixed optics, whereas up to four additional cameras are connected via the MultiCam unit with a motorized tube-lens system.By adding or removing beam splitters, any number of cameras can be used simultaneously.An outline of the complete setup is shown in Figure 1a and the detailed optical setup of a single camera with a motorized lens is shown in Figure 1b.For brevity, only one camera is shown in full.
The objective lenses used were Nikon Plan Apo λ 60×/0.95 and Plan Fluor 100×/0.90.Both objectives have a correction collar for the cover glass thickness.Positioning of the objective is achieved either via the manual focus knob on the microscope or by a piezo positioner (Physik Instrumente PI GmbH, P-725.2CD,Karlsruhe, Germany).The motorized lens systems are custom-built by Cairn Research Ltd.The coordinate along the optical axis is defined as z with its origin at the interface between the glass substrate and the sample, if not stated otherwise.The actual position within the sample is given without accent, while objective positions are indicated with a tilde, z obj , having the same origin as z.The motorized lenses allow for changes to be made to the focal plane position, z f ocus , without moving the objective.For the lens position, s axis = 0 mm, all cameras observe the same focal plane.
Research Ltd (Faversham, United Kingdom).One camera (Andor Zyla 5.5 sCMOS, 2560 × 2160 px , 6.5 μm/px) is connected to the backport of the microscope with fixed optics, whereas up to four additional cameras are connected via the MultiCam unit with a motorized tube-lens system.By adding or removing beam splitters, any number of cameras can be used simultaneously.An outline of the complete setup is shown in Figure 1(a) and the detailed optical setup of a single camera with a motorized lens is shown in Figure 1 The lateral or transverse magnification, M T , is the ratio of the lateral size of the recorded image to the corresponding dimensions of the physical object, whereas the axial or longitudinal magnification, M L , is the ratio of the axial distance between two planes in the image space to the corresponding distance in the object space.Both magnifications are linked via Equation (1) for a thin lens [50].For off-the-shelf microscopes, M T is given on the objectives (60× resp.100×).In our µPTV setup, however, deviations may occur due to the additional motorized lens systems.Therefore, we determined the transverse magnification for different lens positions, s axis , with a 1/100 mm scale printed on a glass slide and the image-processing program, ImageJ.The longitudinal magnification can either be calculated from (1) or measured by performing a z-scan of objects with a known axial distance.
Fluorescent monodisperse tracer particles (FSDG003, Bangs Laboratories Inc., Fishers, IN, USA) with a 0.52 µm diameter were used to visualize the flow.The observed three-dimensional image of an ideal, dot-like object is known as the point-spread function (PSF).The shape of the PSF accounts for refraction and aberration effects and is characteristic for the respective optical system [51].It can be measured performing a z-scan on a tracer particle with a fixed position in a sample.With known optical properties, a PSF can be simulated using the well-established Gibson-Lanni model [52], which assumes an ideal point-like object.Its limitations, however, are its computation speed.Recently, a fast and yet accurate implementation was proposed [53].Table 1 gives an overview of the required optical input parameters.The design properties, marked with superscript *, are set by the objective manufacturer.Only the design cover-glass thickness, t * g , can be adjusted for objectives with a correction collar.As pointed out in previous work, an asymmetric PSF caused by spherical aberration is beneficial due to the unambiguous occurrence of ring structures on only one side of the focal plane [44].On the opposite side, the particle image blurs out and vanishes quickly with an increasing distance to the focal plane.Considering that the refractive index of polymer solutions is typically above 1.3, the mismatch in the refractive indices of the immersion and sample inevitably leads to an asymmetric PSF [52].Therefore, we used air-immersion objectives only.The mismatch also results in a displacement of the position of the focal plane [44,54], hence the differentiation between z and z obj .Due to the non-design conditions, even if a tracer particle is in focus, the observed image is not razor-sharp.Therefore, we defined the "best focus" of a PSF z f ocus as the lateral plane at which the axial intensity profile has its maximum.Figure 2 shows a vertical cross-section of a simulated PSF (a) and the corresponding axial intensity profile (b).The dashed line indicates the best focus.Although the particle position in the sample is set to be z p = 54 µm, the observed position is z obj = 33.5 µm.
Colloids Interfaces 2019, 3 FOR PEER REVIEW 5 out in previous work, an asymmetric PSF caused by spherical aberration is beneficial due to the unambiguous occurrence of ring structures on only one side of the focal plane [44].On the opposite side, the particle image blurs out and vanishes quickly with an increasing distance to the focal plane.
Considering that the refractive index of polymer solutions is typically above 1.3, the mismatch in the refractive indices of the immersion and sample inevitably leads to an asymmetric PSF [52].Therefore, we used air-immersion objectives only.The mismatch also results in a displacement of the position of the focal plane [44,54], hence the differentiation between  and ̃ .Due to the non-design conditions, even if a tracer particle is in focus, the observed image is not razor-sharp.Therefore, we defined the "best focus" of a PSF ̃ as the lateral plane at which the axial intensity profile has its maximum.Figure 2 shows a vertical cross-section of a simulated PSF (a) and the corresponding axial intensity profile (b).The dashed line indicates the best focus.Although the particle position in the sample is set to be  = 54 μm, the observed position is ̃ = 33.5 μm.The outmost ring is of the highest intensity and therefore is most suitable for tracking [48,49,55].From the PSF-simulation results, we obtained the outmost ring from the maximum of the vertical intensity profiles for each value of the radius,  , and smoothed the numerical errors using a polynomial fit of 4th degree.The dash-dotted line in Figure 2(a) shows the fit.The ring detection in the experimental data was performed using our GPU-enhanced detection algorithm presented in [44].
For calibration purposes, a sample with a known refractive index and known vertical position of tracer particles is needed.In Cavadini et al. 2018, we proposed a stack of transparent tape strips with tracer particles in between the layers [44].Similar calibration samples were prepared with tesafilm crystal clear ( = 1.472,Δ = 54 ± 3 μm; tesa SE, 57315, Norderstedt, Germany).This is, however, limited to one sample refractive index.To cover a wider range of  we prepared calibration samples by depositing tracer particles on two glass slides with spacers of a known thickness between them.The gap was filled with different fluids with varied refractive indices.The fluids used are air ( = 1.000, [56]), bidistilled water ( = 1.336, [57]; Carl Roth GmbH, 3478.1,Karlsruhe, Germany), and refractometer calibration oil ( = 1.557;Bellingham + Stanley Ltd., 90-235, Tunbridge Wells, United Kingdom).Unlike the tesa stacks, multiple particle positions could only be realized in separate calibration samples.The height of all individual layers was measured with a The outmost ring is of the highest intensity and therefore is most suitable for tracking [48,49,55].From the PSF-simulation results, we obtained the outmost ring from the maximum of the vertical intensity profiles for each value of the radius, r, and smoothed the numerical errors using a polynomial fit of 4th degree.The dash-dotted line in Figure 2a shows the fit.The ring detection in the experimental data was performed using our GPU-enhanced detection algorithm presented in [44].
For calibration purposes, a sample with a known refractive index and known vertical position of tracer particles is needed.In Cavadini et al. 2018, we proposed a stack of transparent tape strips with tracer particles in between the layers [44].Similar calibration samples were prepared with tesafilm crystal clear (n s = 1.472, ∆z = 54 ± 3 µm; tesa SE, 57315, Norderstedt, Germany).This is, however, limited to one sample refractive index.To cover a wider range of n s we prepared calibration samples by depositing tracer particles on two glass slides with spacers of a known thickness between them.The gap was filled with different fluids with varied refractive indices.The fluids used are air (n s = 1.000, [56]), bidistilled water (n s = 1.336, [57]; Carl Roth GmbH, 3478.1,Karlsruhe, Germany), and refractometer calibration oil (n s = 1.557;Bellingham + Stanley Ltd., 90-235, Tunbridge Wells, United Kingdom).Unlike the tesa stacks, multiple particle positions could only be realized in separate calibration samples.The height of all individual layers was measured with a digital dial gauge (Mitutoyo Europe GmbH, 543-561D, Neuss, Germany).Figure 3 shows schematic drawings of the different calibration sample architectures.
Colloids Interfaces 2019, 3 FOR PEER REVIEW 6 digital dial gauge (Mitutoyo Europe GmbH, 543-561D, Neuss, Germany).Figure 3 shows schematic drawings of the different calibration sample architectures.The calibration experiments as well as the PSF simulations yield a correlation between the objective position, ̃ , and the outmost ring radius,  , for a known and fixed particle position.In drying experiments, however, the objective position is known while the particle positions need to be determined from the detected ring sizes.The conversion from either the experimental calibration results or the simulations is done as follows: First, the ring radius for multiple known particle positions with otherwise constant optical properties (i.e.,  ,  ,  * ) is combined to one dataset as shown in Figure 4(a) as black solid lines.Second, the intersections with constant objective positions are calculated (black circles and dotted red lines, respectively).Third, the intersections for each value of ̃ are fitted as isolines in a plot  over  as shown in Figure 4(b).With a known objective position, measured ring-sizes, and otherwise constant optical properties, the particle positions can be determined from the polynomial fits of the isolines.For a comparison with other work, the particle distance to the focal plane, Δ, is calculated with: where  accounts for the focal displacement due to the refractive index mismatch, which will be discussed in Section 3.1.The calibration experiments as well as the PSF simulations yield a correlation between the objective position, z obj , and the outmost ring radius, r ring , for a known and fixed particle position.In drying experiments, however, the objective position is known while the particle positions need to be determined from the detected ring sizes.The conversion from either the experimental calibration results or the simulations is done as follows: First, the ring radius for multiple known particle positions with otherwise constant optical properties (i.e., n s , t g , t * g ) is combined to one dataset as shown in Figure 4a as black solid lines.Second, the intersections with constant objective positions are calculated (black circles and dotted red lines, respectively).Third, the intersections for each value of z obj are fitted as isolines in a plot z p over r ring as shown in Figure 4b.With a known objective position, measured ring-sizes, and otherwise constant optical properties, the particle positions can be determined from the polynomial fits of the isolines.For a comparison with other work, the particle distance to the focal plane, ∆z, is calculated with: where m d accounts for the focal displacement due to the refractive index mismatch, which will be discussed in Section 3.1.
Colloids Interfaces 2019, 3 FOR PEER REVIEW 6 digital dial gauge (Mitutoyo Europe GmbH, 543-561D, Neuss, Germany).Figure 3 shows schematic drawings of the different calibration sample architectures.The calibration experiments as well as the PSF simulations yield a correlation between the objective position, ̃ , and the outmost ring radius,  , for a known and fixed particle position.In drying experiments, however, the objective position is known while the particle positions need to be determined from the detected ring sizes.The conversion from either the experimental calibration results or the simulations is done as follows: First, the ring radius for multiple known particle positions with otherwise constant optical properties (i.e.,  ,  ,  * ) is combined to one dataset as shown in Figure 4(a) as black solid lines.Second, the intersections with constant objective positions are calculated (black circles and dotted red lines, respectively).Third, the intersections for each value of ̃ are fitted as isolines in a plot  over  as shown in Figure 4(b).With a known objective position, measured ring-sizes, and otherwise constant optical properties, the particle positions can be determined from the polynomial fits of the isolines.For a comparison with other work, the particle distance to the focal plane, Δ, is calculated with: where  accounts for the focal displacement due to the refractive index mismatch, which will be discussed in Section 3.1.2), including the experimental fit from [44].
A simple drying experiment as proposed in [4] was performed to demonstrate a quantitative flow field reconstruction.A solution of poly(vinyl acetate) (PVAc, Carl Roth, 9154.1) and methanol (MeOH, Carl Roth, 4627.1) with tracer particles and an initial solvent load of X 0 = 2 g MeOH /g PVAc was blade-coated onto a transparent glass substrate (t g = 144 µm, T = 20 • C) and observed with µPTV from below.Inhomogeneous drying conditions were induced by partially covering the drying film (see Figure 1b).A Marangoni flow from the covered area towards the uncovered region was expected [4].The 60× objective was used with its correction collar set to t * g = 110 µm.The objective position was set to z obj = 45 µm with the piezo positioner and two cameras with s axis,A = 0 mm and s axis,B = 122 mm recorded 10 s of particle movements with 50 frames per second (fps).The delay between the coating and the start of the recording was approximately 5 s.
The refractive index of the sample is required for quantitative analysis of the recordings.In [58], a similar drying experiment with the same material system was performed.Instead of a partial cover, drying was controlled with lateral forced convection in a flow channel with 0.2 m/s air velocity.Simultaneously, the change in the solvent load within the film was measured over time and the film height by inverse Raman spectroscopy.The data were used to estimate the sample refractive index by applying the mixing rule from [59]: with ϕ being the volume fraction.The drying rate due to the forced air convection with 0.2 m/s is of a similar magnitude as the drying rate due to free convection [60].This implies a similar decrease in the solvent concentration in the uncovered area of our partially covered drying experiment.The covered part, however, has a significantly lower drying rate.For a conservative estimate, a solvent load between X t=0 s = 2.0 g MeOH /g PVAc and X t=35 s = 1.0 g MeOH /g PVAc was used to calculate the sample refractive index.The required material properties are listed in Table 2. 1.32843 (20 • C) [63]   1 As specified by the manufacturer.
The diffraction rings in the recorded series of images were detected using the GPU-enhanced algorithm proposed in [44].Linking of individual particles to trajectories was performed using an algorithm from [64].To reduce experimental noise, the trajectories were smoothed by applying a Savitzky-Golay filter [65] with a window length and polynomial degree of 7 and 1, respectively.

Focal Displacement Calibration
Initially, the focal displacement due to a mismatch between the immersion and sample refractive index was investigated.Calibration samples with different refractive indices were prepared and z-scans were carried out.A series of simulations with z p = 0-200 µm was performed for each sample using the optical properties of the objective and the measured substrate thickness as input parameters.Figure 5 shows a comparison of the best focus.The simulations (lines) show a very good agreement with the experimental results (markers) over a wide range of n s , as do the data from our previous work [44].The impact of the sample refractive index on the displacement can be determined from the slope, m d , of the data.The following fit is based on simulations in the range of n s = 1.0-1.7:

Experimental Calibration of Motorized Lens System
The multifocal setup introduces motorized lens systems uncommon in standard of-the-shelf microscopes.To compare the experimental ring sizes with the calculated PSFs, the transverse magnification was determined for different positions of the motorized lens.The experimental results indicate that even for  = 0 mm, the measured  is slightly smaller than the nominal one imprinted on the objectives (Figure 6(a)).Furthermore, it can be clearly seen that the optical system is not telecentric, since the transverse magnification decreases with an increasing axial lens position.

Experimental Calibration of Motorized Lens System
The multifocal setup introduces motorized lens systems uncommon in standard of-the-shelf microscopes.To compare the experimental ring sizes with the calculated PSFs, the transverse magnification was determined for different positions of the motorized lens.The experimental results indicate that even for s axis = 0 mm, the measured M T is slightly smaller than the nominal one imprinted on the objectives (Figure 6a).Furthermore, it can be clearly seen that the optical system is not telecentric, since the transverse magnification decreases with an increasing axial lens position.

Experimental Calibration of Motorized Lens System
The multifocal setup introduces motorized lens systems uncommon in standard of-the-shelf microscopes.To compare the experimental ring sizes with the calculated PSFs, the transverse magnification was determined for different positions of the motorized lens.The experimental results indicate that even for  = 0 mm, the measured  is slightly smaller than the nominal one imprinted on the objectives (Figure 6(a)).Furthermore, it can be clearly seen that the optical system is not telecentric, since the transverse magnification decreases with an increasing axial lens position.The main purpose of the motorized lens system, however, is to change the focal plane position for each camera individually.Therefore, multiple calibration samples were scanned with different lens positions using the 60× objective.The best focus of each particle layer is depicted in Figure 7. Different color shades indicate individual samples.Data for each particle layer are connected with dashed gray lines for better readability.Each sample has an arbitrary origin of z f ocus , which was kept constant for all lens positions.Therefore, only the slope of each dataset is of relevance.As shown, it is virtually identical for different particle layers within one sample as well as for different sample media, indicating that neither influences the focal shift due to the motorized lens.
The main purpose of the motorized lens system, however, is to change the focal plane position for each camera individually.Therefore, multiple calibration samples were scanned with different lens positions using the 60× objective.The best focus of each particle layer is depicted in Figure 7. Different color shades indicate individual samples.Data for each particle layer are connected with dashed gray lines for better readability.Each sample has an arbitrary origin of ̃ , which was kept constant for all lens positions.Therefore, only the slope of each dataset is of relevance.As shown, it is virtually identical for different particle layers within one sample as well as for different sample media, indicating that neither influences the focal shift due to the motorized lens.An analytic solution can be derived from Equation (1).Considering that Δ̃ is an axial distance in the object space and Δ is an axial displacement in the image space, the quotient is equivalent to the axial magnification,  (Equation ( 5)).The negative sign needs to be introduced since an increase in ̃ decreases the distance between the objective lens and the observed position [50].Inserting the linear relation from Figure 6 (Equation ( 6)) and Equation ( 1) in ( 5), the subsequent integration yields Equation (7).
By shifting the origin of ̃ of the experimental datasets from Figure 7, the data coincide as depicted in Figure 6(b).The dataset for the 100× objective was derived accordingly.The integration constant in (7) is fitted for each objective using the sum of least squares.For the 100× objective, the theory and experiment are in excellent agreement.For 60×, however, the theoretic solution deviates An analytic solution can be derived from Equation (1).Considering that ∆ z f ocus is an axial distance in the object space and ∆s axis is an axial displacement in the image space, the quotient is equivalent to the axial magnification, M L (Equation ( 5)).The negative sign needs to be introduced since an increase in z f ocus decreases the distance between the objective lens and the observed position [50].Inserting the linear relation from Figure 6 (Equation ( 6)) and Equation ( 1) in ( 5), the subsequent integration yields Equation (7).
By shifting the origin of z f ocus of the experimental datasets from Figure 7, the data coincide as depicted in Figure 6b.The dataset for the 100× objective was derived accordingly.The integration constant in ( 7) is fitted for each objective using the sum of least squares.For the 100× objective, the theory and experiment are in excellent agreement.For 60×, however, the theoretic solution deviates significantly from the experimental results.The reason for the deviation is most likely Equation (1), which is rigorously valid only for single thin lenses and not for complex lens systems as in microscope objectives.A linear fit is better suited to match the experimental data.

Diffraction-Ring Size Calibration for Off-Focus Particle Positions
Up to this point, only the best focus of the particles was considered.The axial position of particles is, however, determined from diffraction-ring sizes occurring whenever a particle is not in focus.With known transverse magnification, the radii of experimental PSFs derived from z-scans can be converted from image pixels to micrometers.A comparison with a simulated PSF is shown as a vertical cross-section in Figure 8. Overall, it shows a very good agreement despite the fact that the Gibson-Lanni simulation model assumes an ideal point source and the tracer particle has a diameter well above the pixel resolution limit of the cameras.With M T (s axis = 0 mm) = 57.5, the tracer particle diameter is d p = 4.6 px and is clearly not an ideal point.As seen in Figure 8, this has no apparent effect on the PSF.

Diffraction-Ring Size Calibration for Off-Focus Particle Positions
Up to this point, only the best focus of the particles was considered.The axial position of particles is, however, determined from diffraction-ring sizes occurring whenever a particle is not in focus.With known transverse magnification, the radii of experimental PSFs derived from -scans can be converted from image pixels to micrometers.A comparison with a simulated PSF is shown as a vertical cross-section in Figure 8. Overall, it shows a very good agreement despite the fact that the Gibson-Lanni simulation model assumes an ideal point source and the tracer particle has a diameter well above the pixel resolution limit of the cameras.With  ( = 0 mm) = 57.5, the tracer particle diameter is  = 4.6 px and is clearly not an ideal point.As seen in Figure 8, this has no apparent effect on the PSF.For quantitative validation of the simulation model regarding the ring sizes, a -scan on a tesa stack was performed.The best focus and rings were detected via the vertical intensity profile and GPU-enhanced detection algorithm, respectively.Figure 9 shows the experimental results (markers) as well as the simulated ring sizes (red lines).The detection algorithm yields a high number of falsepositive rings (blue crosses) for very small radii only.The deposition of the particle layers for the calibration samples was done from aqueous solution.The dry layers show a high amount of agglomerates, which tend to have a much higher intensity, resulting in a less pronounced blur-out.Drying experiments with polymer solutions do not suffer from tracer-particle agglomeration.Both the false positives as well as the detected agglomerates were sorted out manually, but are shown in Figure 9 for full disclosure.The quality of the detection algorithm, however, is not in the scope of this article.
Since experimental -scans have an arbitrary origin, the detected best focus positions (dashed horizontal lines) were used to vertically shift the experimental results by calculating the sum of least squares with the best focus positions derived from simulations (dotted horizontal lines).Neither the distances between the individual particle layers nor the rings were fitted in any way.Figure 9 shows For quantitative validation of the simulation model regarding the ring sizes, a z-scan on a tesa stack was performed.The best focus and rings were detected via the vertical intensity profile and GPU-enhanced detection algorithm, respectively.Figure 9 shows the experimental results (markers) as well as the simulated ring sizes (red lines).The detection algorithm yields a high number of false-positive rings (blue crosses) for very small radii only.The deposition of the particle layers for the calibration samples was done from aqueous solution.The dry layers show a high amount of agglomerates, which tend to have a much higher intensity, resulting in a less pronounced blur-out.Drying experiments with polymer solutions do not suffer from tracer-particle agglomeration.Both the false positives as well as the detected agglomerates were sorted out manually, but are shown in Figure 9 for full disclosure.The quality of the detection algorithm, however, is not in the scope of this article.
Since experimental z-scans have an arbitrary origin, the detected best focus positions (dashed horizontal lines) were used to vertically shift the experimental results by calculating the sum of least squares with the best focus positions derived from simulations (dotted horizontal lines).Neither the distances between the individual particle layers nor the rings were fitted in any way.Figure 9 shows an excellent agreement between the experiments and simulations regarding the ring sizes for multiple settings of the cover-glass correction-collar (t * g ).The actual substrate thickness was t g = 144 µm.For a design value much larger than the actual value (see Figure 9c), the diffraction rings appear above the respective focal plane for the lower two particle layers and below for the upper two layers.This results in an area where a detected ring size is ambiguous regarding its corresponding particle position (i.e., the intersection of the red lines).For quantitative measurements, this ambiguity is unwanted and needs to be avoided.detection is impossible.Considering a larger absolute slope, as observed for  = 110 μm, a ring with a constant diameter would have a larger distance to the corresponding focal plane if compared to  * = 150 μm .Thus, if a large vertical observation volume is desired, a large absolute slope is beneficial.Overall, the correction-collar setting is a trade-off between the vertical extent of the observed volume and the error tolerance of ring size detection.We found that  * = 110 μm was most suitable for our application.In [44], we reported data on ring-size calibration, plotted using the distance between the particle position and the focal plane.The simulated data in Figure 4(c) were derived using matching optical properties.A good agreement between the current simulations and our old data can be seen, but only for a single, large  position.Similar ring-size data can be found in [48,49,55,66].Speidel et al. (2003) reported a linear relation between the ring-size and the distance to the focal plane, but they provided the ring radius in arbitrary units only [48].The other authors reported non-linear relations [49,55,66].Qualitatively, their data and shape of the plot compares well with our findings.However, a quantitative comparison with PSF simulations was not possible due to the lack of optical properties of the respective setups.Comparing the remaining two correction-collar settings, a change in the ring data's slope is apparent (Figure 9a,b).Considering a smaller absolute slope, as observed for t * g = 150 µm, this setting would be beneficial if large errors in the detected ring radii are to be expected.A deviating ring size would have a less pronounced effect on the derived vertical particle position.The signal-to-noise ratio of the recorded rings decreases with an increasing ring size up to a point where ring detection is impossible.Considering a larger absolute slope, as observed for t * g = 110 µm, a ring with a constant diameter would have a larger distance to the corresponding focal plane if compared to t * g = 150 µm.Thus, if a large vertical observation volume is desired, a large absolute slope is beneficial.Overall, the correction-collar setting is a trade-off between the vertical extent of the observed volume and the error tolerance of ring size detection.We found that t * g = 110 µm was most suitable for our application.In [44], we reported data on ring-size calibration, plotted using the distance between the particle position and the focal plane.The simulated data in Figure 4c were derived using matching optical properties.A good agreement between the current simulations and our old data can be seen, but only for a single, large z f ocus position.Similar ring-size data can be found in [48,49,55,66].Speidel et al. (2003) reported a linear relation between the ring-size and the distance to the focal plane, but they provided the ring radius in arbitrary units only [48].The other authors reported non-linear relations [49,55,66].Qualitatively, their data and shape of the plot compares well with our findings.However, a quantitative comparison with PSF simulations was not possible due to the lack of optical properties of the respective setups.

Flow Field of Partially Covered Drying Experiment
Recordings of both cameras from the drying experiment are provided in the Supplementary Materials.The dataset was evaluated as follows: Using the objective and lens positions and the linear fit in Figure 6b, the calculated focal planes were at z obj,A, f ocus = 45.0 µm and z obj,B, f ocus = 8.3 µm for the two cameras, respectively.The estimated sample refractive index was n s = 1.374 ± 0.011 during the observation period, including variations due to a vertical concentration gradient.With these input parameters, PSF simulations were performed and polynomial fit functions were derived as described in Figure 4a,b.Figure 10 shows the resulting lookup data, accounting for the uncertainty in n s and the focal planes of the two cameras.The transverse magnifications are M T,A = 57.5 and M T,B = 50.8,respectively.The gray area indicates the vertical field of view.It results from the minimal and maximal observed ring radius in the recordings.It is noteworthy that the significant concentration range the film undergoes within the observed timeframe results in only a small uncertainty of n s .This is important, since n s has a major impact on z p calibration (see Figure 5).The range of uncertainty of n s is shown with black lines in Figure 10.The maximal error in the vertical particle position for the given experiment resulting from variations in n s are ∆z p,A,max = ±0.7 µm and ∆z p,B,max = ±0.5 µm.∆z p increases with increasing r ring and increasing z obj .
, , , , for the two cameras, respectively.The estimated sample refractive index was  = 1.374 ± 0.011 during the observation period, including variations due to a vertical concentration gradient.With these input parameters, PSF simulations were performed and polynomial fit functions were derived as described in Figure 4(a,b).Figure 10 shows the resulting lookup data, accounting for the uncertainty in  and the focal planes of the two cameras.The transverse magnifications are  , = 57.5 and  , = 50.8,respectively.The gray area indicates the vertical field of view.It results from the minimal and maximal observed ring radius in the recordings.It is noteworthy that the significant concentration range the film undergoes within the observed timeframe results in only a small uncertainty of  .This is important, since  has a major impact on  calibration (see Figure 5).The range of uncertainty of  is shown with black lines in Figure 10.The maximal error in the vertical particle position for the given experiment resulting from variations in  are Δ , , = ±0.7 μm and Δ , , = ±0.5 μm.Δ increases with increasing  and increasing ̃ .The resulting three-dimensional particle trajectories for both observation volumes are shown in Figure 11.The black arrows at the tip of each trajectory indicate the flow direction and velocity.It can be clearly seen that the dominant current is from the covered area towards the uncovered area.The colors indicate the temporal occurrence of particles.In this diagram, only the relative velocity of two trajectories with a similar length can be estimated by comparing the color gradient.A relatively smaller gradient on a trajectory of a similar length allows a conclusion on a higher velocity.The resulting three-dimensional particle trajectories for both observation volumes are shown in Figure 11.The black arrows at the tip of each trajectory indicate the flow direction and velocity.It can be clearly seen that the dominant current is from the covered area towards the uncovered area.The colors indicate the temporal occurrence of particles.In this diagram, only the relative velocity of two trajectories with a similar length can be estimated by comparing the color gradient.A relatively smaller gradient on a trajectory of a similar length allows a conclusion on a higher velocity.
For quantitative assessment of the velocity, Savitzky-Golay filtering using the first derivate was applied to each trajectory.The results were averaged in slices of ∆z = 2 µm and for each cartesian velocity component.Figure 12 shows the vertical velocity distribution.The highest mean velocity occurs in the positive x direction, from the covered to uncovered area.In (a), a vertical x-velocity gradient can be seen.This is to be expected for Marangoni convection, which is an interfacial flow.The viscosity of the polymer solution and the non-slip condition at the substrate results in the observed gradient.The in-plane velocity, u y , perpendicular to x is shown in (b).The mean velocity is small in the lower observation volume and scattered in the upper volume.The broad fluctuations at around z = 70 µm indicate that averaging over the complete xy-domain might not be well suited if several different flow domains occur.Furthermore, it should be considered that tracer particles only provide Lagrangian information on the flow field.Considering the two upmost data points in (b), they indicate a negative mean y-velocity with moderate fluctuations.Considering the small number of averaged particles in these slices, as depicted in (c), the reason could be an uneven distribution or insufficient number of tracer particles.This may lead to a misinterpretation of the obtained data.An in-depth analysis of the flow regimes, however, is not in the scope of this article.The highest observed total velocity from quantitative µPTV is u = 84 µm/s.For quantitative assessment of the velocity, Savitzky-Golay filtering using the first derivate was applied to each trajectory.The results were averaged in slices of Δ = 2 μm and for each cartesian velocity component.Figure 12 shows the vertical velocity distribution.The highest mean velocity occurs in the positive  direction, from the covered to uncovered area.In (a), a vertical -velocity gradient can be seen.This is to be expected for Marangoni convection, which is an interfacial flow.The viscosity of the polymer solution and the non-slip condition at the substrate results in the observed gradient.The in-plane velocity,  , perpendicular to  is shown in (b).The mean velocity is small in the lower observation volume and scattered in the upper volume.The broad fluctuations at around  = 70 μm indicate that averaging over the complete -domain might not be well suited if several different flow domains occur.Furthermore, it should be considered that tracer particles only provide Lagrangian information on the flow field.Considering the two upmost data points in (b), they indicate a negative mean y-velocity with moderate fluctuations.Considering the small number of averaged particles in these slices, as depicted in (c), the reason could be an uneven distribution or insufficient number of tracer particles.This may lead to a misinterpretation of the obtained data.An in-depth analysis of the flow regimes, however, is not in the scope of this article.The highest observed total velocity from quantitative µPTV is  = 84 μm/s.

Discussion
It has been pointed out that despite more than a century of research on Marangoni convection, there is still a need for more quantitative experimental data.Past research efforts have been focused on convection patterns in pure liquids, omitting the complexity of superimposed solution drying.While often unwanted in functional thin films, surface deformation due to Marangoni flows may lead to homogeneous inkjet-printed patterns.The recently established measurement technique, µPTV, enables transient three-dimensional flow field measurements within thin films and small printed structures without obstructing the space above the sample.This facilitates the implementation of controlled drying boundary conditions.The calibration of such a setup is strongly dependent on optical parameters, such as the refractive index of the sample or the substrate glass thickness.We were able to demonstrate that our experimental calibration data are in very good agreement with simulated data from an existing, but recently accelerated, mathematical model, which relies only on the available data of the optical setup.Consequently, we derived polynomial functions from the simulation results, which can be used to match experimentally observed diffraction-ring sizes with vertical positions in the sample.In past research, experimental correlations between the ring diameter and the distance to the focal plane have been reported [48,49,55].Due to a lack of data on the respective optical setups, a quantitative comparison could not be performed.Furthermore, we are

Discussion
It has been pointed out that despite more than a century of research on Marangoni convection, there is still a need for more quantitative experimental data.Past research efforts have been focused on convection patterns in pure liquids, omitting the complexity of superimposed solution drying.While often unwanted in functional thin films, surface deformation due to Marangoni flows may lead to homogeneous inkjet-printed patterns.The recently established measurement technique, µPTV, enables transient three-dimensional flow field measurements within thin films and small printed structures without obstructing the space above the sample.This facilitates the implementation of controlled drying boundary conditions.The calibration of such a setup is strongly dependent on optical parameters, such as the refractive index of the sample or the substrate glass thickness.We were able to demonstrate that our experimental calibration data are in very good agreement with simulated data from an existing, but recently accelerated, mathematical model, which relies only on the available data of the optical setup.Consequently, we derived polynomial functions from the simulation results, which can be used to match experimentally observed diffraction-ring sizes with vertical positions in the sample.In past research, experimental correlations between the ring diameter and the distance to the focal plane have been reported [48,49,55].Due to a lack of data on the respective optical setups, a quantitative comparison could not be performed.Furthermore, we are the first to report a deliberate mismatch in the cover-glass correction collar settings to increase the axial field-of-view.The dependency of the lookup data on the sample refractive index is significant.Despite obvious concentration changes while drying the PVAc-MeOH solution at near ambient conditions, the refractive index within the observation period changes only slightly.This may vary for different material systems and drying conditions.
In previous work, PIV measurements resulted in two-dimensional qualitative flow field data where either only size and form of the convective cells were analyzed or velocities were derived only for a selected few representative particles.With a simple drying experiment, we demonstrated that µPTV measurements grant access to transient, three-dimensional microscopic flow fields, resulting in quantitative velocity data over a significant observation volume.As expected for a partially covered PVAc-MeOH film, the dominant flow occurs in the x-direction from the covered towards the uncovered area due to Marangoni convection (see Figure 12a).The x-velocity profile in the lower observation volume closer to the substrate appears to be linear and can be fitted accordingly as shown with a dashed blue line.Such a linear profile strongly resembles the flow profile of a planar Couette-flow, where a Newtonian fluid undergoes a shear-driven flow between a horizontally moving upper wall and a stationary lower wall.The resulting shear rate is ∆u x /∆z = 1.3 s −1 , which is very low and justifies the assumption of Newtonian behavior for the PVAc-MeOH solution.Instead of a moving upper wall, however, the horizontal Marangoni flow from the covered towards the uncovered area acts as the driving force.
The velocities in the upper observation volume strongly deviate from a Couette-flow profile.Especially in the y-direction in the film plane, but perpendicular to the dominant x-flow, fluctuations increase drastically.This indicates a secondary instability, which needs to be investigated further in the near future.Errors induced by tracer particle sedimentation, inertia, and Brownian motion are discussed in detail in Appendix A. The sedimentation velocity is less than 1 µm/h while the observation period is 10 s.Therefore, it is safe to assume that this has no noticeable impact on the results.However, Brownian motion may result in a velocity fluctuation of close to 1 µm/s.Averaging multiple particle velocities significantly mitigates these undirected fluctuations.An axial extension of the observation volume would be beneficial and could be achieved by employing a high-power laser illumination.

Conclusions
Micro particle tracking velocimetry (µPTV) enables microfluidic flow-field measurements with only one optical access through a transparent substrate.The multifocal approach as well as the diffraction ring analysis allows for the detection of line-of-site velocities, resulting in a quantitative three-dimensional measurement technique.In this article, we have presented an in-depth calibration and evaluation routine for the line-of-site particle positions dependent on the optical parameters of the setup and the sample.Our experimental calibration data were in good agreement with a well-established simulation model.The unobstructed space above the sample makes µPTV ideal for the investigation of surface-tension driven convection in thin films and small printed structures, since the drying conditions can be adapted.This, as well as the evaluation routine, was successfully demonstrated on a simple drying experiment.In a partially covered thin film of poly(vinyl acetate)-methanol drying at near ambient conditions, a dominant flow from the covered towards the uncovered area was measured.The highest observed velocity with µPTV was u = 84 µm/s.Fluctuations perpendicular to the dominant flow were detected and demand further investigations.Overall, we consider µPTV to be a useful addition to the existing portfolio of microfluidic measurement techniques.
with k B and T being the Boltzmann constant and the film temperature, respectively.The mean squared displacement due to diffusion in one Cartesian direction, s = {x, y, z}, is: with ∆t being the time between two frames.This results in a particle displacement of s 2 = 0.019 µm.
With ∆t = 0.02 s, the absolute diffusion velocity in one dimension is less than 1 µm/s.By comparing the displacement due to diffusion with the displacement due to the flow field in one Cartesian direction, ∆s = u s ∆t, the relative error due to Brownian motion can be assessed by: For a single particle with a velocity, u s < of 10, 50, or 80 µm/s, the resulting error is 9.4, 4.7, or 1.2%, respectively.The error after averaging over multiple particles, N, is ε B / √ N for unidirectional statistical errors [69].We set N = 7 to be the window length of the applied Savitzky-Golay filter.Consequently, for a particle trajectory derived from at least seven particle positions with a velocity of 10, 50, or 80 µm/s, the resulting maximal error due to Brownian motion reduces to 3.5, 1.8, or 0.4%, respectively.The component velocities in Figure 12 were, in most cases, averaged over or more particles, which reduces the error even further.Considering the total velocity in all three dimensions, Equation (A4) would change to s 2 = 6 D p ∆t and the total diffusion velocity would be 1.6 µm/s.

Figure 1 .
Figure 1. Outline of micro particle tracking velocimetry (µPTV).The objective lens is attached to a piezo positioner.The signal is distributed to up to five cameras using beam splitters.Motorized lens systems precede four out of five cameras.(a) The complete setup.(b) Complete optical setup of a single camera.For reasons of clarity, only one camera system is depicted in full.

Figure 3 .
Figure 3. Schematic drawings of the calibration samples used for detecting tracer particles with known positions and a known sample refractive index.(Left) Tesa stack; (Right) Glass slides with spacers and calibration fluid.

Figure 4 .Figure 3 .
Figure 4. Conversion steps to obtain particle-position fit functions from simulations: (a) Ring sizes for various particle positions ( ) and intersections for constant objective positions.(b) Particle positions as a polynomial function of the ring sizes for various objective positions derived from intersections in

Figure 3 .
Figure 3. Schematic drawings of the calibration samples used for detecting tracer particles with known positions and a known sample refractive index.(Left) Tesa stack; (Right) Glass slides with spacers and calibration fluid.

Figure 4 .
Figure 4. Conversion steps to obtain particle-position fit functions from simulations: (a) Ring sizes for various particle positions ( ) and intersections for constant objective positions.(b) Particle positions as a polynomial function of the ring sizes for various objective positions derived from intersections in

Figure 4 .
Figure 4. Conversion steps to obtain particle-position fit functions from simulations: (a) Ring sizes for various particle positions (z p ) and intersections for constant objective positions.(b) Particle positions as a polynomial function of the ring sizes for various objective positions derived from intersections in (a).(c) Particle distance to the focal plane depending on the ring sizes and focal plane position derived from (b) with Equation (2), including the experimental fit from[44].

Figure 5 .
Figure 5. Focal displacement due to mismatch in the refractive indices of immersion (air,  = 1.000) and sample media ( ).Points indicate experimental values from the calibration samples; lines indicate simulations.

Figure 6 .Figure 5 .
Figure 6.Experimental calibration of the magnification.(a) Dependency of the transverse magnification ( ) on the position of the motorized lens ( ) for two objectives with different nominal magnification.Experimental values and linear fit.(b) Axial focal displacement as a function

Figure 5 .
Figure 5. Focal displacement due to mismatch in the refractive indices of immersion (air,  = 1.000) and sample media ( ).Points indicate experimental values from the calibration samples; lines indicate simulations.

Figure 6 .Figure 6 .
Figure 6.Experimental calibration of the magnification.(a) Dependency of the transverse magnification ( ) on the position of the motorized lens ( ) for two objectives with different nominal magnification.Experimental values and linear fit.(b) Axial focal displacement as a function Figure 6.Experimental calibration of the magnification.(a) Dependency of the transverse magnification (M T ) on the position of the motorized lens (s axis ) for two objectives with different nominal magnification.Experimental values and linear fit.(b) Axial focal displacement as a function of the motorized lens position.Theoretical trend for thin lenses (dotted line), experimental results (points), and linear fit (dashed/dash dotted lines).

Figure 7 .
Figure 7. Axial focal displacement as a function of the motorized lens position derived from the -scans of calibration samples with a 60× objective.Different color shades indicate individual samples.Data for each particle layer are connected with dashed gray lines for better readability.Each sample has an arbitrary origin of ̃ , which remains constant within individual datasets.

Figure 7 .
Figure 7. Axial focal displacement as a function of the motorized lens position derived from the z-scans of calibration samples with a 60× objective.Different color shades indicate individual samples.Data for each particle layer are connected with dashed gray lines for better readability.Each sample has an arbitrary origin of z f ocus , which remains constant within individual datasets.

Figure 8 .
Figure 8. Vertical cross-section of the PSF from a single tracer particle.Parameters: z p = 54 µm, n s = 1.472,N A = 0.95, t * g = 150 µm, t g = 144 µm.(a) Experimental result from a z -scan on a tesa calibration sample.Cross-section image stitched with ImageJ.(b) Simulation with fast Gibson-Lanni model implementation [53].Dashed line indicates the best focus at z obj,exp = 34 µm and z obj,sim = 33.5 µm, respectively.

Figure 9 .
Figure 9.Comparison of experimental and simulated diffraction ring sizes in a tesa stack for different cover-glass correction-collar settings,  * .Experimental data were shifted vertically to match best with the simulated focal plane positions.Ring sizes were not fitted in any way.

Figure 9 .
Figure 9.Comparison of experimental and simulated diffraction ring sizes in a tesa stack for different cover-glass correction-collar settings, t * g .Experimental data were shifted vertically to match best with the simulated focal plane positions.Ring sizes were not fitted in any way.

Figure 11 .
Figure 11.Detected particle trajectories in PVAc-MeOH during a drying experiment with partially covered film.Tracks with six or less particle positions observed in consecutive frames were omitted for reasons of clarity.Black arrows at the end of each trajectory indicate the flow direction smoothed with the Savitzky-Golay filter.

Figure 11 .Figure 12 .
Figure 11.Detected particle trajectories in PVAc-MeOH during a drying experiment with partially covered film.Tracks with six or less particle positions observed in consecutive frames were omitted for reasons of clarity.Black arrows at the end of each trajectory indicate the flow direction smoothed with the Savitzky-Golay filter.Colloids Interfaces 2019, 3 FOR PEER REVIEW 14

Figure 12 .
Figure 12.Cartesian velocity components from a partially covered drying experiment.The velocities of multiple particles were averaged in slices of ∆z = 2 µm.Error bars indicate the standard deviation.(a) Mean x -velocity component as a function of the vertical position.Dashed blue line indicates a linear fit resembling a Couette-flow profile.(b,d) Mean y -and z -velocity component, respectively, as a function of the vertical position.(c) Number of particles, N p , averaged in each slice.