Emulsion Flow Analysis of a Sensor Probe for Sustainable Machine Operation

Working fluids possess several applications in manufacturing processes, for instance lubricants in metals machining. Typical metal working fluids are formulated as oil-in-water emulsions. The maintenance of the physical stability of the working fluid during operation is a key factor for the sustainability of the relevant process. Therefore, continuous control of the working fluids stability and performance during machine operation is an essential tool for maintenance of the process performance. Turbidity measurement (TM) is a process technique for emulsion stability and quality assessment, where light transmission and absorption of an emulsion system is analyzed. However, for in-process measurement and fluid quality detection during the machine operation by TM, it is necessary to implement a transmission inline sensor into the working fluid flow line. The continuous flow measurement may cause problems for long-term sensor operation regarding, e.g., biofouling of the sensor optical glasses or erroneous measurements due to emulsion droplets segregation effects. In the present investigation, computational fluid dynamic (CFD) simulations have been adapted to obtain the emulsion flow conditions within a typical TM sensor probe, thereby allowing an assessment of the adhesion probability of microorganisms as well as droplet segregation effects. The simulation results indicate some temporal changes of the dispersed phase concentration in the detected emulsion flow. Due to droplet segregation in the emulsion, the flow velocity needs to exceed a certain value for reliable operation. It is shown here that in this flow regime microbiological attachments on the probe surfaces may be sufficiently avoided. A minimum critical flow velocity is derived to avoid biomolecule adhesion and thus durable operation of the sensor.


Introduction
Emulsions as dispersed multiphase systems are used as metal working fluids (MWFs) in material processing operations for lubrication and heat removal purposes [1].For a typical oil-in-water formulation of an MWF emulsion, the dispersed oil concentration amounts to 2-10 v/v% where the droplets have mean sizes of 0.1-2.0µm [2].MWF typically consist of different (mineral) oil mixtures and chemical additives to stabilize the emulsion, inhibit corrosion, kill bacteria via biocides, and avoid foam formation.These ingredients increase the performance of the MWF.The use of MWF in machining operation decreases mechanical stresses caused by shearing and friction in the contact zone between the tool and work piece.MWF further washes away chips and fines from the nascent metal Fluids 2017, 2, 9 2 of 10 surface.This way, rewelding of the newly formed surface is a preventative measure and the wetted surface is protected [3].MWF operation cycles within machining processes may last over several weeks and months.Due to metal working machining operation stresses, the MWF composition and morphology may change with time, resulting, e.g., in increased mean droplet size and broader droplet size distributions (DSDs), which may impact the MWF performance and, therefore, the machine operation efficiency and sustainability ( [1,4]).
In-situ and in-process emulsion quality measurements are to be applied for monitoring the change of the physical stability of the working fluid, where turbidimetry measurements (TM) has recently been successfully applied in different emulsion systems ( [5][6][7]).The turbidity spectrum of a MWF emulsion is influenced by the chemical composition, droplet size (distribution), internal structure, and concentration of the dispersed phase ( [5,8]).Measurement of the turbidity of an emulsion system requires contacting the fluid flow and optical access to the fluid.Therefore, during machine operation, sensor probes need to be dipped into the stream of the flowing MWF emulsion.The optical path through the liquid necessary for the measurement principle may be realized by an internal reflection arrangement.Thereby, the emulsion flows through a small channel within the sensor probe where a focused light beam passes the stream twice, forward and backward (see Figure 1).
MWF composition and morphology may change with time, resulting, e.g., in increased mean droplet size and broader droplet size distributions (DSDs), which may impact the MWF performance and, therefore, the machine operation efficiency and sustainability ( [1,4]).
In-situ and in-process emulsion quality measurements are to be applied for monitoring the change of the physical stability of the working fluid, where turbidimetry measurements (TM) has recently been successfully applied in different emulsion systems ( [5][6][7]).The turbidity spectrum of a MWF emulsion is influenced by the chemical composition, droplet size (distribution), internal structure, and concentration of the dispersed phase ( [5,8]).Measurement of the turbidity of an emulsion system requires contacting the fluid flow and optical access to the fluid.Therefore, during machine operation, sensor probes need to be dipped into the stream of the flowing MWF emulsion.The optical path through the liquid necessary for the measurement principle may be realized by an internal reflection arrangement.Thereby, the emulsion flows through a small channel within the sensor probe where a focused light beam passes the stream twice, forward and backward (see Figure 1).
Typical flow velocities of an MWF inside pumping lines of metal working machines are in the range 1 <  < 2.5 m/s [9].Segregation effects for the different emulsion droplet sizes for flow around the sensor probe may impact the measurement performance and hence need to be evaluated.Additionally, due to microorganisms in the water-based emulsion, biological films inside the probe that will interfere with the measurement may form.To avoid adhesion of microorganisms to sensor surfaces, it has been shown that minimum wall shear stresses τwall > 3-10 Pa are necessary ( [10,11]), where wall shear stress levels τwall > 20-100 Pa are found to even completely remove existing bio-films from surfaces ( [12]).Thus, in the present investigation, the multiphase emulsion flow of a working fluid around a TM sensor probe will be analyzed for the derivation of flow regimes of sustainable process conditions.

