1. Introduction
Axial impellers are used, among other things, in processes defined as flow-sensitive operations [
1]. They include the homogenization of miscible liquids, preparation of suspensions of a solid in a liquid, or in heat and mass transfer processes. The term “flow sensitive” suggests that the speed of these processes depends on the liquid flow in the entire volume of the mixing tank. Therefore, knowledge of the characteristics of the liquid flow in the mixer is often the starting point for a quantitative description of these operations.
In our previous work [
2], it was established that the application of an aerodynamic (airfoil) profile for the propeller blade cross-section reduces the mixing power significantly.
Figure 1 presents two of the tested propellers. The use of the NACA4424 airfoil resulted in almost a 25% reduction in mixing power in comparison with the flat cross-section at the same projection area of the blade on a horizontal plane.
However, the reduction in mixing power cannot be the only criterion for assessing the efficiency of new—or only improved—design solutions. The hydrodynamic effects during the operation of the impeller are also important, since they affect the time taken to achieve the required mixing state. In the cited case, to determine the efficiency of the impeller operation, Relationship (1) was used, which determines mixing energy
UE that defines how much energy should be supplied to the system in order to obtain good mixing of the liquid appropriate from a technological point of view.
where
E—energy required to achieve the assumed mixing degree, J;
V—volume of the mixed liquid, m
3;
P—mixing power,
W; 𝜏
m—mixing (homogenization) time, s. Relationship (1) defines the mixing efficiency well, but for its determination, it is necessary to know the mixing time 𝜏
m needed to achieve the required degree of homogeneity. The mixing time can be associated with the circulation time 𝜏
c, which is defined by Equation (2) as the average time needed for each liquid element to make one complete circulation loop in the axial–radial plane.
Information in the literature indicates that three to four circulation loops are required to achieve complete mixing [
3,
4,
5,
6,
7,
8]. This means that the simple Equation (3) can be used to calculate the mixing time 𝜏
m.Information on circulation and homogenization times for individual impellers is most often given on the basis of experimental data [
9,
10]. However, if there are no literature data for a specific impeller, then Relationship (1) is of little use in determining the efficiency of this impeller. While the average circulation time can be calculated, among things, from the relationship [
11,
12], the arbitrarily assumed value of the numerical coefficient in Equation (3) may be a source of error. When the exact mixing time is not known, a better efficiency criterion is Relationship (4), based on the dimensionless pumping number
Flpump and the power number
Po [
13].
where
Vpump—pumping capacity of the impeller, m
3/s;
Fl—flow (pumping) number;
ρ—density, kg/m
3;
D—impeller diameter, m;
T—tank diameter, m;
Po—mixing power number. Pumping capacity
Vpump is the value of a liquid stream thrown out by the impeller from its space, i.e., the area specified by the blades of the impeller during its operation. For axial impellers and a known radial profile of axial velocity
Uz(
r), in the cross-section just below the impeller with radius
RM =
D/2, the values of
Vpump and
Flpump will be expressed by the equations
where
N—rotational frequency of the impeller, s
−1. According to Equation (4), impellers with higher values of
Ep are more effective, because for a given mixing power they generate a larger liquid stream.
For the mixing process, the power number is defined by the relationship
During the rotation of the impeller, its blades are subjected to aerodynamic force resulting from two sources:
Normal force caused by pressure on the body surface;
Shear force due to fluid viscosity, often referred to as skin friction.
The resultant aerodynamic force acting on the body is equal to the forces due to pressure and the sum of shear forces acting over the entire body surface. This aerodynamic force is usually broken down into two components (both components pass through the center of pressure):
Drag force, i.e., the component of the force parallel to the direction of rotation of the impeller blade (depending on the resistance coefficient Cd);
Lifting force, i.e., the component of the force perpendicular to the direction of rotation of the impeller blade (depending on the coefficient of the lifting force Cl).
It can be shown that the mixing power number
Po depends on the blade drag coefficient
Cd [
14], assuming that
,
,
(
F—force, N).
That is, the drag coefficient of the impeller Cd and the power number Po are directly proportional. The above description is simplified, because when the impeller rotates, induced drag is also created. It is generated by the circulation of fluid around the rotor blade, which creates lift. For a blade pitch of 45° and flow completely stalled at the airfoil pressure surface, shape drag is likely to have a greater effect on torque and, consequently, on the mixing power (large projection of airfoil area on the r–z plane). During the experimental determination of the torque, its value is influenced by both the shape drag and the induced drag. However, for an accurate theoretical description of the phenomenon, induced resistance should also be taken into account.
As is known, the power number is a measure of the efficiency of converting the mechanical energy supplied to the mixing tank into kinetic energy of the liquid [
3,
4]. The smaller the value of the power number, the greater the conversion efficiency, because less energy input needs to be provided to obtain the same kinetic energy of the liquid. As can be seen from Equation (8), impellers with more streamlined blade shapes, i.e., with lower values of the drag coefficient
Cd, are characterized by lower power demand. In addition, the application of an aerodynamic blade cross-section makes sense for axial impellers, because impellers of this type require less mixing power than radial impellers.
Axial impellers are mainly used for the preparation of suspensions. The energy supplied to the liquid–solid system generates turbulence on a microscale, i.e., on the Taylor and Kolmogorov scale, as well as eddies on a scale comparable to the scale of disturbing elements, i.e., the Brodkey scale [
15,
16]. The latter lift solid particles from the bottom and disperse them in the entire volume of the liquid. For obvious reasons, the average velocity of the liquid stream should be higher than the falling velocity of the suspension particles. In a study on propellers with profiled blades [
2], no deterioration in hydrodynamic conditions was found compared to the classic propellers; however, it should be noted that the system pumping the liquid to the top of the tank was tested. The direction of flow of the main stream in the mixing tank does not affect the mixing time, while in the preparation of suspensions, the flow of liquid from the area of the mixer towards the bottom is more often used. Therefore, it is important to check how the use of the aerodynamic profile in the construction of the impeller blade affects the hydrodynamics of the liquid in the mixer, i.e., the velocity average components for other types of impellers and a different direction of pumping the liquid from the area of the impeller.
At the same time, based on the knowledge of the size of the liquid main stream and the structure of its turbulence, it is possible to determine individual scales of turbulence [
17] and the rate of energy dissipation [
18,
19,
20,
21,
22,
23,
24]. Knowing the rate of energy dissipation can be useful, among other things, for determining the value of the mass transfer coefficient in liquid–gas systems [
10,
25,
26,
27] or liquid–liquid systems [
28,
29,
30,
31]. Because PBT mixers are often used in industry, they have therefore been the subject of various studies [
32,
33,
34,
35,
36,
37,
38]. However, most of the research on rotors with airfoil cross-section blades concerns marine propellers and wind turbines [
39,
40,
41,
42,
43,
44]. A more extensive discussion of the literature data on PIV measurement methods and the obtained results can be found in
Section 2.1 and
Section 3.2.
The aim of this study is to investigate whether the use of impellers with blades of aerodynamic sections, most likely characterized by lower mixing power, will make it possible to achieve higher mixing efficiency from a hydrodynamic point of view.
2. Materials and Methods
This paper presents research on PBT impellers with six inclined blades with an inclination (attack) angle of
α = 45° (
Figure 2) that pump the liquid towards the bottom. Individual impellers differed with respect to the cross-sections of their blades. Blades with a flat section were used (standard impeller PBT45-6: hub diameter
dh = 0.2·
D, blade width
b = 0.2·
D,
Figure 2) as well as blades with profiled cross-sections characterized by high values of lifting force, developed by the National Advisory Committee for Aeronautics (NACA): symmetrical NACA0021 and asymmetrical NACA4424 in two orientations. The individual cross-sections of the blades are shown in
Figure 3. For blades with airfoils, blade width is equal to chord length (
b =
C).
As shown in
Figure 3, the asymmetrical profile NACA4424 can be positioned in two ways relative to the direction of movement, as shown in
Figure 4, where the arrows show the direction of the blade movement. The forces acting on the blade set as in
Figure 4b,c will be identical, but for the setting in
Figure 4c (negative angle of attack), the drag coefficient values should be more accessible, at least for small angles of attack. Therefore, in subsequent parts of the paper, the flow of the impeller blade set as shown in
Figure 4c is analyzed.
Despite the fact that the profiles that were used were developed many years ago, they are still being studied [
45]. The impellers shown in
Figure 2 with the cross-sections of the blades as shown in
Figure 3 were made using the 3D printing technique.
2.1. Measurements of the Liquid Velocity in the Mixer
Measurements of velocity distributions in the vertical r–z plane were carried out in a transparent vessel with a diameter T = 292 mm with a flat bottom and filled with distilled water to a height H = 300 mm (H ≅ T). Inside the tank, there were four standard baffles with a width of B = 0.1⋅T placed at equal intervals. A PBT45-6 turbine-blade impeller with a diameter of D = 100 mm suspended at a height of Z = 100 mm above the bottom of the tank was tested. As flow markers hollow glass spheres with an average diameter of 10 µm were used. To reduce the possibility of optical distortions, the cylindrical vessel was placed in an additional rectangular glass container also filled with water. Based on the analysis of the results of preliminary measurements, the rotational frequency of the stirrer was limited to two limit values: N = 90 min−1 (Frm = 0.023, Rem = 15,000) and N = 240 min−1 (Frm = 0.163, Rem = 40,000). The direction of rotation of the impeller was selected so that the liquid flowing through the impeller zone was directed down the mixer as is appropriate when mixing suspension under industrial conditions.
The measurements were performed using the
LaVision PIV measurement system with a two-pulse laser with a maximum power of 135 mW and an
ImagePro camera with a resolution of 2048 px × 2048 px with a Nikkor 1.8/50 lens. The lens aperture was closed to the value ensuring the maximum resolution (i.e., the aperture value was 5.6 [
46]).
In the measurements, the plane of the light knife was placed about 2° before the baffle. The diagram of the light knife position is shown in
Figure 5. The measurement field was approximately 300 mm × 300 mm. The frequency of the laser double flashes (pulses) was 2.7 Hz. Thus, the images were recorded for different positions of the blades in relation to the plane of the light knife. For the double laser flash, the intervals were
Dt = 4200 µs for
N = 90 min
−1 and Δ
τ = 1900 µs for
N = 240 min
−1 and were determined for the velocity of 0.65·
Utip and the shift of the tracer particles 10 px. The DaVis 7.2 program uses sub-pixel interpolation, and under these conditions the particle displacement is estimated up to 0.1 pixels. This means that the expected velocity accuracy is σ
u = 0.1/(
M·Δ
τ), where
M is the magnification factor (
M = 6.83 px/mm). The calculated estimated accuracy of the instantaneous velocity measurement is 0.0032 m/s and 0.0077 m/s respectively. For each rotational frequency and blade cross-section, 100 image pairs were taken to average the results [
47]. The proprietary program
DaVis 7.2 was used for data processing. Two-pass data processing was used with the final size of the analyzed field 32 px × 32 px (i.e., approximately 4.7 mm × 4.7 mm) without overlapping. To make comparisons easier, measured liquid velocities were divided by the tangential velocity of the blade tip to obtain a dimensionless value.
2.2. Measurements of Mixing Power
The measurements of the mixing power were carried out in a flat-bottomed tank with a diameter of
T = 400 mm, equipped with four standard baffles (
B = 0.1·
T). The tank was filled with water (
t = 20°) to the height of
H =
T = 400 mm. For the tests impellers with a diameter of
D = 133 mm (
T/
D = 3) were used. This system was geometrically similar to the system with a diameter of
D = 292 mm shown in
Figure 5, in which measurements of velocity distributions were made.
Torque measurements were made with the IKA EURO-ST P CV meter. Measurement data were recorded using the Labworldsoft 4.6 program with the sampling time Δt = 2 s. The values of the torque M were recorded for the set rotational frequency of the impeller N. The rotational frequency was changed within the range from N = 50 min−1 to N = 400 min−1, which corresponds to the values of Reynolds number for the mixing process from 14,750 to 118,000.
The measured values of the rotational frequency
N and the torque
M were used to calculate the mixing power number according to the relationship