Initial values of the numerical model.^{ a} Main properties of the medium grid are shown in Table 2.

## 1. Introduction

Aquaculture is expanding all over the world, and the net cage is becoming prevalent in the aquaculture industry. Good knowledge of the hydrodynamic characteristics of the net cage has important significance to the design of the net cage and the welfare of the fish. An important component of a net cage is the fishing net, where the fish are kept and grown. The fishing net is completely submerged and infinitely flexible. When exposed to current, the fishing net changes its shape to reduce the hydrodynamic force acting on it, and the deformed net in turn affects the flow field around the fishing net. In this way, the flow field and the deformation of the fishing net interact and mutually influence each other. Therefore, it is very necessary to consider the fluid-structure interaction between flow and the fishing net when calculating the hydrodynamic characteristics of the net cage.

To investigate the hydrodynamic characteristics of net cages, much progress has been made by both experimental measurements and numerical simulations. Aarsnes et al. [1] conducted a series of physical model tests to measure the flow-velocity reduction, fishing-net deformation and the drag force acting on net cages. Lader and Enerhaug [2] performed experiments to investigate the forces and geometry of a net cage in uniform flow, and the flow-velocity reduction in the centre of the net cage was given as 20%. It was found that the global forces acting on a flexible fishing net were not proportional to the square of the flow velocity for all currents. Fredriksson et al. [3] conducted a numerical study on mooring tensions of a large fish farm containing 20 net cages while the flow-velocity reduction was estimated using a simplified control volume approach. DeCew et al. [4] investigated the submergence behavior of a fish cage in a single-point mooring system under currents by a numerical model and a reduction of 12.7% of the current velocity was applied to the aft portion of the net chamber. Huang et al. [5] and Xu et al. [6] analyzed the hydrodynamic behaviors of net cages in waves and a reduction factor of 0.85 was adopted in the numerical model to represent the retardation effect on the water particle speed after passing through each previous net.

Based on the above analysis, the reduction in current downstream from the fishing net has been taken into consideration in studying the hydrodynamic force and the effective volume of the net cage. However, only a reduction factor was estimated to represent the shielding effect of the fishing net for most of the numerical works. There has not been a numerical model focusing on the fluid-structure interaction between flow and a flexible net. As an extension of our previous works [7-8], a coupled fluid-structure model is developed based on a joint use of the porous-media fluid model and the lumped-mass mechanical model. Using an appropriate iterative scheme, the fluid-structure interaction between flow and fishing net can be solved.

## 2. Model descriptions

The coupled model includes two numerical models: the porous-media fluid model and the lumped-mass mechanical model. The porous-media fluid model can simulate the flow field around a plane net with no deformation, and the lumped-mass mechanical model can simulate the deformation of fishing net in certain current. The joint use of the above two models is presented to solve the fluid-net interaction problem in current. The steady flow field around fishing net can be obtained by using an appropriate iterative scheme as explained later in Section 2.3. This model has been introduced detailed in our previous studies [9-10].

### 2.1. The porous-media fluid model

The flow field around a plane net is simulated using a commercial computational fluid dynamics (CFD) software FLUENT from ANSYS Inc. The porous-media fluid model is introduced to model the plane net, and the finite volume method is used to solve the governing equations of the numerical model. Thus, a brief review of the porous-media fluid model is given here.

The governing equations of the porous-media fluid model are mainly the Navier-Stokes equations as follows:

where *t* is time; *µ* and *ρ* are the viscosity and the density of the fluid, respectively; *µ*_{t} is the eddy viscosity; *P*=*p*+(2/3)*ρk*, where *p* is pressure; *k* is the turbulent kinetic energy; *u*_{i} and *u*_{j} are the average velocity components, respectively; *g*_{i} is the acceleration due to gravity; *i*, *j*=1, 2, 3 (*x*, *y*, *z*); and *S*_{i} is the source term for the momentum equation.

The porous-media fluid model is a hypothetical model which produces the same water-blocking effect as the fish net by setting the coefficients of the hypothetical porous media. For flow through porous media, a pressure gradient exists:

where a and b are constant coefficients and

Outside the porous media, the source term for the momentum equation is given by *S*_{i}=0. Inside the porous media, *S*_{i} is calculated by the following equation:

where *C*_{ij} is given by

where *C*_{ij} represents material matrices consisting of porous media resistance coefficients, *C*_{n} is the normal resistance coefficient, and *C*_{t} is the tangential resistance coefficient.

For flow through the porous media, the drag force (*F*_{d}) and the lift force (*F*_{l}) acting on it can be expressed as follows:

where *λ* is the thickness of the porous media and *A* is the area of the porous media.

The porous coefficients in Eqs. (6) and (7) can be calculated from the drag and lift forces, while the drag and lift forces acting on a fishing net have strong relationship with the features of the net. So the connection between net features and flow field is built by the porous coefficients in the numerical model. Generally, the drag and lift forces of the plane net are obtained from laboratory experiments. These forces can also be calculated from the Morison equation:

where *C*_{d} and *C*_{l} are coefficients that can be calculated using empirical formulas [1, 12-14].

The porous coefficient *C*_{n} can be calculated from a curve fit between drag forces acting on a plane net and corresponding current velocities using the least squares method when the plane net is oriented normal to the flow. The other coefficient *C*_{t} can be obtained when a plane net is oriented at an attack angle (see Figure 1a). In this case, the coefficient *C*_{t} should be transformed into *C*_{tα} using Eq. (10) [15]. Next, *C*_{tα} can be calculated from a curve fit between lift forces acting on a plane net and corresponding current velocities using the least squares method.

where *α* is the attack angle that is defined to be the angle between the flow direction and the plane net in the horizontal plane; *C*_{tα} is the tangential inertial resistance coefficient for the plane net at an attack angle *α*.

When a plane net is oriented at an arbitrary position (see Figure 1b), the porous coefficients should be transformed using Eq. (11):

where *β* is the angle between the normal vector of the plane net (*z*-direction. It should be noted that the cosines of *α* and *β* can be calculated between the corresponding two vectors.

### 2.2. The lumped-mass mechanical model

The lumped-mass mechanical model is introduced to simplify the fishing net, and the motion equations can be established mainly based on Newton's Second Law. Given an initial net configuration, the motion equations can be solved numerically for each lumped-mass point. Finally, the configuration of the net in current can be simulated. The calculation method for the model has been explained fully in our previous studies [16-17]. Here only a brief outline is described.

In a uniform current, the motion equations of lumped-mass point *i* can be expressed as follows:

where *M*_{i} and Δ*M*_{i} are the mass and added mass of lumped-mass point *i*, *ij* (*j* is the code for the knot at another end of the bar *ij*), *n* is the number of adjacent knots of point *i*,

To consider the direction of the hydrodynamic forces acting on mesh bars, local coordinates (*τ*, *η*, *ξ*) are defined to simplify the procedure (see Figure 2). The origin of the local coordinate is set at the center of a mesh bar, and the *η*-axis lies on the plane including *τ*-axis and flow velocity *τ*-direction can be obtained from the Morison equation as follows:

where *C*_{Dτ} is the drag coefficient in the *τ*-direction, *A*_{τ} is the projected area of the twine normal to the *τ*-direction, *τ*-direction. Similar expressions can be applied to the other two components of the drag force: *F*_{Dη} and *F*_{Dξ}.

The mesh-grouping method [18-19] is used to reduce the computational effort of the net. Each grouped mesh is defined as a plane-net element (see Figure 3). The normal vector of the plane-net element (

where

Because of the deformation of the net, the area of a plane-net element is not constant, thus impacting the solidity of the plane-net element. Consequently, the solidity *S*_{n} must be corrected accordingly based on the varying area of the plane-net element:

where *S*_{n} and *A* are the solidity and area of the undeformed plane-net element, and

### 2.3. The coupled fluid-structure interaction model

Based on the above work, the interaction between flow and fishing net can be simulated by dividing the net into many plane-net elements (see Figure 4). The main concept of the numerical approach is to combine the porous-media fluid model and the lumped-mass mechanical model to simulate the interaction between flow and flexible nets. Taking a single flexible net as example, a calculating flow chart for the numerical approach is shown in Figure 5. A more detailed calculation procedure is given as follows.

The calculation procedure includes four steps. Step 1 is to simulate the flow field around a vertical plane net without considering the net deformation using the porous-media fluid model. Step 2 is to calculate the drag force and the configuration of a flexible net in current using the lumped-mass mechanical model. The given initial flow velocity is the average value exported from the upstream surface of the plane net in last step. In Step 3, according to the configuration of the net, the orientation and the porous coefficients of each plane-net element can be obtained, thus obtaining the flow field around the fishing net. In Step 4, the drag force and the configuration of the net are calculated again. The flow velocity acting on each lumped-mass point is the average value exported from the upstream surface of the corresponding plane-net element in last step.

The terminal criterion of the calculation procedure is defined as *F*_{i} and *F*_{i+1} are the drag forces of two adjacent calculations. If Δ*F*< 0.001, it was determined that the force acting on the fishing net is stable and the configuration of the net is the equilibrium position considering the fluid-structure interaction. Thus the calculation procedure is completed. Otherwise, continue with the next iteration of the procedure and repeat Steps 3 and 4 until the criterion meets the precision requirement. Finally, the fluid-structure interaction problem can be solved and the steady flow field around the fishing net can be obtained.

## 3. Discussions of model parameters

Due to uncertainties in numerical simulation, convergence studies were performed to analyze the influence of different model parameters, such as thickness of the porous media, grid refinement, time-step size, number of net elements, number of iterations, porous coefficients, etc. When varying one parameter, all other parameters were kept at their initial values as shown in Table 1.

Parameter | Value |

Thickness of porous media | 20 mm |

Grid refinement | Grid 2 ^{a} |

Time-step size | 0.04 s |

Number of net elements | 5×5 |

Number of iterations | 2 |

Porous coefficients | Cn =35.3 m^{-1} |

Ct =10.0 m^{-1} |

Grid | Min. cell size (mm) | Max. cell size (mm) | Total cells |

Grid 1 | 5 | 50 | 413,697 |

Grid 2 | 10 | 50 | 117,294 |

Grid 3 | 20 | 50 | 49,807 |

Grid 4 | 10 | 60 | 34,492 |

Grid 5 | 10 | 70 | 21,124 |

The modeled net was a 0.3 m×0.3 m knotless polyethylene (PE) net with 15 meshes in width and 15 meshes in height. The twine diameter was 2.6 mm, and the mesh bar length was 20 mm. Mounted as square meshes, the net solidity ratio was 0.243.

### 3.1. Thickness of the porous media

A study on the thickness of the porous media was performed for a plane net at an attack angle *α*=90° and incoming velocity *u*_{0}=0.226 m/s. The corresponding Reynolds number is 5.88×10^{2} calculated according to the twine diameter of the net. The dimensions of the porous media are the same as that of the fishing net in both width and height, and the porous media has a certain thickness. In order to investigate the influence of the thickness on the numerical results, the porous-media fluid model with different thicknesses, *λ*=5, 10, 20, 30, 40 and 50 mm, were performed.

The thickness of the porous media can affect the distribution of the flow velocity around the fishing net to some extent (see Figure 6). However, the numerical results of the flow velocity are almost the same for the thickness of 5, 10 and 20 mm. It was determined that the thickness of the porous media does not affect the numerical results obviously until the thickness is greater than 20 mm. Therefore, a thickness of 20 mm is adopted hereafter for all the computations considering the total number of cells and computational accuracy. An increase in the thickness of the porous media can effectively reduce the number of cells, thus reducing the computational effort greatly. However, it is suggested that the maximum thickness of the porous media should not be greater than 10% of width or height of the plane net.

### 3.2. Grid refinement

The influence of grid on the numerical results was investigated by repeating the computations with different levels of refinement. Five grids as shown in Table 2 were taken into account. The minimum cell size was adopted inside the porous-media region. A growth rate of 1.2 was set around the net model, and the maximum cell size in the computational domain was set to control the distribution of the grid.

Numerical results of flow velocity obtained with different grids were compared (see Figure 7). It is found that numerical results are free of grid influence for Grid 1 and Grid 2. But for the coarser grid models, the numerical results tend to divergence, especially downstream from the fishing net. Therefore, Grid 2 is adopted hereafter for all the computations.

### 3.3. Time-step size

In order to investigate the influence of the time-step size, the computations were repeated with different time-step sizes: 0.01, 0.02, 0.04, 0.06, 0.08 and 0.10 s. The flow velocity both upstream and downstream from the plane net shows a convergence result for the time-step size smaller than 0.04 s (see Figure 8). Therefore, a time-step size of 0.04 s is adopted for computational effort saving. The computational effort decreases with increasing time-step size. However, it is suggested that the time-step size should not exceed the ratio between the minimum value of the cell size and the incoming velocity. Because the larger the time-step size is, the greater the residual error of each step will be, thus causing a divergence result without meeting the convergence criteria.

### 3.4. Number of plane-net elements

In the numerical model, each mesh is defined as a plane-net element. In order to reduce the computational effort, the mesh-grouping method was used to reduce the number of plane-net elements. The original model of the fishing net was a 15×15 model with 15 meshes in width and 15 meshes in height. Two grouped models were used to simplify the original model. For the original model, nine meshes were grouped into one mesh, thus obtaining a 5×5 model with a total of 25 elements. A 3×3 model with 9 elements can be obtained by grouping 25 meshes into one mesh.

The calculated net configuration of the 5×5 model is almost coincided with that of the original 15×15 model (see Figure 9). The result of the 3×3 model agrees well with that of the original model except for a slight discrepancy at the bottom of the net. There is no significant difference between the results of flow velocity obtained with the grouped models and the original model (see Figure 10). Therefore, the mesh-grouping method can be used to reduce the number of plane-net elements without affecting the calculation accuracy. Because of better agreement on net configuration, the 5×5 model is chosen for further numerical simulations.

### 3.5. Number of iterations

The interaction between flow and fishing net can be solved after several iterations. Taking a single net as example, the calculation precision reached the terminal criterion, Δ*F*< 0.001, within two iterations, and the calculation procedure terminated for incoming velocity of 0.226 m/s. The drag-force difference between the 1st iteration and the 2nd iteration was approximately 0.7%. One more iteration was also calculated, and the result tended to converge.

Comparing the net configurations after different iterations, there is a little discrepancy between the 1st iteration and the 2nd iteration. However, the calculated net configuration of the 3rd iteration is almost coincided with that of the 2nd iteration (see Figure 11). The flow velocity upstream of the flexible net is greater than that of the plane net with no deformation, while the flow velocity downstream from the flexible net is smaller (see Figure 12). It was determined that the deformed net has more obvious shielding effect than the plane net with no deformation. There exists rather obvious discrepancy of the flow-velocity distribution between the 1st iteration and the 2nd iteration (see Figure 12), and the maximum discrepancy is approximately 6.3%. The result of the 2nd iteration is almost the same as that of the 3rd iteration, and the maximum discrepancy is approximately 0.17%. Therefore, two iterations is adopted to model a single net for computational effort saving.

### 3.6. Porous coefficients

A set of porous coefficients, *C*_{n}=35.3 m^{-1} and *C*_{t}=10.0 m^{-1}, can be obtained from the experimental data (see Tables 3 and 4). The corresponding laboratory experiment is presented in our previous study [9]. Porous coefficients are key to the numerical results. Before the numerical simulations, a sensitivity study on the flow field with respect to the porous coefficients was carried out with a plane net at different attack angles.

u (m/s) | 0.058 | 0.113 | 0.170 | 0.226 |

Fd (N) | 0.082 | 0.395 | 0.968 | 1.678 |

u (m/s) | 0.058 | 0.113 | 0.170 | 0.226 |

Fl (N) | 0.032 | 0.131 | 0.321 | 0.534 |

When *C*_{n} is kept at a constant, the coefficient *C*_{t} increases with an increment of 2 from 6 to 12. As shown in Figure 13a, varying *C*_{t} has no effect on the flow field around the plane net for attack angle of 90°. However, when the plane net is positioned at a smaller attack angle; the flow velocity decreases with increasing *C*_{t}. As the attack angle decreases, the influence of variation of *C*_{t} to the flow velocity increases. When *C*_{t} is a constant, the flow velocity decreases linearly with increasing *C*_{n} (see Figure 13b). In addition, the variation tendencies are almost the same for different attack angles. With a same increment of 6, the coefficients *C*_{t} and *C*_{n} have different influence on the flow field around the plane net. The maximum discrepancies in flow velocity are 4.3% for the variation of *C*_{t} at an attack angle of 30° and 2.1% for the variation of *C*_{n} at an attack angle of 90° respectively.

### 3.7. Brief conclusion

In summary, numerical results depends on the different model parameters to some extent. Appropriate values of the thickness of the porous media, the grid refinement and the time-step size are beneficial to reducing the computational effort and keeping high calculation accuracy. The mesh-grouping method can be used to reduce the number of plane-net elements and the computational effort. Two iterations are sufficient to achieve a converged solution of a flexible fishing net. The porous coefficients have a significant effect on the flow velocity around the fishing net and the relationship between each coefficient and the flow velocity is approximately linear. In addition, the influence of the coefficient *C*_{t} to the flow velocity is different with different attack angles. Based on the above, the initial values (see Table 1) are the optimal parameters for the numerical model.

## 4. The advantage of the coupled model

In present numerical model, the fluid-structure interaction (FSI) between flow and fishing net is taken into consideration. Because of the shielding effect of the fishing net, there exist the flow-velocity-reduction regions upstream of and downstream from the fishing net (see Figure 14). According to the numerical simulation, the reduction regions become smaller when considering FSI due to the decrease in height of the fishing net. However, the flow-velocity reduction is larger than that without considering FSI downstream from the fishing net. In the numerical results, the streamlines present the flow direction (see Figure 14). It was determined that the diversion of flow direction around the fishing net is relatively small but more obvious than that without considering FSI.

Comparing the flow-velocity distribution along the *x*-direction, there exists rather obvious discrepancy between the fishing net with and without considering FSI (see Figure 15). The flow velocity upstream of the fishing net is greater than that without considering FSI, while the flow velocity downstream is smaller. The maximum flow-velocity reductions downstream are 15.0% for the FSI model and 11.9% for the model without considering FSI, respectively. It was found that the deformed net has more obvious shielding effect than the net with no deformation. This is due to the deformation of the flexible net. When exposed to current, the flexible net changes its shape and deviates from the original position. If the net is evenly divided into several plane-net elements, the angle between the incoming current and each net element (attack angle) decreases. Thus there is a relative increase in net solidity along the flow direction, and resulting in a more obvious blockage effect.

As shown in Figure 16, the deformation of the fishing net is less serious than that without considering FSI. This is because that the net configuration is calculated based on the actual velocity acting on the flexible net instead of the incoming velocity in present numerical model. The actual velocity is slightly smaller than the incoming velocity due to the blockage effect of the net. As a result, the drag force is a bit smaller at approximately 7-8% than that without considering FSI for different incoming velocities (see Figure 17).

According to the analyses above, there are obvious difference between the flow field, the hydrodynamic force and the deformation of the fishing net with and without considering FSI. For multiple nets, the influence of FSI on the hydrodynamic characteristics of the fishing net will become more obvious. Therefore, it is necessary to take FSI into consideration in studying the interaction between flow and fishing net.

## 5. Conclusions

A numerical model is proposed to simulate the interaction between flow and fishing net based on the joint use of the porous-media fluid model and the lumped-mass mechanical model. Convergence studies are performed to analyze the influence of different model parameters on the numerical results and a set of optimal parameters are obtained for the numerical model. The numerical results of net deformation, drag force and flow velocity around the fishing net agree well with the experimental data. It is indicated that present numerical model can simulate the interaction between flow and fishing net accurately.

Compared with the numerical results without considering FSI, present numerical model is considered to present more accurate results: (i) the flow-velocity reduction is larger downstream from the fishing net, (ii) the diversion of flow direction around the fishing net is more obvious, (iii) the deformation of the fishing net is less serious and (iv) the drag force is a bit smaller at approximately 7-8% for different incoming velocities.

This study forms a foundation for studying the flow field and hydrodynamic characteristics of net cages and offers some valuable information for the aquaculture industry.

## Acknowledgments

This work was financially supported by the National Natural Science Foundation (NSFC) Project Nos. 51239002, 51221961 and 51109022 and the Fundamental Research Funds for the Central Universities (DUT13LK55).