COMBINING ICE CLASS RULES WITH DIRECT CALCULATIONS FOR DESIGN OF ARCTIC LNG VESSEL PROPULSION

ABS TECHNICAL PAPERS 008 COMBINING ICE CLASS RULES WITH DIRECT CALCULATIONS FOR DESIGN OF ARCTIC LNG VESSEL PROPULSION Sing-Kwan Lee American Bureau of Shipping Originally presented at the Gastech 008 Conference held in Bangkok, Thailand, March 10-1, 008 ABSTRACT Ships navigating in ice areas perform quite differently due to ice resistance as compared to their open sea operations. It is a challenge for ship designers to find a design solution which not only optimizes the propulsion performance in open sea but also provides the ship with good ice performance. Traditionally, many ship designers have used rule formulae in Finnish-Swedish Ice Class Rules (FSICR) to estimate the required engine power. While the use of these formulae have proven to be satisfactory thus far, recently increased growth in the activities in oil and gas production in the Arctic area call into question the continued applicability of these Rules for a very obvious reason. The transportation of these oil and gas products requires ice tankers and ice LNG ships to operate in ice conditions much thicker than those previously assumed in FSICR. In such cases, more and more ship designers are turning to direct calculations and model tests for propulsion designs. In this paper, the two most critical design issues in ice ship propulsion will be discussed, namely engine powering and propulsor performance. Comparisons of propulsion designs based on FSICR and on direct calculation are presented for different propulsors. These include fixed pitch, controllable pitch, and ducted propellers for their performance in ice operation. 1. INTRODUCTION Approximate estimates indicate about one third of the world s known, though not yet developed, reserves of natural gas are in Russia. The overwhelming majority of these reserves are in Artic and Subarctic areas. Russia is not alone in having large natural gas reserves located in Arctic areas, but other countries like Canada and the USA share the same situation with having natural gas reserves in harsh, ice-covered sea environments. As a consequence, the LNG ship technology is moving towards Arctic LNG carriers. New developments in ice navigation and the winterization issues are generating a new challenge for shipping and shipbuilding industries all over the world. 1.1 Design Challenges in Ice Propulsion The economics of LNG transportation by ships from the Arctic are dependent on the efficiency of operation in two areas-- the efficiency of the ship when traveling through the difficult ice conditions and its efficiency in open sea transit. Of special note is the fact that on a direct route from the Arctic to Europe or North America the vessel will spend more than 90% of the time in open sea. In other words, if direct transportation operation is adopted whereby the LNG is loaded on a carrier and transported all the way from the Arctic to its destination, an Arctic LNG vessel is not allowed to sacrifice too much of its open sea performance in order to fulfill its good ice performance. However, designing efficient open sea propulsion and designing good ice propulsion are always in contradiction to each other. In the resistance-propulsion point of view, the contradiction is mainly caused by the fact that high resistance is encountered while operating in low ship speed in ice condition. In such an operation, a propulsor has difficulty performing with its full propulsion capability and efficiency, yet the extreme resistance created by the harsh ice Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion 17

ABS TECHNICAL PAPERS 008 conditions demands that a large thrust be generated enabling the ship to break through the ice. There are many innovative designs proposed to solve this design dilemma, for example, DAT (Double Acting Tanker), podded propulsion, and hull air bubble ejection, to name some of them. Indeed, new design experience continues to grow in industry. Research and development activities have been quite active in recent years in the area of new ice propulsion design. In this study, however, we put more focus on the traditional propulsion concepts such as CPP (Controllable Pitch Propeller) and ducted propeller (nozzle propulsor) instead of the edge-cutting design concepts, while trying to explore more of their potentials in ice propulsion. 1. Integrated Design Consideration In ice going ship design, the strengthening of the hull, rudder, propellers, shafts and gears is clearly related to the safety of navigation in ice. In principle, all parts of the hull and the propulsion machinery exposed to ice loads have to be ice-strengthened. Hence, an ice propulsion design cannot be successful with just the use of the propulsion capability, but should include both the thrust and the propulsion power. A comprehensive design should include strength design of the propulsor under ice load. In addition, propeller cavitation and its induced vibration must also be considered. It should be noted that an ice-strengthened blade design makes it much easier to trigger cavitation problems due to its increased thickness as compared to a non-ice propeller design. For blade strength assessment, finite element analysis is usually requested to be performed based on the ice loads and failure criteria proposed in IACS Polar class rules (IACS, 008; Lee, 007). Analyses of cavitation for propeller blade and its induced hull vibration are relied upon more in the direct calculation approach. For more details of the issues refer to (Lee, 006).. ICE RULES FOR SHIP POWERING.1 Power Requirement in Ice Rules Generally, FSICR (Finnish-Swedish Ice Class Rule) is the accepted Rules used for vessels trading in the Northern Baltic Sea in winter. The Rules primarily address matters directly relevant to the capability of ships to advance in ice. Ships are required to have a certain speed, 5 knots, in a brash ice channel in order to ensure the smooth progress of traffic in ice conditions. The required minimum speed, 5 knots, is based on the maximum waiting time, 4 hours, for icebreaker assistance. According to Section in FSICR (00), for ships entitled to ice class IA Super, IA, IB, or IC, the engine output is to be not less than determined by the following formula. ( RCH /1000) P = K (1) C D p / where, K C is a propeller-related constant factor, D p is propeller diameter and R CH is resistance, which can be calculated by the following rule formula. R CH = C 1 + C + C C μ (H F + H M ) (B + C Ø H F ) + C 4 L PAR H F + C5 LT AWL () B L Detailed meanings of the symbols used in the formula can be found in Section.. of FSICR and are not repeated here. When using equation, the following restriction is to apply: LT 0 B 5 where L is length of perpendiculars, B is maximum breadth, and T is draught. In fact, this restriction constrains the use of the engine output formulae for large size ship. For Arctic ships, although IACS Polar Rule for Machinery Requirement, URI, will soon be published (IACS, 008), the powering for Arctic navigation is not included. 18 Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008. Difficulties in the use of FSICR Difficulties in using FSICR for a rational design of propulsion power are learned by shipyards and classification society such as ABS when applying FSICR to large ships (Bezinovic, 005, Kim et al., 005, Lee, 006). For a FPP (Fixed Pitch Propeller) design for an Aframax IA ice class tanker, it is found that the engine power can be largely reduced based on a more accurate direct calculation for propulsor thrust and model test results for ice resistance. To illustrate the general idea of the conservativeness of FSICR, Figure 1 shows a comparison done by an engine manufacturer (MAN B&W, 00, 006) for the non-ice class tankers engine powering and FSICR requested engine powering for ice tankers. As seen, the Rule requests an Aframax size IA ice class tanker to install the engine power with a non-ice VLCC tanker. It not only causes the extensive increase in cost in building the ice tanker but also creates an extremely difficult situation in trying to fit the huge engine into an Aframaxsize engine room. In fact, in determining the ice propulsion powering, ice resistance and propulsor performance are two dominant factors, which are related to the coefficients, R CH and K C in FSICR (see the formula (1)). The coefficients are empirically estimated based on average ice resistance and propeller performance in FSICR. With the continuously, everincreasing introduction of new ice ship designs such as the new hull and bow forms designed to reduce ice resistance as well as more efficient propellers, FSICR ice powering formulae have come into a questionable review on the current validity of the rationale.. Furthermore, in the case of these Arctic ships, it has not been verified that the ice resistance and ice powering formulae in FSICR are applicably suited to meet the demands of thicker multi-year ice propulsion. 1a) Engine power demand for non-ice tankers 1b) Engine power demand for Finnish/Swedish ice classed vessels based on FSICR calculations Figure 1 Ice and non-ice tanker engine power comparison. ICE POWERING DETERMINATION BASED ON DIRECT CALCULATION A comprehensive design for ship ice propulsion involves many different elements such as hull resistance, propulsor design, engine character, shafting/gear/bearing lubrication and friction loss, CP mechanism operation, and sometimes propeller cavitation induced by machinery vibration etc.. Although all of them are more or less influencing the finial ice propulsion performance, ice resistance, propulsor performance and engine Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion 19

