Next Article in Journal
Improvement of Machine Learning-Based Modelling of Container Ship’s Main Particulars with Synthetic Data
Previous Article in Journal
A Borehole Acoustic Calculation Approach with Gas Hydrate Saturation Inversion in Gas Hydrate-Bearing Sediments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Dimensional Wave Interaction with a Rigid Body Floating near the Marginal Ice Zone

1
School of Naval Architecture and Ocean Engineering, Jiangsu University of Science and Technology, Zhenjiang 212100, China
2
State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(2), 272; https://doi.org/10.3390/jmse12020272
Submission received: 21 December 2023 / Revised: 22 January 2024 / Accepted: 22 January 2024 / Published: 1 February 2024
(This article belongs to the Section Ocean Engineering)

Abstract

:
The interaction problem of waves with a body floating near the marginal ice zone is studied, where the marginal ice zone is modeled as an array of multiple uniformly sized floating ice sheets. The linear velocity potential theory is applied for fluid flow, and the thin elastic plate mode is utilized to describe the ice sheet deflection. A hybrid method is used to solve the disturbed velocity potential; i.e., around the floating body, a boundary integral equation is established, while in the domain covered by ice sheets, the velocity potential is expanded into an eigenfunction series, and in the far-field with a free surface, a similar eigenfunction expansion is used to satisfy the radiation condition. The boundary integral equation and the coefficients of the eigenfunction expansions are solved together based on the continuous conditions of pressure and velocity on the interface between the sub-domains. Extensive results for the equivalent Young’s modulus of the ice sheet array and hydrodynamic force on the body are provided, and the effect of individual ice sheet length as well as wave parameters are investigated in detail.

1. Introduction

When the free surface wave propagates into the ice-covered region, the large ice sheet would bend to let the energy pass through. As the wave amplitude is larger than a specific one, the large ice sheet would break into small ones, forming the marginal ice zone (MIZ). To avoid additional ice resistance, a ship may navigate along the edge of the MIZ. In such a situation, the disturbed waves by the ship may be reflected by the MIZ, and exerting on the ship again, which will lead to different features of the hydrodynamic force compared to that in open water.
Sea ice has been an area of research interest, with studies mainly focusing on wave–ice interaction. The research has focused on the propagation of waves through ice sheets. The authors of [1] explored wave dispersion relationships and the viscous attenuation of sea ice cover in wave diffusion to a marginal ice zone. They discussed the dispersion relationship of waves under bending and compression, as well as the viscosity attenuation of ice sheets. The authors of [2] calculated the equivalent modulus of elasticity using empirical equations for the elastic modulus and dispersion relationship of ice. Using this modulus and the dispersion relationship equation based on the linear theory of wave motion under infinite elastic floating plates, the wave propagation velocity in the sea ice region was estimated. Some researchers have shown that the properties of ice sheets have an important influence on wave scattering. Ref. [3] considered the effect of a draft on wave scattering between two floating flexible plates. Ref. [4] studied the scattering of ice-coupled waves by cracks in ice sheets, analyzed the effects of cracks and pressure ridges on scattering, and determined the reason for the minimum reflection coefficient. Ref. [5] studied the propagation of curved gravity waves in sea ice of different thicknesses and analyzed the effects of ice size and draft on wave energy distribution. Ref. [6] studied the effect of plate thickness on wavenumber and critical angle. Ref. [7] studied the interaction of waves with a fixed horizontal plate on a fluid surface, providing solutions for the scattering of incident waves from an infinite plate with a constant finite width and a finite parallelogram plate. Ref. [8] deduced the transport equation for flexural-gravity wave propagating under ice sheets of variable thickness in space.
The problems of wave attenuation and scattering has been addressed by various studies. For example, residual calculus techniques were used by [9] to solve the scattering problem of water waves by the edge of a semi-infinite ice sheet in a finite depth of water. The case where the obliquely incident plane wave is from the open sea region and the complementary problem where the wave is incident from the ice-covered region are both considered and exact solutions are obtained. Ref. [10] proposed a linear Boltzmann equation to simulate wave scattering in marginal ice regions. Ref. [11] proposed a method to address the wave scattering problem of an ice sheet of variable thickness. Ref. [12] presented an asymptotic solution for a single crack in an ice sheet, which shed light on the existing wide spacing approximation. As for wave scattering caused by randomly distributed pack ice, the authors of [13] proposed a diffusion approximation method to decompose the wave action density function into a transmission part and a scattering part, and investigate the effects of various factors on wave scattering. Furthermore, ref. [14] presented a three-dimensional model for wave attenuation caused by ice pack scattering, which takes into account the horizontal dimensions of individual pack ices. Ref. [15] proposed a wave dissipation model to study the relationship between ice thickness and wave frequency. Ref. [16] used a direct phase resolution simulation to study the propagation and scattering of nonlinear ocean waves in broken sea ice. They calculated the local wave spectra of the entire ice field for various wave states and floating ice configurations, and performed least square fitting to extract the spatial attenuation rate as a function of wave frequency. Through these studies, the problem of wave attenuation and scattering has been addressed from various aspects, which facilitated our understanding of wave–ice interactions in marginal ice regions.
Experimental studies on the interaction between waves and ice have been conducted. The authors of [17] experimentally studied the dispersion relationship and the attenuation of waves propagating through ice sheets composed of a mixture of oil and thin ice. Ref. [18] performed an experimental study of the kinematic response between isolated ice floes and waves. The results showed that the scale effect, ice pack size, ice pack azimuth, and surface roughness had no effect on the heave motion. Ref. [19] experimentally studied wave propagation in ice/pancake ice, pancake ice, and crushed ice floes, obtained the wave number and attenuation of several monochromatic waves, and determined the equivalent viscoelasticity required for different ice sheets. Ref. [20] proposed an experimental model for ice floes, which used thin plastic sheets with different material properties and thicknesses to study the effects of period and wave steepness on incident and transmitted waves. Ref. [21] studied the effects of incident wave period and steepness on reflection and transmission, and measured the wave field reflected and transmitted by thin floating elastic plates. Ref. [22] conducted experiments on the attenuation and scattering of surface waves in various ice sheets to study the effects of ice thickness, wave frequency, and flexural stiffness on the attenuation rate. The authors of [23] determined the equivalent modulus of elasticity through the wavenumber and dispersion relationship experimentally. They proposed the empirical relationship between equivalent elasticity in wavelength, floating ice size, and length scale. The authors of [24] conducted experiments to study monochromatic waves of different amplitudes and frequencies propagating through three types of ice sheets: thin ice, crushed ice floes, and grease ice. They analyzed the spatial evolution of wave attenuation and phase velocity variation. The wave dispersion and attenuation under different ice properties and wave characteristics were studied.
In recent years, the study of wave, body, and ice interactions in the Arctic region has attracted significant research interest. The authors of [25] studied the wave-induced motion of a floating body on the water surface between two semi-infinite ice sheets. They divided the fluid domain into sub-regions and solved the unknown coefficients by matching the eigenfunction expansions of each sub-region. They discovered that the added mass, damping coefficient, and excitation force of the body oscillate with the wave number. Ref. [26] investigated the hydrodynamics of a cylinder undergoing large amplitude oscillation in ice-covered water using the multipole expansion method. The study examined the effects of properties of the ice sheet as well as the frequency of oscillation. The authors of [27] studied wave diffraction and radiation of a body in a polynya enclosed by an ice sheet extending to infinity. They utilized the eigenfunction expansion method to solve the problem. Recently, ref. [28] presented a hybrid method to solve the interaction of waves with a three-dimensional floating structure in a polynya. These studies contribute to a better understanding of wave–ice–body interactions, which can be helpful for improving our ability to navigate and operate in ice-covered regions.
There are many large and small pieces of ice in the MIZ, which may alternate the way of wave propagations. This work aims to study the wave interactions with a body floating in the vicinity of the MIZ. The disturbed velocity potential was solved using the hybrid method, which involved the establishment of a boundary integral equation around the floating body. The velocity potential within the domain covered by ice sheets was expanded into an eigenfunction series, while the similar eigenfunction expansion was used to satisfy the radiation condition in the far-field free surface part. The continuous conditions of pressure and velocity on the interface between the sub-domains were used to solve the boundary integral equation and the coefficients of the eigenfunction expansions. The study focuses on the influence of the length of an individual floating ice sheet as well as the wavelength. The effects of floating ice array on the hydrodynamics of the body are studied in detail.

