Numerical Simulation of Wave Interaction with Payloads of Different Postures Using OpenFOAM

<jats:p>A three-dimensional numerical wave tank (NWT) is established with Open Source Field Operation and Manipulation (OpenFOAM) software and waves2foam to investigate wave interaction with payloads with different postures in the process of offshore lifting or lowering. Numerical results of regular wave interaction with a vertically suspending cylinder are presented first for validation by comparison with the published data. A series of simulation experiments are carried out, and the forces and the moments exerted by the regular waves on a fixed suspending cylinder payload and a fixed suspending cuboid payload with different postures are presented. It can be concluded from the results that the rotating rectangular payload (cuboid and cylinder) suffers a drastically changed moment when it is initially vertically placed, and the projection area of payload vertical to the force affects the corresponding force. The simulation results also show how the forces and the moments change with different posture angles. With some certain posture, the suspending payload suffers minimum forces and moments. Parametric study for the cuboid payload is done in the case of normal incidence. The influence of the payload’s size and wave parameters on forces and moments are analyzed. All of the numerical simulation results and conclusions provide the fundamentals for further research and safe control of offshore lifting or lowering.</jats:p>

section_content: With ever-increasing marine exploration and subsea resource exploitation, offshore cranes which are mounted on vessels and carry out lifting/lowering have been widely used in marine operations. While working on the sea, offshore cranes suffer from persistent disturbances induced by ocean waves. During lifting or lowering, the payloads may be subject to large hydrodynamic forces, which could cause payload damages or cable breaks. This would further cause accidents and impair the safety of life and property [1]. 

In order to lift/lower payloads on the sea safely and efficiently, the capability to estimate the hydrodynamic loads on payloads is of vital importance. The hydrodynamic loads on stationary structures in waves have been studied for the safe and cost-effective design of coastal and offshore structures in the past decades. Compared to physical experiments, which need to establish scaled models, numerical modeling is more practical. The numerical models based on potential flow theory and Navier-Stokes (N-S) equations are two main categories for the simulation of wave-structure interactions. 

The potential flow model is applied for wave interaction with large structures where viscous and turbulence effects can be ignored, such as the second-order potential flow theory model [2, 3] and the fully nonlinear potential flow theory model [4]. With the assumption that the flow is inviscid and flow irrotational, it is challenging for the potential flow theory to capture the nonlinear free surface correctly when wave breaking occurs. Computational Fluid Dynamics (CFD) based on Navier-Stokes (N-S) equations is used for highly nonlinear wave-structure interactions in the case of breaking wave impacts and evolution of vortices. Various methods or models have been considered for wave-structure interaction, such as the Institute of Environmental Hydraulics of Cantabria Field Operation and Manipulation (IHFOAM) model, which solves Volume-Averaged Reynolds-Averaged Navier-Stokes equations (VARANS) [5, 6], the multiple-layer σ-coordinate model [7], the Immersed Boundary Method [8], the Smooth Particle Hydrodynamics method [9, 10], and the Constrained Interpolation Profile method [11]. 

OpenFOAM, a free open-source C++ toolbox for the development of customized numerical solver (such as the naoe-FOAM-SJTU solver [12] ) based on CFD, has been applied in coastal and offshore engineering recently. Regular wave interaction with two tandem cylinders is studied with OpenFOAM [13], and an improved model named IHFOAM is used to study wave interaction with porous coastal structures [14, 15]. The performance of OpenFOAM for nonlinear wave interactions with offshore structures is assessed, with up to eighth order harmonics correctly modeled [16]. 

In addition to the normal incident wave interaction with structures, many researchers have also investigated the interaction of oblique waves with stationary structures, such as perforated caissons [17], bridge decks [18], and various other structures [19] [20] [21] [22] [23] [24]. The stationary nature of the structure makes it hard to rotate around different axes, the above oblique papers only focus on the situation of one single posture angle. Compared with stationary structures, the payloads can move with much more freedom while lifting or lowering payloads on the sea. Here, we want to reach a general conclusion when considering different posture angles, and to the authors' knowledge, there has been no previous research about the general postures' study of the payload. 

