Energy Efficiency Analysis of a Deformable Wave Energy Converter Using Fully Coupled Dynamic Simulations

: Deformable wave energy converters have significant potential for application as flexible material that can mitigate structural issues, while how to design the dimensions and choose an optimal deployment location remain unclear. In this paper, fully coupled computational fluid dynamics and computational solid mechanics were used to simulate the dynamic interactions between ocean waves and a deformable wave energy converter. The simulation results showed that the relative length to wave, deployment depth and aspect ratio of the device have significant effects on the energy conversion efficiency. By calculating the energy captured per unit width of the device, the energy efficiency was found to be up to 138%. The optimal energy conversion efficiencies were achieved when the structure length was 0.25, 0.5 or 0.75 of the dominating wavelength and sub-merged at a corresponding suitable depth. The aspect ratio and maximum stress inside the wave energy converter showed a nonlinear trend, with potential optimal points revealed. The simulation approach and results support the future design and optimisation of flexiable wave energy converters or other marine structures with notable deformations.


Introduction
Ocean waves have great development potential as a renewable energy source.They are ubiquitous in the ocean, a major part of the world [1], and the energy density of ocean waves is also higher than that of solar energy and wind energy due to water weight [2].Devices that convert wave energy into a form of energy that humans can utilise are called wave energy converters (WECs).In the process of WEC research and development, there are three major challenges: cost of energy, efficiency, and survivability [3].However, the current WECs are not ready for commercial operation, mainly because the current structures cannot survive in the marine environment in an affordable way, and thus cannot meet the needs of the three challenges mentioned above [4].Adopting flexible materials in WECs is considered a viable technical route that can reduce the impact of wave loads on WECs to improve survivability in harsh sea states and use deformation to transform energy significantly [5].
A wide variety of flexible WECs exist, as shown in Figure 1.The bulge wave device is an example of a device that utilises the feature of flexible material to push airflow, which was discovered by Farney et al. [6].Alam [7] introduced a WEC concept that uses a synthetic seabed carpet to imitate the absorption behaviour of muddy seafloors.The synthetic seabed carpet is expected to respond like film to the action of waves, thus transferring the deformation into energy through pumping [8].A poly-OWC is a device that combines an oscillating water column (OWC) and dielectric elastomer generators (DEGs) [9,10].Point absorbers such as OPT PowerBuoy and Wavebob [11] can also incorporate flexible materials to improve their ability to adapt to harsh sea states [12].

