Maneuvering predictions in the early design phase using CFD generated PMM data

Size: px
Start display at page:

Download "Maneuvering predictions in the early design phase using CFD generated PMM data"

Transcription

1 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 Maneuvering predictions in the early design phase using CFD generated PMM data Claus D. Simonsen, Janne F. Otzen, Christian Klimt, Nikolaj L. Larsen (FORCE Technology, Denmark) and Frederick Stern (IIHR, niversity of Iowa, SA) Abstract The present paper covers the work made in order to be able to perform standard deep water IMO manoeuvring simulations based on a combination of computed and measured hydrodynamic input data. Based on a full set of measured PMM data a reduced test matrix was identified and the standard and 0-0 zigzag and the 35 turning circle manoeuvres was simulated. Based on the reduced test matrix all the static PMM conditions were computed with the RANS code STAR-CCM+ in order to obtain the hydrodynamic forces and moments. The computed static PMM data was subsequently used to replace the corresponding measured static PMM data in the simulations. Comparison between the computed and measured forces and moments showed a quite good agreement. Formal verification and validation, taking numerical and experimental uncertainties into account, was made for a selected condition. The result showed that the X-force is validated, while the Y-force and the yaw moment are not. For the simulated manoeuvres, the results also look promising, since reasonable good agreement was found when comparing the simulated turning circle and zigzag manoeuvres, obtained with measured and computed input data. Introduction Today manoeuvring simulators are used to find the trajectory of the ship in arbitrary manoeuvres. During the ship design phase, focus is mainly on the requirements to the IMO standard manoeuvres (10/10 and 0/0 zigzag and 35 turning circle), which the new ship must comply with. The simulator models are based on integral quantities and the idea is to solve the equations of motion together with large measured sets of hydrodynamic force/moment data. The major advantage of the simulators compared to CFD is their ability to simulate general manoeuvres in real time. But, since they rely on usage of experimental PMM data, which is costly and time consuming to produce, they are usually used late in the design phase, after the final hull form is known and the scale model can be built. Consequently, the manoeuvring characteristics of the ship are evaluated late in the design process. sing CFD instead of the experiment for the data generation would not require the scale model and the simulation could therefore be done earlier in the design process. It is well known that CFD requires much CP time when it comes to producing hydrodynamic force/moment data by sweeping through large test matrices. But, taking into account how much extra useful information the CFD results provide about the flow field, this is acceptable. Maneuvers have traditionally been simulated by empirical methods, for instance as done in Agdrup (008) or by simulators based on mathematical models using experimental rotating arm or PMM data, for instance as done in Otzen (008). The application of CFD in connection with maneuvering has shown promising results in the past. However, much of the work performed so far has primarily focused on individual PMM conditions like static drift and static rudder configurations, e.g. Simonsen and Stern (003), Simonsen and Stern (005) and Simonsen et. al (006) or pure yaw, e.g. Simonsen and Stern (008). The SIMMAN008 manoeuvring workshop, showed the same tendency, except for very few cases where CFD was used to predict the final manoeuvres directly. Only one participant, Hochbaum et. al (008), tried the approach used in this present paper, i.e. CFD based PMM. The approach was somewhat simplified, since effects of free surface and dynamic sinkage and trim were not taken into account, but even so quite good results were obtained for the final manoeuvres. Based on above, the objective of the current work is to be able to perform standard deep water IMO maneuvering simulations in the early design phase purely based on computed input data, i.e. without using experimental data. This also includes investigations on how the numerical PMM test can be reduced to minimize the computational effort, while still

2 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 producing sufficient information to give reliable maneuvering simulations. Approach All the experimental and numerical activities in the present work are focused on the appended KCS containership, which, according to the 4th Maneuvering Committee ITTC Report, is appointed as one of the benchmark ships within the maneuvering community. A full standard test matrix for the PMM test is defined and the PMM test is conducted in FORCE Technology s towing tank in Denmark. Based on this full data set, a mathematical model for the simulator is made and the IMO maneuvers are simulated. A reduction of the full test matrix is then made by removing points from the matrix and corresponding maneuvering simulations are performed to study the influence on the maneuvers, which will be quantified by overshoot angles and turning circle diameters. For each of the simulations the resulting maneuvers will be compared to the initial reference case obtained with the full matrix to quantify the changes. Based on this work, the test matrix for the computational work will be chosen. With the reduced test matrix, CFD computations will be made for all the conditions and the resulting X and Y forces and the yaw moment for each condition will be compared with data from the measurement as a point to point comparison. The computations are performed with the RANS solver Star-CCM+. Finally, all the computed force and moment results will be included in the simulator model in order to rerun the maneuvering simulations based on the numerical data input. The simulated maneuvers will be compared with the results from the previous simulations based on the measured data with full and reduced test matrix. The structure of the paper is as follows. First the model test is presented. This is followed by a presentation of the CFD work. In this connection, point to point comparison of all measured and computed quantities will be made. Next part covers the maneuvering simulations, i.e. the applied mathematical model and simulations based on measured PMM data and the reduction of test matrix. Finally, simulations will be made where the EFD data from the static conditions in the reduced test matrix are replaced with data computed with CFD. The result based on EFD and CFD will be compared to conclude about how well the CFD based input data performs. Data generation based on model testing The model testing is carried out in FORCE Technology s towing tank with a scale model of the KCS container ship using a scale of 1: The towing tank is 40m long, 1m wide and 5.5m deep. Figure 1. PMM model test setup at FORCE Technology. The main particulars and propeller data for the ship are shown in Table 1. The approach speed for the manoeuvres is 0 =1.701 m/s (4.0 knots full scale) corresponding to a Froude number of 0.6. The test is described in detail in Larsen (01), so only a brief summary is given here. Table 1: Hull and propeller data. L pp [m] B [m] T [m] S [m ] C B [-] D P [m] Z [no. of blades] 5 P/D 0.7 [-] The considered experimental test conditions cover a set of tests, which are representative for a 1 st quadrant 4 DOF PMM test and which can be used for assessment of the experimental uncertainty for representative conditions. The PMM testing technique enables various test conditions to be studied individually. The conditions, which are considered in this work, are static rudder, static drift, static drift and rudder, static heel static heel and drift, pure sway, pure yaw yaw and rudder, yaw and drift and yaw and heel. In the static test the model is towed in the same steady condition through the tank. The other tests are dynamic, i.e. the model is oscillating harmonically. In the present paper focus is placed on 3 DOF simulations, which means that only a subset of the conditions is covered, leaving all quantities related to heel and roll out. The contents of the considered test tests can be summarized as follows:

3 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 Static rudder : The model travels straight ahead through the tank with the rudder in a given angle δ. Static drift : The model travels through the tank in oblique flow due to a given drift angle β. Static drift and rudder : Same as Static drift but the rudder is deflected. Pure sway : The model travels through the tank on straight ahead course while it is oscillated from side to side. With u, v and r being the surge velocity, the sway velocity and the yaw rate in the ships local coordinate system, the pure sway motion can also be expressed in terms of the velocities, i.e. u= c (carriage speed), r=0 and v oscillates harmonically. Pure yaw : The model travels through the tank while it performs a pure yaw motion, where it is forced to follow the tangent of the oscillating path. In terms of velocities this means that v=0, while r and u oscillate harmonically. u oscillates, since the carriage speed in the present set-up is constant. Yaw and rudder : Same as Pure yaw but the rudder is deflected. Yaw and drift : The model travels through the tank, while it performs a pure yaw motion as described in Pure yaw. However, a fixed and preset drift angle is overlaid on the motion in order to obtain a drift angle relative to the tangent of the oscillating path. In terms of velocities this means that v 0, but constant, while r and u oscillate harmonically. For all of the above conditions, the tests were conducted according to FORCE s standard PMM testing procedures. This means that the model was constrained in all degrees of freedom except for heave and pitch to account for dynamic sinkage and trim. The PMM setup is shown in Figure 1. During the test the model is equipped with propeller and rudder. The propeller RPMs are kept constant at the model s self propulsion point. The model is shown in Figure and Figure 3, where the sand strips used for turbulence stimulation can also be seen. The PMM test program is carried out for speeds down to corresponding to a Froude number of The test program carried out including the repeat tests for the uncertainty assessment are shown in Table. Table : PMM test programme for the design condition. 0 is nominal approach speed and N 0 is corresponding impeller revolutions. Speed [1/ 0 ] Revs. [1/N 0 ] Drift angle [deg] Rudder angle [deg] Sway velocity [-] Yaw velocity [-] Standard program Static speed Static drift Static rudder Drift & rudder Pure sway Pure yaw Yaw & drift Yaw & rudder Repeated tests (no. of repeats) Static drift (13,1) 0.78, Static rudder (1) Pure sway (4) Pure yaw (1)

4 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 FX measured Tm X (1) 0.5ρ L pp FY measured Tm Y () 0.5ρ L pp M Z measured TmL pp N (3) 0.5ρ Figure. Bow of KCS model with sand strips for turbulence stimulation while in the dynamic cases the measured quantities include inertial contributions due to acceleration of the ship. In this case the inertial contributions are subtracted from the measured quantise to get the hydrodynamic part as FX measured M ( u rv X r Y r ) X (4) 0.5ρ T L m pp FY measured M ( v ru Y r X r ) Y (5) 0.5ρ T L m pp N M Z measured I r M ( X Z 0.5ρ ( v ru) Y Tm L pp ( u rv)) (6) Figure 3. Stern of KCS model with rudder and propeller. During the test the instantaneous operating conditions for the ship like speeds, positions etc. are measured together with the resultant forces. All forces are measured in a coordinate system following the ship, meaning that X-components act in the longitudinal direction of the ship (positive forward) and Y- components perpendicular to this direction (positive starboard). The yaw moment is taken with respect to the mid-ship position at /. All hydrodynamic forces and moments presented in the present work are non-dimensionalized by the data reduction equations shown below. It should be mentioned that for the static conditions, the hydro dynamic forces and moments are equal to the measured quantities, i.e., and are the longitudinal force, the transverse force and the yaw moment, respectively. is the water density. is the ship speed, where u and v are the surge and sway velocities respectively. r is the yaw rate. Finally, the dots above the velocity quantities indicate the corresponding accelerations. See Figure 4. and are the mean draft and the length between perpendiculars. and are the mass and moment of inertia of the model, i.e. of the model itself, the gauges and the ballast weights. and are the X- and Y- distances from the center of gravity of the model to the point, which the model rotates around. Figure 4. Definition of state variables.

