2.2. Design and Analysis of Hubless Propeller
Figure 1 represents a flowchart for the design and analysis of the HRDP. Based on the determined parameters
DH,
DH_stk,
DP, NP_rated, and
TP_rated of the HRDP, the basic design of the propeller with the duct is performed by applying NACA66 airfoils and subsequently evaluating the initial performance of the HRDP through the CFD analysis.
Subsequently, for the optimization design of the propeller using the RSM, the thrust and torque of the propeller are defined as objective functions. The pitch (
P), chord length (
c), and camber (
f) of the propeller are selected as factors that directly influence the objective function [
34,
35]. To standardize the factors, pitch was represented as the pitch diameter ratio (
P/
DP), chord length was expressed as the chord length ratio (
c/
DP), and camber ratio (f/
c) was defined as the ratio between blade camber and chord length ratio.
For the 3D shape design of the hubless propeller, the airfoil of the blade was based on the NACA66 profile. The simplicity and standardization of the NACA66 profile facilitates the ease of design and manufacturing, further contributing to its widespread adoption in ship propeller applications [
36].
To represent changes in the chord length ratio (
c) of the propeller based on each radius ratio, the following equation can be used:
where
r is the radial distance from the center hub of the propeller,
R is the total radius of the propeller,
rR is the non-dimensional radius ratio, and
α and
β are the coefficients used to vary the
rR.
Table 3 shows the initial design values for the propeller. The CFD analysis was performed using the flow simulation tool in SolidWorks 2016.
Figure 2 shows how the boundary conditions for the CFD analysis were configured.
The entire flow field is divided into two distinct regions: a rotating domain that encompasses the propeller, stator, and the inner surface of the duct, and a static domain comprised of other components. To establish the dimensions and locations of these various computational domains, a Cartesian coordinate system is employed. The origin, denoted as O, is centered at the blade’s core, with the positive x-axis pointing in the direction of the free flow.
The inlet of the static domain, characterized by all sides spanning 5Dp in length, is situated 3Dp upstream of the HRDP. In contrast, the outlet is positioned 7Dp downstream from the HRDP. Within the computational domain, the fluid is modeled as water. At the inlet, the velocity is set to 0 m/s, establishing an advance ratio J of 0. The outlet is governed by a static pressure condition.
At the rated speed of 3300 rpm, the rated thrust of the initial propeller is 28.6689 N, and the required torque is 0.373902 Nm. This exceeds the rated thrust in the required specifications by approximately 3 N. A higher rated torque would result in an increased output rating for the SPMSM, consequently increasing the volume of the SPMSM. Therefore, it is necessary to satisfy the rated thrust of 25 N while minimizing the rated torque.
The MINITAB 16, a software package for statistical analysis and data analysis, was used for the RSM analysis. The range for
α is set between 0.4 and 0.7, the range for
β is between 0.5 and 0.8, the range for
P/DP is set between 1.4 and 2.6, and the range for f/
c is set between 0 and 0.15.
Figure 3 represents an optimization plot of the propeller using the RSM. The optimized factors have been selected as follows:
α is 0.55,
β is 0.65,
P/
DP is 1.9974, and
f/
c is 0.2250. The confidence level is 97.677%. The RSM results reveal the following values: the thrust is 25.325 N and the required torque is 0.3 Nm. It can be confirmed that all the specifications have been met.
Figure 4a–c illustrate the cross-sections of the initial blade and the optimized blade on the XY plane as a function of the radius ratio R of the propeller.
Figure 4d shows the 3D shape of the optimized propeller. The CFD analysis for the optimized propeller was conducted using the same boundary conditions as those of the initial propeller in
Figure 2.
Compared to the initial propeller, the optimized propeller exhibits a rated thrust of 25.7295 N, which is approximately a 2.93 N reduction, and a rated torque of 0.302528 Nm, which also decreased by about 0.071% as
Table 4 shows. The optimized propeller, with thrust values close to the target of 25 N and reduced required torque, enables energy savings at the rated speed of 3300 rpm by reducing power consumption to approximately 24.66 W.
Figure 5 illustrates the performance comparison of the thrust and required torque between the initial propeller and optimized propeller according to the speed using CFD analysis. Compared to the initial propeller, it can be confirmed that the thrust of the optimized propeller decreased across all speed ranges, leading to a corresponding reduction in torque. The torque-thrust ratio
TTR can be expressed as follows:
where
QP is the required torque of the propeller and
TP is the thrust of the propeller. A higher value of
TTR indicates a greater required torque to generate any given thrust. In all operating speed ranges in
Figure 5d, the
TTR value of the optimized propeller was reduced to 0.0116~0.0122, compared to the
TTR value of the initial propeller, which ranged from 0.013 to 0.1887. This result is possible to save energy because the required torque is minimized while meeting the required thrust specifications.
Figure 5c shows the mechanical output power of the propeller according to the speed. The mechanical output power of the propeller
PP can be expressed as follows:
where
NP is the speed of the propeller. The reduction ratio
PRR in
Figure 5d can be expressed as follows:
where
Pinitial is the mechanical output power of the initial propeller and
Poptimized is the mechanical output power of the optimized propeller. In
Figure 5d, the
PRR values are all positive, ranging from 0.2146 to 0.2442, which means that power consumption is reduced, so replacing the initial propeller with an optimized propeller can increase the operating time of the unmanned underwater drone.
2.3. Design and Analysis of SPMSM
According to the results in
Table 3, the required torque for achieving a rated thrust of 25.7295 N with the optimized propeller should be greater than or equal to 0.3 Nm.
Table 5 represents the required specifications of the SPMSM to drive the optimized propeller. The objective functions chosen for the use of the RSM are as follows: the winding fill factor of 0.35 or less for ease of manufacturing, efficiency of 90% or higher for power consumption minimization, a torque ripple affecting noise and vibration of 3% or less, and an output torque of at least 0.3 Nm.
When satisfying the rated thrust at the rated speed of the propeller, the selection of the rotor diameter of the SPMSM is based on the propeller diameter in
Figure 1. The diameter of the propeller (
DP) is 52 mm, so the inner diameter of the rotor core was chosen as 54 mm, taking into account a 1 mm epoxy injection space for the coupling between the propeller and the rotor core. To prevent water ingress in underwater conditions, it is necessary to apply epoxy coating to both the rotor and stator. Therefore, a 2 mm airgap length was sufficiently chosen to prevent mechanical interference between the rotor and stator.
The torque per unit of the rotor volume (
TRV) for selecting the rotor size of the SPMSM can be expressed as follows [
37]:
where
kw1 is the fundamental harmonic winding factor,
A is the electric loading, and
B is the magnetic loading. The electric field and magnetic field can be expressed as follows:
where
m is the number of phases,
Tph is the number of turns in series per phase,
I is the RMS phase current,
D is the diameter of the airgap,
ϕ1 is the flux density,
p is the number of pole pair, and
Lstk is the stack length.
The HRDP is used in an unmanned underwater drone powered by DC 16 V, and the SPMSM is supplied with the electrical power through an inverter capable of vector control based on the speed PI controller and current PI controller [
38]. Therefore, it is necessary to determine the back electromotive force (BEMF) based on the upper voltage limit of the SPMSM and select the rated current accordingly.
Table 6 represents the results of the basic design of the SPMSM using Equations (5)–(7).
Figure 6a represents the initial shape of the SPMSM based on the basic design.
Figure 7 and
Table 7 present the performance comparison of the initial SPMSM and the optimized SPMSM. The initial SPMSM satisfies the required specifications with 91.33% at a rated current of 10.5 A and an output torque of 0.31995 Nm. However, the winding fill factor is 0.3926 and the torque ripple is 3.4647%, which does not satisfy the required specifications in
Table 4, so it is necessary to satisfy them through optimization design.
To solve these problems, the optimization design of the SPMSM was performed using the RSM according to the procedure in
Figure 1. The tooth width (TW), slot opening (SO), and magnet arc (MA) of the initial SPMSM are selected as factors that directly influence the objective function in
Figure 6a.
Figure 8 represents an optimization plot of the SPMSM using the RSM. The selected factor values are as follows: TW is 4.3437 mm, SO is 4.1204 mm, and MA is 151.1713 degrees, with a confidence level of 97.677%. The fill factor is 0.348, efficiency is 90.3767%, torque ripple is 2.99%, and the output torque is 0.3004 Nm. It can be confirmed that all the specifications have been met.
To validate the reliability of the RSM results, the selected factors were applied to create the shape in
Figure 6b, followed by performing EF analysis. Compared to the initial SPMSM, the optimized SPMSM demonstrates an efficiency of 91.186%, with a slight reduction in output torque to 0.3143 Nm, which still satisfies the required specifications. Furthermore, the winding fill factor is 0.3493, and the torque ripple is 2.9789%, both of which also conform to the required specifications.