A floating Oscillating Water Column (OWC) wave energy converter, a Backward Bent Duct Buoy (BBDB), was simulated using a stateoftheart, twodimensional, fullynonlinear Numerical Wave Tank (NWT) technique. The hydrodynamic performance of the floating OWC device was evaluated in the time domain. The acceleration potential method, with a fullupdated kernel matrix calculation associated with a mode decomposition scheme, was implemented to obtain accurate estimates of the hydrodynamic force and displacement of a freely floating BBDB. The developed NWT was based on the potential theory and the boundary element method with constant panels on the boundaries. The mixed EulerianLagrangian (MEL) approach was employed to capture the nonlinear free surfaces inside the chamber that interacted with a pneumatic pressure, induced by the timevarying airflow velocity at the air duct. A special viscous damping was applied to the chamber free surface to represent the viscous energy loss due to the BBDB’s shape and motions. The viscous damping coefficient was properly selected using a comparison of the experimental data. The calculated surface elevation, inside and outside the chamber, with a tuned viscous damping correlated reasonably well with the experimental data for various incident wave conditions. The conservation of the total wave energy in the computational domain was confirmed over the entire range of wave frequencies.
INTRODUCTION
Since the 19
^{th}
century, ocean wave energy has often been studied as a key renewable energy resource. Over 1000 patented wave energy takeoff devices have been proposed as energy converters in the last several decades. Recently, as a specific functional structure system, an array of porous circular cylinders has been studied for use as an efficient wave energy takeoff or attenuation system (Park et al., 2010). In general, wave energy takeoff techniques are based on nine ideas (McCormick, 2007). One of the most practical and energy efficient concepts is the Oscillating Water Column (OWC), which is based on a pneumatic power takeoff inside the chamber achieved by using specially designed air turbines, such as the Wells, Impulse, and DennissAuld turbines.
Since Masuda (1971) first proposed a commerciallyavailable OWC device, several commerciallevel fixedtype OWC plants have recently been constructed and operated successfully (Heath et al., 2000). Numerical analyses of fixed OWC systems have been performed by many researchers (Brendmo et al., 1996; Wang et al., 2002; Delaure and Lewis, 2003). However, most of these analyses were founded on linearbased numerical models. Koo and Kim (2010) recently developed a fully nonlinear timedomain model of a landbased OWC system. Their simulation included viscous energy loss from the chamber skirt and pneumatic pressure from oscillatory airflows in the chamber.
For a floatingtype OWC system, Masuda et al. (1987) proposed a special type of OWC, a Backward Bent Duct Buoy (BBDB), which is thought to be one of the most energyefficient OWC devices. The BBDB has the typical characteristics of dynamic behavior due to its unique body shape, such as a reverse timemean drift force. Therefore, many studies on the hydrodynamic behaviors of BBDBs, including reverse drifting, have been conducted, either numerically or experimentally (McCormick and Sheehan, 1992; Hong et al. 2004a; 2004b; Kim et al., 2006; 2007; Nagata et al., 2008; 2009; Imai et al., 2009; Toyota et al., 2008; 2009). Suzuki et al. (2011) performed a numerical investigation to determine the optimal twodimensional (2D) hull design of the BBDB using the eigenfunction expansion method, with the relevant experiment being conducted by Toyota et al. (2010). Using a 2D numerical wave tank (NWT) technique, Koo and Lee (2011) and Koo et al. (2012) calculated the hydrodynamic behavior of a BBDB and chamber free surfaces with different shapedcorners. A proper viscous damping coefficient was applied to their numerical model, which was deducted from the experimental results.
In this study, a stateoftheart 2D fullynonlinear NWT technique was fully described for a simulation of a floating OWC wave energy converter, a BBDB. Using the acceleration potential method associated with pneumatic pressured chamber and viscous energy loss, the hydrodynamic performance of the floating BBDB was simulated in the time domain. The developed NWT technique was based on the potentialfluid theory and the boundary element method with constant panels on the boundaries, on which fully nonlinear freesurface and moving body boundary conditions were applied. A mixed EulerianLagrangian (MEL) scheme was used for a timevarying nonlinear freesurface treatment, along with the RungeKutta fourthorder time integration scheme, as a timemarching approach. To accurately predict the time derivative of the velocity potential for the freely floating BBDB, the mode decomposition scheme combined with the acceleration potential was implemented to calculate the body force and displacement by simultaneously solving the equations of body motion and the velocity potential.
In order to consider the effect of pneumatic pressure acting on the freesurface inside the chamber, the damped free surface condition was applied to the NWT technique, which was first presented by Evans (1982), Sarmento and Falcao (1985) and Falnes and McIver (1985). A linear relation between the chamber pressure and airflow velocity at the air duct, as observed in various experiments, was used to model the OWC chamber (Gato and Falcao, 1988; Suzuki and Arakawa, 2000). The timevarying pneumatic pressure, due to instantaneous airflow velocity interacting with freesurface fluctuation, was numerically modeled in an OWC system. The energy loss due to viscous flow at the entrance of the BBDB hull, which could be amplified by the body motions, was also modeled by imposing an artificial viscous damping coefficient upon the free surface inside the chamber. In the potential flow calculation, the viscous damping coefficient can be obtained by a comparison of the experimenttal data in open chamber conditions, which represents no pneumatic chamber pressure. The difference of chamber free surface elevation between the potentialfluidbased numerical results and the experimental data in open chamber conditions can be interpreted as the viscous damping.
MATHEMATICAL FORMULATION
 Boundary value problem for a floating BBDB