ABS TECHNICAL PAPERS 008 character are the three most critical areas in ice powering determination. They are addressed briefly in the following sections..1 Resistance In addition to the concern of open sea resistance when working on ice propulsion design, ice resistance is another factor to be considered for power determination. Ship resistances in ice and open sea can be estimated by numerical simulation or model test measurement. As CFD simulations have become more and more practical for the study of open sea resistance in recent years (MARNET, 00), it is anticipated that these simulations will play an increasingly more important role in resistance prediction. However, in the current stage, CFD is still state-of-the-art and used as a compromising method with model test results for resistance prediction. For ice resistance simulation, some progress has been made in numerical simulation for level ice resistance (Valato, 001) but the brash ice resistance under waterline is still not certain and needs to use the semi-empirical model of Lindqvist (1989). In general, model tests are recommended by the Guidelines of FSICR (005) for ice resistance assessment, especially when vessel displacement is greater than 70,000 tonnes. As model test uncertainty is always an issue, it should be noted that the difference in ice modeling and measurement procedure may come up with a large discrepancy on ice resistance. To assist designers in performing ice model test, some aspects of the ice model tests such as the geometry of the ice channel, ice friction coefficient, determination of the propulsion power in full scale, and model test documentation are addressed in the Guidelines of FSICR (Appendix 4). Ice resistance can usually be expressed as a quadratic curve. R ice = d 0 + d 1 V S + d V S () where R ice is ice resistance, V S is ship speed and d 0, d 1, and d are curve fitting constants. Resistance tests usually are performed based on experiences and calibrations in different towing tanks. ITTC also recommends some standard procedures for resistance test and performance prediction method (ITTC, 000 and 1999). If open sea resistance for designed speed is obtained, resistance at other speeds can be calculated based on cubic law as follows: R OW = c 1 V S + c V S + c V S (4) where R OW is open sea resistance, V S is ship speed and c 1, c, and c are curve fitting constants, which can be determined based on resistance and ship speed measured in model tests. Figure is the typical outlook of R ice and R ow curves. As seen, ice resistance usually has higher value than open sea resistance under the same ship speed. However, the increasing rate of the ice resistance curve is smaller than the open sea resistance curve. It is also the reason behind the use of a quadratic curve for ice resistance instead of a cubic curve. Figure Typical ice and open sea resistance curves. Propulsor Performance As is well known, propeller thrust and torque are dependent on its geometry and operating condition (rpm and inflow velocity toward propeller). They can be expressed as nondimensional forms as follows: 0 Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 K T = a + (5a) 1 + aj + aj a4j 1 + b J + b J b4j K Q = b + (5b) where K T is thrust coefficient defined as T/ρn D 4 and K Q is torque coefficient as Q/ρn D 5, in which T and Q are thrust and torque. J is advanced coefficient defined as V/nD; ρ, n, V, and D are fluid density, propeller rps, inflow velocity, and diameter. For a specific design, {a 1 a a a 4 } and {b 1 b b b 4 } are constant vectors. In additional to propeller open water model test, K T and K Q curves can also be constructed by CFD simulation (Chen and Stern, 1999, Chen and Lee, 005) or regression formulae (for standard series propellers, see Oosterveld and Ossanen, 1975). In general, equations 5a and 5b can be extended to any higher degree polynomial based on the propeller open water measurement data, however, around the range [0, 1] of J, third order polynomial usually provides sufficient accuracy for K T and K Q values. This can be seen obviously from the typical pattern of K T and K Q curves shown in Figure for both model and full scale Reynolds number R n.. Ship/Engine/Propeller Interaction Propeller absorbing power (see propeller demand curve shown in Fig. 4) is usually limited by engine power curve (see torque limit line in Fig. 4). For an appropriately designed propeller, it should absorb all engine power at MCR rpm when advancing with the design ship speed. For Controllable Pitch Propeller CPP), it can adjust its pitch to keep its original hydrodynamic performance at different ship speed. Thus, CPP can approximately maintain its rpm as the MCR engine speed and absorb the MCR engine power at different ship speed. However, for Fixed Pitch Propeller (FPP), when ship speed becomes slower (operating in ice condition), in order to maintain the rpm as MCR speed, the propeller needs to be supplied with more power from the engine to overcome the higher hydrodynamic torque on the blade compared to the original ship speed condition. This leads to the original propeller demand curve shifting up to curve A as shown in Fig. 4. As seen in Fig. 4, due to the limitation of the engine torque limit curve, when the propeller operates off its design condition (design ship speed along with MCR rpm), the engine can only provide the power along the torque limit curve and accordingly the propeller will rotate slower (see point x of curve A in Fig. 4). At point x, engine power P is equal to propeller absorbing power, i.e., P p = πnk Q ρn D 5. Using equation 5b and definition of J, we have 5 V V V P = πρ D [ b1n + b n + b n + b4 ] (6) D D D Usually, engine power along the torque limit line follows the following relationship with rpm (MAN B&W, 1994): P rpm (7a) Figure Typical propeller performance curves (K T and K Q ) or more precisely n P = P ( ) MCR (7b) nmcr Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion 1