Importantly, the posture of the payload has an impact on the force and moment exerted by the wave; additionally, the force and moment can also change the posture. This paper focuses on studying the influence of different postures of the payloads on wave forces and moments exerted on the payloads; thus, we assume that the payload is fixed without linear motion and rotation. A cylinder payload and a cuboid payload, both fixed and suspended with different postures in regular waves, are investigated, respectively. By carrying out a series of simulations, the influence of the payloads' posture angles relative to the regular waves on the hydrodynamic forces and moments exerted on the payloads are analyzed. It can be concluded from the results that the rotating rectangular payload (cuboid or cylinder) suffers a drastically changed moment when it is initially vertically placed, and the direction of the moment is the same as axis' rotation except for one situation. The projection area of the payload vertical to the force affects the corresponding force. The analysis could provide help for developing control strategies for offshore cranes, such as choosing the appropriate payload posture during water entry, and then using a controller to keep the payload on a certain posture that suffers minimal forces or moments during water entry.
section_name: Introduction

section_name: Numerical Methods

section_content: In order to represent the payload's posture in the wave, two Cartesian frames are defined, as shown in Figure 1a Both air and water are assumed to be incompressible laminar fluid. The motion of the fluid continuum is described with the governing equations, i. e., the Navier-Stokes equations and the continuity equation [16],

where U is the fluid velocity, ρ is the fluid density, p is the fluid pressure, µ is the dynamic viscosity, t is the time, g is the gravity acceleration, and f σ is the surface tension. Only the laminar flow is considered in the study. 

,, θ φ ψ (pitch, roll, and yaw angles), which represent the pose relationship between the payload's body frame ( ) Both air and water are assumed to be incompressible laminar fluid. The motion of the fluid continuum is described with the governing equations, i. e., the Navier-Stokes equations and the continuity equation [16],

where U is the fluid velocity, ρ is the fluid density, p is the fluid pressure, μ is the dynamic viscosity,t is the time, g is the gravity acceleration, and f σ is the surface tension. Only the laminar flow is considered in the study.
section_name: Governing Equations

section_content: The Volume of Fluid (VOF) method is applied for tracking the free surface in OpenFOAM. In the VOF method, a phase function α is defined in each cell, which indicates the quantity of water in the cell. α is 1 if the cell is full of water, and it is 0 in empty cells. On the air-water interface, the value of α is between 0 and 1. The fluid density ρ and the dynamic viscosity µ in each cell are calculated with the equations,

where the subscripts 1 and 2 mean the values of water and air, respectively. The phase function α can be determined by solving an advection equation,

where the last term on the left-hand side is an artificial compression term and U α is the relative compression velocity [25].
section_name: Free Surface Tracking

section_content: The library waves2Foam is used to generate regular waves. The boundary condition and solve procedures are listed below.
section_name: Waves2Foam Library and WaveFoam Solver

section_content: The library waves2Foam is a toolbox for generating and absorbing water waves [26]. Waves are generated at the inlet and absorbed at the outlet. 

The velocities of regular waves are based on the linear Stokes' wave theory,

where u(x, z, t) is the horizontal velocity distribution, A is the wave amplitude, ω is the wave frequency, ϕ = kx -ωt and k is the wave number, h and is the water depth. 

The relaxation zone technique is used to absorb waves at the outlet. The relaxation function is

It is applied into the relaxation zone as follows,

where λ is either U or α. The variation of α R is the same as given in [27], and χ R represents a certain point in the relaxation zone. The definition of χ R is such that it is always 1 at the interface between the nonrelaxed part of the computational domain and the relaxation zone.
section_name: Waves2Foam Library

section_content: 1.
section_name: WaveFoam Solver and Boundary Conditions

section_content: The boundary name is just as shown in Figure 2. At the inlet of the 3D NWT, a specified boundary condition of fluid velocity U is set to waveVelocity, the boundary condition of the indicator phase function α is set to waveAlpha, and the boundary condition of the fluid pressure p is set to zeroGradient. At the top of the NWT, the velocity U is set to pressureInletOutletVelocity, which is a default boundary condition in OpenFOAM, the pressure p is set to totalPressure, and the phase function α is set to inletOutlet. 

For the remaining parts of the NWT and the fixed suspending objects, the boundary conditions are considered as solid walls, where the fluid velocity U is set to a fixed value of zero, the fluid pressure p and the indicator phase function α are set to zeroGradient.
section_name: Boundary conditions

section_content: Solving procedure

The waveFoam solver starts with the preprocessor, which is used to set up wave properties and computational meshes. The meshes of the NWT are generated by using the built-in tool blockMesh and snappyHexMesh. The N-S equations are discretized into a set of algebraic equations by integrating the boundary conditions over the whole solution domain and time domain. The physical parameters of the whole domain like the fluid pressure p and the fluid velocity U, etc., are calculated and updated at each timestep by calling solver waveFoam. 