Multiphase Flow Analysis
Multiphase computational fluid dynamic (M-CFD) simulations are utilized to evaluate the flow patterns within and around typical sensor probes used in turbidity measurements with a rectangular slit channel as shown in Figure 1 [7].Multiphase flow around cylindrical objects has already been intensively investigated ( [13,14]).Liquid droplets tend to segregate in regions with low fluid velocities, leading to an increased local droplet concentration which may cause a decrease of the physical stability of the droplets [15].
For low Reynolds numbers (Re = 5-40), recirculating vortices may arise downstream of a cylinder.As the Reynolds number increases (Re > 40-350), these eddies periodically detach from the rear surface of the cylinder, carrying the droplets away from the probe [16].The authors of [17] have shown that, for low Stokes numbers Typical flow velocities of an MWF inside pumping lines of metal working machines are in the range 1 < u < 2.5 m/s [9].Segregation effects for the different emulsion droplet sizes for flow around the sensor probe may impact the measurement performance and hence need to be evaluated.Additionally, due to microorganisms in the water-based emulsion, biological films inside the probe that will interfere with the measurement may form.To avoid adhesion of microorganisms to sensor surfaces, it has been shown that minimum wall shear stresses τ wall > 3-10 Pa are necessary ( [10,11]), where wall shear stress levels τ wall > 20-100 Pa are found to even completely remove existing bio-films from surfaces ( [12]).Thus, in the present investigation, the multiphase emulsion flow of a working fluid around a TM sensor probe will be analyzed for the derivation of flow regimes of sustainable process conditions.

Multiphase Flow Analysis
Multiphase computational fluid dynamic (M-CFD) simulations are utilized to evaluate the flow patterns within and around typical sensor probes used in turbidity measurements with a rectangular slit channel as shown in Figure 1 [7].Multiphase flow around cylindrical objects has already been intensively investigated ( [13,14]).Liquid droplets tend to segregate in regions with low fluid velocities, leading to an increased local droplet concentration which may cause a decrease of the physical stability of the droplets [15].
For low Reynolds numbers (Re = 5-40), recirculating vortices may arise downstream of a cylinder.As the Reynolds number increases (Re > 40-350), these eddies periodically detach from the rear surface of the cylinder, carrying the droplets away from the probe [16].The authors of [17] have shown that, for low Stokes numbers the droplets will recirculate behind the cylinder and then entrain to the wall.At higher Stokes numbers, this effect vanishes.In this paper, simulations were carried out assuming various mean-sized monodispersed droplets (no distribution).Typical droplet sizes of metalworking fluid emulsions for machine operation are in the order of 100 nm for fresh mixed emulsions up to several microns for aged fluids [18].This leads to Stokes numbers in the order of 5 × 10 −9 for minimal droplet size and inlet velocity to 5 × 10 −3 in the case of maximal droplet size and velocity, thus at minimum three orders of magnitude less than unity.

Measurement Facilities
The numerical flow simulation results are evaluated by experiments in a test flow channel and a typical metal working machine.The fluid velocity in the flow channel varied between 0.04 m/s and 2.5 m/s and in the metal working machine between 0.3 m/s and 3.5 m/s; therefore, the range of Reynolds numbers investigated is Re cyl = u inlet d cyl ν = 260-22,750.The emulsion concentration was varied in the range of 1-10 v/v% in the flow channel, whereas the concentration in the machine was constant at 5 v/v%.The MWF under investigation is a commercial emulsion system with a narrow mono-modal droplet size distribution and a mean droplet diameter of d mean = 150 nm.