5 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 A part of the experimental work also covered assessment of the experimental uncertainty. Following the approach in ITTC 1999a and b and Simonsen 004 the uncertainty assessment, which covers both precision and bias limits, is based on the data reduction equations for the forces and moments listed above. The uncertainties are expressed as X Y N X X B P (7) Y Y B P (8) N N B P (9) B and P are the bias (systematic errors) and precision (random errors) limits, respectively. The bias contributions are found by sensitivity analysis of the data reduction equations and error estimates related to the individual components of the measurement system, The bias limits are assessed based on a study of the measuring system. According to ITTC 1999a they can be estimated on the basis of B r J i1 i i J 1 J B B (10) i k i1 k i1 where is the influence coefficient defined by r i (11) X i is the bias limits in and is the correlated bias limits in and. B ik L 1 ( B ) ( B ) (1) i k where is the number of correlated bias error sources that are common for measurement of variables and. The bias error for each variable in the data reduction equation may consist of a number of bias errors, so in order to calculate the combined bias error the root-sum-square is used J B i ( B i ) k (13) k 1 is the number of the considered variable in the data reduction equation. ik The precision limits are assessed through repeated tests, which are built into the test program. The model has not been dismounted from the carriage during the test, so in order to disturb the system the repeat tests has been mixed with the other test configurations. According to ITTC 1999a the precision limit is estimated from Sr Pr (14) M where is the number of repeats and the factor of is applied for 10. is the standard deviation defined as ½ M ( r ) k r S r (15) M 1 k 1 Here is the value from each repeat test and is the mean value of all the quantities from the repeat tests. is defined as 1 r M M r k k 1 Table 3. ncertainties for static drift, β=8, Fr=0.01. (16) X Y N % The uncertainty assessment procedure is very extensive and is not described here. It has only been done for a few conditions as spot checks of the uncertainty level. It is described in details in Larsen (01). Table 3 and Table 4 show the results of the uncertainty assessment for the static drift case with 8, and static rudder case with 10, respectively at model speed of 1.318m/s, i.e Table 4. ncertainties for static rudder, δ=10, Fr=0.01. X Y N %

6 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 The level of the uncertainties is higher than for a resistance test, but this is to be expected since the PMM test is conducted with a two-gauge system, which typically gives higher uncertainties than a single gauge system. The uncertainty level seems reasonable and the level is comparable to other to PMM tests in Benedetti et al. (006), Yoon (009) and Simonsen (004). With respect to the measured data, examples are shown in Figure 8 to Figure 19 in the section below, in connection with the discussion of the comparison between computed and measured PMM data. Data generation based on CFD computations The CFD computations in the present work are supposed to replace both dynamic and static PMM conditions. However, at present time of writing the dynamic conditions are not completed, so only results for static cases are covered. Computational method The computations are performed with the Reynolds Averaged Navier-Stokes (RANS) solver STAR-CCM+ from CD-adapco. The code solves the RANS and continuity equations on integral form on an unstructured mesh by means of the finite volume technique. Both steady state and transient calculations are considered depending on the considered test condition. For the steady state calculations the temporal discretization is based on a first order Euler difference, while a second order difference is used for transient calculations. Spatial discretization is performed with second order schemes for both convective and viscous terms. The pressure and the velocities are coupled by means of the SIMPLE method. Closure of the Reynolds stress problem is achieved by means of the isotropic blended k-ε/k-ω SST turbulence model with an all Y+ wall treatment, which depending on the Y + value selects the near wall treatment. The rotating propeller is not modelled directly, as this requires very long computational times that are not realistic when computing a large number of cases. Therefore, a prescribed body force propeller is applied instead as is commonly done for this type of simulations. The propeller model is based on a radial varying body force field, which follows the variation of a theoretically derived circulation distribution. It is prescribed by means of the ship speed based advance coefficient and the propeller coefficients and. The model is integrated in the CFD-code, where the axial and tangential force components are given by fb X A X r r 1 r 1 r, fbθ Aθ (17) (1 r ) r r Where the non-dimensional propeller and hub radii are defined by r ( r rh )/ (1 rh ) and rh Rh /Rp. With being the longitudinal extension of the disk divided by the ship length, the coefficients and are defined by CTh Ax Δx 16 A θ V r 1 r h h h h h h (18) KQ 105 (19) ΔxJ π4 3r 1 r The free surface is modelled with the two phase volume of fluid technique (VOF). Further details about the code can be found in STAR-CCM+ user guide (01). Finally, it must be noted that all computations are done in model scale with a model similar to the one used in the model test. The motions of the ship are described in terms of an earth-fixed inertial reference coordinate system and a ship fixed coordinate system. The motion of the ship, i.e. translations and rotations, are described with respect to the inertial frame. The linear and angular velocities plus the forces and moments are described with respect to the ship-fixed frame. In the present calculation, three degrees of freedom are considered and the ship motions in surge, sway and yaw are all prescribed. The motions are applied in the earth-fixed frame by mowing the complete computational domain and hereby introduce grid velocities according to the considered motion. This applies to both static conditions like static drift where the ship is moved with constant speed to obtain the desired drift angle and to dynamic tests like pure yaw motion, which are not covered in this paper. Heave and pitch motions in terms of dynamic sinkage and trim are taken from the model test, so they are not predicted in the CFD solver. The thrust and torque values applied in the prescribed bodyforce model are also taken from the model test. Computational domain The computational mesh is based on the trimmed mesh approach, which means that the mesh is dominated by hexahedral cells in most regions except close to the hull where the hull surface is used to trim the hexahedral cells leaving a set of polyhedral cells. Further, prism layer mesh is used near the hull surface to resolve the boundary layer flow. Concerning boundary conditions on the domain, an inlet condition

7 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 prescribing the inlet velocity and the still water level for the volume fraction is given up-stream of the model and on the sides of the domain. On the outlet boundary a pressure condition is applied in order to have zero gradients for the velocity and volume fraction and prescribed hydrostatic pressure. On the top and bottom of the domain, slip walls are applied so there is no flow through the boundary. The inlet boundary is located 4.0 ship lengths from the ship, Outlet and sides are located 6.4 and 5.0 ship lengths from the ship, respectively. Figure 7. Mesh in bow region. Verification and validation Figure 5. Illustration of domain boundaries. Figure 6 and Figure 7 show the computational mesh in the bow and stern regions, respectively. Refinement zones are applied in the free surface region and on the hull where the wave formation is expected to take place. Figure 6. Mesh in bow region. Concerning verification focus is on mesh sensitivity. To study the sensitivity of the mesh size on the results one condition has been picked to be investigated. Three different grid sizes have been used and the corresponding forces have been evaluated for a static drift condition with 0 degrees drift angle at model speed of 1.318m/s, i.e Further, no running body-force propeller is included in this case. With the unstructured mesh, it is not possible to make a completely systematic grid refinement, since it is not possible to directly control the size of the individual cells. However, all mesh quantities are given as a percentage of a base size, so in order to change the grid as systematically as possible the coarse (No. 3), medium (No. ) and fine (No. 1) grids were made by changing the base size with a constant factor. Based on the three grids, the CFD simulations were conducted and the longitudinal X-force, X, the transverse Y-force, Y, and the yaw moments, N, taken around the Z-axis at mid-ship position were computed. It is difficult to estimate the grid uncertainty for unstructured meshes, so in lieu of an uncertainty method for unstructured meshes and since the present grid is hexahedral dominated, it has been attempted to estimate the grid uncertainty based on Richardson extrapolation. Solutions were obtained on 3 systematically refined grids by means of the refinement ratio,. This gave mesh sizes ranging from approximately 1.5 to 7.0 million cells. Exact grid numbers and computed force and moment coefficients can be seen in Table 5. Forces and moments are nondimensionalized with ½ and ½ respectively.

8 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 Table 5: Non-dim forces and moments for different grids. No. of cells X Y N 1,491, ,068, ,969, The grid uncertainty is estimated according to the procedure given in Stern et al. 001 and Wilson et al Briefly, the procedure is as follows. The changes in solutions between coarse and medium grids, S 3 3 S, and between medium and fine grids, S 1 S1, are used to calculate the convergence ratio R / 1 3. Depending on the value of three conditions can occur: i) 0 1, grid convergence, ii) 0, oscillatory convergence and iii) 1, grid divergence. In condition iii) no uncertainty can be estimated. In condition ii) the uncertainty is estimated by ½( S SL ) (0) where S and S L are the maximum and minimum of the solutions from the considered grids. In condition i) grid convergence occurs and generalized Richardson extrapolation (RE) is used to estimate the grid error (1) RE and the order of accuracy p 1, which are given as (1) 1 RE (1) 1 p r 1 and p ln / 3 1 () ln( r ) r x / x1 x3 / x is the refinement factor, which defines the relation between the considered grids. When and p are know it is (1) RE 1 possible to estimate the grid uncertainty. There are two ways to do this depending on whether the solutions are close to the asymptotic range or not. If the correction factor defined by r 1 C (3) pest r 1 p where p est is the limiting or theoretical accuracy of the applied numerical method, is close to unity, the solutions are close to the asymptotic range. In this case the sign of the error is known, so the numerical error, SN, benchmark, S C and uncertainty, C can be calculated from SN (1) C (4) S C RE 1 SN S (5) (.4(1 C ) 0.1) RE, 1 C C 1 C RE, 1 C (6) If the correction factor is away from unity only the numerical uncertainty is calculated (9.6(1 C ) 1.1) ( 1 C 1) RE 1 RE1,, 1 C 1 C (7) Table 6: % change in non-dim forces and moments between grids. % Change X Y N 3 to to Table 7: Estimated order of accuracy and grid uncertainty. Quantity p % S 1 X Oscill Y N Oscill Table 6 shows the change in quantities between grids. It is seen that all quantities are converging as the grid is refined, but both X and N show oscillatory convergence, i.e. the order of accuracy cannot be estimated. For Y Table 7 shows that an order of accuracy of.3 is achieved. This is close to the theoretical order of accuracy of. Finally, Table 7 also shows the estimated grid uncertainty, which ranges from 0.6 to 1.4 % of the fine grid solution. It should be noted that if the solution for Y is corrected with the

