Issue 
Metall. Res. Technol.
Volume 119, Number 1, 2022



Article Number  109  
Number of page(s)  16  
DOI  https://doi.org/10.1051/metal/2021088  
Published online  20 January 2022 
Regular Article
Experimental and numerical modelling of mass transfer in a refining ladle
^{1}
Sorbonne Université and CNRS, UMR 7190, Institut Jean Le Rond d’Alembert, BP 162, 4, Place Jussieu, 75252 Paris Cedex 05, France
^{2}
ArcelorMittal Maizières Research, BP 30320, Voie Romaine, 57283 MaizièreslèsMetz Cedex, France
^{3}
Institut Universitaire de France, Paris, France
^{*} email: nelson.joubert@sorbonneuniversite.fr
Received:
2
July
2021
Accepted:
8
November
2021
Mass transfer between liquid steel and slag is an important physical phenomenon during secondary metallurgy for prediction of the chemical reaction rate and adjustment of liquid steel composition. We study this phenomenon at ambient temperature with a water experiment and perform Direct Numerical Simulations, aiming to reproduce an argongas bottomblown ladle. First, we measure the evolution of the timeaveraged openeye area as a function of the air flow rate. Both simulation and experiment agree relatively well and are close to other water experiments in the literature. Secondly, the mass transfer of thymol between water and oil is investigated. The experimental results show that two mass transfer regimes can be observed. The regime change coincides with atomization of the oil layer resulting in the continuous formation of oil droplets in the water whenever the air flow rate rises above a critical value. The numerical results for the masstransfer rate or Sherwood number are obtained at small Schmidt numbers and are then extrapolated to the experimental Schmidt number of 1480. A good agreement with experiment is observed although with large error bars. The Sherwood numbers at the two largest simulated flow rates show a steep increase.
Key words: mass transfer / computational fluid dynamics / volumeoffluid / large Schmidt numbers / multiphase flow
© N. Joubert et al., Published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
In numerous natural and industrial processes, chemical reaction is strongly coupled with a fluid flow. In the steel industry, during the secondary metallurgy step, the adjustment of liquid steel composition at high temperature in a ladle involves several chemical reactions between steel and slag [1]. Chemical reactions can be broken down into three stages: transport by convection of reactive species within the liquid phase toward the interface, transport by molecular diffusion through the concentration boundary layer and finally chemical reaction at the interface. The global kinetics is governed by the kinetics of the slowest stage. In the case of a chemical reaction between liquid steel and slag, it is generally accepted that chemical reactions at the interface are very fast. Then, the steps governing the global kinetics are the convective and diffusive transport or mass transfer of reactive species. In order to better adjust the steel composition, it is important to characterize the mass transfer of reactive species between the liquid steel and slag phase. Here, for example, the process to be studied is the desulfurization of liquid steel by chemical reaction with slag under argon gas stirring. Previous studies attempted to model experimentally the process either with ambient temperature ladle [2–4] or hightemperature ladle [3,5,6]. Other studies try to model the process numerically [7–9]. The results showed some change in the mass transfer regime between the two phases when the gas flow rate reaches a particular value. In the literature, it is generally acknowledged that the observed mass transfer regime change is due to the high deformation of the slag layer leading ultimately to its fragmentation into slag droplets in the liquid steel [10–13]. Slag droplets formation would increase the exchange area between the two phases. In order to verify this assumption and the global behavior of the mass transfer, we developed simultaneously an experimental water ladle of the process and its numerical model. Considering constant temperature and fluids properties in the ladle, the gas flow rate is the control parameter of the process. Then we will vary the gas flow rate to measure its influence on the mass transfer. In this study, we start by giving details on the experimental and numerical procedure. Then the result section is divided in two parts. The first one is a brief hydrodynamic characterization of the water experiment and, the second one concern the mass transfer characterization.
2 Model description
The main characteristics of the flow consist of the interaction of a gas bubble plume with a slag layer at the free surface. Given the numbers of parameters of the physical model, many dimensionless number quantities can be established. As a result, the choice of similitude to follow is not trivial and depends on the phenomena of interest. In this study, we are interested in the characterization of the phenomenon of mass transfer between liquid steel and slag. Nevertheless, this mass transfer depends strongly on the fluid flow characteristics.
2.1 Experimental model
If we assume that within the gas flow rate range considered in the industrial case, we are in presence of a bubble plume [14,15], the vertical rising velocity due to the bubble plume in the liquid steel can be approximated by a characteristic velocity scale given by(1)
Now, considering the expression of (1) we can rewrite the Froude number(2)
as a function of the gas flow rate to obtain the expression of the Froude number N(3)
To determine the corresponding gas flow rate between a 180ton industrial ladle and our experiment a Froude number similarity can be used in which the model experiment and the industrial ladle have the same Froude number.(4)
with h_{m} and h_{s} respectively the height of the steel and the slag equivalent phase. il and rl are respectively the indices for the industrial and reduced ladle and Q is the gas flow rate considered at the operating temperature and pressure at the injector height. In this study, we use a cubic ladle filled of water with an oil layer at the top. Both oil and water are at ambient temperature. Figure 1 displays a sketch of the experimental setup of the water ladle and the main geometric parameters of the ladle can be found in Table 1. It is a squaresection transparent ladle made with acrylic glass, with a single bottom centered circular air injection hole. Air flow is injected from a compressed air network and can be varied through an automate controlling the valve opening. In the following, we denote the water, oil, and air phase with the respective indices w, o and a.
As in the industrial process, it is a threephase experiment where the air phase represents the argon gas, water represents the liquid steel, and the oil mixture represents the liquid slag. The fluids choice follows the one used in [2], the oil phase is a 50–50% in volume mixture of cottonseed oil and paraffin oil. All the fluids physical parameters can be found in Table 2. Following [2], to reproduce the behavior of the Sulphur dissolved in liquid steel we have chosen thymol (C_{10}H_{14}O) dissolved in water as chemical tracer in our experiment. In the experiment thymol concentration in water is measured every 10 minutes during 3 hours by analyzing water samples with a refractometer previously calibrated. Measurements can be done with two size of injection diameters d_{inj} = 2.35 mm and d_{inj} = 7.9 mm, with 14.4 l of water with a thymol concentration of 0.9 g/l and 0.486 1 of oil mixture atop of it.
We detail here the general solution of the evolution of the concentration of thymol in the water. The mass transfer equation of thymol in each phase gives(5) (6)where C_{w} (g/l) and C_{o} (g/l) are respectively the concentration of thymol in water and oil bulk, and are respectively the concentration of thymol at the interface on water and oil side, K_{w} (m/s) and K_{o} (m/s) are respectively the global mass transfer coefficient of thymol in water and oil, V_{w} (m^{3}) and V_{o} (m^{3}) are respectively the water and oil volume and J (g/m^{3}s) is the diffusive flux. The partition coefficient of thymol comparing the solubilities of thymol in water and oil at equilibrium is given by(7)
Using mass balance for thymol between water and oil from (5) and (6) we can write(8)
After time integration we can write(9)where C_{w,o} (g/l) is the initial concentration of thymol in the water phase bulk. Using (8) and (9) we can rewrite (5) and (6) (10) (11)with τ_{w} = V_{w}/K_{w}A and τ_{o} = V_{o}/K_{o}A, combining (10) and (11) we can obtain the following differential equation(12)with α = K_{w}/PK_{o} a casespecific constant assuming K_{w} and K_{o} are time independent. β = V_{w}/V_{o}P is a constant depending only on the fluid quantity and nature. Considering that at t = 0 we know that C_{o,0} = 0 g/l, we obtain the following solution for the evolution of concentration of thymol in the oil phase(13)with . Using (9) and (13) we can obtain the evolution of concentration of thymol in the water phase(14)
Following (14), we can see that the concentration of thymol follows an exponential decay from the initial value C_{w,0} to an equilibrium value of C_{w,0} − C_{w,0}/(1 + β). The characteristic concentration decay time is given by α and β need to be small to have a minimum concentration decay time around τ_{c,min} ∼ τ_{w}. It is generally acknowledged that in an industrial ladle the desulphurization reaction of liquid steel with slag has a high PK_{o} value. Then, the steel phase mass transfer is controlled by the steel phase mass transfer resistance. In the case of water experiment, thymol has a high partition coefficient between oil and water (P > 350). Assuming that K_{o} is not too small it leads to a mass transfer controlled by the water phase resistance.
Now, if we take the logarithm of (14) and use the expression of λ we can write(15)
In reference [2], the authors assume that α ≪ 1/β is verified in the mass transfer configuration associated with the desulfurization process. Considering the same assumption in our study, we can rewrite (15) to obtain(16)
In the expression above, the product K_{w}A/V_{w} is called the capacity coefficient. From relation (16), it is possible to estimate directly the capacity coefficient from the measured thymol concentration in water on the righthand side of (16). This formulation is used to determine the global mass transfer coefficient in the following experimental results.
Fig. 1 Schematic of the experimental setup of the water ladle. 
Main geometrical parameters of the experiment.
Physical properties of fluids used in our water experiment.
2.2 Numerical model
To resolve the problem, we solve numerically the partial differential equations with the opensource Free Software library Basilisk. Basilisk allows to perform Direct Numerical Simulation or DNS, meaning there are no turbulence models as the NavierStokes equations are resolved down to the smallest scales. Basilisk has the possibility to adapt the mesh dynamically following a quadtree /octree structure. Typically, the minimum mesh size Δ_{min} can be computed from the maximum level of refinement maxlevel which defines the maximum number of cells in one direction and the size of the domain L_{0} with Δ_{min} = L_{0}/2^{maxlevel}. The numerical schemes used in Basilisk are based on its ancestor the Gerris solver and can be found in [16,17]. The flow is described with the incompressible NavierStokes equations with surface tension term treated by a Continuous Surface Force [18] which can be writtenwith D = (∇ u + (∇ u) ^{T})/2 the deformation tensor. Because we are dealing with multiphase flow the variations of density and viscosity inside the domain are described with the help of the Volume of Fluid method. The volume fraction χ (x, t) is described by(17)
which obeys to an advection equation(18)
Because we are in presence of a threephase flow, we cannot use the implicit declaration of phase two as in a twophase flow. Instead, we declare explicitly three volume fractions corresponding to each phase. To ensure that a cell is not filled with more than one phase the sum of n volume fractions should respect(19)
To do so at each timestep we normalize the volume fractions χ_{k} of the k − th phase with the sum of volume fractions(20)
This way we can follow a one fluid description with variable density and viscosity determined in the domain with arithmetic mean(21)
Surface tension is acting on the interface between two fluids, but here there is more than one possibility of fluid neighbor for one phase. That is why we decompose the physical surface tension into phasespecific surface tension just depending on the phase and not the fluid in contact following [19,20](22)
To take into account the concentration variation we need to consider also the generic scalar transport equation of concentration c (amount of chemical species/unit volume) for an incompressible flow and without sources or sinks(23)
We can separate (23) in two parts:
∇ ⋅ (D ∇ c) describes the molecular diffusion of c;
∇ ⋅ (uc) describes convection of c by the flow.
c has no influence on the velocity field determined by solving the NavierStokes equation and is only used to consider diffusion, therefore it is a passive scalar or tracer. We solve the two parts of (23) in two distinct steps and start by solving the convection part. To do so, we define a kthphasespecific concentration following(24)
T_{k} is confined to a VOFdetermined phase in the same manner as described in [21]. The method in the latter reference ensures that the tracer c cannot cross a VOF interface and limits the unwanted numerical diffusion. Therefore, instead of directly doing the convection of c we do it on T_{k}.
Before considering the diffusion of the tracer c computed from (24), we need to consider the concentration jump in the interfacial cells. Because of the formulations used in the Basilisk code, we could not easily add the equivalent of the partition ratio P at the interface. We use instead a Dirichlet condition imposed on the kth phase side(25)
This is an important assumption made to approximate the diffusion in the configuration of interest. Within this assumption it is not possible to ensure that the absolute diffusive flux of tracer leaving one phase stays constant when entering another phase.
Then we solve the second part of (23) which is the molecular diffusion of c. Once the tracer c has diffused, the value of T_{k} is updated with the new value of c using (24).
Finally, to simplify the model we consider constant phasespecific diffusion coefficients in the water and oil phase. In the air phase we set the diffusion coefficient to D_{a} = 0 m^{2}/s as we neglect its influence on the mass transfer.
To have an idea of the smallest scales to be resolved in the simulation we can use the Kolmogorov length scale given by(26)
Using an approximation of the kinetic energy dissipation in the water experiment configuration we can directly compute (26) from the gas flow rate value with(27)
At the lowest air flow rate considered in our experiment Q = 0.6 l/min it gives η = 165 μm.
In reference [22], the authors postulated a criterion to determine the minimum mesh size needed to resolve the Kolmogorov length scale in a DNS Δ_{min} ≤ 2.1η. As we are interested about mass transfer, we need to consider the Schmidt number Sc = μ/ρD characterizing the ratio of momentum diffusivity and mass diffusivity. In the experiment the Schmidt number of thymol in water is Sc_{th,w} = 1.48 × 10^{3}. The higher is the Schmidt number, the smaller is the concentration boundary layer. To avoid the use of an extremely small mesh size we choose to do simulation with smaller Schmidt number than Sc_{th,w} and consider four tracers in the water phase in the range of Sc_{w} ∈ {1, 4, 10, 40}. We also consider four tracers in the oil phase but at higher Schmidt number than in the water phase because of the difference in the viscosity and density. The range in the oil phase is Sc_{o} ∈ {10, 40, 100, 400}. Then we will extrapolate the numerical Sherwood number values to the experimental Schmidt number. At low value of Schmidt number, the concentration boundary layer thickness is approximately equal to the hydrodynamic boundary layer thickness. If we consider a minimum mesh size of Δ_{min} = 104.4 μm corresponding approximately to a maximum mesh resolution of 2^{11} cells in one direction. At this mesh resolution we can describe the initial oil layer with a maximum of 50 cells. Also, if we consider an oil droplet of 3 mm of diameter with the same mesh resolution we end up with 22 cells in its diameter.
3 Results
In the following, for easier comparison between experiments using other fluids than water and oil we will denote the steel and the slag equivalent phase respectively with the indices m and s. We also have to mention that we have used in the numerical simulation a reference value for the surface tension value of σ_{a/w} = 72 × 10^{−3} N/m instead of the measured value of Table 2.
3.1 Hydrodynamic characterization
3.1.1 Experimental results
In Table 3 we compute the viscosity ratio m = μ_{m}/μ_{s}, the density ratio r = ρ_{m}/ρ_{s} together with the bath height ratio h_{m}/L_{x} and the slag height ratio h_{s}/h_{m}. We can see from Table 3 that the viscosity ratio is close to the liquid steelslag value of the industrial measurement of [23] for almost all the experiments except the one of [24] using motor oil. Now, if we look at the density ratio, we can see that the water experiments present a lower density ratio compared to the liquid metal experiments. If we look now at the bath height ratio, we can see that in all the experiments, it is mostly around 0.75 except for the experiment of [24]. Finally, for the slag height ratio, we can see that in most of the experiments the value is not far from the industrial configuration corresponding to a thin slag layer.
Images of the top surface of the water experiment are recorded with a video camera making possible to follow the deformation of the oil layer under the influence of gas injection. In Figure 2 showing an instantaneous image of the top view of the water experiment, we can identify the contour of the open eye of the oil layer with a yellow line and the bubbles spout contour in orange. Unfortunately, it has not been possible to automatically plot the open eye contour of all the images for each air flow rate. So, we choose to rather plot manually the open eye contour for four frames separated by approximately 1.8s for each air flow rate. Then, these plots have been postprocessed using ImageJ to determine the experimental time average open eye area and its standard deviation.
In Figure 3 we plot the dimensionless open eye area as a function of the Froude number. From Figure 3, first we can note that most of the water experiments with cross symbols show similar evolution of the open eye area as a function of the Froude number. On the other hand, there is more dispersion between the liquid metal experiments. Indeed, we observe for the liquid metal experiments a smaller slope for the evolution of the dimensionless open eye area than for the water experiments. If we look at the industrial results of [23] (filled red symbols in Fig. 3), we can see that the latter are closer to the water experiments than the cold liquid metal experiments of [23,25]. This may be due to the very high value of r and in the case of [22] of h_{s}/h_{m} in the cold liquid metal experiments. As a result, it is likely that the implicit correlation in Figure 3 (28)
the exact form of the function f would still require some investigations. To conclude, the literature shows that the open eye area is very strongly affected by the fluid density ratio and perhaps the ratio of heights h_{s}/h_{m}. Extrapolation from the experimental water/oil configuration to the industrial liquidsteel/slag configuration should take this into account.
Fig. 2 Top view of the instantaneous experimental open eye obtained at Q = 0.6 l/min. The color code is in yellow: open eye contour; orange: spout zone. 
Fig. 3 Ratio of time average open eye area with the square of the bath height as a function of the Froude number for our experiment and the experiments of Krishnapisharody et al. [24], Yonezawa and Schwerdtfeger [23] and Thunman et al. [25]. 
3.1.2 Numerical results
To capture the open eye numerically first we do simulations of the water experiment configuration. Then we postprocess the generated images of the top view of the oil interface with the software ImageJ to get the open eye area. From Figure 4 we can make a clear distinction from the oil interface colored in red and the other phase in blue. Taking advantage of this we could obtain this time the temporal evolution of the open eye area from all the images of the numerical simulations.
In Figure 5, we compare the time evolution of the open eye area obtained for different air flow rates with d_{inj} = 7.9 mm. First, it is important to mention that the data gap observed for Q = 1.5 l/min is due to a format problem of the generated images and is not interpolated from the nonmissing data. We can see that globally the open eye area follows the increase of the air flow rate and that waiting ten seconds of physical time seems to be enough to reach a statistical steady state. If we look at the data for Q < 3.5 l/min we can see an important change in the slope during the transient regime and on the asymptotic value of the open eye area when the air flow rate is increased. Looking now at results for Q ≥ 3.5 l/min we observe mostly difference on the fluctuations between Q = 3.5 l/min and Q = 5.5 l/min. We can note that the convergence of the open eye area toward an asymptotic value is slower for high air flow rate than for low air flow rate.
We can see on Figure 6 that the numerical simulation reproduces qualitatively well the same flow behavior as the one that will be seen later on Figure 7. We can identify a central bubble plume generating an open eye in the oil layer when it reaches the free surface. The oil layer is relatively stable, and no oil droplet is observed at low air flow rate (cf. Fig. 6 top picture). While the free surface is highly perturbated at high air flow rate (cf. Fig. 6 bottom picture) and some oil droplets of various sizes can be identified in the water phase.
Fig. 4 Top view of the instantaneous numerical open eye, obtained at Q = 0.6 l/min. The numerical results are obtained with a minimum grid size of Δx = 527 μm and the oil interface is colored in red. 
Fig. 5 Time evolution of the open eye area for different air flow rates with d_{inj} = 7.9 mm. The results are obtained with a minimum grid size Δ_{min} = 527 μm corresponding to a maxlevel = 9 of adaptive refinement. 
Fig. 6 Front view of the simulation with d_{inj} = 7.9 mm with a minimum grid size Δ_{min} = 527 μm. The air interface is colored in green, and the oil interface is colored in red. 
3.1.3 Comparison of the results
In Figure 8 we plot the timeaveraged numerical and experimental open eye area. First, we can see that the increase of the open eye area with flow rate can be well described by a logarithmic function of the flow rate in both cases. Furthermore, we can note that the numerical results overestimate the experimental open eye area, and that this tendency is more pronounced at low air flow rate than at high air flow rate. We have searched unsuccessfully for an explanation of the discrepancy.
One possibility is that the spatial discretization error is still large. We have compared Δ_{min} = 527 μm and Δ_{min} = 264 μm simulations without noticing an improvement when the grid is refined but it is possible that very thin oil layers form because of the positivity of the spreading coefficient(30)
Since S = 0.0148 N/m > 0, oil spreads in thin layers [26]. These layers could be much thinner than even the smallest grid size used, so that even finer grids would be needed. This would explain the discrepancy between numerics and experiment.
Fig. 7 Front view picture of the water experiment with thymol dissolved in the water obtained with d_{inj}=7.9 mm. 
3.2 Mass transfer characterization
3.2.1 Experimental results
From the time evolution of the concentration of thymol in water (cf. Fig. 9, top picture), we can observe a good fit of the measurement with an exponential decay in agreement with correlation (14). Now, from the concentration measurement we can plot minus the righthand side of (16) in Figure 9 bottom picture. Taking the slope of a linear fit of the curve, we can compute the product of the time average mass transfer coefficient in the water with the interfacial area or capacity coefficient.
From the results of Figure 9 bottom picture, the slope of the linear fit is 6.2 × 10^{−5} so, at this air flow rate it gives a capacity coefficient of .
We repeat this process for different air flow rates and different injection diameters in Figure 10. From our experimental results in Figure 10 first we can observe that the capacity coefficient increases when the air flow rate is increased for both injection diameters. Even if we have less data for the d_{inj} = 2.35 mm case it appears that capacity coefficient varies only a little compared to the d_{inj} = 7.9 mm case. We can also identify two different mass transfer regimes below and above a critical air flow rate Q_{c} ∼ 5.0 l/min for both injection diameters. Secondly, if we compare our results with the one obtained by [2], we observe a similar evolution with a critical air flow rate qualitatively the same as the one we measured with d_{inj} = 7.9 mm. However, the slope of the first mass transfer regime is higher in the results of [2] than in our results. For the second mass transfer regime similar slope is observed between our results and the results of [2].
In reference [2] the authors justify the change of mass transfer regime observed in Figure 10 by the formation of oil droplets of various sizes at high air flow rate. In the top picture of Figure 7 (Q = 1.1 l/min) no oil droplet is visible in the water. At higher air flow rates, oil droplets appear transiently a few minutes after the start of air injection, then disappear probably by flotation to the top oil layer. If we increase the air flow rate further, a lot of dark spots which correspond to oil droplets can be identified (see the bottom picture of Figure 7 for Q = 7.5 l/min). Indeed, oil droplets of various sizes and shapes are present in water both at the start of air injection and after three hours of agitation at this high air flow rate. These oil droplets increase the interfacial area between water and oil resulting in an abrupt increase of mass transfer.
Now, to compare our results with results in the literature we compute the time average Sherwood number which is the ratio of convective mass transfer to the rate of diffusive mass transport. We consider as characteristic length scale L_{c} = L_{x}/4, (the 1/4 factor considers that boundary layers develop on the side of the open eye overs distances shorter than L_{x}) then we can define(31)
Figure 11 shows a comparison of the time average Sherwood number for water experiments (our results with d_{inj} = 7.9 mm and the one of [2]) and liquid metal experiments [3,6]. From the results of Figure 11, we observe for each experiments an abrupt increase of the mass transfer corresponding to the start of the second mass transfer regime. The results of [6] show a higher value of the Sherwood number than the other experiments with a difference of a factor ten with the results of [3]. This is interesting to note from Table 4 that these two experiments are done with similar Schmidt number in the steel equivalent phase. Now for the water experiments, if we look back at Figure 10 we see that there is a small difference in the value of the capacity coefficient between our results and the one of [2]. If we assume that these two experiments follow the same chemical reaction between water and oil, it means that the difference observed in Figure 11 is mainly due to the difference in the size of the ladle. The same goes for the difference in the Sherwood number for the liquid metal experiments of [3,6]. Beside the difference in the Sherwood number value, we can see that the liquid metal experiments also show a smaller slope for the second mass transfer regime than for the water experiments.
Fig. 8 Evolution of the time averaged open eye area of several experiments with its standard deviation and numerical timeaveraged openeye area with its temporal standard deviation as a function of the air flow rate obtained with d_{inj}=7.9 mm. The simulation results are obtained with a minimum grid size of D_{min}=527 mm corresponding to a maxlevel=9 of adaptive refinement. 
Fig. 9 Top: Time evolution of the concentration of thymol in water at Q=2.5 l/min and d_{inj}=7.9 mm. The error bars correspond to the measurement errors (cf. Supplementary Material: Fig. 1). Bottom: Time evolution of minus the righthand side of (16) obtained with Q=2.5 l/min and d_{inj}=7.9 mm. 
Fig. 10 Evolution of the capacity coefficient or product of the time average mass transfer coefficient with the interfacial area as a function of gas flow rate for our experiment and the one of Kim et al. [2]. 
3.2.2 Numerical results
We perform simulations for five distinct air flow rates below and around the mass transfer transition observed experimentally. The simulations are started after t = 7.6 s of physical time with a minimum grid size of Δ_{min} = 527 μm corresponding to a maxlevel = 9 of adaptive refinement. Then when the mass transfer transitory regime is finished, we use the obtained results to restart the simulation with half the minimum grid size corresponding to a maxlevel = n + 1 of adaptive refinement. This allows us to save computational time by not recomputing the mass transfer transitory regime for each mesh resolution. First, we consider a case with a constant air flow rate of Q = 0.6 l/min and various mesh resolutions. On Figure 12 we observe a slow decrease of the normalized concentration of tracer in water. The decrease over the eight seconds of simulation is small (2%) due to the limited simulation time. When the mesh resolution increases from a maxlevel = 9 to maxlevel = 10 of adaptive refinement we can see an important change in the slope of the normalized concentration. We remark that the change of the slope is smaller when increasing the maxlevel of adaptive refinement from maxlevel = 10 to maxlevel = 11. This indicates that at Sc = 10 and Q = 0.6 l/min we start to obtain mesh convergence from a maxlevel = 10 of adaptive refinement corresponding to a minimum mesh size of Δ_{min} = 263.5 μm.
If we now look at the results of Figure 13 obtained at a higher air flow rate of Q = 3.5 l/min, we observe globally the same behavior than in Figure 12. However, even if we have less data for this high air flow rate, in Figure 13 we can observe some small oscillations on the normalized concentration. It appears also that there is a significative difference of the slope of the normalized concentration when the mesh resolution increases from maxlevel = 10 to maxlevel = 11 of adaptive refinement.
From the time evolution of the concentration of tracer in water, we can compute the mass transfer coefficient in the water assuming that the concentration at the interface is set to zero because of the Dirichlet condition. So, we can write(32)
As we do not know the experimental oilwater interfacial area, we consider the wateroil area without gas injection as reference area . From this we can compute the instantaneous numerical Sherwood number using again definition(33)
In Figure 14 we can observe for a higher mesh resolution than a maxlevel = 9 a transitory regime where the numerical Sherwood number increases until it reaches an asymptotic value. We can also observe some small fluctuations around an asymptotic value which decrease while the mesh resolution is increased. If we compare the asymptotic value of the numerical Sherwood number (Sh_{w,num}) for the different mesh resolutions, we can note a higher value of Sh_{w,num} at a low mesh resolution than at a very fine mesh resolution.
Now, if we look at the time evolution of the numerical Sherwood number at Q = 3.5 l/min in Figure 15 we observe that we have not yet reached an asymptotic value of the Sherwood number for a maxlevel = 11 of adaptive refinement. We can also remark that in addition to some small fluctuations of the numerical Sherwood number some more important oscillations from an asymptotic value can be seen in Figure 15 considering a maxlevel < 11. This leads to a higher standard deviation on the fitted value of the asymptotic or time average numerical Sherwood number in Figure 15 than in Figure 14.
We repeat this process for several Schmidt numbers Sc_{w} ∈ {1, 4, 10, 40} and for all the air flow rates Q ∈ {0.6, 2.5, 4.5, 5.5, 6.5} l/min. (Some complementary data on the simulations and their results can be found in the Supplementary Material section and in the author’s PhD thesis [27]).
We also want to verify the assumption made before that the product of coefficients α and β defined above is small. First, we compute the time average numerical mass transfer coefficient from the timeaveraged value of the Sherwood number with(34)
Then the ratio of the time averaged mass transfer coefficients ratio between the water and the oil phase is(35)
To compute (35) we need to determine the value of the Sherwood number in the oil phase Sh_{o}. To do so we proceed in the same way than for the tracer in the water phase by first measuring the time evolution of the concentration of tracer in the oil phase (cf. Fig. 16). We can see from Figure 16 results that at Q = 0.6 l/min the normalized concentration of tracer in the oil phase decreases more rapidly than in the water phase (cf. Fig. 12). This can be explained by a higher diffusion coefficient in the oil phase due to the difference in the viscosity and density with the water phase. We can also note that almost no difference in the slope of the normalized concentration can be seen. Now, if we look at the Sherwood number in the oil phase in Figure 17 we can see that all the curves converge toward the same asymptotic value of for all the mesh resolutions. We can also note more fluctuations on the numerical Sherwood number in the oil phase that what could be observed in the water phase (cf. Fig. 14).
Using the time average numerical Sherwood number in the oil phase (cf. Sect. 8.3 and appendix 11.2 of [27]) and Table 4 we can compute (35) for each air flow rates. Then, we can plot the evolution of the ratio of the mass transfer coefficient as a function of the Froude number in Figure 18. We can see from Figure 18 that when the ratio of Schmidt number Sc_{r} = Sc_{w}/Sc_{o} decreases the ratio of timeaveraged mass transfer coefficients increases but never goes higher than K_{w}/K_{o} = 3. Considering P = 350 corresponding to the partition coefficient of thymol between oil and water, and the water and oil volume of the experiment we can see that the approximation α ≪ 1/β remains valid in the Froude number range considered.
Fig. 11 Time average Sherwood number in the steel equivalent phase as a function of the Froude number for our experiment and the experiment of Kim and Fruehan [2], Lachmund et al. [6] and Ishida et al. [3]. 
Fig. 12 Time evolution of the normalized concentration of tracer in the water phase C_{w}/C_{w,o} with different mesh resolutions at Sc_{w} = 10 and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 
Fig. 13 Time evolution of the normalized concentration of tracer in the water phase C_{w}/C_{w,0} with different mesh resolutions at Sc_{w} = 10 and Q = 3.5 l/min and with d_{inj} = 7.9 mm. 
Fig. 14 Time evolution of the Sherwood number in the water phase Sh_{w,num} with different mesh resolutions at Sc_{w} = 10 and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 
Fig. 15 Time evolution of the Sherwood number in the water phase Sh_{w,num} with different mesh resolutions at Sc_{w} = 10 and Q = 3.5 l/min and with d_{inj} = 7.9 mm. 
Fig. 16 Time evolution of the normalized concentration of tracer in the oil phase C_{o}/C_{o,0} with different mesh resolutions at Sc_{o} = 40 and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 
Fig. 17 Time evolution of the Sherwood number in the oil phase Sh_{o,num} with different mesh resolutions at Sc_{o} = 40 and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 
3.2.3 Extrapolation of numerics to large Schmidt numbers
First, we must keep in mind that to save computational time the numerical results are obtained for a maximum Schmidt number of Sc_{w} = 40. This means that we cannot directly compare our numerical results with the experimental one. Instead, we need to extrapolate our time average numerical Sherwood number to Sc_{w,th} = 1480.
On Figure 19 we plot the numerically obtained Sherwood numbers for several Schmidt numbers at a constant air flow rate of Q = 0.6 l/min. The error bars show the standard deviation of the numerical results from an asymptotic Sherwood number or time average Sherwood number Sh_{w} represented by dashed lines in Figures 14 and 15.
We can see that the three smallest Schmidt number results for the most resolved simulation (Δ_{min} = 131.75 μm) align with a Sc^{5/12} law. This power law is remarkable as it interpolates between the Sc^{1/2} law expected for mass transfer in a free surface flow and the Sc^{1/3} law expected for mass transfer in flow over a rigid (no slip) boundary. We expect the mass transfer behavior to interpolate between the Sc^{1/2} and the Sc^{1/3} law since the oil is much more viscous than water. Thus, relatively large chemical boundary layers in the water “see” a rigid oil phase while very thin boundary layers see a flowing oil phase. The interpolation would result in(36)where the coefficient b depends on the ratio of viscosities of oil and water (or slag and metal) and the coefficient a is obtained by fitting. The derivation of (36) will be published elsewhere and is still being perfected. However, it is likely that (36) can be approximated by a power law with an exponent 5/12 that is the exact average of 1/3 and 1/2.(37)where a can be determined as the numerically obtained value for Sc = 1. This scaling gives extremely good results as shown on Figure 19, as it fits both the three first numerical data points for Sc = 1, 4 and 10 and the experimental one at Sc_{w,th} = 1480. The values of the Sherwood numbers for various air flow rates extrapolated in this manner for Sc_{w,th} = 1480 are gathered in Table 5.
However, these good results could be the effect of two compensating errors: the spatial discretization error and the extrapolation error. A “safer” manner of extrapolating the low Sc results is to consider the two limiting scaling laws as upper and lower bounds so that(38)where a is the value for Sc = 1. Therefore, we plot the extrapolation of the numerical results with the two limiting scaling laws as error bars on Figure 20.
From the limited number of numerical results that we have, we cannot obviously observe both mass transfer regimes seen experimentally (cf. Fig. 20) but only the first regime. However, we notice a steep increase of the Sherwood numbers for the two largest air flow rates, indicative of the beginning of a transition to the second regime.
Finally, we want to verify the assumption made on the Section 3.2.1 concerning the role of the wateroil interfacial area in the mass transfer regime change. As we can see from the numerical results of Figure 21 when the air flow rate increases the wateroil interfacial area decreases. However, above a critical air flow rate Q > 5.5 l/min, we observe that the interfacial area reaches a value close to the interfacial area obtained at the lowest flow rate. Now, if we normalize the wateroil interfacial area with the ratio of the difference of the area of the section of the box with the time average open eye area in Figure 22. We can see that when the air flow rate is increased the normalized interfacial area increases. This is more striking for Q ≥ 5.5 l/min where we can see a steep increase. On the other hand, the values of the normalized interfacial area at low air flow rates are close to one. This indicates that despite having a bigger open eye in the oil layer, the numerous oil droplets in the water increase well the normalized interfacial area at high air flow rate compared to low flow rate. So, from the numerical results the change in the normalized interfacial area plays an important role in the increase of the Sherwood number.
Fig. 18 Evolution of the ratio of time average mass transfer coefficient in the water and oil phase as a function of the Froude number for several ratios Sc_{r} = Sc_{w}/Sc_{o}. The simulation results are obtained with a minimum grid size of Δ_{min} = 263.5 μm corresponding to a maxlevel = 10 of adaptive refinement. 
Fig. 19 Dependency of the time averaged numerical Sherwood number on the Schmidt number for different mesh refinements and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 
Fig. 20 Time averaged experimental Sherwood number (black triangles) and extrapolation at Sc_{w,th} = 1480 of the time averaged numerical Sherwood number values (purple triangles with purple error bars) as a function of the modified Froude number. 
Fig. 21 Time evolution of the numerical wateroil interfacial area for different air flow rates. The results are obtained with a constant minimum grid size of Δ_{min} = 527 μm corresponding to a maxlevel = 9 of adaptive refinement. 
Fig. 22 Time evolution of the ratio of the numerical wateroil interfacial area to the difference of the area of the section of the box with the time average open eye area for different air flow rates. The results are obtained with a constant minimum grid size of Δ_{min} = 527 μm. 
4 Conclusion
In this study, we have investigated the mass transfer between water and oil as a model of the one between liquid steel and slag in an industrial process. To reproduce this complex industrial process, we designed an experiment at ambient temperature using water and oil instead of metal and slag. We used thymol in solution as a chemical tracer transferred from the water to the oil layer. We have also designed a Direct Numerical Simulation of the experimental configuration in addition to this experiment. We have first verified that we could reproduce the hydrodynamic behavior of the industrial process by measuring the evolution of the open eye area as a function of the air flow rate. Comparison with the literature reveals that one observes areas close to those observed in other water experiments. The literature also shows that the open eye area is very strongly affected by parameters such as the fluid density. Moreover, our numerical results for the open eye area give an approximate agreement with the experimental ones, showing the same logarithmic trend and an error limited to 20%.
In the second part of this paper, we measured the mass transfer of thymol from the water to the oil phase using the same experimental configuration. Our experimental results show that two different mass transfer regimes can be identified depending on the air flow rate value. Furthermore, the transition between the two regimes coincides with the critical Froude number for the entrainment of oil droplets into the water. This agrees with the results reported in the literature.
Finally, the mass transfer of the water experiment was simulated. We used an extrapolation of our numerical results at a Schmidt number Sc = 1 to the experimental Schmidt number Sc_{w,th} = 1480. This led to some encouraging results, in agreement with the experiments although with large error bars. However, from the limited numerical results we obtained we could not clearly observe the second mass transfer regime. Nevertheless, at high air flow rate the results show something similar to a transition region between the two regimes.
For the first mass transfer regime the numerics and the experiments are well predicted by the following correlation
This linear law (exponent of unity) should be compared with the correlations found by our analysis of the literature in Figure 11, which involves exponents 1.8, 1.9 and 3.1. A definite disagreement between the various authors is noticed while our experiment (exponent 0.8) agrees with our numerics. The causes of this intriguing disagreement should be the object of future study, including numerical simulation of the other author’s experiments.
Nomenclature
Greek symbols
μ: Dynamic viscosity, [kg/m.s].
ν: Kinematic viscosity, [m^{2}/s].
δ_{s} : Dirac distribution function.
ϵ: Rate of dissipation of turbulent kinetic energy, [m^{2}/s^{3}].
η: Kolmogorov length scale, [m].
Latin symbols
n: Normal to the interface, [m^{−1}].
Q: Gas flow rate at operating condition and injector height, [m^{3}/s].
g: Gravitational acceleration, [m/s^{2}].
C: Concentration of tracer, [g/l].
D: Diffusion coefficient, [m^{2}/s].
K: Global mass transfer coefficient, [m/s].
U_{p} : Velocity of bubble plume, [m/s].
T: Numerical tracer concentration, [g/l].
P: Partition coefficient, [−].
N: Reexpression of the Froude number [−].
Sub/superscripts
Supplementary Material
Figures 1–23. Access here
Acknowledgments
This work was supported by ANRT with convention CIFRE N° 2017/0694 and benefited from access to the HPC resources of CINES under the allocations 2018 A0052B07760 and 2019 − A0072B07760 and the resources of the TGCC under the project 2020225418 granted respectively by GENCI and PRACE. Support by the ERC ADV grant TRUFLOW is acknowledged.
References
 P. Riboud, R. Vasse, Désulfuration de l’acier en poche: synthèse des résultats théoriques et industriels, Revue de Métallurgie. 82, 801–810 (1985) [CrossRef] [EDP Sciences] [Google Scholar]
 S.H. Kim, R.J. Fruehan, Physical modeling of liquid/liquid mass transfer in gas stirred ladles, Metall. Trans. B 18, 381–390 (1987) [CrossRef] [Google Scholar]
 J. Ishida, Effects of stirring by argon gas injection on metallurgical reactions in secondary steelmaking, DenkiSeiko (Electr. Furn. Steel) 52, 2–8 (1981) [CrossRef] [Google Scholar]
 J. Mietz, S. Schneider, F. Oeters, Model experiments on mass transfer in ladle metallurgy, Steel Res. 62, 1–9 (1991) [CrossRef] [Google Scholar]
 M. Hirasawa, K. Mori, M. Sano, A. Hatanaka, Y. Shimatani, Y. Okazaki, Rate of mass transfer between molten slag and metal under gas injection stirring, Trans. Iron Steel Inst. Jpn, 27, 277–282 (1987) [CrossRef] [Google Scholar]
 H. Lachmund, Y. Xie, T. Buhles, W. Pluschkell, Slag emulsification during liquid steel desulphurisation by gas injection into the ladle, Steel Res. Int. 74, 77–85 (2003) [CrossRef] [Google Scholar]
 D. Bothe, M. Koebe, K. Wielage, J. Prüss, H.J. Warnecke, Direct numerical simulation of mass transfer between rising gas bubbles and water, in Bubbly Flows (Springer, 2004), pp. 159–174 [CrossRef] [Google Scholar]
 W. Lou, M. Zhu, Numerical simulation of desulfurization behavior in gasstirred systems based on computation fluid dynamicssimultaneous reaction model (CFDSRM) coupled model, Metall. Mater. Trans. B 45, 1706–1722 (2014) [CrossRef] [Google Scholar]
 W. Lou, M. Zhu, Numerical simulation of slagmetal reactions and desulfurization efficiency in gasstirred ladles with different thermodynamics and kinetics, ISIJ Int. 55, 961–969 (2015) [CrossRef] [Google Scholar]
 R.V. Calabrese, T. Chang, P. Dang, Drop breakup in turbulent stirredtank contactors. Part I: Effect of dispersedphase viscosity, AIChE J. 32, 657–666 (1986) [CrossRef] [Google Scholar]
 R. Calabrese, C. Wang, N. Bryner, Drop breakup in turbulent stirredtank contactors. Part III: Correlations for mean size and drop size distribution, AIChE J. 32, 677–681 (1986) [CrossRef] [Google Scholar]
 C. Wang, R.V. Calabrese, Drop breakup in turbulent stirredtank contactors. Part II: Relative influence of viscosity and interfacial tension, AIChE J. 32, 667–676 (1986) [CrossRef] [Google Scholar]
 M. Iguchi, Y. Sumida, R. Okada, Z. Morita, Evaluation of critical gas flow rate for the entrapment of slag using a water model, ISIJ Int. 34, 164–170 (1994) [CrossRef] [Google Scholar]
 A. Monin, A.M. Yaglom, Statistical fluid mechanics (MIT Press, Cambridge, MA, 1971), vols. 1 and 2, p. 11 [Google Scholar]
 G.R. Hunt, T.S. van den Bremer, Classical plume theory: 19372010 and beyond, IMA J. Appl. Math. 76, 424–448 (2010) [Google Scholar]
 S. Popinet, Gerris: a treebased adaptive solver for the incompressible Euler equations in complex geometries, J. Comput. Phys. 190, 572–600 (2003) [CrossRef] [MathSciNet] [Google Scholar]
 S. Popinet, An accurate adaptive solver for surfacetensiondriven interfacial flows, J. Comput. Phys. 228, 5838–5866 (2009) [CrossRef] [MathSciNet] [Google Scholar]
 J.U. Brackbill, D.B. Kothe, C. Zemach, A continuum method for modeling surface tension, J. Comput. Phys. 100, 335–354 (1992) [Google Scholar]
 D.C. Kurt, A. Smith, F.J. Solis, A projection method for motion of triple junctions by level sets, Interfaces Free Boundaries (2002) [Google Scholar]
 X. Chen, Y. Sun, C. Xue, Y. Yu, G. Hu, Tunable structures of compound droplets formed by collision of immiscible microdroplets, Microfluid. Nanofluid. 21, 109 (2017) [CrossRef] [Google Scholar]
 J. LópezHerrera, A. GañánCalvo, S. Popinet, M. Herrada, Electrokinetic effects in the breakup of electrified jets: a volumeoffluid numerical study, Int. J. Multiphase Flow 71, 14–22 (2015) [CrossRef] [MathSciNet] [Google Scholar]
 S.B. Pope, Turbulent Flows (IOP Publishing, 2001) [Google Scholar]
 K. Yonezawa, K. Schwerdtfeger, Spout eyes formed by an emerging gas plume at the surface of a slagcovered metal melt, Metall. Mater. Trans. B 30, 411–418 (1999) [CrossRef] [Google Scholar]
 K. Krishnapisharody, G.A. Irons, Modeling of slag eye formation over a metal bath due to gas bubbling, Metall. Mater. Trans. B 37, 763–772 (2006) [CrossRef] [Google Scholar]
 M. Thunman, S. Eckert, O. Hennig, J. Björkvall, D. Sichen, Study on the formation of openeye and slag entrainment in gas stirred ladle, Steel Res. Int. 78, 849–856 (2007) [CrossRef] [Google Scholar]
 D. Bonn, J. Eggers, J. Indekeu, J. Meunier, E. Rolley, Wetting and spreading, Rev. Mod. Phys. 81, 739 (2009) [CrossRef] [Google Scholar]
 N. Joubert, Liquidliquid mass transfer characterization applied to metallurgical process, PhD dissertation, Sorbonne Université/Université Pierre et Marie Curie − Paris VI (2021). https://hal.archivesouvertes.fr/tel03227840 [Google Scholar]