Numerical Model of Multiphase Flow
The Eulerian-Eulerian approach has been utilized for the numerical investigation of the multiphase flow around the turbidity sensor, where continuous and dispersed phases are both treated by solving separate momentum, mass, and continuity equations.The governing conservation equations are numerically closed with the Reynolds averaged Navier-Stokes method (RANS).The non-conservative form of the momentum equation in an incompressible two-component multiphase flow [19] is used as where α is the phase fraction, u is the velocity, p is the pressure, g is the gravity constant, Ik c,p is the combined Reynolds and viscous stress, and S includes all momentum source terms ( [20,21]).The combined turbulent stress ν t and viscous stress ν equal the effective viscosity ν eff c = ν + ν t .As shown by [22], the viscosity of emulsions can change significantly, especially for nanodroplets smaller than ca.200 nm and high volumetric fractions.In the present case, however, since the simulations were carried out for droplets with mean sizes of 0.1-10.0µm, and a volumetric fraction of 0.05, which results in a relatively small interaction area between dispersed and continuous phases, the emulsion viscosity was assumed to be constant.This assumption is consistent with the observed experimental values, and predictions by emulsion viscosity models [22].Furthermore, the Eulerian-Eulerian model approach is a two-fluid model and is based on the individual properties of each fluid.Therefore, the viscosity of the emulsion is not pertinent to the model.The ratio of viscosities is not pertinent, either.As the turbulent viscosity also depends on the phase fraction, the equation can be solved using an empirical correlation, e.g., by using a turbulent dispersion response coefficient C t .This approach for a turbulent flow description is governed by the fact that the discrete phase turbulent momentum fluxes are directly linked to their continuous phase counterparts [20].This leads to an equation for the dispersed phase: where the indices p and c correspond to the particulate dispersed and continuous phase, respectively.The value for the response coefficient is usually set to a constant value of 1, which is adequate for phase fractions above approximately 5% according to [23].Note that this way particle turbulent dispersion is modeled while turbophoresis is neglected.The continuity equation is applied to fulfill the assumption of an incompressible flow.Three main forces act on the dispersed droplets and define the source term for the momentum transfer where V is the cell volume.The drag F d , lift F l and virtual mass force F vm are derived from empirical correlations and theoretical models [24].The drag force depends on the relative velocity u rel = u p − u c and the shape of the particle, given by its area towards the flow A and its drag coefficient C D .Multiple correlations identify the value of the coefficient C D as a function of the Reynolds number.
The standard model for rigid spherical particles is an expanded version of the derivation by [25] to ensure the correct limiting behavior in the inertial regime and valid for small emulsion particles with d p < 10 µm in liquid/liquid flows as Re p < 30, Eo < O(10 −3 ) and µ p /µ c >> 1 [26].In addition, the emulsifier builds a layer around the droplets, thus preventing internal circulation and, therefore, changes in C D .The swarm influence may be neglected for volume concentrations below 10 v/v%.M-CFD simulations were carried out with the OpenFOAM [27] 2.1.0solver "twoPhaseEulerFoam" using a hexahedron structured mesh as shown in Figure 2. The region around the cylindrical probe is refined by generating an o-grid type mesh structure that allows perpendicular cell adjustment towards the probe walls (compare Figure 2, left, top).Hence, numerical diffusion can be minimized, and convergence and mesh performance are improved.The first grid point wall distance is a critical factor due to the small dimension of the measurement slit in the probe and varied depending on the flow speed.Low distance correlation [28] was utilized, where the y + wall distance is kept between 1 and 4 for the k-ω-model.The grid structure has been checked for proper resolution and grid size dependency.The mesh is validated regarding to the dispersed mass flow through the slit and the pressure drop.A mesh with about 850 thousand hexahedron cells finally is used that offers reasonable computational times while deviating less than 10% from the results obtained when using 1.8 million cells.
Fluids 2017, 2, 9 4 of 10 adequate for phase fractions above approximately 5% according to [23].Note that this way particle turbulent dispersion is modeled while turbophoresis is neglected.The continuity equation is applied to fulfill the assumption of an incompressible flow.Three main forces act on the dispersed droplets and define the source term for the momentum transfer where V is the cell volume.The drag   , lift   and virtual mass force  vm are derived from empirical correlations and theoretical models [24].The drag force depends on the relative velocity  rel =   −   and the shape of the particle, given by its area towards the flow A and its drag coefficient   .Multiple correlations identify the value of the coefficient   as a function of the Reynolds number.
The standard model for rigid spherical particles is an expanded version of the derivation by [25] to ensure the correct limiting behavior in the inertial regime and valid for small emulsion particles with dp < 10 μm in liquid/liquid flows as Rep < 30, Eo < O(10 −3 ) and μp/μc >> 1 [26].In addition, the emulsifier builds a layer around the droplets, thus preventing internal circulation and, therefore, changes in CD.The swarm influence may be neglected for volume concentrations below 10 v/v%.M-CFD simulations were carried out with the OpenFOAM [27] 2.1.0solver "twoPhaseEulerFoam" using a hexahedron structured mesh as shown in Figure 2. The region around the cylindrical probe is refined by generating an o-grid type mesh structure that allows perpendicular cell adjustment towards the probe walls (compare Figure 2, left, top).Hence, numerical diffusion can be minimized, and convergence and mesh performance are improved.The first grid point wall distance is a critical factor due to the small dimension of the measurement slit in the probe and varied depending on the flow speed.Low distance correlation [28] was utilized, where the y + wall distance is kept between 1 and 4 for the k-ω-model.The grid structure has been checked for proper resolution and grid size dependency.The mesh is validated regarding to the dispersed mass flow through the slit and the pressure drop.A mesh with about 850 thousand hexahedron cells finally is used that offers reasonable computational times while deviating less than 10% from the results obtained when using 1.8 million cells.The physical properties of the dispersed phase (oil) were taken from [29] (a) and (b) and the continuous phase was water.The flow inlet into the channel is characterized by a plug flow with a uniform velocity, a constant particle concentration, and a zero pressure gradient, whereas velocity and concentration at the outlet are represented by zero gradient boundary conditions.