ABS TECHNICAL PAPERS 008 where P = engine power corresponding to rotation speed n P MCR = engine power corresponding to MCR rotation speed n MCR Based on (7b) and (6), if propeller rotation speed n at off design operating condition is given, the engine MCR power can be determined by the following equation. P MCR n n MCR 5 = πρd [ b1n + b V n D V V + b n + b4 ] D D Figure 4 Engine and propeller demanded curve (8) For ship under ice operating condition, the determination of actual propeller rotation speed n in equation (8) is related to the balance between propulsor thrust and ice resistance. Details are given in the following. When the ship operates in ice condition, the propeller thrust generated should be used to overcome both hydrodynamic resistance (or open sea resistance) and ice resistance, i.e., T R OW + R i = R ice, where T is propeller thrust, R OW is open sea resistance, and R i is net ice resistance, and R ice is the total resistance under ice navigation. In other words, T should be large enough to overcome the total resistance, R ice. Accordingly, we have T R ice 4 ρ n D K T d 0 + d1v S + d V (9) S ρn D + V + nd V nd + V nd 4 [ a1 a a ( ) a4 ( ) ] d 0 d1v S d VS + + It should be noted that V and V S in equation (9) represent different physical meanings. The former is the actual velocity toward the propeller and V S is the ship speed. Usually, V can be roughly estimated as V S (1-w), where w is a wake fraction due to hull boundary layer effect. In equation (9), the thrust reduction effect is not considered as it becomes smaller when the ship is in slow speed. Based on equation (9), if total resistance R ice, ship velocity V S, wake fraction w, are given, the actual propeller speed n under ice condition can be calculated..4 Procedure for Power Determination In summary, the overall procedure for determining required engine power is as follows: Perform model test and measure the open sea resistance at design ship speed Construct total resistance curve equation () based on ice resistance model test results Construct propeller K T and K Q curves equations (5a) and (5b) For FPP, use equation (9) to determine the actual propeller rotation speed n; for CPP, n is same as the engine MCR rotation speed but propeller performance curves (K T and K Q ) will be based on those at new turning position (more details for CPP will be illustrated in the Case Study section). Use equation (8) to determine the required engine MCR power. 4. CASE STUDY The aforementioned powering determination procedure is used in this section for an Arctic LNG carrier. Different propulsor designs, namely FPP, CPP, and ducted propeller, are compared for their propulsion powers for this LNG carrier in Baltic ice sea and Arctic sea operations. Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 4.1 Twin Screw LNG Vessel LNG ships have traditionally been a single screw ship driven by a steam turbine. The size of the steam turbine plant is relatively large and the turbine is expensive both to manufacture and to maintain. Therefore, a twin steam turbine design has not been an attractive solution. New and improved techniques have resulted in more efficient reliquefication plants to take care of the LNG boil off and more efficient dual fuel engines. The expansive steam turbine is no longer the only choice for the LNG ship. Model tests for ice resistance and open sea resistance are conducted for this shuttle. The ice thickness range for ice resistance is from 0.6 m to.4 m with ship speed range between 0 to 18 knots. The resistance curves are reproduced in Figure 6. Later, these resistance curves will be used to estimate the requested engine power and propulsor operation conditions such as rpm for FPP and turning angle for CPP under different ice conditions. An Arctic LNG Shuttle with twin screws (Figure 5) and the following design parameters are selected for this propulsion power comparison study. In this study, the selected LNG ship is assumed to be driven by two Diesel engines. Length = 00 m Breadth (water line) = 8 m Draft (max.) = 9 m Cargo capability = 40, 000 m Figure 5 A twin screw Arctic LNG Shuttle (after http://www.hoegh.com) According to the information posted in the internet site (http://www.hoegh.com), the LNG shuttle is designed to transport LNG through areas with severe ice, such as the Russian Arctic. It has the capability to operate independently (without ice breaker escort) in up to.0 meter nominal ice thickness. Figure 6 Ice and open sea resistance curves In Figure 6, Baltic ice classes IA, IB and IC are indicated as the curves corresponding to the ice thickness H ice = 0.6m, 0.8m and 1.0m. Due to the limited data releasing in the public internet site, the ice thickness for the ice resistance curves is only up to H ice =.4m although the ship is designed to operate with ice thickness up to m. It is noted that ice resistance increases rapidly with the ice thickness and compared to the open sea resistance R OW, the net ice resistance R i dominates the total resistance R ice. To have a general idea of percentage for different resistances, the comparison among R OW (open sea resistance), R i = R ice R OW (net ice resistance), and R ice (total resistance) under ship speed 5 knots is summarized as follows: Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 Ice thickness H ice 0.6 m (IC class) Ship speed V S Open sea Resis. R ow Total Resistance R ice Net ice Resistance R i = R ice R OW 5 knots 157.0kN 114.4 kn 1057.4 kn (87% of R ice ).0 m 5 knots 157.0 kn 5000.0 kn 484.0 kn (97% of R ice ) 4. Powering for Baltic Ice Sea Navigation As explained in section, in determining the engine power for ice navigation a sequence of calculations needs to be performed based on the information of ice/open sea resistance, propulsor performance and engine character. To initiate this calculation procedure, propeller diameter and its design rpm need to be determined first. For the selected LNG carrier, a propeller (including FPP, CPP and Ducted propeller) with 6 m diameter and 10 rpm is a reasonable choice. For rotation rpm = 10, one option of the engine can be the following. MCR rpm =10 Cylinder. MCR power in kw number 5 585 6 6990 7 8155 8 90 Table 1 A example engine for MCR rpm = 10 Figure 7 AU-CP Propeller geometry - 4 blades with EAR 55% and 70% 4..1.1 Propeller curves K t and K q In FPP case calculation, although the propeller is a CPP design, propeller turning angle is fixed as θ = 0 o, which is corresponding to the pitch ratio P/D = 1.0. The propeller performance curves, K T, K Q, and efficiency η, can be found from the pervious model test results (Yazaki, 1964, Yazaki and Sugai, 1966). For design pitch ratio P/D = 1.0 under θ = 0 o position, the K T and K Q curves (Figure 8) can be expressed as follows: K t + = 0.47 0.4J 0.16J 0.04J (10) 10K q + J = 0.69 0.46J 0.17J 0.001 (11) 4..1 Fixed Pitch Propeller Design In order to compare the engine power requested due to the performance of FPP and CPP, the same propeller geometry AU-CP 4 blades EAR = 55% is used for both FPP and CPP cases. Figure 7 is an outline plot of the propeller blade with two different EAR 55% and 70%. Figure 8 K t, K q and η for AU-CP propeller (4 blade, EAR = 55%) 4 Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 4..1. Ice resistance for IA, IB and IC classes Under 5 knots speed (FSICR requirement), the ice resistances for IA, IB and IC class can be estimated using the model test results shown in Figure 6. The following table summarizes the results of the total ice resistance curves (based on order curve fitting) and the value of R ice at ship speed V s = 5 knots. IA: H ice = 1.0 m IB: H ice = 0.8 m IC: H ice = 0.6 m Ice resistance curve R ice at 5 knot speed R ice = 1114.8 + 191.1V 157.8 kn S R ice = 846. + 14.7V 1669.55 S kn R ice = 59.4 + 10.9V S + 114.4 kn Table Ice resistance for IA, IB and IC classes IA: 1 157.8 /. = 1078.9 = 0.69[0.47n 7.49n 68.4 + 64.01 ] n IB: 1 1669.55 = 84.775 = 0.69[0.47n 7.49n 68.4 + 64.01 ] n IC: 1 114.4 = 607. = 069[0.47n 7.49n 68.4 + 64.01 ] n (1a) (1b) (1c) By solving the equations (1a), (1b) and (1c), the rpms for the AU-CP FPP under IA, IB and IC ice conditions are obtained as follows: Class 0.5 R ice Actual rpm of FPP IA 1078.9 kn 88.0 IB 84.7 kn 79.0 IC 607. kn 69.0 Table Actual rpm of FPP 4..1. FPP rpm under ice operations To determine the actual rpm for the FPP under ice operation, equation (9) is used for the three scenarios of IA, IB and IC ice resistance. Since a twin screw design is used, the balance of one propeller thrust should be equal to half of the R ice, i.e., 0.5R ice = ρn D 4 K t. Therefore, we have (see also Figure 9). Figure 9 Thrust and R ice curves 4..1.4 Engine power for ice navigation Once the actual propeller rpms under the ice operations are known, the engine power needed to reserve in MCR speed can be calculated based on equation (8). For the MCR rpm = 10, the requested engine powers for IA, IB and IC classes are 11.98, 10.48 and 8.8. The details of the calculations for these engine power values are summarized in the Table 4. It is also interesting to compare these direct calculation results with FSICR calculations to see how much reduction can be obtained based on the current approach. In FSICR, the engine power is calculated based on equation 1. For the FPP twin case, K c is taken as 1.6. The following is a table summarizing the comparison. Class R ch FSICR Power IA 157.8 kn 6.79 IB 1669.55 kn 18.191 IC 114.4 kn 11.85 Table 5 FSICR requirement for engine powers for IA, IB and IC classes Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion 5

ABS TECHNICAL PAPERS 008 Class n J 10Kq (eq. Propeller absorbing n MCR = 10 rpm actual rpm 11) power (n MCR /n) P MCR (eq. 8) P = πρn D 5 K Q IA 88.0 0.5 0.566 8.94 1.4 11.98 IB 79.0 0.8 0.550 6.87 1.667 10.48 IC 69.0 0. 0.50 4.06.185 8.819 Table 4 Engine powers for IA, IB and IC classes Through this comparison, it is found the engine powers requested in FSICR for this FPP case are % of the power from the direct calculation for IA class, 174% for IB class and 18% for IC class. 4..1.5 Open sea performance To ensure this FPP and engine powering combination to be suitable for open sea operation, the propeller performance needs to be checked at the MCR condition in open sea. In this study, the MCR condition assumes that a full speed of the LNG is 17 knots under a propeller rotation speed = 10 rpm. Under the 17 knots speed, the flow velocity in front of the FP propeller is estimated to be 14.57 knots. Here, the wake factor is taken as 0.14 (Manen and Oossanen, 1988) based on the LNG ship block coefficient range, 0.7 ~ 0.8. Accordingly, the advance coefficient J is calculated to be 0.7 for rpm = 10. From Figure 8 (propeller performance curves), we notice the FPP selected is operated in good efficiency (see the η curve). The K t and K q values under the J = 0.7 can be obtained from equations (10) and (11). Based on the K t and K q values calculated, the thrust, T, and propeller absorbing power, P, can be determined. The details are summarized in the following table: rp m 10 J K t K q T = ρn D 4 K t 0.7 0.1 6 0.0 78 614. 4 kn Q = ρn D 5 K q 640. 6 knm P =πn Q 6.84 0 Table 6 Propeller thrust and absorbing power at open sea operation Based on the propeller thrust value, the twin propulsion thrust is 18.48 kn. In open sea operation, the resistance R at 17 knots can be calculated by the following equation. It is about 119.55 kn. R + = 5.44V + 0.59V 0.1V (1) In summary, the selected FP propeller can generate enough thrust to propel the LNG ship with 17 knots speed. The requested engine power is in the level for the ship operating in 5 knots in IB class ice condition (see table 4). Basically, the engine powers (IA, IB and IC classes) requested to reserve in MCR condition are large enough for open sea condition (see Table 4). 4.. Controllable Pitch Propeller Design Unlike FPP, CPP can change its pitch by turning the blade in certain angle Δθ (Figure 10) to obtain more thrust for propulsion. Usually the propeller rpm can be maintained as the MCR speed even at low ship speed in ice navigation. These advantages make CPP more capable and efficient in ice operation. Figure 10 Turning angle Δθ for CPP to adjust its pitch 6 Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 4...1 Blade turning angles at ice operations To determine what θ should be turned for the blade to generate enough thrust to overcome ice resistance, propeller performance (K t and K q ) curves at different turning angles should be obtained first. Figure 11 shows the thrust, torque and efficiency curves for the AU-CP propeller at different blade turning angles (Yazaki, 1964). consideration. For the three Baltic ice classes, IA, IB and IC, their values are given in the following table. In C ice calculation, the ice resistance formulae in Table are used for the ice resistance R ice calculations. The J is calculated based on inflow velocity in front of the propeller which takes into account the wake factor w = 0.14. Class rpm J R ice / C ice IA 1078.9 kn 0.81 IB 10 0.16 84.8 kn 0.17 IC 607. kn 0.158 Table 7 C ice values for IA, IB, and IC ice classes In Figure 1a, the J-C ice points for IA, IB and IC classes are plotted in the J-θ-K t curves. Through interpolating the C ice values into the K t values, the turning angles θ can be calculated. They are -5.09 o for IA class, - 8.7 o for IB class, and -11.48 o for IC class. Based on these θ values, K q values corresponding to the θs are calculated and plotted in Figure 1b. 11a) K t and η at different θ 11b) K q and η at different θ Figure 11 Thrust, torque and efficiency curves for the AU-CP propeller at different θ To interpolate the turning angle value from the above J-θ-K t curves, the operation point (J- C ice ) for ice resistance should be determined first. Here C ice is the ice resistance coefficient non-dimensionalized by ρn D 4 ; here factor is taken for twin screw propulsion Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion 7

