Liquid blending and solid-liquid mixing are important operations in the production processes of various industries such as chemical, pharmaceutical, and food. A typical method for the blending and mixing is agitation of a fluid in a vessel with impellers. Various techniques, such as visualization of mixing process  and velocity measurement by particle image velocimetry    , have been used to clarify the mixing mechanism and improve mixing performance. Theoretical consideration  and numerical simulations of mixing performance  -  were also performed using computational fluid dynamics (CFD), which can be applied to prediction in practical problems. Although agitation is a typical mixing technique, there is a demand for blending and mixing without impellers and agitators for mixing liquids with higher viscosities.
A blade-free planetary mixer is a new technology for mixing liquids without impellers and agitators. This type of mixer can mix highly viscous fluids with less contamination as the flow in the vessel is induced by precession that combines rotation and revolution of the vessel. The fluid mechanics in a precessing cylinder, which is the basic principle of the blade-free planetary mixer, has been studied in connection with the instability of the fluid motion  . However, the design methodology of this type of mixer has not been thoroughly studied in literature, which may be due to the complexity of the flow field of the blade-free planetary mixer even in a laminar flow regime. Only a paper in literature is a numerical simulation of flow field in the blade-free planetary mixer of laminar flow  . It reported that spiral vortical flows were generated in the vessel, and their sizes, shapes, and rotation axes were influenced by the precession rate, which is the ratio of the revolution speed to the rotation speed. Nevertheless, the influence of the precession rate on laminar mixing performance has not been clarified, because it is difficult to evaluate the mixing performance only from the flow field information.
Two approaches to evaluate mixing performance with CFD have been reported in the literature: Eulerian-Lagrangian method      and Eulerian-Eulerian method    . The Eulerian-Lagrangian method was employed in this study, because this method is desirable for advection-dominated problems, and it can provide the trajectory of individual particles in the mixing process. An example of the Eulerian-Lagrangian method is the mixing performance analysis in a precessing sphere  . It was found that the flow mixing was enhanced under weak precession in the turbulent flow regime.
The objective of this study is to investigate the laminar mixing performance by numerical simulation in a blade-free planetary mixer, which is related to flow structures induced by the combined influence of rotation and revolution. The influence of the precession rate on the flow structure and mixing performance was also investigated.
2. Numerical Methods
2.1. Flow Field
Flow structures in a blade-freeplanetary mixer vessel were investigated by three-dimensional numerical simulations  . The simulations were carried out with the SRFPimpleFoam solver in OpenFOAM 4.1. The solver used in this study was an unsteady-state solver for incompressible Newtonian fluid flows including a rotating reference frame in the computational domain. However, the target flow field was the steady state laminar flow of blade-free planetary mixer. Fluid motion in a rotating frame is expressed by Navier-Stokes equations with the Coriolis and centrifugal forces acting as external forces, as well as continuity equation. The governing equations are written as follows:
where t, u, p, ρ, and ν denote time, fluid velocity, pressure, fluid density, and viscosity, respectively. Note that the fluid velocity and position vectors are defined as the relative velocity and position vectors in the rotating frame. The third term on the right-hand side of Equation (1) is the Coriolis force, where Ω is the revolution speed of the rotating frame. The fourth term is the centrifugal force, where r represents the position vector from the revolution axis.
Figure 1 shows the computational configuration of a blade-free planetary mixer. The mixing vessel has cylindrical geometry with an inner diameter of d = 76 mm and a height of h = 77 mm. The rotation axis of the vessel was inclined 45 degrees with respect to the revolution axis. The revolution radius was Rrev = 83 mm at the bottom center of the vessel. The configuration of the mixer was determined based on an actual machine designed for the purpose of experimental investigation. The rotation speed of the vessel was set to nrot = 60 rpm, and the fluid was silicone oil with kinematic viscosity of ν = 1.5 × 10−5 m2/s to keep the flow regime laminar in the range of the computational conditions. Therefore, the Reynolds number was Re ,where ωrot is the angular velocity. Note that the present Reynolds number is smaller than the critical Reynolds number of 1 × 104 in a stirred cylindrical vessel with weak precession  . Simulations were performed with revolution speeds of nrev = 0, 15, 30, 45, and 60 rpm, and the precession rates in the simulation were Γ (= nrev/nrot) = 0, 0.25, 0.5, 0.75, and 1. The rotation and revolution directions were opposite.
The computational domain was in the cylindrical vessel which was filled with the fluid. At the wall boundaries of the vessel, the rotation speed was applied to define the rotational motion of the vessel. The governing equations were discretized using the finite volume method with the quadratic upstream interpolation for convective kinematics (QUICK) scheme for the convective term and with the central differencing scheme for the other terms. The pressure implicit with splitting of
Figure 1. Computational configuration and grid of a planetary centrifugal mixer.
operator (PISO) algorithm was used to solve the governing equations. Independence of the grid convergence was evaluated using fine, medium, and coarse grids with 6.5 × 105, 3.8 × 105, and 4.8 × 104 cells, respectively. Root-mean-square errors for the velocity fields based on the fine grid were 0.1% and 1.8% for the medium grid and the coarse grid, respectively. Consequently, the medium grid shown in Figure 1 was used for numerical simulations in this study.
2.2. Mixing Performance
In this study, the mixing performance of the blade-free planetary mixer was evaluated by the Eulerian-Lagrangian method. The degree of mixing can be estimated from the number density distribution of particles by tracking a sufficient number of the particles in the vessel. Particle trajectories were simulated with the Lagrangian particle tracking method as passive particles after flow in the vessel reached the steady state. In this study, for experimental validation, it was assumed that solid particles gravitated to the bottom before mixing and were transported throughout the vessel by mixing. The particles were spherical with 10 μm diameter, and the particle mass density was equal to the mass density of the fluid in order to evaluate the mixing performance of the flow field itself. Collisions between particles were ignored.
The mixing performance was evaluated by the mixing index, which was calculated based on the particle number density in the vessel  . The entire vessel was divided into N sub-domains. The relative particle number density ρi was defined as the ratio of the particle number density in the ith sub-domain to the maximum particle number density in all sub-domains. The standard deviation of the number density σ over all sub-domains is defined as follows:
where, ρm represents the average particle number density in the entire vessel. The mixing index Im is given by the following equation:
In this study, the number of sub-domains was N = 6000 and the number of supplied particles was 6000 per second.
3. Results and Discussion
3.1. Flow Structure
Velocity fields in a blade-free planetary mixer vessel were computed to understand the mixing mechanism. Figure 2 shows three-dimensional flow structures and cross-sectional velocity fields in the mixer vessel for a rotation speed of nrot = 60 rpm and precession rates Γ (= nrev/nrot) =0 (a), 0.25 (b), 0.5 (c), and 1 (d). The results for Γ = 0.75 are omitted because they are similar to the results for Γ = 1. Streamlines (left) are shown to analyze the three-dimensional flow structures, and velocity magnitude contours of the in-plane velocity components are shown
Figure 2. Flow structures and cross-sectional velocity fields: (a) nrev = 0 rpm (Γ = 0), (b) nrev = 15 rpm (Γ = 0.25), (c) nrev = 30 rpm (Γ = 0.5), and (d) nrev = 60 rpm (Γ = 1). Each column shows the streamlines (left), velocity contours in the x-z plane (middle), and in the y-z plane (right).
in the x-z plane (middle) and in the y-z plane (right). The velocity component was normalized by the wall velocity U0 (= ). At this Reynolds number, a steady flow was observed after a certain transient state. The results from panel (a) show circular streamlines parallel to the bottom and no in-plane velocity component due to the two-dimensional rotating flow. In the precession cases, the streamlines show the generation of vortical flows with axes inclined with respect to the spin axis of the vessel. The velocity fields for the lower precession rate of Γ = 0.25 (b) show strong rotational flow in the vertical cross-section, but the flow structure was nearly two-dimensional. The results for Γ = 0.5 (c) show a spiral axis of the rotating flow and a complex vortical flow structure. Non-axisymmetric rotating flow was observed in the x-z plane, and twin corner vortices occurred near the upper left and lower right corners in the y-z plane. As the precession rate increases further, the flow structure for Γ = 1 (d) became more complex. However, the low-speed region extended around the vortical flow axis. These characteristics of the vortical flow structure with the spiral and inclined axis in the blade-free planetary mixer are similar to the results from others presented in the literature  . The occurrence of the rotating flows in the vertical cross-section of the cylindrical vessel is caused by the inclination of the axis of the rotating flow, which is caused by the Coriolis force due to revolution. It is noted that the rotating flow does not occur in the cylindrical vessel without revolution. With increasing the precession rate, the flow around the axis of the rotating flow is highly magnified due to the instability of the rotating flow  . Therefore, the complex flow pattern due to the Coriolis force is generated in the vessel and results in the enhancement of flow mixing in the blade-free planetary mixer.
3.2. Mixing Performance
The mixing performance was investigated by tracking tracer particles in the vessel and evaluating the particle number density. Figure 3 shows the temporal evolution of particle mixing in the vertical (x-z) plane. The particle positions were calculated based on the velocity field at different precession rates of (a) Γ = 0, (b) 0.25, (c) 0.5, and (d) 1. Particles were supplied near the bottom surface. The degree of mixing converged to a steady state in any cases by t = 200 s after mixing began. In the case of only rotation (a), the particles remained at the bottom of the vessel due to the two-dimensional rotational flow. In the case of Γ = 0.25 (b), some particles moved to the upper region along the vessel walls by the rotating flow induced through the combined effect of rotation and revolution. However, the particles continued to move in the outer area along the vessel walls. This result indicates that although the rotating flow axis was inclined with respect to the vessel spin axis, the flow in the vessel was a two-dimensional rotating flow. As the revolution speed increased, the particles were dispersed throughout the vessel and particle mixing was maximized for Γ = 0.5 (c). The relatively high concentration of particles was observed in the lower part of the vessel due to the supply of particles from the bottom surface. As the revolution speed increased further, unmixed areas were observed for Γ = 1 (d) due to the low-speed rotating flow region in the y-z plane (Figure 2).
To quantify the mixing performance, the mixing index was evaluated from the particle number density at different precession rates. Figure 4 summarizes the mixing index in the steady state as the precession rate varies. As shown in the figure, the precession rate for Γ = 0.5 shows the highest mixing index value. The mixing index at Γ = 0.25 was comparatively lower, though strong rotating flow occurred in the vertical plane. At higher precession rates, the mixing index gradually decreased due to the occurrence of the unmixed region, since the low-speed region extended around the rotating flow axis. These results indicate that the moderate precession enhances the mixing performance by the magnified rotating flow, but too strong precession prevents the formation of the large-scale vortical flow. Therefore, an optimum precession rate, which is around Γ =0.5, is found to enhance the mixing performance of a blade-free planetary mixer. It should be mentioned that the optimum precession rate in a precessing sphere operating in a turbulent flow  was reported at Γ = 0.1. Therefore, the optimum precession rates may be affected by the influence of container shapes, precession angles, Reynolds numbers and the laminar or turbulent flow. Therefore, further investigation is needed to clarify the influence of these parameters on the optimum precession rate.
Figure 3. Temporal evolution of mixing of tracer particles from the bottom of the vessel: (a) nrev = 0 rpm (Γ = 0), (b) nrev = 15 rpm (Γ = 0.25), (c) nrev = 30 rpm (Γ = 0.5), and (d) nrev = 60 rpm (Γ = 1).
Figure 4. Mixing index versus precession rate.
The flow structures and the mixing performance of a blade-free planetary mixer were investigated by numerical simulations. The flow fields in the mixer vessel were determined by solving the Navier-Stokes equation considering the influence of the Coriolis and centrifugal forces in a rotating reference frame. The mixing process was numerically visualized by a particle tracking method, and the mixing index was evaluated from the concentration of particles. In the numerical results, a vortical structure with an axis inclined with respect to the rotation axis was induced by the combined influence of rotation and revolution at lower precession rates. This vortical flow generated a rotating flow in the vertical cross-section of the vessel. As the precession rate increases, vortical flows showed complex three-dimensional flow structures with a spiral rotation axis. It is found from the present numerical simulation that the laminar mixing performance of the blade-free planetary mixer becomes maximum around the precession rate 0.5, which is contributed by the occurrence of rotating flow with an inclination angle to the spin axis and the generation of highly magnified flow around the axis caused by the precession. However, too strong precession prevents the formation of the large-scale vortical flow.
The authors express thanks to Mitsuboshi Co. Ltd. for the help and suggestions during the course of this study.
d vessel diameter [mm]
h vessel height [mm]
Im mixing index
N number of sub-domains
nrev revolution speed [rpm]
nrot rotation speed [rpm]
p pressure [Pa]
Re Reynolds number
Rrev revolution radius [mm]
r position vector from a revolution axis
t time [s]
U0 wall velocity [m/s]
u velocity vector [m/s]
u, v, w velocity components [m/s]
x, y, z Cartesian coordinates
Γ precession rate (=nrev/nrot)
ν kinematic viscosity [m2/s]
ρ fluid density [kg/m3]
ρi particle number density in the ith sub-domain
ρm averaged number density
σ standard deviation of number density
Ω revolution rate [rad/s]
ωrot rotation rate [rad/s]
 Wang, S., Wu, J. and Bong, E.Y. (2013) Reduced IMRs in a Mixing Tank via Agitation Improvement. Chemical Engineering Research and Design, 91, 1009-1017.
 Aubin, J., Le Sauze, N., Bertrand, J., Fletcher, D.F. and Xuereb, C. (2004) PIV Measurements of Flow in an Aerated Tank Stirred by a Down- and an Up-Pumping Axial Flow Impeller. Experimental Thermal and Fluid Science, 28, 447-456.
 Zhu, H., Nienow, A.W., Bujalski, W. and Simmons, M.J.H. (2009) Mixing Studies in a Model Aerated Bioreactor Equipped with and Up- or a Down-Pumping “Elephant Ear” Agitator: Power, Hold-Up and Aerated Flow Field Measurement. Chemical Engineering Research and Design, 87, 307-317.
 Guida, A., Nienow, A.W. and Barigou, M. (2010) The Effect of the Azimuthal Position of the Measurement Plane on the Flow Parameters Determined by PIV within a Stirred Vessel. Chemical Engineering Science, 65, 2454-2463.
 Yu, Z., Cope, R.F., Kar, K.K., Evan, R.C. and Guillaudeu, S.J. (2012) Mixing Performance of the Novel Kar Dynamic Mixer Impeller in Small Laboratory-Scale Systems. Industrial & Engineering Chemistry Research, 51, 15282-15292.
 Ng, K.C. and Ng, E.Y.K. (2013) Laminar Mixing Performances of Baffling, Shaft Eccentricity and Unsteady Mixing in a Cylindrical Vessel. Chemical Engineering Science, 104, 960-974.
 Matos, J., Brito, M.S.C.A., Dias, M.M., Lopes, J.C.B. and Santos, R.J. (2018) Lagrangian Mixing Simulation and Quantification of Scales. Chemical Engineering Science, 192, 199-210.
 Zughbi, H.D. and Rakib, M.A. (2004) Mixing in a Fluid Jet Agitated Tank: Effects of Jet Angle and Elevation and Number of Jets. Chemical Engineering Science, 59, 829-842.
 Cheng, D., Feng, X., Cheng, J. and Yang, C. (2013) Numerical Simulation of Macro-Mixing in Liquid-Liquid Stirred Tanks. Chemical Engineering Science, 101, 272-282.
 Mishra, P. and Ein-Mozaffari, F. (2017) Using Computational Fluid Dynamics to Analyze the Performance of the Maxblend Impeller in Solid-Liquid Mixing Operations. International Journal of Multiphase Flow, 91, 194-207.
 Chergui, N., Lateb, M., Lacroix, E. and Dufresne, L. (2015) CFD Study of Flow Dynamics in a Blade Free Planetary Mixer (BFPM)—A Qualitative Flow Study. Chemical Engineering Research and Design, 102, 100-115.
 Yamagata, T., Matsuzawa, T. and Fujisawa, N. (2017) Numerical Simulation of Flow Structure in Planetary Centrifugal Mixer. Proceedings of the Fourteenth International Conference on Flow Dynamics, Sendai, 1-3 November 2017, 444-445.