2. Numerical Method

In Figure 1, the interaction of waves with a body floating near the semi-infinite floating ice array is depicted. The width and draft of the body are a and b , respectively. The origin O of the cartesian coordinate system is defined at the center of gravity of the body. The x-axis is horizontally to the right along the mean free surface, the z-axis is vertically upward, and y-axis pointing into the paper. The motion of the body is excited by an incident wave which propagates along the x-axis.
When the body was floating on the water surface and in the stationary state, the center of buoyancy of the body was passed through by the z-axis. The density of an inviscid, incompressible, and homogeneous fluid is assumed to be ρ , the depth is H , and its motion is irrotational. The flow of fluid can be described through the velocity potential Φ . There are several individual ice sheets floating on the right extend from x 2 to infinity, as shown in Figure 1, forming a floating ice array. The floating ice array can be considered as a semi-infinite continuous ice sheet with the equivalent elastic modulus E e q , the average thickness h , draught d , density ρ i , Young’s modulus E , Poisson’s ratio ν . The equivalent process of such a model for the floating ice array has been discussed by ref. [23] based on both physical measurements and theoretical approach.
We assume that the amplitude of incident wave is small relative to the wavelength and width of body, so the linearized velocity potential theory can be used. The wave motion is sinusoid in time with the frequency ω ; the total velocity potential can be obtained as [29].
Φ x , z , t = Re α 0 ϕ 0 x , z e i ω t + Re i = 1 3 i ω α i ϕ i x , z e i ω t .
In Equation (1), α 0 is the amplitude of the incident wave potential. The incident potential ϕ I and diffracted potential ϕ D are included in the ϕ 0 . The radiation potential ϕ i ( i = 1 ,   2 ,   3 ) is generated by oscillation of the floating body. Degrees of freedom include translations in the x and z directions, respectively, and rotation about the y-axis. The potential ϕ i in the flow field should satisfy Laplace’s equation:
2 ϕ i = 0 , i = 0 ,   1 ,   2 ,   3 .
The combination of the linearized dynamic and kinematic free surface boundary conditions provides
ϕ i z ω 2 ϕ i = 0 , < x < x 2 ,   z = 0 ,
where g   is the gravitational acceleration. On the vertical plane of the draft of the floating ice array, the impermeable condition is
ϕ i x = 0 , x = x 2 , d z 0 .
In addition, the impermeable conditions of the floating body are
ϕ 0 n = 0   and   ϕ i n = n i ,   i = 1 ,   2 ,   3 ,
where n 1 , n 2 are the x-axis and z-axis components of the unit normal vector n pointing into body, and n 3 = z z n 1 x x n 2 with ( x , z ) as the centre of rotation.
The seabed boundary condition can be written as
ϕ i z = 0 , < x < + , z = H .
The radiation conditions in the far-field are
l i m x ϕ i δ 0 , i ϕ I x k 0 1 ϕ i δ 0 , i ϕ I = 0 ,
l i m x + ϕ i x + κ 0 2 ϕ i = 0 ,
where δ 0 , i is the Kronecker delta function, and δ 0 , i = 1 if i = 0; otherwise, δ 0 , i = 0 . In the dispersion equations of regions Ω 1 and Ω 2 , the purely positive imaginary roots k 0 1 and κ 0 2 will be discussed in Section 3.1 below.
The floating ice array is considered as a continuous ice sheet with an equivalent elastic modulus. Boundary conditions at the surfaces under the ice sheet can be written as [30]
L 4 x 4 m ω 2 + ρ g ϕ i z ρ ω 2 ϕ i = 0 , x x 2 , z = d ,
where L = E e q h 3 / [ 12 1 υ 2 ] is the effective flexural rigidity of the elastic thin plate, and m = h ρ i is the mass per unit area. E e q is the equivalent modulus of ice, and it can be obtained by an empirical formula [2],
E i m = E e q / E ,
where E i m can be obtained by ref. [23]
log 10 E i m = 6.88 1 + 1.281 e I G .
Here, I G consists of two dimensionless length scales
I G = 2.921 h I c 0.012 log 10 I i I c λ o w I i 0.001
with
I c = E h 3 12 ρ w g 1 ν 2 1 / 4 .
Here, h and I i are the thickness and length of an individual plate, respectively. I c is a characteristic length relating to h ; λ o w is the wavelength of the incident wave in open water. E and ν are the intrinsic Young’s modulus and Poisson’s ratio of the plate, respectively. ρ is the density of water, and g is the gravitational acceleration.

3. Solution Procedures

3.1. Solution of Velocity Potential of Each Sub-Domain