ABS TECHNICAL PAPERS 008 Class R ch FSICR Power IA 157.8 kn 4.056 IB 1669.55 16.7 kn IC 114.4 kn 10.156 % the ratio of FSICR power to direct calculation power 79 % 71 % 56 % Table 9 Power comparison between FSICR and direct calculations 1a J-C ice points for IA, IB and IC classes 1b K q values for IA, IB and IC classes Figure 1 CPP turning angle interpolations 4... Engine power for ice navigation The engine powers required for the ice classes, which are equal to the propeller absorbing powers, are calculated based on the K q values obtained, i.e. P = πρn D 5 K q. They are summarized as follows: Class n CPP rpm J θ 10Kq Propeller absorbing power P = πρn D 5 K Q IA -4.18 o 0.87 8.6 IB 10.0 0.16-7.44 o 0.685 6.040 IC - 10.7 o 0.17615.96 Table 8 Engine required powers for IA, IB and IC classes for CPP designs As done in the pervious FPP design case, FSICR formulae are used to check the Rule required engine powers for this CPP design. For a twin CPP propeller, K c value is taken as 1.44 for equation 1. Comparisons between FSICR and the direct calculations are given in the following table. For open sea operation, since the design is same as the FPP case, the engine power for the LNG operates in 17 knots is with the same value as in the table 6, i.e. 6.84, which is 79.% of the IA class engine power in this CPP design. 4.. Ducted Propeller Design It is well-known that a ducted propulsor has better propulsion capability when operating in slow speed compared to an open propeller. This is not only due to the better design of the propeller at low J operation but also due to the contribution of the thrust generated by the duct (see Figure 1). As seen in Figure 1, a typical thrust generated by the duct is more than half of the propeller thrust at low J condition (J = 0.1). However, it also needs to be noted that a duct will become a disadvantage when a vessel navigates in high speed. As seen in Figure 1, thrust over J = 0.6 from the duct becomes negative, i.e. extra resistance is generated by the duct. Later, in this case study it will be shown that the use of the ducted propeller can create the problem whereby the LNG carrier cannot reach a 17-knots speed in open sea operation although it does reduce quite a lot of the engine power required for ice propulsion. 8 Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 Figure 1 Thrust coefficients of duct and propeller for a ducted propeller. Figure 14 Geometry of ducted propeller 19A Ka 4-70. 4... Ducted propeller curve K t and K q 4...1 Ducted propeller 19A Ka-4-70 A fixed pitch ducted propeller selected in this study is a standard design of a 19A nozzle with a Ka 4-70 propeller (4 blade with EAR = 70%). In order to compare the final engine power results with the pervious open propellers under an approximate base, the same propeller diameter, 6 m, as the open propeller is adopted to this ducted propeller. Figure 14 shows the outline of the 19A Ka-4-70 ducted propeller (blade and nozzle profile). The propeller performance curves of this ducted propeller are available from the wellknown model tests done by van Gent and Oosterveld (198). Figure 15 shows the performance curves. It should be noted that the K t curve represents the total thrust from propeller and duct. Similar to the pervious open propeller, the K t and K q curves can be expressed as cubic curves. Their mathematical forms are as follows: K t = 0.56 0.599J + 0.85J 0.91J (14) 10K q = 0.44 0.014J 0.85J 0.116J (15) Figure 15 K t, K q and η for 19A Ka 4-70 ducted propeller (4 blade, EAR = 70%) 4... Operating rpm determination Since this ducted propeller is a fixed pitch design, propeller rotation speed at ice operation will be dropped due to the heavy loading condition encountered. Based on the thrust-ice resistance balance principle (equation 9), the ice operation rpms for IA, IB and IC class can be determined through the similar procedure as the FPP case. The solutions are summarized in the following table. Class 0.5 R ice Actual rpm of FPP IA 1078.9 kn 87.0 IB 84.7 kn 78.0 IC 607. kn 68.5 Table 10 Actual rpm of ducted propeller Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion 9

