Effects of Inclination Angles on the Hydrodynamics of Knotless Net Panels in Currents

: Knotless nets have been extensively used in sea cages and offshore ﬁsh farms; the explicit illustrations of local ﬂow patterns and subjected hydrodynamic loads have practical implications for the safe operations of aquaculture pens in open oceans. However, the effects of inclination angles of knotless nets on local ﬂow patterns have not yet been fully discussed. In this paper, the effects of inclination angles on local time-averaged and instantaneous ﬂow ﬁelds, velocity reductions and force coefﬁcients are investigated using a k - ω shear stress transport-improved delayed detached eddy simulation (SST-IDDES) high-ﬁdelity turbulence model. It is demonstrated that the inclination angles have dominant effects on the time-averaged velocity magnitudes around the net meshes. Secondly, the considerable ﬂow interactions amongst the transverse and longitudinal net meshes, as well as cross-ﬂow effects, are observed but weakened with the inclination angles. Thirdly, the proﬁles of velocity decelerations behind the net panels for 0 ◦ –45 ◦ are not as trivial as those for 45 ◦ –90 ◦ . The associations between the embedded mechanism of subjected loads and pressure ﬁelds are established to analyse the current loads on nets in small inclination angles. The inclination angles are ﬁnally elaborated into the area-averaged velocity reduction and screen force models, then the improved characteristics are validated against the experimental data.