9 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 numerical error, the grid uncertainty can be reduced to.3%. In order to validate the computations against measurements, the measured data shown in Table 8 are used. No specific uncertainty assessment was made for the un-propelled 0 degrees static drift case, so the uncertainty estimates from the propelled 8 degrees drift case is applied, by assuming that the percentage in Table 3 is representative for the un-propelled case. Table 8: Measured values and data uncertainty Quantity D D D % D X Y N Quantity Table 9: Validation investigation. E=D-S V E % D X Y N To finally check if Validation has been obtained it is necessary to check if the absolute comparison error E is smaller than the validation uncertainty calculated as V D, i.e. if the comparison error is within the combined noise from the CFD simulation and the experiments. From Table 9 it is seen that the computation agrees with the measurement within approximately maximum 5%, so the agreement is pretty good taking into account that it is a pretty complex flow field at 0 degrees drift angle. However, in spite of this it is only X that is validated. Taking into account that going from medium to fine grids will only change the results with approximately % while the computational time will be more than doubled, it was judged that the medium mesh was sufficient for the present study, which has much focus on computation of large data sets for practical application. It should also be mentioned that later in the project it was necessary to increase the mesh size a little in order to have a mesh that work for non-zero rudder angles. Therefore, the final mesh used for the production runs contains app million cells, which is a little finer than the initial medium mesh. It was designed to give a Y + of maximum 75.5 and an average over the hull of 3.5. After the verification and validation work, all the production runs for the static maneuvers to be used in the simulator model were made. The figures below show comparisons between CFD and EFD for all the computed cases, which cover static drift, static rudder and the combined static drift and rudder conditions. As mentioned earlier only static CFD data will be used as input for the simulator, as the dynamic computations are not finalized at the present time of writing. Figure 8 to Figure 10 show the computed and measured non-dimensionalized forces and moments, X, Y and N as function of the rudder angle for three different speeds. For X it is seen to be close to 0 and negative for the highest speed, whereas it increases to larger positive values as the speed is reduced. This means that the propeller is over-thrusting at the lower speeds. However, based on the applied constant RPM approach, where the propeller RPMs are kept constant at service speed self-propulsion point, the propeller will deliver too much thrust at the lower speeds giving the observed behavior. It is seen that the CFD computations generally predicts X fairly well. The slightly larger deviations at large negative rudder angles are probably related to the simplified propeller model not being able to capture the rudder propeller interaction. X CFD, Fr=0.60 EFD, Fr=0.60 CFD, Fr=0.01 EFD, Fr=0.01 CFD, Fr=0.156 EFD, Fr=0.156 Figure 8. Computed and measured non-dim X force for static rudder. Y CFD, Fr=0.60 EFD, Fr=0.60 CFD, Fr=0.01 EFD, Fr=0.01 CFD, Fr=0.156 EFD, Fr=0.156 Figure 9. Computed and measured non-dim Y force for static rudder. For Y good agreement between EFD and CFD is also observed. Increased rudder angle leads to increased Y-force which makes sense since the rudder lifts more at higher rudder angles. It is seen that the slope of the curves increases as the speed is reduced and the weak non-linearity occurs at the larger rudder angles.

10 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 N CFD, Fr=0.60 EFD, Fr=0.60 CFD, Fr=0.01 EFD, Fr=0.01 CFD, Fr=0.156 EFD, Fr= Y Drift angle [deg] CFD, Fr=0.60 EFD, Fr=0.60 CFD, Fr=0.01 EFD, Fr=0.01 CFD, Fr=0.156 EFD, Fr=0.156 Figure 10. Computed and measured non-dim yaw moment for static rudder. N also shows increased slope of the curves as the speed is reduced and the weak non-linearity at the larger rudder angles. This behavior is expected as the yaw moment for the static rudder case to a large degree is caused by the rudder lift force. Again good agreement between EFD and CFD is observed. Figure 11 to Figure 13 show the computed and measured non-dimensionalized forces and moments, X, Y and N as function of the drift angle for three different speeds. Again X shows the effect of the constant RPM approach discussed above. Further, it can be noticed that the largest deviations occur at the lowest speed, which presumably is caused by the bodyforce propeller not being able to model the propeller flow in the case where the propeller loading is increasing. However, in spite of the difference, the CFD code generally predicts the force fairly well compared to the measurement, plus it is able to capture variation in X with the drift angle. X CFD, Fr=0.60 EFD, Fr=0.60 CFD, Fr=0.01 EFD, Fr=0.01 CFD, Fr=0.156 EFD, Fr=0.156 Figure 11. Computed and measured non-dim X-force for static drift. Figure 1. Computed and measured non-dim Y-force for static drift Drift angle [deg] N CFD, Fr=0.60 EFD, Fr=0.60 CFD, Fr=0.01 EFD, Fr=0.01 CFD, Fr=0.156 EFD, Fr=0.156 Figure 13. Computed and measured non-dim yaw moment for static drift. Figure 1 and Figure 13 show Y and N for the three speeds. Compared to the static rudder case it is here seen that the non-dimensionalization basically collapses the force and moment curves into one curve, which is often seen for the static drift force and moment. It is also seen that as the drift angle is increased the Y-force becomes non-linear. The yaw moment also shows non-linearity, but not as pronounced as the Y-force. Further, it seems that the non-linearity in N is more pronounced for lower speeds. Concerning the comparison between model test and CFD results, it appears the computation agrees well with the measurement. The final set of results to be shown covers the static drift and rudder results. Figure 14 to Figure 16 show the forces and moments as function of the rudder angle for the two negative drift angles β=-4 and β=- 1, while Figure 17 to Figure 19 show the forces and moments as function of the rudder angle for two corresponding positive drift angles β=4 and β=1.

11 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August X CFD, Fr=0.01, β= 4 EFD, Fr=0.01, β= 4 CFD, Fr=0.156, β= 1 EFD, Fr=0.156, β= X CFD, Fr=0.01, β=+4 EFD, Fr=0.01, β=+4 CFD, Fr=0.156, β=+1 EFD, Fr=0.156, β=+1 Figure 14. Computed and measured non-dim X-force for static drift and rudder. Negative drift angles Figure 17. Computed and measured non-dim X-force for static drift and rudder. Positive drift angles Y CFD, Fr=0.01, β= 4 EFD, Fr=0.01, β= 4 CFD, Fr=0.156, β= 1 EFD, Fr=0.156, β= Y CFD, Fr=0.01, β=+4 EFD, Fr=0.01, β=+4 CFD, Fr=0.156, β=+1 EFD, Fr=0.156, β= Figure 15. Computed and measured non-dim Y-force for static drift and rudder. Negative drift angles N CFD, Fr=0.01, β= 4 EFD, Fr=0.01, β= 4 CFD, Fr=0.156, β= 1 EFD, Fr=0.156, β= 1 Figure 16. Computed and measured non-dim yaw moment for static drift and rudder. Negative drift angles. Only the two lower speeds are covered in this case. The reason is that the test program reflects the conditions the ship will experience during maneuvering and typically the combined drift and rudder angle condition occurs later in the maneuver, where the ship is experiencing a speed loss. Figure 18. Computed and measured non-dim Y-force for static drift and rudder. Positive drift angles. For all the forces and moments the overall agreement between the measured and computed PMM results is fairly good. Though, it generally seems that it is most difficult to capture the measured data at the lowest speed. One of the main reasons for this is believed to be the propeller model. The rudder load is strongly dependent of the flow over the propeller since it is located in the propeller slip stream. On the other hand the propeller flow is strongly dependent of the wake field behind the ship as this generates the propeller inflow field. The simplified propeller does not really reflect the difference in propeller load distribution over the disk occurring from the change in the wake field that occurs during drift conditions and reduced speed due to the prescribed body-force approach. If the load distribution over the propeller disk is not correct, the distribution of the propeller slipstream velocity that hits the rudder will not be correct, which again means that the rudder performance will be influenced. This again will influence X, Y and N and probably lead to larger deviations between EFD and CFD, particularly for lower speeds and large drift angles, which also change the wake flow significantly.

12 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August N CFD, Fr=0.01, β=+4 EFD, Fr=0.01, β=+4 CFD, Fr=0.156, β=+1 EFD, Fr=0.156, β=+1 Figure 19. Computed and measured non-dim yaw moment for static drift and rudder. Positive drift angles. Based on the above comparison the idea of generating PMM data on the basis of CFD looks promising. However, the question is if the small differences between EFD and CFD will influence the final maneuvers after applying the data as input to the maneuvering simulator. This is investigated in the following sections where simulations with both EFD and CFD data will be compared. Mathematical model applied for the manoeuvring simulations A mathematical model was generated to determine a specific set of X, Y and N derivative with the most optimal fit to the measured forces of the KCS. The analysis was performed in a classic way where initially the linear Y and N coefficients, Y v and N v were determined from the pure drift test at nominal speed with drift angle in the linear range between +/-4 degrees. These two coefficients were then fixed and included in the following general static analysis where the total matrix of static test data is applied to solve the total matrix of static coefficients, i.e. the remaining static derivatives are determined in a single combined analysis. The two added mass coefficients, Y v and N v was determined based on a special pure sway analysis which comprised the time series from all the dynamic pure sway tests. The general dynamic analysis was performed similar to the static analysis with a combined analysis for all of the dynamic test types. In the dynamic analysis all the static coefficients were fixed to the values adopted from the static analysis. The final coefficient formulation applied for the simulation of manoeuvres were output from the dynamic analysis (including the static analysis) supplemented with the two added mass coefficients from the pure sway analysis and an empirical value for the added mass coefficient, which was calculated as 6 percent of the non-dimensional mass of the vessel. I.e. based on the model tests the final coefficient formulation becomes: X X X X X X Y Y Y Y o rr r X r X vv Y o Y Y Y Y v u v X r u X vr 4 X Y vvv u rr N N N N N o N N N N rrr vrr vv u v r X 3 v Y v Y r v r u u u Y 3 r Y u Y r Y N vvv rr v rrr vrr u vv uu u X u u X X X v v Y u r v r uv 3 v N 3 r N v N r v v Y Y Y u v u X v r Y Y r Y v X X u X v r u v Y r Y v r u N v u r u N v N r N r uv v r vv vvr u v Y v N N N u r N r v r N N vv r X Y rr r r v X r v v Y v rr Y v r Y 3 u Y v Y v u v N r N v r v u r v r vv vvr u v N r r r v v N rr r v v 3 v N r r v r Y vr N v r N 3 u N v N v v r r r v v v 3 r r v r N vr v r v r r v r v r r v r v r r v r (8) (9) (30) In case of a 4DOF model a corresponding set of derivatives for the roll (K) moment would have been generated. However, this has been omitted due to the focus on 3DOF. Since the vessel is a single screw ship it possesses inherent asymmetry regarding port starboard manoeuvring ability. This asymmetry is reflected in a number of X, Y and N coefficients. The asymmetry is partly due to the propeller side force and partly due to unavoidable asymmetry in the physical model used for the model tests. To account for the scewness of the model, the Y and N forces acting on the model in a bare hull resistance test at drift angle zero is subtracted from

