R E D U C I N G A D D E D RESISTANCE USING A N ANTI-PITCH FOIL By John Spencer Avis B.A.Sc. (Systems Design Engineering), University of Waterloo A THESIS THE SUBMITTED IN P A R T I A L REQUIREMENTS MASTER OF FOR FULFILLMENT THE APPLIED DEGREE OF SCIENCE in THE FACULTY OF GRADUATE MECHANICAL STUDIES ENGINEERING We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA April 1989 © John Spencer Avis, 1989 OF In presenting this thesis in partial fulfilment of the degree at the and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may or by his or her representatives. be permission. Mechanical Engineering The University of British Columbia Vancouver, Canada Date DE-6 (2/88) April 28,1989 granted by It is understood publication of this thesis for financial gain shall not be Department of advanced University of British Columbia, I agree that the Library shall make it freely available for reference department requirements for an the head of my that or copying allowed without my written Abstract In this study, the use of a passive bow anti-pitch foil for the purpose of reducing the vertical motions and added resistance in waves of a small marine vehicle is investigated. The difficulties, limitations and assumptions inherent in studying ships in waves are discussed. A mathematical model is developed incorporating linear strip theory, Joosen added resistance theory and oscillating hydrofoil theory in an attempt to predict the effect of a hypothetical anti-pitch foil on the heave and pitch motions and added resistance in head seas of the 12-metre yacht Canada II. The results of regular head wave towing tank experiments on a 1:8 scale model of Canada II are presented for both the bare hull and for the hull fitted with a model antipitch foil to demonstrate the experimental effect of the foil on vertical motions and added resistance. The results of the theoretical and experimental investigations show that the anti-pitch foil is in fact very effective in reducing the motions and added resistance of the Canada II hull. Experimental results show heave amplitude is reduced by as much as 15%, pitch amplitude by 22% and added resistance by as much as 40%. The heave and pitch motion amplitude responses predicted by the linear strip theory for the hull with no foil, do not agree well with experiment, especially in the region around resonance. It is concluded that the linear strip theory used is not completely adequate in predicting the hydrodynamic coefficients and consequently the heave and pitch motions of Canada II due to non-linearities in her hull form. These non-linearities, in particular the flared topsides and large fore and aft overhangs, violate the linear assumptions of strip theory and result in increased hydrodynamic damping and speed dependent restoring ii coefficients that are not accounted for by theory. The Joosen added resistance theory predictions for the bare hull do agree well with experiment despite the hull non-linearities. The change or percentage reduction in heave and pitch amplitude response due to the addition of the anti-pitch foil is predicted well by the oscillating hydrofoil theory, in spite of the inaccuracies of the strip theory in predicting the actual motions. Also, the Joosen theory predictions of the percentage reduction in added resistance due to the addition of the foil agree very well with experiment. It is concluded that the mathematical model provides a very good conservative estimate of the amount of reduction in heave and pitch motion and added resistance that can be expected from the addition of a passive bow anti-pitch foil. The model could be very useful as a preliminary design tool for naval architects, where no such tool existed before. The findings of this study indicate that the concept of using a passive bow anti-pitch foil to reduce the vertical motions and added resistance in waves is well worth pursuing. in Table of Contents Abstract ii List of Tables viii List of Figures ix Nomenclature xi Acknowledgement xiv 1 Introduction 1 1.1 Motivation 1 1.2 Needs Analysis 3 1.3 Purpose of Research 7 1.4 Scope of Investigation 7 1.4.1 Experimental and Theoretical Constraints 2 Historical Background 8 9 2.1 Phenomenon! of Added Resistance 9 2.1.1 Causes of Added Resistance 9 2.1.2 Mechanisms of Added Resistance in Head Seas 2.2 Theoretical Techniques for the Predicton of Added Resistance 10 11 2.2.1 Havelock Theory 11 2.2.2 Joosen Theory 13 iv 2.2.3 Gerritsma and Beukelman Theory 2.3 Ship Motion Theory 15 17 2.3.1 Coupled Heave-Pitch Equations of Motion 17 2.3.2 Determination of Hydrodynamic Coefficients 18 2.3.3 Assumptions 19 2.3.4 20 Added Mass and Damping Coefficients 2.3.5 Exciting Force and Moment 2.4 Experimental Work on Ship Motions and Added Resistance 21 21 2.4.1 Experimental Techniques 21 2.4.2 Experimental Findings for Ships and Yachts 23 2.5 Reducing Motions and Added Resistance in Waves 26 2.5.1 Methods of Motion Reduction 27 2.5.2 Experimental Findings on Anti-Pitch Foils 28 3 Mathematical Model Formulation 31 3.1 Principal Parameters and Limitations 31 3.2 Frank Close-Fit Motion Program 33 3.2.1 2-D Added Mass and Damping 33 3.2.2 2-D Exciting Force and Moment 34 3.3 Modelling Effect of Anti-Pitch Foil 35 3.3.1 Foil Parameters 35 3.3.2 Anti-Pitch Foil Theory 36 3.3.3 Modification of Coefficients in Equation of Motion 37 3.4 Added Resistance Calculation 38 3.5 Results of Theoretical Calculations 40 3.5.1 Heave Motion 40 v 3.5.2 Pitch Motion 42 3.5.3 Added Resistance 44 4 Experimental Program 47 4.1 Experimental Set-up 47 4.1.1 Facilities 47 4.1.2 Model Parameters 50 4.1.3 Model Foil Design 54 4.2 Test Program 55 4.2.1 Calm Water Tests 55 4.2.2 Regular Wave Tests 58 4.3 Observations and Data Reduction 59 4.3.1 Test Observations 59 4.3.2 Data Reduction 60 4.4 Experimental Results 4.4.1 64 Calm Water Tests 64 4.4.2 Heave Motion 66 4.4.3 Pitch Motion 68 4.4.4 68 Added Resistance 5 Comparison of Theoretical and Experimental Results 6 72 5.1 Ship Motion Comparison For Bare Hull 72 5.2 Comparison of Effect of Anti-Pitch Foil on Motions 77 5.3 Added Resistance Comparison 83 Conclusions and Recommendations 88 6.1 Conclusions 88 vi 6.2 Recommendations Bibliography Appendices A Sample Designer Questionnaire B Model Foil Drawings C Experimental Data List of Tables 1.1 Summary of Designer Questionnaire 4 3.2 Canada II Principal Particulars 32 3.3 Full Scale Anti-Pitch Foil Parameters 35 4.4 Canada II 1:8 Scale Model Particulars 54 4.5 Model Anti-Pitch Foil Parameters 55 4.6 Calm Water Pitch and Heave Natural Frequencies 66 vm List of Figures 1.1 Canada II in Waves 2 2.2 Co-ordinate System 12 3.3 Theoretical Heave Amplitude Response 3.4 Theoretical Pitch Amplitude Response 41 43 3.5 Theoretical Added Resistance 45 4.6 Wave Paddle 48 4.7 Heave Post Arrangement 49 4.8 Wave Probe Set-up 51 4.9 Data Acquisition Equipment on Carriage 52 4.10 Canada II Model 53 4.11 Model Anti-Pitch Foil Fitted For Testing 56 4.12 Underwater Sideview of Model Anti-Pitch Foil 57 4.13 Bow Immersion of Canada II Bare Hull 61 4.14 Bow Immersion of Canada II Hull With Anti-Pitch Foil 62 4.15 Calm Water Model Resistance For Canada II 65 4.16 Experimental Heave Amplitude Response 67 4.17 Experimental Pitch Amplitude Response 69 4.18 Experimental Added Resistance 70 5.19 Heave Amplitude Response For Bare HulhTheory vs. Experiment 73 5.20 Pitch Amplitude Response For Bare Hull: Theory vs. Experiment 75 ix 5.21 Heave Amplitude Response With and Without Foil: Theory vs. Experiment 79 5.22 Heave Amplitude Difference: Theory vs. Experiment 80 5.23 Pitch Amplitude Response With and Without Foil: Theory vs. Experiment 81 5.24 Pitch Amplitude Difference: Theory vs. Experiment 82 5.25 Added Resistance For Bare Hull: Theory vs. Experiment 84 5.26 Added Resistance With and Without Foil: Theory vs. Experiment 85 5.27 Added Resistance Difference: Theory vs. Experiment 87 x Nomenclature A = Aspect ratio of foil Added mass coefficients B = Beam of ship Bij Damping coefficients C Wave celerity — g/u) Cij Restoring coefficients w C(K) = Theodorsen circulation function F = Heave exciting force F = Pitch exciting moment J{») = Bessel function K = Lateral coefficient of accession to inertia 3 5 Moment of inertia for pitch h L M jk = Length of ship = Mass coefficients of ship Planform area of foil = be S V = Forward velocity of ship v = Vertical velocity amplitude = Two-dimensional sectional added mass coefficient b = Foil span m = Sectional beam bij z= Two dimensional sectional damping coefficient za c d(0 - Chord of foil Sectional draft of ship /(O = Sectional exciting force 9 — Gravitational acceleration h = Depth of foil below waterline k = Wave number I = Distance of foil from LCG of ship = Sectional exciting moment Pv = 2-D pressure in phase with velocity Pa = 2-D pressure in phase with acceleration Z = Heave displacement (positive up) Heave amplitude z a Zp-p ct, ee : Experimental peak-to-peak heave = Steady angle of attack of foil (usually zero) = Heave phase angle = Pitch phase angle Ca Wave amplitude Cp-p Experimental wave height, peak-to-peak - = e Pitch angle (positive bow down) Pitch amplitude e a Op-p Angular and translatory displacements = Experimental peak-to-peak pitch XII K = Reduced frequency = A = Wavelength fi = Non-dimensional encounter frequency v — e ( ui c/2V e TTC/X = Position of longitudinal centre of gravity c p = Mass density of water °~aw — Non-dimensional added resistance coefficient 2TTI/X </> = UJ = Wave frequency ui = Frequency of encounter V = Volume of displacement e xiii Acknowledgement I would like to express my thanks to Bruce Kirby for the use of the Canada II model and hull lines drawings; to Gerry Stensgaard, George Roddan and Gary Novelsky of the B.C. Research Ocean Engineering Centre for the use of their facilities and invaluable assistance in conducting the experiments and to Dr. Bruce Dunwoody for his supervision and direction on this project. A special thanks goes to yacht designer Steve Killing who from our first meeting has never ceased to provide me with his enthusiastic support and through it all has been my major source of inspiration. Finally, I would like to acknowledge the financial assistance of the Natural Sciences and Engineering Research Council, Robert Allan Ltd. and the Canadian Institute of Marine Engineers. xiv Chapter 1 Introduction 1.1 Motivation Frustration with the performance in waves of Canada II during the 1987 America's Cup Elimination Trials in Fremantle Western Australia was the primary motivation behind this research. As a member of the sailing crew during the races, the author witnessed first hand the importance of seakeeping on the performance of a sailing yacht in a seaway. While Canada II was a solid performer in calm water conditions, her performance relative to the other yachts dropped drastically during rough wave conditions. The Canada II program was not alone in its lack of understanding of the problem. In fact, though millions of dollars were spent on hull research and development by the 13 challengers and 5 defenders vying for the Cup, very little of this was devoted to the problem of predicting added resistance in waves and subsequently performance in waves. The majority of research dollars were spent optimizing the calm water performance of the hull through computer simulations and towing tank tests. In fact, this has been the general practice in 12-metre design for many years. Yacht designers have relied on calm water performance predictions when optimizing their designs and have been forced to assume that a design that is predicted to perform better than another in calm water will also perform better in waves [37]. It should be noted that the prediction of added resistance and performance of a yacht in waves is an extremely complex and difficult problem which helps explain the lack of 1 Figure 1.1: Canada II in Waves Chapter 1. Introduction 3 useful research i n this area. It should also be noted that until 1987, the America's C u p was held i n Newport Rhode Island where w i n d conditions were predominantly light and performance i n waves was not considered to be a crucial factor. O n the other hand, i n Fremantle during the 1987 C u p trials, the winds were mostly strong (up to 30 knots) and due to the shallow water and short fetch the waves conditions were quite extreme, making performance i n waves a critical factor. A s a result, the interest i n performance prediction i n waves is relatively recent. It was the author's goal to shed more light on the problem of added resistance i n waves and investigate possible solutions that might be useful for future Canadian America's C u p challenges. W h i l e the research is directed toward the problem of added resistance of 12metre yachts, the findings should also be applicable to other types of small marine craft. 1.2 Needs Analysis Before undertaking this research, it was first necessary to establish whether a bona-fide need exists to find ways of solving the problem of reducing added resistance i n waves. In order to determine how important a factor 12-metre designers felt that added resistance was and to ascertain how they incorporated it into their 1987 12-metre designs, the author distributed a questionnaire to most of the yacht designers involed i n the 1987 America's C u p . A sample of the questionnaire appears i n A p p e n d i x A . Response to the questionnaire was good including representatives from the Stars and Stripes, Heart of A m e r i c a , Eagle, New Zealand, Canada II, and U.S.A. syndicates. A summary of the answers to the questionnaire is shown i n Table 1.1. T h e results show that a l l designers felt that reducing added resistance due to waves would lead to a significant increase i n upwind performance and was worth pursuing. In the highly irregular sea conditions encountered i n Fremantle, designers estimated that the added resistance Chapter 1. Introduction Table 1.1: Summary of Designer Questionnaire Heart Stars & Stripes USA Yes Yes To what % of t o t a l resistance? 30-40 15-30 D i d y o u have a design program f o r added r e s i s t ance? Yes poor results Is added resistance significant? Did you perforin t a n k wave t e s t s ? Some Did you use tank test results? No- u n reliable Were c o m p u t e r ship motion studies part of your program? Yes in-conclus ive Did you use any p i t c h reducing techniques i n F r e m a n 11e ? No Some New Zealand Eagle of America Canada II Yes Yes Yes 5-20 20 20 30 No used experience No No No Yes No Yes Yes No No - No- u n reliable - No No No Yes undere s t imatec No No Yes (not s p e c i f i e d ) No modified trim No Steering t e c h n i que No Chapter 1. Introduction 5 accounted for up to 50% of the total upwind resistance, making it an extremely crucial factor. In contrast, when asked about the importance placed on reducing added resistance in the actual design programs, few designers had committed any time to research and development in this area. Some had gone as far as to conduct towing tank wave tests and computer motion and added resistance computer studies, however not one had actually used these results in their final design. All the designers indicated that in the end they had to rely on their previous experience and intuition when addressing the seakeeping aspects of their designs. It is interesting to note that only Stars and Stripes (the eventual Cup winner) stated that they had used a pitch reducing technique to improve seakeeping between the start of the Trials and thefinalCup races (though they would not mention the method!). It is very clear from this polling of state of the art yacht designers that a definite need does indeed exist to gain a clearer understanding of added resistance of yachts in waves and tofindways of reducing it to improve performance in waves. While researching added resistance in waves, it became apparent to the author that added resistance was a problem for many types of marine vehicles other than racing yachts including cargo vessels, naval ships, passenger ferries,fishingvessels, drill ships, etc. In fact the bulk of the research in the area of added resistance and seakeeping has been carried out on cargo ships and naval vessels (not yachts). Some of the problems that arise as a result of inadequate seakeeping of these other forms of marine craft are: • Speed loss due to inadequate horsepower to overcome added resistance. • Voluntary speed reduction due to excessive pitching accelerations. • Increased fuel consumption to maintain speed in waves. Chapter 1. Introduction 6 • Passenger discomfort. • Inability for workers to carry out duties. • Inability for naval vessels to carry out operations (eg. landing aircraft on a carrier in rough seas). In some types of craft, the problems are quite minor such as passenger discomfort on ferries, however for vessels such as cargo ships or drill ships these problems can have major economic and scheduling implications costing millions of dollars. In the case of naval vessels, the inability to carry out operations due to rough water could, in the worst case, result in loss of life. Though, to some, the need to reduce added resistance and improve seakeeping of racing yachts in order to win races may seem somewhat trivial, one can see that throughout the marine industry there is a definite pressing need to address this problem as improved seakeeping and reduced added resistance can lead to reduced costs for cargo carriers, increased passenger comfort, increase in worker productivity, improved fishing operations, and most importantly in naval situations decrease the frequency of life threatening situations. By now the reader may have guessed that there is a clear relationship between the motions of a ship in waves and its added resistance. This relationship will be discussed quantitatively in later chapters; but for now let us state that, in general, a decrease in the motions of a ship subject to a certain sea condition will lead to a decrease in added resistance. Chapter 1. 1.3 Introduction 7 Purpose of Research The preceeding discussion leads to the primary purpose of this investigation which was to show that the use of a bow anti-pitch foil can reduce the motions and added resistance of a marine vehicle in waves. As the author's interest is primarily with racing yachts, the hull form to be investigated will be that of the 12-metre Canada II; however the arguments presented should apply to most conventional small marine vehicles. To achieve the objective of the study, a mathematical model was developed from existing linear ship motion theory, added resistance theory and oscillating hydrofoil theory in an attempt to predict the effect of an anti-pitching foil on Canada II. A model foil was designed and implemented on a 1/8 scale model of Canada II and model experiments were conducted at the B.C. Research Ocean Engineering Centre towing basin. The results of the experiments, with and without the foil, are compared to demonstrate that the foil does in fact reduce the motions and added resistance. The experimental results are also compared to the theoretical predictions in order to determine if the theory could be useful as a preliminary design tool for naval architects considering the use of an anti-pitch foil. 1.4 Scope of Investigation A literature search was conducted on the causes and mechanisms of added resistance in waves. Various methods of motion and resistance reduction were identified and studied leading to the choice of a bow anti-pitch foil as the most promising method for detailed study. Historical theoretical and experimental research into the phenomena of ship motions and added resistance was then studied in order to more precisely define the scope this research should take. Chapter 1. 1.4.1 Introduction 8 Experimental and Theoretical Constraints The study of motions and added resistance becomes extremely complex when considering sailing vessel dynamics due to the coupled interaction between the aerodynamic forces on the sails and the hydrodynamic forces on the hull and appendages [37, 32]. Thus, in order to carry out an experimental and theoretical investigation that would yield meaningful results, it became necessary to simplify the problem drastically in order to isolate the most important aspects of the problem. It is not practical to study all of the conditions the vessel might encounter on the ocean therefore this study will look only at the worst case scenerio under somewhat realistic conditions. This investigation is therefore limited to studying the motions and added resistance of the upright yacht when encountering head seas. The analysis was carried out for only one speed and one waveheight corresponding to realistic values encountered in Fremantle. Only the response of the hull and appendages are studied; the effect of the sails and rigging being neglected. The yacht is restricted to the vertical modes of motion, pitching and heaving, only. Motion and added resistance responses were measured and predicted for regular sine wave conditions only. Finally, practical considerations regarding full scale implementation of a bow anti-pitch foil such as induced yaw moment, fabrication, installation, vibration etc. are not addressed. The validity of these simplifications, limitations and assumptions are addressed in subsequent chapters. Chapter 2 Historical Background 2.1 Phenomenon! of Added Resistance A d d e d r e s i s t a n c e i s t h e e x t r a r e s i s t a n c e , over and above t h e c a l m w a t e r r e s i s t a n c e , t h a t a s h i p e n c o u n t e r s w h e n u n d e r t h e i n f l u e n c e o f w i n d a n d waves. T h u s a s h i p o p e r a t i n g at a c e r t a i n s p e e d i n w i n d a n d waves w i l l u s u a l l y r e q u i r e m o r e p o w e r t h a n t h e s a m e s h i p o p e r a t i n g at t h e s a m e s p e e d i n p e r f e c t l y c a l m w a t e r . 2.1.1 Causes of Added Resistance T h e a d d e d r e s i s t a n c e e x p e r i e n c e d b y a s h i p i s c a u s e d b y t h e f o l l o w i n g f a c t o r s [35]: 1. W a v e r e f l e c t i o n o r d i f f r a c t i o n o f t h e i n c i d e n t w a v e s y s t e m o n t h e s h i p 2. T h e m o t i o n o f t h e s h i p ( r i g i d b o d y m o t i o n s ) . T h e vertical motions, heaving and p i t c h i n g , h a v e t h e l a r g e s t effect. 3. S i d e s l i p o r d r i f t angle o f t h e s h i p d u e t o t h e w i n d o r w a v e s y s t e m sideforce e x e r t e d on the ship. 4. W i n d a g e o n t h e s u p e r s t r u c t u r e a n d h u l l a n d a u x i l i a r y p r o p u l s i o n d e v i c e ( s u c h as sails). 5. M o t i o n o f t h e c o n t r o l surfaces. R u d d e r a n d o t h e r m o v a b l e a p p e n d a g e s w i l l i n f l u e n c e the m o t i o n of the ship. 9 Chapter 2. Historical Background 10 6. Change i n propeller efficiency because the propellor is operating in a flow due to the motion of the ship. fluctuating O n a sailing yacht the efficiency of the sail driving force may be affected for the same reason. A s one can see, the added resistance experienced by a ship can be the result of many interacting components. This investigation deals only with components (1) and (2) with the motions being restricted to heave and pitch in head seas. Although the other components may be important they are difficult to study and separate. Also, previous researchers have identified items (1) and (2) as the largest contributing factors in the head sea case and the head sea case has been shown to produce the m a x i m u m added resistance [2, 35]. 2.1.2 Mechanisms of Added Resistance in Head Seas Several researchers have examined the mechanisms that lead to the added resistance of a ship when it is heaving and pitching in head seas. Vossers [38] and Hanaoka et al [15] suggested that the following three mechanisms are responsible: 1. Interference between incident waves and waves generated by the ship when heaving and pitching. Sometimes called the drift force. 2. D a m p i n g force associated with the forced heaving and pitching i n calm water; for example under the exitation of some external mechanism such as an oscillator. 3. A component due to the wave system experiencing some wave reflection from the ship's hull. This is sometimes called the diffraction force. A l l three of these mechanisms involve energy transer from the ship to the water. This energy is transmitted to the water and remains i n the waves generated by the ship (except for a negligible amount consumed by viscous friction). In principle, these three wave Chapter 2. Historical Background 11 components can be superimposed. However, the forces developed by these mechanisms will be proportional to the waveheight squared and are therefore nonlinear. Therefore, in reality, these components are nonlinear and interrelated and cannot be separated. Of the three mechanisms, the drift force has been shown to have the largest contribution to added resistance followed by the damping force and finally the least contribution is from the diffraction force [35, 39]. 2.2 Theoretical Techniques for the Predicton of Added Resistance Many formulations have been put forward over the years for predicting added resistance in waves. These follow one of three methods: • Hull Pressure Methods • Momentum and Energy Methods • Radiated Energy Methods In this section, three formulations are considered starting with Havelock [16], followed by Joosen [17] and culminating with the work of Gerritsma and Beukelman [10]. Before developing these theories it is necessary to first fix the coordinate system of a ship in waves for consistancy. Figure 2.2 identifies the variables and sign convention associated with a ship heaving and pitching in regular waves, z is heave displacement, positive up, and 9 is pitch angular displacement, positive bow down. 2.2.1 Havelock Theory One of the first attempts to calculate added resistance of a ship in waves was carried out by Havelock [16] in 1932. He obtained a simple and easy to use formula for the mean added resistance due to a heaving and pitching ship in head seas at any given frequency Chapter 2. Historical Background Figure 2.2: Co-ordinate System 12 Chapter 2. Historical 13 Background of encounter. Havelock's method involved determining the mean value of the longitudinal component of the pressure forces integrated over the wetted part of the oscillating ship's hull. The water pressure was assumed to be the undisturbed pressure of the incident wave, implying the use of the Froude-Kryloff hypothesis. The Havelock expression is: (F z a sin e a zF + M8 a a sin e ) eM (2.1) where k is the wave number, F and M are amplitudes of excitation force and moment a a in heave and pitch and z and 9 are the corresponding heave and pitch amplitudes. e p a and eeM a r e a z the phase angles between the exciting function and the response for heave and pitch respectively. Through use of the Froude-Kryloff hypothesis, Havelock avoided the difficult problem of evaluating the complicated diffracted waves. Thus, this solution was considered a first approximation only. Even though Havelock neglects diffraction effects and pitch-heave cross coupling, the theory still provides some valuable insights into added resistance. Firstly, we see that both the motion amplitudes and their phases are important factors in calculating added resistance. Also, it is possible to note that the maximum added resistance is likely to occur in the region of pitch and heave resonance. Because of this direct pitch and heave motion relationship to added resistance found in the Havelock equation, it can now be stated categorically that hulls with poor motion characteristics will suffer from higher added resistance or conversely an improvement in motion characteristics will lead to less added resistance. This validates the general statement to this effect made in the introduction. 2.2.2 Joosen Theory The foundation of Joosen's theory [17] developed in 1966 was to consider a control volume around the ship and then derive an energy or momentum balance. This approach was Chapter 2. Historical 14 Background first suggested by Maruo [23] in his analysis of the drifting force of a body in waves. The basis of the solution is to subdivide the velocity potential into the incoming wavefield, the diffraction of the wave system by the body and that due to the oscillating body in a regular wave system. The resultant boundary value problems were solved by Maruo. Joosen then expanded Maruo's expression into an asymptotic series with respect to the slenderness ratio L/B and kept only first order terms. Joosen's simple expression for added resistance is given by: R*. = ^{NA + N,0l) (2.2) where N and Ne are the damping coefficients of heave and pitch. This slender body z approximation is valid for short waves only and produced a speed independent added resistance. The theory was extended by Joosen to account for forward speed by substituting the wave encounter frequency w for the wave frequency w in equation 2.2 (though not e entirely mathematically consistant). This equation has been shown to be equivalent to Havelock's equation 2.1 and suffers again from neglecting the pitch-heave cross coupling and the diffraction effect [35]. However, a fundamental concept is inherent in its formulation as well. Equation 2.2 shows that the added resistance can be regarded as a result of the radiated damping waves. That is, the energy loss associated with the damping characteristics of the ship is related to the work that must be done to keep the ship motions in a constant phase relationship with the forcing functions. This energy is supplied by the ship and is dissipated in the damping waves radiating from the hull. Therefore, the added resistance of a ship in waves can be found by computing the total energy content of these radiated damping waves. Chapter 2. 2.2.3 Historical 15 Background Gerritsma and Beukelman Theory The concept of the direct relationship of the added resistance to the energy contained in the damping waves was developed further by Gerritsma and Beukelman [10]. Their approach was to determine the energy radiated from the hull during one wave encounter period and then equate this to the added work being done by the ship during the same period. The energy radiated during one wave encounter period T is given by: e P = [ 'f b(x)V (x,t)dxdt Jo Jo (2.3) 2 z where b(x) is the velocity dependent damping coefficient of the ship at any longitudinal position. b(x) = N{x) - V(dm(x)/dx) (2.4) where m(x) and N(x) are the zero speed sectional added mass and damping coefficients respectively and V is ship velocity. V z in equation 2.3 is the vertical velocity of the ship section relative to the disturbed water surface elevation and is given by: V (x,t) z = z-x 6 b (2.5) + V6-(* £* is the effective vertical water displacement for a cross section. Here Gerritsma and Beukelman have made a dynamic correction to the Froude-Kryloff hypothesis in order to take into account the presence of the ship's hull. This correction takes the form of a modification of the incident wave amplitude by the hull surface and is given by: where £ is the undisturbed wave height, y is the half-width of the waterline and y and w zy, are the offsets of the hull. b Chapter 2. Historical 16 Background If we assume that V is a harmonic function of time of frequency u/ it can be expressed z e as: V = V z cos{u; t + e) za (2.7) e and the time dependence of equation 2.3 can be integrated to yield: The added work of the ship done by the towing force R can be given as [15]: aw P = Raw(V + c)T = R X e (2.9) aw where c is the wave celerity and A is the wavelength. From equations 2.8 and 2.9 we obtain the final form of the Gerritsma and Beukelman equation: R aw = (2.10) ^~ f b(x)V (x)dx 2uj Jo L 2 z a e This equation requires an accurate knowledge of the zero speed added mass and damping coefficient distribution as well as the pitching and heaving motions in order to calculate the added resistance. However the reader will note the elegant and simple derivation which does not require solving any hydrodynamic boundary conditions. This approach attempts to correct for the Froude-Kryloff hypothesis and the pitchheave cross coupling is inherent in its formulation since V is a function of both z and 6. z However other simplifying assumptions are still present. From equation 2.10, it is clear that the added resistance is assumed to vary as waveheight squared since V is proporza tional to waveheight. Also, in head seas, transverse symmetry of all physical processes is assumed. The Gerritsma and Beukelman theory overcomes some of the drawbacks of previous theories in that it attempts to account for the diffraction mechanism of added resistance in very short waves, which Havelock and Joosen neglect. Chapter 2. 2.3 Historical 17 Background Ship Motion Theory It is evident from all of the added resistance theories just discussed that an accurate knowledge of ship heaving and pitching motion amplitude and phases as well as added mass and damping coefficients is necessary in order to calculate added resistance. In this section, we will briefly summarize the classical linear hydrodynamic ship motion theory for the coupled heave-pitch equations of motion which allows us to predict these quantities. 2.3.1 Coupled Heave-Pitch Equations of Motion The first ship motion theory suitable for numerical computations which had an adequate accuracy for engineering applications was the well known strip theory developed by Korvin-Kroukovsky and Jacobs in 1957 [20]. This theory was extended to predict motions in sway, yaw, and roll as well as pitch and heave by Salvesen et al in 1970 [26]. In doing so, improvements to the basic pitch-heave formulation of Korvin-Kroukovsky and Jacobs were realized. Under the assumptions that the motion responses are linear and harmonic, the six linear coupled equations of motion can be written in the form: 6 £ [{M jk + A )rj ik k + B j k V k + C ] jkVk = Fj e "'* 1 j = !,-•-, 6 (2.11) k=i where n are the translatory and angular displacements. k Under the further assumption of ship lateral symmetry, the six coupled equations of motion reduce to two sets of equations, one set of three coupled equations for pitch, heave and surge and another set for yaw, sway and roll. That is, for a ship with lateral symmetry surge, heave and pitch are not coupled to yaw, sway and roll. If one assumes that the ship has a long slender hull form, it can be shown [26] that the hydrodynamic forces associated with the surge motion are very much smaller than the Chapter 2. Historical 18 Background forces associated with heave and pitch. Thus, under these assumptions, the three coupled equations of motion for pitch, heave and surge reduce to two coupled equations for heave and pitch. The heave-pitch coupled equations can be written as: (M + A )z 33 A z 53 + B z 53 + B z 33 + C z + C z 53 33 +B 9 + A 9 3S 35 + C 9 35 +B 6 + C 9 + {h + A )9 55 55 S5 = Fe ° t (2.12) = Fe iUet (2.13) iw 3 z where we have returned to the notation of Figure 2.2 for heave and pitch displacement and the subscript 3 denoting heave and 5 denoting pitch. F and F are the amplitudes 3 5 of the exciting force and moment. 2.3.2 Determination of Hydrodynamic Coefficients The relationships for the added mass and damping coefficients, A,*, and Bjk, and the amplitude of the exciting force and moment, F and F , involves the general hydrodynamic 3 5 problem of solving for the total velocity potential of the fluid. The velocity potential must satisfy the Laplace equation subject to a free surface and hull boundary condition as well as a suitable radiation condition. The method of Salvesen et al [26] is to first linearize the hull and free surface boundary conditions. With the boundary conditions linearized appropriately, it is possible to solve for the velocity potential. The next step is to obtain the hydrodynamic forces and moments acting on the hull using Bernoulli's equation. By solving Bernoulli's equation, an expression can be obtained for the hull added mass and damping coefficients as well as the exciting force and moment. It is possible to reduce this three dimensional potential flow problem to a sum of two dimensional problems, with the result that we need only solve for the two dimensional added mass and damping coefficients and the two dimensional exciting force and moment along the length of the hull. Chapter 2. Historical Background 19 A complete mathematical derivation of the hydrodynamic coefficients and exciting force and moment appears in the Appendix of Salvesen et al [26] if the reader wishes a more detailed treatment. 2.3.3 Assumptions In the context of this study, it is more important to clearly state the assumptions that are made in the derivation rather than examine the derivation itself as these assumptions are very significant in the application of the theory. First of all, in order to make it a potential flow problem, it is assumed that all viscous effects can be disregarded hence the only damping considered is the energy loss associated with creating surface waves. This is valid for heave and pitch motion as the viscous damping is very small in this case [35]. To have a potentialflowproblem, one must also assume incompressible, irrotational flow. In order to linearize the potentialflowproblem, it is necessary to assume that the wave resistance perturbation potential and all its derivatives are small enough to be ignored in the formulation of the motion problem. Physically, this means the free surface waves created by the ship advancing at constant speed in calm water are assumed to have no effect on the motions. This appears to be a reasonable assumption for fine slender hull forms. Finally, in order to reduce the three dimensional problem to a summation of two dimensional problems, it is necessary to assume that the wave frequency is relatively high. This means that the waves created by the ship's oscillations should have a wave length of the order of the ship beam rather than the ship length. This is a critical assumption in that in reality the maximum motion responses are in the fairly low frequency long wavelength range in which the theory appears to be invalid. However, the pitch and heave motions in the low frequency range are dominated by the hydrostatic restoring forces so that these Chapter 2. Historical 20 Background apparent inaccuracies in the hydrodynamic coefficients in the low frequency range only have a minor effect on the predicted motion characteristics. 2.3.4 Added Mass and Damping Coefficients The speed and frequency dependent added mass and damping coefficients derived by Salvesen, Tuck and Faltinsen [26] are: A33 = a d£ (2.14) £33 = J b d£ (2.15) A = 35 A53 J 33 33 ~ j{t-ic)a d(-— B 33 2 r V = - Ht-tc)a dt + 33 B35 = ~ j(l B53 = -fit- Ass = Bss = 33 -ic)bzzdZ —B (2.16) (2.17) 33 + VAaa tc)b d(-VA 33 33 (2.18) (2.19) [(S-tc) a dt + ~A u>l 33 (2.20) ({i-ic) b di + ^B 33 (2.21) 2 33 J 2 33 The reason for the underlined term in equation 2.21 is explained later in Chapter 3. The hydrostatic restoring coefficients are assumed to be independent of frequency and forward speed and are given by: C33 = pgj b(i)di -pgJ(t-t )b(t)d( C7 35 = C C55 = P9J(t-L) b(0dt B3 = (2.22) e 2 (2.23) (2.24) Chapter 2. 2.3.5 Historical 21 Background Exciting Force and Moment The final form of the exciting force and moment expressions derived by Salvesen et al including an empirical assumption suggested by Korvin-Kroukovsky [20] are: F F J f{t)e»«-M J m{()e «-^ 3 = Re{C 5 = Re{( = [P9b(0-^ea a a d^} (2.25) d(} (2.26) ik where fit) (0 2.4 = - ( * - & ) fit) 3 3 -ib )]e- ^ d 3 3 (2.27) k V - — w(a; a e 33 ib ) -«*(0«(0* e (2.28) 33 Experimental Work on Ship Motions and Added Resistance A great deal of experimental work has been carried out in an effort to validate ship motion theoretical models and their assumptions. This section summarizes the experimental techniques used and the conclusions made from experimental work with respect to the added resistance and motions of ships and yachts. 2.4.1 Experimental Techniques Two different methods of measuring model scale added resistance and motions of ships are commonly used in towing tanks. These are [28, 35]: 1. Constant Thrust Method 2. Constant Velocity Method In the constant thrust method, the model is towed by a constant weight and the resultant speed of the model is measured. The average speed is then recorded over an Chapter 2. Historical Background 22 integer number of waves. This enables one to determine the speed loss (as compared to calm water) associated with the constant towing force. Using this method, the model is restricted in roll, yaw, and sway but is free to heave, pitch and surge. In the constant velocity method, the model is firmly attached to the resistance dynanometer and no speed variations are allowed. The model is again allowed to heave and pitch but is restricted to surge in addition to being restricted to yaw, sway and roll. The quantities measured in both methods of testing include resistance, ship speed, wave frequency, wave amplitude, heave and pitch amplitudes and heave and pitch phase. Motion and added resistance tests can be carried out in regular sine waves or in irregular wave spectra. For regular sinusoidal wave tests, in order to obtain a complete frequency response spectrum, the model should be tested for at least ten wavelengths for each speed under consideration. At a particular wavelength the tests should be carried out at a number of different waveheights. Heave and pitch motion response spectra can be obtained directly from the measured heave and pitch motions found at each encounter frequency. The motion spectra are usually presented in non-dimensional form by dividing heave and pitch by wave amplitude and wave slope respectively. Added resistance is presented in the form of the non-dimensional added resistance coefficient given by: a where R aw ' w = pgGVB'L) ( 2 2 9 ) is the average resistance in waves less the calm water resistance at the cor- responding speed and is the measured wave amplitude. Note that the assumption that the added resistance is proportional to wave amplitude squared is inherent in this expression. Chapter 2. Historical 23 Background The regular wave motion responses and added resistance are plotted against the nondimensional frequency of encounter given by: (2.30) where the dimensional frequency of encounter w is related to the wave frequency u> by: e (2.31) This results in completely non-dimensional plots. To obtain the predicted ship response in irregular waves, the measured regular wave motion and added resistance response spectra are applied to a given wave spectrum using the principle of superposition (assuming this is valid). Conversely, the ship response spectrum can be measured directly for one particular irregular wave spectrum by running the model through that spectrum. 2.4.2 Experimental Findings for Ships and Yachts Most experimental research on motions and added resistance has been carried out on cargo ships and naval vessels which tend to have long slender slab-sided hull forms. Tests conducted on these hull forms have resulted in the following findings and conclusions: • Principle of linear superposition is valid for motion and added resistance tests [13]. • Ship motions are proportional to wave amplitude [35]. • Added resistance is proportional to wave amplitude squared [35]. • Scaling of added resistance can be carried out under Froude's Law [2]. • Strip theory of computing added mass and damping is valid [39]. Chapter 2. Historical Background 24 Researchers studying the motions and added resistance of yachts have tried to extend these findings for regular hull forms to include yacht hull forms, however yacht hulls tend to be very non-linear in form causing deviations from the assumptions which are made in strip theory calculations for cargo ship motions and added resistance [8]. These non-linearities include [19]: • Yachts are rarely slab-sided. • L/B and L/D ratios are low in comparison to ships. • Heeled shapes with appendages are difficult to model analytically. • Due to long overhangs and heel, the dynamic waterline may be significant especially in terms of damping and restoring force [12]. In addition, when testing yachts for motions and added resistance there are extra difficulties that are not present when testing ships including [19]: • yachts operate at large heel and leeway angles. • yachts operate in oblique waves. • the vortices shed by a yacht hull as a consequence of sideforce generation may influence the motion response. • sails may influence the damping and the sail driving forces may in turn be affected by the motions due to deformation of the rig and the sails. With these difficulties and limitations in mind, researchers were still able to make some useful conclusions regarding motions and added resistance of yachts. Chapter 2. Historical Background 25 Research on Yachts In 1967 Spens et al [32] laid the groundwork for all future studies on yachts in waves. They studied the motions and added resistance of yachts in regular oblique and head waves resulting in the following findings: • Response in oblique waves is in general agreement with response in head seas of corresponding frequency. • Change in leeway angle due to waves is negligible. Thus they concluded that tests in head waves could be used to analyse a yacht's performance in oblique waves. Gerritsma [8] expanded on the work of Spens, arriving at the following conclusions: • Yacht motions are proportional to wave amplitude. • The waveheight squared law is not always valid and the relationship between added resistance and waveheight can be frequency dependent. • Unlike cargo ship hulls, for yachts the restoring force and moment in pitch and heave increase with forward speed due to its own wave system. • Present linear strip theory methods are not quantitatively accurate when applied to yacht hull forms. • Motions and resistance for 20 degrees heel do not differ substantially from the upright condition values. Pedrick [25] in 1974 analysed the results of yacht tests in oblique waves determining that force testing in oblique waves can be used to make predictions of motions and speed Chapter 2. Historical Background 26 to windward accounting for both added resistance and side force effects in any specified sea condition. A recent study was carried out by Klaka and Penrose [19] in 1987 for the Kookaburra syndicate. They determined through model tests that: • Winged keels have little effect on the motions and added resistance of 12-metre hull forms. • Heel angle had little effect on motions but caused a significant decrease in added resistance (approx. 50%). • Upright yacht testing alone is not adequate for prediction of speed loss due to waves. These final two conclusions disagree with the results of Gerritsma [8] which found that heel angle had only a small effect on added resistance (though for a very different hull form). However, a very recent study just completed by Gerritsma [11] tends to agree with Klaka and Penrose that for heavy displacement yachts, heel angle can affect added resistance. 2.5 Reducing Motions and Added Resistance in Waves To reduce the motions of a vessel in a seaway, motion stabilization devices and methods are often adopted. Using these methods, the amplitude, rate and accelerations of the motions as well as the associated dynamic effects such as deck wetness and slamming can be reduced considerably. Since added resistance is directly related to ship motions, one might expect that motion stabilization may also be effective in reducing added resistance in a seaway. With this in mind, this section will examine the various methods of motion reduction and the experimental findings in this area. Chapter 2. 2.5.1 Historical Background 27 Methods of Motion Reduction With reference to the equations of motion (equations 2.12 and 2.13 in section 2.3.1) one can see that there are three mechanisms by which one can reduce forced motions [2]. 1. By increasing the damping coefficient, (damping stabilization) 2. By reducing the natural frequency of the ship to shift the region of resonance, (tuning stabilization) 3. Reducing the total exciting force or moment by applying a stabilizing moment that is opposite in phase to the wave exciting moment, (equilibrium stabilization) Some of the devices and methods that are used for motion reduction and their associated mechanisms are: • Changing mass distribution (tuning) • Movement of weight (equilibrium) • Passive and active tanks or other vibration absorber (tuning and equilibrium) • Anti-motionfinsor foils (damping) • Gyroscopic stabilizers (equilibrium) • Passive stabilizers such as bilge keels for roll (damping) Most motion reduction techniques have been developed for roll stabilization such as anti-roll tanks, bilge keels, anti-rollfinsetc. This is because the forces and moments involved with rolling are comparatively small. Also, normal vessels have inherently low damping properties with respect to roll resulting in a sharply tuned resonant region. Chapter 2. Historical Background 28 The problem of reducing pitch motions is quite different from that of roll stabilization. The fact that pitch is already heavily damped by the vessel's hull means that very large forces are required for any further increase in damping. The pitch exciting moment is also comparatively large requiring very large forces in order to counteract it. It is quite impractical to generate these large forces by the means of internal devices generally associated with roll stabilization such as vibration absorbers, moving weights or gyros. Thus most attempts at pitch stabilization have involved the use of external devices capable of generating and sustaining the large forces required. The external devices that have been investigated for pitch stabilization include fixed bow foils, movable stern fins and, to a lesser degree, fixed stern fins and controllable bow foils [14]. 2.5.2 Experimental Findings on Anti-Pitch Foils Many studies have been done on the use of anti-pitch foils on large cargo ships and naval vessels such as aircraft carriers or destroyers. Most of this work has been concerned only with reducing the pitching motions of the ship so that it can maintain speed or maintain its ability to operate in a seaway, with little attention being paid to the reduction of added resistance. Abkowitz [1] showed in 1957 that the pitching motion of ships could be reduced to onehalf normal values by the use of reasonable sized anti-pitching fins at the bow, particularly at speeds in the vicinity of synchronism or in irregular seas with synchronous components. Spens in 1962 [31] found that: • Fixed fins were effective at the bow but not at the stern. • Controlled fins are not necessary at the bow. • No change in phase relationship of pitch to waves appears possible. Chapter 2. Historical Background 29 Stefun ixi 1959 [33] carried out experiments to optimize the parameters of fixed bow fin design. He found that: • Fins of approximately 2% of the waterplane area of the hull were most effective for their size. • Fins of aspect ratio 2 showed the best results. • The use of tip fences causes an additional reduction in pitch of 5% over those without fences. Full scale tests have been carried out using anti-pitching fins with mixed results. Although in several cases significant pitch reduction was achieved, severe vibrations caused their removal [3, 14]. Other findings of researchers in this area include: • Deeper submergence of fins, greater fin span and tip fences can help reduce structural loadings caused by vorticity effect. • If fins remain a minimum distance below the water surface (depending on hull form) an air cavity will not form and vibration will not occur. • Maximum pitch reduction can be achieved by fins located in the forward 10% of the ship Stefun [33] was the only one to measure added resistance as well as motions during his anti-pitch fin tests and he found that although the fins caused a 10% to 15% increase in the calm water resistance, this was more than offset by a decrease in motion-induced resistance due to the effect of the anti-pitching fins. Conolly [3] also noticed in full scale trials that the RPM of a ship equipped with foils tended to be lower than that for a normal ship when travelling at the same speed in identical seas. Chapter 2. Historical Background 30 No-one has previously researched the effect of an external appendage such as a bow anti-pitch foil on the motions and added resistance of a yacht hull. The idea of wings on rudders of yachts as an anti-pitch device was mentioned briefly by Cox and Whitaker [4] but has not yet been pursued. The only research that has been carried out on reducing motions and added resistance of yachts in waves was conducted by Spens [32] and Gerritsma and Beukelman [12] and dealt only with changing mass distribution. Both found that decreasing the longitudinal radius of gyration of the hull will improve the seakeeping performance of the yacht. This agrees with the the usual practice of yachtsmen to concentrate the weight of their yachts in the centre of the boat whenever possible. Gerritsma showed further that a light displacement yacht tends to show better seakeeping qualities then a heavy displacement yacht of similar hullform [11, 12]. The results of the research on anti-pitch foils on ships led the author to believe that even though anti-pitching fins did not seem practical for large ships due to structural and vibrational problems, a properly placed and sized anti-pitch foil could be useful for reducing the added resistance and motions of yachts and other small marine craft. Chapter 3 Mathematical Model Formulation From the historical work on ship motions and added resistance theory presented in Chapter 2, a mathematical model was developed to predict the heave and pitch motions and added resistance of the Canada II hull form. This theory was then extended to include the effect of a bow anti-pitch foil. In this chapter, the equations used in the numerical calculations are presented and the results of the predicted motions and added resistance for Canada II with and without the anti-pitch foil are compared. 3.1 Principal Parameters and Limitations The principal particulars for the hull form of Canada II are listed in Table 3.2. The offsets for Canada II were obtained by digitizing the lines drawings which were kindly supplied by the designer Bruce Kirby. It should be noted at this point that Canada IPs hull form, as with most yacht hull forms (see section 2.4.2) deviates from the linear assumptions associated with linear motion theory. The hull is not slab-sided, has a relatively low L/B ratio and possesses long fore and aft overhangs not accounted for in linear theory. In addition, the theory neglects the effect of deck wetness and forefoot emergence on the motions and added resistance. Thus, the validity of the use of the linear motion theory may be suspect if one were expecting accurate quantitative results. However since the object of this theoretical investigation is to detect the relative difference of the hull with the anti-pitch foil and 31 Chapter 3. Mathematical Model Formulation 32 Table 3.2: Canada II Principal Particulars Length of Waterline 45 ft. Maximum Beam on Waterline 11.36 ft. Draft to Waterline (Hull) 4.96 ft. Draft to Waterline (Keel) 9.08 ft. Displacement 59,360 lbs .(S.W.) Length/Beam Ratio 3.96 Length/Draft Ratio 9.07 Beam/Draft Ratio 2.29 Waterplane Area 356.9 sq. ft. Wetted Surface Area 695.0 sq. ft. Vertical Center of Gravity 3.67 ft. below waterline Longitudinal Center of Gravity 2.04 ft. aft of midships without, it was hoped that the linear theory used along with added resistance theory would be sufficient in a qualitative context and therefore would still be useful as a design tool. As shown by Scragg et al [27] of the Stars and Stripes syndicate, even though their wave resistance studies showed disappointing quantitative results, the theory was still very useful in showing the relative difference between candidate designs. In order to make direct comparisons to later experimental results the theoretical scope was constrained to the following: • Zero heel angle only. • Head seas only. • One speed and one waveheight (8 knots and 2.5 ft. full scale) • Neglect effect of sails and rigging. Chapter 3. 3.2 Mathematical Model 33 Formulation Frank Close-Fit Motion Program To predict the heave and pitch motions of Canada II, the Frank Close-Fit Computer Program [7] was used. This program is based on the theory of Salvesen et al [26] presented in section 2.3. The coupled equations for heave and pitch are repeated here as: {M + A )z 33 A z 53 + B z 53 + B z 33 + C z 53 + C z + AJ + B 6 + C 6 = F e <- + (h + A )9 + B 6 + C d = Fe 33 3 55 35 55 35 55 iw t 3 iWet z (3.32) (3.33) with the added mass and damping coefficients being identical to those given in section 2.3.4 except for the expression for the pitch damping coefficient B . The Frank program 5B does not include the underlined forward speed term shown in the expression for f? in B5 equation 2.21 of section 2.3.4. This forward speed correction has been included in this study to be consistant with the Salvesen theory and with the conclusions of Smith [30]. 3.2.1 2-D Added Mass and Damping The Frank program solves for the 2-D sectional added mass and damping coefficients a 33 and 633 using the close-fit method developed by Frank [7] which is a form of boundary element method. The geometric shape of the section is mathematically represented by a given number of offset points with straight line segments between the points. The 2-D velocity potential is then obtained for a distribution of source singularities over the submerged surface of the cylinder with constant strength over each of the segments. These source singularities are such that each one of them satisfies the 2-D Laplace equation, the linearized free surface condition and the infinity radiation conditions. The cylinder wall boundary condition is finally satisfied by solving for the appropriate strength for each of the source segments. This gives the velocity potential for the boundary value Chapter 3. Mathematical Model 34 Formulation problem. By using the linearized Bernoulli equation, the hydrodynamic pressures on the cylinder wall can be obtained. Integration of these pressures over the immersed portion of the cylinder yields the hydrodynamic force and therfore the added mass and damping coefficients. This method is very accurate, however it breaks down at higher frequencies. Thus there is a limiting frequency at which the added mass and damping can no longer be calculated for a section. For a rectangular cylinder with beam b and draft d it can be shown that the first frequency for which the method fails is [7]: /ti = [(7rL/6)coth(7rti/6)] (3.34) 1/2 where p,\ is non-dimensional with respect to ship length L. Thus the program only computes motion results for non-dimensional encounter frequencies less than that given by equation 3.34. pi is directly proportional to ship beam b and for b —* oo the lower limit of pi is \JL/d. Therefore, as a general rule, the close-fit method is applicable for all encounter frequencies; fi 3.2.2 e (3.35) sjlTd < 2-D Exciting Force and Moment The complex 2-D sectional exciting force and moment /(£) and m(£) are also computed by the close-fit method. The resulting expressions differ from those given in equations 2.27 and 2.28 in section 2.3.5 as those equations use an empirical assumption. The Frank program expressions for sectional exciting force and moment are: = m(0 = ) (3.36) + iPv(v)]je ^dr J^ { ~~lPM 2 k l ] - L)f(C) - 2i-£ / /x e Jo [ () Pa V + i ( )} Pv V e k<M dv (3.37) Chapter 3. Mathematical Model 35 Formulation where (77, £) are the sectional co-ordinates with the £-axis positive upward and the 77-axis at the undisturbed free surface. Here p and p are the 2-D hydrodynamic pressures in a v phase with the acceleration and the velocity respectively. This method is more correct than that of equations 2.27 and 2.28 since the vertical variation in the pressures is included very accurately. Equations 2.25 and 2.26 are then used to find F and F5. 3 The geometric parameters for Canada II were input to the Frank Close-Fit Computer Program for a speed of 8 knots. The resultant coefficients and exciting force and moments were then input to a Lotus spreadsheet and equation 2.21 was used to correct the pitch damping coefficient £? for forward speed. The equations of motion (equations 3.32 55 and 3.33) were then solved for pitch and heave. 3.3 3.3.1 Modelling Effect of Anti-Pitch Foil Foil Parameters From the findings of Stefun [33] and others listed in section 2.5.2, a set of near optimal foil parameters was developed for a hypothetical bow anti-pitch foil fitted to the Canada II hull. The relevant parameters of the hypothetical foil are: Table 3.3: Full Scale Anti-Pitch Foil Parameters Foil Span 4.0 ft. Foil Chord 2.0 ft. Depth Below Free Surface 6.0 ft. Initial Incidence Angle 0.0 degrees Distance From LCG of Ship 24.5 ft. Aspect Ratio 2.0 Planform Area 8.0 sq. ft. Chapter 3. Mathematical 3.3.2 Model 36 Formulation Anti-Pitch Foil Theory The problem of a horizontal hydrofoil undergoing oscillatory motion is analogous to the nutter problem of a wing in aerodynamics. If we assume potential flow, it is essentially a problem of solving for the unsteady lift and moment experienced by the foil. To investigate the effect of the hypothetical foil on Canada II, the method of Spens [32] which is based on the work of Kaplan [18] was used. By this method, it is assumed that the foil is a flat plate of large aspect ratio in a potential flow field undergoing linear harmonic oscillations. The theoretical (unsteady motion) hydrodynamic lift L on a hydrofoil of large aspect ratio and its moment M about the center of gravity of the system are, when it is undergoing oscillatory motion in still water [18]: I '1 pirSV 2 V 1 + i TTSV fl 2 M P = i + Note that 4 U -9 z .6 C(«) + Kp-Sc[-V6 '1 (3.38) (3.39) C(K) V V - z + W\ is the steady state damping force of the foil per unit vertical velocity by aerodynamic theory. Although limited to fairly large aspect ratios it may still be considered to be applicable to the case of aspect ratio 2. The term Kp^Sc is the virtual or added mass of the flat plate (the actual mass of the foil is insignificant). When evaluating the lift and moment of the foil attached far from the LCG of the system, the quarter chord c/4 of the foil can be neglected as insignificant when compared to the distance / to the LCG. Equations 3.38 and 3.39 can therefore be reduced to: pirbcV ^ 2 2 L M = = -IL [~V9 -z + ie]+ Kp^-hc 2 [-V6 -z + IB (3.40) (3.41) Chapter 3. Mathematical Model 37 Formulation In these equations C(n) is the Theodorsen circulation function for a two dimensional oscillating hydrofoil evaluated from tables in Luke and Dengler [22] and K is the coefficient of accession to inertia from Munk [24]. 3.3.3 Modification of Coefficients in Equation of Motion By separating the terms associated with z,9 and their derivatives in equations 3.40 and 3.41 one can obtain expressions that must be added to the coefficients in the heavepitch equations of motion. The increments that must be added due to the foil are: = Kp-bc 2 (3.42) 55 = [AA )l (3.43) AA35 = AB = AB = AB l AB 35 = -AB l AB 53 = -AB l (3.48) AC35 = AB V (3.49) AC 55 = -AB Vl (3.50) AC33 = AC AA33 AA 33 55 2 33 AA = 53 -(AA33)/ (3.44) (3.45) 2 33 33 - (3.46) AA V (3.47) 33 + 33 33 33 53 AA Vl = 0 33 (3.51) These account for the effect of the system oscillating in still water. A change in lift and moment due to the wave orbital motion acting on the foil as the system moves through regular head waves must also be included. These result in increments to the right hand Chapter 3. Mathematical Model 38 Formulation side of the equations of motion as follows [31]: AL = TrSV*\w \(l W AM P ffi 0 - ia.)e**' *> j + J o { y \ ~ = —f(AL„,) t7(«) + itf ( l + ^ ) J X H | (3.52) where 1*1 = <. ( ^ ) \w \ is the ratio of the vertical component of the orbital velocity of the wave to the velocity 0 along the x-axis. The increments to the motion coefficients and the exciting force and moment were calculated for Canada II with an anti-pitch foil, for the foil parameters stated in Table 3.3 for a speed of 8 knots into head seas. These calculations were integrated into the existing Lotus spreadsheet and the modified equations of motion were solved for heave and pitch for the hull with foil. 3.4 Added Resistance Calculation After reviewing the literature on added resistance theories, it was decided to use a modified version of the Joosen equation presented in section 2.2.2. Joosen extended his work on equation 2.2 to account for pitch-heave cross coupling [17]. It should be restated that this theory assumes a slender body approximation for the drift forces and is valid only for short waves. Also, the wave diffraction effect is neglected; therefore the theory also breaks down for very short waves where the diffraction effect becomes dominant. Since in this investigation we are concerned mostly with the resonant region, these assumptions should not present a problem. The final form of the modified Joosen equation is [35]: <r aw = E + E + E 1 2 3 (3.53) Chapter 3. Mathematical Model 39 Formulation where E = C B7s(zJCa) E = C {2irL/X) BT (d /kC ) E = 1 (3.54) 2 0 2 2 3 0 5 -2C {2TT 0 a (3.55) 2 a L/X)B7 (z /Ca){e a/Ha) cose 5 a ° = nFWf)' ' c - v/i e = (3.56) (3 57) (3.58) \e -ee\ z The damping coefficients appearing in the above equations are given by: Bis = (l^VKs/L) " Jb U)d( (3.59) 1 33 5s7 = 5^ = (ll^){glLfl 2 {\l^){glLfl J{(-( )b {()di (3.60) j ^ - i ^ ^ d i (3.61) c 3 33 The Joosen method was chosen over the Gerritsma and Beukelman method (see section 2.2.3) for its ease of implementation. The Gerritsma and Beukelman method requires an accurate knowledge of the distribution of the sectional added mass and damping while the Joosen method requires only thefinaltotal damping coefficients of the hull, which are more readily available. Again, since only a relative measure is required in this study, it was decided that though the Gerritsma and Beukelman method might be more accurate (due to its more accurate handling of diffraction effects and cross coupling) its increased complexity did not warrant its use. It is recommended that if more accurate predictions of added resistance are desired in future investigations, one should consider using the Gerritsma and Beukelman method. The predicted added resistance for Canada II with and without the foil was calculated using equation 3.53 using the spreadsheet solutions to the equations of motion and the calculated values of the damping coefficients. Chapter 3. 3.5 Mathematical Model 40 Formulation Results of Theoretical Calculations In this section, the predicted heave motion and pitch motion responses and added resistance resulting from the theoretical calculations are presented for Canada II with and without the bow anti-pitch foil. 3.5.1 Heave Motion The calculated heave response for Canada II is shown in Figure 3.3. Examining the bare hull results first, one notices that the heave response follows the wave amplitude at low frequencies of encounter as one would expect, then begins to rise rapidly as p increases e into the resonant region. Heave resonance for the bare hull is predicted to occur at p = 3.2 which corresponds to a wavelength A of 77 ft. or 1A1L for Canada II at a speed e of 8 knots. At resonance, the nondimensional heave response peaks at 1.27 indicating that heave is quite heavily damped. As p increases beyond resonance the heave response e drops rapidly. Unfortunately, the close-fit method for calculating added mass and damping breaks down for Canada IPs hull form at fi = 3.6 (see section 3.2.1). Hence, no calculations e can be made beyond this frequency. This corresponds to A = 65 ft. or 1.45L which is still a relatively long wave, however it is safe to assume that the heave response will decrease monotonically beyond the cutoff frequency since at very high frequencies the heave response will tend toward zero. If one examines the heave response of Canada II with the anti-pitch foil attached, one finds that the response is quite different from that for the bare hull. For p < 2.7 the e heave response with the foil is predicted to be slightly higher. Above p = 2.7, near the e resonant region, the response with foil does not rise as rapidly as the bare hull, reaching a peak of only 1.16 at /z = 3.0. The resonant frequency of p = 3.0 is slightly lower than e e Chapter 3. Mathematical Model Formulation 41 Theoretical Heave Response 1.4 1.2 CD T5 1- "75. E D CP > D 0.8- CD _C "6 c o "to c a> E 0.6- 'TD 1 1 co z 0.4- 0.2 Legend • Strip T h e o r y - n o foil • Strip T h e o r y - f o i l I 1 I 2 _ ™ ! 3 I 4 ~ N o n - d i m e n s i o n a l encounter frequency Figure 3.3: Theoretical Heave Amplitude Response Chapter 3. Mathematical Model 42 Formulation that of the bare hull at fi = 3.2. Above resonance the response can again be expected to e drop monotonically. In the resonant region and beyond the predicted decrease in heave response as a result of the foil is up to 11% which is quite significant. Though at first one might be surprised to find that for p < 2.7 the heave response e with the foil is actually higher than for the bare hull, it can be explained as follows. If one examines the changes in the coefficients of motion from the foil (see section 3.3.3) one sees that AA , 33 the added mass of the foil, is a constant while AB , 33 the damping due to the foil, is a function of C(«) which is a function of frequency and vertical velocity. At the lower frequencies of encounter the magnitude of C(K) is quite small with the result that the added mass term dominates. This results in a greater heave response since extra added mass reduces the damping ratio. As p increases beyond 2.7, C(K) increases and e the damping term starts to dominate the heave response causing the response to cross over the bare hull response. The slightly lower resonant frequency of the hull with foil can be attributed to the increase in total added mass due to the foil. Since the restoring coefficient is assumed to be constant, the natural frequency would be expected to decrease. 3.5.2 Pitch Motion Referring to Figure 3.4, one sees the predicted pitch response. At first glance one notices immediately a significant reduction in the non-dimensional pitch amplitude for the hull with foil compared to the bare hull throughout the frequency spectrum. As with heave, the pitch response is predicted to rise for the low encounter frequencies up to the resonant frequency. However, unlike heave, the rise is much more gradual and results in a less pronounced resonant hump. Again, we see a slight decrease in the resonant frequency with the addition of the foil, shifting from p = 3.0 for the bare hull to p = 2.8 for the e hull with foil. e Chapter 3. Mathematical Model Formulation 43 Theoretical Pitch Response 1.4 1.2- 0.8 0.6 0.4 0.2 Legend • Strip T h e o r y — n o foil • Strip T h e o r y — f o i l 1 2 3 4 Non-dimensional encounter frequency Figure 3.4: Theoretical Pitch Amplitude Response Chapter 3. Mathematical Model 44 Formulation For the bare hull, the peak non-dimensional pitch amplitude at resonance is 1.33. With the addition of the foil this peak drops to 1.14 at resonance representing over a 15% predicted decrease in pitch amplitude. Even at the lower encounter frequencies the foil is predicted to decrease pitch amplitudes by 10% implying that, in theory at least, the anti-pitch foil is expected to live up to its name. Added Resistance 3.5.3 The predicted added resistance of Canada II, with and without anti-pitch foil, is shown in Figure 3.5. For the bare hull the non-dimensional added resistance coefficient tr , is au predicted to be zero for very long waves (low encounter frequency) then rises rapidly to a peak of a aw = 3.4 at a resonant frequency of u = 3.4. Due to the break down of e the close-fit method for the ship motions when p, > 3.6, the added resistance cannot be e calculate beyond this frequency either, since added resistance is calculated directly from the heave and pitch motions. Note, however, that the added resistance drops quickly between p = 3.4 and p = 3.6 and can be expected to continue to drop. Note also e e that the added resistance curve has a narrow and steep resonant peak which corresponds closely to the heave and pitch resonant frequencies, as expected. Turning now to the added resistance response of the hull with foil attached, one sees that the resonant frequency is essentially the same and a aw is predicted to drop significantly at this point from 3.4 for the bare hull to 2.35 for the hull with foil. This represents a 31% predicted reduction in added resistance due to the addition of the foil. If one examines the modified Joosen equation (equation 3.53) closely, this large reduction is predicted. The additional damping provided by the foil does not affect the values of the purely hydrodynamic damping coefficients B , B 5 and B in the expressions for 33 5 35 Ei, E and E . However, the effect of the foil is to reduce the predicted heave and pitch 2 3 motion amplitude significantly as we have seen. Since the expressions for Ei,E 2 and E 3 Chapter 3. Mathematical Model Formulation Theoretical Added Resistance Response Spectrum 4 3.5H Non-dimensional encounter frequency Figure 3.5: Theoretical Added Resistance Chapter 3. Mathematical Model Formulation 46 are second order in motion amplitudes the reduced heave and pitch amplitudes translate into an even higher percentage reduction in added resistance. It should be noted that in the calculation of added resistance with the foil, the extra hydrodynamic drag of the foil due to the lift it creates has not been included. However it is expected that this extra drag is a small percentage of the total added resistance (especially at resonance). Chapter 4 Experimental Program In this chapter, the anti-pitch foil experimental program carried out by the author is described including a description of the test facilties, model preparation, method of testing, test program, observations and results. 4.1 4.1.1 Experimental Set-up Facilities All experiments were conducted at the B.C. Research Ocean Engineering Centre towing tank. The tank measures 67.0m x 3.7m x 2.4m and the hydraulic drive towing carriage is capable of speeds up to 4.5 m/s. A hydraullically operated, computer controlled flap type wavemaker is installed at one end of the tank and is capable of producing regular or irregular waves up to a height of 0.25m or just over 10 ins. for wave frequencies between 0.3 and 1.5 Hertz (see Figure 4.6). Models are attached to the carriage and towed by a single heave post arrangement as shown in Figure 4.7. Resistance is measured by a shear beam load cell dynamometer attached at the bottom of the heave post and is accurate to ±0.05 lbs. Average speed of the carriage for the run is measured using a digital encoder wheel. Pitch angle is measured directly from the pivot point of the heave post by an angular deflection potentiometer with an accuracy of ±0.05 degrees. Heave displacement is measured using a rolling wheel potentiometer fixed to the carriage and resting against the free-to-heave post. Heave accuracy is ±0.1 inches. 47 Chapter 4. Experimental Program Figure 4 . 6 : Wave Paddle 48 Figure 4.7: Heave Post Arrangement Chapter 4. Experimental Program 50 Waveheight is measured using a Keiseko electromechanical water level meter with an accuracy of ±0.01 inches. The wave probe is mounted on a cantilever bridge that extends to near the midbreadth of the tank located near the midpoint of the length of the tank. Figure 4.8 shows the waveprobe set-up. The data acquisition system consists mainly of a ST41B signal conditioner and a Digital Electronics MINC 11 mini-computer with analogue to digital module. The signal conditioner provides the required excitation voltage to the measurement transducers as well as amplification and filtering of the returned signal. The MINC 11 computer contains an analogue to digital conversion unit for storing the test data in a binary digital format. Figure 4.9 shows the data acquisition equipment on the towing carriage. 4.1.2 Model Parameters A wooden model of Canada II shown in Figure 4.10 was kindly made available by the designer Bruce Kirby. The model was constucted to 1:8 scale according to a lines plan supplied by the designer. The model was fitted with a lead ballast keel for proper location of the VCG. The principal dimensions for the model are given in the Table 4.4. It was attempted to place the heave post attachment point at the VCG of the model; however this was not possible as the VCG is actually in the keel. As a result, the heave post pivot point was placed 1 inch above the VCG which was as close as was practical. It was felt that this would not significantly affect the motions of the model. Turbulence stimulation was provided on the hull by cylindrical studs of 0.125 inch diameter and 0.1 inch high spaced 1 inch apart on station 2. On the keel, a similar row of studs were installed 0.75 inches aft of the leading edge [34]. Chapter 4. Experimental Program Figure 4.8: Wave Probe Set-up 51 Chapter 4. Experimental Program Figure 4.9: Data Acquisition Equipment on Carriage 52 Figure 4.10: Canada II Model Chapter 4. Experimental Program 54 Table 4.4: Canada II 1:8 Scale Model Particulars Length of Waterline 5.6 ft. Maximum Beam on Waterline 1.42 ft. Draft to Waterline (Hull) 0.62 ft. Draft to Waterline (Keel) 1.14 ft. Displacement 113.04 lbs.(F.W.) Length/Beam Ratio 3.96 Length/Draft Ratio 9.07 Beam/Draft Ratio 2.29 Waterplane Area 5.58 sq. ft. Wetted Surface Area 10.86 sq. ft. Vertical Center of Gravity 5.5 inches below waterline Longitudinal Center of Gravity 3.06 inches aft of midships 4.1.3 Model Foil Design A 1:8 scale model anti-pitch foil was designed and built corresponding to the full scale foil parameters used in the theoretical investigation (see Table 3.3). A NACA 0012 foil section of rectangular planform was used for the main foil because of its high lift to drag ratio and relatively high stall angle [32]. Originally, it was planned to attach the foil to the hull using a single non-surface piercing strut attached just below the waterline, aft of the forward perpendicular. However, it was determined that with this arrangement it would be too cumbersome to take the foil on and off and would require a through-hole in the model, which was undesirable. Instead, a double surface piercing strut arrangement was designed that could easily be attached and removed with the model in the water and required only three screwholes in the deck of the model. The model foil assembly drawings are shown in Appendix B. Note that the struts are foil shaped to reduce drag. Also, tip fences were attached to the end of the foils to reduce end effects as suggested by Stefun [33]. Figure 4.11 shows the Canada II model as it was for testing with the foil attached. An Chapter 4. Experimental 55 Program underwater sideview of the model with foil attached appears in Figure 4.12. The model foil parameters are listed in the Table 4.5. Table 4.5: Model Anti-Pitch Foil Parameters Foil Span 6.0 inches Foil Chord 3.0 inches Depth Below Free Surface 9.0 inches Initial Incidence Angle 0.0 degrees Distance From LCG of Ship 3.06 ft. Aspect Ratio 2.0 Planform Area 18.0 sq. inches 4.2 Test Program The model test program consisted of calm water resistance tests, calm water natural frequency tests and regular sinusoidal head wave tests. 4.2.1 Calm Water Tests In order to obtain the experimental added resistance it is necessary to first determine the calm water resistance since: Raw — Rwave Realm Thus calm water resistance tests were conducted for the Canada II model for both the bare hull and the hull with anti-pitch foil attached. In addition to determining the calm water resistance at the target speed of 8 knots, it was of interest to see the effect of the foil on the calm water resistance at all speeds. Therefore the tests were done for model speeds corresponding to 6-9 knots full scale. In order to obtain the natural frequencies in heave and pitch for the Canada II hull, forced oscillation tests were conducted on the bare hull. This involved forcing the model Chapter 4. Experimental Program Figure 4.11: Model Anti-Pitch Foil Fitted For Testing 56 Chapter 4. Experimental Program Figure 4.12: Underwater Sideview of Model Anti-Pitch Foil 57 Chapter 4. Experimental 58 Program to a large heave or pitch displacement then releasing it and measuring the period of oscillation. These tests were done at zero speed and at a model speed corresponding to 8 knots in order to determine if forward speed affected the heave and pitch natural frequencies. 4.2.2 Regular W a v e Tests Before carrying out regular wave tests, the wave paddle had to be calibrated so that a reasonably constant waveheight could be produced at any frequency in the frequency range. All regular wave tests were carried out at a constant model speed corresponding to 8 knots full scale and at a constant waveheight corresponding to 2.5 feet full scale. These were chosen by the author as realistic conditions encountered in Fremantle. All tests were done using the constant velocity method of wave testing described in section 2.4.1. Using this method the model is free to pitch and heave but restricted in surge, sway, yaw, and roll. Regular wave tests were completed for the bare hull and the hull with foil for 15 regular spaced non-dimensional encounter frequency intervals of 0.25 between p = 1.25 e and p = 4.75. Calibration of all equipment was carried out before each day of testing. e All runs were videotaped for later examination of motions. The procedure for the each run of the wave tests was as follows: 1. Take transducer zero values. 2. Start wavemaker. 3. Allow 3-4 waves to pass waveprobe then begin measuring waveheight. 4. Start model just before waves get to it. 5. When model is at speed, begin data acquisition. Chapter 4. Experimental 59 Program 6. Stop waveheight measurements before model passes waveprobe. 7. Stop data acquisition before decelerating model. This test procedure is designed to maximize run length and minimize the starting response transient damping time [6]. 4.3 Observations and Data Reduction 4.3.1 Test Observations Several difficulties surfaced during the course of testing. Firstly, it proved to be very difficult to maintain a constant waveheight over all of the runs. The waveheight seemed very sensitive to small changes in the tank depth. As a result, waveheights were found to vary over the course of testing by about ±0.3 inches from the target waveheight of 3.75 inches. It was observed that at very low frequencies, the waves did not have a sinusoidal waveform. This was probably due to the small depth of the wave paddle flap. As a result, no data could be gathered for non-dimensional encounter frequencies of less than 1.25 which corresponds to a full scale wavelength of over 300 ft. This a relatively long wave therefore well out of the resonant range and should not limit the usefulness of the data. At certain frequencies, it was observed that the water behind the wave paddle would slosh around violently. When this occurred, water would escape around the sides of the wave flap into the main tank and adversly affect the wave form at these frequencies. This problem was most apparent at the wave paddle frequencies corresponding to p = 4.25. e During the calm water natural frequency tests, it was observed that the hull tended to damp out oscillations after only 1 to 2 cycles indicating very heavy damping. Several interesting observations were noted on the models motions during the wave Chapter 4. Experimental Program 60 tests from the video footage. Firstly, it was visually clear that there was a reduction in pitch motions with the addition of the foil especially at resonance where motion amplitudes were greatest. It also appeared that the phase of the pitch motion with respect to the waves was different when the foil was present than when it was not. Any difference in heave motions was not clearly apparent from the video. At p = 3.0 the pitching motion was observed to be especially large and violent, e especially for the bare hull. In fact, the bare hull frequently submerged its bow into the crest of the next wave at this frequency causing green water to be shipped over the deck. Figure 4.13 is a photograph of the bare hull at p — 3.0. Notice the bow underwater and the splash rail that had to be fitted to divert the water running over the deck. It was also observed at these frequncies that the flat aft counter tended to slap down on the water. Bow submergence was not so evident in the tests of the model with the foil attached however did occur from time to time. Figure 4.14 is a photograph of a test at p — 3.0 with the foil attached. 4.3.2 Data Reduction As mentioned before, the quantities measured during testing in waves were model speed, waveheight, pitch angle, heave displacement and resistance force. Pitch, heave and resistance data was collected on three channels at a sampling frequency of 40 samples per second for an average run time of 15 seconds. This data was multiplexed and stored as one binary data file per run. Speed was recorded directly off the carriage drive console and was assumed to be constant over the length of the run for which data was gathered. Wave height was measured separately, again at 40 samples per second for at least ten wave cycles before the model passed the wave probe location. The raw data files for heave, pitch and resistance were then demultiplexed into separate files. These files, along with the waveheight datafile, were then visually edited using a Chapter 4. Experimental Program Figure 4.13: Bow Immersion of Canada II Bare Hull 61 Chapter 4. Experimental Program Figure 4.14: Bow Immersion of Canada II H u l l W i t h A n t i - P i t c h Foil 62 Chapter 4. Experimental 63 Program program developed by the author. Each file was inspected and a window of at least 10 cycles were visually selected that most closely resembled steady state harmonic oscillation. From the selected window for resistance, the average value of the resistance over an integer number of oscillations within the window was calculated to find R w a v e for each run. From the selected windows for heave, pitch and waveheight the average peak-to-peak value was calculated over the number of steady state oscillations within the window giving Zp-p, 9 -p and £ _ for each run. p p p Added resistance for each run was then calculated as: where Realm is measured from the calm water resistance results at the model speed of the wave run. The non-dimensional added resistance coefficient t r is then calculated using: a u ) 4R <Tntn ~ (4.63) pgQ-pB'L Note that the added resistance is assumed to vary as waveheight squared over the range of waveheights of the tests. The average heave peak-to-peak value for each run was non-dimensionalized by the waveheight: Non-dimensional heave amplitude = (4.64) CP-P and the average pitch peak-to-peak value for each run is nondimensionafized by wave slope: Non-dimensional pitch amplitude = -— ~ p p (4.65) (U /9KP-P 2 where OJ is the wave frequency measured from the wave data. The results were tabulated for a total of 30 runs, 15 for bare hull and a corresponding 15 for hull with foil and appear in full in Appendix C. The added resistance coefficient, Chapter 4. Experimental Program 64 non-dimensional heave amplitude and non-dimensional pitch amplitude were then plotted against non-dimensional frequency of encounter (see section 2.4.1). 4.4 4.4.1 Experimental Results Calm Water Tests Figure 4.15 shows a plot of the calm water resistance for the Canada II model with and without the anti-pitch foil for model speeds corresponding to 6 to 9 knots full scale speed. Notice that the the foil causes a significant increase in the calm water resistance at all speeds tested, as would be expected. At the speed of interest for the wave tests, around 4.8 ft/sec or 8 knots full scale, the increase in resistance due to the foil is 25%. It is thought that this rather larger than expected increase is due mostly to the two surface piercing struts rather than the foil itself, though misalignment of the foil may also have contributed. The wave resistance of the struts may be quite significant. The original idea of a single non-surface piercing strut would probably result in a smaller increase in the calm water resistance due to the foil assembly. This difference in calm water resistance does not affect the added resistance as the calm water effect is subtracted out as shown by equation 4.62. However in terms of the overall performance of the boat, it is desirable to have the lowest calm water resistance possible. Table 4.6 presents the results of the heave and pitch natural frequency tests for the bare hull. In both types of motion, one sees that the non-dimensional natural frequency, p , increases with an increase in speed from 0 to 8 knots. This shows that the natural e frequencies of heave and pitch are strongly influenced by speed. It is probable that the restoring force and moment increase with forward speed due to the ship's own wave system as a result of the non-linear form of the hull. Linear motion theory assumes that the restoring coefficients in the equation of motion are not speed dependent and this has Chapter 4. Experimental Program 65 Experimental Calm Water Resistance 0.5H -I 3.5 1 1 — 1 5 Model s p e e d ( f t / s e c ) 4 4.5 1 5.5 ! *~ Figure 4.15: C a l m Water M o d e l Resistance For C a n a d a II 6 Chapter 4. Experimental 66 Program Table 4.6: Calm Water Pitch and Heave Natural Frequencies Heave Pitch 0 Knots 8 knots 2.80 3.15 3.16 3.55 been shown to be a valid assumption for regular cargo ship hull forms [8]. However it appears that some sort of speed correction to the restoring coefficients is necessary for yacht hull forms. 4.4.2 Heave M o t i o n Figure 4.16 shows the non-dimensional heave amplitude plotted against the non-dimensional frequency of encounter for the Canada II model without and with the foil. The resonant frequency in heave is at p e = 3.16 from Table 4.6. For the bare hull the heave amplitude rises slowly from 1.0 to a maximum of 1.12 at u = 3.01 then begins to drop e steeply to 0.14 at p = 4.7. e The heave amplitude response of the hull with foil behaves quite differently. It rises quite rapidly from 1.0 at p = 1.0 to a peak of 1.13 at p = 2.25 which is well below the e e resonant frequency. It then starts to drop and crosses the bare hull response at fi = 2.8. e It then continues to drop steeply to 0.1 at p = 4.7. e Before the crossover frequency the heave response of the hull with foil is up to 5% higher than that of the bare hull. Above the crossover frequency of 2.8 the response of the hull with foil is 15-20% lower than the bare hull response. Chapter 4. Experimental Program 67 Experimental Heave Response 1.4-1 1.2 0> TJ 3 H "D. E o a> > D 0.8 QJ JC "o c o 'v) 0.6- c E ~o I c O 0.4- 0.2 Legend • E x p e r i m e n t - n o foil • Experiment—foil 1 2 3 4 Non-dimensional encounter frequency Figure 4.16: Experimental Heave Amplitude Response Chapter 4. Experimental 68 Program Pitch Motion 4.4.3 Figure 4.17 shows a plot of non-dimensional pitch amplitude from the experiments. The resonant frequency according to the calm water natural frequency tests is p = 3.55. e Looking at the bare hull results, one sees that the pitch response at the lowest tested encounter frequency is just over 1.0. The response then unexpectedly peaks at 1.22 for p = 1.5 however then drops again and stays relatively level at around 1.1 up to ii = 2.75. e The pitch amplitude response then drops rapidly from 1.13 at p = 2.75 to 0.125 at e p = e 4.75. For the hull with foil attached, one sees that the pitch amplitude response again is just over 1.0 for the lowest tested frequency. It then peaks at a level of 1.14 for fi = 1.5. e For p > 1.5 the pitch amplitude drops steadily to a value of 0.1 for p = 4.68. e e Notice that throughout the frequency range tested, the pitch amplitude response is lower when the foil is fitted than for the bare hull. At lower frequencies this reduction is in the order of 5% while in the resonant region near p = 2.75 and above the reduction e in pitch amplitude due to the foil is higher than 20%. 4.4.4 Added Resistance Figure 4.18 shows the measured added resistance from the experiments. For the bare hull at low frequencies, <x, is very small. At p = 2.25 the added resistance begins to rise au e rapidly to a peak of cr = 3.75 at p aw e = 3.66. Above this frequency o a w drops rapidly but then hits a secondary local peak of 3.07 at p = 4.21. After this, the added resistance e again decreases with increasing fi . e Looking at the results of the tests of the hull with foil, one sees that at the lowest frequency o~ is very similar to that of the bare hull. As encounter frequency increases to aw fi = 2.25, the added resistance of the hull with foil is slightly less than that of the bare e Chapter 4. Experimental Program 69 Experimental Pitch Response 1.4 1.2 TJ 1 "Q. E D -g 0.8 D C o £ 0.6 E TJ o 0.4- Legend B Experiment—no foil 0.2 • Experiment—foil 1 2 -r 3 T 4 Non-dimensional encounter frequency Figure 4.17: Experimental Pitch Amplitude Response Chapter 4. Experimental Program 70 Experimental Added Resistance 4 I 0 I I 1 1 1 1 2 3 4 5 N o n - d i m e n s i o n a l encounter frequency Figure 4.18: Experimental Added Resistance Chapter 4. Experimental hull. At p = 2.25, cr e aw 71 Program begins to rise steeply to a much lower peak than the bare hull of 2.3 at ti = 3.5. It then drops and again rises unexpectedly to a second peak of 2.3 at e f i — 4.25. Above tt = 4.25 the added resistance drops rapidly with increasing f i . In the e e e region between p = 3.25 and tt = 4.0 the reduction in added resistance due to the foil e e is up to 37%. Notice that the region of maximum added resistance is in the same region of resonance for pitch and heave as one would have expected. The appearance of a second peak in the added resistance results at p = 4.25 might e be connected to the wave form at this frequency. Recall from section 4.4 that near this frequency it was observed that the wave form was adversly affected by the sloshing water behind the wave flap spilling around the flap into the main tank. This might be the cause of the extra resistance on the model. Chapter 5 Comparison of Theoretical and Experimental Results In this chapter, the theoretical results calculated in Chapter 3 and the experimental results presented in Chapter 4 are compared and discussed. This is done to determine the applicability of the mathematical model to the problem of predicting the effect of an anti-pitch foil on the motions and added resistance of Canada II. The validity of the linear strip theory in predicting the heave and pitch motions of the bare hull is discussed first. This is followed by a discussion of the effectiveness of the anti-pitch foil theory in predicting the reduction in motions. Finally, the added resistance results are compared to determine if the Joosen theory gives a useful estimate of the reduction in added resistance that can be expected from the foil. 5.1 Ship Motion Comparison For Bare Hull Before dealing with the anti-pitch foil at all, it is useful to examine if the strip theory model for ship motions agrees with the experimental results for the Canada II hull itself. Comparing the heave motion results for the bare hull, shown in Figure 5.19, one finds that for /x < 2.0 the theory matches the experimental results very well. For u > 3.6 the e e theory can also be said to match experiment reasonably, if one extrapolates the theoretical results to continue to drop monotonically. In the mid-frequency range between u — 2.0 e and a = 3.6, however, the theory overpredicts the heave amplitude response significantly. e Strip theory predicts a pronounced resonant hump at u = 3.2 while the experimental e resonant hump is barely noticable. 72 Chapter 5. Comparison of Theoretical and Experimental Results 73 Heave Response For Bare Hull Theory vs. Experiment 1.4-| I 0 I 1 1 2 1 3 1 4 r 5 N o n - d i m e n s i o n a l encounter frequency Figure 5.19: Heave Amplitude Response For Bare Hull.Theory vs. Experiment Chapter 5. Comparison of Theoretical and Experimental Results 74 This discrepency indicates to the author that the Canada II hull is more heavily damped than the theory is capable of predicting, resulting in much lower heave amplitude values at resonance. It is thought that this extra damping is a result of the non-linearities in the Canada II hull form. In calculating the damping coefficient, B33, in the equation of motion, strip theory assumes (among other things) that the sections are slab-sided at the waterline and the hull does not extend out beyond the fore and aft perpendiculars. The Canada II hull form violates these assumptions. Firstly, there is considerable flare to the topsides of the hull, leading to a dramatic change in the section shapes as the hull heaves. In addition, Canada II has extremely long fore and aft overhangs that are not accounted for by the theory. When the hull is heaving, these overhangs become immersed causing large fluctuations in the damping distribution over the hull. The effect of these non-linearities is small at low frequencies as the heave motion with respect to the wave is also small; therefore not as much of the overhangs or flared topsides would be immersed. However, in the resonant frequency range where the heave motions, and therefore the vertical velocities of the hull, are at their greatest, a large portion of the flared topsides and overhangs would become immersed. Recall from the calm water heave natural frequency test, that the model was observed to be very heavily damped in heave. Comparing now the pitch motion of the bare hull, shown in Figure 5.20; with the exception of one point at p = 1.5, the theory tends to overpredict the pitch amplitude e response throughout the frequency spectrum, especially at the resonant frequency and above. The strip theory calculation predicts a steadily rising response in the low frequency range up to resonance, while the experiment shows the pitch amplitude response to be relatively flat in this region. It appears from this comparison that strip theory has severely Chapter 5. Comparison of Theoretical and Experimental Results 75 Pitch Response For Bare Hull Theory vs. Experiment 1.4 1.2- 0.8- 0.6 0.4- Legend • E x p e r i m e n t - n o foil • Strip T h e o r y - n o foil 0.2 — 1 2 3 4 5 N o n - d i m e n s i o n a l encounter frequency Figure 5.20: Pitch Amplitude Response For Bare Hull: Theory vs. Experiment Chapter 5. Comparison of Theoretical and Experimental Results 76 underpredicted the pitch damping, even more so than for heave. Once again, the large overhangs of the Canada II hull are suspected of being responsible for the difference in pitch motion response as the strip theory does not include their effect on the pitch damping coefficient # 5 5 . When the hull pitches, the fore and aft overhangs will alternately immerse, changing the underwater shape of the hull and moving the distribution of buoyancy fore and aft. The damping and added mass distribution will also be continually changing with more of it being found in the ends of the ship, well away from the LCG. When the pitch motion amplitudes are large, and the overhangs are immersed even further, the added damping will be increased substantially resulting in the pitching motion being damped out. This proposal of increased hydrodynamic damping due to the overhangs can be justified by the fact that damping waves were observed to radiate out from the overhangs when they were immersed during the wave tests and the calm water natural frequency test. As has been discussed earlier, these radiating waves indicate additional damping. Note also that at resonance, the bow was observed to fully immerse and water was shipped on deck. This would contribute even further damping at resonance. This helps to explain the experimental results. Recall again that the calm water natural frequency tests showed that pitch was very heavily damped. In fact, it was difficult to measure the natural pitch period as the pitch motion was damped out within 1-2 cycles. The calm water natural frequency tests also show another interesting result in that the natural frequency of heave and pitch increased with forward speed. This indicates that the restoring coefficients C 3 3 and C 5 5 are speed dependent. This increase is unexpected and not accounted for by strip theory, which assumes that the restoring coefficients are speed and frequency independent. The added masses, A and A , increase with speed as shown 3 3 55 in section 2.3.4; therefore according to strip theory the natural frequency should decrease with speed. The increase in the natural frequncy shown by experiment is probably due Chapter 5. Comparison of Theoretical and Experimental Results 77 to the ship's own wave system effect on the restoring moment and force [12]. The conclusion of this comparison and discussion of experimental and strip theory results for ship motion is that the strip theory does not seem to be completely adequate as used for predicting the pitch and heave motions of a non-linear type hull form such as Canada II. It appears that in order to overcome the effect of the non-linearities of the hull form, some sort of method is needed that calculates the added mass and damping and exciting force amplitude distributions for the immersed portion of the hull at any instant in time and then time steps through the oscillations to calculate the heave and pitch motions. A method of this sort was developed by Parissis [32] however is extremely time-consuming and not generally used. It also appears that some kind of velocity effect needs to be incorporated into the restoring coefficients. Although the strip theory results are not very encouraging quantitatively, they do provide us with a first approximation of the heave and pitch motion amplitude response spectrums. 5.2 Comparison of Effect of Anti-Pitch Foil on Motions In this section, the usefulness of the anti-pitch foil mathematical model in predicting the change in the heave and pitch motion due to the foil is discussed. Note that the emphasis here has shifted from predicting the amplitudes of the actual motions themselves, to predicting the change in motion due to the foil. As shown in section 3.3.3, the effect of the foil according to the theory, is to change the coefficients in the equation of motion. Thus, even though the strip theory does not give good agreement with experiment on motion amplitudes, it is hoped that the relative effect on the motions due to the foil can still be predicted reasonably using the anti-pitch foil theory. This would be a very useful preliminary design tool for a naval architect Chapter 5. Comparison of Theoretical and Experimental Results 78 considering the use of an anti-pitch foil. Figure 5.21 shows the experimental and theoretical results for heave of both the bare hull and the hull with foil. In looking at the effect of the foil, the mathematical model agrees with the experiment in predicting that the heave amplitude response for the hull with foil will be higher than that for the bare hull for p < 2.8. The crossover point e agrees quite well. Visually, it appears that above the crossover frequency, the magnitude of the decrease in heave amplitude between the bare hull and the hull with foil is similar for theory and experiment. Looking at this reduction due to the foil more closely, Figure 5.22 shows the percentage difference in heave amplitude response between the bare hull case and when the foil is attached for both experiment and theory. This figure shows that the theoretical prediction of the reduction in heave motion due to the foil agrees reasonably with experiment. Even though the strip theory does not give a good estimate of the heave motion itself, the effect of the foil on the heave amplitude response is predicted reasonably accurately. A comparison of the experimental and theoretical results for pitch motion of both the bare hull and the hull with foil attached is shown in Figure 5.23. Notice that the theory matches the experiment in predicting that the pitch amplitude decreases over the entire frequency range as a result of the addition of the foil. A similar plot to Figure 5.22 is shown for pitch in Figuire 5.24. Though not as good a match as the heave results, the theory again shows reasonable agreement with experiment concerning the amount of reduction in pitch amplitude that can be expected by adding the anti-pitch foil. In fact, near the resonant region there is actually more reduction in pitch shown by experiment, indicating that the theory gives a conservative estimate of the reduction due to the foil which could be useful as a preliminary design tool. The result of this analysis is that the mathematical model gives reasonable estimates of the effect of the foil in reducing the heave and pitch motions, in spite of the fact the Chapter 5. Comparison of Theoretical and Experimental Results 79 Experimental a n d Theoretical Heave Response 1.4 I Figure 5.21: Heave Amplitude Response With and Without Foil: Theory vs. Experiment Chapter 5. Comparison of Theoretical and Experimental Results Heave Response Difference Theory vs. Experiment 20-1 i Figure 5.22: Heave Amplitude Difference: Theory vs. Experiment 80 Chapter 5. Comparison of Theoretical and Experimental Results 81 Experimental a n d Theoretical Pitch R e s p o n s e Legend • E x p e r i m e n t - n o foil • Experiment—foil # Strip_Theory_^no_foij O Strip T h e o r y - f o i l 1 2 3 4 5 N o n - d i m e n s i o n a l encounter frequency Figure 5.23: Pitch Amplitude Response With and Without Foil: Theory vs. Experiment Chapter 5. Comparison of Theoretical and Experimental Results 82 Pitch Response Difference Theory vs. Experiment 20 10 T -ioH g p - • 2 I 0 -30 H -40 Legend • Strip Theory • Experimental -50 0 1 2 3 4 N o n - d i m e n s i o n a l encounter frequency Figure 5.24: Pitch Amplitude Difference: Theory vs. Experiment Chapter 5. Comparison of Theoretical and Experimental Results 83 strip theory predictions of actual motions are inadequate. 5.3 Added Resistance Comparison Figure 5.25 shows a comparison of the theoretical and experimental added resistance response for the bare hull. At p > 2.2, the theory matches the experiment fairly well, e though in the experimental results there is a small shift to the right in the response spectrum. This shift is probably due to the change in restoring coefficients due to forward speed, as discussed earlier, and the resultant effect of increasing the heave and pitch natural frequencies. This would cause a shift to the right in the experimental added resistance response as its natural frequency would also be increased. At resonance, the theory underpredicts the peak value of cr aw . One might have ex- pected that the experimental value of o- at resonance would be less than the theoretically aw calculated value since both heave and pitch amplitudes at resonance are much lower in the experiment than predicted by strip theory. However, as noted in section 4.3.1 the bow of the model was frequently observed to immerse under each wave at resonance causing a large amount of green water to be shipped on deck. Energy is expended in carrying and accelerating this water on deck to the velocity of the model and then dumping it again. This would show up as an increase in the resistance force on the model when the bow is submerged which explains the higher values of a from experiment. aw To examine how the added resistance theory predicts the effect of the foil on the added resistance refer to Figure 5.26 which shows the experimental and theoretical results for the bare hull and the hull with foil. From the figure, the difference between the bare hull and hull with foil cases appears to be similar in theory and experiment. Notice that the theory tends to be more accurate in predicting the height of the resonant peak in the foil case. This is probably because with the foil attached the pitch motions were reduced to er 5. Comparison of Theoretical and Experimental Results Added Resistance Response For Bare Hull Theory vs. Experiment 3.5- Legend • E x p e r i m e n t - n o foil • J o o s e n - n o foil 2.5 2- 1.5- 1- 0.5 •e- 1 2 3 4 N o n - d i m e n s i o n a l encounter frequency Figure 5.25: Added Resistance For Bare Hull: Theory vs. Experiment Chapter 5. Comparison of Theoretical and Experimental Results Experimental a n d Theoretical Added Resistance Response N o n - d i m e n s i o n a l encounter frequency Figure 5.26: Added Resistance With and Without Foil: Theory vs. Experiment 85 Chapter 5. Comparison of Theoretical and Experimental Results 86 the point that little bow submergence was observed. Figure 5.27 shows the percentage difference between the added resistance with the foil attached and that of the bare hull for theory and experiment. Notice that the added resistance is usually reduced by a larger amount in the experiments than the theory predicts. However the theory tends to give a reasonably conservative envelope of the amount of reduction in added resistance that can be expected from the foil. Again, this could be very useful to naval architects when trying to determine the effect of an anti-pitch foil on the performance of a ship in waves. It should be noted that the larger discrepency shown in the added resistance difference results as compared to the heave and pitch difference results may be due to variance in the waveheights that was observed during the experiments. This variance should not affect the motion responses as they have been shown to be proportional to waveheight for yacht hull forms as stated in section 2.4.2. However the variance in waveheight may affect the added resistance results adversely as added resistance has been shown to violate the waveheight squared law for some yacht hull forms. Chapter 5. 87 Comparison of Theoretical and Experimental Results Added Resistance Response Difference Theory vs. Experiment 1 20-1 10- l 0 1 1 1 1 [ 1 2 3 4 5 N o n - d i m e n s i o n a l encounter frequency Figure 5.27: Added Resistance Difference: Theory vs. Experiment Chapter 6 Conclusions and Recommendations 6.1 Conclusions This research has shown through a theoretical and experimental investigation that a passive bow anti-pitch foil is very effective in substantially reducing the heave and pitch motions and added resistance in head seas of the Canada II hull. Certain characteristics of the Canada II hull such as the flared topsides and extremely long overhangs, violate the linear assumptions made in ship motion theory. These nonlinearities lead to additional hydrodynamic damping that is not accounted for by the linear theory. As a result, ship motion theory overpredicts the response amplitudes of the bare hull for both heave and pitch, especially in the region of resonance. Experiments show an increase in the heave and pitch natural frequencies with increased speed indicating speed dependent restoring coefficients. Linear theory on the other hand, assumes speed independent restoring coefficients. It appears that the restoring forces are dependent on the the dynamic waterline resulting from the ship's own wave system. It is therefore concluded that linear strip theory is not completely adequate for accurately predicting the hydrodynamic coefficients and consequently the heave and pitch motions of the Canada II bare hull. However, strip theory does provide us with a good first approximation of these motion responses. The percentage reduction in heave and pitch motions due to the addition of the antipitch foil is predicted reasonably well by the oscillating hydrofoil theory. The 88 relative Chapter 6. Conclusions and Recommendations 89 effect of the foil is in good agreement with experiment in spite of the fact that the actual heave and pitch amplitudes, as predicted by strip theory, are not in adequate quantitative agreement. The Joosen theory gives a good prediction of the added resistance of the bare hull despite the limitations of strip theory. The percentage reduction in added resistance due to the addition of the foil is also predicted well by the Joosen theory. Experimentalfindingsfor added resistance may have been adversely affected by problems with the wave making apparatus. Difficulties in maintaining constant wave height and bad waveforms at certain frequencies may have resulted in experimental error. It is concluded that the mathematical model is very effective in predicting the percentage difference in motion amplitudes and added resistance due to the addition of a bow anti-pitch foil. For both heave and pitch motions and added resistance, the model tends to provide a conservative envelope of the expected reduction due to the foil. The mathematical model therefore provides a very useful preliminary design tool for naval architects, where no such tool existed before. 6.2 Recommendations This study was deliberately constrained in order to obtain useful conclusions. Therefore, recommendations for future study mostly involve the freeing of these constraints. The effect of heel angle on added resistance should be investigated as previous researchers have obtained conflicting results, especially in recent studies. Heel angle may also affect the efficiency of the anti-pitch foil with respect to reducing the heave and pitch motions. The motion and added resistance response spectra for regular waves should be applied Chapter 6. Conclusions and Recommendations 90 to actual coastal sea spectra to determine the predicted effectiveness of the foil in real random sea conditions. Irregular wave model tests would also be useful in this respect. The effect of yaw angle on the foil should be studied to determine if the strut and foil produce yaw induced moments that might adversely affect steering and directional stability of the vessel. It is recommended that work be done to optimize the strut and foil design in order to reduce the rather large increase in calm water resistance that was observed in this study. It is thought that one non-surface piercing strut rather than two surface piercing struts would go far in solving this problem. Any decrease in calm water resistance will improve overall performance which is a primary goal of the naval architect. In terms of future theoretical work, the use of the Gerritsma and Beukelman technique of calculating added resistance should be considered as it may provide more accurate results (at the expense of increased complexity). A theoretical investigation into the non-linear effect of long fore and aft overhangs on the hydrodynamic coefficients and the dynamic waterline effect on the restoring coefficients is recommended in order to improve the quantitative predictions of heave and pitch amplitudes. It would be useful to extend the theoretical results to higher frequencies. Use of the finite element method may overcome the problem of breaking down at high frequencies that is inherent in the formulation of boundary element methods. This study has provided the first step towards the possible use of bow anti-pitch foils to reduce motions and added resistance of small marine vehicles, in particular racing yachts. Of course there are several practical problems to consider before attempting a full scale trial of an anti-pitch foil. For instance, the effect of the weight of the system near the bow, the structural integrity of the hull and foil assembly and the effect of the foil on steering and directional stability, to name a few. In addition, the performance in Chapter 6. Conclusions and Recommendations 91 calm water can not be forgotten as the foil can only adversely affect this. Perhaps a form of retractable foil that is only deployed in rough water would be useful in this respect. In summary, it appears from the findings of this study that the concept of using a passive bow anti-pitch foil on racing yachts and other small marine vehicles to reduce vertical motions and added resistance in waves is well worth pursuing. Bibliography [1] M.A. Abkowitz. The effect of anti-pitching fins on ship motions. tions, SNAME Transac- pages 210-252, 1959. [2] Rameswar Bhattacharyya. Dynamics of Marine Vehicles. John Wiley and Sons, 1978. [3] J.E. Conolly and G.J. Goodrich. Sea trials of anti-pitching fins. tions, R.I.N.A. Transac- pages 87-100, 1969. [4] G.L. Cox and R. G. Whitaker. Some miscellaneous aspects of modern 12-metre design. Conference on Yachting Technology 1987, pages 5-11, 1987. [5] Frank DeBord Jr., A.C. Dovell, and Roderick Y. Edwards Jr. Extrapolation of sailing yacht model data to full scale. Technical Report, Offshore Technology Corporation, 1985. [6] A.C. Dovell. Model tests of 12-metre Canada I with keels a,b and c. Technical Report OTC-85-252, Offshore Technology Corporation, January 1986. [7] W. Frank and N. Salvesen. The Frank close-fit computer program. Technical Report 3289, NSRDC, 1970. [8] lr. Gerritsma. Course keeping qualities and motions in waves of a sailing yacht. In Procedeeings of the Third AIAA Symposium on the Aero/Hydronautics pages 10-29. American Institute of Aeronautics and Astronautics, 1971. 92 of Sailing, 93 Bibliography [9] Ir. Gerritsma and W. Beukelman. Analysis of the modified strip theory for the calculation of ship motions and wave bending moments. Progress, International Shipbuilding 14(156):319-337, 1967. [10] Ir. Gerritsma and W. Beukelman. Analysis of the resistance increase in waves of a fast cargo ship. International Shipbuilding Progress, 19(217):285-293, 1972. [11] Ir. J. Gerritsma and J. A. Keuning. Performance of light and heavy-displacement sailing yachts in waves. Marine Technology, 26(l):14-22, 1989. [12] Ir. J. Gerritsma and G. Moeyes. The seakeeping performance and steering properties of sailing yachts. Symposium of Yacht Architecture, pages 107-156. HISWA, 1973. [13] Ir. J. Gerritsma, Ir. J.J. Van Den Bosh, and W. Beukelman. Propulsion in regular and irregular waves. International Shipbuilding Progress, 8(82), 1962. [14] A. Giddings and R. Wermter. A survey of ship motion stabilization. on Naval Hydrodynamics, Fifth Symposium 1964. [15] T. Hanaoka and et al. Researches on seakeeping qualities of ships in Japan. The , Society of Naval Architects of Japan, 60th Anniversary [16] T.H. Havelock. Drifting force on a ship among waves. Series, 8, 1963. Philosophical Magazine, 33, 1942. [17] W.P.A. Joosen. Added resistance of ships in waves. In Proceedings posium on Naval Hydrodynamics, [18] Paul Kaplan. Motion. A Hydrodynamic of the 6th Sym- 1966. Theory for the Forces on Hydrofoils PhD thesis, Stevens Institute of Technology, 1955. in Unsteady 94 Bibliography [19] K. Klaka and J.D. Penrose. Performance prediction of sailing yachts in waves. ference on Yachting Technology 1987, Con- pages 34-39, 1987. [20] B.V. Korvin-Kroukovsky and Winnefred R. Jacobs. Pitching and heaving motions of a ship in regular waves. SNAME Transactions, pages 590-632, 1957. [21] Theodore Loukakis and Paul Sclavounos. Some extentions of the classical approach to strip theory of ship motions, including the calculation of mean added forces and moments. Journal of Ship Research, 22(1):1—19, 1978. [22] Y.L. Luke and M.A. Dengler. Tables of the Theodorsen circulation function for generalized motion. [23] H. Journal of the Aeronautical Sciences, Maruo. The drift of a body floating on waves. 1951. Journal of Ship Research, 4(3), 1960. [24] M.M. Munk. Aerodynamic Theory. Durand Reprinting Committee, California Insti- tute of Technology, 1943. [25] David R. Pedrick. The performance of sailing yachts in oblique seas. nology, Marine Tech- pages 383-392, October 1974. [26] Nils Salvesen, E.O. Tuck, and Odd Faltinsen. Ship motions and sea loads. SNAME Transactions, pages 250-287, 1970. [27] C. Scragg, B. Chance, J. Talcott, and D. Wyatt. The analysis of wave resistance in the design of the twelve meter yacht Stars and Yacht Symposium, pages Stripes. In Eighth Chesapeake Sailing 109-121. SNAME, 1987. [28] O.J. Sibul. Constant thrust versus constant velocity method for resistance measurements in waves. Technical Report NA-71-1, College of Engineering, University of 95 Bibliography California, 1971. [29] O.J. Sibul. Measurememts and calculations of ship resistance in waves. Technical Report NA-71-2, College of Engineering, University of California, 1971. [30] W. E. Smith. Computation of pitch and heave motions for arbitrary ship forms. International Shipbuilding Progress, 14:267-291, 1967. [31] P. G. Spens. Research on the reduction of pitching motions of ships by controllable fins. Technical Report 913, David Taylor Model Basin, 1962. [32] P.G. Spens, Pierre DeSaix, and P.W. Brown. Some further experimental studies of the sailing yacht. SNAME Transactions, pages 79-111, 1967. [33] G.P. Stefun. Model experiments with fixed bow anti-pitching fins. Research, 3(2), Journal of Ship 1959. [34] G.N. Stensgaard and G.A. Roddan. Model tests of Canada II. Technical Report 6-02-588-1, B.C. Research Ocean Engineering Centre, December 1986. [35] Jorgen Strom-Tejsen, Hugh Yeh, and David Moran. Added resistance in waves. SNAME Transactions, pages 109-143, 1973. [36] Matao Takagi, Yukiharu Nekado, Masahiro Matsui, and Masaru Mokunaka. Activated anti-pitching fins. Japan Shipbuilding and Marine Engineering, pages 31-40, 1967. [37] R.F. Van Hemmen. 12-metre design: State of the art in 1986. Marine Technology, 23(4):320-337, 1986. [38] G Vossers. Fundementals of the behaviour of ships in waves. Technical Report 115a, Netherlands Ship Model Basin, 1961. Bibliography 96 [39] P.A. Wilson. A review of the methods of calculation of added resistance for ships in a seaway. Journal of Wind Engineering, 20:187-199, 1985. Appendix A Sample Designer Questionnaire What do you f e e l are the most Important boat parameters and f a c t o r s c o n t r i b u t e to the added r e s i s t a n c e of a yacht due t o waves? that In approximate order o f importance; 1) 2) 3) 4) 5) S a i l i n g length Volume d i s t r i b u t i o n Weight d i s t r i b u t i o n Freeboard Hull flare T h e s e a r e a l l areas where c o r r e c t c h a r a c t e r i s t i c s w i l l reduce added r e s i s t a n c e due to waves; most of them are s e l f - e v i d e n t , f o r i n s t a n c e g r e a t e r length w i l l g e n e r a l l y r e s u l t i n l e s s motion, hence l e s s added r e s i s t a n c e . Do you t h i n k a r e d u c t i o n i n added r e s i s t a n c e due t o waves would r e s u l t i n a s i g n i f i c a n t i n c r e a s e i n upwind performance or i s i t t o o small a f a c t o r ? In a s i g n i f i c a n t seaway, a d d e d r e s i s t a n c e due t o waves i s c e r t a i n l y important. O b v i o u s l y , t h e c o n v e r s e i s a l s o t r u e ; i n no w a v e s , t h e r e i s no a d d e d r e s i s t a n c e . 97 Appendix A. Sample Designer Questionnaire What would you estimate the added r e s i s t a n c e due t o waves t o be o f t h e t o t a l r e s i s t a n c e of a yacht t r a v e l l i n g upwind a t 8 knots i n t o a Fremantle Seaway? (%) Roughly 201. Did you have a segment o f your d e s i g n program aimed r e d u c i n g added r e s i s t a n c e o f the yacht due t o waves? No. Money a n d t i m e constraints prevented it. s p e c i f i c a l l y at 99 Appendix A. Sample Designer Questionnaire How important a f a c t o r was added r e s i s t a n c e due to waves i n the development of your 12-metre design? What aspects of your design s p e c i f i c a l l y addressed this problem? Not t o o i m p o r t a n t . U n a b l e t o t e s t , we r e l i e d on f l a t w a t e r p e r f o r m a n c e and t r i e d t o e n s u r e s u f f i c i e n t f l a r e and f r e e b o a r d to keep most o f t h e waves from coming a b o a r d . Did you perform towing tank waves program? tests during your 12-metre design No. If so, to what extent did you actually use the results of these tests in your design process? N/A 100 Appendix A. Sample Designer Questionnaire Do you feel the relative performance of different models in towing tank tests in calm water can be extrapolated to rough water with any confidence? No. But i n a h i g h l y r e s t r i c t e d c l a s s l i k e Twelve Meters, where weight d i s t r i b u t i o n i s p r e t t y w e l l c o n t r o l l e d and length tends to get stereotyped, there are not too many c o n t r o l l a b l e f a c t o r s l e f t to a l t e r r e l a t i v e performance dh waves. Offshore designs are probably a more f e r t i l e ( i f l e s s w e l l financed) f i e l d f o r advances i n the r e d u c t i o n of added r e s i s t a n c e due to waves. Did you carry out computer ship motion prediction simulations as part of your design program? No. Time and money a g a i n . I f s o , what inconclusive? N/A programs did you use and were the results useful or Appendix A. 101 Sample Designer Questionnaire Did you use o r contemplate using d u r i n g s a i l i n g i n Fremantle? any sort of p i t c h r e d u c i n g technique Yes. Buddy was c a r e f u l to s t e e r around the waves ( l a t e r i n the r e g a t t a , when both the wave h e i g h t and the crew's a b i l i t y had i n creased s i g n i f i c a n t l y over the e a T l y r o u n d s ) , Do you have any papers on the s u b j e c t o f yacht performance i n waves? I f so, p l e a s e l i s t the t i t l e s and authors and i n what p u b l i c a t i o n I would find i t . The Performance of S a i l i n g Yachts i n Oblique Seas by David P e d r i c k Marine Technology, October 1974 S u r f i n g : Motions o f a V e s s e l Running Letcher, J r . i n Large Waves by John s. The T h i r d Chesapeake S a i l i n g Yacht Symposium (1977) Appendix B Model Foil Drawings 102 fAooEL F O I L f ^ S E ^ f c i - W Appendix C Experimental Data Results of experiments on Canada 2 with and without foil RUN NUMBER HAVE FREQ HZ HAVE FREfi RAD/SEC HAVE HEIGHT INS HAVE LENGTH FT CALC. MEASURED MODEL MODEL ENCOUNTER ENCOUNTER VELOCITY VELOCITY FREQUENCY FREQUENCY H/S FT/S RAD/S HI RAD/S SHIP VELOCITY WTS 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 0.349 0.41 0.453 0.505 0.552 0.597 0.638 0.682 0.722 0.758 0.799 0.837 0.871 0.907 0.939 2.19 2.58 2.85 3.17 3.47 3.75 4.01 4.29 4.54 4.76 5.02 5.26 5.47 5.70 5.90 3.767 3.742 3.755 3.795 3.839 3.862 3.746 3.777 3.747 3.886 3.834 4.045 4.212 3.747 4.011 42.08 30.49 24.97 20.10 16.82 14.38 12.59 11.02 9.83 8.92 8.03 7.32 6.76 6.23 5.81 1.46 1.46 1.46 1.46 1.46 1.46 1.46 1.46 1.46 1.46 1.46 1.46 1.46 1.46 1.46 4.79 4.79 4.79 4.79 4.79 4.79 4.79 4.79 4.79 4.79 4.79 4.79 4.79 4.79 4.79 0.463 0.567 0.645 0.743 0.837 0.930 1.018 1.117 1.209 1.295 1.396 1.492 1.580 1.676 1.763 0.475 0.574 0.659 0.751 0.843 0.928 1.023 1.119 1.21 1.36 1.393 1.484 1.567 1.661 1.754 2.98 3.61 4.14 4.72 5.30 5.83 6.43 7.03 7.60 8.55 8.75 9.32 9.85 10.44 11.02 8.02 8.02 8.02 8.02 8.02 8.02 8.02 8.02 8.02 8.02 8.02 8.02 8.02 8.02 B.02 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 0.35 0.405 0.455 0.503 0.551 0.6 0.642 0.681 0.721 0.759 0.798 0.836 0.874 0.909 0.942 2.20 2.54 2.86 3.16 3.46 3.77 4.03 4.28 4.53 4.77 5.01 5.25 5.49 5.71 5.92 3.8 3.8 3.803 3.873 3.9 4.086 3.938 3.906 3.851 4.027 3.883 4.103 4.293 3.707 3.897 41.84 31.24 24.75 20.26 16.88 14.24 12.43 11.05 9.86 8.90 8.05 7.33 6.71 (.20 5.78 1.44 1.47 1.41 1.49 1.46 1.43 1.44 1.46 1.46 1.44 1.44 1.47 1.46 1.45 1.46 4.72 4.82 4.62 4.89 4.79 4.69 4.72 4.79 4.79 4.72 4.72 4.82 4.79 4.76 4.79 0.463 0.559 0.642 0.744 0.835 0.929 1.022 1.114 1.207 1.290 1.385 1.494 1.588 1.676 1.771 0.463 0.563 0.654 0.745 0.838 0.927 1.016 1.121 1.195 1.295 1.381 1.474 1.574 1.654 1.742 2.91 3.54 4.11 4.68 5.27 5.82 6.38 7.04 7.51 B.14 8.68 9.26 9.89 10.39 10.95 7.91 8.07 7.74 8.18 8.02 7.85 7.91 8.02 8.02 7.91 7.91 8.07 8.02 7.96 8.02 104 Appendix C. Experimental Data 105 PITC MO* MEASURED BIN M K P - f ) PHKP-P) PITCHAHP/ EMC. FREQ. EMC. FREO. MEASURED MA VE SLOPE HZ BAD/SEC K6REES RADIANS «*VE MEASURED EMC. rtEQ. HZ RAD/SEC MM M l MEAVE(P-f) HEAVE/ EMC. FRED. MAVEHEI6HT IMS. 0.4(3 0.571 0.66 0.75 0.84 0.93 1.023 1.119 1.21 1.303 1.395 1.485 1.564 1.664 1.753 2.909 3.588 4.147 4.712 5.278 5.843 6.428 7.031 7.603 8.187 8.765 9.331 9.827 10.455 11.014 1.245 1.536 1.775 2.017 2.259 2.501 2.751 3.010 3.254 3.504 3.752 3.994 4.206 4.475 4.715 2.973 4.473 5.072 6.219 7.5BB 8.886 9.931 11.069 10.84 10.736 8.472 6.967 6.177 3.84 2.587 0.052 0.078 0.089 0.109 0.132 0.155 0.173 0.193 0.189 0.187 0.148 0.122 0.108 0.067 0.045 1.107 1.215 1.124 1.098 1.108 1.103 1.113 1.076 0.948 0.821 0.591 0.420 0.330 0.213 0.125 0.481 0.571 0.661 0.754 0.848 0.94 1.039 1.121 1.209 1.304 1.395 I.48S 1.56 1.667 1.756 3.022 3*588 4.153 4.738 5.328 5.906 6.528 7.043 7.596 8.193 8.765 9.331 9.802 10.474 11.033 1.294 1.536 1.778 2.028 2.281 2.528 2.794 3.015 3.252 3.507 3.752 3.994 4.196 4.483 4.723 3.696 3.842 3.9 3.977 4.034 3.945 3.9 4.218 4.046 4.017 3.059 2.284 1.85 0.964 0.569 0.981 1.027 1.039 1.048 1.051 1.021 1.041 1.117 1.080 1.034 0.798 0.565 0.439 0.257 0.142 0.463 0.565 0.647 0.752 0.838 0.923 1.015 1.116 1.192 1.288 1.381 1.473 1.567 1.656 1.739 2.909 3.550 4.065 4.725 5.265 5.799 6.377 7.012 7.490 1.093 8.677 9.255 9.846 10.405 10.926 1.245 1.520 1.740 2.023 2.254 2.482 2.730 3.002 3.206 3.464 3.714 3.962 4.215 4.454 4.677 3.011 4.157 4.891 5.882 7.039 7.998 8.516 1.787 1.168 7.764 (.114 5.108 4.356 2.839 2.07 0.053 0.073 0.085 0.103 0.123 0.140 0.149 0.153 0.143 0.136 0.107 0.089 0.076 0.050 0.036 1.105 1.139 1.061 1.025 1.016 0.929 0.896 0.829 0.697 0.572 0.422 0.304 0.227 0.158 0.102 0.483 0.568 0.648 0.753 0.838 0.924 1.015 1.118 1.195 1.288 1.379 1.474 1.561 1.651 1.728 3.035 3.569 4.072 4.731 5.265 5.806 6.377 7.025 7.508 8.093 8.665 9.261 9.808 10.374 10.857 1.299 1.528 1.743 2.025 2.254 2.485 2.730 3.007 3.214 3.464 3.709 3.964 4.198 4.440 4.648 3.81 4.096 4.201 4.283 4.395 4.329 4.152 4.001 3.515 3.119 2.144 1.608 1.265 0.675 0.425 1.003 1.078 1.105 1.106 1.127 1.059 1.054 1.024 0.913 0.775 0.552 0.392 0.295 0.1B2 0.109 Appendix C. Experimental Data 106 CALH HON. DIM. HATER TOTAL ADDED ENC. FREQ. RESIST. RESIST. RESIST. LBS LBS LBS S16AA 1.28 1.43 1.641 1.54 1.77 1.43 1.43 1.646 2.02 2.27 1.43 1.43 1.783 1.777 2.50 2.75 3.01 3.25 3.66 3.75 3.99 4.21 4.47 4.72 1.43 1.43 1.43 1.43 1.43 1.43 1.43 1.43 1.43 1.43 1.865 2.079 0.435 0.649 2.433 2.789 3.529 3.256 3.177 3.449 2.715 1.003 1.359 2.099 1.826 1.747 2.019 1.285 3.351 2.881 3.070 2.469 2.541 1.111 1.863 1.25 1.51 1.74 1.84 1.964 2.014 0.224 0.174 1.76 1.63 1.855 0.225 0.419 0.325 0.420 2.00 2.25 2.49 2.73 3.01 3.21 3.48 3.71 3.96 4.23 4.45 4.69 1.9 1.81 1.7 1.74 1.81 1.81 1.74 1.74 1.84 1.81 1.77 1.81 2.224 2.102 2.056 2.25 2.662 2.66 3.104 1.675 0.211 0.216 0.245 0.353 0.347 0.401 0.416 0.469 0.661 0.635 0.787 1.248 1.897 2.612 3.750 0.51 0.852 0.85 1.364 0.583 0.518 0.575 0.887 1.507 1.546 2.269 2.937 2.931 3.384 2.838 1.197 1.091 1.574 1.068 2.142 1.748 2.304 2.097 2.683 0.873 1.551 0.324 0.292 0.356
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Reducing added resistance using an anti-pitch foil
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Reducing added resistance using an anti-pitch foil Avis, John Spencer 1989
pdf
Page Metadata
Item Metadata
Title | Reducing added resistance using an anti-pitch foil |
Creator |
Avis, John Spencer |
Publisher | University of British Columbia |
Date Issued | 1989 |
Description | In this study, the use of a passive bow anti-pitch foil for the purpose of reducing the vertical motions and added resistance in waves of a small marine vehicle is investigated. The difficulties, limitations and assumptions inherent in studying ships in waves are discussed. A mathematical model is developed incorporating linear strip theory, Joosen added resistance theory and oscillating hydrofoil theory in an attempt to predict the effect of a hypothetical anti-pitch foil on the heave and pitch motions and added resistance in head seas of the 12-metre yacht Canada II. The results of regular head wave towing tank experiments on a 1:8 scale model of Canada II are presented for both the bare hull and for the hull fitted with a model anti-pitch foil to demonstrate the experimental effect of the foil on vertical motions and added resistance. The results of the theoretical and experimental investigations show that the anti-pitch foil is in fact very effective in reducing the motions and added resistance of the Canada II hull. Experimental results show heave amplitude is reduced by as much as 15%, pitch amplitude by 22% and added resistance by as much as 40%. The heave and pitch motion amplitude responses predicted by the linear strip theory for the hull with no foil, do not agree well with experiment, especially in the region around resonance. It is concluded that the linear strip theory used is not completely adequate in predicting the hydrodynamic coefficients and consequently the heave and pitch motions of Canada II due to non-linearities in her hull form. These non-linearities, in particular the flared topsides and large fore and aft overhangs, violate the linear assumptions of strip theory and result in increased hydrodynamic damping and speed dependent restoring coefficients that are not accounted for by theory. The Joosen added resistance theory predictions for the bare hull do agree well with experiment despite the hull non-linearities. The change or percentage reduction in heave and pitch amplitude response due to the addition of the anti-pitch foil is predicted well by the oscillating hydrofoil theory, in spite of the inaccuracies of the strip theory in predicting the actual motions. Also, the Joosen theory predictions of the percentage reduction in added resistance due to the addition of the foil agree very well with experiment. It is concluded that the mathematical model provides a very good conservative estimate of the amount of reduction in heave and pitch motion and added resistance that can be expected from the addition of a passive bow anti-pitch foil. The model could be very useful as a preliminary design tool for naval architects, where no such tool existed before. The findings of this study indicate that the concept of using a passive bow anti-pitch foil to reduce the vertical motions and added resistance in waves is well worth pursuing. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2010-08-25 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0080679 |
URI | http://hdl.handle.net/2429/27758 |
Degree |
Master of Applied Science - MASc |
Program |
Mechanical Engineering |
Affiliation |
Applied Science, Faculty of Mechanical Engineering, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1989_A7 A94.pdf [ 14.74MB ]
- Metadata
- JSON: 831-1.0080679.json
- JSON-LD: 831-1.0080679-ld.json
- RDF/XML (Pretty): 831-1.0080679-rdf.xml
- RDF/JSON: 831-1.0080679-rdf.json
- Turtle: 831-1.0080679-turtle.txt
- N-Triples: 831-1.0080679-rdf-ntriples.txt
- Original Record: 831-1.0080679-source.json
- Full Text
- 831-1.0080679-fulltext.txt
- Citation
- 831-1.0080679.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0080679/manifest