228
Oceans 2024, 5 2 absorbers such as OPT PowerBuoy and Wavebob [11] can also incorporate flexible materials to improve their ability to adapt to harsh sea states [12].For this wide variety of WEC devices, many scholars have evaluated their energy capture capabilities.Zheng et al. [13] built an elastic disk-shaped WEC model based on linear potential flow theory and the eigenfunction matching method.The wave energy absorption was evaluated by taking into account the power take-off system of the device.Similar structures were investigated experimentally in [14].Another case of assessing WEC energy absorption by mathematical modelling is [15].In this study, the energy capture capability of a clam-shaped WEC device deployed in floating mode was evaluated by considering different sea states and mooring modes.Benites et al. [16] used an Open-FOAM-based model to evaluate the relationship between the energy absorption capacity of an oscillating wave surge converter (OWSC) and wavelength.
For the simulation of ocean structures such as WECs, the boundary element method (BEM) [11,17,18] and computational fluid dynamics (CFD) are widely used.The former usually requires less computational resources and has a faster computational speed [19].However, the calculation accuracy under nonlinear waves is not reliable [20].From the perspective of the long-term development of WECs, the BEM does not meet the demand of WECs in the design process [21].Although the latter solves the problem of calculation accuracy under strong nonlinear states and obtains satisfactory results [22][23][24], it does not trigger the interaction between the deformation of WECs and the fluid.
The functionality of flexible WECs involves complex structural deformations that are based on fully coupled fluid-structure interaction (FSI), i.e., any deformation of the structure triggers a response of the wave flow and vice versa; in the WEC field, a modelling method with such capability has not been reported in the literature, which has been pointed out as a crucial gap by the recent review of Collins et al. [5].This motivated the present work to develop a tool to analyse the FSI and energy efficiency of deformable WECs.
Aiming to bridge the modelling gap of flexible WECs, a fully coupled CFD+CSM approach is adopted in this paper.This approach can simulate nonlinear wave conditions, complex structural deformations, and their transient interactions, which has been For this wide variety of WEC devices, many scholars have evaluated their energy capture capabilities.Zheng et al. [13] built an elastic disk-shaped WEC model based on linear potential flow theory and the eigenfunction matching method.The wave energy absorption was evaluated by taking into account the power take-off system of the device.Similar structures were investigated experimentally in [14].Another case of assessing WEC energy absorption by mathematical modelling is [15].In this study, the energy capture capability of a clam-shaped WEC device deployed in floating mode was evaluated by considering different sea states and mooring modes.Benites et al. [16] used an OpenFOAMbased model to evaluate the relationship between the energy absorption capacity of an oscillating wave surge converter (OWSC) and wavelength.
For the simulation of ocean structures such as WECs, the boundary element method (BEM) [11,17,18] and computational fluid dynamics (CFD) are widely used.The former usually requires less computational resources and has a faster computational speed [19].However, the calculation accuracy under nonlinear waves is not reliable [20].From the perspective of the long-term development of WECs, the BEM does not meet the demand of WECs in the design process [21].Although the latter solves the problem of calculation accuracy under strong nonlinear states and obtains satisfactory results [22][23][24], it does not trigger the interaction between the deformation of WECs and the fluid.
The functionality of flexible WECs involves complex structural deformations that are based on fully coupled fluid-structure interaction (FSI), i.e., any deformation of the structure triggers a response of the wave flow and vice versa; in the WEC field, a modelling method with such capability has not been reported in the literature, which has been pointed out as a crucial gap by the recent review of Collins et al. [5].This motivated the present work to develop a tool to analyse the FSI and energy efficiency of deformable WECs.
Aiming to bridge the modelling gap of flexible WECs, a fully coupled CFD+CSM approach is adopted in this paper.This approach can simulate nonlinear wave conditions, complex structural deformations, and their transient interactions, which has been successfully applied to the hydroelasticity of sea ice and breakwaters [25][26][27].On the basis of open source CFD software OpenFOAM (Version: Foam extend 4.0) [28], "solids4foam" is adopted in this research to realise fluid-structure interaction (FSI) simulations [29].Wave generation and absorption are achieved by the "Waves2foam" toolbox using the Relaxation Zone Method (RZM) [30].In order to determine the location of the interface between the gas phase and the liquid phase, the Volume of Fluid (VOF) [31] method is adopted to model multiphase flows.
The rest of the paper is organised as follows: in Section 2, the numerical theories and practicalities of building a fully coupled CFD+CSM model are introduced to simulate the dynamic response between ocean waves and the deformable structure.Afterwards, a series of model verification work is reported.Then, in Section 3, by changing the length of the flexible structure, key parameters are investigated to analyse their impacts on energy efficiency, such as the device length and the deployment depth.Wave-induced stress of the device is also investigated.Finally, Section 4 summarises the paper with its key conclusions, limitations, and recommendations for future work.

Computational Modelling
A Numerical Wave Tank (NWT) is the generic term used to describe numerical simulators that can model nonlinear free surface waves, hydrodynamic forces, and floating body motions [32].To calculate the dynamic behaviour of a flexible structure used for energy conversion under a wave load, an NWT is established that combines the functions of the FSI, wave generation and absorption, and the VOF method in this research (as shown in Figure 2).Even though this is only a two-dimensional model, it primarily focuses on the impact of incoming waves propagating along the landwards direction on the flexible structure.This study will investigate the energy efficiency of a unit width of WEC, where the two-dimensional assumption stands.Given that this dimensional simplification can significantly reduce computational costs, numerous simulation cases can be analysed, allowing for a comprehensive study of the influential parameters.The case setup and solution process are presented as follows.
successfully applied to the hydroelasticity of sea ice and breakwaters [25][26][27].On the basis of open source CFD software OpenFOAM (Version: Foam extend 4.0) [28], "solids4foam" is adopted in this research to realise fluid-structure interaction (FSI) simulations [29].Wave generation and absorption are achieved by the "Waves2foam" toolbox using the Relaxation Zone Method (RZM) [30].In order to determine the location of the interface between the gas phase and the liquid phase, the Volume of Fluid (VOF) [31] method is adopted to model multiphase flows.
The rest of the paper is organised as follows: in Section 2, the numerical theories and practicalities of building a fully coupled CFD+CSM model are introduced to simulate the dynamic response between ocean waves and the deformable structure.Afterwards, a series of model verification work is reported.Then, in Section 3, by changing the length of the flexible structure, key parameters are investigated to analyse their impacts on energy efficiency, such as the device length and the deployment depth.Wave-induced stress of the device is also investigated.Finally, Section 4 summarises the paper with its key conclusions, limitations, and recommendations for future work.