13 those measured in the PMM series with the propeller working. The generated model is strictly speaking only applicable for the given approach speed. In case of different approach speeds the model ought to be reevaluated to find the corresponding coefficients. Test matrix reduction and simulation results based on EFD For standard PMM tests a full standard test program (STP) similar to the programme presented in Table is usually applied. When running just tank tests there has been no significant motivation for minimizing the test programme since the main costs are associated with the initial setup of the model and there has rather been a tendency to utilize the length of the entire towing tank to pick up as many points as possible. This approach is suitable for tank tests, but in connection with numerical model tests where the calculation time is the most expensive cost, there is a need to minimize the number of points/calculations. To evaluate a minimum test programme (MTP), a set of simulated manoeuvres based on MTP were compared to similar manoeuvres based on the STP. To establish whether the MTP was acceptable or not a set of criteria or tolerances were defined. If the initially tested MTP did not meet the defined tolerances the number of measuring points was to be increased until the limits were fulfilled. The predefined acceptable tolerances between STP and MTP are presented in Table 10. In order to meet the criteria a number of iterations were performed increasing the number of data-points in the MTP. Different approaches were tested for example adding points mainly focusing on static points and correspondingly adding points mainly focusing on 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 dynamic points. The final minimum test program comprised the total matrix of dynamic test points from the STP while the static part only contained 58 data points of the total 8 in the STP. In fact the final MTP did not meet all the tolerances, but it was accepted as close enough and robust enough to fulfill the purpose. The criteria values, i.e. differences between final MTP and STP maneuvering parameter, are listed in Table 11. The final minimum test matrix is presented in Table 1. Table 10: Acceptable tolerances between results of MTP and STP simulated manoeuvres, where zz is zigzag, OS is overshoot angle, ITA is initial turning ability and tc is turning circle. 10/10 zz 1st OS [deg] +/ /10 zz nd OS [deg] 0.8 ITA [L pp ] 0.0 0/0 zz 1st OS [deg] 0.8 0/0 zz nd OS [deg] tc advance [L pp ] tc tactical diameter [L pp ] tc steady drift [deg] 1.0 Table 11: Difference between results of MTP and STP simulated manoeuvres. Port Starboard 10/10 zz 1st OS [deg] /10 zz nd OS [deg] ITA [L pp ] /0 zz 1st OS [deg] /0 zz nd OS [deg] tc advance [L pp ] tc tactical diameter [L pp ] tc steady drift [deg] Table 1. Minimum PMM test programme for the design condition. The dynamic matrix is not included as this is similar to the STP. Speed [1/ 0 ] Revs. [1/N 0 ] Drift angle [deg] Rudder angle [deg] Static speed 0.6, 0.78, Static drift , -, -0.5, 1, , -4, 8, , -8, -4, 0, 4, 8, 1, 16 0 Static rudder , -0, 10, , -10, 0, , -10, 0, 35 Drift & rudder , -10, 0, 10, 0, 30, , -30, -0, -10, 0, 10, , -10, 0, 10, 0, 30, , -30, -0, -10, 0, 10, 0

14 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 Simulation results based on CFD data As the next step in line towards a full model test matrix of pure numerical model tests the static part of the found MTP was generated using CFD. The purpose was to substitute the static derivative found from the model tests with static derivatives evaluated from numerical model tests. The numerical data was evaluated just as if they have been recorded in the towing tank, i.e. the calculated hydrodynamic forces are faired by the same combination of independent variables as presented in equation 8 to 30. The dynamic coefficients are not reevaluated based on the new static coefficients; to make it simple the static coefficients from the MTP is just replace by the new CFD static coefficients. The predicted static hydrodynamic derivatives are presented in Table 13. Simulations were made using both sets of hydrodynamic derivatives presented in Table 11. The main values for 35/-35 degree turning circles and 10/10 and 0/0 starboard and port zigzag manoeuvres are presented in Table 14 and Table 15. The simulation results are plotted in Figure 0 to Figure 3. Table 13: Predicted hydrodynamic derivatives based (multiplied by 10 5 ) 10 5 MTP CFD X X X V X VV X D X DD X D X DD X VD X VVD X VDD Y Y V Y VV Y VIVI Y V Y D Y DD Y DDD Y D Y DD Y DDD Y DV Y DVV Y DIVI Y VDD N VV N VIVI N V N D N DD N DDD N D N DD N DDD N DV N DVV N DIVI N VDD Table 14: MTP simulated manoeuvres. Port Starboard 10/10 zz 1st OS [deg] /10 zz nd OS [deg] ITA [L pp ] /0 zz 1st OS [deg] /0 zz nd OS [deg] tc advance [L pp ] tc tactical diameter [L pp ] tc steady drift [deg] Table 15: CFD simulated manoeuvres. Port Starboard 10/10 zz 1st OS [deg] /10 zz nd OS [deg] ITA [L pp ] /0 zz 1st OS [deg] /0 zz nd OS [deg] tc advance [L pp ] tc tactical diameter [L pp ] tc steady drift [deg] Advance [m] Transfer [m] Figure degree turning circle to port and starboard, where blue is the model with static CFD and the red is pure MTP. N N V

15 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 Speed [kn] Time [s] Heading and rudder angle [deg] Time [s] Drift angle [deg] Time [s] Heading and rudder angle [deg] Time [s] Figure 1. Speed and drift angle in a 35/-35 degree turning circle, where blue is the model with static CFD, red is pure MTP, dashed is -35 deg and solid is 35 deg. Heading and rudder angle [deg] Heading and rudder angle [deg] Time [s] Time [s] Figure. 10/10 port and starboard zigzag, where blue is the model with static CFD and the red is pure MTP. Figure 3. 0/0 port and starboard zigzag, where blue is the model with static CFD and the red is pure MTP. The simulations based on the data from CFD indicate a more course stable ship causing a longer initial turning ability and larger turning circles. The differences are apparently most severe for the manoeuvres to port side, which might be due to the simple propeller model applied for the CFD calculations. Considering the hydrodynamic derivatives presented in Table 11 these also reveals that the largest deviations are seen among the derivatives dependent on the rudder. Appling a simple propeller model will affect the flow field around the rudder due to the simplification of the flow field from the propeller. Instead of applying a real propeller, the propeller is modeled using an actuator disc. For normal CFD calculations (resistance/self propulsion) this approach seems reasonable since scewness of the model due to the propeller is of minor interest, but when simulation manoeuvres at high drift and rudder angles the port/starboard asymmetry and correspond change of flow field down to the rudder appears to have a larger influence on the manoeuvres than initially assumed as discussed earlier in the CFD result section. To investigate the influence of asymmetry simulations were made for the CFD predicted model using modified Y 0 and N 0, applying at first the values from the original pure experimental MTP model and

16 9th Symposium on Naval Hydrodynamics othenburg, Sweden, 6-31 August 01 afterwards the values from the static CFD calculations. It turned out that when applying the values from pure static numeric calculation, the port starboard asymmetry leveled out even more than previously revealing for example an advance/l pp of 3.9 to port and 3.37 to starboard against the previous 3. and 3.45 respectively, i.e. the model is quite sensitive to the applied values of Y 0 and N 0. Figure 1 showing the speed and the drift angle during the turning circles also clearly shows that the asymmetry for the CFD model is not as pronounced as for the MTP model. This could be a consequence of the propeller side force, which is missing in the present propeller model. Concluding remarks In the present work the KCS container ship is investigated during maneuvering by means of EFD and CFD in order to see if CFD results can be used as input for the maneuvering simulator instead of measured data, as this will allow pure numerical evaluation of the ships maneuvering characteristics in the early design phase where no physical scale model is available. The experimental study was conducted in FORCE Technology s towing tank in Denmark and the CFD study including verification was conducted using the RANS code STAR-CCM+. The experiment covered a 4DOF PMM test, but only 3DOF results were used to simulate the IMO standard maneuvers in this paper. The results from the model test cover the hydrodynamic X-force, Y-force and yaw moment acting on the fully appended ship during the captive PMM motion including experimental uncertainty for selected conditions. Based on a mathematical model of the maneuvering ship, the results from the full test program were first used to simulate the standard IMO maneuvers, i.e. turning circle and and 0-0 zigzag. Afterwards the simulation results were used as basis for a reduction of the test matrix in order to reduce the amount of input data required to get as close as possible to the initial maneuvers. The idea of reducing the test matrix is to reduce the number of CFD computations required to perform the maneuvering simulations with as reasonable accuracy. After reducing the test matrix, all the static conditions in the matrix were computed with the CFD code and the results were compared with the measured data. After ensuring that the results looked reasonable the CFD data were analyzed and fed into the maneuvering simulator by replacing the measured values in the mathematical model and the standard maneuvers were simulated again to see how well the simulated maneuvers based on the reduced experimental test matrix compared with the maneuvers obtained with the CFD based input data. With respect to the results obtained, the experimental work showed that the experimental uncertainty of the forces and moment for the selected conditions all were within 5.5%, which is comparable to the level obtained in previous PMM measurements with the 5415 frigate and the KVLCC tanker. Based on the comparison between standard IMO maneuvers obtained with the full test matrix and a gradually reduced matrix it was possible to reduce the number of static PMM conditions with approximately 30%, while still obtaining reasonable simulation results. It should be noted that the reduction involved experience of the test engineer and it was not done by systematically removing points. The systematic approach is difficult to do since the maneuvering model easily becomes unstable and diverges. Also, it cannot be concluded that removing 30% of the points for all ships will work, as this was only done for one ship. Concerning the CFD work it was initiated with a verification of the numerical model. It was assumed that the numerical error primarily comes from the discretization so a grid study on 3 systematically refined grids was made. The grid uncertainty was estimated assuming that the hexahedral dominated mesh could be refined systematically in spite of being unstructured in some regions. The verification showed grid uncertainties within.5%, which is reasonable. However when doing the validations against measured data, taking the experimental uncertainty into account, the comparison errors became larger than the validation uncertainty leading to validation of X alone. However, generally the difference between measurement and computations was quite small giving comparison errors of 5% of the measured value. After computing all the conditions in the reduced test matrix, comparison between measurement and computations was made for all data points and as indicated initially by the verification conditions quite good agreement was generally obtained. There are conditions at higher drift and rudder angles, where differences start to occur, but this may be expected as the flow starts to separate and therefore is more difficult to compute accurately. Further, the simplified and prescribed propeller model also has limitations, as it cannot capture the change in propeller load distribution over the disk as the inflow changes. This for instance happen at large drift angles. Consequently, the rudder, which is located in the propeller slip stream, experiences an inflow that is not correct leading to less accurate rudder loads. However, in spite of the deviations the forces and moments predicted by CFD look promising. The maneuvering simulations shown in the paper were initially made for the full standard test matrix in order to find reference values for IMO maneuvers. Subsequent reduction of the test matrix showed that the

