Froude Krylov Excitation Force

Transcription

1 13.4 Design Principles for Ocean Vehicles Prof. A.H. echet Spring 005 Froude Krylov Ecitation Force 1. Radiation and Diffraction Potentials he total potential is a linear superposition of the incident, diffraction, and radiation potentials, i t φ = ( φ + φ + φ ). (1) I D R e ω he radiation potential is comprised of si components due to the motions in the si directions, φ where = 13456,,,,,. Each function φ is the potential resulting from a unit motion in th direction for a ody floating in a quiescent fluid. he resulting ody oundary condition follows from lecture 15: φ = iωn ; ( = 1,, 3) φ = iω( r n ) 3; ( = 4, 5, 6) r = (, y, ) () (3) (4) n = n ( = 1,, 3) = ( n, n, n ) (5) y In order to meet all the oundary conditions we must have waves that radiate away from the ody. hus φ e i as ±. For the diffraction prolem we now that the derivative of the total potential (here the incident potential plus the diffraction potential without consideration of the radiation potential) normal to φ the ody surface is ero on the ody: n = 0 on S B, where φ = φi + φd. version 1.0 updated 3/9/ , aht

2 φi φ D = ; on S B (6) We have so far taled primarily aout the incident potential. he formulation of the incident potential is straight forward from the oundary value prolem (BVP) setup in lecture 15. here eist several viale forms of this potential function each are essentially a phase shifted version of another. he diffraction potential can also e found in the same fashion using the BVP for the diffraction potential with the appropriate oundary condition on the ody. his potential can e approimated for a long wave condition. his long wave approimation assumes that the incident wavelength is very long compared to the ody diameter and thus the induced velocity field from the incident waves on the structure can e assumed constant over the ody and approimated y the following equation: i φ φ φ φd φ φ φ (7) ω y I I I Further eplanation of this approimation can e found in Newman (p. 301). Ultimately, if we assume the ody to e sufficiently small as not to affect the pressure field due to an incident wave, then we can diffraction effects can e completely ignored. his assumption comes from the Froude-Krylov hypothesis and assures a resulting ecitation force equivalent to the froude-rylov force: φi F () t = ρ nds (8) t version 1.0 updated 3/9/ , aht

3 . Vertical Froude-Krylov Force on a Single Hull Vessel B Deep water incident wave potential is: aω i( t ) I e Re ie ω φ = (9) he force in the vertical direction is found from the incident potential using eq. 8 along the ottom of the vessel. Here the normal in the -direction, n, is negative: n = 1, so the force per unit length in the -direction is B / iaω i( ωt ) F = Re ρ iω e e d (10) B/ ρoω = Re a e e e e (11) iωt ib/ ib/ ω iωt = Re ρ a e e sin( B/ ) (1) i i e e Recall that sin =. i Using the vertical velocity we can rewrite the force in terms of the velocity. wt () = Re aω e ie ω (13) i( t ) wt () = Re aω e e ω (14) i( t ) (15) i t w ( 0 0 t) Reaω e ω =, =, = version 1.0 updated 3/9/ , aht

4 Now we can write the force in the vertical direction as a function of the vertical (heave) acceleration, ρ F Re e = sin( B/ ) w (0, 0, t). (16) Let s loo at the case where ω 0 the wavenumer, following simplifications can e made: 0 = ω / g, also goes to ero and the t e 1 (17) sin( B/ ) B/ (18) to yield a simplified heave force. ω iωt F Re ρ a (1 )( B/ ) e (19) ω iωt Re ρ g ab 1 e (0) g If we loo at the case where ω 0 and consider the heave restoring coefficient, C33 = ρ g B, i( t ) and the free surface elevation, η( t, ) = Re ae ω we can rewrite this force as { η =, } F Re C33 ( 0 t) (1) 3. Horiontal Froude-Krylov Force on a Single Hull Vessel he horiontal force on the vessel aove can e found in a similar fashion to the vertical force. F = n S ds () B t I ρ φ i a 0 ω = Re ρ i ω e d e e (3) i( ωt B/ ) i( ωt+ B/ ) aω iωt = Re iρ 1 e e sin( B/ ) (4) version 1.0 updated 3/9/ , aht

5 As frequency approaches ero similar simplifications can e made lie aove for the vertical force: aω iωt F () t Re iρ ( K) e B/ (5) ut () = Rea ω e e ω (6) i( t ) ut () = Reia ω e e ω (7) i( t ) { ρ =, =, } F () ( 0 0 ) t Re Bu t (8) Where ρb = ρ, and is the vessel volume such that we are left with the surge force F ρ u (9) F C η + ρ w (30) Multi Hulled Vessel -B/ B/ Again, let s mae a few asic assumptions: ( / λ << 1), ( B/ λ 1), ( a< ), and ( ). Let s loo at the force in the -direction: F ρ u ( = B/, = 0, t) + ρ u ( = B/, = 0, t) (31) η(, t) = acos( ωt ) (3) ut (,, ) = aω e sin( ωt ) (33) (34) version 1.0 updated 3/9/ , aht

6 { } F ρ aω ωt+ B/ + ωt B/ ( ) sin( ) sin( ) (35) ρ ( aω )cos( B/ )sin( ωt) (36) Note that when B/ = π/ (or B = λ/ ) then F ( t ) = Multi Hulled Vessel with additional pontoon -B/ B/ c Use the same assumptions from aove to find the -force adusted for the additional pontoon etween the two hulls. F a B/ t + cp( = B/ + /, = 0, t) cp( = B/ /, = 0, t) ρ ( ω )cos( )sin( ω ) he last two terms are the adustment to the force for the addition of the pontoon, δ F () t. Pressure is found from the incident potential: p( t,, ) = ρgae cos( ωt ). δf = ρ g asin( ωt)sin ( B ) (37) For B >> using g = ω / we get a force: { } F ( t) ρ a ω sin( ωt) cos( B/ ) + δ/ sin( B/ ) (38) version 1.0 updated 3/9/ , aht