Based on Green’s identity, the boundary integral equation in each sub-domain Ω 1 , Ω 2 , and Ω 3 can be written as
S 1 ϕ i 1 δ 0 , i ϕ I ψ m 1 n ϕ i 1 δ 0 , i ϕ I n ψ m 1 d S = 0   in   Ω 1 ,
S 2 ϕ i 2 ψ m 2 n ϕ i 2 n ψ m 2 d S = 0   in   Ω 2 ,
and
α p ϕ i 3 p = s 3 G p , q ϕ i 3 q n q G p , q n q ϕ i 3 q   d S q   in   Ω 3 .
It should be noted that the boundary of region Ω 3 includes the body surface S 0 , the free surface S F , the vertical interface 1 and 2 , and the seabed S B . Thus, S = S 0 + S F + 1 + 2 + S B .
In Equations (14)–(16), the velocity potential ϕ i is expressed by the eigenfunctions as
ϕ i 1 = δ 0 , i ϕ I + m = 0 R i , m ψ m 1 ,
ϕ i 2 = m = 2 T i , m ψ m 2 ,
and
ϕ i 3 x j , z = m = 0 C i , m j ψ m 3 ,   j = 1 ,   2 ,
where
ϕ I = g i ω e k 0 1 ( x x 1 ) cos k 0 1 z + H cos k 0 1 H ,
ψ m 1 = e k 0 1 x 1 x cos k m 1 z + H cos k m 1 H ,
ψ m 2 = e κ 0 2 x 2 x cos κ m 2 h + z cos κ m 2 h d ,
ψ m 3 = cos k m 3 H + z cos k m 3 H .
The unknown coefficient T i . m needs to be determined. k m 1 and k m 3 satisfy the dispersion equation for free surface
k m i tan k m i H = ω 2 g ,   ( i = 1 ,   3 ) ,
where k 0 1 and k 0 3 are purely positive imaginary roots and k m 1   and   k m 3 m = 1 ,   2 ,   . . . are the purely positive real roots. It should be noted that κ m 2 is the roots of the dispersion equation for the ice sheet [2]
κ m 2 tan κ m 2 H d = ρ ω 2 L κ m 2 4 + ρ g m 2 ω 2 .
Similar to the method of [25], Equations (14) and (15) can be rewritten as
H 0 ϕ i 1 ψ m 1 x ϕ i 1 x ψ m 1   d z = δ 0 , i H 0 ϕ I ψ m 1 x ϕ I x ψ m 1   d z   x = x 1
H d ϕ i 2 ψ m 2 x ψ m 2 ϕ i 2 x   d z + L ρ ω 2 3 ψ m 2 x 2 z 2 ϕ i 2 z x 4 ψ m 2 x 3 z ϕ i 2 z z = d = 0   x = x 2 .
Similar to that in [31], Equation (16) can be rewritten as
α p ϕ i 3 p = S 0 G p , q ϕ i 3 q n q G p , q n q ϕ i 3 q   d S q + S F G p , q ϕ i 3 q n q G p , q n q ϕ i 3 q   d S q + m = 0 C i , m 1 H 0 G p , q ξ ψ m 3 ζ ξ = x 1 d ζ m = 0 C i , m 2 H 0 G p , q ξ ψ m 3 ζ ξ = x 2 d ζ + H d G p , q ϕ i 3 q ξ ξ = x 2 d ζ .
where p ( x , z ) is the field point in the field; q ( ξ , ζ ) is the source point. Then, the simple source function is defined as
G p , q = ln 1 / r 1 + ln 1 / r 2 ,
where
r 1 = x ξ 2 + z ζ 2 ,
r 2 = x ξ 2 + z + ζ + 2 H 2 .
In order to solve the velocity potential of each fluid sub-domain, we use the matching method in [31]. According to the continuous condition on the interface
ϕ i j = ϕ i 3 ,   ϕ i j x = ϕ i 3 x   x = x j , j = 1,2 ,
Equations (14) and (15) can be rewritten as
H 0 m = 0 C i , m 1 ψ m 3 ψ m 1 x ϕ i 1 x ψ m 1   d z = δ 0 , i H 0 ϕ I ψ m 2 x ϕ I x ψ m 2   d z   x = x 1
and
H d m = 0 C i , m 2 ψ m 3 ψ m 2 x ϕ i 2 x ψ m 2   d z + L ρ ω 2 3 ψ m 2 x 2 z 2 ϕ i 2 z x 4 ψ m 2 x 3 z ϕ i 2 z z = d = 0 x = x 2 .
When p is located on the vertical interface x = x j ( j = 1 ,   2 ) , Equation (28) can be converted to
π C i , m j cos 2 k m 3 H H 2 + sin 2 k m 3 H 4 k m 3 = H 0 ψ m 3 z   d z S 0 G p , q ϕ i 3 q n q G p , q n q ϕ i 3 q   d S q + H 0 ψ m 3 z   d z S F G p , q ϕ i 3 q n q G p , q n q ϕ i 3 q   d S q + m = 0 C i , m 1 F F m , m x , x 1 m = 0 C i , m 2 F F m , m x , x 2 + m = 2 T i , m G G m , m 2 x , x 2 ,
where
F F m , m x , ξ = H 0 ψ m 3 z   d z H 0 G p , q n q   ψ m 3 ζ   d ζ ,
G G m , m 2 x , ξ = H 0 ψ m 3 z   d z H d G p , q ϕ i 2 q n q   d ζ .

3.2. Solution of Hydrodynamic Force Acting on the Body

Once the velocity potentials have been solved, the pressure on the body can be obtained through the linear Bernoulli equation. Integrating the dynamic pressure over the mean wetted body surface, the hydrodynamic force acting on the floating body can be calculated out. The hydrodynamic force can be divided into three parts; i.e., the wave exciting force from the scattering potential, and the added mass and damping coefficient from the radiation potential. Additionally, due to changes in buoyancy during body oscillation, there may also be a hydrostatic force. Complex motion amplitudes can be determined using linear equations [32]
k = 1 3 ω 2 m j k + μ j k + i ω λ j k + C j k   α k = f E , j α 0 ,   j = 1 ,   2 ,   3 ,
where k = 1 is the sway mode of the body, k = 2 is the heave mode of the body, and k = 3 is the roll mode of the body. The body mass is m j k , and the hydrostatic restoring coefficient is C j k . μ j k is the added mass coefficient, and λ j k is the damping coefficient. They can be computed through the following equation
μ j k i λ j k ω = ρ S 0 ϕ k n j d S .
In Equation (38), f E , j is the exciting force of the wave with unit amplitude. It can be obtained by
f E , j = i ω ρ S 0 ϕ 0 x , z n j d S

4. Numerical Results

In the following sections, we use the density of water ρ = 1025   k g / m 3 , the acceleration due to gravity g = 9.8   m / s 2 , and the width of the body a = 10   m as the basic parameters, and the results will be expressed in the dimensionless form.

4.1. The Wave Dispersion Relation Beneath the Floating Ice Array