VERIFICATION AND VALIDATION OF RESISTANCE AND PROPULSION COMPUTATION

VERIFICATION AND VALIDATION OF RESISTANCE AND PROPULSION COMPUTATION VERIFICATION AND VALIDATION OF RESISTANCE AND PROPULSION COMPUTATION G. Deng, A. Leroyer, E. Guilmineau, P. Queutey, M. Visonneau & J. Wackers (ECN-LHEEA,CNRS, France) A. del Toro Llorens (Spanish Institution

More information

Propeller Loads of Large Commercial Vessels at Crash Stop

Propeller Loads of Large Commercial Vessels at Crash Stop Second International Symposium on Marine Propulsors smp 11, Hamburg, Germany, June 2011 Propeller Loads of Large Commercial Vessels at Crash Stop J.W. Hur, H. Lee, B.J. Chang 1 1 Hyundai Heavy Industries,

More information

CFD Based Hull Hydrodynamic Forces for Simulation of Ship Manoeuvres

CFD Based Hull Hydrodynamic Forces for Simulation of Ship Manoeuvres International Journal on Marine avigation and Safety of Sea Transportation Volume 3 umber 1 March 9 Based Hull Hydrodynamic Forces for Simulation of Ship Manoeuvres T. Tabaczek, T. Gornicz & J. Kulczyk

More information

The Calculations of Propeller Induced Velocity by RANS and Momentum Theory

The Calculations of Propeller Induced Velocity by RANS and Momentum Theory J. Marine Sci. Appl. (2012) 11: 164-168 DOI: 10.1007/s11804-012-1118-1 The Calculations of Propeller Induced Velocity by RANS and Momentum Theory Qiuxin Gao *, Wei Jin and Dracos Vassalos Department of

More information

Manoeuvring Committee Report & Recommendations

Manoeuvring Committee Report & Recommendations Manoeuvring Committee Report & Recommendations Andrés Cura Hochbaum 26 th International Towing Tank Conference, 28 th August 4 th September 2011, Rio de Janeiro Members & Meetings Tasks Report 1. Introduction

More information

Confined water effects on the viscous flow around a tanker with propeller and rudder

Confined water effects on the viscous flow around a tanker with propeller and rudder International Shipbuilding Progress 60 (2013) 309 343 309 DOI 10.3233/ISP-130101 IOS Press Confined water effects on the viscous flow around a tanker with propeller and rudder L. Zou and L. Larsson Chalmers

More information

developed at "Dunarea de Jos" University of Galati Presented by: BENZOHRA Abdelmalek

developed at Dunarea de Jos University of Galati Presented by: BENZOHRA Abdelmalek Master Thesis presented in partial fulfillment of the requirements for the double degree: Advanced Master in Naval Architecture conferred by University of Liege "Master of Sciences in Applied Mechanics,

More information

NUMERICAL SIMULATION OF THE SHIP BOTTOM INTERACTION OF DTC CONTAINER CARRIER FOR DIFFERENT KEEL CLEARANCE IN PURE SWAY MOTION

NUMERICAL SIMULATION OF THE SHIP BOTTOM INTERACTION OF DTC CONTAINER CARRIER FOR DIFFERENT KEEL CLEARANCE IN PURE SWAY MOTION 4th MASHCON, Hamburg - Uliczka et al. (eds) - 2016 Bundesanstalt für Wasserbau ISBN 978-3-939230-38-0 (Online) DOI: 10.18451/978-3-939230-38-0_9 NUMERICAL SIMULATION OF THE SHIP BOTTOM INTERACTION OF DTC

More information

Trajectory Tracking of a Near-Surface Torpedo using Numerical Methods

Trajectory Tracking of a Near-Surface Torpedo using Numerical Methods ISSN (Print) : 2347-671 An ISO 3297: 27 Certified Organization Vol.4, Special Issue 12, September 215 Trajectory Tracking of a Near-Surface Torpedo using Numerical Methods Anties K. Martin, Anubhav C.A.,

More information

Circular motion tests and uncertainty analysis for ship maneuverability

Circular motion tests and uncertainty analysis for ship maneuverability DOI 1.17/s773-9-65-2 ORIGINAL ARTICLE Circular motion tests and uncertainty analysis for ship maneuverability Michio Ueno Æ Yasuo Yoshimura Æ Yoshiaki Tsukada Æ Hideki Miyazaki Received: 13 October 2 /

More information

Numerical Study of the Roll Decay of Intact and Damaged Ships by Q. Gao and D. Vassalos

Numerical Study of the Roll Decay of Intact and Damaged Ships by Q. Gao and D. Vassalos Session 7 Stability of Damaged Ships Numerical Simulation of Progressive Flooding and Capsize Numerical Study of the Roll Decay of Intact and Damaged Ships by Q. Gao and D. Vassalos Qiuxin Gao and Dracos

More information

Calculation of hydrodynamic manoeuvring coefficients using viscous-flow calculations

Calculation of hydrodynamic manoeuvring coefficients using viscous-flow calculations Calculation of hydrodynamic manoeuvring coefficients using viscous-flow calculations Serge Toxopeus Maritime Research Institute Netherlands (MARIN), Wageningen, The Netherlands Delft University of Technology,

More information

ITTC Recommended Procedures and Guidelines

ITTC Recommended Procedures and Guidelines 7.5-06 - Page of 2 Table of Contents Captive Model Test Procedure... 2. PURPOSE OF PROCEDURE... 2 2. DESCRIPTION OF PROCEDURE... 3 2. Preparation... 3 2.. Selection of model dimensions.... 3 2..2 Model

More information

ITTC Recommended Procedures and Guidelines Testing and Extrapolation Methods Propulsion, Performance Propulsion Test

ITTC Recommended Procedures and Guidelines Testing and Extrapolation Methods Propulsion, Performance Propulsion Test 7.5- Page 1 of 13 Table of Contents... 2 1. PURPOSE OF PROCEDURE... 2 2. PARAMETERS... 2 2.1 Data Reduction Equations... 2 2.2 Definition of Variables... 3 3. DESCRIPTION OF PROCEDURE... 3 3.1 Model and

More information

Calculation of the Flow around the KVLCC2M Tanker

Calculation of the Flow around the KVLCC2M Tanker Calculation of the Flow around the KVLCC2M Tanker L. Eça 1, M. Hoekstra 2 and S.L. Toxopeus 2 1 Instituto Superior Técnico, Portugal 2 Maritime Research Institute, Netherlands SUMMARY The flow around the

More information

ITTC Recommended Procedures Testing and Extrapolation Methods Resistance Resistance Test

ITTC Recommended Procedures Testing and Extrapolation Methods Resistance Resistance Test -0- Page 1 of 11 CONTENTS 1. PURPOSE OF PROCEDURE. PARAMETERS.1 Data Reduction Equations. Definition of ariables 3. DESCRIPTION OF PROCEDURE 3.1 Model and Installation 3.1.1 Model 3.1. Test condition 3.1.3

More information

ZIG-ZAG MANEUVER SIMULATION BY CFD FOR A TANKER LIKE VESSEL

ZIG-ZAG MANEUVER SIMULATION BY CFD FOR A TANKER LIKE VESSEL V International Conference on Computational Methods in Marine Engineering MARINE 2013 B. Brinkmann and P. Wriggers (Eds) ZIG-ZAG MANEUVER SIMULATION BY CFD FOR A TANKER LIKE VESSEL G. DUBBIOSO, D. DURANTE,

More information

A Study on Effects of Blade Pitch on the Hydrodynamic Performances of a Propeller by Using CFD

A Study on Effects of Blade Pitch on the Hydrodynamic Performances of a Propeller by Using CFD Journal of Shipping and Ocean Engineering 8 (2018) 36-42 doi 10.17265/2159-5879/2018.01.005 D DAVID PUBLISHING A Study on Effects of Blade Pitch on the Hydrodynamic Performances of a Propeller by Using

More information

A Study of Resistance of High-Speed Catamarans and Scale Effects on Form Factor Sarawuth Srinakaew Dominic J. Taunton Dominic A.

A Study of Resistance of High-Speed Catamarans and Scale Effects on Form Factor Sarawuth Srinakaew Dominic J. Taunton Dominic A. A Study of Resistance of High-Speed Catamarans and Scale Effects on Form Factor Sarawuth Srinakaew, Dominic J. Taunton, Dominic A. Hudson Fluid Structure Interactions Group, University of Southampton ss17e111@soton.ac.uk

More information

ITTC Recommended Procedures and Guidelines

ITTC Recommended Procedures and Guidelines Page 1 of 9 CONTENTS Model Test Experiments... 2 1. PURPOSE OF PROCEDURE... 2 2. PARAMETERS... 2 2.1 Model Parameters... 3 2.2 Environmental Parameters... 3 2.3 Operation of Thrusters... 3 2.3.1 Thruster-Current

More information

Deliverable D.6.1. Application of CFD tools to the development of a novel propulsion concept

Deliverable D.6.1. Application of CFD tools to the development of a novel propulsion concept TRIple Energy Saving by Use of CRP, CLT and PODded Propulsion Grant Agreement Number: 265809 Call identifier: FP7-SST-2010-RTD-1 Theme SST.2010.1.1-2.: Energy efficiency of ships WP 1 Deliverable D.6.1

More information

Comparison of Thruster Axis Tilting versus Nozzle Tilting on the Propeller-Hull Interactions for a Drillship at DP-Conditions

Comparison of Thruster Axis Tilting versus Nozzle Tilting on the Propeller-Hull Interactions for a Drillship at DP-Conditions DYNAMIC POSITIONING CONFERENCE October 12-13, 2010 Thrusters Comparison of Thruster Axis Tilting versus Nozzle Tilting on the Propeller-Hull Interactions for a Drillship at DP-Conditions Michael Palm,

More information

Full scale thruster performance and load determination based on numerical simulations

Full scale thruster performance and load determination based on numerical simulations Third International Symposium on Marine Propulsors smp 13, Launceston, Tasmania, Australia, May 213 Full thruster performance and load determination based on numerical simulations Norbert Bulten 1, Rik

More information

On the advanced extrapolation method for a new type of podded propulsor via CFD simulations and model measurements

On the advanced extrapolation method for a new type of podded propulsor via CFD simulations and model measurements Fifth International Symposium on Marine Propulsors smp 17, Espoo, Finland, June 2017 On the advanced extrapolation method for a new type of podded propulsor via CFD simulations and model measurements Tomi

More information

FACTORS OF SAFETY FOR RICHARDSON EXTRAPOLATION

FACTORS OF SAFETY FOR RICHARDSON EXTRAPOLATION FACTORS OF SAFETY FOR RICHARDSON ETRAPOLATION by Tao ing 1 and Fred Stern Sponsored by The Office of Naval Research Grant N00014-01-1-0073 IIHR Technical Report No. 476 IIHR Hydroscience & Engineering

More information

ITTC Recommended Procedures and Guidelines

ITTC Recommended Procedures and Guidelines Page of Table of Contents Waterjet Propulsion Test and Extrapolation... PURPOSE OF PROCEDURE.... PARAMETERS.... Nomenclature... 3. DESCRIPTION OF PROCEDURE... 3 3. Model and installation... 3 3.. Resistance

More information

CFD PREDICTION AND VALIDATION OF SHIP-BANK INTERACTION IN A CANAL

CFD PREDICTION AND VALIDATION OF SHIP-BANK INTERACTION IN A CANAL CFD PREDICTION AND VALIDATION OF SHIP-BANK INTERACTION IN A CANAL L Zou, Chalmers University of Technology, Sweden L Larsson, Chalmers University of Technology, Sweden G Delefortrie, Flanders Hydraulics

More information

ADVANCES IN FULL-SCALE WAKE-FIELD PREDICTIONS AND THE IMPLICATIONS FOR THE PROPELLER DESIGN

ADVANCES IN FULL-SCALE WAKE-FIELD PREDICTIONS AND THE IMPLICATIONS FOR THE PROPELLER DESIGN ADVANCES IN FULL-SCALE WAKE-FIELD PREDICTIONS AND THE IMPLICATIONS FOR THE PROPELLER DESIGN Gert-Jan Zondervan*, Bram Starke Maritime Research Institute Netherlands PO Box 28, 67 AA Wageningen, The Netherlands

More information

CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH

CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH 82 CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH The coefficient of lift, drag and power for wind turbine rotor is optimized using an iterative approach. The coefficient

More information

Numerical Analysis of Unsteady Open Water Characteristics of Surface Piercing Propeller

Numerical Analysis of Unsteady Open Water Characteristics of Surface Piercing Propeller Third International Symposium on Marine Propulsors smp 13, Launceston, Tasmania, Australia, May 2013 Numerical Analysis of Unsteady Open Water Characteristics of Surface Piercing Propeller Kohei Himei

More information

Estimating Maneuvering and Seakeeping Characteristics with Neural Networks

Estimating Maneuvering and Seakeeping Characteristics with Neural Networks 070131-071 1 Estimating Maneuvering and Seakeeping Characteristics with Neural Networks Paulo Triunfante Martins 1, Victor Lobo 2, Member, IEEE Maneuvering and seakeeping are two very important naval architecture

More information

Offshore Hydromechanics Module 1

Offshore Hydromechanics Module 1 Offshore Hydromechanics Module 1 Dr. ir. Pepijn de Jong 6. Real Flows part 2 Introduction Topics of Module 1 Problems of interest Chapter 1 Hydrostatics Chapter 2 Floating stability Chapter 2 Constant

More information

Influence of Flow Transition on Open and Ducted Propeller Characteristics

Influence of Flow Transition on Open and Ducted Propeller Characteristics Fourth International Symposium on Marine Propulsors smp 15, Austin, Texas, USA, June 2015 Influence of Flow Transition on Open and Ducted Propeller Characteristics Anirban Bhattacharyya 1, Jan Clemens

More information

Propeller Analysis Using RANS/BEM Coupling Accounting for Blade Blockage

Propeller Analysis Using RANS/BEM Coupling Accounting for Blade Blockage DRDC-RDDC-2015-N005 Fourth International Symposium on Marine Propulsors smp 15, Austin, Texas, USA, June 2015 Propeller Analysis Using RANS/BEM Coupling Accounting for Blade Blockage David Hally 1 1 Defence

More information

RANS Computational Fluid Dynamics Predictions of Pitch and Heave Ship Motions in Head Seas

RANS Computational Fluid Dynamics Predictions of Pitch and Heave Ship Motions in Head Seas Journal of Ship Research, Vol. 49, No. 2, June 2005, pp. 80 97 RANS Computational Fluid Dynamics Predictions of Pitch and Heave Ship Motions in Head Seas Gabriel David Weymouth, Robert Vance Wilson, and

More information

ITTC Propeller Benchmark

ITTC Propeller Benchmark ITTC Propeller Benchmark PPTC PPTC, Propeller P1727 - VP1304 Report 4488 Potsdam, March April 2016 2016 Schiffbau-Versuchsanstalt Potsdam GmbH, Marquardter Chaussee 100, 14469 Potsdam Tel. +49 331 56712-0,

More information

Experimental and Numerical Analysis of the Roll Decay Motion for a Patrol Boat

Experimental and Numerical Analysis of the Roll Decay Motion for a Patrol Boat Experimental and Numerical Analysis of the Roll Decay Motion for a Patrol Boat R. Broglia, B. Bouscasse, A. Di Mascio and C. Lugni INSEAN Italian Ship Model Basin Rome, Italy P. Atsavapranee NSWCCD Naval

More information

Numerical Predictions of Hydrodynamic Forces and Squat of Ships in Confined Waters

Numerical Predictions of Hydrodynamic Forces and Squat of Ships in Confined Waters Numerical Predictions of Hydrodynamic Forces and Squat of Ships in Confined Waters *Y. Liu¹, L. Zou 1,2, Z. J. Zou 1,2, T.C. Lu 1 and J.X. Liu 3,4 1 School of Naval Architecture, Ocean & Civil Engineering,

More information

ITTC Propeller Benchmark

ITTC Propeller Benchmark ITTC Propeller Benchmark Tip Rake PPTC, Propeller P1727 - P1727 Report 4487 Potsdam, March April 2016 2016 Schiffbau-Versuchsanstalt Potsdam GmbH, Marquardter Chaussee 100, 14469 Potsdam Tel. +49 331 56712-0,

More information

Reliability assessment of ship powering performance extrapolations using Monte Carlo methods

Reliability assessment of ship powering performance extrapolations using Monte Carlo methods Third International Symposium on Marine Propulsors smp 13, Launceston, Tasmania, Australia, May 2013 Reliability assessment of ship powering performance extrapolations using Monte Carlo methods Iwan M.

More information

Research on Prediction of Ship Manoeuvrability

Research on Prediction of Ship Manoeuvrability Journal of Shipping and Ocean Engineering 8 (08 30-35 doi 0.765/59-5879/08.0.004 D DAVID PUBLISHING Research on Prediction of Ship Manoeuvrability CUI Jian, WU Zixin and CHEN Weimin Shanghai Ship and Shipping

More information

Numerical Investigation of the Impact of SES-Waterjet Interactions and Flow Non-uniformity on Pump Performance

Numerical Investigation of the Impact of SES-Waterjet Interactions and Flow Non-uniformity on Pump Performance 11 th International Conference on Fast Sea Transportation FAST 2011, Honolulu, Hawaii, USA, September 2011 Numerical Investigation of the Impact of SES-Waterjet Interactions and Flow Non-uniformity on

More information

ITTC Recommended Procedures

ITTC Recommended Procedures 7.5-0 -03-0. Page of 6 CONTENTS PURPOSE OF PROCEDURE EXAMPLE FOR OPEN WATER TEST. Test Design. Measurement System and Procedure.3 Uncertainty Analysis.3. ias Limit.3.. Propeller Geometry.3.. Speed of advance

More information

Motions and Resistance of a Ship in Regular Following Waves

Motions and Resistance of a Ship in Regular Following Waves Reprinted: 01-11-2000 Revised: 03-10-2007 Website: www.shipmotions.nl Report 440, September 1976, Delft University of Technology, Ship Hydromechanics Laboratory, Mekelweg 2, 2628 CD Delft, The Netherlands.

More information

SEAKEEPING AND MANEUVERING Prof. Dr. S. Beji 2

SEAKEEPING AND MANEUVERING Prof. Dr. S. Beji 2 SEAKEEPING AND MANEUVERING Prof. Dr. S. Beji 2 Ship Motions Ship motions in a seaway are very complicated but can be broken down into 6-degrees of freedom motions relative to 3 mutually perpendicular axes

More information

MASTER THESIS PRESENTATION

MASTER THESIS PRESENTATION MASTER THESIS PRESENTATION TURBULENT FLOW SEPARATION AROUND A ROV BODY Presented by: KOFFI Danoh Kouassi Ange Gatien Outline 1- Introduction 2- Numerical modeling of the ROV 3- Choice of the turbulence

More information

Computation for the Backward Facing Step Test Case with an Open Source Code

Computation for the Backward Facing Step Test Case with an Open Source Code Computation for the Backward Facing Step Test Case with an Open Source Code G.B. Deng Equipe de Modélisation Numérique Laboratoire de Mécanique des Fluides Ecole Centrale de Nantes 1 Rue de la Noë, 44321

More information

CFD Analysis of Scale Effects on Conventional and Tip-Modified Propellers

CFD Analysis of Scale Effects on Conventional and Tip-Modified Propellers Fifth International Symposium on Marine Propulsors smp 17, Espoo, Finland, June 2017 CFD Analysis of Scale Effects on Conventional and Tip-Modified Propellers Keun Woo Shin 1, Poul Andersen 2 1 Propeller

More information

Numerical Simulation of a Complete Francis Turbine including unsteady rotor/stator interactions

Numerical Simulation of a Complete Francis Turbine including unsteady rotor/stator interactions Numerical Simulation of a Complete Francis Turbine including unsteady rotor/stator interactions Ruprecht, A., Heitele, M., Helmrich, T. Institute for Fluid Mechanics and Hydraulic Machinery University

More information

Uncertainty assessment for free-running model cases at the IIHR wave basin

Uncertainty assessment for free-running model cases at the IIHR wave basin University of Iowa Iowa Research Online Theses and Dissertations Summer 016 Uncertainty assessment for free-running model cases at the IIHR wave basin Michael John Bottiglieri University of Iowa Copyright

More information

ITTC Recommended Procedures Testing and Extrapolation Methods Manoeuvrability Captive Model Test Procedure

ITTC Recommended Procedures Testing and Extrapolation Methods Manoeuvrability Captive Model Test Procedure Page of 7 CONTENTS. PURPOSE OF PROCEDURE. DESCRIPTION OF PROCEDURE. Preparation.. Selection of model dimensions... Model inspection. The model should be inspected, prior to launching and testing, for:..3

More information

Alexander Phillips 1, Stephen R. Turnock 1, Maaten Furlong 2 School of Engineering Sciences, University of Southampton

Alexander Phillips 1, Stephen R. Turnock 1, Maaten Furlong 2 School of Engineering Sciences, University of Southampton Evaluation of manoeuvring coefficients of a self-propelled ship using a Blade Element Momentum propeller model coupled to a Reynolds averaged Navier Stokes flow solver Abstract 1 Alexander Phillips 1,

More information

Performance Prediction of the Francis-99 Hydroturbine with Comparison to Experiment. Chad Custer, PhD Yuvraj Dewan Artem Ivashchenko

Performance Prediction of the Francis-99 Hydroturbine with Comparison to Experiment. Chad Custer, PhD Yuvraj Dewan Artem Ivashchenko Performance Prediction of the Francis-99 Hydroturbine with Comparison to Experiment Chad Custer, PhD Yuvraj Dewan Artem Ivashchenko Unrestricted Siemens AG 2017 Realize innovation. Agenda Introduction

More information

Measurement of speed loss due to waves

Measurement of speed loss due to waves Third International Symposium on Marine Propulsors smp 13, Launceston, Tasmania, Australia, May 213 Measurement of speed loss due to waves Sverre Steen 1 and Zhenju Chuang 1 1 Department of Marine Technology,

More information

Local Velocity Field Measurements around the KCS Model (SRI M.S.No.631) in the SRI 400m Towing Tank

Local Velocity Field Measurements around the KCS Model (SRI M.S.No.631) in the SRI 400m Towing Tank Local Velocity Field Measurements around the KCS Model (SRI M.S.No.631) in the SRI 400m Towing Tank 1. Introduction The present report describes the local velocity field measurements around a model of

More information

This thesis is submitted for evaluation at Ålesund University College.

This thesis is submitted for evaluation at Ålesund University College. MASTER THESIS TITLE: Scale effect on the wake field of a Single Screw Ship CANDIDATE NAME: Yue Ding Candidate No. 2110 DATE: COURSE CODE: COURSE TITLE: RESTRICTION: 29/05/2015 IP501909 MSc thesis, discipline

More information

AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract

AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo

More information

MARINE PROPULSION SYSTEM DYNAMICS DURING SHIP MANOEUVRES

MARINE PROPULSION SYSTEM DYNAMICS DURING SHIP MANOEUVRES MARINE PROPULSION SYSTEM DYNAMICS DURING SHIP MANOEUVRES Michele Viviani, Department of Naval Architecture and Marine Engineering (DINAV), Genoa University, Italy Marco Altosole, Department of Naval Architecture

More information

Strathprints Institutional Repository

Strathprints Institutional Repository Strathprints Institutional Repository Tezdogan, Tahsin and Incecik, Atilla and Turan, Osman (2016) A numerical investigation of the squat and resistance of ships advancing through a canal using CFD. Journal

More information

Wake fraction and thrust deduction during ship astern manoeuvres

Wake fraction and thrust deduction during ship astern manoeuvres Wake fraction and thrust deduction during ship astern manoeuvres J. Artyszuk Maritime University of Szczecin, Poland Abstract A relatively small amount of data concerning the behaviour of propulsion coefficients,

More information

Is My CFD Mesh Adequate? A Quantitative Answer

Is My CFD Mesh Adequate? A Quantitative Answer Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011

More information

The Numerical Investigation on Hydrodynamic Performance of Twisted Rudder during Self-propulsion

The Numerical Investigation on Hydrodynamic Performance of Twisted Rudder during Self-propulsion The Numerical Investigation on Hydrodynamic Performance of Twisted Rudder during Self-propulsion Cong Liu *, Jianhua Wang,Decheng Wan State Key Laboratory of Ocean Engineering, School of Naval Architecture,

More information

Department of Aerospace and Ocean Engineering Graduate Study Specialization in Ocean Engineering. Written Preliminary Examination Information

Department of Aerospace and Ocean Engineering Graduate Study Specialization in Ocean Engineering. Written Preliminary Examination Information Department of Aerospace and Ocean Engineering Graduate Study Specialization in Ocean Engineering Written Preliminary Examination Information Faculty: Professors W. Neu, O. Hughes, A. Brown, M. Allen Test

More information

Study of the hydrodynamic flow around a 70m sailing boat for powering, wave pattern and propeller efficiency prediction

Study of the hydrodynamic flow around a 70m sailing boat for powering, wave pattern and propeller efficiency prediction Study of the hydrodynamic flow around a 70m sailing boat for powering, wave pattern and propeller efficiency prediction Romain Baudson Supervisor: Prof. Dario Boote (UNIGE) External Reviewer: Prof. Pierre

More information

OMAE OMAE

OMAE OMAE Proceedings Proceedings of the ASME of ASME 2010 29th International Conference on on Ocean, Offshore and Arctic Engineering OMAE2010 June 6-11, 2010, Shanghai, China June 6-11, 2010, Shanghai, China OMAE2010-20373

More information

Potsdam Propeller Test Case (PPTC) Test Case Description

Potsdam Propeller Test Case (PPTC) Test Case Description Second International Symposium on Marine Propulsors smp 11, Hamburg, Germany, June 2011 Workshop: Propeller performance Potsdam Propeller Test Case (PPTC) Test Case Description Ulf Barkmann 1, Hans-Jürgen

More information

Transactions on the Built Environment vol 24, 1997 WIT Press, ISSN

Transactions on the Built Environment vol 24, 1997 WIT Press,  ISSN Comparison of model test with ship sea trial results for a given vessel series C.Behrendt & T.Kucharski Institute of Marine Plant Operation, Maritime University of Szczecin, 70-500 Szczecin, Poland Abstract

More information

Sailing Performance and Maneuverability of a Traditional Ryukyuan Tribute Ship

Sailing Performance and Maneuverability of a Traditional Ryukyuan Tribute Ship sia Navigation Conference 009 ailing Performance and Maneuverability of a Traditional Ryukyuan Tribute hip by Yutaka MUYM (Kanazawa Institute of Technology) Hikaru YGI (Tokai University ) Yutaka TERO (Tokai

More information

Study on Strut Effect on Turning Characteristics of LNG Carrier

Study on Strut Effect on Turning Characteristics of LNG Carrier Study on Strut Effect on Turning Characteristics of LNG Carrier Mohd Amirul Hamdi Bin Mohd Alwi, a Jaswar Koto, a, b,* and Zulkarnain, b a) Department of Aeronautics, Automotive and Ocean Engineering,Mechanical

More information

PLEASURE VESSEL VIBRATION AND NOISE FINITE ELEMENT ANALYSIS

PLEASURE VESSEL VIBRATION AND NOISE FINITE ELEMENT ANALYSIS PLEASURE VESSEL VIBRATION AND NOISE FINITE ELEMENT ANALYSIS 1 Macchiavello, Sergio *, 2 Tonelli, Angelo 1 D Appolonia S.p.A., Italy, 2 Rina Services S.p.A., Italy KEYWORDS pleasure vessel, vibration analysis,

More information

Large Eddy Simulation of Crashback in Marine Propellers

Large Eddy Simulation of Crashback in Marine Propellers 26th Symposium on Naval Hydrodynamics Rome, Italy, 17 22 September 2006 Large Eddy Simulation of Crashback in Marine Propellers Martin Vyšohlíd, Krishnan Mahesh (Aerospace Engineering & Mechanics, University

More information

Investigation of ship maneuvering with hydrodynamic effects between ship and bank

Investigation of ship maneuvering with hydrodynamic effects between ship and bank Journal of Mechanical Science and Technology Journal of Mechanical Science and Technology (8) 3~36 www.springerlink.com/content/738-494 Investigation of ship maneuvering with hydrodynamic effects between

More information

Discretization error analysis with unfavorable meshes A case study

Discretization error analysis with unfavorable meshes A case study Discretization error analysis with unfavorable meshes A case study Denis F. Hinz and Mario Turiso Kamstrup A/S, Denmark ASME 2016 V&V Symposium May 16 20, 2016, Las Vegas, Nevada Agenda Who we are Background

More information

A Mathematical Model for the Maneuvering Simulation of a. Propelled SPAR Vessel. Soroosh Mohammadafzali. A thesis submitted to the

A Mathematical Model for the Maneuvering Simulation of a. Propelled SPAR Vessel. Soroosh Mohammadafzali. A thesis submitted to the A Mathematical Model for the Maneuvering Simulation of a Propelled SPAR Vessel By Soroosh Mohammadafzali A thesis submitted to the School of Graduate Studies in partial fulfillment of the requirements

More information

Simulation of Aeroelastic System with Aerodynamic Nonlinearity

Simulation of Aeroelastic System with Aerodynamic Nonlinearity Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School

More information

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force

More information

CFD Analysis of a Propeller Flow and Cavitation

CFD Analysis of a Propeller Flow and Cavitation CFD Analysis of a Propeller Flow and Cavitation S. Ramakrishna* GVP College of Engineering (A) Visakhapatnam V. Ramakrishna Scientist D N S T L Visakhapatnam A. Ramakrishna Professor A U C E (A) Visakhapatnam

More information

ENGR 4011 Resistance & Propulsion of Ships Assignment 4: 2017

ENGR 4011 Resistance & Propulsion of Ships Assignment 4: 2017 Question 1a. Values of forward speed, propeller thrust and torque measured during a propeller open water performance test are presented in the table below. The model propeller was 0.21 meters in diameter

More information

Journal of Ship Research

Journal of Ship Research Journal of Ship Research, Vol. 53, No. 4, December 009, pp. 179 198 Journal of Ship Research Model- and Full-Scale URANS Simulations of Athena Resistance, Powering, Seakeeping, and 5415 Maneuvering Shanti

More information

Final Exam TTK4190 Guidance and Control

Final Exam TTK4190 Guidance and Control Trondheim Department of engineering Cybernetics Contact person: Professor Thor I. Fossen Phone: 73 59 43 61 Cell: 91 89 73 61 Email: tif@itk.ntnu.no Final Exam TTK4190 Guidance and Control Friday May 15,

More information

Studying Fluid Mechanics for Better Understanding Simulation Technologies. A Special Approach for Maritime Universities

Studying Fluid Mechanics for Better Understanding Simulation Technologies. A Special Approach for Maritime Universities Studying Fluid Mechanics for Better Understanding Simulation Technologies. A Special Approach for Maritime Universities DUMITRU DINU Mechanical Engineer Department Constanta Maritime University 104 Mircea

More information

The effect of top tension on VIV model analysis of a vertical flexible riser

The effect of top tension on VIV model analysis of a vertical flexible riser The Second Conference of Global Chinese Scholars on Hydrodynamics The effect of top tension on VIV model analysis of a vertical flexible riser Muyu Duan 1,2, Bowen Fu 1, Decheng Wan 1* 1 State Key Laboratory

More information

Computation of Unsteady Flows With Moving Grids

Computation of Unsteady Flows With Moving Grids Computation of Unsteady Flows With Moving Grids Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de Unsteady Flows With Moving Boundaries, I Unsteady flows

More information

Simple Estimation of Wave Added Resistance from Experiments in Transient and Irregular Water Waves

Simple Estimation of Wave Added Resistance from Experiments in Transient and Irregular Water Waves Simple Estimation of Wave Added Resistance from Experiments in Transient and Irregular Water Waves by Tsugukiyo Hirayama*, Member Xuefeng Wang*, Member Summary Experiments in transient water waves are

More information

Performance and wake development behind two in-line and offset model wind turbines "Blind test" experiments and calculations

Performance and wake development behind two in-line and offset model wind turbines Blind test experiments and calculations Journal of Physics: Conference Series OPEN ACCESS Performance and wake development behind two in-line and offset model wind turbines "Blind test" experiments and calculations To cite this article: Lars

More information

EFFECT OF FORCED ROTATING VANELESS DIFFUSERS ON CENTRIFUGAL COMPRESSOR STAGE PERFORMANCE

EFFECT OF FORCED ROTATING VANELESS DIFFUSERS ON CENTRIFUGAL COMPRESSOR STAGE PERFORMANCE Journal of Engineering Science and Technology Vol. 6, No. 5 (2011) 558-574 School of Engineering, Taylor s University EFFECT OF FORCED ROTATING VANELESS DIFFUSERS ON CENTRIFUGAL COMPRESSOR STAGE PERFORMANCE

More information

DETERMINING OPTIMUM GEOMETRY FOR A BOW THRUSTER PROPELLER

DETERMINING OPTIMUM GEOMETRY FOR A BOW THRUSTER PROPELLER DEERMINING OPIMUM GEOMERY FOR A BOW HRUSER PROPELLER Y.H. OZDEMIR, S. BAYRAKAR,. YILMAZ, M.GUNER Yildiz echnical University, Dept. of Naval Architecture and Marine Engineering, 34349 Besiktas, Istanbul,

More information

Validation of Computational Fluid-Structure Interaction Analysis Methods to Determine Hydrodynamic Coefficients of a BOP Stack

Validation of Computational Fluid-Structure Interaction Analysis Methods to Determine Hydrodynamic Coefficients of a BOP Stack Validation of Computational Fluid-Structure Interaction Analysis Methods to Determine Hydrodynamic Coefficients of a BOP Stack The MIT Faculty has made this article openly available. Please share how this

More information

Numerical Analysis of Azimuth Propulsor Performance in Seaways: Influence of Oblique Inflow and Free Surface

Numerical Analysis of Azimuth Propulsor Performance in Seaways: Influence of Oblique Inflow and Free Surface Journal of Marine Science and Engineering Article Numerical Analysis of Azimuth Propulsor Performance in Seaways: Influence of Oblique Inflow and Free Surface Nabila Berchiche, Vladimir I. Krasilnikov

More information

On the transient modelling of impinging jets heat transfer. A practical approach

On the transient modelling of impinging jets heat transfer. A practical approach Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,

More information

Investigation of Flow Profile in Open Channels using CFD

Investigation of Flow Profile in Open Channels using CFD Investigation of Flow Profile in Open Channels using CFD B. K. Gandhi 1, H.K. Verma 2 and Boby Abraham 3 Abstract Accuracy of the efficiency measurement of a hydro-electric generating unit depends on the

More information

The Simulation of Wraparound Fins Aerodynamic Characteristics

The Simulation of Wraparound Fins Aerodynamic Characteristics The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:

More information

The... of a particle is defined as its change in position in some time interval.

The... of a particle is defined as its change in position in some time interval. Distance is the. of a path followed by a particle. Distance is a quantity. The... of a particle is defined as its change in position in some time interval. Displacement is a.. quantity. The... of a particle

More information

Draft Tube calculations using OpenFOAM-1.5dev and validation with FLINDT data

Draft Tube calculations using OpenFOAM-1.5dev and validation with FLINDT data 6th OpenFOAM Workshop, June 13-16 2011, PennState University, USA Draft Tube calculations using OpenFOAM-1.5dev and validation with FLINDT data C. Devals, Y. Zhang and F.Guibault École Polytechnique de

More information

Hydrodynamic Forces due to Orbital Stokes 5 th Order Waves on Subsea Pipelines Resting on Porous Seabed

Hydrodynamic Forces due to Orbital Stokes 5 th Order Waves on Subsea Pipelines Resting on Porous Seabed Hydrodynamic Forces due to Orbital Stokes 5 th Order Waves on Subsea Pipelines Resting on Porous Seabed Annelise Karreman Dr Jeremy Leggoe School of Mechanical and Chemical Engineering W/Prof Liang Cheng

More information

PRINCIPLES OF FLIGHT

PRINCIPLES OF FLIGHT 1 Considering a positive cambered aerofoil, the pitching moment when Cl=0 is: A infinite B positive (nose-up). C negative (nose-down). D equal to zero. 2 The angle between the aeroplane longitudinal axis

More information

MODELLING THE DYNAMICS OF SHIPS WITH DIFFERENT PROPULSION SYSTEMS FOR CONTROL PURPOSE

MODELLING THE DYNAMICS OF SHIPS WITH DIFFERENT PROPULSION SYSTEMS FOR CONTROL PURPOSE POLISH MARITIME RESEARCH 1(89) 2016 Vol. 23; pp. 31-36 10.1515/pomr-2016-0005 MODELLING THE DYNAMICS OF SHIPS WITH DIFFERENT PROPULSION SYSTEMS FOR CONTROL PURPOSE Witold Gierusz, Assoc. Prof. Gdynia Maritime

More information

Trajectory tracking & Path-following control

Trajectory tracking & Path-following control Cooperative Control of Multiple Robotic Vehicles: Theory and Practice Trajectory tracking & Path-following control EECI Graduate School on Control Supélec, Feb. 21-25, 2011 A word about T Tracking and

More information

Numerical Prediction Of Torque On Guide Vanes In A Reversible Pump-Turbine

Numerical Prediction Of Torque On Guide Vanes In A Reversible Pump-Turbine Journal of Multidisciplinary Engineering Science and Technology (JMEST) ISSN: 3159 Vol. 2 Issue 6, June - 215 Numerical Prediction Of Torque On Guide Vanes In A Reversible Pump-Turbine Turbine and pump

More information

Finite Element Modelling with Plastic Hinges

Finite Element Modelling with Plastic Hinges 01/02/2016 Marco Donà Finite Element Modelling with Plastic Hinges 1 Plastic hinge approach A plastic hinge represents a concentrated post-yield behaviour in one or more degrees of freedom. Hinges only

More information