1. Introduction
The propulsion test of ship models on actual seas is a test performed using a large-scale model in a natural water environment. The test site is selected in an offshore area where the wind, waves, and currents are similar to the navigation conditions of a real ship. This ensures that the performance of the ship propulsion system can be predicted more accurately. Since the test object is a large-scale model, such tests can also effectively alleviate the scale effect in the towing-tank tests, so as to improve the accuracy of extrapolation calculation. It also significantly saves manpower, materials, and financial resources compared to tests using a full-scale ship, so that repeated tests can be carried out. In addition, it is difficult to modify the full real ship once it has been built and launched. Hence, the method has broad development prospects. However, as this kind of test technology is still in the exploratory stage and there is insufficient comparative research data, it is necessary to evaluate the uncertainty of the test to determine the reliability of its scheme and the validity of the measurement data.
Uncertainty indicates the degree of trust in the measurement results. For the ship model test, the International Towing Tank Conference (ITTC) initially followed the American Institute of Aeronautics and Astronautics (AIAA) standard (1995) to analyze the uncertainty from the viewpoint of bias and precision errors [
1,
2]. Stern et al. provided a summary of the AIAA standard for the experimental uncertainty assessment methodology [
3]. They assert that the benefits of uncertainty analysis far outweigh the time it takes to evaluate measurement results. Subsequently, concomitant with the widespread application of the “Guide to the Expression of Uncertainty in Measurement (GUM)” around the world [
4], the uncertainty Analysis Committee of ITTC also voted in to recommend the use of this method to evaluate relevant experiments [
5]. In order to compare the characteristics of the old and new methods, Delen et al. used two ITTC (2002 and 2014) procedures based on different standards to analyze the uncertainty of ship model resistance tests in towing tanks [
6]. With the gradual improvement of the uncertainty analysis method of resistance tests, Nikolov et al. studied the uncertainty of tow-tank testing of high-speed planning craft, in order to generate more data to incorporate into the evaluation process [
7].
The test sites used in the above studies were all in towing tanks, and most of them evaluated resistance tests. With the diversification of the experimental environment and the complexity of the computational model, the traditional rules based on the Taylor series method (TSM) cannot fully meet the present evaluation needs [
8]. Recently, several scholars have used the Monte Carlo method (MCM) to study the uncertainty of ship tests, especially for trials on actual seas. For example, Insel analyzed the navigation test of a group of sister ships at sea by utilizing the ITTC standard speed/powering trial analysis procedure with Monte Carlo simulations, with particular focus on the influence of the Beaufort scale and other uncertainty factors [
9]. Kamal introduced the MCM into the extrapolation prediction of dynamic ship performances [
10]. By referring to the ITTC1978 procedure, the uncertainty caused by the error source in the resistance and propulsion test was studied. Samples were taken from the existing sea trial and towing tank test, and the uncertainty of the system was calculated. Woodward analyzed the uncertainty of the maneuverability test of the ship model, showed the transfer process of the uncertainty component by the MCM, and determined the error source [
11,
12]. Subsequently, Vrijdag compared the advantages and disadvantages of the linear method and the MCM in ship performance prediction and analyzed the factors that significantly influence ship speed and bollard pull predictions [
13]. Gavrilin analyzed the maneuvering sea trial results of the ship “Gunnerus” and proposed a method to evaluate the uncertainty of test results caused by environmental impact using a simulation model based on the MCM [
14]. Aldous detailed and categorized the relevant sources of uncertainty in performance measurements and proposed a method to quantify the overall uncertainty on a ship performance indicator based on the MCM [
15]. A sensitivity analysis was also conducted on the sources of uncertainty to highlight the relative importance of each source. Seo studied the speed and power performance results of real ships using the MCM [
16]. Based on the correction method of real ship data published by the ISO, the Monte Carlo results were approximated to the normal distribution, and the uncertain sources under actual sea conditions were analyzed.
This present study was conducted with the objective of using prototype ships to perform the tests. However, for ship model propulsion tests on real seas, there is no relevant uncertainty evaluation data. Therefore, in this study, according to the methodology recommended by the ISO, the uncertainty of large-scale ship model propulsion tests on actual seas was evaluated and analyzed based on the Monte Carlo numerical simulation method. The uncertainty of the environment and propulsion efficiency were also simulated and calculated using an adaptive program.
The aim is to quantify the uncertainty caused by various factors in the test through the analysis of a ship model propulsion test on actual seas, so as to obtain the influence of different error sources on the measurement accuracy. The main part of this study is to analyze the interference of the test environment and discuss the influence of wind wave and current. At the same time, the uncertainty variation trend caused by wind and waves of different degrees was simulated through the program, which provided data reference for similar experiments and clarified the direction of further work.
3. Results
Similar to the towing tank test, uncertainty analysis of the propulsion test on an actual sea requires consideration of the geometric structure uncertainty of the ship model, installation of equipment, and measurement of related quantities, where the Taylor series method can be used in this part. In addition, because the experiment was performed under actual sea conditions, it is necessary to consider the interference of wind, wave, and current in the environment. The fluctuation of these factors in the entire test process suggests non-linear characteristics, and the related correction formula is complex. Thus, in this study, the MCM was used to evaluate the uncertainty.
3.1. Uncertainty of Ship Model
The uncertainty of the ship model conditions can be divided into the following parts [
18]:
Before each sea test, the floating state of the ship model requires adjustments in the wharf, and the relative error of displacement is generally required to be less than 2%. In this test, the estimated error of the displacement of the ship model was 1000 kg, and the relative deviation was about 1%.
The trim condition of the ship model is an important parameter for the test, but it does not consider the final range of uncertainty evaluation. Generally, the trim is less than 0.1% of the length between the vertical lines, and the error of the bow draught within a range of ± 0.1 m can be considered to meet the test requirements.
During ship navigation, the hull roughness or growth of marine life will significantly increase the resistance of the ship. Presently, there is no accurate correction method in this respect. Hence, the hull and propeller require a careful check before the test. The ship model in this experiment spent a short time at sea, and it was observed that the surface of the propeller was relatively smooth when the propeller was dismantled in the middle of the test. Thus, the influence of roughness was ignored.
3.2. Uncertainty of Shafting Measuring Instrument
In this test, the torque and thrust of the shafting were measured by a special self-propulsion instrument in which the torque measuring device was a T40B (HBM) digital torque sensor with a measuring range of 500 Nm, and the measuring accuracy was 0.25 Nm. The measuring accuracy of the force sensor used in the thrust was 2 N. The shafting motor was a Siemens motor, which offers a high precision and a rotational speed precision of 0.5 rpm. The standard uncertainty of the parameters in the shafting can thus be expressed as follows, the coverage factor k corresponding to the instrument is 3:
3.3. Uncertainty Caused by Environment
Based on the results collected from the instrument, the environmental parameters in the two tests are as shown in
Table 2. In fact, the relative wind speed is different with different ship speed. However, due to the low accuracy of the anemometer used in this test, which could also be due to variability of the wind, the difference between the measured data is small. So we use the average relative wind speed over the entire time period of each group as the final speed. Also because of this setting, we can compare the uncertainty of different speeds at the same relative wind speed.
3.3.1. Uncertainty of Wind Resistance
The resistance increase due to relative wind is calculated using:
where
RAA is the wind resistance;
ρA is the mass density of air,
ρA = 1.226 kg/m
3; and
AXV is the projection area of the hull above the waterline. Using the three-dimensional graph of the model, we estimated
AXV = 8.7 m
2.
VWRef and
VG represent the relative wind speed and the ship speed, respectively;
CAA is the wind resistance coefficient, the relative wind direction is represented by
ψWRef, and the ship is sailing upwind when the degree is zero. Based on the suggestion, in this study, the data of a portable bulk carrier was selected to calculate the wind resistance coefficient [
18,
19].
Through the analysis of the data obtained by the wind speed and direction meter and by referring to the relevant literature, it was assumed that the wind speed in the parameters is in agreement with the normal distribution, accuracy is 0.5 m/s, wind direction obeys the uniform distribution, and precision is 10°. The input distribution of the two tests is presented in
Table 3.
The wind resistance of the ship model under various working conditions was simulated using a Monte Carlo program, and the expanded uncertainties of the wind resistance of the two groups of propulsion tests are presented in
Table 4. Therefore, the relative expanded uncertainty of wind resistance in the propulsion test is extremely large.
Notably, in order to facilitate calculation of the total uncertainty and the statistics of the data, the distribution function of the wind resistance was approximated to a normal distribution, which is not completely equal to the actual simulation results. For example, the wind resistance distribution of the ship model when the propeller speed is 250 r/min in the second test is shown in
Figure 3.
It can be seen from the figure that the probability density function of wind resistance shows a slight asymmetry. In the 95% probability inclusion interval, the distance between the two endpoints and the mean is close. Comparing the frequency distribution histogram with the normal distribution curve using the same mean and standard deviation, it is observed that the figure deviation is small and the inclusion intervals are similar. Therefore, it can be considered that this approximation is reasonable, and the corresponding wind resistance can be expressed as follows:
3.3.2. Uncertainty of Resistance Increase Due to Wave
In this propulsion test, the sea conditions in the sailing sea area were relatively mild. Based on the real-time measurement results of the wave altimeter, the waves encountered by the ship were high-frequency short waves. In this case, the drag increase in hull motion caused by waves can be ignored, and the additional resistance is mainly caused by the wave reflection of the hull on the waterline. Therefore, in this study, the STWAVE-1 method recommended by the ISO was used to estimate the wave resistance, as shown in the following formula:
Among them, ρS is the sea water density = 1026 kg/m3; g = 9.80 m/s2; HW1/3 represents the significant wave height of 0.14 m; B is the width of the ship model, which is 4.04 m; and LBWL is the distance from the bow waterline to 95% of the width of the ship, with a measured value of 2.51 m. The vertical acceleration at the bow measured by the sensor was 0.18 m/s2 < 0.05 g = 0.49 m/s2, which means the trim and pitch conform to the conditions of the formula. In the wave drag increase correction, only the area of the wave direction angle within the heading angle ± 45 is considered, and if the wave direction is outside this area, the wave resistance increase is not considered. In this experiment, only the wave direction angle of the second group test meets the requirement of correction, thus the wave resistance increase of the first group test was zero.
Based on the measurement accuracy of the wave altimeter, the distribution of
HW1/3 in the input is defined as N (0.14, 0.02). Other parameters including the sea water density and the uncertainty of the ship model value were very small, thus negligible here. The frequency distribution histogram of the simulation is shown in
Figure 4.
Comparing the histogram with the normal distribution curve, it can be observed that the peak value of the histogram shifts to the left, which may be related to the non-negativity of wave resistance. Considering that their figure deviation is small and the inclusion interval is similar, the distribution function of wave resistance can be approximated as a normal distribution, which can be expressed as follows in 95% confidence probability:
3.3.3. Uncertainty of Water
The temperature and density of the water will affect its viscosity, thus affecting the resistance of the ship model. The actual seawater temperature of the test differs from the standard temperature, the prediction and calculation in the propulsion test is usually based on the water temperature of 15 °C, and the viscosity corresponding to 1026 kg/m
3 is 1.1892 × 10
−6 m
2/s. The sea water temperature in this test was 16.8 °C, thus it requires correction. Since the effect of roughness is not considered, the effect of sea water temperature and density on the resistance can be expressed by the following formula:
where
CF is the friction resistance coefficient of the actual water temperature and water density;
CF0 and
CT0 are the friction resistance coefficient and total resistance coefficient of the reference water temperature and water density, respectively;
RF is the friction resistance of actual water temperature and water density;
RF0 and
RT0 are the friction resistance and total resistance of the ship model under reference water temperature and water density, which are calculated from the data obtained from the resistance test;
S is wet surface area;
VS is the speed of ship model; and
ρ and
ρ0 are the sea water density under actual and reference water temperatures and salt content, respectively. The parameters in the formula were calculated using the ITTC1978 formula, as shown in
Table 5.
Among them, the input water temperature is normally distributed, and the interpolation function of water density and viscosity with respect to temperature change is shown in the following formula:
The tests under various working conditions were simulated respectively, and the expanded uncertainty of resistance affected by water temperature and density of the two groups of tests is shown in
Table 6.
It can be seen from the table that the uncertainty of resistance caused by sea water temperature and density increases with the rotational speed, but the relative expansion uncertainty of resistance itself is maintained at roughly 10.8%, which is caused by the measurement accuracy of the thermometer. The value of this uncertainty is large, and it is difficult to reduce from the instrument’s measurement. However, considering that the resistance caused by seawater is small, it has little effect on the total resistance of the ship model.
Figure 5 is the frequency distribution histogram of a working condition, from which we can see that the resistance distribution frequency diagram caused by water temperature and density is consistent with the normal distribution curve.
Furthermore, the distance between the two endpoints and the mean is almost equal in the 95% probability range. Therefore, the related resistance can be expressed as in the following formula:
3.3.4. Influence of Shallow Water Effect on Speed
The influence of the shallow water effect can be expressed by the following formula (Lackenby shallow-water correction):
where
AM is the waterline area of the ship model,
h is the water depth of the test area,
VS is the speed of the ship model, and
ΔV is the influence of the shallow water effect on the speed. The site of this experiment was an area with a large water depth selected through chart selection as there was no professional instrument. Hence, the influence of shallow water effect was not considered. The uncertainty of shallow water effect is mainly determined by the measurement uncertainty of velocity and water depth.
3.3.5. Influence of Current on Speed
This propulsion test adopted the “MOM” test method recommended by ISO, which assumes that the velocity change of the current flow is parabolic with time. In the test, we carried out two double runs tests under the same rotational speed, that is, four times for each working condition. The speed of the ship model can be expressed as follows:
where
VS represents the relative speed of the ship model to the water and
VG1, VG2, VG3, and
VG4 are the four ground velocities at each power measured in the test.
The second group of test data was selected to modify the speed of the ship model. All the velocities in this experiment were measured by the Global Positioning System (GPS). Based on its measurement accuracy, the uncertainty of current on velocity was calculated by Taylor formula, as shown in
Table 7.
It can be seen that as only the systematic error of speed was considered, the expanded uncertainty at each speed is only related to the measurement accuracy of the GPS. Therefore, with increasing ship model speed, the influence of ocean current on speed uncertainty will decrease.
3.4. Uncertainty of Propulsion Test
One purpose of this self-propulsion test was to judge the propulsion performance of the ship. To express the overall hydrodynamic performance of the propulsion system in ship propulsion, the effective power
PE of the ship was compared with the received power
PDB of the propeller, which is called propulsion efficiency η
D:
RX =
R + Δ
R, where resistance
RX is obtained by adding the resistance value measured by the large-scale ship model resistance test performed in the same period to the correction value in the propulsion test. The equation Δ
R =
RAA +
RAWL +
RAS expresses the correction and includes changes in resistance caused by wind, wave, and sea water characteristics. Here, we collate the uncertainty calculation results of the previous sections, select the second group of propulsion test data, and summarize the distribution of the parameters in the propulsion coefficient in
Table 8. As this section only evaluates the uncertainty of the propulsion test, the uncertainty of the resistance test is not considered. On the other hand, the resistance value can also be converted from the resistance test data of the towing tank.
The calculation results substituted by the Monte Carlo program are shown in
Table 9.
Figure 6 shows the distribution frequency of the propulsion efficiency at 200 rpm.
It can be seen from the table that when the propeller speed is 200–300 rpm, the uncertainty of propulsion efficiency is between 4.8% and 11.4%, and the uncertainty decreases significantly with the increase in speed. It should be noted that as it is a single test, the uncertainty evaluation of the propulsion test in this study mainly used the evaluation method of type B uncertainty, and other data fluctuations caused by the environment (random error) were not considered, which will lead to the results in the final evaluation being smaller than the actual value.
The proportion of uncertainty of each part of the propulsion coefficient is shown in
Figure 7. It can be seen from the figure that the effect of resistance correction is the most important as it reached 98%, and the sum of wave resistance and wind resistance accounts for 99% of the uncertainty of resistance correction.
4. Discussion
4.1. Ship Model and Shafting Measurement
The uncertainty related to the conditions of the ship model itself was a small amount in this test. Considering that the interference of these factors to the test was constant throughout the test, their influence was ignored in the evaluation. The Taylor series method was used to analyze the shafting. The uncertainty of shafting conduction has been proven to be very large in a previous ship test [
9,
16]. However, because the broken-shaft-type self-propulsion in this experiment can directly measure the thrust and torque received by the propeller, the uncertainty of this part was found to be very small. In this case, we can better understand the effects of the environment on the experiment.
4.2. Wind
The relative expanded uncertainty of wind resistance in this propulsion test was very large. On the one hand, it was caused by the measurement accuracy of the instrument. After many simulations, it was found that the selection of measuring instruments with higher accuracy of wind speed improved the test accuracy better. On the other hand, the measurement uncertainty of wind resistance has little relationship with the numerical value of wind resistance itself. This means when the wind resistance of the ship model itself increases, the relative measurement uncertainty can be effectively reduced.
Because the velocity of the ship model has a direct influence on the wind resistance, we numerically simulate the uncertainty of the wind resistance at different speeds. It can be seen from
Figure 8 that for the ship model in this test, the relative uncertainty of wind resistance can be controlled below 35% when the speed is above 7 m/s, and the relative uncertainty of wind resistance can be controlled below 25% when the speed of the ship model in this test reaches above 8 m/s. In the case of downwind, when the speed increases to a certain extent, the action direction of the wind resistance will change, and the wind resistance of the ship model will first decrease and then increase with the increase of the speed. Therefore, the uncertainty of the second test will reach a peak when the speed is close to the wind speed.
From the final evaluation of the propulsion coefficient, the wind resistance significantly influenced the test, accounting for about 40%. Therefore, it is important to select the appropriate wind level of the weather and accurately measure the stroke information during the trial. Higher relative velocity will result in greater wind resistance, which is conducive to improving the accuracy of the trial. For this, special attention should be paid to the test under downwind conditions.
4.3. Wave
The relative expanded uncertainty of wave resistance was approximately 56.2%, which shows that the measurement uncertainty of wave resistance is very large. Specifically, it was difficult to measure the wave height of 0.14 m with the wave altimeter in this experiment, mainly because the test site was offshore and the sea conditions on the test day were very mild. Through the numerical simulation of the uncertainty of wave resistance under different wave heights, the curve of
Figure 9 can be obtained. It can be seen that for the equipment and ship model in this test, the uncertainty decreases significantly when the wave height increases from 0.1 m to 0.3 m.
Therefore, for the ship model test on actual seas, choosing a more significant wave environment is beneficial to improve the measurement accuracy. However, when the wave height reaches a certain level, the promotion will be slow, and the conversion of the actual working conditions, the safety of the test and other factors need to be taken into account. For example, in this experiment, the environment needs to be close to the condition of still water, so the wave height is required. The best way for this case is to improve the measurement accuracy of the instrument (wave altimeter). Through calculation, when the measurement accuracy of wave height is doubled, the uncertainty will be reduced to about half of the original.
In addition, since the STWAVE-1 method only considers the influence of meaningful wave height, which may be different from the actual wave resistance experienced by the ship model, especially in strong sea conditions. Therefore the wave resistance estimation method other than the STWAVE-1 method can be adopted to obtain more accurate results by comprehensively considering wave direction, period and other factors.
In the final evaluation, it was found that wave had the greatest effect on this test, which is consistent with our perception in the actual test process. Owing to the requirements of the propulsion test, it is important to choose a relatively calm test site, although this is difficult. Moreover, when we analyzed the data of the wave altimeter, we found that the actual wave height might be inaccurate. This is not only because its true value is exceedingly small, but also because our data can only represent the wave information of a certain point near the ship model. Therefore, it is necessary to obtain real-time and multi-point wave information in future research.
4.4. Other Parameters
Finally, the uncertainty caused by water temperature and speed is also insignificant, which is mainly caused by the accuracy of the instrument itself. Notably, the data in this study are a demonstration of the entire test process. However, being limited by the number of tests and measurement mechanism, it cannot effectively express the fluctuation of the data in the test process. Hence, the data is insufficient for studying the effect of environmental changes. In other words, although the uncertainty of the current is very small in this study, it does not mean that the impact of the current on the test can be ignored. This is partly because currents are difficult to measure with our equipment and may require testing in somewhat more pronounced environmental conditions. On the other hand, we did not consider the influence of random error in this study, so the influence of current in the test was not fully shown. In fact, during the test on seas, we could feel that the effect of the current on the ship model is significant and variable.
According to the results of the evaluation, the effect of wind and wave is the most important factor in this propulsion test. Thus, it is important to improve the test accuracy of wind speed, direction, and wave measurement instruments. In this experiment, the measurement accuracy of wind direction has much room for improvement, and it is necessary to set up two devices to collect wave direction and water depth, respectively. On the other hand, when the speed of the ship model becomes higher, the relative uncertainty of parameters such as wind resistance and speed will be significantly reduced. Therefore, if the conditions permit, we can consider increasing the speed of the ship model to obtain a higher test accuracy.
It can be seen that even though the uncertainty of the environmental parameters in this experiment is very large, the total uncertainty of the propulsion coefficient is relatively small. This is because the wind and wave conditions selected in this experiment are mild, so the corresponding correction value is relatively small. If the experiment is carried out under the condition of a severe environment, the uncertainty of environmental parameters will be smaller, but the correction value will be larger. Therefore, before the test, the research team should estimate the uncertainty according to the test content and the environment at that time, and adjust the test instrument or working condition according to the initial test data. The improvement measures in this study should be helpful to the measurement of most tests on actual seas.
5. Conclusions
In this study, we mainly used the self-adaption program based on the Monte Carlo method to analyze the uncertainty of a ship model propulsion test on a real sea. By comparing the histogram simulated by the Monte Carlo method with the probability density function curve of the normal distribution, only a little deviation was observed between them. Hence, the uncertainty components of each part can be expressed by the results of an approximate normal distribution. Thus, the influence of the test environment was quantified, and the composition of the total uncertainty was analyzed. From the analysis in this study, it can be concluded that the wind and waves in an environment are the two factors with the greatest influence on the test accuracy.
In the subsequent design process of similar tests on actual seas, the experimenter should focus on these two aspects and select appropriate measuring instruments and an environment according to the test content. In addition, an important research work is to quantify the uncertainty of current, and we believe that its influence cannot be ignored.
Owing to the limitations of instruments and test time, this study did not take into account the impact of data fluctuations in the test process. In fact, parameters such as wind, waves and currents in the marine environment are constantly changing, and it may be inaccurate to replace these values with only a few data points or averages. Therefore, the next goal will be to undertake an in-depth study of the sources of uncertainty in the environment by improving the experimental scheme and analysis methods. Furthermore, uncertainty analysis of relevant towing-tank tests in controlled conditions will also be carried out, so as to reflect the interference factors in the actual environment directly.