In this section, the effects of wavenumber under the floating ice arrays which consisted of individual ice sheets of different sizes are discussed. The ice thickness is h = 0.1 , the draft of ice is d = 0.09 , and the quality of the ice is m = 0.09 . Six individual ice sheet lengths are considered; i.e., I i = 4 ,   2 ,   1 ,   0.6 ,   0.3 ,   0.15 .
Figure 2 shows the wavenumber in the water covered by individual ice sheets with different lengths, against σ = ω 2 / g a . It can be found that the wavenumber k increases with σ . For a given σ , the wavenumber k decreases with the increase of individual ice sheet length I i . Depending on the individual ice sheet length, the wavenumber may become larger or smaller than that for open water. In fact, when the wavelength of the incident wave is shorter than the size of the individual ice sheet, complementary waves emerge beneath it. This leads to an increase in the wavenumber for smaller-sized individual ice sheets as the wavelength decreases. This phenomenon is observed in ref. [11].

4.2. Equivalent Flexural Stiffness of the Floating Ice Array

The flexural stiffness L of floating ice arrays is given as
L = E e q h 3 / 12 1 υ 2 ,
where E e q is the equivalent modulus of elasticity. E e q can be obtained by Equation (10) [2]. h is the average thickness of the floating ice array and υ is the Poisson’s ration. It may be noticed that E e q = E when the floating ice array becomes the continuous ice sheet.
The effect of individual ice sheet length on the flexural stiffness is investigated according to Equation (41). The average thickness of the floating ice array is h = 0.1 , Poisson’s ratio is υ = 0.3 , and the sizes of individual ice sheets are I i = 4 ,   2 ,   1 ,   0.6 ,   0.3 ,   0.15 . The flexural stiffnesses of floating ice arrays with different I i is shown in Figure 3. It can be found that I i has a significant effect on the flexural stiffness. Because the values have different orders of magnitude, they are shown in different sub-figures. From the figure, it is observed that the flexural stiffness for I i = 4 ,   2 decreases with the increase of wavelength. However, for small-sized ice sheets, the flexural stiffness of the floating ice array increases with the wavelength, namely for I i = 1 ,   0.6 ,   0.3 ,   0.15 .
The variation of flexural stiffness with wavelength is mainly due to the change in the equivalent elastic modulus E e q , which is defined as
E e q = E i m E ,
as given in Equation ( 10 ) . The experience value E i m is determined by empirical coefficient I G in Equation ( 11 ) , in which the terms of ( h / I c ) and l o g 10 ( I i / I c ) are mainly included. Here, the average thickness of the floating ice array is taken to be h = 0.1 , the characteristic length I c is the fixed value which is solved by Equation ( 13 ) , then ( h / I c ) is also a fixed one. Thus, the main variable is the size of individual ice sheet I i . When the length of individual ice sheet I i is more than the characteristic length I c , the ratio ( I i / I c ) > 1 . At this moment, the correction factor I G directly influences the trend of variation in flexural stiffness. Thus, the flexural stiffness decreases with the increase of the wavelength, as shown in Figure 3a,b. While for a small I i , the flexural stiffness increases with the wavelength.
Then, the effect of wavelength on flexural stiffness is considered. As shown in Figure 4, the flexural stiffness L increases with the length of individual ice sheet I i in all cases, as it can be expected.

4.3. Wave Transmission and Reflection by the Floating Ice Arrays

Wave transmission and reflection by the floating ice array are considered in this section. The average thickness of a floating ice array is h = 0.1 , the draft of ice is d = 0.09 , and the mass per unit area is m = 0.09 . The reflection coefficients K R and transmission coefficients K T in case of different sizes of ice sheets are shown in Figure 5 and Figure 6. When σ < 0.8 , K R rises quickly, and K T decreases correspondingly. When 0.8 < σ < 1.8 , K R shows a fluctuation, and for I i = 1 ,   0.6 ,   0.3 ,   0.15 , K R is nearly indistinguishable. Within the same frequency range, we can observe the similar fluctuation of K T . It is found that the transmission coefficient in the case of small-sized individual ice sheets is greater than that of large-sized individual ice sheets from Figure 6. This is because for a small-sized individual ice sheet, the flexural stiffness of the floating ice array is small.

4.4. The Features of Hydrodynamic Forces on the Body near the Floating Ice Array

We investigate the effect of changes in the size of an individual ice sheet on the hydrodynamic forces of a square box floating near the MIZ, and analyze the degree of influence on the hydrodynamic force of the square box by changing the mean thickness and draft of the floating ice array.

4.4.1. Verification of Hydrodynamic of a Body Floating near the MIZ

