The motion responses of the hybrid energy platform were calculated using AQWA, a commercial hydrodynamic software developed by ANSYS, based on the boundary element method (BEM) and linear potential flow theory (
ANSYS, 2021).
Fig. 3 shows the overall numerical analysis procedure. First, frequency-domain analysis was performed using the AQWA-LINE module to calculate the hydrodynamic coefficients that account for the interactions between the FOWT and the WECs under incident wave frequencies. These coefficients were used as input data for the time-domain analysis module AQWA-NAUT, which calculates the external forces acting on the hybrid platform. In AQWA-NAUT, hinge constraints were applied to connect the FOWT and four WECs. Time-domain simulations were conducted considering the nonlinear Froude–Krylov forces and nonlinear restoring forces resulting from changes in the wetted surface of the floating bodies.
2.2.1 Mathematical formulation
The three-dimensional frequency-domain commercial software AQWA-LINE was used to obtain the hydrodynamic coefficients of the hybrid energy platform. This software is a numerical hydrodynamic analysis program based on the BEM and linear potential flow theory. The flow can be described by a velocity potential assuming that the fluid field in the computational domain is inviscid, irrotational, and incompressible, and the time-domain velocity potential is expressed as
Eq. (1). Under these assumptions, the governing equation in the computational domain becomes the Laplace equation, as shown in
Eq. (2):
where
X→(X,Y,Z) denotes the position vector on the body surface;
t is the time;
aw is the incident wave amplitude;
φ is the complex velocity potential; and
ω is the wave frequency. The complex velocity potential is composed of the sum of the incident, diffraction, and radiation velocity potentials, as expressed in
Eq. (3), based on the principle of linear superposition:
where
φ1 and
φd are the first-order incident and diffraction velocity potentials, respectively.
φrj denotes the radiation velocity potential caused by unit amplitude motion in the
jth direction.
φ1 can be obtained using the theoretical expression shown in
Eq. (4).
where
g represents gravitational acceleration;
k is the wavenumber;
h is the water depth; and
α is the incident wave angle. To calculate
φd and
φr Green’s theorem and Green’s function were applied to
Eq. (2), resulting in the following boundary integral equation:
where S represents the wetted surface of the body, and n→ is the normal vector on the surface. G is the Green’s function, which is expressed as follows:
where ξ(ξ, η, ζ) represents the source position in the source distribution method, and J0 is the Bessel function of the first kind:
φd and
φr can be obtained by solving the diffraction and radiation problems, respectively, by substituting the following boundary conditions into
Eq. (5).
where v→ is the velocity vector of the floating body.
The Froude–Krylov force and diffraction force acting on the body due to the incident wave are calculated using φ1 and φd. In addition, φr is used to calculate the radiation force caused by the motion of the body, which is represented by the added mass and radiation damping coefficients. Using these, the equation of motion of the body is formulated as follows:
where M is the mass matrix of the hybrid platform. Ma and C are the added mass matrix and radiation damping matrix, respectively, which consider multi-body interactions. Khys represents the stiffness matrix of the platform, and Km denotes the stiffness matrix of the mooring lines. F is the external force vector acting on the platform, consisting of the sum of the Froude–Krylov force, diffraction force, and PTO system force.
The time-domain numerical analysis of the hybrid energy platform was performed using AQWA-NAUT. The equations of motion of the platform are expressed in the time domain using Cummins’ equation, as shown in
Eq. (10).
where
R denotes the velocity impulse function matrix;
F(
t) represents the external force. When
Eq. (10) is defined as the equations of motion for the FOWT and the WEC, respectively, they can be expressed as
Eqs. (11) and
(12).
where
Fm and
Km represent the force and stiffness exerted by the mooring lines, respectively, and
FEX is the wave excitation force, which is the sum of the Froude–Krylov force and diffraction force acting on the floating body.
cPTO is the PTO damping coefficient acting on the WEC, and
FFOWT,con and
FWEC,con denote the constraint forces acting on the FOWT and WEC because of the hinge constraints. Detailed mathematical formulations can be found in the AQWA theory manual (
ANSYS, 2021).