Computational Modelling
A Numerical Wave Tank (NWT) is the generic term used to describe numerical simulators that can model nonlinear free surface waves, hydrodynamic forces, and floating body motions [32].To calculate the dynamic behaviour of a flexible structure used for energy conversion under a wave load, an NWT is established that combines the functions of the FSI, wave generation and absorption, and the VOF method in this research (as shown in Figure 2).Even though this is only a two-dimensional model, it primarily focuses on the impact of incoming waves propagating along the landwards direction on the flexible structure.This study will investigate the energy efficiency of a unit width of WEC, where the two-dimensional assumption stands.Given that this dimensional simplification can significantly reduce computational costs, numerous simulation cases can be analysed, allowing for a comprehensive study of the influential parameters.The case setup and solution process are presented as follows.

Computational Domain and Boundary Conditions
In an NWT, the boundary conditions include the inlet, outlet, wall, open boundary, etc.In the fluid subdomain, the top boundary applies the static pressure boundary to represent the atmospheric conditions, and the bottom is determined as the no-slip boundary to simulate the seabed.The left and right boundary are the inlet and outlet, respectively.

Computational Domain and Boundary Conditions
In an NWT, the boundary conditions include the inlet, outlet, wall, open boundary, etc.In the fluid subdomain, the top boundary applies the static pressure boundary to represent the atmospheric conditions, and the bottom is determined as the no-slip boundary to simulate the seabed.The left and right boundary are the inlet and outlet, respectively.The x-axis is parallel to the seabed and points from the inlet to outlet.The y-axis is positive upwards.
In the solid subdomain, except for the fixed boundary on the left, the rest of the boundary can move freely under external loads.The information exchange between the fluid subdomain and the solid subdomain also takes place within these boundaries.
The relaxation zone method is the main approach to generating and absorbing waves.Relaxation zones are provided on both the left and right sides of the NWT.This method introduces the spatial weighting factor χ and the equation is as follows: ϕ can represent either the fluid velocity u or γ. γ represents the change in spatial variation of any fluid property, such as dynamic molecular viscosity, density and so on.
In the wave generation zone, χ gradually changes from 0 to 1 along the x-axis, and in the wave absorption zone, it gradually changes from 1 to 0. The above-mentioned wave generation and absorption are achieved by the open-source wave toolbox, waves2foam [30].

Computational Method
The finite volume method is a method used to discretise the computation domain into small cells in order to calculate and obtain the final numerical solution in the space domain.The finite volume method is widely used for solving fluids.Cardiff et al. developed a solver named solids4foam to solve the solid subdomain problem through FVM [29], which means that the calculation of FSI can be realised only by using FVM.Following the completion of the solution process, time and space will be discretised.In space, the computing domain is discretised into hexahedral elements in the form of mesh.In Openfoam, the front and back patches of a 2D model do not impose any boundary conditions, so they can be considered as "empty".Therefore, hexahedral elements actually solve problems in 2D.In time, continuous time will be divided by a certain time step.The time step depends on the Courant number (Co) and it can be obtained using the following equation: where ∆t is the step size, u is the normal velocity, and ∆x is the distance between the centres of two adjacent cells.∆t is set to 0.001 s in this NWT, and the Co will be significantly less than 1.

Fluid Solution
The Navier-Stokes (N-S) equations, based on conservation of mass and momentum, is the most commonly used governing equation in CFD.Because of different assumptions about fluids, the N-S equations are expressed differently depending on the situation.

∇ • v =0
(3) Incompressible, isothermal, Newtonian flow can be described by Equation (3) and Equation ( 4), where ∇ is the Hamiltonian, v is the velocity vector, P is pressure, ρ is the density, and τ is the viscous stress in which µ is the dynamic viscosity, and g is the gravitational acceleration.τ can be described by the following equation: The interface between water and air is the free surface, and it can be captured by the VOF method.A passive scalar field α is introduced to determine whether a cell should be water or air by describing the volume fraction of air and water in the cell.When α = 1 for this cell, the fluid inside is water, and vice versa, it is air.When α is between 0 and 1, it means that the cell contains a mixture of water and air.The nature of the VOF method could cause the free surface to blur, which requires the mesh refinement near the free surface area; this could also be addressed by using isoAdvector or the ghost fluid method [33].