Based on the numerical calculation method of a floating body in two-dimensional MIZ, the square box is used for numerical calculations in the study. Ref. [25] solved the hydrodynamic problem of the interglacial channel between two semi-infinite floating ice sheets, and the analytical solution of the floating square box with width a and draft b was obtained. Thus, we used the same parameters, the width of the square box a = 1 , draft b = 0.5 , rotation center ( x , z ) = ( 0 , b / 2 ) , the depth of water H = 10 , the distance between the square box and floating ice array x = 5 , the mean ice thickness h = 0.1 , the draft of ice d = 0.09 , the quality of the ice m = 0.09 , and the flexural stiffness of the continuous ice sheet L = 4.5582 , to conduct a comparative verification.
The following Figure 7, Figure 8 and Figure 9 show the added mass, damping coefficients, and wave exciting force of the square box, respectively, which is floating near the MIZ. Additionally, the far-field damping coefficient and far-field wave exciting force are shown in Figure 8 and Figure 9, respectively, which are calculated by the far-field formula [25]
λ k j = ρ ω ( Q 0 ( 1 ) C g ( 1 ) R k , 0 * R j , 0 + Q 0 ( 2 ) C g ( 2 ) T k , 0 * T j , 0 ) , ( k , j = 1 ,   2 ,   3 )
and
f E , j = 2 i ρ g R j , 0 C g ( 1 ) , ( j = 1 ,   2 ,   3 )
where
Q 0 ( 1 ) = ω g ,
Q 0 ( 2 ) = ρ ω [ L ( κ 0 ( 2 ) ) 4 + ρ g ] [ L ( κ 0 ( 2 ) ) 4 + ρ g m ω 2 ] 2 ,
and C g ( 1 ) is the group velocity of the propagating wave in the open water which can be given as
C g ( 1 ) = g k 0 ( 1 ) t a n h ( k 0 ( 1 ) H ) ,
C g ( 2 ) is the group velocity of a propagating wave under the ice sheet in subregion Ω 2 which can be given as
C g ( 2 ) = d ω d ( i κ 0 ( 2 ) ) = i ω 2 κ 0 ( 2 ) ( 1 + 2 κ 0 ( 2 ) ( H     d ) sin [ 2 κ 0 ( 2 ) ( H     d ) ] ) + 2 L ( κ 0 ( 2 ) ) 3 ω L ( κ 0 2 ) 4   +   ρ g     m ω 2 L ( κ 0 ( 2 ) ) 4   +   ρ g L ( κ 0 ( 2 ) ) 4   +   ρ g     m ω 2 .
Figure 7. The verification of added mass [25]. ( a ) sway; ( b ) heave; ( c ) roll; ( d ) sway–roll.
Figure 7. The verification of added mass [25]. ( a ) sway; ( b ) heave; ( c ) roll; ( d ) sway–roll.
Jmse 12 00272 g007
Figure 8. The verification of damping coefficient [25]. (a) sway; (b) heave; (c) roll; (d) sway–roll.
Figure 8. The verification of damping coefficient [25]. (a) sway; (b) heave; (c) roll; (d) sway–roll.
Jmse 12 00272 g008
Figure 9. The verification of wave exciting force [25]. ( a )  sway; ( b )  heave; ( c )  roll.
Figure 9. The verification of wave exciting force [25]. ( a )  sway; ( b )  heave; ( c )  roll.
Jmse 12 00272 g009aJmse 12 00272 g009b
It should be noticed that the coupling items μ j k and λ j k are both zero, because of the symmetry of the floating body with respect to the x-axis; thus, zero terms are not given in the following figures. As shown in Figure 7, Figure 8 and Figure 9, we have discovered that the results of mesh size 1 M = 70 , N 0 = 90 , N F = 180 , mesh size 2 M = 100 , N 0 = 180 , N F = 360 , and the case of [25] basically coincide when the frequency σ is lower than 0.6 , where M is the number of truncation terms in the vertical interfaces, N 0 is the number of body surface meshes, and N F is the number of free surfaces meshes. For the small wavenumber, the thickness of the ice sheet is relatively small compared to the wavelength, which can result in the negligible effect of ice thickness on hydrodynamics. Then, whether it is two ice sheets or a single ice sheet, the result is the same. However, when the frequency σ is greater than 0.6 , the oscillation trends of the curves are basically the same, but the amplitudes of mesh size 1 and mesh size 2 are only about half of the result of [25]. This is because the square box is only floating near the floating ice array in this study as shown in Figure 1, so a portion of the incident wave is reflected by the left side of the square box to an infinite distance on the left, and the other part of the incident wave passes through the square box and is reflected by the floating ice array to the right of the square box, generating hydrodynamic effects. While in the study of [25], where the square box is floating between two semi-infinite ice sheets, the incident wave is reflected by the left side of square box to the semi-infinite ice sheet on the left. Similarly, the waves passing through the square box are reflected by the right ice sheet to the right side of the square box; the ice sheets on both sides have hydrodynamic effects on the square box. However, only one side of the ice sheet has a hydrodynamic impact on the square box in our study. For the far-field calculation results, they are basically consistent with the results of the two types of meshes. So, the amplitude of our numerical results are roughly half of that. This difference can be observed in Figure 8 and Figure 9. Due to the fact that the positions of the rigid body in the studies are the same, the propagation distances of waves between the floating body and the ice sheet are also the same, resulting in corresponding to the maximum and minimum values at equal frequencies. In addition, due to the interaction of wave–ice sheet–body, the rigid body is subjected to wave forces opposite to the acceleration in the sway direction, resulting in positive and negative values of the added mass in sway and sway–roll modes, which is similar to [33]. This can be seen in Figure 7a,d.
Based on the above analysis of Ren’s and our results, the similarities of oscillation trends and amplitude of the curve are obtained, and the rationality of the numerical method in our study is proved. In following sections, mesh size 1 is adopted for the calculations to improve computational efficiency.

4.4.2. The Influence of the Sizes of Individual Ice Sheets

The hydrodynamic forces of a square box near floating ice arrays in the case of the sizes of individual ice sheets I i = 4 ,   2 ,   1 ,   0.6 ,   0.3 ,   0.15 under the action of incident waves are calculated. The width and draft of the square box are a = 1 and b = 0.5 , respectively, and the parameters of the individual floating ice sheets are: average thickness h = 0.1 , the draft d = 0.09 , the mass m = 0.09 . It is noted that the distances between vertical interfaces and origin point x 1 = x 2 = 5 . The results for the continuous ice sheet and open water are also provided for comparison.
Figure 10, Figure 11 and Figure 12 show the added mass, damping coefficient, and wave exciting force of the square box floating near the floating ice arrays, respectively. The incident wave passes through the square box, with some being reflected by the floating ice array or the continuous ice sheet to the right wall of the square box, and the other passing through the ice cover, resulting in oscillation of the added mass and damping coefficient as the wave frequency increases, and the extreme values of the oscillation curve is affected by the individual size of the floating ice array.
As shown in Figure 10 and Figure 11, when σ < 0.5 , the results for all cases are very similar. In fact, when the wavelength of the incident wave is relatively longer to the average ice thickness, the hydrodynamic forces on the square box are very small. When 0.5 < σ < 4 , the hydrodynamic forces of a continuous ice sheet are significantly larger than that in the case of another floating ice sheet. This means that the interaction between the square box and waves which are reflected by the floating ice array begins to strengthen. Due to the equivalent elastic modulus of I i = 4 ,   2 being close to that of the continuous ice sheet, the results for floating ice arrays consisting of individual ice sheet I i = 4 ,   2 are also close to that in the case of the continuous ice sheet, while the hydrodynamic forces caused by floating ice arrays of individual ice sheet I i = 1 ,   0.6 ,   0.3 ,   0.15 are closer to that in the case of open water. As wave frequency increases, the wavelength is much smaller than the length of the individual ice sheet; this leads to extreme values in the case of floating ice arrays which consisted of individual ice sheet I i = 4 ,   2 ,   1 ,   0.6 ,   0.3 ,   0.15 , exhibiting a gradual increase when σ > 4 . Moreover, at a high frequency, the shorter the length of an individual ice sheet, the greater variation in extreme values. This phenomenon can be observed in Figure 10 and Figure 11.

4.4.3. The Influence of the Individual Ice Thickness

The results for the floating ice arrays for two groups with (1) the size of individual ice sheet I i = 4 and average thickness of ice h = 0.1 ,   0.15 and (2) I i = 0.6 , h = 0.1 ,   0.15 are presented, respectively. The main parameters such as the width of the square box a = 1 , the draft of square box b = 0.5 , and the depth of water H = 10 are adopted in the following calculation. Meanwhile, the results of that in open water are also provided.
As shown in Figure 13, Figure 14 and Figure 15, the added mass, damping coefficient, and wave exciting force on the square box can be categorized as two sets of comparative results, respectively. When σ < 0.8 , it is difficult to distinguish the differences between the curves. It indicates that the impact on hydrodynamic forces is minimal at long wavelengths. When 0.8 < σ < 4 , the hydrodynamic force caused by floating ice arrays with the average ice thickness h = 0.15 is more significant than that caused by that with h = 0.1 . This is because the wavelength decreases and the influence of the ice thickness begins to become stronger. When σ > 4 , for the group (1), the extreme values of hydrodynamic forces with h = 0.1 ,   0.15 both begin to decrease. On the other hand, as the frequency increases, the extreme values with h = 0.1 ,   0.15 begin to rise. This is because when the wavelength is small enough, the effect of the length of an individual ice sheet becomes vital. It is not surprising that the large-sized individual ice sheet will be seldom affected relative to the small-sized one. In addition, as shown in Figure 13 and Figure 14, the trend of rapid change of the added mass and damping coefficient for h = 0.15 is so remarkable, compared with that of h = 0.1 .
In fact, the thickness of the floating ice sheet shall be related to its stiffness. When h increases, it will impede short waves to penetrate the floating ice arrays. Then, it results in greater reflection and scattering waves. Moreover, as the thickness of ice increases, so does its mass, which will directly enhance its wave exciting force, and greater resistances are excited for body motion. Thus, it enforces more complex interactions between the ice array and floating bodies.