Introduction
The expanding progress of moving mariculture cages from coastal areas to open oceans is inevitable nowadays. Diverse from closed or semienclosed sites, offshore aquacultures are ordinarily located in fully exposed regions. It is well known that open oceans are substantially isolated from lands or islands, and combined with the harsh sea states, hint at the potential threat to the safety of aquaculture infrastructures. Rigid steel frames are utilised so as to form a constant cultural volume and shield the structures from extreme waves and currents, while flexible nets, are however justified, fragile enough that the subjected excessive hydrodynamic loads can induce fatigues or breakages of twines. Thus, the safety issues of nets in the excitation of energetic waves and currents are emphasised continuously in offshore fish farming. Knotless nets with square meshes have been widely utilised in deep-water sea cages and offshore aquaculture farms presently, for instance, braided nylon fabric nets ( Figure 1a) and welded copper rigid nets (Figure 1b). The clarifications of flow distributions around nets, including time-averaged velocity properties, instantaneous turbulence fields and momentum reductions at the downstream vicinity, are theoretically and practically significant [1]. The predicted accuracies of subjected hydrodynamic loads, as well as the structural responses of sea cages, can be improved with a better knowledge of flow fields [2,3]. Further practical insights into flow fields are given as the correlations with fluid exchange rates and spreading patterns of nutrients across fish farms [4]. For example, the clarification of local flow fields can be used to assess the transportation of dissolved oxygen across net cages. Fish swimming behaviours are also closely associated with flow velocities inside cages. The design and operational guidelines of fishing nets or aquaculture farms require a better knowledge of hydrodynamic characteristics with various environmental scenarios and structural parameters. Practical implications, such as the improvement of trawl mouth expansions and the reduction in resistances of fishing gears, as well as biofouling, are closely concentrated with the hydrodynamic characteristics of nets. In summary, the flow around and wake behind nets, along with the subjected wave or current loads around nets, are of great significance for scientific research and engineering applications in aquaculture. The hydrodynamic characteristics of nets in currents, including local flow distributions and induced loads, are ordinarily investigated through numerical methods and physical experiments. With regard to flow fields, Løland [5] derived the velocity attenuation patterns behind a 2D net panel theoretically and qualitatively. Taking wake-interacting features behind twines into consideration, the predicted formulae for velocity reductions behind a full-scale net panel were initially given. Cha et al. [6] obtained a near-field wake behind twines of chain-link woven copper alloy nets using particle image velocimetry (PIV) measurements, thereby confirming the significant effects of the unique structural shape and the weaving patterns on wakes. Bi et al. [7] analysed the flow velocity reductions behind fabric net panels with various net solidities using PIV and acoustic Doppler velocimetry (ADV), presenting the dominant relevance with net solidities and inclination angles. These findings keep in accordance with the numerical study by Endresen et al. [8] and the experimental test of knotless nets by Føre et al. [9]. The turbulence characteristics and mass transport in the near-wake region of the net panel are assessed through the PIV-planar laser-induced fluorescence (PIV-PLIF) method [10]; herein, the turbulent diffusivity is found to theoretically be in reasonable agreement with the Taylor diffusivity.
In addition to experimental techniques, numerical methods, including fluid-structure interaction (FSI) and computational fluid dynamics (CFD) methodologies, are also capable of analysing flow patterns. Patursson et al. [11], Bi et al. [12] and Yang et al. [13] studied downstream velocity reductions of fixed net panels through the equivalence of porous media, and the accuracies of the numerical method were validated against the published experimental data. Chen and Christensen [14] and Cheng et al. [15] incorporated porous media modelling with the lumped mass method and Morison force model, enabling a two-way coupling scheme and expanding its applicability to flexible nets and cages. Later, a novel coupling algorithm comprising the Lagrangian approach for net panels and the Eulerian fluid domain was proposed in Martin et al. [16], then flow reductions behind nets, cage-to-cage effects and full-scale net hydrodynamics in steel-frame fish farms [17,18] were reported. Inspired by cruciform configurations in net meshes [19], the Reynolds-averaged Navier-Stokes (RANS)/large eddy simulation (LES) studies in Bi et al. [20], Zou et al. [21], Tang et al. [22] and Dutta et al. [23] were conducted to provide velocity profiles around cylinders, knots and cruciforms, then the impacts of mesh orientations, Reynolds numbers and so-called blockage ratios were also addressed. Benefiting from the advancements of computational abilities, a unique perspective for hydrodynamics opened up through the structural connections amongst circular cylinders, spheres, cruciforms and nets. Wang et al. [1] modelled net panels with 6 × 6 meshes using LES, then the distributions of turbulence kinetic energies and vorticity patterns over nets were firstly given. Furthermore, the lattice Boltzmann method was implemented and validated for the flow around net planes by Tu et al. [24]. The wake interactions behind planes were also justified through velocity profiles in this study, which are consistent with Løland [5] and Wang et al. [1].
On the other hand, the hydrodynamic loads on net panels are considerably correlated with materials, mesh shapes, structural parameters and biofouling, along with angles of attack (AOAs). In terms of metal net panels, Tsukrov et al. [25] and Cha et al. [6] measured hydrodynamic coefficients of woven, welded and chain-linked copper nets through water flume and PIV tests, respectively. Apart from solidity ratios and Reynolds numbers, the influences of weaving styles, surface roughness, knots and AOAs on the hydrodynamic coefficients of various metal nets were also discussed in detail. As for fabric nets, a comparative study between rectangular and cylindrical net panels was carried out to justify the close relevance of Reynolds numbers, solidities, inclinations and hydrodynamic loads [26]. With the same geometrical parameters, it was substantiated that the drag coefficient of a net bar is higher than a circular cylinder, and the drag coefficient of rectangular panels is also larger than that of cylindrical nets. Lader et al. [19] focused on the hydrodynamic properties of cruciform structures constituting net meshes, indicating that the sphere in a cruciform configuration can experience two to three times larger size than that without spheres. Following the concept of cylindrical cruciform elements forming nets, the drag coefficient of net panels was deduced from knots and twines, then the relevance of Reynolds numbers and solidity ratios were discussed through the experimental data [27]. For the long-term duration of aquaculture operations, the effects of biofouling were considered by Bi et al. [28], where it was found that accumulation can result in over 10 times larger loads on clean nets. With respect to net materials, hydrodynamic loads on the nylon and polyethylene net panels with twisted or knitted styles were examined through flume experiments [29,30], wherein the influences of monofilament or multifilament processing on the hydrodynamics were compared initially.
To the author's best knowledge, the hydrodynamics of nets have been investigated for decades, whereas the effects of inclination angles of rigid and fabric nets on local flow patterns have not been fully discussed. Secondly, the embedded mechanism of subjected hydrodynamic loads has not been addressed by the flow distributions. Thirdly, Kristiansen and Faltinsen [31] proposed a novel screen-type force model predicting the hydrodynamic drags and lifts, represented by the Fourier series and bounded at the drag or lift coefficients of θ = 0 • , 45 • and 90 • . The predicted drag coefficient equals 0 for the θ = 90 • , i.e., the incoming flow keeps parallel to the net panel. This estimation is, however, opposed to the physical observations of hydrodynamic drags of knotless net panels by Tang et al. [30] and Tang et al. [32], hinting at the downside of existing formulae.
Provided the adoption of high-fidelity CFD methods, it is feasible to introduce a wide range of inclination angles into the fitted formulae of the wake reduction model [1] and the screen-type force model [31], leading to further advancements in accuracies. In this research, the k-ω shear stress transport-improved delayed detached eddy simulation (SST-IDDES) turbulence model, enabling scale resolution at a lower computational cost than LES, was adopted to resolve the turbulence fields around net twines, while the far-field regions were modelled using RANS. The numerical framework, as well as the convergences of domains and grids, were successfully validated and justified in the author's former published work [2]. The effects of the inclination angles of smooth rigid net panels on time-averaged and instantaneous flow fields, downstream momentum reductions and hydrodynamic coefficients were investigated in pure currents; the inclination angles were subsequently introduced into the wake reduction model and the screen force model. Finally, the updated screen force model was applied and validated in an FSI framework for simulating fullscale rigid and fabric nets in viscous fluids, including the approximations of subjected hydrodynamic loads and momentum reductions.

