Numerical Study on Hydrodynamic Performance of a New Semi-Submersible Aquaculture Platform

: To address the limitations of traditional aquaculture cages, such as poor performance in offshore environments, challenges in energy supply, and the inability to integrate modern farming facilities, a new semi-submersible aquaculture platform has been developed and put into operation. The hydrodynamic performance of the platform is analyzed using a fully coupled numerical model. The numerical model incorporates the main structure of the platform, the net system, and the mooring system, and is based on the linear potential ﬂow theory and the Morison equation. The natural period of the platform is found to be far away from the common wave period. The mooring system for the platform in a speciﬁc sea area is proven to be safe. The net system signiﬁcantly affects the hydrodynamic performance and alleviates the dynamic response of the platform. This study provides valuable insights for researchers and engineers in the design and optimization of this type of semi-submersible aquaculture platform in the future.


Introduction
As the global population and economy continue to grow, human demand for aquaculture products increases rapidly, and natural fishery resources decline drastically.Aquaculture farms play an increasingly important role in avoiding overfishing of aquatic products.However, aquaculture farms in inshore and harbor areas are facing limited development and receiving increased criticism due to the lack of production sites and the adverse impact on local ecosystems [1].In contrast, moving aquaculture farms to the offshore seas is a popular and inevitable trend due to the broader development space, better water quality, and stronger environmental purification capacity in more exposed sea regions [2].Nevertheless, aquaculture facilities in open seas will endure much more severe environmental conditions.The hydrodynamic performance of aquaculture facilities is a crucial factor that affects their stability, safety, and economy.Evaluating the hydrodynamic performance of these structures has important theoretical and practical significance for their design, optimization, and management.Therefore, the design and development of aquaculture facilities with excellent hydrodynamic performance can significantly support the advancement of offshore aquaculture technology.
Surface floating gravity-type cages have been widely used in fish farming in the past decades due to their large breeding capacity and simple daily maintenance compared with other types of cages.Researchers have widely investigated the hydrodynamic characteristics of this kind of cage and its components, such as the floating collar, the net system, and the mooring system, through numerical simulations and physical models.Fredriksson et al. [3] studied the dynamic characteristics of a central spar fish cage deployed in the Gulf of Marine using physical and numerical models.Lader and Enerhaug [4] investigated the dependency between the forces and the geometry of a flexible circular net cage with different sizes of bottom weights in a flume tank by measuring the global forces and net deformation.Zhao et al. [5] developed a numerical model to analyze the motion response and the mooring forces of a gravity cage in a combined wave-current flow and found that the numerical results agreed well with the experimental data.Huang et al. [6] utilized a previously validated numerical model to examine the effects of waves and currents on the tension of mooring lines and volume reduction of two types of gravity net cages: a double-cage system and a single-cage system.Dong et al. [7] proposed a numerical model based on statistical approaches and Fourier transforms to analyze the motion response and mooring line tension of gravity cages in irregular waves.Li et al. [8] studied the floater motion and the net volume reduction of a gravity net cage in waves and currents using the finite element method and the buoyancy distribution method.Li et al. [9] investigated the nonlinear hydroelastic response of a deep-water gravity fish cage in irregular waves using the finite element method and the modal superposition method.Shainee et al. [10] proposed a self-submersible single-point-mooring cage concept and investigated its submergence characteristics in regular waves with following currents using the numerical model Aqua-FE.Shainee et al. [11] further investigated its submergence characteristics in random waves.Cifuentes and Kim [12] studied the hydrodynamic response of a cage system with a flexible surface collar, net, ballast, and mooring lines in regular/irregular waves and currents using a Morison-force model and validated it with experimental data.Shen et al. [13] investigated the mooring loads of a realistic aquaculture fish farm system with a floating collar, a flexible net cage, and a sinker tube using a reliable numerical tool and validated it with experimental data.Liu et al. [14] numerically studied the mooring forces of a fish cage array with six gravity cages in currents and waves using the finite element method.Liu et al. [15] introduced a new net cage with a shielding device and analyzed its hydrodynamic characteristics using a numerical method based on computational fluid dynamics.
The gravity cages with floating collars and flexible nets mentioned above suffer from large cage deformation at offshore sites with high sea states.It is urgent to develop aquaculture facilities in typical offshore environments to overcome this.Various innovative concepts of aquaculture facilities have been developed recently by combining technologies from the offshore industry with fish farming.Correspondingly, researchers have conducted a series of numerical and physical experimental studies to evaluate and enhance the hydrodynamic performance of these structures.
Ocean Farming AS, a subsidiary of the SalMar Group, developed the world's first offshore fish farm, Ocean Farm 1, which is dodecagonal in shape.It has a height of 68 m, a diameter of 110 m, and a volume of 250,000 m 3 [16].Dou [17] conducted hydrodynamic analysis of Ocean Farm 1 in both frequency and time domains using the HydroD V4.10 and SIMA V4.1 software packages, respectively.Jin et al. [18] modeled Ocean Farm 1 using SIMO/RIFLEX/SIMA V4.1 software packages and investigated its motion responses under different combined conditions of irregular waves and currents.Yu et al. [19] assessed the structural strength of Ocean Farm 1 under supply vessel collisions using the nonlinear finite element program USFOS.Zhao et al. [20] presented a semi-submersible offshore fish farm similar to Ocean Farm 1 and conducted a series of physical model experiments to investigate its motion response and mooring line tension.Liu et al. [21] analyzed the hydrodynamic characteristics of an aquaculture platform similar to Ocean Farm 1 in uniform currents using the porous media combined with a rigid wall based on the SST k-omega turbulence model.
NSK Ship Design for a Norwegian fishing company designed a novel vessel-shaped fish farm, Havfarm 1, which consists of a vessel-shaped floater, five flexible cages under the floater, and a single point mooring system, with a total length of 385 m [22].Li et al. [23,24] analyzed the global responses of Havfarm 1 and the tensions in the mooring lines under various wave and current conditions using the coupled SIMO-RIFLEX V4.1 computer program.Li et al. [25] proposed an integrated method to predict the misalignment angle between the heading of the vessel-shaped offshore fish farm and the incoming currents under combined waves and currents, to facilitate the design of the dynamic positioning system used for the fish farm.Li et al. [26] proposed an optimization methodology for mooring system design and designed the mooring system of the vessel-shaped offshore fish farm accordingly.Ma et al. [27] presented a single-point moored vessel-shaped floating aquaculture platform with two different steel floating frames, nets, and weight system pipes and analyzed its motion, mooring force, and net reaction force in regular waves using the boundary element method and the lumped mass method.Ma et al. [28] depicted a hinged multi-body floating aquaculture platform with square floating frames, nets, and weight system pipes and examined the effects of the hinge joint rotational stiffness on its dynamic responses under different wave steepnesses using the boundary element method and the lumped mass method.
Ocean University of China and Hubei Ocean Engineering Equipment Research Institute co-designed China's first fully submersible deep-sea fish farm, ShenLan 1.The farm features a rigid frame with eight vertical columns, two cylindrical liquid-gas tanks, and ultra-high-strength polyethylene square nets.The full-scale diameter and height are 60.44 m and 34.45 m, respectively [29].Wang et al. [30,31] used the open-source CFD toolbox REEF3D to analyze the effects of the wave parameters and structural variations on the response motions and mooring forces of ShenLan 1. Miao et al. [32] presented a semi-submersible offshore fish farm similar to ShenLan 1 with a pontoon, main frame system, net system, and upper building and developed a hybrid method that combined the direct time domain  [33] investigated the mooring forces and motion response of Dehai Farm 1 in waves and currents through a series of physical model tests.Chu and Wang [34] presented a novel offshore fish farm consisting of a deep draught floating spar and an octagonal net steel cage (named COSPAR) and performed the hydrodynamic response analysis of free-floating COSPAR and the coupled analysis involving COSPAR and mooring lines using the software Design Modeler R2012 and Aqwa R2017 owned by Ansys, Inc. (Canonsburg, PA, USA) [35,36].
Guangzhou Institute of Energy Conversion (GIEC) of the Chinese Academy of Sciences designed a semi-submersible aquaculture platform by combining Sharp Eagle WEC technology and offshore aqua-culture.This innovative design led to the successful deployment of the world's first semi-submersible wave energy aquaculture platform named "Penghu" (Figure 1) in Zhuhai, China, in 2019 [37].The invention belongs to the team of the authors of this paper.By utilizing this platform, the challenges encountered by traditional aquaculture cages in offshore environments, including poor performance, limited energy supply, and inability to accommodate modern farming equipment, can be effectively addressed.
GIEC continues to develop a series of semi-submersible aquaculture platforms to further promote industrial applications.Systematic numerical or experimental analyses are essential for the design of these offshore structures to assess their performance in complicated environmental conditions.However, open publications on the detailed hydrodynamic and structural analysis of these concepts are rare, to the authors' knowledge.
In this paper, we aim to fill the gap in the understanding of the hydrodynamic performance of a semi-submersible aquaculture platform, which is newly designed by GIEC and put into operation.While Yue et al. [38] developed a numerical model to analyze the structural response of this type of platform, their model did not consider the mooring system and nets and only evaluated the platform under pure wave action.In contrast, our study proposes a comprehensive numerical model that takes into account the semi-submersible main structure, net system, and mooring system and further analyzes the mooring line tension and motion response of the platform under wave, wind, and current conditions.
In this paper, we aim to fill the gap in the understanding of the hydrody formance of a semi-submersible aquaculture platform, which is newly design and put into operation.While Yue et al. [38] developed a numerical model to structural response of this type of platform, their model did not consider th system and nets and only evaluated the platform under pure wave action.In c study proposes a comprehensive numerical model that takes into account th mersible main structure, net system, and mooring system and further analyze ing line tension and motion response of the platform under wave, wind, and c ditions.
The rest of this paper is organized as follows: Section 2 provides a descri new aquaculture platform.Section 3 presents a coupled numerical model of th under the combined action of waves, winds, and currents.Section 4 validates ical results and analyzes the hydrodynamic performance of the platform in reg and the dynamic responses of the platform with and without nets in waves, currents.Section 5 summarizes the main conclusions and the prospects for fut  The rest of this paper is organized as follows: Section 2 provides a description of the new aquaculture platform.Section 3 presents a coupled numerical model of the platform under the combined action of waves, winds, and currents.Section 4 validates the numerical results and analyzes the hydrodynamic performance of the platform in regular waves and the dynamic responses of the platform with and without nets in waves, winds, and currents.Section 5 summarizes the main conclusions and the prospects for future study.In this paper, we aim to fill the gap in the understanding of the hydrodynam formance of a semi-submersible aquaculture platform, which is newly designed b and put into operation.While Yue et al. [38] developed a numerical model to ana structural response of this type of platform, their model did not consider the m system and nets and only evaluated the platform under pure wave action.In cont study proposes a comprehensive numerical model that takes into account the se mersible main structure, net system, and mooring system and further analyzes th ing line tension and motion response of the platform under wave, wind, and curr ditions.

Description of the Platform
The rest of this paper is organized as follows: Section 2 provides a descriptio new aquaculture platform.Section 3 presents a coupled numerical model of the p under the combined action of waves, winds, and currents.Section 4 validates the ical results and analyzes the hydrodynamic performance of the platform in regula and the dynamic responses of the platform with and without nets in waves, win currents.Section 5 summarizes the main conclusions and the prospects for future  The platform is equipped with seawater desalination and sewage treatmen storage space, restaurants, and other facilities.These additions transform it into a marine industrial integration center that encompasses green intelligent equipmen culture, fishing, sightseeing, and tourism.An intelligent fishery breeding system integrated into the platform to enable modern fishery production functions, in The platform is equipped with seawater desalination and sewage treatment plants, storage space, restaurants, and other facilities.These additions transform it into a modern marine industrial integration center that encompasses green intelligent equipment, aquaculture, fishing, sightseeing, and tourism.An intelligent fishery breeding system can be integrated into the platform to enable modern fishery production functions, including automatic feeding, fish swarm monitoring, and water quality monitoring.Furthermore, the platform is equipped with photovoltaic power generation equipment, utilizing clean energy and ensuring its green and environmentally friendly operation.With a generation and storage capacity of 100 kilowatts, the platform achieves energy self-sufficiency [39].Shape-Eagle wave energy converters are not installed on Pu Sheng 1 due to economic considerations.A more modern and luxurious agriculture and fishing tourism platform named Mintou No.1 with Shape-Eagle wave energy converters had a soft opening in 2023 [40].

Main Parameters of the Platform
The platform has a total length of 86 m, a width of 30 m, and a height of 18 m.It consists of a steel semi-submersible main structure and soft nets that form a breeding water volume of approximately 30,000 cubic meters.Both the operational and survival drafts are 10 m.Table 1 lists the detailed parameters of the platform.In the hydrodynamic analysis of the platform, a three-dimensional Cartesian coordinate system is adopted, with its origin located on the mean water surface and intersecting with the middle longitudinal and transverse sections.The x-axis points forward to the bow of the platform, the y-axis points forward to the port side of the platform, and the z-axis points upward vertically.The mass properties of the platform are presented in Table 2.The hydrodynamic model of the platform is depicted in Figure 3.Only the structure below the waterline is considered in the hydrodynamic calculation.

Environment Parameters
The platform has been in operation at the Longqi Bay National Intelligent Marine Ranch, located in Ledong County in southwestern Hainan since 24 March 2022 (see Figure

Environment Parameters
The platform has been in operation at the Longqi Bay National Intelligent Marine Ranch, located in Ledong County in southwestern Hainan since 24 March 2022 (see Figure 4).

Environment Parameters
The platform has been in operation at the Longqi Bay National Intelligent Ranch, located in Ledong County in southwestern Hainan since 24 March 2022 (se 4).As shown in Figure 5, the longitude and latitude coordinates of the platform tion are 18°19′29.839″N and 108°50′13.633″E, respectively.Table 3 presents the environmental conditions with a return period of 20 yea project site.The site has a mean water depth of 20 m.The wave is expressed

Environment Parameters
The platform has been in operation at the Longqi Bay National Intelligent Ranch, located in Ledong County in southwestern Hainan since 24 March 2022 (se 4).As shown in Figure 5, the longitude and latitude coordinates of the platform tion are 18°19′29.839″N and 108°50′13.633″E, respectively.Table 3 presents the environmental conditions with a return period of 20 yea project site.The site has a mean water depth of 20 m.The wave is expressed Table 3 presents the environmental conditions with a return period of 20 years at the project site.The site has a mean water depth of 20 m.The wave is expressed by the JONSWAP spectrum.The wind speed is referenced at 10 m above the mean sea level, and the variation of the current velocity with water depth is neglected.The wind is expressed by the NPD spectrum, and the current is assumed to be constant.

Numerical Model
The motion equation for the time domain coupling analysis of the platform is expressed as follows: where ..
x, and x represent the acceleration, velocity, and displacement of the platform, respectively; t is the time; M is the mass matrix; µ ∞ is the infinite-frequency added-mass matrix; h(t − τ) is the wave-radiation-retardation kernel matrix; B is the additional damping matrix; C is the hydrostatic stiffness matrix; F Wave is the total wave excitation force on the platform; F Mooring is the mooring system force on the platform; F Current is the current force on the platform; F Wind is the wind force on the platform; and F Net is the net system force on the platform.

Potential Flow Model
For the main structure of the platform, the linear wave theory is adopted.The fluid is assumed to be inviscid and the motion is irrotational.The flow field's total velocity potential ϕ consists of the incident wave velocity potential ϕ i , the diffraction potential ϕ d , and the radiation potential ϕ r .
The governing equations and boundary conditions are given as follows [21,28]: Laplace equation: Linearized free surface condition: Structure surface condition: Sea bed condition: ∇ϕ → 0 Infinity condition: where n is the unit normal vector; ω is the wave frequency; g is the gravitational acceleration; and k is the wave number.
The boundary element method is used to calculate wave loads on the main structure.Based on the linearized Bernoulli equation, the distribution of hydrodynamic pressure p is calculated by: where ρ is the density of water with ρ = 1025 kg/m 3 .The forces and moments of the six-degree-of-freedom can be calculated by integrating the wave pressure acting on the outer surface of the structure: The first-order equation of motion for the structure in the frequency domain can be established as follows: where m ij is the generalized mass matrix; c ij is the generalized hydrostatic restoring stiffness; µ ij is the frequency-dependent added mass; λ ij is the frequency-dependent radiation damping; b ij is the additional damping; f f k i is the Froude-Krylov force; and This study considers only the linear additional damping resulting from the fluid viscosity.Due to the lack of relevant model test data for the platform, the additional linear damping in heave, roll, and pitch directions are assumed to be 5%, 8%, and 8% of the critical damping, respectively.
The added mass matrix, radiation damping matrix, and the first-order wave force are obtained using the three-dimensional Green function method, which is based on the wave diffraction and radiation theory.The motion response characteristics can be quantified by normalizing the response amplitude to the input amplitude.This normalization process yields the response amplitude operator (RAO).The RAOs for the motion of the structure are determined by solving the frequency domain motion Equation (10).
In addition to the first-order wave force f the structure is also subjected to the second-order wave force, which includes the mean wave drift force, difference-frequency wave force, and sum-frequency wave force.In this study, the second-order wave force is calculated by the middle-field method proposed by Chen [41].
The hydrodynamic coefficients, including the added mass, the radiation damping, the first-order wave force, and the second-order wave force, are then converted into the motion Equation (1) to conduct the time domain analysis, which integrates the main structure, mooring system, and net system.

Wind Force
In accordance with the rules of CCS [42], the wind force F Wind acting on the platform should be calculated using the following formula: where P represents the wind pressure with P = 0.613 × 10 −3 U 2 kPa and U represents the wind speed; S denotes the orthographic projection area of the wind-exposed component on the platform, which is 630 m 2 along the x-axis direction (longitudinal) and 271.5 m 2 along the y-axis direction (transverse) in this study; C h stands for the height coefficient and is 1 in this study, while C s represents the shape coefficient and is also 1 in this study.

Current Force
In accordance with the rules of CCS [42], the current force F Current acting on the platform should be calculated using the following formula: where C * D is the drag coefficient with C * D = 1.5 in this study; V * denotes the current velocity; and A * stands for the projected area of the component on the plane perpendicular to the current velocity, which is 609 m 2 along the x-axis direction (longitudinal) and 219 m 2 along the y-axis direction (transverse) in this study.

Net System
Empirical formulas are usually employed for the rapid prediction of hydrodynamic forces on the nets to calculate the loads.The Morison equation is used to compute the hydrodynamic loads on each unit, and the loads on the nets are obtained by summing up the calculated forces on all units.The equation used for this calculation is presented below: where C M is the inertia coefficient; V is the volume of water displacement; C D represents the drag coefficient; u is the flow velocity; . u is the flow acceleration; and A represents the cross-sectional area of the structure perpendicular to the flow direction.

Hydrodynamic Coefficients
The nets of the project are made of nylon cables.As shown in Table 4 and Figure 6, the diameter of a single nylon cable is 5 mm, and the density is 6.5 × 10 −6 t/m.The length of the nylon cable is the same as the distance between the two ends of the opening area of the platform.The spacing between two adjacent nylon cables is 50 mm.The forces on the nets are calculated using the Morison Equation (13).According to the guideline issued by the CCS [43], the drag coefficient in the Morison equation can be estimated using the following formula: where ' 90  =− ,  represents the angle of attack that is the angle between the net and the flow direction on the horizontal plane, and when the flow direction is perpendicular to the nets, the angle of attack is 0; n S represents the solidity ratio of nets, which is the ratio of the project area of the nets to the whole contour area(Figure 6).In this study, the solidity ratio of nets is Sn = 0.22, and the corresponding drag coefficient of the whole net is =0.305The forces on the nets are calculated using the Morison Equation (13).According to the guideline issued by the CCS [43], the drag coefficient in the Morison equation can be estimated using the following formula: where α = 90 − α, α represents the angle of attack that is the angle between the net and the flow direction on the horizontal plane, and when the flow direction is perpendicular to the nets, the angle of attack is 0; S n represents the solidity ratio of nets, which is the ratio of the project area of the nets to the whole contour area (Figure 6).
In this study, the solidity ratio of nets is Sn = 0.22, and the corresponding drag coefficient of the whole net is C d = 0.305.The inertia coefficient is C m = 1.0.According to the principle of force equivalence, the drag coefficient of the single net line is C d = 1.525, and the inertia coefficient of the single net line is C m = 1.0.

Equivalent Simulation
As shown in Figure 7, the entire nets are divided into multiple subparts in the vertical direction, and the cables within each subpart are assumed to form a new single cable.For instance, the cables within the dashed box are replaced by a new cable positioned in the middle of the dashed box.The equivalent substitution follows two principles: (1) The weight of the new cable equals the total weight of the original cables; (2) The hydrodynamic force of the new cable equals the total hydrodynamic forces of the original cables.The same approach is used for equivalent substitution in the horizontal direction [44].
where d represents the diameter of a single cable; and n is the number of cab form the equivalent cable.In this study, 95 cables in the vertical direction are consolidated into 1 cable, and 60 cables in the horizontal direction are consolidated into 1 equiv The relevant equivalent coefficients of nets based on Equations ( 15)-( 17) are p Table 5.Based on the Morison equation, the equivalent diameter D, equivalent drag coefficient C D , and equivalent inertia coefficient C M of the new cable are calculated as follows [44,45]: where d represents the diameter of a single cable; and n is the number of cables used to form the equivalent cable.
In this study, 95 cables in the vertical direction are consolidated into 1 equivalent cable, and 60 cables in the horizontal direction are consolidated into 1 equivalent cable.The relevant equivalent coefficients of nets based on Equations ( 15)-( 17) are presented in Table 5.

Mooring System
The mooring lines are simulated using the finite element model and the lumped-mass method.Each line is divided into several line segments with a node at each end.The segments only model the axial and torsional behaviors of the line.Other properties (mass, weight, buoyancy, etc.) are lumped at the nodes [46].The hydrodynamic loads on the mooring lines are calculated using Morison's equation.

Composition of the Mooring System
Previous studies [23][24][25] have predominantly focused on single-point mooring systems in the context of vessel-shaped fish farms.While the single-point mooring system provides good environment direction-following capability, it requires a large sea area due to extensive platform movements.In order to reduce the range of motion, this paper adopts a multipoint mooring system for the presented semi-submersible platform.As Figure 8 shows, the system arranges eight mooring lines symmetrically into four groups, with two lines in each group.At the project site, the angle between the geographical north direction, and the position direction of the X-axis is 135 • .Line 1 makes an angle of 45 • with the positive X-axis direction, while line 2 makes an angle of 37.5 • .The horizontal distance between the fairlead and the anchor measures 240 m. l.Sci.2023, 13, x FOR PEER REVIEW to extensive platform movements.In order to reduce the range of adopts a multi-point mooring system for the presented semi-subme Figure 8 shows, the system arranges eight mooring lines symmetrical with two lines in each group.At the project site, the angle between the direction, and the position direction of the X-axis is 135°.Line 1 ma with the positive X-axis direction, while line 2 makes an angle of 37 distance between the fairlead and the anchor measures 240 m.Eight identical mooring lines are utilized.The normal drag, axial f and axial inertia coefficients of each line are 2.6, 1.4, 1.0, and 0.5, resp tailed parameters are presented in Table 6.Eight identical mooring lines are utilized.The normal drag, axial force, normal inertia, and axial inertia coefficients of each line are 2.6, 1.4, 1.0, and 0.5, respectively.Other detailed parameters are presented in Table 6.Each mooring line connects one end to the fairlead on the platform and the other end to the anchor on the seabed.The coordinates of fairleads and anchors are provided in Table 7.The fairlead coordinates are based on the local Cartesian coordinate system, as described in Section 2.1.The anchor coordinates are based on the global Cartesian coordinate system, with the origin located at the mean water level.The X-axis points forward to the bow of the platform, the Y-axis points toward the port of the platform, and the Z-axis points vertically upward.

Loading Cases for Mooring Analysis
In the mooring analysis, the wave direction is defined as follows: 0 • represents the forward propagation of the wave along the x-axis of the platform's coordinate system, 90 • represents the forward propagation of the wave along the y-axis of the platform's coordinate system, and other wave directions follow the same principle.The direction definition for winds and currents is identical to that of waves, as depicted in Figure 9.
Appl.Sci.2023, 13 In the mooring analysis, the wave direction is defined as follows: 0° represe forward propagation of the wave along the x-axis of the platform's coordinate syste represents the forward propagation of the wave along the y-axis of the platform's c nate system, and other wave directions follow the same principle.The direction defi for winds and currents is identical to that of waves, as depicted in Figure 9. Based on the environmental parameters in Section 2.2 and the symmetric ar ment of the platform and mooring system, Table 8 presents the loading cases for su conditions, following the CCS rules [42].As the environmental conditions are more Based on the environmental parameters in Section 2.2 and the symmetric arrangement of the platform and mooring system, Table 8 presents the loading cases for survival conditions, following the CCS rules [42].As the environmental conditions are more severe for survival conditions than for operation conditions, this study only considers loading cases for survival conditions.

Verification of Numerical Results
Hydrodynamic coefficients, including the added mass, damping coefficient, the firstorder wave force, and the second-order wave force, are obtained using the software package HydroStar v8.00 [47].This software is based on the three-dimensional Green function potential flow theory in a finite water depth.
The project site has a water depth of only 20 m, which classifies it as a shallowwater area.Therefore, the effects of the second-order differential frequency forces become important.It is necessary to calculate the full quadratic transfer function (QTF).
In this study, the wave frequency ranges from 0.05 rad/s to 2.0 rad/s, with an interval of 0.05 rad/s, resulting in a total of 40 wave frequencies.The wave frequency range covers the most common wave frequencies observed in our study area, and the interval ensures a sufficient resolution for the analysis.The wave direction ranges from 0 • to 345 • , with an interval of 15 • , giving a total of 24 wave directions.The wave direction range represents the full circle of possible wave directions, and the interval is consistent with the wind direction data we collected from the meteorological station.Basically, the same selection of the range and interval of the wave frequency and wave direction for the common frequency domain analysis can be found in [38].
The mesh convergence analysis is conducted to examine the effects of mesh size on the numerical results of hydrodynamic coefficients.Three different mesh sizes of 1.0 m, 1.5 m, and 2.0 m are selected for the analysis, with respective total mesh numbers of 4846, 2221, and 1022, respectively.A quadrilateral mesh is utilized, and the mesh size is decreased in areas with complex surface shapes.Figures 10-12 present the RAO results of surge, heave, and pitch under three different sizes.Based on these results, a mesh size of 1.0 m is adopted to balance the accuracy and efficiency of the numerical results in this study.The mesh model of the platform is illustrated in Figure 13.
the numerical results of hydrodynamic coefficients.Three different mesh sizes of 1.0 m, 1.5 m, and 2.0 m are selected for the analysis, with respective total mesh numbers of 4846, 2221, and 1022, respectively.A quadrilateral mesh is utilized, and the mesh size is decreased in areas with complex surface shapes.Figures 10-12 present the RAO results of surge, heave, and pitch under three different sizes.Based on these results, a mesh size of 1.0 m is adopted to balance the accuracy and efficiency of the numerical results in this study.The mesh model of the platform is illustrated in Figure 13.The middle-field method, proposed by Chen [41], has been implemented in the software HydroStar v8.00 [47].The mean drift forces are calculated using both the far-field method and the middle-field method.Figures 14-16 compare the mean drift forces of The middle-field method, proposed by Chen [41], has been implemented in the software HydroStar v8.00 [47].The mean drift forces are calculated using both the far-field method and the middle-field method.Figures 14-16 compare the mean drift forces of surge, sway, and yaw under different wave incident angles obtained using these two methods.The results indicate that the mean drift forces calculated using the middle-field method agree with those obtained using the far-field method, demonstrating that the middle-field method is a reliable approach to calculating the full quadratic transfer functions in the analysis.The maximum value of the mean drift force in the sway direction is significantly larger than that in the surge direction.This can be attributed to the larger area subjected to forces in the sway direction compared with the surge area.The peak values of the mean drift force in the sway direction occur at approximately 1.6 rad/s (corresponding to a wave period of 3.93 s), 1.15 rad/s (corresponding to a wave period of 5.46 s), and 0.35 rad/s (corresponding to a wave period of 17.95 s, with a value of 14,085 N/m).In addition, the peak value of the mean drift force in the sway direction is achieved at a wave frequency of approximately 1.75 rad/s (corresponding to a wave period of 3.59 s) and 0.35 rad/s (corresponding to a wave period of 17.95 s, with a value of 18,436 N/m).Furthermore, the peak value of the mean drift force in the yaw direction occurs at a wave frequency of around 1.55 rad/s (corresponding to a wave period of 4.05 s) and 1.85 rad/s (corresponding to a wave period of 3.39 s).surge, sway, and yaw under different wave incident angles obtained using these two methods.The results indicate that the mean drift forces calculated using the middle-field method agree with those obtained using the far-field method, demonstrating that the mid dle-field method is a reliable approach to calculating the full quadratic transfer functions in the analysis.The maximum value of the mean drift force in the sway direction is signif icantly larger than that in the surge direction.This can be attributed to the larger area subjected to forces in the sway direction compared with the surge area.The peak values of the mean drift force in the sway direction occur at approximately 1.6 rad/s (correspond ing to a wave period of 3.93 s), 1.15 rad/s (corresponding to a wave period of 5.46 s), and 0.35 rad/s (corresponding to a wave period of 17.95 s, with a value of 14,085 N/m).In addition, the peak value of the mean drift force in the sway direction is achieved at a wave frequency of approximately 1.75 rad/s (corresponding to a wave period of 3.59 s) and 0.35 rad/s (corresponding to a wave period of 17.95 s, with a value of 18,436 N/m).Furthermore, the peak value of the mean drift force in the yaw direction occurs at a wave fre quency of around 1.55 rad/s (corresponding to a wave period of 4.05 s) and 1.85 rad/s (corresponding to a wave period of 3.39 s).For the platform without nets in a free-floating state, natural periods of he roll,and pitch motions are calculated using the software HydroStar.The calculated nat periods are as follows: 16.908 s for heave, 16.801 s for roll, and 18.208 s for pitch.Nota these periods are far from the wave periods typically existing in the East and South C seas.This indicates that the designed platform is capable of effectively mitigating st tural resonance caused by waves.

Hydrodynamic Performance in Regular Waves
This study employs the software OrcaFlex 11.0a to conduct a coupled analysis o platform's dynamic responses.OrcaFlex is a software for dynamic analysis of marine tems, which can simulate various types of structures.A fully nonlinear numerical m is constructed in OrcaFlex.
Figures 17-19 illustrate the surge, heave, and pitch RAOs of the moored platf with and without nets in regular waves with a heading angle of 180°.The wave hei considered are 1 m and 3.5 m, with the periods ranging from 5 s to 30 s.It is evident the surge, heave, and pitch RAOs of the platform are quite small within the wave pe of 5-15 s, where the majority of energy is concentrated, regardless of the presence of n For the platform without nets in a free-floating state, natural periods of heave, roll, and pitch motions are calculated using the software HydroStar.The calculated natural periods are as follows: 16.908 s for heave, 16.801 s for roll, and 18.208 s for pitch.Notably, these periods are far from the wave periods typically existing in the East and South China seas.This indicates that the designed platform is capable of effectively mitigating structural resonance caused by waves.

Hydrodynamic Performance in Regular Waves
This study employs the software OrcaFlex 11.0a to conduct a coupled analysis of the platform's dynamic responses.OrcaFlex is a software for dynamic analysis of marine systems, which can simulate various types of structures.A fully nonlinear numerical model is constructed in OrcaFlex.
Figures 17-19 illustrate the surge, heave, and pitch RAOs of the moored platform with and without nets in regular waves with a heading angle of 180 • .The wave heights considered are 1 m and 3.5 m, with the periods ranging from 5 s to 30 s.It is evident that the surge, heave, and pitch RAOs of the platform are quite small within the wave period of 5-15 s, where the majority of energy is concentrated, regardless of the presence of nets.The corresponding surge RAOs are less than 1 m, the heave RAOs are significantly smaller than 1 m, and the pitch RAOs are less than 1 • .These results demonstrate the excellent hydrodynamic performance of the new platform.For wave periods larger than 15 s, under the conditions of H = 1 m, the impact of nets on the surge, heave, and pitch RAOs is minimal for most periods.However, the presence of nets reduces the RAOs near the natural period of heave and pitch, which may be attributed to the increased damping effect induced by the nets.Similar effects of nets on the peak RAOs of heave and pitch can be found in Figure 9 of [48].For wave periods larger than 15 s, under the conditions of H = 3.5 m, the effects of nets on the RAOs are complex.The existence of nets reduces the values of surge and heave RAOs for most periods, but the opposite results are observed for the roll RAOs.Additionally, incident waves with larger wave height correspond to smaller surge, heave, and pitch RAOs of the platform for wave periods larger than 15 s.

Dynamic Response in Waves, Winds, and Currents
For each loading case, five random waves with different seed numbers are generated, and the simulation duration is set to three hours.The calculations for determining the response of the platform's motion and the line tension were conducted using a time step of 0.01 s.In the analysis, the results of the first 1000 s are truncated for reliability purposes.The tension of mooring lines and the motions of six-degree-of-freedom of the platform are calculated in accordance with the rules [42].
The safety factor is defined as the ratio of the minimum breaking load on the mooring line to the maximum tension of the mooring line.For the dynamic analysis method, the safety factor of the mooring line is 1.67 for intact survival conditions and 1.25 for damaged survival conditions.The permissible mean and maximum offsets are generally specified based on the water depth of operating areas, environmental conditions, clearance requirements, and limitations of the system [42].

Mooring Line Tension
Table 9 presents the statistics of the mooring line tension (the average value of the results obtained from five different seed numbers) under various loading cases for survival conditions.The statistics of mooring line tension for operation conditions are omitted to reduce the length of this paper.For intact survival conditions, the maximum tension occurs in Line 1, which is 2460 kN when the wave, wind, and current approach simultaneously from 90 • (LC5).For damaged survival conditions, the maximum tension occurs in Line 2, which is 3337 kN.It also occurs when the wave, wind, and current come from 90 • simultaneously (LC10).The minimum safety factor of all mooring lines for all intact survival conditions is 1.81 in LC5, which exceeds the CCS requirement of 1.67.Similarly, the minimum safety factor of all mooring lines for all damaged survival conditions is 1.33 in LC10, exceeding the CCS requirement of 1.25.For all the loading cases, the safety factors of the mooring lines meet the design rule requirements.Consequently, the mooring system designed for this semi-submersible aquaculture platform can be deemed reliable.For both the intact survival conditions and the damaged survival conditions, the most critical loading case occurs when the wave, wind, and current have an incident angle of 90 • .The reason is that the platform has a larger lateral area.When the incident wave, wind, and current are transverse, the platform is subject to greater forces.When one of the two mooring lines in the same group is damaged, the tension of the remaining mooring line increases significantly.For example, when comparing the conditions LC5 and LC10, the maximum line tension of Line 2 increases from 1262 kN to 3337 kN, representing a significant increase of 164%.Additionally, the average line tension of Line 2 increases from 285 kN to 937 kN, which corresponds to a substantial increase of 228%.Three typical loading cases (LC1, LC3, and LC5), were considered to analyze the net effects on the mooring line tension.Table 10 presents the statistics of the tension of the weakest mooring line under these conditions.It is evident that for all these loading cases, there are significant differences among the maximum, minimum and mean mooring line tensions, which is also verified by the standard deviation.For instance, the minimum and mean tensions of Line 1 for LC5 without nets are 2 kN and 420 kN, respectively.However, the maximum tension of Line 1 is 2471 kN, representing an increase of 488% compared with the mean tension.Moreover, it can be observed that the average and minimum mooring line tensions with nets are larger than those without nets, while the maximum mooring line tension with nets is smaller than that without nets.This suggests that the presence of nets enhances the safety of the entire mooring system, as it is determined by the maximum value of the mooring line.Li et al. [24] also investigated the net effects on the mooring line tension and found that the presence of nets can alleviate the maximum mooring line tension in pure wave conditions.However, when the current coming from the same direction as the wave is included, the different solidity ratio of nets can either increase or decrease the maximum mooring line tension compared with the case without nets.This implies that, for the platform in this study, changes in the solidity ratio of nets and current may also increase the maximum mooring line tension, contrary to the results shown in Table 10.Figures 20-22 illustrate the segments of time histories of the mooring line tension where the maximum line tension occurs.Under the conditions LC1 and LC3, the mooring line tension with nets is generally larger than that without nets.However, the maximum mooring line tension exists in the case with nets, which is consistent with the results in Table 10.Under the conditions LC5, both the maximum mooring line tension and the variation range of mooring line tension with nets are larger than those without nets.The minimum amplitude of surge motion is observed at the condition LC5, as the loads at a heading angle of 90 • have little impact on surge motion.However, the maximum amplitude of surge motion is observed at the condition LC4, rather than LC1.This is likely due to the asymmetry in the longitudinal direction.Similar results for the amplitude of surge motion can be found for damaged survival conditions.The amplitude of sway motion increases gradually with the incident angle gradually varying from 180 • (LC1) to 90 • (LC5), for both intact and damaged survival conditions.The sway motion is minimally impacted by the condition of the mooring lines.Under the conditions LC1, LC2, and LC3, the amplitude of heave motion is almost the same, measuring 0.88 m for LC1, 0.87 m for LC2, and 0.85 m for LC3.Under the condition LC4, the amplitude of heave motion increases to 1.34 m, a 58% increase compared with the condition LC3.The maximum amplitude of heave motion, reaching 2.22 m, is observed at the condition LC5.This indicates that the loads of 90 • heading angle have the largest effects on the motion of heave, which is consistent with the variation of the mooring line tension of Line 1.The differences in heave motion between intact and damaged survival conditions are not significant.It is found that the amplitude of roll motion increases as the incident angle gradually changes from 180 • (LC1) to 90 • (LC5), with the largest amplitude of 8.96 • occurring at the condition LC5.When one of the mooring lines is damaged, the change in roll motion is not noticeable compared with the intact survival condition.Similar to the variation observed in surge motion amplitude, the maximum amplitude of pitch motion is observed at the condition LC4, rather than the condition LC1.The maximum value of pitch motion is 3.41 • .Under the conditions LC1 and LC5, the motion of the yaw is quite small.However, the amplitude of yaw motion increases obviously under the condition LC3 for both the intact survival and damaged survival conditions.
The typical loading case of LC3 is further considered to analyze the effects of nets on the six-degree-of-freedom motion of the platform.Table 12 presents the statistics of the motion of six-degree-of-freedom under the condition LC3.It is found that the range of roll motion without nets increases significantly compared with the case with nets.The amplitude varies from 4.55 • to 6.59 • , with an increase of 45%.Compared with the case with nets, the range of pitch motion without nets increases slightly, and the amplitude changes from 5.45 • to 5.83 • , with an increase of 7%.Compared with the case with nets, the range of yaw motion without nets also increases slightly by 7%.The range of heave motion without nets significantly increased from 1.66 m to 2.27 m compared with the case with nets, with an increase of 36%.Lei et al. [48] investigated the effects of nets on the motion of a floating offshore wind turbine with a steel fish farming cage (FOWT-SFFC) and found that the presence of nets can also alleviate the motion of FOWT-SFFC under certain conditions.Figure 23 shows the partial results of the whole-time histories of the heave of the platform under the condition LC3. Figure 24 shows the planar motion of the platform in the whole-time histories under the condition LC3.From Figure 23, it is observed that the variation range of heave motion with nets is apparently smaller than that without nets.From Figure 24, it is observed that the range of the whole plane motion becomes larger significantly without nets, and both the amplitude ranges of the sway and the surge obviously increase, compared with those with nets.These findings imply that the presence of nets alleviates the motion response of the platform.

Conclusions
In this paper, a numerical study of the hydrodynamic performance of a newly designed semi-submersible aquaculture platform, developed by GIEC and put into operation, is conducted.The present study proposes a numerical model that considers the semisubmersible main structure, the net system, and the mooring system to accurately capture the dynamic response of the platform under wave, wind, and current conditions.The wave force is calculated based on the linear potential theory and the boundary element method.The wind force and current force are calculated according to the rules of CCS.The net system is modeled based on the Morison equation and the equivalent simulation.The mooring system is modeled based on the lumped-mass method and the finite element method.Hydrodynamic coefficients and wave forces are calculated by the software Hy-droStar.The fully coupled motion equation is solved by the software OrcaFlex.
The main findings of this study are summarized as follows: Firstly, the reliability of the middle-field method to calculate the second wave force is confirmed through compar-

Conclusions
In this paper, a numerical study of the hydrodynamic performance of a newly designed semi-submersible aquaculture platform, developed by GIEC and put into operation, is conducted.The present study proposes a numerical model that considers the semisubmersible main structure, the net system, and the mooring system to accurately capture the dynamic response of the platform under wave, wind, and current conditions.The wave force is calculated based on the linear potential theory and the boundary element method.The wind force and current force are calculated according to the rules of CCS.The net system is modeled based on the Morison equation and the equivalent simulation.The mooring system is modeled based on the lumped-mass method and the finite element method.Hydrodynamic coefficients and wave forces are calculated by the software Hy-droStar.The fully coupled motion equation is solved by the software OrcaFlex.
The main findings of this study are summarized as follows: Firstly, the reliability of the middle-field method to calculate the second wave force is confirmed through compar-

Conclusions
In this paper, a numerical study of the hydrodynamic performance of a newly designed semi-submersible aquaculture platform, developed by GIEC and put into operation, is conducted.The present study proposes a numerical model that considers the semi-submersible main structure, the net system, and the mooring system to accurately capture the dynamic response of the platform under wave, wind, and current conditions.The wave force is calculated based on the linear potential theory and the boundary element method.The wind force and current force are calculated according to the rules of CCS.The net system is modeled based on the Morison equation and the equivalent simulation.The mooring system is modeled based on the lumped-mass method and the finite element method.
Hydrodynamic coefficients and wave forces are calculated by the software HydroStar.The fully coupled motion equation is solved by the software OrcaFlex.
The main findings of this study are summarized as follows: Firstly, the reliability of the middle-field method to calculate the second wave force is confirmed through comparing the results of the far-field method and the middle-field method for the mean drift wave force.Secondly, the analysis of the hydrodynamic performance of the platform in regular waves reveals that the natural period of the platform is far away from the common wave period in the specific sea area where the project is located.Thirdly, when subjected to a combination of waves, winds, and currents, the maximum tension of the mooring line is found to be 2460 kN with a safety factor of 1.81 for intact survival conditions and 3337 kN with a safety factor of 1.33 for damaged survival conditions.These parameters meet the requirements guided by CCS, indicating the reliability of the designed mooring system.Furthermore, the motion of six-degree-of-freedom is deemed acceptable for the platform's operation in the specific sea area.Lastly, the net system is found to have noticeable effects on the hydrodynamic performance of the platform.Unter certain conditions, the maximum mooring line tension with nets is smaller than that without nets.The range of roll motion without nets is increased by 45% compared with the case with nets, and the range of heave motion without nets is increased by 36% compared with the case with nets.Therefore, it is recommended that the net system should be carefully considered during the design stage to ensure the platform's ability to maintain safety in severe environments.
This study provides valuable insights for researchers and engineers to design and optimize this kind of semi-submersible aquaculture platform in the future.However, it is important to acknowledge several limitations of this study and identify future research directions.Firstly, the linearized potential flow models utilized in this study do not incorporate the nonlinearity caused by breaking waves.To accurately capture these nonlinear effects, advanced numerical methods such as computational fluid dynamics (CFD) simulations should be employed in future studies.Secondly, this study does not consider the installation of Sharp Eagle wave energy converters on the platform.Investigating the effects of these converters on the hydrodynamic performance of the platform should be a priority for future research.Lastly, the mooring system modeled in this study is relatively simple.Future investigations should focus on exploring more complex mooring systems to better understand their impact on the platform's performance.By acknowledging these limitations and suggesting future research directions, this study lays the foundation for continued exploration and advancement for this kind of platform.

Data Availability Statement:
The data used to support the findings of this study are available from the corresponding author upon request.The data are not publicly available due to privacy.
and indirect time domain to analyze the motion responses of the cage.Physical model tests in the wave basin verified the numerical model.South China Sea Fisheries Research Institute and Zhuhai Desai Marine Fishery Technology Co., Ltd.(Zhuhai, China) co-designed a semi-submersible fish farm named Dehai Farm 1.It consists of the truss-work system and ballast tanks with dimensions of 91.3 m in length, 27.6 m in width, and 10.3 m in height.Huang et al.

Figure 2
Figure 2 displays the newly designed semi-submersible aquaculture named Pu Sheng 1, created by GIEC and constructed by Guangzhou Shipya tional Co., Ltd.(Guangzhou, China) for Hainan Pusheng Marine Science and Development Co., Ltd.(Ledong, China) This platform has earned a classifica cate from the China Classification Society (CCS), making it the first marin equipment in South China to do so.

Figure 2
Figure 2 displays the newly designed semi-submersible aquaculture platform, named Pu Sheng 1, created by GIEC and constructed by Guangzhou Shipyard International Co., Ltd.(Guangzhou, China) for Hainan Pusheng Marine Science and Technology Development Co., Ltd.(Ledong, China) This platform has earned a classification certificate from the China Classification Society (CCS), making it the first marine breeding equipment in South China to do so.

Figure 2
Figure 2 displays the newly designed semi-submersible aquaculture p named Pu Sheng 1, created by GIEC and constructed by Guangzhou Shipyard tional Co., Ltd.(Guangzhou, China) for Hainan Pusheng Marine Science and Tec Development Co., Ltd.(Ledong, China) This platform has earned a classification cate from the China Classification Society (CCS), making it the first marine b equipment in South China to do so.

Figure 4 .
Figure 4. Pu Sheng 1 in the operation condition.

Figure 5 .
Figure 5.The location of the project site.

Figure 4 .
Figure 4. Pu Sheng 1 in the operation condition.As shown in Figure 5, the longitude and latitude coordinates of the platform's location are 18 • 19 29.839 N and 108 • 50 13.633E, respectively.

Figure 4 .
Figure 4. Pu Sheng 1 in the operation condition.

Figure 5 .
Figure 5.The location of the project site.

Figure 5 .
Figure 5.The location of the project site.

Figure 6 .
Figure 6.Division of flow field region.
the principle of force equivalence, the drag coefficient of the single net line is ' =1.525 d C , and the inertia coefficient of the single net line is ' =1.0 C .

Figure 6 .
Figure 6.Division of flow field region.

Figure 7 .
Figure 7.The schematic diagram of net equivalent substitution.

Figure 7 .
Figure 7.The schematic diagram of net equivalent substitution.

Figure 9 .
Figure 9. Definition of the directions of waves, winds, and currents.

Figure 9 .
Figure 9. Definition of the directions of waves, winds, and currents.

Figure 12 .
Figure 12.RAO results of pitch for different mesh sizes at the 180° heading angle.

Figure 12 .Figure 12 .
Figure 12.RAO results of pitch for different mesh sizes at the 180° heading angle.Figure 12. RAO results of pitch for different mesh sizes at the 180 • heading angle.

Figure 16 .
Figure 16.The mean drift force of yaw at the 90 • heading angle.

Figure 19 .
Figure 19.Pitch RAOs for the moored platform with/without nets.

Figure 23 .
Figure 23.Time history result of the motion of the heave of the platform under the condition LC3.

Figure 24 .
Figure 24.Time history result of the motion of the plane under the condition LC3.

Figure 23 .Figure 23 .Figure 24 .
Figure 23.Time history result of the motion of the heave of the platform under the condition LC3.

Figure 24 .
Figure 24.Time history result of the motion of the plane under the condition LC3.

Table 1 .
Detailed parameters of the platform.

Table 4 .
Detailed parameters of the nets.

Table 4 .
Detailed parameters of the nets.

Table 5 .
Equivalent coefficients of nets.

Table 6 .
Detailed parameters of mooring lines.

Table 7 .
Coordinates of fairleads and anchors.
, x FOR PEER REVIEW 3.3.2.Loading Cases for Mooring Analysis

Table 8 .
Loading cases of survival conditions.

Table 9 .
The mooring line tension under survival conditions.
* SF represents the safe factor.

Table 10 .
Statistics of mooring line tension under typical conditions.Std.Dev.represents the Standard Deviation of the mean value of the parameter. *

Table 12 .
Statistics of the six-degree-of-freedom motion under the condition LC3.