where λ is either U or α. The variation of R α is the same as given in [27], and R χ represents a certain point in the relaxation zone. The definition of R χ is such that it is always 1 at the interface between the nonrelaxed part of the computational domain and the relaxation zone.
section_name: 2.

section_content: The boundary name is just as shown in Figure 2. At the inlet of the 3D NWT, a specified boundary condition of fluid velocity U is set to waveVelocity, the boundary condition of the indicator phase function α is set to waveAlpha, and the boundary condition of the fluid pressure p is set to zeroGradient. At the top of the NWT, the velocity U is set to pressureInletOutletVelocity, which is a default boundary condition in OpenFOAM, the pressure p is set to totalPressure, and the phase function α is set to inletOutlet. 

For the remaining parts of the NWT and the fixed suspending objects, the boundary conditions are considered as solid walls, where the fluid velocity U is set to a fixed value of zero, the fluid pressure p and the indicator phase function α are set to zeroGradient.
section_name: WaveFoam Solver and Boundary Conditions 1. Boundary conditions

section_content: The waveFoam solver starts with the preprocessor, which is used to set up wave properties and computational meshes. The meshes of the NWT are generated by using the built-in tool blockMesh and snappyHexMesh. The N-S equations are discretized into a set of algebraic equations by integrating the boundary conditions over the whole solution domain and time domain. The physical parameters of the whole domain like the fluid pressure p and the fluid velocity U, etc., are calculated and updated at each timestep by calling solver waveFoam.
section_name: Solving procedure

section_content: A 3D numerical wave tank (NWT) is established with the above numerical methods of OpenFOAM and waves2foam. To validate the 3D NWT model, we compare the numerical results of wave interaction with a fixed and vertically suspended cylinder payload with the published data [16].
section_name: Comparison Against Published Data

section_content: A 3D numerical wave tank (NWT) is established, as shown in Figure 3. Its geometry has the outer dimensions 15 and the relaxation zone of 1. 5L, where L is the wavelength. A cylinder whose radius 0. 125 a m = is stationary and vertical is suspended in the tank, leaving a 1mm gap beneath to the bed of the tank. The length of the cylinder is1m. The cylinder is located at 7. 5m from the paddles in the center of the tank. A wave gauge WG2 is placed 2mm in front of the upstream stagnation point of the cylinder to monitor the wave field around the cylinder, and a wave gauge WG1 is placed 0. 77m from the inlet to monitor the wave elevation.
section_name: Numerical Wave Tank

section_content: A 3D numerical wave tank (NWT) is established, as shown in Figure 3. Its geometry has the outer dimensions 15 m × 4 m × 1. 2 m with the water depth h = 0. 505 m and the relaxation zone of 1. 5 L, where L is the wavelength. A cylinder whose radius a = 0. 125 m is stationary and vertical is suspended in the tank, leaving a 1 mm gap beneath to the bed of the tank. The length of the cylinder is 1 m. The cylinder is located at 7. 5 m from the paddles in the center of the tank. A wave gauge WG2 is placed 2 mm in front of the upstream stagnation point of the cylinder to monitor the wave field around the cylinder, and a wave gauge WG1 is placed 0. 77 m from the inlet to monitor the wave elevation. Two regular wave cases [16] R1 and R2 are reproduced with our 3D NWT. The wave parameters are shown in Table 1, where h is the water depth, k is the wavenumber, A is the wave amplitude, and T is the wave period. The mesh resolution in the computational domain affects the numerical solution. The built-in mesh generator blockMesh is used to generate meshes of hexahedral cells, then snappyHexMesh in OpenFOAM is used to generate the cylinder. The mesh consists of multilevel grids, as shown in Figure 4. In the areas around the payload, the grid cells have a resolution of x Δ in the horizontal direction and z Δ in the vertical direction, which are measured by the cells per wavelength and wave height. 

Three different time steps are used here for the convergence examination. As shown in Figure 5, three cases are set to a fixed time-step, and the results are convergent. For each time-step, one, two, and three inner iterations (nOuterCorrectors in OpenFOAM) are used for convergence examination, the result is the same as Figure 5. For each inner iteration, the PIMPLE algorithm is called three times (nCorrectors in OpenFOAM). For the remaining cases in the paper, one nOuterCorrectors and three nCorrectors are used, the simulation time is 18 s, and the fixed time step is set to 0. 005 s, the courant numbers during the simulation are all less than 0. 1. 