Net Panel Models
As revealed by Wang et al. [1,2] and Tang et al. [27], the dependencies of hydrodynamic behaviours with net solidities and Reynolds numbers are almost identical for rigid and fabric nets. The detailed weaving patterns of twines and knots cannot be accurately captured or modelled during numerical simulations, otherwise the accuracies can be sacrificed considerably. Thus, only smooth rigid knotless net panels with square meshes are used as the targeted models in CFD modellings ( Figure 2). The investigated model scale was chosen as 6 × 6 meshes, achieving the mesh-number convergence of flow fields and drag force coefficients in the author's published work [1]. Under the premise of saving computing resources, the high-resolution CFD simulation of small-scale net panels can reflect flow fields and hydrodynamic characteristics around full-scale nets to the greatest extent, while ignoring net deformations in currents. The total of 65 computational cases were set to analyse the effects of inclination angles and structural and inflow properties on the hydrodynamic characteristics of net panels, as shown in Table 1. Following Føre et al. [9], the discrete sets of 0 • , 22.5 • , 45 • , 67.5 • and 90 • were utilised to represent a broad scope of inclination angles. According to the marine fishery practice, the broadly adopted case C1 was selected to analyse the local flow properties, of which structural and inflow parameters were the medians across the investigated ranges. Concerning wake attenuation behind the net panels and subjected hydrodynamic loads, the structural parameters of the twines kept constant, while the integrated influences of incoming velocities U 0 and angles of attack θ were examined in cases C1-C5. C1, C6-C9 and C10-C13 were considered for the integrated effects of bar diameters d and bar lengths l, along with θ, while keeping U 0 unaltered.

Incoming Bar Bar Reynolds Solidity Velocities Diameters Lengths
Number