In order to simulate a freely floating BBDB associated with a pneumatic chamber and viscous damping, the mixed boundary value problem has to be solved using proper boundary conditions. The boundaries on the computational domain are the free surfaces inside and outside of the chamber, the body surface, the incident wave boundary, the rigid sea bottom, and the rigid endwall. Using the velocity potential (
φ
) and a continuity equation that satisfied the entire fluid domain, the Laplace equation (Eq. (1)) was chosen as a governing equation under the assumption of inviscid, incompressible, and irrotational flows. However, a special damping can be applied on the free surface inside the chamber to represent the viscous energy loss due to body shape and motion. Using the Green function (
G
), the Laplace equation can be transformed into a boundary integral equation (Eq. (2)) that is solved with the corresponding boundary conditions on the respective boundaries:
where
G_{ij}
(
x_{i}
,
z_{i}
,
x_{j}
,
z_{j}
)= (1/2
π
)ln
R
_{1}
for a 2D problem, the solid angle
α
is 0.5 for this study, and
R
_{1}
is the distance between the source (
x_{j}
,
z_{j}
) and the field points (
x_{i}
,
z_{i}
) on the boundary of the entire fluid domain.
The rigid sea bottom and the vertical endwall of the computational domain can be described by applying a nopenetration condition so that the water particle velocity in the normal direction is zero on the rigid boundary.
The moving body boundary condition can be expressed by the fact that the water particle velocity adjacent to the body surface is the same as the body velocity in the normal direction:
where
V_{B}
is the body velocity at the center of gravity. When the body is stationary (
V_{B}
= 0), Eq. (4) becomes the rigid boundary condition. For incident waves, an analytic linear incident wave profile (Eq. (5)) was fed on the left vertical input boundary.
where
A
,
ω
,
k
, and h are the wave amplitude, frequency, wave number, and water depth, respectively. The linear incident wave profile can be simply replaced by Stokes’ 2
^{nd}
order wave profile in the case of higher wave steepness. Fully nonlinear dynamic and kinematic freesurface boundary conditions are described as:
where
η
is the freesurface elevation and
P_{a}
is the air pressure on the freesurface, which is set to zero (atmospheric pressure) outside the chamber. However, the pneumatic pressure should be imposed inside the chamber when the air duct is installed.
The mixed EulerianLagrangian (MEL) approach was adopted to capture timevarying nonlinear freesurfaces, by which the freesurface boundary conditions can be modified using the total time derivative
Since the node on the free surface was designed to follow the water particle velocity (
, materialnode approach and
is node velocity), the fully nonlinear freesurface boundary conditions were transformed in the Lagrangian frame:
where
is the node location (
x
,
z
). A RungeKutta 4
^{th}
order (RK4) scheme was used to integrate the timedependent free surface boundary conditions (socalled timemarching) in order to obtain the velocity potential (
φ
) and surface elevation (
η
) at every time step.
During the nonlinear timedomain simulation, nodes on the floating body and the free surfaces must be updated and rearranged using a regridding scheme to avoid the numerical instability caused by the local accumulation of material nodes on the surface. The regridding process was carried out at every time step for a robust timedomain simulation. An artificial damping zone, as a numerical beach, was placed near the end of the fluid domain to absorb the transmitted waves and to prevent the reflected waves at the end wall. The damping coefficients were applied to both the dynamic and kinematic free surface boundary conditions (Eq. (8)).
where
μ
_{1}
and
μ
_{2 }
are damping coefficients, whose values are properly selected depending on the wave conditions.
A frontal damping zone was also installed on the free surface in front of the incident wave boundary to prevent a rereflection at the wave maker, which enables long time simulation of surfacepiercing bodies. This special damping scheme has to be designed so as to suppress only the reflected waves from the floating body, while preserving the original incident waves. In this regard, the adopted damping term should be applied to the difference between the total waves and incident waves. Then, the free surface boundary conditions can be expressed as:
where
μ
_{1f}
and
μ
_{2f}
are the frontal damping coefficients for absorbing the reflection waves from the body, and the reference values of (
∂φ
/
∂n
)* and
η
* were computed using the same computational conditions in the absence of bodies. For the case of moderate nonlinear incident waves, proper analytic solutions (e.g., secondorder Stokes waves) can be used in the practical applications, as determined by Tanizawa and Naito (1997). For the efficient absorption of wave energy, a spatial ramp function was applied to the selected damping coefficients (
μ
_{1}
,
μ
_{2}
and
μ
_{1f}
,
μ
_{2f}
) over the damping zone of two wavelengths (2λ). The direction of the ramp to the target value in the frontal damping is opposite to the case of rear damping zone. In the present study, open sea conditions can be realized in the computational domain by using both the frontal damping and the numerical beach.
A temporal ramp function at the input boundary was applied during the first two wave periods (2T) to prevent the impulselike behavior of the wave maker and to reduce the corresponding unnecessary transient waves. A Chebyshev fivepoint smoothing scheme was used along the free surface during the time integration to avoid the nonphysical, sawtooth numerical instability on a highly nonlinear free surface. The smoothing scheme was applied at every fifth time step for a stable time simulation. More detailed explanations of the general numerical schemes and formulations in the NWT technique are given in a previous study (Koo and Kim, 2004).
 OWC chamber modeling associated with power takeoff and viscous damping