4.4.4. The Influence of the Ice Array Draft

The hydrodynamic forces of the square box with a draft considered or not for floating ice arrays are computed and discussed in this section. The calculation parameters are set up as that the width of square box a = 1 , the draft of square box b = 0.5 , rotation center ( x , z ) = ( 0 , b / 2 ) , the depth of water H = 10 , the distance between the original point and vertical interfaces x 1 = x 2 = 5 , the average thickness of ice h = 0.1 , the mass of floating ice array m = 0.09 , and the length of individual ice I i = 0.6 ,   4 . When the ice sheet draft is included, it should be d = 0.09 . Otherwise, d = 0 .
As shown in Figure 16, Figure 17 and Figure 18, the added mass, damping coefficient, and wave exciting force of a square box with I i = 0.6 ,   4 under different drafts are presented. When σ < 0.8 , it is undistinguished for the results of floating ice arrays of lengths of individual ice I i = 0.6 ,   4 with d = 0 and d = 0.09 . When 0.8 < σ < 5 , the peaks of floating ice arrays of I i = 4 with d = 0.09 begin to increase. In the case of d = 0 , the variation of peak values of added mass and the damping coefficient is seemingly small. For I i = 0.6 , the difference in the results of d = 0 and d = 0.09 is not obvious. When σ > 5 , the results of I i = 4 with d = 0 ,   0.09 begin to decrease. On the contrary, for I i = 0.6 , the results begin to increase.
It indicates that, for the floating square box near the floating ice arrays, the draft influence on the hydrodynamic forces will be different under different length scales of individual ice sheets. When near that consisting of a large-sized individual one, it has a greater impact on forces with the ice draft than that with non-draft. However, for floating ice arrays consisting of small-sized ones, the ice draft’s influence is relatively small.

4.4.5. The Impact of an Individual Ice Sheet on Wave Elevation

We further investigate the total wave elevation during the entire field. The computational parameters are adopted as the same as those in Section 4.4.2 but only for I i = 4 ,   0.6 , continuous ice sheet, and open water. Figure 19 shows the wave elevations at four different frequencies σ = 0.22 ,   2.86 ,   4.96 ,   6.34 . [ 5 , 0.5 ) presents the free surface on the left side of rigid body; the free surface between the rigid body and floating ice sheet lies ( 0.5 ,   5 ] . Surely, the [ 0.5 ,   0.5 ] represents the position of the rigid body.
In Figure 19a, the wave elevations on the left show the same values and decline, gradually approaching x = 0.5 . This is because the wavelength is much longer than the width a , completely penetrating through the rigid body which does not reflect the incident wave. For the right wave elevations, the results are less than the left ones. The wave elevations initially decrease, then slowly rise near the right side of the rigid body, and reach a maximum value at x = 5 . In addition, due to the reflection of waves by the ice sheet, the results in the case of an ice sheet and open water can be distinguished. For a greater σ , the wavelength becomes short, the waveform of the incident wave can be fully displayed, and some incident waves can be reflected by rigid body. As shown in Figure 19b, the waveform of the incident wave can be fully displayed. Due to the obstruction of the rigid body, the wave elevation near the left side of the rigid body is higher and the right elevation significantly decreases. As the frequency increases, the obstruction effect of the rigid body becomes stronger. It can be observed that the wave elevation in case of I i = 0.6 on the left side of the rigid body increases with frequency in Figure 19c,d, while the situation on the right side is the opposite. However, for the cases of an ice sheet and I i = 4 , the results seem to have no significant change. This is because a small-sized floating ice array has smaller flexural stiffness, resulting in fewer reflected waves than a large-sized one. This can be well observed from Figure 6.

5. Conclusions

A hybrid method was used to solve the problem of a floating body near the MIZ. The fluid domain is divided into three sub-domains. Around the floating body, a boundary integral equation is established, while in the sub-domain covered by ice sheets the velocity potential is expanded into an eigenfunction series, and in the far-field sub-domain with a free surface, a similar eigenfunction expansion is used to satisfy the radiation condition. The boundary integral equation and the coefficients of the eigenfunction expansions are solved together based on the continuous conditions of pressure and velocity on the interface between the sub-domains. The wave dispersion relation beneath the ice array and flexural stiffness and equivalent elastic modulus of the floating ice arrays are studied. In addition, the hydrodynamic impacts of floating ice arrays comprising different individual ice lengths on a floating square box are investigated. Furthermore, hydrodynamic forces resulting from variations in average ice thickness and the draft of floating ice arrays of different sizes are discussed. This study yields the following conclusions:
(1)
The flexural stiffness of floating ice arrays consisting of large-sized individual ice sheets decrease with the increase of wavelength, while for that of small-sized individual ice sheets, they gradually increase.
(2)
The equivalent elastic modulus of the floating array is minimally impacted by the wavelength of the incident wave. However, the change in the equivalent elastic modulus of the floating ice array mainly depends on the length of the individual ice sheet in the arrays.
(3)
Due to the reflection of waves by the floating ice arrays, the hydrodynamic forces of the floating body will oscillate, and the amplitude of the oscillation is affected by the length of the individual ice sheet, approaching the result of a continuous ice sheet as the wave frequency increases.
(4)
As the average thickness of ice increases, its stiffness also increases. It can weaken the ability of smaller waves to penetrate the ice surface. Then, it can induce that the hydrodynamic effect in the case of a thicker individual ice sheet for a large-sized or small-sized floating ice array is significant. Plus, at a high frequency, the changes in thickness are more sensitive to hydrodynamic effects.
(5)
Considering the different drafts of the floating ice arrays, the hydrodynamic forces of a square box in the case of a large-sized floating ice array with a draft are larger than that in the case of without a draft. However, for the small-sized floating ice array, the results in the case of a draft or no draft have no differences.
(6)
In the high-frequency stage, compared to large-sized floating ice arrays, small-sized floating ice arrays have a weaker ability to reflect waves, which further leads to an increase in the hydrodynamic coefficient of the rigid body.

