1. Introduction
The list of variables used in the paper is presented in
Table 1.
The field experiments [
1] of a self-propelled remotely operated river push barge model in 1:20 geometrical scale equipped with the bow steering system consisting of two rotating cylinders (rotors) showed the significant enhancement of her motion control due to the rotor-generated hydrodynamic force.
The push barge model during the field experiments is presented in
Figure 1.
The model tests followed by the numerical simulation confirmed that the hydrodynamic steering force generated by the bow system was dependent on rotors dimensions, their rotational speed, water inflow velocity, and shape of the barge bow [
2,
3].
The future application of the bow steering system for autonomous river barges is related to the development of an algorithm for the rotor angular velocity control based on the reliable model of the generated hydrodynamic force. The implementation of this algorithm in the ship-handling model along with the automation of main on-board systems and their control integration in a single communication channel is the first step of transition from level zero to level one of autonomy, e.g., remote vessel controlled by a human operator [
4].
The program of the tests presented in the paper is adapted to the rotor models in a larger geometrical scale than the scale used in the field experiments, e.g., 1:10. The velocity scale according to the Froude principle of similarity is
. The assumed rotational velocity is in the range of 0–550 rpm and the Reynolds number related to the flow is in the range 2 × 10
4–1.23 × 10
5 The examined cylinder aspect ratios are chosen as 0.5, 1, and 1.5. The round endplates are installed on the cylinder to reduce the strength of the expected tip vortices affecting the hydrodynamic force [
5].
This paper presents a method developed for the measurements of the hydrodynamic force generated on a rotating cylinder under the water current and the experimental test setup constructed for the model tests. The verification and validation methods of the hydrodynamic force computed using the CFD (Computation Fluid Dynamics) FLOW-3D program are proposed.
The experimental test setup presented in the paper is a unique construction consisting of a bearing mounted platform with the rotor drive and sensors measuring the components of the hydrodynamic force, installed in the flume tank equipped with an advanced wave generator and system for measurements of the flow field parameters.
The available experimental results for the prediction of the rotor-generated flow field are mainly based on the experiments carried out in wind tunnels or carried out in the water without the free surface effect. The published experimental and numerical data related to both the rotor aerodynamics and hydrodynamics present a different combination of the significant influencing parameters.
The endplate diameter to cylinder diameter ratio D
e/D and the cylinder aspect ratio (length to diameter ratio) L/D were studied by Thouault et al. [
5]. Their general conclusions from CFD simulations are in line with the results presented in this study. The flow separates on each endplate edge and the impact of tip vortices on the flow field generation is significant at high rotational rates, especially for low-aspect-ratio cylinders with small endplates. The obtained lift coefficient (C
L) was equal to 4.71 for the cylinder with L/D = 1.3, D
e/D = 2, α = 2, and Re = 7.1 × 10
4. The free surface effect was not included in that study.
The unsteady flow past a rotating circular cylinder at flow related Reynolds numbers Re = 1 × 10
3 and Re = 1 × 10
4, and cylinder rotation rate ratio α (tangential velocity to free steam velocity ratio) between 0.5 and 3, was investigated both theoretically and experimentally by Badr et al. [
6]. This experimental investigation based on flow visualization showed the difference between numerical and experimental results especially at the highest rotation rates, when the three-dimensional effects become significant.
Karabelas et al. [
7] presented lift, drag, and hydrodynamic moment coefficients obtained from numerical modeling at Re = 5 × 10
5, Re = 1 × 10
6, and Re = 5 × 10
6.
Chen and Rheem [
8] presented the experiments carried out for different aspect ratios of cylinders for Re = 1 × 10
5 and α in the range between 0 and 8.
Zhou et al. (2015) [
9] presented the experimental study on the water flow past a non-rotating circular cylinder at Re = 7 × 10
3–8 × 10
4, considering the different groove patterns on the cylinder surface.
Tokumaru & Dimotakis [
10] described the influence of forced rotary oscillations added to the steady rotation on the wake separation and the mean lift force for the cylinder of the large aspect ratio equal to 18.7. They presented the lift coefficient values obtained from measurements and compared with other published results at Re = 3.8 × 10
3–1.2 × 10
5 and cylinder aspect ratios 4.7–18.8. For α = 1, the lift coefficient values were between 0.8 and 1.2, and for α = 2 they were between 3 and 4, which can be compared with the results of measurements presented in this paper for the cylinder with aspect ratio 1.5 when the free surface effect at small rotational speed can be neglected.
Wong et al. [
11] presented the results of experimental investigations of flow-induced rotor vibrations at flow dependent Reynolds numbers in the range of 1.1 × 10
3–6.3 × 10
3 and noticed the vortex shedding cease and decrease of the oscillation amplitude at α = 2. The observed phenomena were significantly different from the comparative CFD predictions and other experimental results at lower Reynolds numbers described by Bourguet & Jacono [
12].
Carstensen et al. [
13] presented the experimental test setup for the measurements of the hydrodynamic force on the rotating spar buoy model. The results included lift and drag coefficients at Re in the range of 3.2 × 10
4 to 9.9 × 10
4 and α between 0–4. In these experiments, the tested cylinder was divided into three parts. Only the middle part with diameter 0.160 m and aspect ratio 2.5 was rotating. The free surface effect was limited due to the not rotating upper part of the cylinder. The bottom tip vortices were limited by not rotating the lower part of the cylinder. The results included lift force coefficients at Re in the range of 1 × 10
4 to 6 × 10
4 and α from 0 to 5. The lift coefficient for α = 2 was equal to 4.
In the experiment presented by Wang et al. [
14], the tested cylinder was installed under the plate situated on the water surface, eliminating the free surface effect and tip vortices on the top of the cylinder. For the cylinder of 0.1 m diameter and aspect ratio 1, at Reynolds numbers in the range of from Re = 2 × 10
5 to 1 × 10
6 and α = 2, the lift coefficient was equal to 3.5.
Mobini & Niazi [
15] presented the lift coefficient at different Reynolds numbers computed using the Large Eddy Simulation method. For Re = 6 × 10
4 and α = 1, C
L was equal to 0.75, compared to C
L = 1.15 presented by Karabelas [
7]. For Re = 1 × 10
4 and α = 1.5, C
L was equal to 1.95, and for α = 2, it was equal to 2.75.
The research presented in this paper demonstrates the sensors used in the experiment and accuracy of measurements determined by comparing the obtained results with available published data. CFD simulation of rotor generated hydrodynamic force was carried out and validated for the same ranges of parameters as the experiment.
2. Materials and Methods
2.1. Hydrodynamic Lift Force Measurements
The measurements presented in the paper were carried out in the wave flume of the hydraulic laboratory of the Institute of Hydro-Engineering of Polish Academy of Sciences in Gdansk. The wave flume was 64.0 m long, 0.6 m wide, and 1.40 m high with side walls made of 0.018 m thick glass sheets allowing for visualization of the measurements. The scheme of the measurement area is presented in
Figure 2.
The flume is equipped with the FLYGT pp4640 power turbines, which, depending on the water depth, can generate water flows with velocity from 0.1 m/s to 1.0 m/s in both directions [
16,
17]. The water current generation is electronically controlled and ensures the stability of the mean flow velocity with deviations no greater than 2%. The flow is regulated by the honeycomb, which is a flow straightening structure used to avoid large eddies that can be generated by pump and to minimize the turbulence effects of the flow [
18]. The honeycomb is composed of pipes of 0.2 m length and 0.035 m diameter each, vertically positioned along the flow axis. The water circulation is provided by return pipe with the diameter of 0.5 m located under the wave flume.
The force measuring part of the experimental test setup is placed on the supporting frame made of aluminum flat bars and aluminum profiles (
Figure 3).
The construction of the frame allows for changing the immersion of the tested cylinder with adjusting screws. The main elements of the force measuring experimental test setup are BLDC (BrushLess Direct-Current) drive motors and replaceable cylinders.
The BG75X50SI drive is integrated with a planetary gear with a 1:4 ratio to ensure the appropriate system dynamics at rotational speeds of the cylinder in the range up to 1000 rpm and torque up to 5 Nm.
There are thermal and overload protections implemented in the control system to protect the drive in the event of uncontrolled dynamic states. The drive is controlled by the operator console equipped with a seven-segment display of the rotational speed, potentiometer regulating the drive voltage level, and a two-position switch enabling the change of the stator rotation direction.
The tested cylindrical rotors were made of PA6 polyamide, which is easy for technological processing, with high mechanical strength, low moisture absorption, low roughness, and optimal rotating mass. Each rotor was equipped with rounded end plates.
The entire drive system along with the measuring devices has been integrated with a bearing-mounted platform free to move in two directions along x and y axes.
In different solutions of force measuring devices, the rotor drive was mounted on a three-component force measuring balance [
16], or two force gauges rigidly fixed to the drive shaft at the top and bottom of the rotating part of the cylinder [
13]. The free movement of the measuring section, relative to the drive shaft, allows in this solution to measure its displacement using the strain gauges on the force gauges [
13].
The free surface effect in the known experimental setups is eliminated by the plate situated on the water surface [
16] or the non-rotating upper part of the cylinder crossing the free surface [
16].
The scheme of the measurement area used in the presented experiment is presented in
Figure 4. The rotating cylinder was located 25 m upstream from the water pump and its vertical position was in the center of the water column. The still water depth was equal to h = 0.262 m. The location of the rotor, rotor drive, and supporting frame in a wave flume and its dimensions are presented in
Figure 4.
2.2. Measuring Sensors and Procedure
The main objective of experiments was to measure the lift force FY generated on the rotating cylinder in the presence of the stationary water current.
The water current velocity was recorded by two Acoustic Doppler Velocimeter probes located 5 m upstream and 5 m downstream of the rotating cylinder (
Figure 2). The Nortek instruments measure the velocity of water by utilizing the Doppler Effect [
19]. The ADV Vectrino, 3D Lab Probe measures the acoustic signal reflection from the seeding particles suspended in water with an operating acoustic frequency equal to 10 MHz, giving the water particles velocities with an error of ±0.001 m/s [
19,
20]. It requires a great amount of seeding material to provide its proper concentration in the whole volume of water.
The aim of those measurements was to confirm the stationary conditions of flow. Additionally, the water current velocities were determined indirectly by the flow rate of the pump.
The forces generated by the rotor were recorded using the KMM20 force sensors equipped with threaded fasteners capable to measure both the tension and compression. Their operation is based on a strain gauge bridge. The small size sensors with a wide measuring range of 200 N nominal rated forces and maintaining the accuracy class of 0.5 have proven themselves very well in this application. An additional advantage of the sensors was the built-in signal amplifiers, which significantly simplified data registration. The voltage readings of the transducer were incorporated into the DHI Filter Cabinet Type 155/IF controlled by the National Instruments DAQmx with the NI 6229 device and converted to forces. Linearity tolerance of the transducers was 0.25% of the reading.
The investigated cylinder diameters were as follows: D1 = 0.05 m, D2 = 0.10 m, and D3 = 0.15 m. The endplates diameters to the cylinders diameters ratios De/D for all cylinders were equal to 1.1. The roughness coefficients of cylinders are as follows: Ra: 2.5 μm, 1.6 μm, and 4 μm respectively. The roughness of the cylinders was measured using a HOMMEL-ETAMIC W5 device equipped with a head for measurements on cylindrical surfaces. The measurement was performed five times in different places of each rotor. The mean value of the measurements was determined and used as the roughness coefficient in simulations.
For each cylinder, five current velocities were tested: 0.4 m/s, 0.5 m/s, 0.6 m/s, 0.73 m/s, and 0.82 m/s. The rotational speed of the cylinders was in the range 0–550 rpm with the step of 50 rpm (
Table 2).
The cylinder rotated counterclockwise. The recording duration of a single measurement was 1 min for each rotational speed and 5-min gap between the subsequent rotor rotational speeds was made.
3. Results of Measurements
The measured lift force F
Y in the adopted coordinate system (
Figure 4) is presented in
Figure 5. For low current velocities, the increase of absolute force values is constant through the whole rpm range, especially for the lowest considered cylinder diameter 0.05 m. The force is higher for higher cylinder diameters.
This increase of force for rotational speed up to 200–250 rpm follows the vortex shedding process, which creates the von Karman street in the wake of the cylinder [
8]. Beyond certain rpm, this process deteriorates limiting further increase of the lift force. This effect can be seen in
Figure 5 for the highest considered rotational speeds when the absolute value of the force seems to achieve its maximum. That means a further increase of the rpm would not lead to an increase of the force. The detailed investigation of this problem was carried out by Karabelas et al. [
7].
The measurement error for the force was 1.2 N. The error related to non-linearity (% F.S.) was 0.25% of 200 N, hysteresis (% F.S.) was 0.25% of 200N, and the creep error 30 min (% F.S.) was 0.1% of 200 N.
The lift force coefficient C
L (1) related to α (2) was analyzed with respect to Re (3) for all the tested flow conditions and cylinder diameters used in the wave flume experiment.
where: C
L is the lift coefficient, F
Y is the lift force obtained from the experiments, L is the height of the cylinder, and D is the cylinder diameter.
where: α is the cylinder rotation rate, r is the angular velocity of the cylinder [rad/s], D is the cylinder diameter [m], and V is the free stream velocity [m/s].
where: V is the free stream velocity [m/s], D is the cylinder diameter [m], and ν is the kinematic viscosity of water: ν = 10
−6 m
2/s.
Reynolds numbers for the tested flow velocities and cylinder diameters are presented in
Table 3.
The lift coefficients for α in the range 0–2 for the cylinders with D
1 and D
2 diameters and in the range 0–4 for the cylinder with D
3 diameter are presented in
Figure 6.
The comparative results are available for α no greater than 2 when the free surface effect at small rotational speeds can be neglected. The results obtained by different authors [
7,
10,
15] are presented in
Table 4.
In the present study, the measured values of CL are in the ranges of values presented in the literature.