The time history of horizontal force x F on the cylinder payload with three different mesh schemes for the regular waves are shown in Figure 6. From this grid convergence examination, it can be seen that the results of Mesh 2 and Mesh 3 are convergent and Mesh 2 uses much less time; thus, the intermediate Mesh 2 is selected in this paper. For the Mesh 2 scheme, multilevel grids are used just as Figure 4 shows: in total, 470 cells in the x-direction, 125 cells in the y-direction, and 100 cells in the z-direction. The mesh around the inlet, outlet, and object is dense, and the rest transitions smoothly. The mesh around the cylinder and free surface is uniform: 110 cells per wavelength and 110 cells in total are set in the x-direction, 30 cells per wave height, and 60 cells in total are set in the z-direction and 60 cells in total in 1 m are set in the y-direction. 

The simulations are run on purchased Dell T7920 workstation with Intel Xeon (R)E5 2699v4 CPU, 128GB RAM, and 44 cores. The comparison of the computation cost, the total cell numbers, the number of cores, and the simulated time under three different mesh schemes are illustrated in Table 2. Two regular wave cases [16] R1 and R2 are reproduced with our 3D NWT. The wave parameters are shown in Table 1, where h is the water depth, k is the wavenumber, A is the wave amplitude, and T is the wave period. The mesh resolution in the computational domain affects the numerical solution. The built-in mesh generator blockMesh is used to generate meshes of hexahedral cells, then snappyHexMesh in OpenFOAM is used to generate the cylinder. The mesh consists of multilevel grids, as shown in Figure 4. In the areas around the payload, the grid cells have a resolution of ∆x in the horizontal direction and ∆z in the vertical direction, which are measured by the cells per wavelength and wave height. 

Three different time steps are used here for the convergence examination. As shown in Figure 5, three cases are set to a fixed time-step, and the results are convergent. For each time-step, one, two, and three inner iterations (nOuterCorrectors in OpenFOAM) are used for convergence examination, the result is the same as Figure 5. For each inner iteration, the PIMPLE algorithm is called three times (nCorrectors in OpenFOAM). For the remaining cases in the paper, one nOuterCorrectors and three nCorrectors are used, the simulation time is 18 s, and the fixed time step is set to 0. 005 s, the courant numbers during the simulation are all less than 0. 1. 

The time history of horizontal force F x on the cylinder payload with three different mesh schemes for the regular waves are shown in Figure 6. From this grid convergence examination, it can be seen that the results of Mesh 2 and Mesh 3 are convergent and Mesh 2 uses much less time; thus, the intermediate Mesh 2 is selected in this paper. For the Mesh 2 scheme, multilevel grids are used just as Figure 4 shows: in total, 470 cells in the x-direction, 125 cells in the y-direction, and 100 cells in the z-direction. The mesh around the inlet, outlet, and object is dense, and the rest transitions smoothly. The mesh around the cylinder and free surface is uniform: 110 cells per wavelength and 110 cells in total are set in the x-direction, 30 cells per wave height, and 60 cells in total are set in the z-direction and 60 cells in total in 1 m are set in the y-direction. 

The simulations are run on purchased Dell T7920 workstation with Intel Xeon (R)E5 2699v4 CPU, 128GB RAM, and 44 cores. The comparison of the computation cost, the total cell numbers, the number of cores, and the simulated time under three different mesh schemes are illustrated in Table 2. L is the wavelength. H is the wave height. L is the wavelength. H is the wave height.
section_name: Numerical Wave Tank

section_content: Before the simulation, the surface elevation at WG1 for the wave R1 is compared with theory results, the result is as Figure 7 shows, the surface elevation agrees well with the theory. The free surface elevation and horizontal force are compared with published data. The time histories of the free surface elevation at WG2, and the corresponding amplitude spectra obtained by applying the FFT algorithm to the time histories are shown in Figure 8. The surface elevation is normalized by the wave amplitude A, and the time is normalized by the wave period t. The time series of the horizontal force on the cylinder and the corresponding amplitude spectra are presented in Figure 9. The force is normalized by 0. 5 gAS ρ, where ρ is density of the water, and S is the cross-sectional area of the payload in the water perpendicular to the wave propagation direction. It can be seen that the results obtained with our NWT model match with the published data [16]. It is validated that our present 3D NWT numerical model can be used to calculate the wave load exerted on the payload with a reasonable degree of accuracy.
section_name: Comparison with the Published Data