Structural Solution
The St. Venant-Kirchhoff constitutive material model is widely used in engineering, which can be combined with the law of momentum conservation to simulate the solid subdomain.The following description of the governing equations was originally developed by Tuković et al. [34].
Assuming that the solid is elastic and compressible, the law of conservation of momentum can be used to obtain the following equation: ρ solid is the density of solid material, u is the displacement vector, F = I + (∇u)T is the deformation gradient tensor, I is the second-order identity tensor, and Σ is the second Piola-Kirchhoff stress tensor, which is related to the Cauchy stress tensor σ by the following expression: The stress/strain relationship is dictated as follows: where the Green-Lagrange strain tensor is as follows: G and Λ are Lamé's coefficients and they are related to Young's modulus E and Poisson's ratio v.
2.1.5.Fluid-Structure Interaction Tuković et al. [34] proposed a partitioned FSI solver to couple the fluid and solid subdomains, which has been implemented in the OpenFOAM framework.Information about the fluid subdomains and solid subdomains can be extracted for the fluid-solid interfaces.First, the velocity and pressure of the fluid in the fluid subdomain are obtained, and the load is applied to the fluid-solid interface.Then, the displacement in the solid domain is solved, and the moving speed of the fluid-solid interface is obtained.When this process meets the convergence condition, it enters the next time step, and if not, the iteration continues.The whole calculation process is shown in Figure 3.The matching of kinematic and dynamic conditions achieved in fluid-solid interfaces follows the following equations: The matching of kinematic and dynamic conditions achieved in fluid-solid interfaces follows the following equations: where v is the velocity vector, u is the displacement vector, n is the normal vector at the interface, and σ is the stress tensor.

Data Acquisition and Processing
The displacement data were used to analyse the efficiency of the energy conversion, and the vertical displacement data of the flexible structure could be obtained with a probe in OpenFOAM.The probes are arranged at a distance of 0.2*L, 0.4*L, 0.6*L, 0.8*L and 1*L from the left edge of the flexible structure in the horizontal direction, and were deployed at the middle position of the flexible structure in the vertical direction.In Figure 4, the positions of the probes are marked with a red dot, where L and h represent the length and thickness of the flexible structure respectively.
where v is the velocity vector, u is the displacement vector, n is the normal vector at the interface, and σ is the stress tensor.

Data Acquisition and Processing
The displacement data were used to analyse the efficiency of the energy conversion, and the vertical displacement data of the flexible structure could be obtained with a probe in OpenFOAM.The probes are arranged at a distance of 0.2*L, 0.4*L, 0.6*L, 0.8*L and 1*L from the left edge of the flexible structure in the horizontal direction, and were deployed at the middle position of the flexible structure in the vertical direction.In Figure 4, the positions of the probes are marked with a red dot, where L and h represent the length and thickness of the flexible structure respectively.In this research, the motion pattern of flexible structures in most cases are not sinusoidal.Therefore, it is necessary to obtain its amplitude in the frequency domain first by using the Fast Fourier Transform (FFT).
Figure 5 shows the displacement data obtained by the above probe located at point E (sea state: wave period 1.1 s, wave height 0.04 m, wave depth 1.5 m; deployment conditions of the flexible structure: L = 1.5112, d = 0.18) and the time required from wave action on the flexible structure to the full development of vibration is about 7 s.Therefore, it would be perfectly feasible to take the data after 9.5 s when analysing the amplitude.
In order to obtain the vibration amplitude of flexible structures more quickly and accurately, the Fast Fourier Transform (FFT) is applied to data processing.After FFT processing, the amplitude-frequency graph as shown in Figure 6 below can be obtained.Its X-coordinate has been normalised to the ratio of the wave frequency.From this figure, four peaks can be clearly discovered, which are roughly distributed in the position of f/f0 of 0.33, 1 and 2. When f/f0 is 1, the amplitude obtained by FFT is 0.0196 m, which is very close to the value of subtracting two peaks in the time domain, which also shows the reliability of this method.In this research, the motion pattern of flexible structures in most cases are not sinusoidal.Therefore, it is necessary to obtain its amplitude in the frequency domain first by using the Fast Fourier Transform (FFT).
Figure 5 shows the displacement data obtained by the above probe located at point E (sea state: wave period 1.1 s, wave height 0.04 m, wave depth 1.5 m; deployment conditions of the flexible structure: L = 1.5112, d = 0.18) and the time required from wave action on the flexible structure to the full development of vibration is about 7 s.Therefore, it would be perfectly feasible to take the data after 9.5 s when analysing the amplitude.The deformable WECs apply Dielectric Elastomer Generators (DEGs) to fo power take-off (PTO) system that generates electricity from structural deformatio this work, the energy harnessed was calculated as the deformation energy per unit of the device.The energy efficiency of the WEC was calculated by per unit width o In order to obtain the vibration amplitude of flexible structures more quickly and accurately, the Fast Fourier Transform (FFT) is applied to data processing.After FFT processing, the amplitude-frequency graph as shown in Figure 6 below can be obtained.Its X-coordinate has been normalised to the ratio of the wave frequency.From this figure, four peaks can be clearly discovered, which are roughly distributed in the position of f /f 0 of 0.33, 1 and 2. When f /f 0 is 1, the amplitude obtained by FFT is 0.0196 m, which is very close to the value of subtracting two peaks in the time domain, which also shows the reliability of this method.The deformable WECs apply Dielectric Elastomer Generators (DEGs) to fo power take-off (PTO) system that generates electricity from structural deformation this work, the energy harnessed was calculated as the deformation energy per unit w of the device.The energy efficiency of the WEC was calculated by per unit width o nessed energy divided by per unit width of wave energy, using the following equat The energy conversion efficiency of the device is obtained by averaging the φF tained at each probe point.In addition, f/f0 of 0.33 should not be involved in The deformable WECs apply Dielectric Elastomer Generators (DEGs) to form a power take-off (PTO) system that generates electricity from structural deformations.In this work, the energy harnessed was calculated as the deformation energy per unit width of the device.The energy efficiency of the WEC was calculated by per unit width of harnessed energy divided by per unit width of wave energy, using the following equation: The energy conversion efficiency of the device is obtained by averaging the φ FFT obtained at each probe point.In addition, f /f 0 of 0.33 should not be involved in the calculation.This is an unstable vibration.When the calculation time is extended, the vibration at this frequency will disappear.