Author Contributions

Conceptualization, Z.L. and Y.S.; methodology, B.W.; software, B.W.; validation, Z.L., Y.S. and B.W.; formal analysis, B.W.; investigation, B.W.; resources, B.W.; data curation, Y.S; writing—original draft preparation, B.W.; writing—review and editing, Z.L.; visualization, Y.S.; supervision, Z.L.; project administration, Y.S.; funding acquisition, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (52071162, 52101315, 12372242).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, A.K.; Holt, B.; Vachon, P.W. Wave propagation in the marginal ice zone: Model predictions and comparisons with buoy and synthetic aperture radar data. J. Geophys. Res. Oceans 1991, 96, 4605–4621. [Google Scholar] [CrossRef]
  2. Sakai, S.; Hanai, K. Empirical formula of dispersion relation of waves in sea ice. In Ice in the Environment: Proceedings of the 16th IAHR International Symposium on Ice, Dunedin, New Zealand, 2–6 December 2002; International Association of Hydraulic Engineering and Research: Singapore, 2002; Volume 2. [Google Scholar]
  3. Williams, T.D.; Squire, V.A. The effect of submergence on wave scattering across a transition between two floating flexible plates. Wave Motion 2008, 45, 361–379. [Google Scholar] [CrossRef]
  4. Vaughan, G.L.; Williams, T.D.; Squire, V.A. Perfect transmission and asymptotic solutions for reflection of ice-coupled waves by inhomogeneities. Wave Motion 2007, 44, 371–384. [Google Scholar] [CrossRef]
  5. Squire, V.A.; Williams, T.D. Wave propagation across sea-ice thickness changes. Ocean Model. 2008, 21, 1–11. [Google Scholar] [CrossRef]
  6. Lin, Q.; Lu, D.Q. Hydroelastic interaction between obliquely incident waves and a semi-infinite elastic plate on a two-layer fluid. Appl. Ocean Res. 2013, 43, 71–79. [Google Scholar] [CrossRef]
  7. Porter, R. Surface wave interaction with rigid plates lying on water. Wave Motion 2016, 66, 118–131. [Google Scholar] [CrossRef]
  8. Mosig, J.E.M.; Montiel, F.; Squire, V.A. A transport equation for flexural-gravity wave propagation under a sea ice cover of variable thickness. Wave Motion 2019, 88, 153–166. [Google Scholar] [CrossRef]
  9. Linton, C.M.; Chung, H. Reflection and transmission at the ocean/sea-ice boundary. Wave Motion 2003, 38, 43–52. [Google Scholar] [CrossRef]
  10. Meylan, M.H.; Masson, D. A linear Boltzmann equation to model wave scattering in the marginal ice zone. Ocean Model. 2006, 11, 417–427. [Google Scholar] [CrossRef]
  11. Smith, M.J.A.; Meylan, M.H. Wave scattering by an ice floe of variable thickness. Cold Reg. Sci. Technol. 2011, 67, 24–30. [Google Scholar] [CrossRef]
  12. Bennetts, L.G.; Squire, V.A. Model sensitivity analysis of scattering-induced attenuation of ice-coupled waves. Ocean Model. 2012, 45–46, 1–13. [Google Scholar] [CrossRef]
  13. Zhao, X.; Shen, H. A diffusion approximation for ocean wave scatterings by randomly distributed ice floes. Ocean Model. 2016, 107, 21–27. [Google Scholar] [CrossRef]
  14. Meylan, M.H.; Horvat, C.; Bitz, C.M.; Bennetts, L.G. A floe size dependent scattering model in two- and three-dimensions for wave attenuation by ice floes. Ocean Model. 2021, 161, 101779. [Google Scholar] [CrossRef]
  15. Yu, J.; Rogers, W.E.; Wang, D.W. A new method for parameterization of wave dissipation by sea ice. Cold Reg. Sci. Technol. 2022, 199, 103582. [Google Scholar] [CrossRef]
  16. Xu, B.; Guyenne, P. Nonlinear simulation of wave group attenuation due to scattering in broken floe fields. Ocean Model. 2023, 181, 102139. [Google Scholar] [CrossRef]
  17. Wang, R.; Shen, H.H. Experimental study on surface wave propagating through a grease–pancake ice mixture. Cold Reg. Sci. Technol. 2010, 61, 90–96. [Google Scholar] [CrossRef]
  18. McGovern, D.J.; Bai, W. Experimental study of wave-driven impact of sea ice floes on a circular cylinder. Cold Reg. Sci. Technol. 2014, 108, 36–48. [Google Scholar] [CrossRef]
  19. Zhao, X.; Shen, H.H. Wave propagation in frazil/pancake, pancake, and fragmented ice covers. Cold Reg. Sci. Technol. 2015, 113, 71–80. [Google Scholar] [CrossRef]
  20. Bennetts, L.G.; Alberello, A.; Meylan, M.H.; Cavaliere, C.; Babanin, A.V.; Toffoli, A. An idealised experimental model of ocean surface wave transmission by an ice floe. Ocean Model. 2015, 96, 85–92. [Google Scholar] [CrossRef]
  21. Nelli, F.; Bennetts, L.G.; Skene, D.M.; Monty, J.P.; Lee, J.H.; Meylan, M.H.; Toffoli, A. Reflection and transmission of regular water waves by a thin, floating plate. Wave Motion 2017, 70, 209–221. [Google Scholar] [CrossRef]
  22. Yiew, L.J.; Parra, S.M.; Wang, D.; Sree, D.K.K.; Babanin, A.V.; Law, A.W.K. Wave attenuation and dispersion due to floating ice covers. Appl. Ocean Res. 2019, 87, 256–263. [Google Scholar] [CrossRef]
  23. Cheng, S.; Tsarau, A.; Evers, K.U.; Shen, H. Floe Size Effect on Gravity Wave Propagation Through Ice Covers. J. Geophys. Res. Oceans 2019, 124, 320–334. [Google Scholar] [CrossRef]
  24. Parra, S.M.; Sree, D.K.K.; Wang, D.; Rogers, E.; Lee, J.H.; Collins, C.O.; Law, A.W.-K.; Babanin, A.V. Experimental study on surface wave modifications by different ice covers. Cold Reg. Sci. Technol. 2020, 174, 103042. [Google Scholar] [CrossRef]
  25. Ren, K.; Wu, G.X.; Thomas, G.A. Wave excited motion of a body floating on water confined between two semi-infinite ice sheets. Phys. Fluids 2016, 28, 127101. [Google Scholar] [CrossRef]
  26. Li, Z.F.; Shi, Y.Y.; Wu, G.X. Large amplitude motions of a submerged circular cylinder in water with an ice cover. Eur. J. Mech. B/Fluids 2017, 65, 141–159. [Google Scholar] [CrossRef]
  27. Ren, K.; Wu, G.X.; Ji, C.Y. Wave diffraction and radiation by a vertical circular cylinder standing in a three-dimensional polynya. J. Fluids Struct. 2018, 82, 287–307. [Google Scholar] [CrossRef]
  28. Li, Z.F.; Shi, Y.Y.; Wu, G.X. A hybrid method for linearized wave radiation and diffraction problem by a three dimensional floating structure in a polynya. J. Comput. Phys. 2020, 412, 109445. [Google Scholar] [CrossRef]
  29. Wehausen, J.V.; Laitone, E.V. Surface Waves. In Fluid Dynamics/Strömungsmechanik; Truesdell, C., Ed.; Springer: Berlin/Heidelberg, Germany, 1960; pp. 446–778. [Google Scholar] [CrossRef]
  30. Squire, V.A.; Dugan, J.P.; Wadhams, P.; Rottier, P.J.; Liu, A.K. Of Ocean Waves and Sea Ice. Annu. Rev. Fluid Mech. 1995, 27, 115–168. [Google Scholar] [CrossRef]
  31. Li, Z.F.; Shi, Y.Y.; Wu, G.X. Interaction of waves with a body floating on polynya between two semi-infinite ice sheets. J. Fluids Struct. 2018, 78, 86–108. [Google Scholar] [CrossRef]
  32. Newman, J.N. Marine Hydrodynamics; The MIT Press: Cambridge, MA, USA, 1977. [Google Scholar] [CrossRef]
  33. Li, Z.F.; Shi, Y.Y.; Wu, G.X. Interaction of wave with a body floating on a wide polynya. Phys. Fluids 2017, 29, 097104. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of floating ice array and floating body.
Figure 1. Schematic diagram of floating ice array and floating body.
Jmse 12 00272 g001
Figure 2. Wavenumber for individual ice sheets with different lengths.
Figure 2. Wavenumber for individual ice sheets with different lengths.
Jmse 12 00272 g002
Figure 3. The effect of individual ice sheet length on flexural stiffness. a   I i = 4 ; b   I i = 2 ; c   I i = 1 ; d   I i = 0.6 ; e   I i = 0.3 ; f   I i = 0.15 .
Figure 3. The effect of individual ice sheet length on flexural stiffness. a   I i = 4 ; b   I i = 2 ; c   I i = 1 ; d   I i = 0.6 ; e   I i = 0.3 ; f   I i = 0.15 .
Jmse 12 00272 g003aJmse 12 00272 g003b
Figure 4. The effect of wavelength on flexural stiffness. a   λ o w = 100 ; b   λ o w = 50 ; c   λ o w = 10 ; d   λ o w = 1 .
Figure 4. The effect of wavelength on flexural stiffness. a   λ o w = 100 ; b   λ o w = 50 ; c   λ o w = 10 ; d   λ o w = 1 .
Jmse 12 00272 g004
Figure 5. Reflection coefficient K R against σ .
Figure 5. Reflection coefficient K R against σ .
Jmse 12 00272 g005
Figure 6. Transmission coefficient K T against σ .
Figure 6. Transmission coefficient K T against σ .
Jmse 12 00272 g006
Figure 10. The added mass against σ .  ( a ) sway; ( b ) heave; c roll; d sway–roll.
Figure 10. The added mass against σ .  ( a ) sway; ( b ) heave; c roll; d sway–roll.
Jmse 12 00272 g010
Figure 11. The damping coefficient against σ . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Figure 11. The damping coefficient against σ . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Jmse 12 00272 g011aJmse 12 00272 g011b
Figure 12. The wave exciting force against σ . ( a ) sway; ( b ) heave; c roll.
Figure 12. The wave exciting force against σ . ( a ) sway; ( b ) heave; c roll.
Jmse 12 00272 g012
Figure 13. The added mass against σ in case of ice thickness h = 0.1 ,   0.15 . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Figure 13. The added mass against σ in case of ice thickness h = 0.1 ,   0.15 . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Jmse 12 00272 g013
Figure 14. The damping coefficient against σ in case of ice thickness h = 0.1 ,   0.15 . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Figure 14. The damping coefficient against σ in case of ice thickness h = 0.1 ,   0.15 . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Jmse 12 00272 g014
Figure 15. The wave exciting force against σ in case of ice thickness h = 0.1 ,   0.15 . ( a ) sway; ( b ) heave; c roll.
Figure 15. The wave exciting force against σ in case of ice thickness h = 0.1 ,   0.15 . ( a ) sway; ( b ) heave; c roll.
Jmse 12 00272 g015aJmse 12 00272 g015b
Figure 16. The added mass against σ in case of draft d = 0 ,   0.09 . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Figure 16. The added mass against σ in case of draft d = 0 ,   0.09 . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Jmse 12 00272 g016
Figure 17. The damping coefficient against σ in case of draft d = 0 ,   0.09 . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Figure 17. The damping coefficient against σ in case of draft d = 0 ,   0.09 . ( a ) sway; ( b ) heave; c roll; d sway–roll.
Jmse 12 00272 g017aJmse 12 00272 g017b
Figure 18. The exciting force against σ in case of draft d = 0 ,   0.09 . ( a ) sway; ( b ) heave; c roll.
Figure 18. The exciting force against σ in case of draft d = 0 ,   0.09 . ( a ) sway; ( b ) heave; c roll.
Jmse 12 00272 g018
Figure 19. The wave elevation distribution to the entire fluid field at different σ . a   σ = 0.22 ; b   σ = 2.86 ; c   σ = 4.94 ; d   σ = 6.34 .
Figure 19. The wave elevation distribution to the entire fluid field at different σ . a   σ = 0.22 ; b   σ = 2.86 ; c   σ = 4.94 ; d   σ = 6.34 .
Jmse 12 00272 g019aJmse 12 00272 g019b
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wan, B.; Shi, Y.; Li, Z. Two-Dimensional Wave Interaction with a Rigid Body Floating near the Marginal Ice Zone. J. Mar. Sci. Eng. 2024, 12, 272. https://doi.org/10.3390/jmse12020272

AMA Style

Wan B, Shi Y, Li Z. Two-Dimensional Wave Interaction with a Rigid Body Floating near the Marginal Ice Zone. Journal of Marine Science and Engineering. 2024; 12(2):272. https://doi.org/10.3390/jmse12020272

Chicago/Turabian Style

Wan, Bingbing, Yuyun Shi, and Zhifu Li. 2024. "Two-Dimensional Wave Interaction with a Rigid Body Floating near the Marginal Ice Zone" Journal of Marine Science and Engineering 12, no. 2: 272. https://doi.org/10.3390/jmse12020272

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop