Issue 
Mechanics & Industry
Volume 21, Number 5, 2020
Scientific challenges and industrial applications in mechanical engineering



Article Number  517  
Number of page(s)  9  
DOI  https://doi.org/10.1051/meca/2020063  
Published online  13 August 2020 
Regular Article
Particleladen flow around an obstacle in a square pipe: experiments and modeling
^{1}
AixMarseille Univ, CNRS, IUSTI, Marseille and CEA, DES, IRESNE, DTN,
Cadarache,
France
^{2}
CEA, DES, IRESNE, DTN,
Cadarache,
France
^{3}
AixMarseille Univ, CNRS, IUSTI, Marseille and Université de Paris, CNRS, Matière et Systèmes Complexes (MSC) UMR 7057,
Paris,
France
^{4}
AixMarseille Univ, CNRS, IUSTI,
Marseille,
France
^{*} email: laurence.bergougnoux@univamu.fr
Received:
8
November
2019
Accepted:
3
June
2020
Particle trapping and deposition around an obstacle occur in many natural and industrial situations and in particular in the nuclear industry. In the steam generator of a nuclear power plant, the progressive obstruction of the flow due to particle deposition reduces the efficiency and can induce tube cracking leading to breaking and damage. The steam generator then loses its role as a safety barrier of the nuclear power plant. From a fundamental standpoint, dilute and concentrated particulate flows have received a growing attention in the last decade. In this study, we investigate the transport of solid particles around obstacles in a confined flow. Experiments were performed in a simplified configuration by considering a laminar flow in a vertical tube. An obstacle was inserted at the middle height of the tube and neutrallybuoyant particles were injected at different locations along the tube. We have investigated first the trajectories of individual particles using particle tracking (PT). Then, the particle trajectories were modeled by using the BoussinesqBassetOseen equation with a flow velocity field either measured using particle image velocimetry (PIV) or calculated by the Code_Saturne software in order to account for the threedimensional (3D) character of the obstacle wake. This paper presents a comparison between the experimental observations and the predictions of the modeling for an obstacle consisting of a rectangular step at a Reynolds number of ≈100 and evidences the importance of accounting for the 3D complex nature of the flow.
Key words: Particleladen flow / obstacle / wakes
© O. Ait Oucheggou et al., published by EDP Sciences 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://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
About 11% of the world electricity is generated by 450 nuclear power reactors. Most of these nuclear power plants (65% of them) are Pressurized Water Reactors (PWR). The Steam Generator (SG) is a crucial component of the PWR where the heat exchange between the primary and secondary circuits produces steam which is used to drive a turbine linked to an alternator to produce electricity. The tube bundle in the SG comprises 3300 to 6000 tubes supported at different heights by internal structures called Tube Support Plates (TPS). Clogging, i.e. the progressive obstruction mainly due by deposits of iron oxides, occurs in the quatrefoil holes between the tubes and the TPS designed for the circulation of the fluid. This phenomenon can have important consequences on the safety of the PWR. Therefore, the understanding of the origin of this mechanism [1,2], especially that related to thermal, chemical, hydrodynamic, and physical issues, requires particular attention. Investigations inside the PWRs are limited because of the difficulty to access to specific location of this phenomenon. The difficulty of undertaking experiments and instrumentation under thermohydraulic (pressure ≃6 MPa, temperature ≃280°C) and chemical representative conditions in the SG justifies to simplify the problem and to perform model experiments. The objective of the present study is thus to examine the flow of particles around obstacles (which mimic the TPSs) and in particular to study their preferential locations of trapping. To simplify further the problem, we have chosen to focus on the hydrodynamic effects and to neglect other phenomena such as thermal and chemical interactions.
From a fundamental standpoint, the flow around a bluff body is controlled by several parameters: the shape of the obstacle, the environment in which the body is located (a confined pipe [3] or a free surface channel [4]), and the flow Reynolds number, Re_{f}. Pressure losses caused by a sudden change in the flow section give rise to recirculation zones. This is the case of a backwardfacing step flow which has been extensively studied and used as a test case in numerical simulations. Experimental investigations of laminar, transitional, and turbulent flows with regions of separation behind a twodimensional backwardfacing step were carried out by Armaly et al. [5] with an air flow in a pipe with a blockage ratio of 0.5. The length of the wake behind the step was measured, as well as the lengths and locations of the other recirculatingflow regions. One of the rare studies focusing on the upstream region of an obstacle has been conducted by Hattori et al. [6]. This numerical study shows that the size of the recirculation zone upstream of the forwardfacing step becomes smaller with increasing Re_{f}. Another result is that the obstacle blockage in the pipe influences the sizes of the recirculation zones.
Dilute and concentrated particulate flows have received a growing attention in the last decade. For example, the study conducted by Matas et al. [7] on dilute suspensions of neutrally buoyant particles in Poiseuille flow extended that of Segrè and Silberberg [8] and showed that particles accumulate on a narrow annulus which moves toward the wall as Re_{f} increases. The same authors [9] investigated the influence of suspended particles on the transition to turbulence in a pipe flow. Their results indicate that particle sizes and concentration influence this transition: small particles cause a delay in transition to larger critical Reynolds number while, for larger particles at low concentrations, the opposite effect is observed.
In the case of suspension flows in obstructed pipes, the studies are mostly focused on the downstream side of the obstacle, although the area affected by clogging is the upstream region of the obstacles in the practical case of steam generators. The behavior of monodisperse and neutrally buoyant particles in the wake of a cylinder has been examined by Haddadi et al. [10] at moderate Re_{f} < 30. Their experimental observationsin a microfluidic device indicated the formation of a particledepleted region in the wake of the obstacle. Using simulations based on a latticeBoltzmann method, they investigated the behavior of a single particle placed at different initial locations. They found that, when released inside the wake region, the particles move towards the boundary of the wake until reaching a limitcycle trajectory. When many particles are introduced into the wake, because of particle hydrodynamic interactions, particles can move from one lobe of the wake to the other, or completely exit the wake by escaping into the free stream. Recently, the same authors [11] examined numerically the trapping of particles in the wake of obstacles having different geometries by varying the size of particles and the density ratio between the particles and the fluid. They showed that the vortex depopulation kinetics depends upon the shape of the obstacle as well as upon the density ratio of the particle and the fluid and the size ratio between the particles and the obstacle.
The particle trapping or deposition phenomenon in particleladen flow inside confined geometries is still poorly understood. In this paper, an experimental investigation of the behavior of spherical particles around an obstacle in laminar flow is presented. Experimental particle trajectories are compared to predictions obtained by using the equation of motion for a small rigid sphere at low Reynolds numbers [12,13].
2 Materials and methods
2.1 Dimensional analysis
Before presenting the experimental setup, the different physical and geometrical quantities of the study are discussed. The square pipe is characterized by its hydraulic diameter, d_{h}, defined as the ratio of four times the crosssectional area of the flow and the wetted perimeter of this crosssection. The fluid flowing through the tube is characterized by its viscosity, μ, its density, ρ_{f}, and the inlet velocity, u_{0}. The particles are spherical with a radius, a_{p}, and a density, ρ_{p}. An additional quantity is the gravitational acceleration, g. We thus have seven independent physical quantities. Using the Buckingham π theorem, dimensionless analysis indicates that the present problem is characterized by four dimensionless parameters: the velocity ratio, W = v_{St}∕u_{0}, where is the Stokes velocity of particles, the Stokes number, , characterizing the behavior of particles suspended in a fluid flow (a particle with a low Stokes number follows fluid streamlines while a particle with a large Stokes number is dominated by its inertia), the density ratio, R = ρ_{f}∕(ρ_{p} + ρ_{f}∕2), and the size ratio, P = a_{p}∕d_{h}, which is always rather small. By combining these four parameters, two Reynolds numbers can be obtained: the Reynolds number of the particles, , and the Reynolds number of the fluid, flow, Re_{f} = u_{0}ρ_{f}d_{h}∕μ = Re_{p}∕(PW).
2.2 Experimental setup
Experiments are performed in a transparent Plexiglass^{®} pipe having an inner square crosssection of 2.5 ×2.5 cm^{2} with d_{h} = 2.5 cm, and a filled height of ≈100 cm obstructed by a rectangular obstacle inserted at the middle height of the tube, as sketched in Figure 1. All the experiments are performed at the ambient room temperature of 23 ± 1 °C. The fluid is a mixture of water and saccharose (13.5 %) with a density ρ_{f} = 1049 kg.m^{−3} and a viscosity μ = 1.5 mPa.s at 23± 1 °C. The upwardflow in the pipe is provided by a pump (SIMATEC/MVPUMPSYSTEM) covering a range of flow up to Re_{f} = 1000. The inlet fluid velocity is fixed at u_{0} = 0.54 cm.s^{−1}, i.e. Re_{f} = 94. Spherical particles made of polystyrene having a radius a_{p} = 250 ± 5 μm and density ρ_{p} = 1050 kg.m^{−3} are used.
A digital camera (JAI SP12000MCXP4FXT, 4096 × 700 pixels, 120 fps) placed in front of the tube is focused on the obstacle and the particles. The flow characterization is performed by a nonintrusive method, particle image velocimetry (PIV) using the Matlab ™PIV software DPIVsoft [14]. The fluid is seeded with hollow particles used as fluid tracers (Dantec Measurement Technology) having a diameter ≈ 15 μm and a density ≈ 1400 kg.m^{−3}. For all PIV measurements, the spatial resolution is 0.13 × 0.13 cm^{2}, corresponding to an interrogation window of 32 × 32 pixels with an overlap of 50%. The error in the measurements of the timeaveraged velocity fields was estimated to remain below 10%.
Particle tracking (PT) is used to track the motion of particles and infer their velocities. A neon light illuminates the region around the obstacle and the digital camera records the particle motions at 15 frames/s. The particle tracking scheme uses the plugin MTrackJ of the ImageJ software [15]. The experimental error in the particle center measurement is ≈ 1∕5 of the particle diameter. Since the camera is at the same location in the PIV and PT measurements, the coordinates of the particle can be accurately positioned inside the velocityvector map of the flow as a function of time.
Fig. 1 Schematic view of the experimental setup: the pipe has an inner square crosssection of 2.5 × 2.5 cm^{2} with d_{h} = 2.5 cm, and a filled height of ≈100 cm obstructed by a rectangular obstacle inserted at the middle height of the tube. 
2.3 Simulation with Code_Saturne
Code_Saturne is a free, opensource software developed and released by EDF to solve computational fluid dynamics (CFD) applications. It solves the NavierStokes equations for 2D, 2Daxisymmetric, and 3D flows, steady or unsteady, laminar or turbulent, incompressible or weakly dilatable, isothermal or not, with scalar transport if required. In this study, a 3D steady simulation of the flow is performed. A structured hexahedral mesh is used with cell sizes of 1 × 1× 1 mm^{3} along the tube dimensions. The mean flow is in the zdirection (opposite to the gravity direction) with an inlet uniform velocity u_{0} = 0.54 cm.s^{−1}. Computational convergence is verified by following the velocities of the flow (u_{x}, u_{y}, u_{z}) throughout the iterations. Convergence is obtained after 200 iterations.
2.4 Modeling
We follow the same approach as that developed by Bergougnoux et al. [16] in their study of particle sedimentation in a cellular flow field at low Stokes numbers (<0.1). Bergougnoux et al. performed experiments with Poly(methyl methacrylate) (PMMA) spherical particles sedimenting in a cellular flow. A model inspired by that of Maxey [17] was developed by considering the additional effect of the particle history force but also the modification of the drag force of the particle when the particle Reynolds Re_{p} increases. This modeling is based on the resolution of the BoussinesqBassetOseen equation (BBO) including Faxén corrections to consider the influence of the particle size [12,13]. The BBO equation predicts the acceleration dv ∕dt and the velocity v of the particle knowing its initial position and the velocity field u = (u_{x}, u_{y}, u_{z}) of the fluid: (1)
where and , are respectively the particle mass and the mass of the displaced fluid. In this equation (1), d∕dt denotes the time derivative in a frame moving withthe particle, while D∕Dt the derivative following the fluid element. The forces included in the righthand side of equation (1) are the buoyancy force, the Stokes drag,the fluid acceleration, and the added mass effect. The history BoussinesqBasset force has been neglected in the present study. The terms in are the Faxén corrections.
We can consider that the flow field u is steady and write the dimensionless form (*) of equation (1), using u_{0} as the velocity scale and d_{h} as the length scale, (2)
For the selected combination of particles and fluid, we have W = 0.019, St = 0.003, R = 0.666, and P = 0.011. The Faxén corrections, proportional to P^{2}, are negligible and the dragbuoyancy term, inversely proportional to St, largely prevails over the other terms as also found by Bergougnoux et al. [16].
3 Results and discussion
3.1 Flow characterization
3.1.1 Experimental characterization
The flow characterization is obtained experimentally using PIV in a plane (x^{*}, z^{*}) by considering 500 pairs of snapshots. From each pair of snapshots, one instantaneous velocity field can be determined. In this steadystate regime, the timeaveraged velocity field , in the center section (y^{*} = 0.5) is presented in Figure 2.
In this Figure 2, recirculation zones are formed upstream and downstream from the obstacle which is a classical result found in the literature. There are some discontinuity zones in the velocity field before reaching the obstacle, due to a lack of seeding in these regions. The fluid is accelerating when passing above the obstacle (u^{*} ≈ 4u_{0}) due to flow rate conservation. In the recirculation zones, the fluid velocity is weak: particles can be trapped in these areas. In the wake of the obstacle, the main recirculation zone is clearly visible.
Fig. 2 Velocity field () measured by PIV in the (x^{*}, z^{*}) plane at y^{*} = 0.5. 
3.1.2 Numerical characterization
Because of the complexity of the flow around the obstacle, access to the third velocity component of the flow, u_{y}, is essential. The 3D flow characterization is obtained using Code_Saturne. The velocity field () computed in the same plane (y^{*} = 0.5) as that of the experimental result is plotted in Figure 3 and shows a decent agreement with the experimental velocity field of Figure 2.
Figure 4 presents sliced views of the flow field. It exhibits the Poiseuille velocity profile in the first two planes before the obstacle and the two last planes towards the outlet. The flow complexity and the importance of the velocity field () in different planes (x^{*}, y^{*}) are reported in Figure 5. This last figure shows that the component is not negligible. Clearly, there is a three dimensional structure of the flow as evidenced for instance on the top of the obstacle for graph (d). The numerical 3D flow obtained by the Code_Saturne will then be used for solving equation (2) and comparing to the experiments.
Fig. 3 Velocity field calculated by Code_Saturne in the (x^{*}, z^{*}) plane, y^{*} = 0.5. 
Fig. 4 Slices of velocity calculated by Code_Saturne in different (x^{*}, y^{*}) planes. 
Fig. 5 Velocity field calculated by Code_Saturne in different (x^{*}, y^{*}) planes at (a) z^{*} = −0.78, (b) z^{*} = −0.22, (c) z^{*} = −0.02, (d) z^{*} = 0.42, (e) z^{*} = 1.46, (f) z^{*} = 2.10, (g) z^{*} = 2.50, (h) z^{*} = 3.22 and (i) z^{*} = 5.62. 
3.2 Particle trajectories
3.2.1 Experimental trajectories
Figure 6 presents three experimental trajectories obtained by PT as well as the results of the modeling with equation (2) when using the velocity field given by PIV measurements. The particles are experimentally injected at different times. The first particle (□, . ) released in the free stream passes over the obstacle and remains away from the recirculation regions. The comparison with the modeling shows a good agreement. Conversely, for particles injected in the downstream wake (△, ) or upstream of the obstacle (○, ), the model cannot capture properly the trajectories because of the 3D character of the flow in these recirculation zones which is not given by the 2D PIV flow field.
Fig. 6 Experimental trajectories of three particles (□, △, ○) and modeling results using the PIV flow field (. , , ). 
3.2.2 Modeling using the numerical 3D flow field
In this section, we now use the 3D flow calculated by Code_Saturne to solve equation (2). The new predictions are compared with the experimental trajectories presented in Figure 6.
Figure 7 presents the results for the particle released in the free stream. It is mostly important for the numerical modeling to fix the coordinate y^{*} at the initial time of the computation. We have taken y^{*} = 0.35 () from t^{*} = 0 and y^{*} = 0.32 () from t^{*} ≈ 0.5. The predictions with P = 0, R = 0, W = 0.019, and St = 0.003 are in excellent agreement with the experimental data. We have also calculated the first part of the trajectory with P = 0.011 and R = 0.666 (). The comparison with the trajectory using P = 0 and R = 0 () shows that these two parameters do not affect the particle trajectory. Therefore, the model is simplified by neglecting these terms in equation (2). We should also note that we computed the trajectory of a fluid particle having the same initial conditions and that the particle does not deviate significantly from being a fluid tracer.
We proceed similarly for the second particle injected in the wake shown in Figure 8. As soon as the model can no longer follow the experimental trajectory (due to the complexity of the fluid in the wake of the obstacle), we adjusted the value of the y^{*} component. Good agreement between the predictions and the experimental data is found. The trajectory of a fluid particle has been also computed ( ) with same initial positions. Clearly, the particle does not simply follow the flow lines and is not a simple fluid tracer.
The trajectory of the third particle injected upstream of the obstacle is shown in Figure 9. To be able to follow the particle, we had to change three times the plane y^{*} (, , ). Experimentally, the particle turns once when arriving in the recirculation zone upstream of the obstacle and then leaves by bypassing it. Figure 9d shows that the particle, in addition to a change of position in the observation plane (x^{*}, z^{*}), also moves in the direction y^{*}. To follow the particle and to predict correctly its trajectory, it is essential to have a measure of its injection plane. Again this particle is not a simple fluid tracer, see the ( ) curve.
Fig. 7 (a) Experimental trajectory (□) and modeling prediction (, ) of particle 1; Temporal evolution of particle position (b) z^{*}, (c) x^{*} and (d) y^{*}; Particle velocity (e) , (f) and (g) . 
Fig. 8 (a) Experimental trajectory (△) and modeling prediction (, ) of particle 2; Temporal evolution of particle position (b) z^{*}, (c) x^{*} and (d) y^{*}; Particle velocity (e) , (f) and (g) . 
Fig. 9 (a) Experimental trajectory (○) and modeling prediction (, ) of particle 3; Temporal evolution of particle position (b) z^{*}, (c) x^{*} and (d) y^{*}; Particle velocity (e) , (f) and (g) . 
4 Conclusion
In this study, we have examined experimentally the dynamics of quasineutrallybuoyant particles around an obstacle. We have also used a model for the Lagrangian motion of particles at low Stokes number, St = 3 10^{−3}, and for Re_{f} ≈ 100. The predicted trajectories have been compared to the experimental observations. The 3D character of the flow makes the problem complex, in particular in the recirculation zones upstream and downstream of the obstacle. Decent agreement between experiments and modeling can be achieved when the 3D nature of the flow is accounted for in a very fine manner.
Nomenclature
d_{h}: Hydraulic diameter of tube
u, v: Fluid and particle velocity
v_{St}: Particle Stokes velocity
Re_{f}: Reynolds number of the flow
Acknowledgements
The authors thank P. Cervetti, S. Noel, and F. Ratouchniak for technical assistance in the experimental setup realization. The work was supported in part by CEA. The thesis of O. Ait Oucheggou benefits from a joint funding by CEA and AixMarseille Université (contract AMUCNRSCEA V4322.001).
References
 G. Yang, V. Pointeau, E. Tevissen, A. Chagnes, A review on clogging of recirculating steam generators in pressurized water reactors, Prog. Nucl. Energy 97, 182–196 (2017) [Google Scholar]
 H. Rummens, J.T. Rogers, C.W. Turner, The thermal hydraulics of tube support fouling in nuclear steam generators, Nucl. Technol. 148, 268–286 (2004) [Google Scholar]
 L. Klotz, S. GoujonDurand, J. Rokicki, J.E. Wesfreid, Experimental investigation of flow behind a cube for moderate Reynolds numbers, J. Fluid Mech. 750, 73–98 (2014) [Google Scholar]
 J.C.R. Hunt, C.J. Abell, J.A. Peterka, H. Woo, Kinematical studies of the flows around free or surfacemounted obstacles; applying topology to flow visualization, J. Fluid Mech. 86, 179–200 (1978) [Google Scholar]
 B.F. Armaly, F. Durst, J.C.F. Pereira, B. Schönung, Experimental and theoretical investigation of backwardfacing step flow, J. Fluid Mech. 127, 473–496 (1983) [Google Scholar]
 H. Hattori, Y. Nagano, Investigation of turbulent boundary layer over forwardfacing step via direct numerical simulation, Int. J. Heat Fluid Flow 31, 284–294 (2010) [Google Scholar]
 J.P. Matas, J. F. Morris, É. Guazzelli, Inertial migration of rigid spherical particles in Poiseuille flow, J. Fluid Mech. 515, 171–195 (2004) [Google Scholar]
 G. Segré, A. Silberberg, Behaviour of macroscopic rigid spheres in Poiseuille flow Part 2. Experimental results and interpretation, J. Fluid Mech. 14, 136–157 (1962) [Google Scholar]
 J.P. Matas, J. F. Morris, É. Guazzelli, Influence of particles on the transition to turbulence in pipe flow, Philos. Trans. Roy. Soc. A 361, 911–919 (2003) [Google Scholar]
 H. Haddadi, S. ShojaeiZadeh, K. Connington, J.F. Morris, Suspension flow past a cylinder: particle interactions with recirculating wakes, J. Fluid Mech. 760, R2 (2014) [Google Scholar]
 H. Haddadi, S. ShojaeiZadeh, J.F. Morris, LatticeBoltzmann simulation of inertial particleladen flow around an obstacle, Phys. Rev. Fluids 1, 024201 (2016) [Google Scholar]
 R. Gatignol, The Faxén formulae for a rigid particle in an unsteady nonuniform Stokes flow, J. Mécaniq. Théor. Appl. 2, 143–160 (1983) [Google Scholar]
 M.R. Maxey, J.J. Riley, Equation of motion for a small rigid sphere in a nonuniform flow, Phys. Fluids 26, 883–889 (1983) [Google Scholar]
 P. Meunier, T. Leweke, Analysis and treatment of errors due to high velocity gradients in particle image velocimetry, Exp. Fluids 35, 408–421 (2003) [Google Scholar]
 E. Meijering, O. Dzyubachyk, I. Smal, Methods for cell and particle tracking, Methods Enzymol. 504, 183–200 (2012) [PubMed] [Google Scholar]
 L. Bergougnoux, G. Bouchet, D. Lopez, É. Guazzelli, The motion of solid spherical particles falling in a cellular flow field at low Stokes number, Phys. Fluids 26, 093302 (2014) [Google Scholar]
 M.R. Maxey, The motion of a small spherical particles in a cellular flow field, Phys. Fluids 30, 1915–1928 (1987) [Google Scholar]