Results and Discussion
Figure 3 illustrates the time averaged mean flow field approaching the probe and through the slit channel in the probe by streamlines.The splitting behavior of the flow is obvious.The flow around the cylinder dominates the flow structure and subsequently that behind the small slit.The emulsion is partly sucked into the slit and around the probe.
Fluids 2017, 2, 9 5 of 10 uniform velocity, a constant particle concentration, and a zero pressure gradient, whereas velocity and concentration at the outlet are represented by zero gradient boundary conditions.

Results and Discussion
Figure 3 illustrates the time averaged mean flow field approaching the probe and through the slit channel in the probe by streamlines.The splitting behavior of the flow is obvious.The flow around the cylinder dominates the flow structure and subsequently that behind the small slit.The emulsion is partly sucked into the slit and around the probe.Thus, high gradients occur at the slit inlet and the velocity pattern is close to a block profile.This pattern cannot develop over the probe slit length of 4.5 mm, resulting in a high velocity gradient throughout the entire slit.Due to vortex shedding in the flow field behind the cylinder, the mean droplet volume fraction inside the sensor slit obtains a time-dependent behavior with a variation of about ∆α ≈ 1%.
Figure 4 illustrates a certain droplet concentration fluctuation in the probe channel for a free stream velocity of 0.03 m/s ( cyl ≈ 195).The droplet concentration in the initial phase  < 5 s converges and develops into a time-dependent state for  > 5 s to be described by a sinusoidal function.The cylinder of the probe defines the droplet movement because of the detaching eddies.This temporal effect correlates with the Strouhal number  =  •   inlet (7) where the frequency f of drop concentration fluctuations is derived from the simulation with the characteristic length L as the diameter of the cylinder.The calculated Strouhal number  calc = 0.227 for the concentration fluctuation is related to the Strouhal number of a cylinder  cyl = 0.2 for the onset of the vortex shedding process at a Reynolds number  cyl ≈ 140-260 [17].Thus, high gradients occur at the slit inlet and the velocity pattern is close to a block profile.This pattern cannot develop the probe slit length of 4.5 mm, resulting in a high velocity gradient throughout the entire slit.Due to vortex shedding in the flow field behind the cylinder, the mean droplet volume fraction inside the sensor slit obtains a time-dependent behavior with a variation of about ∆α ≈ 1%.
Figure 4 illustrates a certain droplet concentration fluctuation in the probe channel for a free stream velocity of 0.03 m/s (Re cyl ≈ 195).The droplet concentration in the initial phase t < 5 s converges and develops into a time-dependent state for t > 5 s to be described by a sinusoidal function.The cylinder of the probe defines the droplet movement because of the detaching eddies.This temporal effect correlates with the Strouhal number where the frequency f of drop concentration fluctuations is derived from the simulation with the characteristic length L as the diameter of the cylinder.The calculated Strouhal number Str calc = 0.227 for the concentration fluctuation is related to the Strouhal number of a cylinder Str cyl = 0.2 for the onset of the vortex shedding process at a Reynolds number Re cyl ≈ 140-260 [17].
Fluids 2017, 2, 9 5 of 10 uniform velocity, a constant particle concentration, and a zero pressure gradient, whereas velocity and concentration at the outlet are represented by zero gradient boundary conditions.