The timevarying pneumatic pressure caused by the oscillating water column can be described as the rate of change of air volume inside the chamber (Kim and Iwata, 1991). The rate of air volume change is also directly related to the relative spatialmean vertical velocity between the body and free surface inside the chamber. The pneumatic pressure in the chamber at each time step is expressed as (Koo and Lee, 2011; Koo et al., 2012):
where
C_{dm}
is an equivalent linear pneumatic coefficient determined by the air ductchamber area ratio, Δ
V
/Δ
t
(
V_{t}

V_{tΔt}
)Δ
t
is the timedifferential change of air volume in the chamber (
V_{t}
is the volume at the current time step, while
V_{tΔt}
is the volume at the previous time step),
A_{d}
is the sectional area of the air duct, and
U_{d}
(
t
) is a timevarying airflow velocity at the air duct under the assumption of incompressible air in the case of a relatively large air opening. Sarmento and Falcao (1985) pointed out the springlike effects, due to air compressibility in the OWC chamber, for a relatively large air volume compared with a small nozzle outlet, which was not considered in the present analysis. Therefore, the case of this study can represent the largediameter air turbine system, such as Wells turbine and Impulse turbine. Koo and Kim (2010) also studied pneumatic pressure formulation and its application to a fixed OWC chamber.
In addition to pneumatic pressure, a pressure drop induced by viscous energy loss occurs inside the chamber because a viscousflow phenomenon, such as vortex shedding, arises at the corner of the BBDB. When the pistonlike movement of chamber surface elevation occurs, as in a case of resonance frequencies, and the ensuing high velocity flows are generated, the magnitude of the viscous energy loss greatly increases. In this regard, the difference in the chamber surface elevation between the potentialflowbased numerical solutions and the experimental results can be interpreted as viscous energy loss (work done). A typical entry pressure drop
where
K_{L}
is the loss coefficient and
V
is the flow velocity, can be linearized, so that the flow velocity is proportional to the vertical velocity of the water column
Then, the pressure drop in the chamber can be described as (Koo and Kim, 2010):
where
denotes a modified energy loss coefficient,
represents the timevarying spatialmean surface velocity inside the chamber, and
ν
is an equivalent linear viscous damping coefficient.
Generally, the pneumatic pressure and viscous energy loss are approximately proportional to the square of the airflow velocity and water column velocity, respectively. Hence, Eqs. (10) and (11), which have an equivalent linear coefficient, can be modified to equations with a quadratic coefficient. The quadratic pneumatic (
P_{ac}
) and viscousloss pressure drop (
P_{v}
) can be applied straightforwardly to the present nonlinear timedomain simulation:
where
D_{dm}
and
ν_{q}
denote the quadratic pneumatic and viscousloss coefficients, respectively. Therefore, the final dynamic freesurface boundary conditions inside the chamber (including an equivalent linear or quadratic pneumatic pressure and viscous pressure drop) can be expressed, respectively, as:
 Energy conservation of the BBDB system
The BBDB system contains five energyflux components: input, reflection, transmission, energy extraction by airflow, and energy loss by viscosity. The transferred wave energy flux into the energy converter is proportional to the square of the wave height and the wave group velocity. Outside the chamber, the energyflux components per unit area are expressed as:
where
C_{g}
is the wave group velocity, and
H_{i}
,
H_{R}
, and
H_{T}
represent the incident, reflected, and transmitted wave heights, respectively.
E_{i}
,
E_{R}
, and
E_{T}
also denote the respective energy components. Inside the chamber, the available pneumatic power in Watts (
AP
) and the rate of energy dissipation (
DR
) based on the linear (or quadratic) coefficients can be simplified as (Koo and Kim, 2010):
where
B_{dm}
is the chamber width, and
V
_{0}
is the amplitude of the relative spatialmean vertical velocity of the free surface. Therefore, the energy conservation of the total energy flux in the BBDB system is given by:
 Acceleration potential method for a freely floating BBDB