Cite this article as: O. Ait Oucheggou, V. Pointeau, G. Ricciardi, É. Guazzelli, L. Bergougnoux, Particleladen flow around an obstacle in a square pipe: Experiments and modeling, Mechanics & Industry 21, 517 (2020)
All Figures
Fig. 1 Schematic view of the experimental setup: the pipe has an inner square crosssection of 2.5 × 2.5 cm^{2} with d_{h} = 2.5 cm, and a filled height of ≈100 cm obstructed by a rectangular obstacle inserted at the middle height of the tube. 

In the text 
Fig. 2 Velocity field () measured by PIV in the (x^{*}, z^{*}) plane at y^{*} = 0.5. 

In the text 
Fig. 3 Velocity field calculated by Code_Saturne in the (x^{*}, z^{*}) plane, y^{*} = 0.5. 

In the text 
Fig. 4 Slices of velocity calculated by Code_Saturne in different (x^{*}, y^{*}) planes. 

In the text 
Fig. 5 Velocity field calculated by Code_Saturne in different (x^{*}, y^{*}) planes at (a) z^{*} = −0.78, (b) z^{*} = −0.22, (c) z^{*} = −0.02, (d) z^{*} = 0.42, (e) z^{*} = 1.46, (f) z^{*} = 2.10, (g) z^{*} = 2.50, (h) z^{*} = 3.22 and (i) z^{*} = 5.62. 

In the text 
Fig. 6 Experimental trajectories of three particles (□, △, ○) and modeling results using the PIV flow field (. , , ). 

In the text 
Fig. 7 (a) Experimental trajectory (□) and modeling prediction (, ) of particle 1; Temporal evolution of particle position (b) z^{*}, (c) x^{*} and (d) y^{*}; Particle velocity (e) , (f) and (g) . 

In the text 
Fig. 8 (a) Experimental trajectory (△) and modeling prediction (, ) of particle 2; Temporal evolution of particle position (b) z^{*}, (c) x^{*} and (d) y^{*}; Particle velocity (e) , (f) and (g) . 

In the text 
Fig. 9 (a) Experimental trajectory (○) and modeling prediction (, ) of particle 3; Temporal evolution of particle position (b) z^{*}, (c) x^{*} and (d) y^{*}; Particle velocity (e) , (f) and (g) . 

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.