Viscous Fluid Dynamics Solver
k-ω SST-IDDES model was utilised to model the viscous fluids around the net panels and the wake behind them. As a hybrid methodology of RANS/LES for enabling a scaleresolving strategy, its original intention was to reduce the considerable computational cost of LES modelling, that is, RANS is used to model Reynolds stress or turbulent kinetic energy on the wall and far-field, while the LES method is capable of resolving the Reynolds stress or turbulent kinetic energy for the large-scale eddy motions on the near-wall regions. In contrast with the strategy that adopts the RANS model globally [20,21], the turbulence disturbance at downstream wakes as well as the flow separations from near-field boundary layers of cylindrical twines can be detected more precisely. As the basis of the IDDES method, the three-dimensional continuity and Navier-Stokes (N-S) equations are written in the convective form as where U represents the velocity vector; ∇p and ρ denote the pressure gradient and the density of viscous fluids, respectively; g is the gravitational acceleration vector.
In the RANS part, the eddy viscosity in the diffusion term was enclosed using k-ω SST turbulence model [33]. The improving characters at accuracies and efficiencies of simulating flow separations from cylinders were dominant against the standard k-ω model. As for the IDDES method, Gritskevich et al. [34] enhanced the performance of the existing k-ω SST delayed detached eddy simulation (DDES) model without sacrificing the turbulence resolving in the vicinity of flow separations. IDDES method is derived by replacing the RANS length scale of the turbulent kinetic energy k equation with the IDDES length scale, while keeping the specific dissipation rate ω unaltered. The governing equations of k and ω in the SST-IDDES model read as [35] ∂ρk ∂t where F 1 and F 2 denote the SST blending functions, which read as follows: where d w represents the distance to the nearest wall. The turbulence production term P k and IDDES length scale in Equation (2) are given below: where C DES1 = 0.78, C DES2 = 0.61 and ∆ is described as where h max denotes the length scale of computational grids, equalling the maximum length of all grid edges. Moreover, the empiric blending function embedded in Equation (14) is The detailed procedures for calculating f dt , f b and f e can be found in the study by Gritskevich et al. [34]. All the cases were simulated using the open-source CFD toolbox OpenFOAM v1806. The fluid properties of CFD simulations are listed in Table 2. The transient solver pisoFoam accounting for pressure-velocity decoupling was used to solve the N-S equation. Considering balancing accuracies and stabilities, the convective terms were spatially discretised by the second-order linear-upwind stabilised transport (LUST) scheme with low-level numerical dissipation. The diffusive terms were treated through the second-order central difference scheme, while the temporal term advancement was processed with the second-order backward scheme. The time steps ∆t within unsteady simulations were determined as 0.005 d/U 0 , keeping the Courant-Friedrichs-Lewy (CFL) condition fluctuating at around 0.3-0.5. The pressure gradient was solved with the generalised geometric-algebraic multigrid (GAMG) method, while the momentum equation was solved through the Gauss-Seidel algorithm with additional smoothing for improved convergence.

Computational Domain, Grids and Boundary Conditions
As depicted in Figure 3a, the computational domain was determined as 16.0 l × 14.4 l × 14.4 l (length × width × height). The coordinate origin coincides with the geometrical centre of net panels, placed at 1/3 of the length and in the middle of the width and height directions. The inclination angle in this paper is defined as the intersection angle between the net plane and the gravity direction, as shown in Figure 3b. Verified by Wang et al. [2], this configuration ensures that the water flow can fully develop before reaching the panel, and minimises the irrelevant influence from side boundaries of the computational domain on the hydrodynamics of nets. The left boundary in the domain was specified as the inlet condition, following the isotropic initialisations of k and ω. The right-side boundary was assumed as the zero-pressure outlet condition, where U 0 , k and ω respect the zero-gradient Neumann boundary conditions. The top and side boundaries were appointed as symmetry planes, that is, the normal components of the velocity and the normal gradients of all other variables were zero. The bottom plane and the net twines were both assumed as non-slip walls, but only the turbulent viscosity close to the bottom surface was modelled by wall functions.
The computational grids were generated using SnappyHexMesh in OpenFOAM. Owing to the meshing difficulties around net twines and numerous cruciform substructures, the unstructured hexahedral body-fitted grids are more capable of preserving an acceptable grid quality. Overview of the computational grids around the C1 case with the θ = 22.5 • as well as three refinement zones close to net twines were also set to capture turbulence fluctuations around the structure. A total of 13-15 boundary layer grids attached to the near-field region were generated with the expansion rate 1.1, satisfying that the first layer is inside the laminar sublayer of the boundary layer, i.e., y + < 1. According to the grid independence tests of flow fields and hydrodynamic coefficients [1,2], it 2qs determined that the grid size in the far-field area was 3.80 mm, while the most refined grids in the boundary layer remained 1.01 × 10 −5 m. The resulting grids had 40-50 M points.