Cite this article as: Nelson Joubert, Pascal Gardin, Stéphane Popinet, Stéphane Zaleski, Experimental and numerical modelling of mass transfer in a refining ladle, Metall. Res. Technol. 119, 109 (2022)
All Tables
All Figures
Fig. 1 Schematic of the experimental setup of the water ladle. 

In the text 
Fig. 2 Top view of the instantaneous experimental open eye obtained at Q = 0.6 l/min. The color code is in yellow: open eye contour; orange: spout zone. 

In the text 
Fig. 3 Ratio of time average open eye area with the square of the bath height as a function of the Froude number for our experiment and the experiments of Krishnapisharody et al. [24], Yonezawa and Schwerdtfeger [23] and Thunman et al. [25]. 

In the text 
Fig. 4 Top view of the instantaneous numerical open eye, obtained at Q = 0.6 l/min. The numerical results are obtained with a minimum grid size of Δx = 527 μm and the oil interface is colored in red. 

In the text 
Fig. 5 Time evolution of the open eye area for different air flow rates with d_{inj} = 7.9 mm. The results are obtained with a minimum grid size Δ_{min} = 527 μm corresponding to a maxlevel = 9 of adaptive refinement. 

In the text 
Fig. 6 Front view of the simulation with d_{inj} = 7.9 mm with a minimum grid size Δ_{min} = 527 μm. The air interface is colored in green, and the oil interface is colored in red. 