The time derivative of the velocity potential
φ_{t}
is a critical value when obtaining the accurate pressure and force acting on a freely floating body, as determined by Koo and Kim (2004). The
φ_{t}
on the bodysurface node has to be obtained by simultaneously solving the fluid particle and body motion equations. Thus, the use of the acceleration potential method for obtaining
φ_{t}
is the most accurate and consistent way to predict the motion of a freely floating body. The wave force on a floating body can be calculated by integrating Bernoulli’s pressure over the instantaneous wetted body surface. Including the gravitational and external restoring spring forces, if any, the total force in the ith direction for a 2D body can be calculated as follows:
where
K
and
C
represent the additional horizontal spring and damping coefficients, respectively.
x
is the surge,
z
is the heave,
is the horizontal body velocity,
W
is the weight of the body (=
mg
),
δ_{ij}
denotes the Kronecker delta function (horizontal,
i
= 1; vertical,
i
= 2), and
S_{B}
is the wetted body surface.
In order to implement the acceleration potential method, we used the mode decomposition method introduced by Vinje and Brevig (1981) to obtain the time derivative of the velocity potential directly from the acceleration field. In a 2D problem, the acceleration field can be decomposed into four modes, corresponding to three unit accelerations for surgeheavepitch and the acceleration due to the velocity field. Hence, each mode can be obtained by solving the respective boundary integral equation in the acceleration field. Using these four modes and the equation of body motion, the body acceleration can be determined.
The time derivative of the velocity potential
φ_{t}
on a floating body for a 2D problem is given by:
where
a_{i}
is the
i
th mode component of generalized body acceleration (1 = surge, 2 = heave, 3 = pitch).
To solve the boundary integral equation for each mode in the acceleration field, the body boundary condition is given as:
where,
n
_{1}
n
_{2}
and
n
_{3}
denote the unit normal vector in the respective direction (1 = surge, 2 = heave, 3 = pitch), and
q_{B}
represents the contribution of the velocity field to the acceleration field, which can be written for a 2D simulation as follows (Koo and Kim, 2004):
where
n
and
s
are normal and tangential unit vectors of the body surface, respectively,
k_{n}
= −1/
ρ
* is the local curvature along the
s
direction,
x_{P}
,
z_{P}
denote the position of point P from a body’s center of gravity, and
υ_{X,Z}
and
ω
_{0}
are the translational and angular velocities of a body, respectively.
The freesurface boundary conditions for the respective free surface areas of the BBDB system in the acceleration field are also given as:
where
ϕ_{i}
for
i
= 1~3 modes is zero for the free surface boundary conditions in the acceleration field. The input boundary and other rigid boundary conditions are given as:
After solving the boundary integral equation in the acceleration field for the aforementioned conditions,
ϕ_{i}
for
i
= 1~4 on the body surface (Eq. (18)) is obtained. To calculate the generalized acceleration (
a_{i}
) in Eq. (18) for each mode (surge, heave, and pitch), the force equation should be combined with Newton’s second law. Including the gravitational force and external force (e.g., horizontal spring), the total force and moment in the ith direction can be calculated as follows:
where
,
m
is the body mass,
I
is mass moment of inertia.
From Eqs. (23) to (25), three unknown accelerations (
a_{i}
) of each mode can finally be determined. Using the RungeKutta Nystrom 4
^{th}
order integration method, the body acceleration, velocities and displacements are then consecutively determined. The calculated translational and rotational body displacement is used to update the body geometry for the next time step.
In the present study, a freely floating BBDB with no additional mooring lines was used to evaluate the hydrodynamic performance. Therefore, neither horizontal spring nor damping coefficients was applied to Eq. (23). In the fully nonlinear simulation, the instantaneous position of the BBDB was updated at every time step and the body drifts during the time simulation.
EXPERIMENTAL SETUP
In order to verify the calculated results from the developed numerical model, an experiment was conducted independently in the 2D wave tank at the University of Ulsan (Koo and Lee, 2011). The wave tank dimensions were 0.5 × 0.4 × 35
m
(width × depth × length). A flaptype wave maker and a sloped beach were positioned at both ends of the tank. Six wave probes were installed to measure surface elevations: three in the weather side, two in the lee side, and one in the middle of the chamber. The measured signals of surface elevations were amplified by the AMP and gathered into the DAQ board.
Table 1
shows the specifications of the BBDB model used in Koo and Lee (2011). Since the width of BBDB model is 0.48 m compared with the tank width (0.5
m
), there is no side wall effect. In order to avoid friction during the body motions in a 2D wave tank, four bearings were mounted on each corner of the model along the wall of the wave tank, by which the BBDB model occupied the whole width of the tank and a threedegreeoffreedom motion can be allowed. According to the tank size and the performance of the wave maker, incident waves of periods from 0.6
s
to 1.4
s
with a height of 0.01
m
were generated.
Since the wave maker used in the present experiment cannot control the reflected waves from the floating body, the rereflection may occur when the duration of wave generation is long. Therefore, the steadystate waves with no contamination from the rereflected waves are only used to compare the numerical results.
Specifications of the BBDB model.
Specifications of the BBDB model.
NUMERICAL RESULTS AND DISCUSSION
A schematic diagram and coordinate system for the present numerical model is shown in
Fig. 1
. The flowchart of computational algorithm is also shown in
Fig. 2
. The dimensions of the BBDD are the same as those in the experimental model (
Table 1
). Each side of the free surface (weather and lee side) is four times the incident wavelength ( 4
λ
), including the artificial damping zone of the two wavelengths ( 2
λ
), to dampen the designated waves sufficiently at the end of the fluid domain. A flat rigid sea bottom was set at a depth of 0.4 m. In order to prevent a rereflection from the wave maker, a special damping scheme on the free surface, called the frontal damping zone (Koo and Kim, 2004) of two wavelengths ( 2
λ
), was also installed near the wave maker zone. From the convergence test for various node numbers and time step intervals (
Fig. 3
), 20 nodes (or more) per incident wavelength and a time step of
dt
=
T
/ 64 (or less) interval were used to produce the numerical results.
Overview of computational domain.
Flowchart of the computational algorithm.
Convergence test for node number and time segment (comparison of surface elevation inside the chamber).
Initially, the calculated wave elevations of the fixed BBDB were compared with the experimental data (
Fig. 4
) to determine the viscous energy loss due to the body shape. The linear and nonlinear results were also compared to evaluate the effect of nonlinear interactions between wave and body. A significant difference was observed between the numerical and experimental results near the resonance frequencies, where the incident wavelength was approximately six times greater than the total body length (
L
). Some discrepancy between numerical results and experimental data was also found at longer wavelength (
λ
/
L
> 8 ), which may be due to a measurement error at the tank experiment and the effect of wave reflection at the endwall of the physical wave tank installed with wave absorbing material. In general, it is difficult to fully control the wave reflection at the downstream of the tank when the incident wave is longer than a certain limitation of the tank capacity. The wave tank used in this study does not have a sufficient capacity of wave absorbing system at the downstream.
Since a small incident wave steepness was applied to the entire frequency range ( 0.00243 <
H
/
λ
< 0.0131 , H = 0.01
m
), there was little difference between the linear and nonlinear results. When the incident wave becomes higher (H = 0.02
m
and 0.04
m
) in the nonlinear calculation, the chamber surface elevation does not change much compared to the case of H = 0.01
m
, which implies that the effect of wave nonlinearity is not significant in the fixed body. Since this study is mainly focused on the development of numerical tool to simulate a floating OWC system, BBDB and to validate the numerical results compared with the experiment, the numerical simulation with large waves and the violent body motions will be studied as the next research topic.
The discrepancy between the numerical and experimental results in the resonance frequency was caused by viscous loss. Energy loss due to vortex generation may occur at the sharp corners of the BBDB.
In order to compensate for the viscous energy loss with the artificial damping in the fixed BBDB system, a small viscous damping coefficient (
ν
= 0.09) was applied to the freesurface boundary inside the chamber (
Fig. 5
). The damping coefficient chosen matched the experimental results in the open chamber condition, where the difference between the numerical results and experimental data was caused by viscous energy loss. The viscous damping coefficient was divided by the water density (
ρ
= 1000).
Fig. 4
shows a comparison of surface elevations with applications of viscous damping on the free surface. Due to the relatively narrow chamber, compared to wavelength, the pumping mode of the free surface elevation was observed over the whole range of frequencies. The simulated elevations with a chosen viscous damping coefficient were in agreement with the experimental values. The applied damping coefficient only affected the resonance frequency region, with the other regions being only minimally affected.
A comparison of surface elevations at the weather and lee sides is shown in
Fig. 6
. All data were measured at 1
m
ahead of the body. The calculated elevations, with a tuned viscous coefficient (
ν
= 0.09), were in good agreement with the experimental
Comparison of open chamber elevations for the experimental, linear and nonlinear calculations of the fixed body case (H = 0.01, 0.02, 0.04 m).
Comparison of open chamber elevations of the experimental and nonlinear calculations for the case of a fixed body (H = 0.01 m).
results. It is also observed that a slight difference between the cases of viscous damping and no damping, which implies that the effect of a small viscous damping coefficient (
ν
= 0.09) on the free surface elevation outside the chamber is not significant. Since a reflected wave cannot be controlled in the experiment, especially in the case of a surface piercing body, a slight deviation can be found between the two results in the weather side. The linear and nonlinear calculations were similar, because of the small incident wave steepness (similar to the effect seen in
Fig. 4
).
Total energy conservation of the fixed BBDB system is compared in
Fig. 7
. For the open chamber conditions (
C_{dm}
= 0), the sum of the three energy fluxes of reflection, transmission, and viscous energy loss was equal to the incident wave energy flux. In the case of the pneumatic chamber (
C_{dm}
= 0.7), the sum of all of the energy components, including the pneumatic energy due to air flow, was conserved over the whole range of frequencies. Comparing the respective energy components showed that the total energy of the fixed BBDB system was conserved in the computational domain during the time simulation. A small discrepancy was observed in the shorter wave length, which may be due to the numerical measurement of large reflected waves in front of BBDB, with an insufficient number of probes on the free surface.
Comparison of wave elevations for the experimental, linear and nonlinear calculations with a viscous damping coefficient ((ν= 0.09) for the fixed body case (H = 0.01 m).
Comparison of energy components with a viscous damping coefficient (ν= 0.09) and energy extraction (C_{dm}= 0.7) in the fixed body case.
For the floating BBDB system,
Fig. 8
shows the comparison of chamber surface elevations for the experimental and numerical calculations. With a tuned linear equivalent viscous damping coefficient (
ν
= 0.52), the nonlinear calculation was in agreement with the experimental data, which was also determined by Koo and Lee (2011). With a quadratic viscous damping coefficient (
ν_{q}
= 13), the calculated results also matched the cases of the linear damping coefficient and experimental data. Practically, there is no fundamental difference between the use of a linear viscous damping and a quadratic viscous damping during the present simulation, because a proper damping coefficient can be obtained from the experiment. In this study, the comparison of the damping effect between the linear and quadratic coefficients was evaluated and the characteristics of the results were investigated. More detailed studies on the effect of viscous damping type, its proper selection from various experimental data, and the general application of the damping coefficients should be conducted as the next research topics.
The proper viscous damping can be obtained from a comparison of the experimental data with the open chamber conditions, where the viscous energy loss only exists without pneumatic pressure. The viscous damping coefficients shown in the figures were divided by water density ((
ρ
= 1000). The incident wave height of 0.01
m
was used for the linear calculation, whereas the height of 0.005
m
was chosen for nonlinear calculation, so as to prevent numerical instability attributed to the extreme nonlinear body motions near the resonance frequencies. Compared with the fixed BBDB case, the deviation between the experimental and numerical results without viscous damping is significant near the resonance frequencies (
λ
/
L
>≈ 4.5 ~ 5.0 ). Thus, it can be inferred that viscous energy loss may greatly increase due to the motion of the BBDB. In particular, the pitch of the BBDB could intensify the vortex generation at the corner of the body.
Comparison of open chamber surface elevations for experimental and numerical calculations with either a linear (ν= 0.52) or a quadratic ( ν_{q} = 13.0) viscous damping coefficient in the case of a freely floating BBDB.
Fig. 9
shows a comparison of the surface elevations at the weather and lee sides of the BBDB with the application of a linear viscous damping coefficient inside the chamber free surface. The calculated elevations agree reasonably well with those from the experimental data. Due to the difficulty in controlling the reflected waves in the floating body experiment, the measured elevations deviate from the numerical results at some frequencies, especially in the long wave region. Since a relatively large motion of the BBDB magnifies the radiated waves, the nonlinear calculation results with updated body and freesurface motions showed greater accuracy than the linear results calculated with a meanpositioned body and free surface. The magnitude of elevation at the lee side is smaller near the resonance frequency region (
λ
/
L
≈ 4.5 ~ 5.0 ) than at the other frequencies. This phenomenon can be explained by the fact that a greater portion of incident wave energy is damped by viscosityinduced loss due to high vertical flow velocities, as well as large body motions. Thus, the transmitted wave energy becomes small. As mentioned earlier, a slight discrepancy of the total energy between the input and total summated energy may be due to an insufficient number of wave probes during the numerical measurement on the free surface including numerical aliasing error.
Comparison of wave elevations outside the camber in the case of a freely floating BBDB with a linear equivalent damping coefficient (ν= 0.52) (Koo and Lee, 2011).
Comparison of energy components of the freely floating BBDB with viscous damping coefficient ((ν = 0.52) and open chamber (C_{dm} = 0.0) and pneumatic chamber conditions (C_{dm} = 0.7).
The respective energy components of the floating BBDB for both chamber conditions are compared in
Fig. 10
. In the open chamber conditions (left graph), as pointed out by Koo and Lee (2011), the relative energy loss by viscosity near the resonance frequencies is up to 80% of the incident wave energy, which indicates that a large amount of energy could be lost at the sharp corners of the body by the resonant pitchinduced vortices. In the pneumatic chamber conditions (right graph), the relative energy loss can be reduced by 50% of the total energy. This can be explained by the fact that the motion of the BBDB decreesed due to the effect of the pneumatic pressured chamber, resulting in a diminished vortex generation. In both chamber conditions, the total energy in the computational domain was conserved over the whole range of frequencies. Comparing with
Fig. 7
, the magnitude of the pneumatic energy in the floating case (about 25% of incident energy) was greater than that in the fixed case (about 10%). Thus, the body motion (heave and pitch) may amplify the relative vertical water column velocity, thus increasing the resultant airflow velocity.
A comparison of the available wave power with various pneumatic coefficients (
C_{dm}
) is shown in
Fig. 11
. As the coefficient increases, indicating that the ratio of chamber surface area to air duct area increases, the airflow velocity at the nozzle outlet increases until the air opening is relatively large, under the assumption of incompressible air. The maximum available pneumatic power can be found near the resonance frequencies.
Comparison of available wave power with various extraction coefficients.
CONCLUSIONS
The hydrodynamic performance of a floating OWC device was evaluated in the time domain. An acceleration potential method using fullupdated matrix calculation, associated with a mode decomposition scheme, was implemented to obtain the hydrodynamic force and displacement of a freely floating BBDB. The developed NWT technique was based on the potential theory and the boundary element method with constant panels on the boundaries, on which fully nonlinear freesurface and moving body boundary conditions associated with OWC power takeoff were applied. The modeling of the pneumatic chamber and viscous damping was proportional to the square of the airflow velocity and water column velocity, respectively. However, the numerical model can also be simplified to a linear relation. Open sea conditions can be realized in the computational domain by using the frontal damping and the numerical beach. Therefore, the NWT technique was better than the experiment in controlling the reflected waves from the wave maker and floating body.
A significant difference of chamber surface elevation was observed between the numerical and experimental results near the resonance frequencies, which may be attributed to the viscous energy loss. The energy loss by vortex generation may increase at the sharp corners of the BBDB. The numerical results under the open chamber conditions were adjusted by adding a proper viscous damping coefficient, obtained from a comparison of the experimental data. The calculated surface elevations affected by pneumatic pressure correlated reasonably well with the experimental values.
The total energy summation of the BBDB system in the computational domain was found to be conserved over the whole range of frequencies. In the case of a floating BBDB with an open chamber, the relative energy loss due to viscosity near the resonance frequencies was up to 80% of the incident wave energy, which indicates that a large amount of energy could be lost at the sharp corners of the body by the resonant pitchinduced vortices.
Since a relatively large motion of the BBDB magnified the radiated waves, the nonlinear results that were calculated using the updated freesurface and body motions were found to be more accurate than the linear ones with a meanpositioned free surface and body.
The developed NWT model for a floating BBDB device can be used to analyze the hydrodynamic performance of a floatingtype OWC system. Using a parametric study for various conditions compared with the corresponding experiments, it is possible to predict the motion characteristics and determine the proper BBDB shape for various environmental conditions.
Brendmo A.
,
Falnes J.
,
Lillebekken P.M.
1996
Linear modelling of oscillating water columns including viscous loss.
Applied Ocean Research
18
(2)
65 
76
DOI : 10.1016/01411187(96)000119
Delaure E.M.C.
,
Lewis A.
2003
3D hydrodynamic modelling of fixed oscillating water column wave power plant by a boundary element methods.
Ocean Engineering
30
(3)
309 
330
DOI : 10.1016/S00298018(02)00032X
Evans D.V.
1982
Wave power absorption by systems of oscillating surfacepressure distributions.
Journal of Fluid Mechanics
114
481 
499
DOI : 10.1017/S0022112082000263
Falnes J.
,
McIver P.
1985
Surface wave interactions with systems of oscillating bodies and pressure distributions.
Applied Ocean Research
7
(4)
225 
234
Gato L.M.C.
,
Falcão A.F.O.
1988
Aerodynamics of the wells turbine.
International Journal of Mechanical Science
30
383 
395
Heath T.
,
Whittaker T.J.T.
,
Boake C.B.
2000
The design, construction and operation of the LIMPET wave energy converter.
Proceedings of 4th European Wave Energy Conference
Aalborg, Denmark
46 December
Hong D.C.
,
Hong S.Y.
,
Hong S.W.
2004a
Numerical study of the motions and drift force of a floating OWC device.
Ocean Engineering
31
(2)
139 
164
DOI : 10.1016/S00298018(03)001185
Hong D.C.
,
Hong S.Y.
,
Hong S.W.
2004b
Numerical study on the reverse drift force of floating BBDB wave energy absorbers.
Ocean Engineering
31
(10)
1257 
1294
DOI : 10.1016/j.oceaneng.2003.12.007
Imai Y.
,
Toyota K.
,
Nagata S.
,
Setoguchi T.
,
Matsunaga N.
,
Manago Y.
,
Shimozono T.
2009
Experimental study on negative drift force acting on a floating OWCtype wave energy converter “Backward Bent Duct Buoy”.
Proceedings of 19th International Offshore and Polar Engineering Conference, ISOPE
Osaka, Japan
2126 June
331 
338
Kim D.
,
Iwata K.
1991
Dynamic behavior of tautly moored semisubmerged structure with pressurized airchamber and resulting wave transformation.
Coastal Engineering in Japan
34
(2)
223 
242
Kim J.H.
,
Lew J.H.
,
Hong D.C.
,
Hong S.W.
2006
A study on motion and wave drift force of a BBDB type OWC wave energy device.
Journal of Ocean Engineering and Technology
20
(2)
22 
28
Kim J.H.
,
Lew J.H.
,
Hong D.C.
,
Kim Y.S.
,
Choi H.S.
,
Hong S.W.
2007
A study on motion of a BBDB type OWC wave energy device considering pneumatic damping coefficients in the duct.
Proceedings of 17th International Offshore and Polar Engineering Conference, ISOPE
Lisbon, Portugal
16 July
483 
488
Koo W.C
,
Kim M.H.
2010
A nonlinear timedomain simulation of a landbased oscillating water column.
Journal of Waterway, Port, Coastal and Ocean Engineering
136
(5)
276 
285
DOI : 10.1061/(ASCE)WW.19435460.0000051
Koo W.C.
,
Lee K.R.
2011
Numerical and experimental analysis of Backward Bent Duct Buoy (BBDB) wave energy converter.
Proceedings of the 21st International Offshore and Polar Engineering Conference, ISOPE
Hawaii, USA
June 2011
655 
660
Koo W.C.
,
Kim S.J.
,
Kim M.H.
2012
Numerical analysis of a floating wave energy converter, Backward Bent Duct Buoy (BBDB) with different shapedcorners.
Offshore and Arctic Engineering
Proceedings of 31st International Conference on Ocean
Rio de Janeiro, Brazil
16 July
OMAE83666
669 
674
Masuda Y.
1971
Wave activated generator.
Proceedings of International Colloquium on Exposition of the Oceans
Bordeaux, France
Masuda Y.
,
Yamazaki T.
,
Outa Y.
,
McCormick M.E.
1987
Study of backward bent duct buoy.
Proceedings of Oceans ’87, IEEE
384 
389
McCormick M.E.
2007
Ocean wave energycConversion.
Dover Publication
Nagata S.
,
Toyota K.
,
Imai Y.
,
Setoguchi T.
2008
Experimental study on hydrodynamic forces acting on a floating wave energy converter ‘Backward Bent Duct Buoy’.
Proceedings of 18th International Offshore and Polar Engineering Conference, ISOPE
Vancouver, BC, Canada
611 July
366 
373
Nagata S.
,
Toyota K.
,
Imai Y.
,
Setoguchi T.
2009
Numerical simulation for evaluating of primary energy conversion of floating OWCtype wave energy converter.
Proceedings of 19th International Offshore and Polar Engineering Conference, ISOPE
Osaka, Japan
300 
307
Park M.
,
Koo W.
,
Choi Y.
2010
Hydrodynamic interaction with an array of porous circular cylinders.
International Journal of Naval Architecture and Ocean Engineering
2
(3)
146 
154
DOI : 10.3744/JNAOE.2010.2.3.146
Sarmento A.
,
Falcão A.F.O.
1985
Wave generation by an oscillating surfacepressure and its application in waveenergy extraction.
Journal of Fluid Mechanics
150
467 
485
Suzuki M.
,
Arakawa C.
2000
Guide vanes effects of wells turbine for wave power generator.
International Journal of Offshore and Polar Engineering
10
(2)
153 
159
Suzuki M.
,
Kuboki T.
,
Nagata S.
,
Setoguchi T.
2011
Numerical investigation of 2D optimal profile of backwardbent duct type wave energy converter.
Journal of Offshore Mechanics and Arctic Engineering
133
041602 
1
DOI : 10.1115/1.4003519
Tanizawa K.
,
Naito S.
1997
A study on parametric roll motions by fully nonlinear numerical wave tank.
Proceedings of 7th International Offshore and Polar Engineering Conference, ISOPE
Honolulu, USA
2530 May
69 
75
Toyota K.
,
Nagata S.
,
Imai Y.
,
Setoguchi T.
2008
Effects of hull shape on primary conversion characteristics of a floating OWC ‘Backward Bent Duct Buoy’.
Journal of Fluid Science and Technology
3
(3)
458 
465
DOI : 10.1299/jfst.3.458
Toyota K.
,
Nagata S.
,
Imai Y.
,
Setoguchi T.
2009
Research for evaluating performance of OWCtype wave energy converter ‘Backward Bent Duct Buoy’.
Proceedings of 8th European Wave and Tidal Energy Conference
Uppsala, Sweden
710 September
901 
913
Toyota K.
,
Nagata S.
,
Imai Y.
,
Oda J.
,
Setoguchi T.
2010
Primary energy conversion characteristics of a floating OWC ‘Backward Bent Duct Buoy’.
Proceedings of 20th International Offshore and Polar Engineering Conference, ISOPE
Beijing, China
2026 June
850 
855
Wang D.J.
,
Katory M.
,
Li Y.S.
2002
Analytical and experimental investigation on the hydrodynamic performance of onshore wavepower devices.
Ocean Engineering
29
(8)
871 
885
DOI : 10.1016/S00298018(01)000580