Verifications
In this section, the rationality and accuracy of the above model are verified following four steps: (a) determine the key parameters of the model; (b) certify fluid solution model simplification to laminar flow is reasonable; (c) determine if the mesh setup can ensure the accuracy of the wave generator and solid sub-domain displacement.

Fluid and Structural Parameters
Tables 1-3 show the parameters required for the calculation.The conditions of wave height and wavelength meet the requirements of Stokes second-order wave theory, and the values of H/gT 2 and d/gT 2 are 1.69 × 10 −3 and 4.25 × 10 −2 , respectively.The parameters of solids are similar to those of rubber because most of the current WECs based on flexible structures are made of rubber.The model-scale wave condition was chosen based on previous validation against experiments, including (a) wave transmission and reflection; (b) structural deformations in waves, which was reported in detail in [26].However, the results will also be applicable on a larger scale following Froude scaling.Although the scaling will involve a mismatch in Reynolds numbers, the influence is expected to be minimal as inertial force is much more significant than viscous force in this application.

Laminar Flow
Compared with laminar flow, the turbulence model is suitable for cases where inertia is dominant in fluid flow.Due to the relatively complex computational mechanism of turbulence, many scholars have proposed a lot of different computational models.Therefore, for one new case, whether to adopt a turbulence model and what kind of turbulence model to adopt should be considered before running the simulation program.
The identification of which model should be used depends on a dimensionless quantity called the Reynolds number.The flexible structure is arranged in the ocean in a way similar to a swept plate.Generally, for the "swept plate" problem, when the Reynolds number is less than 5 × 10 5 , it is considered that there is no turbulent component in the fluid.In other words, the laminar flow model can be used in this model to reduce the computational cost.It has become common sense that the critical Reynolds number for laminar flow to become turbulent is between 3.2 × 10 5 and 3 × 10 6 .The equation to calculate Reynolds number is as follows: Re is the Reynolds number, ρ is the fluid density, v is the characteristic velocity, d is the characteristic length, and µ is the dynamic viscosity.In the case of the flexible structure, the characteristic length can be thought of as the length of the plate.Near the plate, the velocity of the fluid will not exceed 0.18 m/s, so the characteristic velocity can be thought of as 0.2 m/s.Even with relatively loose calculation parameters (Re = 2.7 × 10 5 ), it is also difficult to reach the Reynolds number of 3.2 × 10 5 .Therefore, it can be considered that the laminar flow model is consistent with the actual situation.In addition, the Keulegan-Carpenter number is around 0.2, which is smaller than 1, indicating that the flow does not separate from the plate [36]; thus, no turbulence effect is expected in the present study.However, we note that this is an ideal case and turbulence can be ubiquitous in real oceans.