In the text 
Fig. 7 Front view picture of the water experiment with thymol dissolved in the water obtained with d_{inj}=7.9 mm. 

In the text 
Fig. 8 Evolution of the time averaged open eye area of several experiments with its standard deviation and numerical timeaveraged openeye area with its temporal standard deviation as a function of the air flow rate obtained with d_{inj}=7.9 mm. The simulation results are obtained with a minimum grid size of D_{min}=527 mm corresponding to a maxlevel=9 of adaptive refinement. 

In the text 
Fig. 9 Top: Time evolution of the concentration of thymol in water at Q=2.5 l/min and d_{inj}=7.9 mm. The error bars correspond to the measurement errors (cf. Supplementary Material: Fig. 1). Bottom: Time evolution of minus the righthand side of (16) obtained with Q=2.5 l/min and d_{inj}=7.9 mm. 

In the text 
Fig. 10 Evolution of the capacity coefficient or product of the time average mass transfer coefficient with the interfacial area as a function of gas flow rate for our experiment and the one of Kim et al. [2]. 

In the text 
Fig. 11 Time average Sherwood number in the steel equivalent phase as a function of the Froude number for our experiment and the experiment of Kim and Fruehan [2], Lachmund et al. [6] and Ishida et al. [3]. 

In the text 
Fig. 12 Time evolution of the normalized concentration of tracer in the water phase C_{w}/C_{w,o} with different mesh resolutions at Sc_{w} = 10 and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 

In the text 
Fig. 13 Time evolution of the normalized concentration of tracer in the water phase C_{w}/C_{w,0} with different mesh resolutions at Sc_{w} = 10 and Q = 3.5 l/min and with d_{inj} = 7.9 mm. 

In the text 
Fig. 14 Time evolution of the Sherwood number in the water phase Sh_{w,num} with different mesh resolutions at Sc_{w} = 10 and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 

In the text 
Fig. 15 Time evolution of the Sherwood number in the water phase Sh_{w,num} with different mesh resolutions at Sc_{w} = 10 and Q = 3.5 l/min and with d_{inj} = 7.9 mm. 

In the text 
Fig. 16 Time evolution of the normalized concentration of tracer in the oil phase C_{o}/C_{o,0} with different mesh resolutions at Sc_{o} = 40 and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 

In the text 
Fig. 17 Time evolution of the Sherwood number in the oil phase Sh_{o,num} with different mesh resolutions at Sc_{o} = 40 and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 

In the text 
Fig. 18 Evolution of the ratio of time average mass transfer coefficient in the water and oil phase as a function of the Froude number for several ratios Sc_{r} = Sc_{w}/Sc_{o}. The simulation results are obtained with a minimum grid size of Δ_{min} = 263.5 μm corresponding to a maxlevel = 10 of adaptive refinement. 

In the text 
Fig. 19 Dependency of the time averaged numerical Sherwood number on the Schmidt number for different mesh refinements and Q = 0.6 l/min and with d_{inj} = 7.9 mm. 

In the text 
Fig. 20 Time averaged experimental Sherwood number (black triangles) and extrapolation at Sc_{w,th} = 1480 of the time averaged numerical Sherwood number values (purple triangles with purple error bars) as a function of the modified Froude number. 

In the text 
Fig. 21 Time evolution of the numerical wateroil interfacial area for different air flow rates. The results are obtained with a constant minimum grid size of Δ_{min} = 527 μm corresponding to a maxlevel = 9 of adaptive refinement. 

In the text 
Fig. 22 Time evolution of the ratio of the numerical wateroil interfacial area to the difference of the area of the section of the box with the time average open eye area for different air flow rates. The results are obtained with a constant minimum grid size of Δ_{min} = 527 μm. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.