Results and Discussion
Figure 3 illustrates the time averaged mean flow field approaching the probe and through the slit channel in the probe by streamlines.The splitting behavior of the flow is obvious.The flow around the cylinder dominates the flow structure and subsequently that behind the small slit.The emulsion is partly sucked into the slit and around the probe.Thus, high gradients occur at the slit inlet and the velocity pattern is close to a block profile.This pattern cannot develop over the probe slit length of 4.5 mm, resulting in a high velocity gradient throughout the entire slit.Due to vortex shedding in the flow field behind the cylinder, the mean droplet volume fraction inside the sensor slit obtains a time-dependent behavior with a variation of about ∆α ≈ 1%.The normalized time and volume averaged droplet concentration in the probe slit for various mean droplet sizes in the range 0.1-10 µm are illustrated in Figure 5.For each mean droplet size, a separate simulation run has been performed.An increase of the mean droplet concentration with increasing droplet size is found with a maximum of the concentration variations for d p ≈ 7 µm.The normalized time and volume averaged droplet concentration in the probe slit for various mean droplet sizes in the range 0.1-10 µ m are illustrated in Figure 5.For each mean droplet size, a separate simulation run has been performed.An increase of the mean droplet concentration with increasing droplet size is found with a maximum of the concentration variations for   ≈ 7 μm.During continuous MWF operation not only bacteria, but also fungi and yeast may appear that affects the machine operation and performance as well as the measurement reliability and accuracy [30].Even though in MWF operation the settling of monocultures of microorganisms is reported [12], their influence on each other is scarcely evaluated.Hence, microbes can act differently under certain circumstances while travelling through the sensor probe slit.Various types of microbes yield critical shear stress values to overcome adhesion ( [10,11,31]).Thus, the criteria of minimum wall shear stress can predict contamination of surfaces.
The calculated wall shear stress distribution in the probe slit channel of the sensor at a fluid velocity of 0.03 m/s is illustrated in Figure 6 as color plot on the slit surfaces.The cylindrical probe configuration can be seen from these results.The flow is separated to a flow around and through the slit.The fluid velocity can be kept high at the inner slit walls, where at the edges of the sensor the velocity gradients are very low.Even at high inlet velocities, the wall shear stress only increases slightly in these areas.Round edges might improve this part of the sensor probe to avoid bacterial contamination and strong differences in the maximum and minimum stresses.During continuous MWF operation not only bacteria, but also fungi and yeast may appear that affects the machine operation and performance as well as the measurement reliability and accuracy [30].Even though in MWF operation the settling of monocultures of microorganisms is reported [12], their influence on each other is scarcely evaluated.Hence, microbes can act differently under certain circumstances while travelling through the sensor probe slit.Various types of microbes yield critical shear stress values to overcome adhesion ( [10,11,31]).Thus, the criteria of minimum wall shear stress can predict contamination of surfaces.
The calculated wall shear stress distribution in the probe slit channel of the sensor at a fluid velocity of 0.03 m/s is illustrated in Figure 6 as color plot on the slit surfaces.The cylindrical probe configuration can be seen from these results.The flow is separated to a flow around and through the slit.The fluid velocity can be kept high at the inner slit walls, where at the edges of the sensor the velocity gradients are very low.Even at high inlet velocities, the wall shear stress only increases slightly in these areas.Round edges might improve this part of the sensor probe to avoid bacterial contamination and strong differences in the maximum and minimum stresses.The normalized time and volume averaged droplet concentration in the probe slit for various mean droplet sizes in the range 0.1-10 µ m are illustrated in Figure 5.For each mean droplet size, a separate simulation run has been performed.An increase of the mean droplet concentration with increasing droplet size is found with a maximum of the concentration variations for   ≈ 7 μm.During continuous MWF operation not only bacteria, but also fungi and yeast may appear that affects the machine operation and performance as well as the measurement reliability and accuracy [30].Even though in MWF operation the settling of monocultures of microorganisms is reported [12], their influence on each other is scarcely evaluated.Hence, microbes can act differently under certain circumstances while travelling through the sensor probe slit.Various types of microbes yield critical shear stress values to overcome adhesion ( [10,11,31]).Thus, the criteria of minimum wall shear stress can predict contamination of surfaces.
The calculated wall shear stress distribution in the probe slit channel of the sensor at a fluid velocity of 0.03 m/s is illustrated in Figure 6 as color plot on the slit surfaces.The cylindrical probe configuration can be seen from these results.The flow is separated to a flow around and through the slit.The fluid velocity can be kept high at the inner slit walls, where at the edges of the sensor the velocity gradients are very low.Even at high inlet velocities, the wall shear stress only increases slightly in these areas.Round edges might improve this part of the sensor probe to avoid bacterial contamination and strong differences in the maximum and minimum stresses.Figure 7 illustrates the calculated occurring minimum and maximum wall shear stress τ wall inside the slit of the sensor, which increases exponentially with increasing inlet velocities u inlet .The minimal wall shear stress τ wall,min is typically located next to the edges of the slit, and the maximal wall shear stress τ wall,max is located at the verge of the slit.Thus, microbiological adhesion to the walls will more likely take place at the middle and edges of the slit channel.Even at high inlet velocities, the wall shear stress stays low in these areas.Therefore, short slits are superior to longer ones.Flow velocities in the range of 1-3 m/s of MWF emulsions in metal working machines will lead to wall shear stress of τ wall = 10-200 Pa (see Figure 7).For the examined sensor geometry, a correlation for the minimum and maximum wall shear stress dependent on the inlet velocity is ln(τ max ) = 1.152• ln(u inlet ) + 4.798 and ln(τ min ) = 1.402• ln(u inlet ) + 2.504. ( Fluids 2017, 2, 9 7 of 10 wall shear stress τ wall,max is located at the verge of the slit.Thus, microbiological adhesion to the walls will more likely take place at the middle and edges of the slit channel.Even at high inlet velocities, the wall shear stress stays low in these areas.Therefore, short slits are superior to longer ones.Flow velocities in the range of 1-3 m/s of MWF emulsions in metal working machines will lead to wall shear stress of τ wall = 10-200 Pa (see Figure 7).For the examined sensor geometry, a correlation for the minimum and maximum wall shear stress dependent on the inlet velocity is An increase of the fluid velocity directly leads to higher velocity gradients and therefore wall shear stresses.The calculated minimum wall shear stress extending the literature values of 3 Pa for avoiding, respectively, the resolution of biomolecules from surfaces is achieved at mean inlet velocities  crit > 0.4 m/s.This threshold should be exceeded to avoid the biological impact on the turbidity measurements.
In order to study the impact of the MWF multiphase flow conditions around the sensor on the measurements, MWF turbidity spectra were detected in a test channel as well as in a machine during operation.The influence of the flow velocity and the droplet concentration on the detected spectra were investigated.Figure 8 illustrates the standard deviation σ between the turbidity spectra measured at a MWF test channel ( =0 ) and the in-machine operation turbidity spectra measured τ( = inlet ) for a concentration range of 1-10 /%: where m is the number of turbidity measurements in the spectra τ(λ).The standard deviation σ for all concentrations and velocities is about 10 −1 for the test channel and less than 10 −1 for the measurement in the machine.The lower standard deviation for the machine is caused by the somewhat aged MWF compared to the fresh MWF used in the flow channel.Figure 8 illustrates the spectra for MWF emulsion with 5 /%.The spectra do not change with flow velocity, indicating that there is no significant change of the concentration inside the slit for velocities up to 3.5 m/s in the metal working machine.The experimental data validate the CFD simulations, as no segregation effects have been detected in the velocity range typically occurring in metal working machines.An increase of the fluid velocity directly leads to higher velocity gradients and therefore wall shear stresses.The calculated minimum wall shear stress extending the literature values of 3 Pa for avoiding, respectively, the resolution of biomolecules from surfaces is achieved at mean inlet velocities u crit > 0.4 m/s.This threshold should be exceeded to avoid the biological impact on the turbidity measurements.
In order to study the impact of the MWF multiphase flow conditions around the sensor on the measurements, MWF turbidity spectra were detected in a test channel as well as in a machine during operation.The influence of the flow velocity and the droplet concentration on the detected spectra were investigated.Figure 8 illustrates the standard deviation σ between the turbidity spectra measured at a MWF test channel τ(λ u=0 ) and the in-machine operation turbidity spectra measured τ λ u=u inlet for a concentration range of 1-10 v/v%: where m is the number of turbidity measurements in the spectra τ(λ).The standard deviation σ for all concentrations and velocities is about 10 −1 for the test channel and less than 10 −1 for the measurement in the machine.The lower standard deviation for the machine is caused by the somewhat aged MWF compared to the fresh MWF used in the flow channel.Figure 8 illustrates the spectra for MWF emulsion with 5 v/v%.The spectra do not change with flow velocity, indicating that there is no significant change of the concentration inside the slit for velocities up to 3.5 m/s in the metal working machine.The experimental data validate the CFD simulations, as no segregation effects have been detected in the velocity range typically occurring in metal working machines.

Figure 1 .
Figure 1.Measurement arrangement of inline light transmission sensor probe with cylindrical shaft in a flow channel of machining process.

Figure 1 .
Figure 1.Measurement arrangement of inline light transmission sensor probe with cylindrical shaft in a flow channel of machining process.

Figure 2 .
Figure 2. Setup of the sensor probe arrangement and display of the main features of the mesh.The physical properties of the dispersed phase (oil) were taken from [29] (a) and (b) and the continuous phase was water.The flow inlet into the channel is characterized by a plug flow with a

Figure 2 .
Figure 2. Setup of the sensor probe arrangement and display of the main features of the mesh.

Figure 4 .
Figure 4. Time-dependent droplet concentration inside the sensor probe slit for a droplet size of 150 nm.

Figure 4
illustrates a certain droplet concentration fluctuation in the probe channel for a free stream velocity of 0.03 m/s ( cyl ≈ 195).The droplet concentration in the initial phase  < 5 s converges and develops into a time-dependent state for  > 5 s to be described by a sinusoidal function.The cylinder of the probe defines the droplet movement because of the detaching eddies.This temporal effect correlates with the Strouhal number =  •   inlet(7)where the frequency f of drop concentration fluctuations is derived from the simulation with the characteristic length L as the diameter of the cylinder.The calculated Strouhal number  calc = 0.227 for the concentration fluctuation is related to the Strouhal number of a cylinder  cyl = 0.2 for the onset of the vortex shedding process at a Reynolds number  cyl ≈ 140-260[17].

Figure 4 .
Figure 4. Time-dependent droplet concentration inside the sensor probe slit for a droplet size of 150 nm.

Figure 4 .
Figure 4. Time-dependent droplet concentration inside the sensor probe slit for a droplet size of 150 nm.

Figure 5 .
Figure 5. Dependency of the droplet concentration in the sensor probe slit on the droplet size.

Figure 6 .
Figure 6.Mean wall shear stress distribution inside the probe slit for an inlet velocity of  inlet = 0.03 m/s and a drop size of 150 nm.

Figure 7
Figure 7 illustrates the calculated occurring minimum and maximum wall shear stress τ wall inside the slit of the sensor, which increases exponentially with increasing inlet velocities  inlet .The minimal wall shear stress τ wall,min is typically located next to the edges of the slit, and the maximal

Figure 5 .
Figure 5. Dependency of the droplet concentration in the sensor probe slit on the droplet size.

Figure 5 .
Figure 5. Dependency of the droplet concentration in the sensor probe slit on the droplet size.

Figure 6 .
Figure 6.Mean wall shear stress distribution inside the probe slit for an inlet velocity of  inlet = 0.03 m/s and a drop size of 150 nm.

Figure 7
Figure 7 illustrates the calculated occurring minimum and maximum wall shear stress τ wall inside the slit of the sensor, which increases exponentially with increasing inlet velocities  inlet .The minimal wall shear stress τ wall,min is typically located next to the edges of the slit, and the maximal

Figure 6 .
Figure 6.Mean wall shear stress distribution inside the probe slit for an inlet velocity of u inlet = 0.03 m/s and a drop size of 150 nm.

8 )Figure 7 .
Figure 7. Maximal and minimal wall shear stress within the sensor probe in dependence of the inlet velocity  inlet and a drop size of 150 nm.

Figure 7 .
Figure 7. Maximal and minimal wall shear stress within the sensor probe in dependence of the inlet velocity u inlet and a drop size of 150 nm.