ABS TECHNICAL PAPERS 008 To crosscheck the solutions, a plot for the ice resistances (IA, IB, and IC classes) and the ducted propeller thrust curve is drawn (Figure 16). The intersection points of the thrust curve and the ice resistance lines provide the actual rpm values. As seen, the rpm values of the intersection points are close to the values in the Table 10. Similar to the FPP case, a comparison of these direct calculations with FSICR calculation is preformed. In FSICR, engine power is calculated based on equation 1. For a FPP twin screw, K c is taken as 1.6. For a ducted propeller, a further 70% reduction from the open propeller result can be applied according to the Guidelines for the Application of FSICR. The following is a table summarizing the comparison. Class R ch IA 157.8 kn IB 1669.55 kn IC 114.4 kn FSICR Power Open propeller 6.79 18.191 11.85 Power for ducted propeller 70% power of open propeller 18.71 1.7 7.90 Table 1 FSICR requirement for engine powers for IA, IB and IC classes Figure 16 Thrust and R ice curves 4...4 Engine power for ice navigation Based on the actual rpms in table 10, the engine powers needed to reserve for MCR condition can be determined based on the equation (8). The details are summarized in the following table. Class n actual rpm J 10Kq (eq. 15) Propeller absorbing power P = πρn D 5 K Q n MCR = 10 rpm (n MCR /n) P MCR (eq. 8) IA 87 0.55 0.419 6.97 1.7 8.764 IB 78 0.85 0.41 4.544 1.71 7.770 IC 68.5 0.4 0.404.010. 6.68 Table11 Engine powers for IA, IB and IC classes From Table 11 and Table 1, it is found that FSICR requests 1% engine power of the value from direct calculation for IA class; 164% for IB class; and 118% for IC class. 4...5 Open sea performance Although the ducted propeller has good performance in ice operations, since its performance is usually bad in fast speed, a crosscheck of its open sea performance is necessary. Assuming that 17 knots is the design full speed and 10 rpm engine is used same as before, the thrust and the engine power can be calculated based on the approach used in the open FPP case. Table 1 summarizes the details of the calculations. rp m J K t K q T = ρn D 4 K t Q = ρn D 5 K q P =πn Q 10 0.7 0.0 0.0 7. 555. 5.81 0 Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 8 66 8 kn 0 kn m Table 1 Propeller thrust and absorbing power at open sea operation As known earlier, at 17 knots speed, the ship resistance is 119.55 kn. From Table 1, the twin screw propulsion from this ducted propeller design is 675.6 kn which is almost 50% smaller than the requested thrust. In other words, the LNG ship can not reach the 17 knots speed by this propulsor design. It should be noted that the problem is not from the engine power but mainly due to the available thrust at the full speed operation. In fact, the engine requested power in this case is smaller than the open FPP case. To estimate the ship speed for this ducted propeller design, the thrust curve for 10 rpm rotation speed based on K t curve (.0 equation 14 a twin crew propulsion) and the open sea resistance curve (equation 1) are plotted in a same diagram (Figure 17). At the intersection point of the two curves, the propeller thrust becomes equal to the open sea resistance and the ship speed at the point is 1.5 knots. For this speed, the engine power requested is about 8.94 for each ducted propeller. Details of the calculations are summarized in the following table. Rpm V s Thrust per propeller 10 1.5 knots Torque per propeller 60 kn 785.8 kn-m Power per propeller 8.94 Table 14 Thrust, torque and power for open sea operation Figure 17 Propeller thrust and open sea resistance curves 4. Powering for Arctic Sea Navigation To illustrate how the ice powering rapidly increases with the ice resistance, the powering for Arctic ice, which is with thicker ice thickness than Baltic ice class, is also studied. As the ice resistances in thick ice conditions are extremely high and if the FSICR ship speed 5 knots is used, it may come up that the ice powering is irrationally high to select a reasonable size engine to fit into an engine room. Hence, in determining the Arctic ice powering, a slower speed.5 knots is used as the operating speed for the Arctic ship design. Basically, all the calculations are the same as that were previously performed for the FPP, the CPP and the ducted propeller. The final results for the FPP, the CPP and the ducted propeller for the ice thickness range from 1. m to.4 m are summarized in the following tables. Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion 1