section_content: Before the simulation, the surface elevation at WG1 for the wave R1 is compared with theory results, the result is as Figure 7 shows, the surface elevation agrees well with the theory. The free surface elevation and horizontal force are compared with published data. The time histories of the free surface elevation at WG2, and the corresponding amplitude spectra obtained by applying the FFT algorithm to the time histories are shown in Figure 8. The surface elevation is normalized by the wave amplitude A, and the time is normalized by the wave period t. The time series of the horizontal force on the cylinder and the corresponding amplitude spectra are presented in Figure 9. The force is normalized by 0. 5ρgAS, where ρ is density of the water, and S is the cross-sectional area of the payload in the water perpendicular to the wave propagation direction. It can be seen that the results obtained with our NWT model match with the published data [16]. It is validated that our present 3D NWT numerical model can be used to calculate the wave load exerted on the payload with a reasonable degree of accuracy.
section_name: Comparison with the Published Data

section_content: The 3D NWT established in Section 3 is applied in a series of simulation experiments in this section. Our study focuses on the influence of different postures of the payloads on wave forces and moments exerted on the payloads while suspending in the sea. We assume that the payloads are fixed and stationary while suspended in the sea without considering their translational and rotational motions caused by wave forces and rotational moments. A cylinder payload and a cuboid payload with different postures that are fixed and suspended in the regular wave R1 are simulated respectively.,,
section_name: Numerical Results

section_content: The 3D NWT established in Section 3 is applied in a series of simulation experiments in this section. Our study focuses on the influence of different postures of the payloads on wave forces and moments exerted on the payloads while suspending in the sea. We assume that the payloads are fixed and stationary while suspended in the sea without considering their translational and rotational motions caused by wave forces and rotational moments. A cylinder payload and a cuboid payload with different postures that are fixed and suspended in the regular wave R1 are simulated respectively. The same cylinder payload in Section 3 is used here. In this subsection, in addition to the vertical suspension, several postures of the cylinder in the NWT are considered. The posture of the cylinder payload in the 3D NWT is shown in Figure 10. The wave condition is the same as the regular wave R1. The same cylinder payload in Section 3 is used here. In this subsection, in addition to the vertical suspension, several postures of the cylinder in the NWT are considered. The posture of the cylinder payload in the 3D NWT is shown in Figure 10. The wave condition is the same as the regular wave R1. To validate the 3D NWT, the horizontal wave force F x is compared with results computed by Morison's equation [27]

where C m is the added mass coefficient (C m = 1. 15), C d is the drag coefficient (C d = 1), V is the volume of the payload in the water, S is the cross-sectional area of the payload in the water perpendicular to the wave propagation direction, and. 

u is the horizontal acceleration. 

As shown in Figure 11, the normalized first-harmonic forces and moments are obtained by applying the FFT algorithm to the time histories. The first-harmonic forces are normalized by 0. 5ρgAS and the first-harmonic rotational moments are normalized by ρgdAS where d is the draft of the cylinder. In the latter study, the same normalization method is used for the forces and moments. 

where m C is the added mass coefficient ( =1. 15

), V is the volume of the payload in the water, S is the cross-sectional area of the payload in the water perpendicular to the wave propagation direction, and u  is the horizontal acceleration. 

As shown in Figure 11, the normalized first-harmonic forces and moments are obtained by applying the FFT algorithm to the time histories. The first-harmonic forces are normalized by 0. 5 gAS ρ and the first-harmonic rotational moments are normalized by gdAS ρ where d is the draft of the cylinder. In the latter study, the same normalization method is used for the forces and moments. It is obvious that y F, the lateral force, and, y z M M the rotational moments about x-axis and zaxis, are much less than the others, and can be neglected regardless of pitch angle θ. The numerical results match the physical phenomena and can be explained easily with the force analysis. In the case In Figure 11a, it can be seen that the horizontal force F x obtained with the 3D NWT matches with that computed by Morison's equation. The horizontal force F x decreases with the pitch angle and F z increases with the pitch angle. This could be explained by the decrease of the projection area of the cylinder on the surface yoz and increase of the projection area of the cylinder on the surface xoy. 