Data Analysis
The time-averaged drag coefficient C D and the root-mean-square lift coefficient C L [11] were utilised to analyse the hydrodynamic drags F D and lifts F L on net panels with the outline area A. They are written as follows:

Effects of Inclination Angles on Local Flow Around Net Panels
The influences of the inclination angles of net panels on local flow fields are illustrated through two slicing planes (P1 and P2) across the net panels in the normal and tangential directions ( Figure 4). For case C1 normal to the incoming flow, the nondimensional mean velocity magnitudes (|U|/U 0 ) simulated by the IDDES method keep in agreement with the LES observations by Wang et al. [1] (Figure 5). For the other cases, the inclination angles have dominant effects on the time-averaged velocity magnitudes in the tangential plane across the net panel, and the extent is obviously more pronounced than twine diameters and lengths. More specifically, there are considerable flow interactions amongst the transverse and longitudinal net meshes for the net panel perpendicular to the incoming current. The interacting effects are nevertheless weakened with larger AOAs, with the areas of local flow acceleration zones around twines presenting a shrinking tendency, while the low-speed regions are of larger dominances gradually. Until the net plane and the incoming flow are in parallel, the central regions of net meshes are fully occupied with downstream wakes, of which velocity magnitudes close to knots and net bars are almost decreased to zero.  In comparison to slice P1, the turbulence intensities of the flow fields at slice P2 are more significant, because viscous fluids are separating from cylindrical twines, resulting in an unsteady vortex shedding phenomenon behind nets in slice P2. Thus, the necessity of analysing instantaneous vorticity results in slice P2 has been emphasised to obtain a higher resolution of flow patterns. The modulus fields of the nondimensional instantaneous vorticity (|ω|d/U 0 ) at slice P2 of the net panels with θ = 0 • , 45 • and 90 • are depicted in Figure 6. For the case with θ = 0 • , it is demonstrated that the vortex street behind the twine at the net boundaries shifts towards the outer edge in a marginal manner (see the deviated angles δ in Figure 6). This has the implications that the flow patterns around nets show similar characteristics to those around bluff bodies, as well as the so-called cross-flow properties for the downstream wakes. The IDDES results have also been observed in the published numerical studies [1,24], in addition to theoretical deductions [5]. However, one should notice that the deviations of downstream streamlines towards outer edges are not prominent for the case with θ = 45 • . It can be attributed that the longitudinal layout of twines caused by the increase in inclination angles weakens the cross-flow. When θ = 90 • , that is, the net bars perpendicular to P2 are entirely maintained in a tandem and longitudinal configuration with a spacing of 11.4d, the vortex shedding cannot be seen, as in Figure 6c. Further comparisons were given by Mittal et al. [36] that the vortex separations behind the tandem cylinders are considerable at a mere 5.5d spacing. A possible explanation for the difference is that the flow interactions induced by twines and cruciforms in vicinities still affect local flow fields, even in the case of small inclination angles.

Effects of Inclination Angles on Velocity Reductions Behind Net Panels
In Figure 7, the very far-field wake regions of flowing past net panels are analysed through the normalised time-averaged velocity modulus extracted at x/d = 20 behind the net panel. The trivial influences from the outlet boundary are avoided, as stated by Wang et al. [1]. For the case with 0 • -45 • , the velocity around the mesh centre can restore the initial momentum, while the velocity reductions are primarily concentrated in the downstream area where twines are projected. Thus, the overall differences in velocity distributions are not significant. With regard to the case with θ = 45 • , the majority of wake regions are concentrated at the bottom of contours. This is attributed to the fact that the bottom vicinities remain much closer to the net panel, forming a dominant negative pressure zone and implying considerable momentum reductions. For θ = 45 • -90 • , the velocity attenuation regions are more concentrated in the area where the panel is projected to x/d = 20 in the incoming flow direction. This trend is also justified quantificationally in the physical measurement of velocity reductions behind a net panel by Patursson et al. [11]. Furthermore, the flow velocities around the mesh centre cannot be restored yet due to the intense momentum reduction. The discrepancies of distributing patterns of velocity contours for the cases with θ = 45 • -90 • are more significant against the cases with θ = 0 • -45 • . Therefore, it is of great necessity to carry out the comparative discussions between the velocity decreasing patterns for θ = 45 • -90 • and θ = 0 • -45 • in further investigations. Based on the flow field data using k-ω SST-IDDES method, the area-averaged velocity reduction (AAVR) model suggested by Wang et al. [1] is utilised to quantify the effects of inclination angles and other factors, such as Re and S n , on the momentum attenuation behind net panels. The definition of AAVR is written as where A denotes the area of the net panel projected to the rear plane of x/d = 20, (U/U 0 ) cell represents the nondimensional mean velocities at the dA dominated by the targeted cell centre and vicinity. Thus, the velocity reduction is estimated by combining the weight of each computational cell. The estimation quality can be improved compared to the commonly applied velocity reduction factor, assuming a constant value at a certain point in wake regions. The AAVRs of net panels perpendicular to currents with Re and S n were fitted in Wang et al. [1], then the velocity reduction estimated formulae for rigid smooth and fabric nets were proposed with the validations against experimental data. The results of the AAVRs of all cases in this paper are summarised in Table 3. It can be concluded that the AAVRs across θ = 0 • -45 • are below 5% with the suitable range of Re and S n ; thus, referring to the AAVR-predicted formulae [1] for nets normal to currents is acceptable, which is in agreement with the varying patterns of velocity contours with AOA ( Figure 7). In contrast, the AAVRs of the cases of θ = 45 • -90 • rapidly reduce to approximately 0.75, indicating the necessity of performing the nonlinear fittings between AAVRs and inclination angles in this scope. On the basis of the data in Table 3, nonlinear fitting is accomplished using the downhill simplex method [37]. As a result, the mean absolute percentage error (MAPE) between the Fourier series fitted formulae and the data in Table 3 remains 1.74%.  On account of the lack of experimental data on velocity reductions behind welded rigid net panels with various inclination angles, only the validations of AAVR models for fabric net panels were executed. Bi et al. [7] measured the mean velocities at the central wake regions of net panels with θ = 0 • -60 • so that the validations of the novel AAVR model could be conducted to show the improved features. Table 4 lists the validation results as well as the relative divergences AAVR−Exp. Exp.
(%). It is demonstrated that the relative divergences are below 5% within U 0 ≤ 0.226 m/s; hence, the predicted accuracies of the improved AAVR models are verified across the cases of a wider range of inclination angles.

Effects of Inclination Angles on Hydrodynamic Coefficients
As mentioned above, the accuracy of the simulated hydrodynamic coefficients of net panels deserves further advancements, especially at small AOAs. In order to improve the accuracies of existing models, the influences of wide-range inclination angles (0 • -90 • ) on C D and C L were studied following the high-resolution k-ω SST-IDDES method. As shown in Figure 8, the distributions of simulated C D and C L with AOAs over all cases are compatible with the laws revealed by Kristiansen and Faltinsen [31]. Nevertheless, the C D and C L for the case with θ = 90 • can be obtained as 0.08 and 0, respectively. As opposed to the typical formulae in Kristiansen and Faltinsen [31], this has the implication that C D of the cases with small inclination angles cannot be overlooked. This statement can be examined through the time-averaged pressure in slice P1 around the net panel with θ = 90 • (Figure 9). The pressure jump over the twines contributes to the hydrodynamic drag, which can be identified as the significant connection between hydrodynamic forces and flow fields. The simulated C L remaining around 0 is substantiated due to the highly symmetrical features of nets in the direction of lifts.  With respect to the corrections of existing screen force models, the nonlinear fittings amongst C D , C L and the inclination angles were accomplished using the downhill simplex method [37]. Then, the estimated hydrodynamic coefficient formulae of the novel screen force model (NSFM) were proposed based on the CFD data and bounded by C D and C L at θ = 0 • and 45 • , as Herein, C D,0 • , C D,45 • and C L,45 • for fabric nets were separately calculated using the polynomial-fitted and power-fitted formulae [2], as in Equation (22). The applicability and accuracies of the polynomial-fitted formula applied for C D as well as the power-fitted formula for C L were elaborated and justified in [2]. C D,0 • = a 1 U 0 + a 2 d + a 3 l + a 4 U 0 d + a 5 U 0 l + a 6 dl + a 7 U 2 0 + a 8 d 2 + a 9 l 2 + a 10 , (0.1 m/s < U 0 < 1 m/s, 0.5 mm < d < 3 mm, 0.01 m < l < 0.03 m, 50 < Re < 2970, 3.31% < S n < 51.0%) with the coefficients in Table 5: Table 5. Coefficients for the polynomial-and power-fitted formulae of C D and C L [2].  [31] to use two harmonics to characterise C D and C L , but the goodness of nonlinear fitting is insufficient if truncated after two coefficients. Thus, one more harmonic was introduced by Martin et al. [16] and Yao et al. [38] to evaluate drags and lifts, respectively. As opposed to the previous patterns, the harmonics and constants were revised extensively in the newly proposed formulae (Equations (20) and (21)) and truncated by C D,0 • , C D,45 • and C L,45 • . The divided-format function bounded by 45 • was also established simultaneously to fit the simulated data better, resulting in the MAPEs of C D and C L achieving 4.74% and 7.01% separately. It can be seen from Figure 8 that the fitting curves captured the majority of simulated data; thus, the effectivenesses of the nonlinear fittings were demonstrated from both qualitative and quantitative perspectives.