ABS TECHNICAL PAPERS 008 H ice Total R ice (kn) Propeller rpm Propeller absorbing n MCR = 10 rpm power () (n MCR /n) P MCR () (m) 1. 1994.69 80.5 7.465 1.6 1.09 1.4 456.41 88.6 10.14 1. 1.489 1.6 905.6 95.7 1.88 1.14 14.69 1.8 400. 10..0 950. 111.0. 4450. 117.6.4 5000. 14. Required propeller rpm > 10; this cause the engine being overloaded to generate the thrusts for.5 knots propulsion in 1.8 m upper ice thickness.. Table 15 Engine power for.5 knot Arctic sea navigation FPP design H ice Total R ice (kn) CPP turning n MCR = 10 rpm; P MCR () (m) θ 1. 1994.69-7.9 o 7.001 1.4 456.41-5.0 o 9.99 1.6 905.6 -.4 o 1.685 1.8 400. 0.5 o 16.554.0 950..9 o 1.946. 4450. 6.9 o 7.964.4 5000. 10. o 5.08 Table 16 Engine power for.5 knot Arctic sea navigation CPP design H ice (m) Total R ice (kn) Propeller rpm Propeller absorbing power () n MCR = 10 rpm (n MCR /n) P MCR () 1. 1994.69 77.98 4.78 1.71 8.177 1.4 456.41 85.77 6.78 1.41 8.99 1.6 905.6 9.74 8.081 1.1 9.778 1.8 400. 99.8 10.1 1.04 10.50.0 950. 107.11 Required propeller rpm > 10; this cause the engine being overloaded to generate the thrusts for.5 knots. 4450. 11. propulsion for.0 m upper ice thickness.4 5000. 119.7 Table 17 Engine power for.5 knot Arctic sea navigation Ducted propeller design According to these case studies, it has been found that the ducted propeller is the most energy-saving design for ice operations for Arctic sea propulsion in the high ice resistance (H ice = 1.4 1.8 m ). However, due to the limitation of its fixed pitch design, the ducted propeller cannot operate in the condition with ice thickness thicker than m. For the CPP design, although the power-saving capability is not as good as the ducted propeller, the propulsion capability in ice is unlimited for any ice thickness providing that enough power can be supplied. As expected, the FPP design has the worst performance in ice navigation as it needs the largest amount of power among the other designs for the same ice thickness, and it becomes unable to operate in the ice thickness equal to 1.8 m and upwards no matter if the engine power continues to increase. However, it should be noted that for Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 the open propeller design, both CPP and FPP, the open sea performance is much better than the ducted propeller. 5. CONCLUDING REMARKS Driven by the large reserves of oil and natural gas and their exploitation in Arctic and subarctic areas, ship propulsion in ice-covered sea becomes a critical concern. The design dilemma of ice propulsion performance against open-water performance makes ice going ship design more challenging. In this paper, commonly used propulsors, namely FPP, CPP, and the ducted propeller, were investigated to compare their performance in ice propulsion conditions. A detailed procedure based on direct calculation for determining propulsion thrust and powering was documented. To illustrate how to process the calculation procedure to obtain the required engine MCR power, an Arctic LNG carrier with ice resistance information available was used as an illustrative example. The overall results of powering for the FPP, the CPP and the ducted propeller for the LNG vessel are summarized again in the following table to highlight the findings in this study. Ship speed V s 17 kts 5. kts.5 kts H ice (m) Open IB IC IA sea 0.6 0.8 1.0 1. 1.4 1.6 1.8.0..4 Resistance (kn) 119.6 114.4 1669.5 157.8 1994.7 456.4 905. 400.0 950.0 4450.0 5000.0 P MCR 6.8 8.8 10.5 1.00 1.1 1.5 14.6 x x x x FPP P io 4.0 6. 8.9 7.5 10.1 1.9 x x x x NA Power rpm 69 79 88 80 89 96 10 111 118 14 () P FSICR NA 11. 18. 6.7.8.5 41.8 5.9 66. 79. 94. P MCR 6.8 4.0 6.0 8.6 7.0 9.4 1.7 16.6 1.9 8.0 5. CPP P io 4.0 6.0 8.6 7.0 9.4 1.7 16.6 1.9 8.0 5. NA Power θ -11.5 o -8. o -5.1 o -6.9 o -.9 o -1.1 o.1 o 5.6 o 8.8 o 1. o () P FSICR NA 10. 16.4 4.1 1.4 9. 7.6 47.6 59.6 71. 84.9 P MCR x 6.7 7.8 8.8 8. 9.0 9.8 10.5 x x x Nozzle P io.0 4.5 6.4 4.8 6.4 8.1 10.1 x x x NA Power rpm 69 78 87 78 86 9 100 107 11 10 () P FSICR NA 7.9 1.7 18.7 16.7.8 9. 7.0 46. 55.4 66.0 Table 18 Comparison of powering for FPP, CPP, and Ducted propeller based on FSICR and direct calculation In the Table 18, both the direct calculations and FSICR calculations are included. P MCR represents the final engine power required to be installed; P io is the power required under different ice operation conditions, for instance, ice operation conditions of IC class are ship speed 5 knots and 0.6 m first year ice; P FSICR is the power calculated based on the FSICR. The FSICR calculations are extended to higher ice condition beyond the FSICR classes, in which ice thickness is 1.0m. In principle, FSICR can be applied if ice resistance and propeller size (diameter) and type are known (see equation 1). The rpms provided in the table for the FPP and the nozzle (the ducted propeller) are the propeller rotation speed at ice operation. For the value θ, it is the turning position of the CPP blade at ice propulsion. As well-known, the CPP can maintain its MCR rotation speed even in high load condition under ice operation. The symbol x means the propeller design cannot be applied. For example, x in the table for P MCR of the ducted propeller (nozzle) means the open sea speed 17 knots cannot be fulfilled due to the shortage of the thrust generated by the ducted propeller. And x in the FPP for P io means no 10 rpm engine can be used to fulfill the ice propulsion for the ice conditions H ice > 1.8 m, due to the rpms requested are higher than 10 in that ice condition. From Table 18, the concluding remarks of this study are summarized as follows: Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

Design dilemma of ice propulsion against open sea propulsion is obvious in the FPP design. As seen, despite the good performance in open sea for the FPP (the propeller fulfills the full speed and is operated in high efficiency), its ice propulsion is not satisfactory. First, it requires quite high power to reserve in the engine in order to perform the ice propulsion at low ship speed and low propeller rpm. Second, for more severe ice conditions (H ice > 1.8m), no matter how large the engine power provided, the FPP still cannot generate enough thrust to propel the ship at the requested speed (.5 knots). The ducted propeller has the best performance for ice propulsion in severe ice conditions (see H ice = 1.4 ~ 1.6m) but its open sea operation is the worst. If the CPP design is further applied to the ducted propeller, its engine power P MCR will be the value of P io and becomes much smaller compared to open CPP for ice operation (see P io for Nozzle and P io for CPP in the table). In fact, there are more advantages associated with a ducted CPP in Arctic ice propulsion. For instance, in the strength point of view, the use of a duct to protect propeller blades will make the CPP mechanism safer in ice operation. The CPP has the best overall propulsion performance among the designs. It fulfills the open sea operation well and its ice performance is also good. As seen in Table 18, its required engine power is the smallest among the other designs for the H ice 0.6m ~ 1,m range. The biggest advantage of CPP is noted when the ice condition becomes more severe. The CPP design can still generate enough thrust (by changing its pitch through turning blade position) to propel the ship to overcome the large ice resistance. However, it should be noted that the strength of CPP is a serious design concern as the large ice impact can easily damage the CPP mechanism. Compared to the direct calculation for ice propulsion power determination, FSICR is quite conservative. Although FSICR did consider the differences of propeller designs through using different K c values for FPP and CPP (single, twin and triple propellers), it seems the Rule is still quite conservative for engine powering. Even more, FSICR cannot identify the limitation of a FPP design for high ice resistance condition. At high resistance, in some FPP designs such as the one selected, no matter how great the amount of engine power that can be provided to a propeller, the propeller still cannot generate enough thrust for ice propulsion. REFERENCES Bezinovic, D. 005, Design of an AFRAMAX tanker Baltic ice class IA acc. to new Swedish-Finnish rules 00 and model testing results basic hydrodynamics problem Marine Transportation and Explotation of Ocean and Coastal Resources, 005. Chen, H.C. and Lee, S.K., Time-Domain Simulation of Four-Quadrant Propeller Flows by a Chimera Moving Grid Approach, Proc. COE VI Conference, Baltimore, MD, Oct., 005. Chen, B. and Stern, F., Computational Fluid Dynamics of Four-Quardrant Marine- Porpulsor Flow, Journal of Ship Research, Vol. 4, No. 4, 1999. Finnish Maritime Administration and Swedish Maritime Administration, Guidelines for the Application of the Finnish-Swedish Ice Class Rule, 005. Finnish Maritime Administration, Finnish- Swedish Ice Class Rule, FMA Bulletin No., 00 (www.fma.fi). IACS, 008, URI Machinery requirement for Polar class ships 4 Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion

ABS TECHNICAL PAPERS 008 ITTC, Testing and Extrapolation Methods Resistance Resistance Test, ITTC Recommended Procedures, 00. ITTC, Performance, Propulsion 1978 ITTC Performance Prediction Method, ITTC Recommended Procedures, 1999. Kim, H.S., Ha, M.I., Ahn, D., Kim, S.H., Park, J.W., 005, Development of 115K tanker adopting Baltic Ice Class IA, the 15 th ISOPE, Seoul, Korea, June 16-4, 005. Lee, S.K., 007, Engineering practice on ice propeller strength assessment based on IACS polar ice rule URI, pp. 17-145, PRADS 007, Oct.1-5, 007, Houston, Texas, USA. Lee, S.K., 006, Rational approach to integrate the design of propulsion power and propeller strength for ice ships, IceTech 006, July 16-19, Banff, Canada. Lindqvist, G., A Straightforward Method for Calculation of Ice Resistance of Ships, Proc. 10 th Int. Conference on Port and Ocean Engineering under Arctic Condition (POAC), Lulea, Sweden, 1989. MAN B&W, S70MC Project Guide, 1994 MAN B&W, Propulsion Trends in Tankers, 00. MAN B&W 'Diesel Facts', pp 14-15, 006 Feb. Manen, J.D. van and Oossanen, P. van, Principles of Naval Architecture, Vol. II, Chapter 6, Propulsion, 1988. MARNET-CFD Thematic Network, Best Practice Guidelines for Marine Applications of Computation Fluid Dynamics, 00. Oosterveld, M.W.C and Ossannen, P. van, Further computer-analyzed Data of the Wageningen B-Screw Series, ISP,, July, 1975 Valanto, P., The Resistance of Ship in Level Ice, SNAME Transaction, vol. 109, 001 van Gent, W. and van Oossanen, P, 198, Ducted propeller systems and energy saving, NSMB publication no. 756a. Yazaki, A. and Sugai, N., 1966, Further model testa on four-blade controllable-pitch propellers, Ship Research Institute, Paper no. 16 Yazaki, A., 1964, Model testa on four-blade controllable-pitch propellers, Ship Research Institute, Paper no. 1 Combining Ice-Class Rules with Direct Calculations for Design of Arctic LNG Vessel Propulsion 5