It is obvious that F y, the lateral force, and M y, M z the rotational moments about x-axis and z-axis, are much less than the others, and can be neglected regardless of pitch angle θ. The numerical results match the physical phenomena and can be explained easily with the force analysis. In the case of ψ = 0 •, φ = 0 • the rotational moment M y exerted on the cylinder about y-axis depends on the horizontal force F x and the vertical force F z. M y increases with the increase of θ from 0 • to 60 •, then it decreases with the increase of θ from 60 • to 90 •. The maximum moment with θ = 60 • is 10 times larger than that of vertical suspension.
section_name: Numerical Results

section_content: Roll angle φ = (0

In Figure 11b, it can be seen that the horizontal force F x obtained from the 3D NWT matches with that computed with Morison's equation only for a limited range near to φ = 0 •. The reason is that the two coefficients C m, C d change with the increase of the roll angle φ. The values of C m, C d at φ = 0 • no longer work with the increase of φ. 

It can be seen that the lateral force F y, the rotational moments M y, M z, about the y-axis and z-axis are not zero but small values. Both the horizontal force F x and the vertical force F z increase with the roll angle φ, but F x decreases when the angle is 75 •. The projection area of the cylinder on the surface xoy increases with the increase of θ from 0 • to 90 •, and the vertical force exerted on the cylinder also increases. The increase is quicker with the roll angel φ from 60 • to 90 •. The rotational moment M x increases with the increase of φ from 0 • to 60 •, and then decreases with the increase of φ from 60 • to 90 •. 

From the above simulations of the two cases, we can see that the changes of force could be explained by the change of the corresponding projection area. Additionally, the moment around a certain axis changes drastically with the change of angle. For example, when the cylinder rotates around the y-axis (pitch angle), the max value of M y is 10 times larger than the initial value. For the vertical cylinder payload, there is no angle where all the forces and moments are minimal, but the initial posture could be an optimal selection.
section_name: 2.

section_content: A cuboid payload is fixed and suspended in the 3D NWT as shown in Figure 1. The size of the cuboid is 1m × 0. 5m × 0. 5m, and the draft is d = 0. 25m. The cuboid's posture in the 3D NWT is represented by the three Euler angles. A series of simulations are done with different postures of the suspending cuboid in the 3D NWT.
section_name: Case 2: A Cuboid Payload Fixed Suspending in the 3D NWT

section_content: Yaw angle ψ = (0 •, 15

The normalized forces and moments on the cuboid payload versus its yaw angle ψ are shown in Figure 12a. With the normal incident regular waves, the horizontal force F x decreases with the increase of ψ from 0 • to 90 •. The projection area on the surface yoz is the biggest when ψ = 0 •. The projection area on the surface yoz decreases with the increase of ψ from 0 • to 90 •. It can be seen only F x changes drastically with the angle, which is not the same as the results of cylinder where the moment changes drastically.
section_name: 1.

section_content: Pitch angle θ = (0 •, 15

The normalized forces and moments on the payload versus the pitch angle θ are shown in Figure 12b. It is obvious that both the lateral force F y and the rotational moment M x and M z are near to zero no matter the pitch angle. The horizontal force F x, the vertical force F z. and the rotational moment M y are all symmetrical around 45 •. This could be easily explained by the change in the projection area. When the pitch angle θ increases from 0 • to 45 •, the horizontal force F x decreases. Additionally, it increases when θ increases from 45 • to 90 •. At θ = 0 •, 90 •, the regular waves which are normally incident to the cuboid's face with the largest surface area exert the maximum horizontal force on the cuboid. It can be seen that there is also not a drastically changed moment with the change of the angle.
section_name: 2.

section_content: Roll angle φ = (0

The normalized forces and moments on the cuboid payload versus the roll angle φ are shown in Figure 12c. It can be seen that compared with M x, the changes of other forces and moments are small, and M x increases very quickly with the increase of φ from 0 • to 30 •, and decreases with φ from 30 • to 90 •. The phenomenon of a drastically changed moment is similar to the results of the cylinder. For the two cases of the cylinder, the moment changes drastically with the angle. For the three cases of the cuboid payload, when the cuboid payload rotates around the x-axis (roll angle), the change of force and moment is similar to the cylinder cases. The changes of force could be explained by the change of the corresponding projection area, and the moment around a certain axis changes drastically with the change of angle. However, results when rotating around the z-axis (yaw angle) and y-axis (pitch angle) show no drastically changed moment. 

To show that the difference could be brought by the initial posture, we plot the normalized forces and moments on the cylinder payload versus its yaw angle ψ when the roll angle φ = 90 •, just as Figure 13 shows. It can be seen that the result is similar to Figure 12a, b. When the roll angle φ = 90 •, the cylinder is horizontally placed, its length side along the y-axis. When this happens, the phenomenon of a drastically changed moment disappears. 

For further study, we also plot the result when the vertically placed cuboid payload rotates around the y-axis and z-axis. Just as Figure 14 shows, the result is similar to Figure 12c. When the cuboid is vertically placed and the roll angle φ = 90 •, its long side along the z-axis. However, in Figure 14b, there is an exception, the drastically changed moment is not around the z-axis but the y-axis, other results are all as expected. 

All three above figures show that the drastically changed moment is brought about by the initial posture. The moment around a certain axis changes drastically with the change of angle when the payload is vertically placed (which means the long side of the payload is vertical to the water surface) such as in Figure 14, and this phenomenon could happen when the horizontally placed payload changes to the vertical posture, such as in Figure 12c. For the two cases of the cylinder, the moment changes drastically with the angle. For the three cases of the cuboid payload, when the cuboid payload rotates around the x-axis (roll angle), the change of force and moment is similar to the cylinder cases. The changes of force could be explained by the change of the corresponding projection area, and the moment around a certain axis changes drastically with the change of angle. However, results when rotating around the z-axis (yaw angle) and y-axis (pitch angle) show no drastically changed moment. 

To show that the difference could be brought by the initial posture, we plot the normalized forces and moments on the cylinder payload versus its yaw angle ψ when the roll angle 0 =90 φ, just as Figure 13 shows. It can be seen that the result is similar to Figure 12a, b. When the roll angle

the cylinder is horizontally placed, its length side along the y-axis. When this happens, the phenomenon of a drastically changed moment disappears. 

For further study, we also plot the result when the vertically placed cuboid payload rotates around the y-axis and z-axis. Just as Figure 14 shows, the result is similar to Figure 12c. When the cuboid is vertically placed and the roll angle

, its long side along the z-axis. However, in Figure 14b, there is an exception, the drastically changed moment is not around the z-axis but the y-axis, other results are all as expected. All three above figures show that the drastically changed moment is brought about by the initial posture. The moment around a certain axis changes drastically with the change of angle when the payload is vertically placed (which means the long side of the payload is vertical to the water surface) such as in Figure 14, and this phenomenon could happen when the horizontally placed payload changes to the vertical posture, such as in Figure 12c.
section_name: 3.

section_content: Yaw ψ and roll φ concurrently change from 0 • to 90 •, pitch θ = 0 •. 

The force and moment exerted on the cuboid versus the yaw angle ψ and the roll angle φ is shown in Figures 15 and 16. For the horizontal force F x, it decreases with the yaw angle ψ regardless of the roll angle φ. When the roll angle is 0 •, the cuboid is horizontally placed, with its long side vertical to the wave direction when the yaw angle is 0 •. When the roll angle increases from 0 • to 90 •, the long side gradually changes to the vertically placed position; thus, the gradient along the yaw angle decreases with the increase of the roll angle. For the lateral force F y, the result is symmetrical about yaw angle and roll angle. 

The results of M x could also be explained by the conclusion raised above. When the yaw angle is 0 •, the roll angle increases to 90 •, the cuboid changes from horizontally placed to be vertically placed, and the M x changes drastically with the roll angle. When the yaw angle is 90 •, the cuboid payload could not change to be vertically placed with the change of roll angle, and the phenomenon of a drastically changed moment disappears. When the yaw angle change from 0 • to 90 •, the phenomenon gradually disappears. 

The result of M y and M z demonstrate the exceptional condition in Figure 14b. When the roll angle is 0 •, the cuboid is horizontally placed, and there is no drastically changed moment. When the roll angle is 90 •, the cuboid is vertically placed and the moment M y, instead of M z, drastically changes with changes with the yaw angle. There is also a transition when the roll angle increases from 0 0 to 0 90

. The amplitude of z M is much less than the others, and can be neglected.
section_name: 4.

section_content: The parameter studies are done to analyze the effects of the cuboid's size and wave parameters on the forces and moments exerted on the cuboid payload. Here, in order to focus on the effects of cuboid's size and wave parameters, no posture angles are considered and the cuboid is horizontally placed.
section_name: Parameter Studies

section_content: Cuboid's size effects on forces and moments

The normalized forces and moments on the cuboid payload with different drafts, lengths, and widths are shown in Figure 17a-c
section_name: 1.

section_content: The normalized forces and moments on the cuboid payload with different drafts, lengths, and widths are shown in Figure 18a, b. It can be seen that both the horizontal force F x and the vertical force F z increase with the increase of the wave amplitude and wavelength. The change of other forces and moments could be neglected.
section_name: Wave's parameters effects on forces and moments

section_content: In order to investigate regular wave interaction with a fixed suspending payload with different postures, a three-dimension NWT based on OpenFOAM and waves2foam is established. Regular wave interaction with a vertically suspended cylinder is investigated. The free surface elevation, horizontal wave force, as well as the corresponding amplitude spectra obtained by the FFT algorithm, are compared with the theory result and the results reported in [16] for validation. Then, the representation of the payload's posture in the regular wave is given. The forces and moments exerted on a suspended cylinder and a suspended cuboid with different postures are investigated separately. Finally, parameter studies in the case of payload's size wave parameters are considered. 

It can be concluded that the moment around a certain axis changes drastically with the change of the same angle when the payload is initially vertically placed (which means the long side of the payload is vertical to the water surface). For example, the moment around the y-axis could change drastically when rotating around the y-axis. This phenomenon could also happen when the horizontally placed payload (which means the long side parallel to the sea level) changes to the vertical posture. There is an exception: when rotated around the z-axis, the drastically changed moment is not around the z-axis but the y-axis. Therefore, for the rectangular shape payload, it is better to keep the payload horizontally placed to prevent the drastic change of the moment. Additionally, the projection area of the payload vertical to the direction of force affects the corresponding force. It is better to keep the short side vertical to the incident direction of the wave; thus, a minimal horizontal force can be obtained. Through the simulations, some certain posture of the payload with the minimum forces and moments can be reached. It can guide the design of control strategies for the safe operation of offshore cranes, such as keeping the payload to a certain posture that suffers minimal force and moment or changing the controller weight of some forces and moments under specific circumstances. From the above parameter simulations, we can see that the horizontal force F x and the rotational moment M y exerted on the cuboid payload increase with the increase of its draft d, and its length l. The vertical force F z and the rotational moment M x increase with its width B. The increase of the wave amplitude A and wave length L cause the increase of the horizontal force F x and the vertical force F z.
section_name: Conclusions

section_content: In order to investigate regular wave interaction with a fixed suspending payload with different postures, a three-dimension NWT based on OpenFOAM and waves2foam is established. Regular wave interaction with a vertically suspended cylinder is investigated. The free surface elevation, horizontal wave force, as well as the corresponding amplitude spectra obtained by the FFT algorithm, are compared with the theory result and the results reported in [16] for validation. Then, the representation of the payload's posture in the regular wave is given. The forces and moments exerted on a suspended cylinder and a suspended cuboid with different postures are investigated separately. Finally, parameter studies in the case of payload's size wave parameters are considered. 

It can be concluded that the moment around a certain axis changes drastically with the change of the same angle when the payload is initially vertically placed (which means the long side of the payload is vertical to the water surface). For example, the moment around the y-axis could change drastically when rotating around the y-axis. This phenomenon could also happen when the horizontally placed payload (which means the long side parallel to the sea level) changes to the vertical posture. There is an exception: when rotated around the z-axis, the drastically changed moment is not around the z-axis but the y-axis. Therefore, for the rectangular shape payload, it is better to keep the payload horizontally placed to prevent the drastic change of the moment. Additionally, the projection area of the payload vertical to the direction of force affects the corresponding force. It is better to keep the short side vertical to the incident direction of the wave; thus, a minimal horizontal force can be obtained. Through the simulations, some certain posture of the payload with the minimum forces and moments can be reached. It can guide the design of control strategies for the safe operation of offshore cranes, such as keeping the payload to a certain posture that suffers minimal force and moment or changing the controller weight of some forces and moments under specific circumstances.
section_name: Conclusions

section_content: Funding: This research was funded by the Joint Fund of the National Nature Science Foundation of China and Shandong Province, grant number No. U1706228.

section_name: Conflicts of Interest:

section_content: The authors declare no conflicts of interest.
section_name: Conflicts of Interest: