CC BY
6
-□ □-
Надано ршення задачi про дифрак-цю хвилювання малог амплтуди, яке набиае тд довшьним кутом на нерухоме судно в умовах мшководдя. Потенщал швидкостей дифрагованого хвильового руху визначений методом зрощуваних асимптотичних розкладань. Виконаш розрахунки амплтуд хвилювання у зав-даних точках навколо судна, наведем приклади хвильових полiв
Ключовi слова: дифракщя лтшних хвиль на судт, метод зрощуваних
асимптотичних розвинень, мшководдя □-□
Представлено решение задачи о дифракции волнения малой амплитуды, набегающего под произвольным углом на неподвижное судно в условиях мелководья. Потенциал скоростей дифрагированного волнового движения определен методом сращиваемых асимптотических разложений. Выполнены расчеты амплитуд волнения в заданных точках вокруг судна, приведены примеры волновых полей
Ключевые слова: дифракция линейных волн на судне, метод сращиваемых асимптотических разложений, мелководье
-□ □-
UDC 629.12:532.059.041
|DOI: 10.15587/1729-4061.2017.993121
DETERMINING THE CHARACTERISTICS OF DIFFRACTED WAVES OF SMALL AMPLITUDE AROUND A VESSEL IN SHALLOW WATER
N. Ef remova
PhD, Associate Professor* E-mail nina290264@gmail.com A. Nilva PhD
Department of seaworthiness Marine Engineering Bureau Tenista str., 15, Odessa, Ukraine, 65009 E-mail: alexnilva@gmail.com N. Kotovskaya Senior lecture* E-mail: kotov-nata@rambler.ru M . D r y h a Senior lecture* E-mail: ridermarina@rambler.ru *J. L. Vorobyov Department of the Theory and Design of the Ship
Odessa National Maritime University Mechnikova str., 34, Odessa, Ukraine, 65029
1. Introduction
In most of the ports of the Black and Azov Seas, large-capacity tankers and bulk carriers, as well as some passenger vessels, have to be additionally loaded, unloaded or fuelled outside of the ports protected by water basins.
The water areas of ports in Ukraine are shallow-water. The ships that are docked on outer anchorage are under the influence of the rough sea and, at the same time, they prevent the propagation of waves. A wave field around the ship at the port roads is essentially three-dimensional; it represents a superposition of incoming and diffracted waves. The height of waves in specific points of the field depends on the location of points around the ship.
The complexity of evaluating the interaction between a ship and such waves is associated with the non-stationarity of the region, which is filled with fluid, and the non-linearity of boundary condition at the free surface of the fluid. The solution of the diffraction problem in a general form is complicated. Solving the particular problems requires refinement and specialization of boundary conditions.
With respect to ships, a diffraction problem is tackled rarely, since overcoming the difficulties associated with the practical evaluation of characteristics of a wave field around the body of a vessel requires special processing of the com-
puting aspects of the obtained theoretical results. At the same time, for determining the hydrodynamic forces acting on a ship in pitching, as shown in articles [1, 2], a solution of a simpler task of radiation will suffice.
Characteristics of the waves around the vessel, which is anchored at the open port roads, are to be considered when planning the operation of auxiliary ships. Auxiliary vessels are tugs, bunkering vessels, pilot and harbor boats, as well as oil/garbage collector boats and boom crafts. The work of these vessels is related to safe navigation (transfer of people and cargo from one ship to another) and environmental safety (elimination of oil and petroleum spills).
The height of waves imposes restrictions on the operation of auxiliary ships. Determining the transformation of waves on outer anchorage is necessary for their safe operation, which is why this is a relevant task.
2. Literature review and problem statement
In recent years, most of the existing solutions of diffraction problem have focused on assessing the impact of waves on the stationary or movable offshore structures.
Article [3] examined diffraction and refraction of waves in a fluid with defined depth. The problem is solved by
©
the finite element method using discrete non-local (DNL) boundary conditions. The examined objects are a channel of rectangular cross section, circular cylindrical island with a parabolic bottom around it. The characteristics of wave fields are presented.
A diffraction problem is important, in particular, for the equipment that employs the energy of ocean waves to generate electricity - Oscillating Water Column (OWC). Paper [4] presented a first order analytical solution of the problem on diffraction of ocean waves on the hollow-body vertical cylinder in the ocean with defined depth. In article [5], authors defined a wave field for the same object, created by a swinging cylinder, and solved a combined problem on diffraction-radiation.
Article [6] investigated numerical simulation (CFD) of the interaction between waves and a vertical cylinder. The Reynolds equation of the averaged turbulent fluid flow is solved (Reynolds Averaged Navier-Stokes, RANS). They modeled regular and irregular waves of small and finite amplitude (second-order Stokes) in the numerical testing pool. The results of determining the wave fields and wave forces are presented.
Paper [7] is devoted to the experimental study of influence of waves on a floating cylinder. In a small-sized wave flume, ocean waves are simulated. They examine their interaction with the floating cylinder, docked to the shore. The characteristics of waves and motion of a floating body are presented.
Article [8] describes numerical modeling of the interaction between nonlinear waves and a system of two vertical circular cylinders. One of them is rigidly fixed to the bottom, the second floats. Authors identified characteristics of wave forces and moments, as well as displacement of the floating cylinder. In connection with the operation of tension leg platforms (TLPs), the problems are solved on the impact of waves on the groups of vertical cylinders.
Paper [9] described the solution of potential problems of first and second orders by the method of finite elements. Vertical motions of vertical circular cylinders are examined. Wave fields and wave loads are determined for a single cylinder and the groups of two and four cylinders.
Article [10] demonstrates interaction between nonlinear waves and a vertical cylinder and the group of four cylinders. Authors employed the method of finite differences (FDM), methods of finite elements (Weakly Nonlinear and Weakly Dispersive FEM, Fully Nonlinear and Weakly Dispersive FEM). They defined the elevation of a free surface and coefficients of hydrodynamic forces at propagation and collision of separate waves. Results of the calculations are presented, as well as a comparison to experimental data.
Worth noting here is a complete analytic solution, presented in paper [11], of second order nonlinear diffraction problem for the two-dimensional stationary rectangular cylinders at the free surface of the fluid with defined depth. Authors determined the magnitudes of vertical and horizontal forces of first and second orders. The likelihood of solution is confirmed by a comparison with experiments and calculations of other authors.
Article [12] presented a solution of the problem on diffraction of monochromatic and bichromatic waves on a stationary horizontal cylindrical body that crosses free surface. The depth of the fluid is infinite, the incident waves are lateral or scant. Diffraction potential is used to determine the forces acting on the floating oil tank (a body with semi-el-
liptical waterlines in the bow part, rectangular frames in the middle part and prismatic stern).
An analysis of the above research results reveals that wave fields are defined around the objects of simplified forms.
Paper [13] considers a two-dimensional nonlinear potential problem on the fluctuations of cramped contour in the fluid of limited depth. Boundary conditions on the contour and the free surface of fluid are non-linear. Nonlinear forces are determined with an accuracy to the second order. Calculations are performed for different cramped contours. The influence of change in the relative depth on the value of nonlinear forces is investigated.
Article [14] addresses a three-dimensional potential problem on the fluctuations of a vessel in the liquid wtih limited depth and its solution by a numerical method. The influence of change in the depth on the values of attached masses and damping coefficients is examined; results of calculations of the given magnitudes for different types of vessels are presented.
Paper [15] defined the potential of radiation and diffraction potential in the pitching of a vessel moving in the significantly shallow water. The method used is the matched asymptotic expansion method (MAEM). The potentials of radiation for longitudinally-horizontal, vertical, and pitching fluctuations were determined at motion in the quiet water. The expression for the components of diffraction potential at vertical pitching is presented.
Article [16], by employing the improved matching method, defined the wave profiles around the ship that moves in deep water. The region, which is filled with fluid, is divided into a near-field and a far-field, which applies a radiation condition while the boundary condition at free surface is considered to be linear. The vessel is replaced with a system of features. Characteristics of these features and the potential of speeds of the near field are determined by a double technique. The first is the use of conditions on free surface and on the body in the near-field. The second is using the continuity of speed potentials and its derivative by a normal in the transition through the interface (matching surface).
For ships, a diffraction problem is solved when determining the wave loads in deep water or for determining the hy-drodynamic forces in pitching in the shallows; wave profiles around the ship are not investigated.
The waves around the ship significantly differ from the waves at a significant distance from the ship. A wave field is transformed around the large tonnage vessels at the outer anchorage. This must be considered when assessing the safety of operation of auxillary (relatively small) vessels.
Thus, the study of wave diffraction processes at diffraction of regular waves of small amplitude that is incident under an arbitrary angle on the hull of a large tonnage ship, which floats without a run on the outer anchorage, is a relevant issue. Considering that under conditions of limited depth of the anchorage the occurrence of long waves that can cause noticeable pitching of a large tonnage ship is highly unlikely, the latter is considered stationary.
A choice of linear wave theory for solving the diffraction problem is predetermined by the following considerations.
First, the waves that come from deepwater areas of the open sea is transformed in a complex way in the shallows -the waves transform from three-dimensional into the two-dimensional ones. Big waves are reduced, while small waves increase in length and height [17].
Second, in the case of moderate waves, calculations of seaworthiness of a vessel using the linear theory yield fully acceptable results [2].
3. The aim and tasks of the study
The aim of present work is to determine the characteristics of a wave field that occurs around the ship, which floats idling at outer anchorage and that is exposed to incident waves of small amplitude. A presence of these characteristics provides for the enhancement of safety of operation of auxiliary vessels at outer anchorage. These data should also be used when carrying out search-and-rescue operations and elimination of oil spills in open water area.
Taking into account the accepted assumptions used to determine the characteristics of a wave field around the ship, the following tasks are to be solved:
- to determine the potential of speeds of perturbed fluid motion caused by the diffraction of waves of small amplitude that are incident under an arbitrary angle to a stationary vessel under conditions of shallow water;
- to define in the specified points around the ship the amplitudes of waves caused by the incident waves and their diffraction on the vessel as an obstacle in wave propagation.
4. Materials and methods of examining the diffraction of regular waves of small amplitude on a stationary ship in shallow waters
4. 1. Statement of diffraction problem
Let us consider the interaction between a stationary vessel of length L floating in shallow water of depth H with the wave of small amplitude, which is incident under an arbitrary angle. Denote length of the wave l, wave amplitude as a, wave propagation speed as c. The fluid is considered perfect, heavy, incompressible, its traffic - potential.
We shall introduce a rectangular Oxyz coordinate system, associated with the vessel. The Oxz plane coincides with the non-perturbed free surface of the fluid. The Ox axis is directed towards the bow; Oy - towards the right board; Oz is vertically upward. The location of the ship relative to the incident wave is defined by the course angle b between the Ox asis and a vector of wave's speed.
Perturbed fluid motion is described by the potential of velocities OE (x,y,z,t). The region of its determining E is limited by the bottom of water reservoir D, wetted surface of a ship S and free surface of the fluid S (part of the plane z=0 outside the waterline of a stationary ship). Given the linearity of the problem, we shall represent the potential OE (x, y, z, t) in the form of a sum:
O1
(x,y,z,t) = O* (x,y,z,t) + Od (x,y,z,t),
(1)
, ag ch [a0 (z - H)]
O =——-r-^sin (a0xcosP + a0ysinb-Gt), (2)
o ch (a0H)
where g=9.81 m/s2; a is the frequency of the incident wave; a0 is the wavenumber (form frequency), which in shallow waters is defined as the only real positive root of transcendental equation:
ag th(aH) = o2. (3)
Let us represent formula (2) in the form:
O* (x,y,z,t) = (x,y,z) cos(ot) + O* (x,y,z) sin(ot); (4)
(x,y.z)l = °s (x,y.z)J
ag ch [a0 (z - H)] [sin 1, ,
= ±— K j \ [(aoxcosp + aoysinp). (5)
o ch (a 0H) [cosJv '
Potential Od (x,y,z,t) of speeds of diffracted wave motion is a solution of the following boundary problem:
9x2
at2"
9y2 + 9z2
d_ !9z
Od (x,y,z,t) = 0, (x,y,z) e E; Od (x,y,0,t ) = 0, (x,y)eS;
aod (x,y,z,t) ao* (x,y,z,t) . .
---- =-----, (x,y,x)eS;
3N 9N V '
dod(x^,H,t)= 0 ( ) . r
lim [gradod (x,y,z,t)] = 0, (x,y,z)eE, r = J x2 + y2
(6)
(7)
(8)
(9) (10)
where N is the outward normal to the surface S.
Similarly (4) Od (x,y,z,t) is represented in the form of a sum:
Od (x,y,z,t) =
= Odd (x,y,z) cos (ot) + Od (x,y, z) sin (ot).
(11)
where O* (x,y,z,t) is the potential of speeds of the system of regular incident waves that the vessel is exposed to; Od (x,y,z,t) is the potential of velocities of wave motion of the fluid, caused by the diffraction of incident wave, on the hull of the ship.
The potential of incident waves' speeds in the coordinate system related to a ship is written in the following form [1]:
Further on, index d in the expressions of diffraction potentials is omitted to simplify the record.
Article [18] stated a general hydrodynamic problem on the potential of perturbed fluid velocities at motion of a ship in waves in a restricted fairway. A stepwise linearization of the given problem was carried out. Solution of the corresponding linear problems by the MAEM method is described.
The above boundary problems differ from the corresponding problems on the potential of radiation [18] only by the form of boundary condition on the vessel hull. Therefore, to determine the potential of velocities of the diffracted motion of fluid, we used the method similar to the one employed in [18].
Consider a structure of the normal derivative on the ship hull. We assume that a ship is an extended body, whose transverse dimensions are small compared to its length while the longitudinal distances at which noticeable changes in the
form of the hull occur, are finite. Given this, we shall refine the type of a normal derivative from the potential on the wetted surface S. It follows from (8) and (11):
9®(x,y,z,t) 9N =
d®c (x,y,z)
=----(
9N
s (at)-
90„
x,y,z)
9N
-sin (a
t).
(12)
According to [18], in the problems on radiation, the expressions for boundary conditions on the ship hull in the case of its vertical fluctuations contain a multiplier cos(N,z), and in the case of cross-horizontal fluctuations - cos(N,y). By analogy, we shall group the constituents of normal derivatives of the potentials that contain these multipliers. According to formulas (5) and (12):
aoc (x,y,z) = BEV 9N C
(x,y,z)
9N
= BE
BO
-bsod;
(13)
(14)
BCV1 _ aga0cos ( N,z)
xsh [a 0 (z - H)]
ach (a0H)
sina0 (xcosb + ysinb),
cosa0 (xcosb + ysinb);
(15)
BOD 1agapcos (N,y) sinb %
xch [a0 (z - H)]
ach (a0H)
cosa0 (xcosb + ysinb),
sina0 (xcosb + ysinb).
(16)
The linearity of a boundary problem for the diffraction potential, as well as an analysis of the structure of the normal derivative from the potential on the wetted surface of a ship indicate that this potential is advisable to represent in the form of a sum:
®(x,y,z,t) = £®, (x,y,z,t).
(17)
problem, uniformly applicable over the entire region, which is filled with fluid.
4. 2. 1. Solving the boundary problems in external zone
As shown in [18], at ship fluctuations by the sinusoidal law, a function of radiation contains not only sine, but also cosine, component, and vice versa. That is why we shall represent potentials 0, (x,y,z,t) in the form of a sum:
0, (x,y,z,t) =
= 0° (x,y,z) cos (ot) + ®s (x,y,z) sin (at). (18)
In the external zone (that is region Eo), we state the boundary problems for amplitude functions Of (x,y,z) and ®s (x,y,z). The problems include conditions of harmony, the boundary conditions at the free surface of fluid So, at the bottom of reservoir D, condition of attenuation of diffracted waves at an infinite distance from the ship. The potentials should also satisfy the fundamental radiations. Boundary conditions at the wetted surface of a ship are not formulated because this surface does not belong in the external zone. It is known only that at diffraction on the ship hull of longitudinal waves, the potential of perturbed velocities is continuous, and its normal derivative is discontinued while passing through the diametral plane of the ship. At diffraction of transverse waves, the normal derivative form the potential is continuous, while the potential is discontinued when passing through the diametral plane. Thus, the bondary conditions of the following form are formulated:
d |±f, x < L/2, 0 < z < H;
—O (x, ±0,z) = 1 c
9y ' [0, \x\ > L/2, 0 <z< H;
—O (x,±0,z) = 0, -¥< x <¥, 0 < z < H; dy
O2 (x, ±0,z) =
l±pc, \x\ < L/2, 0<z< H; 10, \x\ > L/2, 0<z< H;
O2 (x, ±0,z) = 0, -¥< x <¥, 0 < z < H;
—O3 (x, ±0,z) = 0, -¥< x <¥, 0 < z < H; dy
Consequently, the diffracted potentials 0, (x,y,z,t) are determined by the MAEM method, similar to [18].
4. 2. Determining the potentials of diffracted wave motion
According to the procedure for appying the MAEM method, we shall conditionally split the region filled with fluid into the following zones: external, where (y/L)=O(1), and internal, where (y/L)=O(e), e<<1.
Boundary transition s^0 at y and z, fixed in the external zone, transforms the hull into a segment A={-L/2<x<L/2, y=z=0}. Region E is transformed into region E0 (a layer of fluid 0<z<H with a cut-out segment A). Free surface S turns into plane S0 (plane z=0 with a cut-out segment A).
In each zone, separate boundary problems are stated, consequently their solutions asymptotically converge on the boundaries of zones, forming an approximate solution of the
dy °3(x, ±M =
|±fs, |x|<L/2, 0<z<H; ]0, Ixl >L/2, 0<z< H;
O4 (x, ±0,z) = 0, -¥< x <¥, 0 < z < H;
1±ps, Ixl < L/2, 0<z< H;
O4 (x, ±0,z) = ■
0, \x\ > L/2, 0 <z< H.
(19)
(20) (21)
(22)
(23)
(24)
(25)
(26)
Note also that (H/L) = o(ea).
In this case, solutions of the boundary problems are applicable both for the moderate shallow waters (boundary transition a^0 is executed, then:
EV
B
S
OD
B
S
(H/L) = °(1)
and for the considerable shallow waters (then a^1 and (H/L) = O(e)).
We shall employ the Fourier method to the boundary problems for potentials OÇ (x,y,z) and O (x,y,z) in each of the regions y>0 and y<0. Decompose the potentials segment [0; H] along the full orthogonal system of functions [18]:
Zo (z) = N-2 ■ ch [a 0 (z - H)]; Zm (z) = Nm2 ■ cos [am (z - H)] ;
(27)
No = 1 0 2
Nm = 2
1 +
sh (2a0H) 2a H
sin (2amH) + 2a H
, m = 1,2,...
(28)
a1 < a2 <... is the sequence of real positive roots of transcendental equation:
agTg (aH) = -c
(29)
Note that the system of functions (27) is the same that was applied in solving the problems in [4, 5, 11] taking into account the transformations of coordinates.
By using a condition of harmony, the boundary conditions in rectangle 0<z<H, -L/2<x<L/2, as well as the principle of radiation, we shall receive a solution of the problems in the external zone. Then:
- for components of the potentials even along y:
O }=2HZ0 <z)(N0 mû jf {£} z0 (iw-
-SÏZm (Z)i K0
oc
-R) H f ((3 Zm ^
jf (x, ç)L
= + 2H Z0 (z) J J0 (a0R) Jf^jZ0 (i)didX;
(30)
(31)
O2
- for components of the potentials odd along y:
L
=2HZ0 (z)dy (N0 M)i jp; Q Z0 wx-
.yx (z)(K0
)J
fpc(X s)l lp. (X s)j
Zm (ç)dçdX; (32)
O2
oc
=, ¿z (z)dy -J (a 0r ) î jpc
Z0 (ç)dçdx, (33)
where R = ^(x -X)2 + y2, X, Z are the variables of integration along the length of a ship and depth, respectively; Jo, No, K0 are the Bessel, Neumann, and Macdonald functions of zero order of real argument, respectively.
4. 2. 2. Solving the boundary problems in internal zone
Let us perform a boundary transition s^0 when y and z are fixed in the external zone. Introduce the "stretched" coordinates Y=y/e, Z=z/e. With accuracy to low O(e), the motion in the internal zone can be considered as two-dimensional. The potential of speeds of the diffracted wave movement - a harmonic function Y and Z - should satisfy the following differential system:
+iZ22]O(Y,Z)=a (Y,Z) e E (x); (34)
dZ + 0^,0) = 0, k = Y >|b(x); (35)
9O(Y,Z)
9N
= B(Y,Z), (Y,Z)eL(x),
(36)
(37)
where region E(x) is the band 0<Z<h with a cut-out in the form of a frame contour L(x), the coordinates of points of which are (Y, Z); b(x) is the width of the frame along a waterline; h is the depth of water.
The magnitudes b(x) and h are taken in a linear scale of the internal zone.
The given boundary problem does not state boundary conditions at an infinite distance from the ship Y^<», that is, on the outer border of the internal zone. These conditions are set at the stage of converging the solutions of boundary problems of the internal and external zones at their border. Therefore, the solution can be found with accuracy to some arbitrary additive function. The merging technique accepted in the present paper, according to [18], allows us to find a solution without determining this function, as for merging the solutions we apply the asymptotics of speed potential at the outer border of the internal zone.
The asymptotics of speed potentials ®(Y, Z) are written in the form:
- for the components even along Y:
'(Y,Z)Q (I0 ) ch [l (Z - h)] sin (^Y); (38)
- for the components odd along Y: Ood(Y,Z)~±-2v(I0)P(l)ch[I0(Z-h)]cos(I0Y), (39) where
V(l° )= 210h + sh21 0h. (40)
P and Q are the real functions k, defined by formulas:
P (k )1 Q (k1 )J
[[T1 (Y) sin (k1Y) + T2 (Y) cos(k1Y)]}
= 2 J e-k>i0 (Y)
-b 2
where fci(Y) is the equation of contour L(x);
[T1 (Y) cos (k1Y)-T2 (Y)sin (k1Y)][
2
H
b
T = k1A+(Y) + B+ (Y) 1 +
dÇ0 (Y)
dY
T2 = kA+(Y)
dÇ0 (Y).
dY
(42)
(43)
1
0(Y,0) = 0, |Y| > ^b(x).
(44)
For the components of speed potentials odd along Y in accordance with boundary conditions (19), (20), (23), (24), not the potential itself is needed but its derivative along Y. We shall obtain from (38):
—OEV (Y,Z) ~ ± 9Y 1 1
Q (l ) ch[l (Z - h)] cos (l0Y).
(45)
>Q (a
—OEV (±0,z ) = ±-9Y 1 1
for the components odd along y :
' ch [a„ (z -H)] V (a0H); (46)
OOD (±0,z) = + ^ch \aa (z - H)] V^H).
(47)
Therefore,
fc 1 = |Q: ix:^["• (z- H)] V(a.H);
(48)
1 (x, a0 )1
2 k (x, a0 )l
ch[a0 (z - H)] V(a0H).
(49)
A+(Y) is the value of speed potential on contour L(x);
B+(Y) is the value of a normal derivative from the potential on contour L(x).
The value of the normal derivative B+(Y) is defined by the boundary condition. Potential A+(Y) on the contour is unknown as it is actually a solution of the problem. In accordance with the practice of applying the Kochin functions in wave problems [18], instead of A+(Y), formulas (42) and (43) are substituted with the value of the potential at infinite frequency. This potential is the solution to the boundary problem in the internal zone, which consists of equations (34), (36), (37), and instead of condition (35), we accept:
The resulting formulas for the components of speed potential of the diffracted wave motion are derived at substituting the functions fc,s and pc,s from (47) and (48) into (30)-(33).
It should be noted that functions Qc,s and Pc,s, which are icluded in (48) and (49), do not depend on the vertical coordinate. Consider internal integrals in (30)-(33). For the components of potentials that contain functions Z0 (<;), upon substitutions and transformations, we obtain:
J ^ch [a. (ç-H)] Z. (ç*. XMI N3.
0 2 2
where it is written to simplify the record
H , ,r , NT, 2a0H + sh(2a0H) I0 = Jch2[a0(ç-H)]dç = 0 1 0 '.
n ^^n
(50)
(51)
For the components of potentials that contain function Zm (?), we receive:
H ch [«0 (S-H)] Zm (,)d,=Xia^0!) Nmi, 0 2 2
where it is written to simplify the record
Im = Hch [a0 (ç - H)]cos [am (ç - H)]dç = 0,
(53)
4. 2. 3. Merging the solutions. Formulas for a speed potential
In order to merge the solutions, we use a method of boundary merging - "the inner boundary of the external boundary equals the outer boundary of the internal boundary" [19].
In chapter 4. 2. 1, the solution of problem on the diffraction of oblique waves on the hull of a stationary elongated ship in the external zone was derived on the assumption that in band -¥<x<¥, 0<z<H the conditions (19)-(26) are satisfied. The value of functions fc,s and pc,s will be determined from the solution of problem in the internal zone. Let us return to (38) and (45) to the external variables y=eY and z=eZ. In this case, loh proceeds to «0H, V(lo) - to V(«oH), while the potentials and their derivatives at y=±0 will take the form:
- for the components even along y:
because the functions of system (27) are orthogonal. Therefore, other terms in formulas (31) and (32) are converted to zero.
When computing functions Qc,s and Pc,s, we use the corresponding formulas for normal derivatives of potentials at the wetted surface of a vessel, in particular: BCV for Qc; BEV - for Qs; BOD - for Pc; B°?D - for Ps.
Note also that for the components of potentials odd along y in the resulting formulas, instead of Bessel function J0 («0R) and Neumann function N0 («0R) of zero order of real argument, we use first order Bessel function J1 («0R) and Neumann function N1 («0R) [20], respectively.
Solution of the problem on the cosine component of the potential is given in article [21].
4. 4. Determining the characteristics of a wave field around a ship
A theoretical solution of the diffraction problem, derived above, is used for determining the amplitudes of wave in the assigned points around the hull of a ship, which floats idling in the shallow waters.
The equation of a wave profile is written in the following form:
zb = gdt [oE (x-y,0,t)].
(54)
Considering (1), (2), (4), (11), (17), formula (54) takes the form:
zB = rccos (at) + rssin (at), (55)
where the amplitude functions are equal to, respectively:
rc = acos (0) + ^ + q2 +q3 +q4;
rs =asin (0)_qi _q2 _q3 _q4.
The components of diffracted waves qc,s, i=1, 2, 3, 4, taking into account the substitution z=0 and function evenness ch («0H) are determined by formulas:
qc =a2g0E(«oH)J Qc(X,«o)No(«oR)dX; (57)
2g L
q 2 = y
a«(
2g
L
E («oH) }pc (X,«o) ^^ dX; (58)
q3 =ar"~E(«oH) J Qs(X,«o)Jo(«oR)dX; (59)
2g L
qc = y a«gLE («oH) JPs (X, «o )j1%R)dX; (60)
_L 2 L 2
qs = -OfE («oH) J Qc (X,«o) Jo («oR) dX; (61)
2g L
q 2 = _ y
a«(
2g
L
■E(«oH)JPc (X,«o )J1^dX; (62)
q3 = —(E(«(H) J Qs (X,«o)No («(R)dX; (63)
2g L
q 4 = y O^E («(H)JPs (X,«o) ^ dX,
2g L
R
where it is written to simplify the record
E («(H) = VM) ch (_«oH) = 1
j«-^ + 2th («0H)
ch2 («0H)
The amplitude of waves is determined by formula:
5. Results of examining a wave field around a ship
Using the formulas given above, we calculated wave amplitudes in the assigned points of observation around a ship. These points form a grid, coordinates of the points are within |x|<L, |y|<L points inside the contour of the waterline of a ship are excluded. A step by abscissas and ordinates of the grid is 0.0125L. Such an arragement of points is chosen according to the results of numerical experiments; it provides a satisfactory accuracy of subsequent calculations at the smallest volume of source data.
In the grid's nodes, we determine the cosine and sine components of rc ans rs (56), which are consequently employed for the calculation of wave amplitudes (66).
As the object of research we chose a wave field around a bulk carrier of "Zoya Kosmodemyanskaya" type. The hull of this ship is a typical representative of the class of large tonnage bulkers and tankers. Its features are a tuber-like bow, transom stern, a long cylindrical insert with vertical sides, a flat bottom. Principal dimensions of the vessel: length L=201.6 m, width B=19.8 m, draught T=11.73 m.
Calculations of waves' amplitudes were performed for the following characteristics:
- relative water depth H/T=1.1, 1.3, 1.5, 2.0;
- course angles of incident waves p=90°, 120°, 135°, 150 °, 180 °;
- relative wavelength X/L=0.5, 0.6, 0.7, 0.8, 0.9, 1.0.
Using results of the calculations for all combinations of
water depth, length of waves and course angles of waves, we obtained relative wave amplitudes r = r / a. The points of observation are around the hull of the ship in a square field |x|<L, |y|<L. We built lines of equal amplitudes at step Ar = 0.1.
As an example, Fig. 1, 2 show the wave fields around the hull of the ship.
(64)
(65)
(66)
When performing calculations, each frame contour L(x) is given in the form of sets of points. The points divide the contour into fairly small elements, each of which is considered as a straight-line segment.
Functions Qc,s and Pc,s are computed for each frame contour L(x), that is for fixed x. Calculation formulas are transformed considering the replacement of variables in (41) at the stage of merging the solutions and trnasition from the integration to the summation by points.
It should be noted that the employed calculation techniques are not related to any special form of frame contrours.
Fig. 1. Wave field around the hull of the vessel H/T=1.1, VL=0, p=120°
The distributions of relative wave amplitudes are shown around the hull of the ship in a square box |x|/L<1, |y|/L<1
at relative depth of water area H/T=1.1, relative length of waves X/L=0.5, and course angles of incident waave p=120 ° and p=150 respectively.
The following patterns are established: 1. At transverse and oblique wave near the ship from the side of incidence, the standing waves appear. Wave antinodes are at a distance from each other that roughly equals a half the wavelength.
Fig
2. Wave field around the hull of the vessel H/T=1.1, X/L= P=150°
2. The arrangement of antinodes depends on the course angle of waves, in particular:
- at p=90°, the antinodes are located in the region of a middle frame, perpendicular to the incident waves' velocity vector;
- at oblique angles (p=120 135°, 150°), the antinodes are shifted to the stern, the angles between the antinodes and an incident waves' velocity vector are blunt; the shorter the incident wave the larger the shift of antinodes towards the stern. _
3. Increase in the relative wave ordinates r = r/a depends on the course angle of waves and wavelengths. The largest values r = 1.45 are at p=90° and X/L=0.5. The longer the waves and the larger the course angle of waves, the less wave ordinates are. Thus, at p=90° and X/L=1.0 r = 1.12.
4. Among the antinodes, there are zones of reduction in wave coordinates - we observe valuer = 0.8^0.9. It should be noted that r < 1.0 does not mean reducing the level of water, but the moderation of fluctuations in a given area.
5. At the waves incident from the bow (P=180 °), an increase in the wave ordinates of the order r = 1.3 is observed before the ship.
6. The magnitudes of reduction in the wave ordinates in a region of shadow are _approximately r = 0.4 at X/L=0.5, r = 0.6 at X/L=0.8 and r = 0.8 at VL=1.0.
6. Discussion of results of determining a wave field around a ship
An analysis of the calculated wave fields reveals that at transverse and oblique waves near a ship, there occur the standing waves from the side of incidence. In this case:
- the slope of waves in this region increases by about twice compared with the region of hydrodynamic shadow or the region at a considerable distance from the hull of the vessel;
- the arrangement of antinodes in the examined range of lengths and course angles of waves changes little at depth of the water area;
- the depth influences the magnitude of growth and reduction in the wave ordinates, in particular: the lower the depth, the larger the magnitude r and the bigger the difference between the maximum and minimum relative amplitudes is.
It was also established that in the region of a hydrody-namic shadow, the standing waves are missing. In this case:
- location and length of the region of shadow is mainly determined by the magnitude of course angle of waves and the wave length and almost does not depend on the depth of water area;
- the shorter the incident waves the less the depth of water area, the more pronounced is the weakening of fluctuations in the region of a hydrodynamic shadow, but the zone of maximum attenuation is also narrower;
- the magnitude of waves' attenuation in the region of a shadow mostly depends on the wavelengths and course angle of waves, and to a lesser extent on the depth.
The longer the incident waves, the smaller are the changes in the wave field caused by a vessel.
The results obtained are preliminary, as they were received using the linear theory. Their refinement is implied at the next stage of study that will address solution of the problem on the diffraction of oblique waves of ultimate amplitude on the hull of a stationary elongated vessel under conditions of shallow water.
7. Conclusions
1. The present work reports a solution by the MAEM method of the problem on diffraction of oblique waves of small amplitude on a vessel under conditions of significant shallow water. By applying this method, we for the first time obtained the formulas for the potential of speeds of diffracted wave motion of fluid in the assigned points near the hull of the ship. The vessel is considered to be elongated and motionless. A shape of the frame contours of the vessel may take U-shaped, V-shaped, tuber-shaped form. The methodogy of solution and the results obtained represent a further development of the MAEM method for the study and solution of boundary problems on the hydrodynamics of a ship.
2. We obtained the equation of a wave profile for a combination of incident waves and the diffracted waves. The calculations are performed of the amplitudes of waves in the preset points of observation around the ship. The variable parameters are the depth of water area, wave length and a course angle of waves. We analyzed the transformation of the incident waves on the hull of the ship and around it. Examples are given of the distributions of relative wave amplitudes around the hull of the ship.
References
1. Haskind, M. D. Gidrodinamicheskaya teoriya kachki korablya [Text] / M. D. Haskind. - Moscow: Nauka, 1973. - 327 p.
2. Morekhodnost' sudov i sredstv okeanotekhniki. Metody ocenki [Text]: monografiya / I. K. Borodaj (Ed.); FGUP «Krylovskij gosu-dar-stvennyj nauchnyj centr». - Sankt-Peterburg, 2013. - 256 p.
3. Bonet Chaple, R. P. Refraction and diffraction of water waves using finite elements with a DNL boundary condition [Text] / R. P. Bonet Chaple // Ocean Engineering. - 2013. - Vol. 63. - P. 77-89. doi: 10.1016/j.oceaneng.2012.10.006
4. Zhu, S.-P. Diffraction of ocean waves around a hollow cylindrical shell structure [Text] / S.-P. Zhu, L. Mitchell // Wave Motion. -2009. - Vol. 46, Issue 1. - P. 78-88. doi: 10.1016/j.wavemoti.2008.09.001
5. Zhu, S.-P. Combined diffraction and radiation of ocean waves around an OWC device [Text] / S.-P. Zhu, L. Mitchell // Journal of Applied Mathematics and Computing. - 2010. - Vol. 36, Issue 1-2. - P. 401-416. doi: 10.1007/s12190-010-0410-y
6. Aggarwal, A. Irregular Wave Forces on a Large Vertical Circular Cylinder [Text] / A. Aggarwal, M. A. Chella, A. Kamath, H. Bihs, O. A. Arntsen // Energy Procedia. - 2016. - Vol. 94. - P. 504-516. doi: 10.1016/j.egypro.2016.09.223
7. Song, J. Dynamic characteristics between waves and a floating cylindrical body connected to a tension-leg mooring cable placed in a simulated offshore environment [Text] / J. Song, S.-H. So, H.-C. Lim // International Journal of Naval Architecture and Ocean Engineering. - 2016. - Vol. 8, Issue 4. - P. 375-385. doi: 10.1016/j.ijnaoe.2016.05.003
8. Bai, W. Simulation of second-order wave interaction with fixed and floating structures in time domain [Text] / W. Bai, B. Teng // Ocean Engineering. - 2013. - Vol. 74. - P. 168-177. doi: 10.1016/j.oceaneng.2013.07.014
9. Wang, C. Z. Second-order wave radiation by multiple cylinders in time domain through the finite element method [Text] /
C. Z. Wang, S. Mitra, B. C. Khoo // Ocean Systems Engineering. - 2011. - Vol. 1, Issue 4. - P. 317-336. doi: 10.12989/ ose.2011.1.4.317
10. Zhong, Z. Modeling fully nonlinear shallow-water waves and their interactions with cylindrical structures [Text] / Z. Zhong, K. H. Wang // Computers & Fluids. - 2009. - Vol. 38, Issue 5. - P. 1018-1025. doi: 10.1016/j.compfluid.2008.01.032
11. Goren, O. Second-order wave diffraction by horizontal rectangular barriers [Text] / O. Goren, S. M. Calisal // Canadian Journal of Civil Engineering. - 2011. - Vol. 38, Issue 5. - P. 546-555. doi: 10.1139/l11-027
12. Newman, J. N. Second-Order Diffraction in Short Waves [Text] / J. N. Newman // Journal of Ship Research. - 2005. - Vol. 49, Issue 4. - P. 263-273.
13. Semenova, V. Yu. Opredelenie nelinejnyh sil vtorogo poryadka, voznikayushchih pri poperechnoj kachke kontura na tihoj vode v usloviyah melkovod'ya [Text] / V. Yu. Semenova, So Chzho Tu // Morskie intellektual'nye tekhnologii. - 2012. - Issue 2. -P. 22-26.
14. Semyonova, V. Yu. O vliyanii melkovod'ya na znacheniya prisoedinennyh mass i koehfficientov dempfirovaniya sudna [Text] / V. Yu. Semyonova, Tan Htun Aung // Morskie intellektual'nye tekhnologii. - 2011. - Issue 1. - P. 10-14.
15. Kostrov, I. S. Gidrodinamicheskaya prodol'noj kachki sudna, dvizhushchegosya na znachitel'nom melkovod'e [Text] / I. S. Ko-strov // Visnyk ODMU. - 2009. - Issue 27. - P. 49-59.
16. Du, S. X. An improved matching method to solve the diffraction-radiation problem of a ship moving in waves [Text] / S. X. Du,
D. A. Hudson, W. G. Price, P. Temarel // Proceedings of the Institution of Mechanical Engineers, Part M: Journal of Engineering for the Maritime Environment. - 2004. - Vol. 218, Issue 3. - P. 139-152. doi: 10.1243/1475090041737930
17. Lopatuhin, L. I. Vetrovoe volnenie [Text]: uch. pos. / L. I. Lopatuhin; Sankt-Peterburgskij gosudarstvennyj universitet. - 2-e izd., dop. - Sankt-Peterburg: VVM, 2012. - 165 p.
18. Vorob'ev, Yu. L. Gidrodinamika sudna v stesnennom farvatere [Text] / Yu. L. Vorob'ev. - Sankt-Peterburg: Sudostroenie, 1992. - 224 p.
19. Van Dyke, M. Perturbation methods in fluid mechanics [Text] / M. Van Dyke. - London: Academic Press, 1964. - 229 p.
20. Kotlyar, Ya. M. Metody matematicheskoj fiziki i zadachi gidro-aehrodinamiki [Text] / Ya. M. Kotlyar. - Moscow: Vysshaya shkola, 1991. - 208 p.
21. Efremova, N. V. Gidrodinamicheskaya zadacha o difrakcii volneniya, nabegayushchego na sudno v melkovodnoj akvatorii [Text] / N. V. Efremova // Visnyk ODMU. - 1998. - Issue 1. - P. 80-85.