Mesh Convergence Study for Wave Modelling
The accuracy of wave modelling is closely related to the number of cells, and it is common practice to obtain the optimal number of cells by grid sensitivity analysis.Ref. [37] quantified the relationship between waves and cells using cells per wavelength (CPW) and cells per wave height (CPH).The discretisation error can be minimised when CPW and CPH reach 50 and 20, respectively.
In this study, CPW was set at 75, which was in line with the value recommended by the above literature.In order to explore the influence of different CPH on the calculation accuracy of liquid level height, the liquid level height located downstream was collected.The Root-Mean Squared Error (RMSE) calculated according to Equation ( 2) is shown in Table 4. Increasing CPH from 3 to 12 had no significant impact on the calculation results.The results shown in Figure 7 coincide with the theoretical data when CPH is set to 6 and 12.The accuracy of wave modelling is closely related to the number of cells, and common practice to obtain the optimal number of cells by grid sensitivity analysis [37] quantified the relationship between waves and cells using cells per wavele (CPW) and cells per wave height (CPH).The discretisation error can be minimised w CPW and CPH reach 50 and 20, respectively.
In this study, CPW was set at 75, which was in line with the value recommende the above literature.In order to explore the influence of different CPH on the calcul accuracy of liquid level height, the liquid level height located downstream was colle The Root-Mean Squared Error (RMSE) calculated according to Equation ( 2) is show Table 4. Increasing CPH from 3 to 12 had no significant impact on the calculation re The results shown in Figure 7 coincide with the theoretical data when CPH is set to 6 12.

Mesh Convergence Study for Structural Modelling
The solid subdomains were spatially dispersed according to the four partition specifications of 38 × 3, 38 × 5, 76 × 5 and 76 × 8, respectively, in order to find the appropriate mesh parameters to avoid their great influence on the calculation results.In these four different specifications, the number of cells is 114,190,380 and 608.As shown in Table 5, there is little difference in the energy conversion efficiency obtained by different discrete settings, and the difference between the maximum and minimum values is only 0.75%.Since the flow of liquid from the surface to the seabed tends to decrease gradually, the vertical mesh density can reduce the computational burden by decreasing gradually.The details of cells are shown in Figure 8.

Energy conversion efficiency
137.93% 138.49% 137.77% 138.52% Since the flow of liquid from the surface to the seabed tends to decrease gradually, the vertical mesh density can reduce the computational burden by decreasing gradually.The details of cells are shown in Figure 8.

Results and Discussion
Existing studies have shown that the efficiency of energy conversion is mainly related to the depth and length of the flexible structure [38].In order to facilitate the development of the research, these two studied parameters were normalised.The length of the flexible structure and its depth in water are divided by the wavelength and height of the wave, respectively.The specific parameters are shown in Table 6.
Since the length of the flexible structure is related to the wavelength, the optimal AR will be obtained by changing the thickness of the flexible structure.For ease of calculation, the flexible structure will be set to 1m in length and the d/H is set to 5. Figure 9 provides an example of the WEC deformation for one wave period.It can be seen that the CFD and CSM solutions evolve together in the time domain, which demonstrates the fully coupled capability of the present model.It shows that the WEC deforms following the dominating wave shape and the end of the plate swings with the movement of the water.At the same time, because of the presence of WEC, the outgoing wave changes significantly compared with the incident wave.In the subfigures, the water phase is beneath the free surface in red, while the air phase is above the free surface in blue.
It takes about 10 h to compute this example (13 wave cycles), using four performance cores of I7-12700K with a memory frequency of 5200 MHz.In particular, it should be pointed out that the model has very low memory requirements, and we only needed 4G RAM to meet the computing requirements.

Effect of Deployment Depth and Length of Flexible Structure on Energy Conversion Efficiency
Figure 10 shows the energy conversion efficiency of flexible structures at different deployment depths and different lengths.It can be seen in the figure that when L/λ = 0.25 and d/H = 3, the energy conversion efficiency of the flexible structure reaches the highest value of 138.45%.Therefore, it can be considered that the best energy conversion efficiency can be obtained when the total length of the flexible structure is one-quarter of the wavelength and the deployment depth is three times the wave amplitude.When the relationship between the total length of the flexible structure and the wavelength is L/λ = 0.50 and L/λ = 0.75, the energy conversion efficiency of the device can also achieve a good effect (at d/H = 5 and d/A = 10, the energy conversion efficiency reaches 124.66% and 135.83%, respectively).
In terms of deployment depth, it seems that a larger L/λ requires a greater deployment depth.Moreover, in the case of insufficient deployment depth, the energy conversion efficiency will decline rapidly and significantly.This phenomenon is like a "wall", and when the deployment depth passes this "wall", the energy conversion efficiency will be significantly improved.
The interesting thing is, that when the value of L/λ is near 0.25, 0.5 and 0.75, the influence of the deployment depth of the flexible structure on the energy conversion efficiency will be weakened, which is very conducive to the deployment and application of this device, and it is a very challenging work to accurately control the deployment depth of the flexible structure.be significantly improved.
The interesting thing is, that when the value of L/λ is near 0.25, 0.5 and 0.75, the influence of the deployment depth of the flexible structure on the energy conversion efficiency will be weakened, which is very conducive to the deployment and application of this device, and it is a very challenging work to accurately control the deployment depth of the flexible structure.

Relationship between AR and Energy Conversion Efficiency of Flexible Structure
Figure 11 shows the relationship between the AR and energy conversion efficiency under the sea state when the wave height, wave depth and period is 0.02 m, 1.5 m and 1.1 s.The flexible structure was deployed at a depth of 0.14 m below the water level, with a total length of 1m.The result shows that if the AR is too small, it will make the deformation of the flexible structure very difficult, and its motion pattern will be like that of a rigid body.As the AR decreases, the efficiency of energy conversion increases significantly.At the same time, with the increase in AR, the maximum stress of flexible structures increases first and then decreases slowly.The influence of the AR value on motion amplitude is also reflected in Figure 12, which is consistent with the conclusion of Figure 11.When the AR

Relationship between AR and Energy Conversion Efficiency of Flexible Structure
Figure 11 shows the relationship between the AR and energy conversion efficiency under the sea state when the wave height, wave depth and period is 0.02 m, 1.5 m and 1.1 s.The flexible structure was deployed at a depth of 0.14 m below the water level, with a total length of 1m.The result shows that if the AR is too small, it will make the deformation of the flexible structure very difficult, and its motion pattern will be like that of a rigid body.As the AR decreases, the efficiency of energy conversion increases significantly.At the same time, with the increase in AR, the maximum stress of flexible structures increases first and then decreases slowly.The influence of the AR value on motion amplitude is also reflected in Figure 12, which is consistent with the conclusion of Figure 11.When the AR is 100 and 150, the limit range of motion is significantly greater than when the AR is 50.The motion amplitude of AR100 and AR150 also does not seem to differ much, which is in line with the content of the energy conversion efficiency in Figure 11.
Oceans 2024, 5 14 is 100 and 150, the limit range of motion is significantly greater than when the AR is 50.The motion amplitude of AR100 and AR150 also does not seem to differ much, which is in line with the content of the energy conversion efficiency in Figure 11.The appearance of the maximum stress is related to the bending degree of the flexible structure, and the greater the bending degree, the greater the stress on the upper and lower planes.Figure 13 shows the stress distribution of a flexible structure under a wave load, when the flexible structure is subjected to its maximum stress at work.These findings support the design and operation of deformable WECs by enabling the prediction of cyclic loading for fatigue life as well as extreme loading for ultimate The appearance of the maximum stress is related to the bending degree of the flexible structure, and the greater the bending degree, the greater the stress on the upper and lower planes.Figure 13 shows the stress distribution of a flexible structure under a wave load, when the flexible structure is subjected to its maximum stress at work. 5 15 strength; this can further inform the inspection schedule and maintenance strategies [39][40][41].Thus, the present work may potentially improve the affordability and safety of the systems.When considering the energy efficiency and survivability of the device in relation to varying parameters such as device length and deployment depth, it should also be noted that costs are always an important trade-off.Therefore, the design should fully consider the deployment location's wave conditions and choose the most cost-effective approach.

Conclusions
In this paper, a fully coupled CFD+CSM method was used to simulate the interaction between a deformable structure with nonlinear waves.On this basis, the energy conversion efficiency of a submerged deformable wave energy converter was studied with varying parameters of total length, deployment depth and AR.It was found that the total length of the flexible structure plays a crucial role in the energy conversion efficiency and when the total length of the flexible structure is 0.25, 0.5 and 0.75 times the wavelength, optimal energy conversion efficiency can be obtained.At the same time, the optimal deployment depth corresponding to these three lengths also increases with the structural length.
When AR increases, the plate structure is more prone to deformation, so the energy conversion efficiency shows a gradually increasing trend, and this trend is no longer obvious after reaching a certain limit.However, the improvement of energy conversion efficiency does not lead to a continuous increase in the load burden of flexible structures.Its maximum stress will reach a maximum value and then slowly decrease.The present work These findings support the design and operation of deformable WECs by enabling the prediction of cyclic loading for fatigue life as well as extreme loading for ultimate strength; this can further inform the inspection schedule and maintenance strategies [39][40][41].Thus, the present work may potentially improve the affordability and safety of the systems.When considering the energy efficiency and survivability of the device in relation to varying parameters such as device length and deployment depth, it should also be noted that costs are always an important trade-off.Therefore, the design should fully consider the deployment location's wave conditions and choose the most cost-effective approach.

Conclusions
In this paper, a fully coupled CFD+CSM method was used to simulate the interaction between a deformable structure with nonlinear waves.On this basis, the energy conversion efficiency of a submerged deformable wave energy converter was studied with varying parameters of total length, deployment depth and AR.It was found that the total length of the flexible structure plays a crucial role in the energy conversion efficiency and when the total length of the flexible structure is 0.25, 0.5 and 0.75 times the wavelength, optimal energy conversion efficiency can be obtained.At the same time, the optimal deployment depth corresponding to these three lengths also increases with the structural length.
When AR increases, the plate structure is more prone to deformation, so the energy conversion efficiency shows a gradually increasing trend, and this trend is no longer obvious after reaching a certain limit.However, the improvement of energy conversion efficiency does not lead to a continuous increase in the load burden of flexible structures.Its maximum stress will reach a maximum value and then slowly decrease.The present work demonstrates a hydroelastic design methodology for a deformable wave energy converter, with a range of empirical rules that are useful for the industry.
Since the current CFD+CSM method has the computational power of two-way full coupling, the deformation of the structure and the change in the flow field can be clearly displayed, which can be used to assist the design work that is needed for marine structures with large deformations.Thus, this work tackles the technical challenge of modelling fully coupled fluid-structure interactions, which are essential for deformable WECs.Also, it facilitates the modelling of other offshore and coastal industries.
However, this model still has its limitations.The current simulation is performed in 2D, which excludes wave scattering near the geometry and potentially influences the energy efficiency in 3D.Advanced incorporation of DEGs as the PTO system should be developed in future work, e.g., mechanical damping; the present study also only simulated a single device and did not model the cases of array configurations, which are of interest [42][43][44].Although regular wave conditions were used to establish relationships between energy efficiency and incident wave conditions, it is noted that realistic sea waves will be irregular and optimal performance might be achieved through a real-time PTO controlling strategy [45].Based on the CFD+CSM approach, it is suggested in future work to explore the option of containing multiple materials in the CSM part [29]; the device geometry could be more complex and optimised in 3D [46], and mooring analysis could be added to enhance the practicalities [47].

Figure 2 .
Figure 2. Schematic of the case: a flat flexible structure arranged under water is subjected to an incident wave and produces elastic deformation.L and d are the total length and deployment depth of the device (All units of length are meters).

Figure 2 .
Figure 2. Schematic of the case: a flat flexible structure arranged under water is subjected to an incident wave and produces elastic deformation.L and d are the total length and deployment depth of the device (All units of length are meters).

Figure 3 .
Figure 3. Solution algorithm for strong two-way coupling (reset the figure in [35]).The matching of kinematic and dynamic conditions achieved in fluid-solid interfaces follows the following equations:

Figure 4 .
Figure 4.The location of displacement probes (A, B, C, D, E).

Figure 4 .
Figure 4.The location of displacement probes (A, B, C, D, E).

Oceans 2024, 5 Figure 5 .
Figure 5. Example of vertical vibration displacement at a point of the flexible structure (The d in the red box is used for analysis).

Figure 5 .
Figure 5. Example of vertical vibration displacement at a point of the flexible structure (The data in the red box is used for analysis).

Figure 5 .
Figure 5. Example of vertical vibration displacement at a point of the flexible structure (The d in the red box is used for analysis).

Figure 6 .
Figure 6.An example of FFT processing to show the displacement amplitude at different f/f0.

Figure 6 .
Figure 6.An example of FFT processing to show the displacement amplitude at different f /f 0 .

Figure 7 .
Figure 7. Liquid level data at CPH 6 and 12 (compared with data obtained through Stokes se order wave theory).

Figure 7 .
Figure 7. Liquid level data at CPH 6 and 12 (compared with data obtained through Stokes secondorder wave theory).

Figure 8 .
Figure 8. Combined fluid and solid mesh (red part for solid subdomain and blue part for fluid subdomain).

Figure 9 .
Figure 9. Periodic interaction between the flexible WEC and incoming waves (C063).The applied wave condition in this figure is T = 1.1 s and H = 0.03 m.In the subfigures, the water phase is beneath the free surface in red, while the air phase is above the free surface in blue.

Figure 10 .
Figure 10.Heat map of energy conversion efficiency.

Figure 10 .
Figure 10.Heat map of energy conversion efficiency.

Figure 11 .
Figure 11.The relationship between conversion efficiency, maximum stress and AR.

Figure 11 .
Figure 11.The relationship between conversion efficiency, maximum stress and AR.

Figure 11 .
Figure 11.The relationship between conversion efficiency, maximum stress and AR.

Figure 12 .
Figure 12.Upward limit (solid line) and downward limit (dash line) of the flexible structure's deformations with different ARs.

Figure 12 .
Figure 12.Upward limit (solid line) and downward limit (dash line) of the flexible structure's deformations with different ARs.

Table 1 .
Basic parameters of wave.

Table 2 .
Basic parameters of fluid.

Table 3 .
Basic parameters of solid.

Table 5 .
The energy conversion efficiency with different numbers of cells.