Applications and Validations of the Novel Screen Force Model
The newly proposed formulae of the screen force model derived from 6 × 6 meshes were applied to the full-scale net hydrodynamics solver in a fluid-structure coupling framework. Then, the simulated hydrodynamic loads on and momentum loss behind the full-scale net panels were compared against the published experimental data, showing the improved character of the NSFM.

The Fluid-Structure Coupling Method
Following the concept raised by Martin et al. [16], the net panel can be characterised by Lagrangian points coupled in an Eulerian fluid domain, modelling the physical fluidstructure interaction through an alternative fulfilment of the boundary conditions at the interface. A source term F net was incorporated into the right-hand side of Equation (1), describing the rapid pressure jump across the net panel, whereby the momentum losses of viscous fluids due to the presence of nets were identified and then estimated, observing the momentum conservation law of fluids. As an instance in Figure 10, it should be noted that one of the essential hypotheses is that momentum losses are assumed to be distributed isotropically around nets regardless of the presence of twines and knots. The uniformly distributed Lagrangian points are required, then the distance between the two ambient points is correlated with the cell size of the surrounding Eulerian grid. This dramatically improves the computational efficiency without evaluating the net geometrical details. This strategy can be naturally joined with screen force models, and similar simplifications were also utilised in previous studies of fluid-structure interactions [11,12,14,15]. The isotropic kernel D around x i was introduced to account for the F net and the dominated vicinity, which is depicted by two averaged split triangles of one computational cell. Each split triangle is related to a specified Lagrangian point x L . Implemented in the open-source CFD toolbox REEF3D, the detailed coupling procedure as well as the sketch of isotropic kernels can be found in the studies by Martin et al. [16] and Martin et al. [17]. In brief, F net is established by relating subjected external forces with the covered region, which can be written as where L i denotes the number of Lagrangian points in the dominated vicinity of the isotropic kernel D around x i [39]: For the fixed nets, the external forces at the Lagrangian point x L dominated by the triangular screen element (with the area of A s ) are the integral of the hydrodynamic forces f h , the gravity vector f g and the buoyancy f b , as

Validations of Hydrodynamic Forces on and Velocity Reductions behind Fixed Net Panels
In this section, the accuracy of subjected current forces as well as momentum loss behind fixed net panels are validated against the published experimental results. The original formulae within a classical screen force model (CSFM) derived by Kristiansen and Faltinsen [31] and refined by Martin et al. [16] are also utilised for further discussions. The net cases used in the validation part are illustrated in Table 6. Three sets of experimental data from Patursson et al. [11], Zhan et al. [26] as well as Føre et al. [9] were adopted to examine the performance and applicability of NSFM for various net panels. In the study by Patursson et al. [11], the subjected drag/lift forces and downstream velocity losses of a 1 m × 1 m knotless nylon net panel with S n = 0.184 were measured in a 37 m (length) × 3.66 m (width) × 2.44 m (depth) towing tank, and the net panel was fixed with the rigid frame beneath the free surface of 0.732 m. The velocity reductions were observed at the centreline 2.5 m behind the net panel. A wide range of net inclination angles and free-stream velocities was considered, especially for the net case forming a small AOA (N8) with the free-stream direction. In contrast, the planar net panels used in Zhan et al. [26] were configured with square diamond meshes, towed in a 204 m (length) × 6 m (width) × 3 m (depth) tank. With regard to the experimen-tal setup by Føre et al. [9], a larger net panel attached with the rigid frame of 1.215 × 0.985 m was utilised in a 10.5 m width and 5.6 m depth towing tank, and the evaluation point of the downstream velocities was located behind the net panel with a horizontal distance of 0.715 m. The lengths of numerical tanks across all net cases were reduced to improve computational efficiencies and discretised uniformly with a cell size of 0.05 m, based on the verification studies by Martin et al. [16]. The comparisons of the predicted hydrodynamic results by CSFM and NSFM against experimental measurements are illustrated in Table 7. As for hydrodynamic loads of N1-N3 with varying U 0 , it can be seen that the NSFM performs better than CSFM with relative differences less than 6%. However, the velocity reduction of N3 with the minimum U 0 is underestimated substantially (36.30%), implying that the predicted ability of NSFM for the case with marginal turbulence fluctuations also needs further improvements. For N4-N8 with the larger inclinations, the hydrodynamic drags, lifts along with velocity losses are exhibited with more significant relative discrepancies between the NSFM and experimental results, but the NSFM in the majority of cases still performs overwhelmingly better than CSFM. The average error magnitude is 17.82% for the new formula versus 23.48% for the previous CSFM. A possible explanation was claimed by You et al. [40] that so-called dual effects, described as the various relationships between hydrodynamic parameters and net solidities with small inclination angles, are dominant for the larger discrepancies. For the cases by Zhan et al. [26] and Føre et al. [9], the improved properties of the NSFM formulae in estimating hydrodynamic drags were demonstrated with less than 5% deviation. In conclusion, the NSFM of hydrodynamic drags for the cases of θ ≤ 45 • is justified, whereas the abilities for θ > 45 • and the decelerations of fluid while passing net panels desire more advancements in terms of small AOAs.

Conclusions
In this paper, the effects of a wide range of inclination angles of knotless net panels on local time-averaged and instantaneous flow fields, momentum decelerations of fluids, and hydrodynamic force coefficients are comprehensively investigated using a k-ω shear stress transport-improved delayed detached eddy simulation (SST-IDDES) high-fidelity turbulence model. The conclusions are drawn as follows: (1) The inclination angles have dominant effects on the time-averaged velocity magnitudes around net meshes. More specifically, considerable flow interactions amongst the transverse and longitudinal net meshes for the net panel perpendicular to the incoming current are observed. The flow interacting effects amongst net meshes are nonetheless weakened at larger inclination angles. It can be seen that local flow acceleration zones around twines present a shrinking tendency, while low-speed regions are of larger dominance gradually. (2) From the transient vorticity fields, the so-called cross-flow properties at downstream wakes proposed by Løland [5] can be validated from our numerical results. Nevertheless, the longitudinal layout of twines induced by the increase in the inclination angles of nets has weakened the cross-flow properties. The accuracies and applicabilities of the novel formulae with a broad scope of AOAs substantiate that the relative discrepancies of all validated cases are within 5% through the comparisons against experimental data. (5) The distributions of simulated C D and C L with inclination angles over all cases are compatible with the laws revealed by Kristiansen and Faltinsen [31]. Sensible differences can nevertheless be detected for the nets parallel to currents. This work has further elaborated Fourier-series formulae of screen force models, and the predicted hydrodynamics results by the improved and the previous formulae have been compared with the experimentally measured data. The improved features for hydrodynamic drags are justified, whereas the abilities for hydrodynamic lifts and the decelerations of fluid while passing net panels desire more advancements at small inclination angles.