The Effect of Longitudinal Rails on an Air Cavity Stepped Planing Hull

: The use of ventilated hulls is rapidly expanding. However, experimental and numerical analyses are still very limited, particularly for high-speed vessels and for stepped planing hulls. In this work, the authors present a comparison between towing tank tests and CFD analyses carried out on a single-stepped planing hull provided with forced ventilation on the bottom. The boat has identical geometries to those presented by the authors in other works, but with the addition of longitudinal rails. In particular, the study addresses the effect of the rails on the bottom of the hull, in terms of drag, and the wetted surface assessment. The computational methodology is based on URANS equation with multiphase models for high-resolution interface capture between air and water. The tests have been performed varying seven velocities and six airﬂow rates and the no-air injection condition. Compared to ﬂat-bottomed hulls, a higher incidence of numerical ventilation and air–water mixing effects was observed. At the same time, no major differences were noted in terms of the ability to drag the ﬂow aft at low speeds. Results in terms of drag reduction, wetted surface, and its shape are discussed. experimental and numerical simulation. The maximum difference between the two curves is to the maximum velocity.


Introduction
The drag reduction is the main issue in limiting fuel consumption or enhancing the performance of High Speed Vessels (HSV). Traditionally the global drag is divided into two principal components: wave-making drag and viscous drag [1]. Under the name of viscous drag are included all the effects due to the fluid viscosity, i.e., all the effects that would be zero if the fluid is inviscid. The viscous part affects several phenomena, including the so-called form effect, the friction effect, the roughness, and the flat plate friction. According to the ship typology and its Froude number, the resistance can change from about 40% to 95% of the global resistance [2]. Over time, there have been numerous attempts to reduce the drag's viscous component, including the use of multihull, the Surface Effect Ships, the use of foils (to reduce the wet surface), or the use of hull ventilation [3]. From the first experiments to today, there have been many efforts to exploit the effect of the ventilation for drag reduction. Hull ventilation can be addressed both in a natural or forced manner. The presence of an air pocket can also have a lift role or only a lubricating effect. The shape or the mass of air can be in the form of cushion [4], layer [5], bubbles, or microbubbles [6].
In order to exploit natural ventilation of planing hulls, wedge boats and stepped bottom were theorized and studied, both experimentally and numerically, since the forerunners' work of Savitsky [7]. In fact, the planing hulls behave in a similar way to an

Methods
The chapter describes the geometric characteristics of the hull and the position of the step for the injection of the air and the geometric dimensions of the rails. In order to have a general idea of the line of the hull, a lines plan is proposed and the trend of the transversal section area is reported for the sake of completeness. The second part of this chapter introduces the simulation settings. The mesh and the refinement zones are presented and a complete description of the sequence of simulations conducted allows us to understand the approach used during this campaign of investigation. Particular attention concerns the methodology used for the wall-treatment inside the commercial CFD software and the initial boundary conditions.

The Model
As reported in [33], the original yacht is 18 m in length and whose principal dimensions are reported by Cucinotta et al. [22].
The bottom part of the hull was modified in order to have a step for injection of the air. A series of rails from the step to the transom was added. The model is fully described by [22] and for the sake of readability, the main dimensions are reported in Table 1.  Figure 1 shows the linesplane, with the presence of the step and its position with respect to the transom. The step is located 8 m away from the transom (1.33 m for the scale model) and the rails starting from this step until the end of the hull. The presence of the step produces a discontinuity in the transversal section area of about 0.3 m 2 .

The Model
As reported in [33], the original yacht is 18 m in length and whose principal dimensions are reported by Cucinotta et al. [22].
The bottom part of the hull was modified in order to have a step for injection of the air. A series of rails from the step to the transom was added. The model is fully described by [22] and for the sake of readability, the main dimensions are reported in Table 1.    The experimental campaign was carried out in the University of Naples facilities and the results are reported in [22]. Figures 2 and 3 and show, respectively, the bottom of the model and a picture from the tests in the towing tank. The experimental campaign was carried out in the University of Naples facilities and the results are reported in [22]. Figures 2 and 3 and show, respectively, the bottom of the model and a picture from the tests in the towing tank.

Simulation Settings
The approach used for solving the fluid dynamics around the hull is the Computational Fluid Dynamic (CFD in this paper). The main goal of simulation is to solve the pressure and velocity field around the hull, the volume fraction (the interface between the air and the water), and, finally, the relative position of the boat in terms of trim and sinkage. In order to capture all characteristics of the fluid field in terms of pressure and velocity, a URANSE (Unsteady RANS) approach was chosen. This approach allows to use the Navier-Stokes equations to solve all the unknowns of fluid. The finite volume method is implemented in the commercial software STAR CCM+ [34] used for this campaign of simulations. A crucial parameter is the time-step and the discretization order of the scheme used (in this simulation 2nd order) for the time-marching solution. The ITTC [35] suggests using a time-step in a range given by a function of the length of the hull and velocity of the boat, according to the formula ∆ = 0.01 ÷ 0.05 The experimental campaign was carried out in the University of Naples facilities and the results are reported in [22]. Figures 2 and 3 and show, respectively, the bottom of the model and a picture from the tests in the towing tank.

Simulation Settings
The approach used for solving the fluid dynamics around the hull is the Computational Fluid Dynamic (CFD in this paper). The main goal of simulation is to solve the pressure and velocity field around the hull, the volume fraction (the interface between the air and the water), and, finally, the relative position of the boat in terms of trim and sinkage. In order to capture all characteristics of the fluid field in terms of pressure and velocity, a URANSE (Unsteady RANS) approach was chosen. This approach allows to use the Navier-Stokes equations to solve all the unknowns of fluid. The finite volume method is implemented in the commercial software STAR CCM+ [34] used for this campaign of simulations. A crucial parameter is the time-step and the discretization order of the scheme used (in this simulation 2nd order) for the time-marching solution. The ITTC [35] suggests using a time-step in a range given by a function of the length of the hull and velocity of the boat, according to the formula ∆ = 0.01 ÷ 0.05

Simulation Settings
The approach used for solving the fluid dynamics around the hull is the Computational Fluid Dynamic (CFD in this paper). The main goal of simulation is to solve the pressure and velocity field around the hull, the volume fraction (the interface between the air and the water), and, finally, the relative position of the boat in terms of trim and sinkage. In order to capture all characteristics of the fluid field in terms of pressure and velocity, a URANSE (Unsteady RANS) approach was chosen. This approach allows to use the Navier-Stokes equations to solve all the unknowns of fluid. The finite volume method is implemented in the commercial software STAR CCM+ [34] used for this campaign of simulations. A crucial parameter is the time-step and the discretization order of the scheme used (in this simulation 2nd order) for the time-marching solution. The ITTC [35] suggests using a time-step in a range given by a function of the length of the hull and velocity of the boat, according to the formula ∆t = 0.01 ÷ 0.05 l v in which l is the length and v the velocity of the boat. A complete uncertainty analysis that comprises the evaluation of time step size, Courant number, grid dimension, and convergence ratio was carried out according to [35]. The detailed results are reported in [31].
The complexity of the simulation is due to the interaction between water and air in two different ways (interface of the sea, injection of the air in the bottom of the hull). Both fluids are treated as incompressible.
The multiphase simulation, with the VoF (Volume of Fluid) scheme, is used and thanks to the High-resolution interface capturing scheme (HRIC), a clear interface is captured during the entire duration of the simulation. A surface tension equal to 0.072 N/m between the two fluids is imposed. All the solver settings are the same as the ones used by Cucinotta et al. [31]. However, for the sake of clarity, Table 2 reports all the settings used during the simulation. In addition to the multiphase problem, the simulation also solves the motion of the boat. The approach used is the overset mesh. It helps to keep the quality of mesh very high. In order to use this procedure, the model was divided into two different regions: the background region (the fixed region in the space) and the overset region. The latter is free to move along the Z-axis (axis perpendicular to initial free-surface) and it is free to rotate itself along the Y-axis. The method used for solving the motion of the boat is the DFBI (Dynamic Fluid Body Interaction). It allows to solve the flow around a rigid body and simultaneously the motion of the rigid body caused by external forces and the forces induced by the flow (viscous forces, pressure forces). The equations of motion are the classic rigid body equations. For each step, the Navier-Stokes equations are solved considering the volume fraction under the hull by means VoF model. The viscous effects and the pressure field are integrated over the vessel in order to obtain the value of forces and moments. In function of these quantities, the solver applies the rigid body equations considering as center of motion the center of gravity of the boat. Finally, considering all the effects solved thanks to Navier-Stokes equations and Vof method (turbulence, air under the hull, waves, pressure field around the hull, viscous effects), the position, velocity, and acceleration of the boat is updated.
The geometry was defined with the commercial CAD software and all lines of the hull were defined with the higher precision possible (NURBS modeling). The ITTC [35] defines a series of guidelines in order to do a virtual towing tank. All these dimensions are shown in Figure 4.
The geometry was defined with the commercial CAD software and all lines of the hull were defined with the higher precision possible (NURBS modeling). The ITTC [35] defines a series of guidelines in order to do a virtual towing tank. All these dimensions are shown in Figure 4.  Figure 5 shows the background block with all the boundary conditions defined inside the solver. Each boundary condition is a field function that involves the velocity and pressure component of the fluid. The inlet has a Velocity Field function, it considers the position of the free surface and the velocity of the boat. The free surface is an unknown quantity that is solved each time step. The hidden surface is defined as symmetry in order to simulate only half of the domain. The outlet surface has a pressure boundary condition. It helps to update the free surface also in the outlet. In the overset region, the external surfaces are the ones linked with the background region. A zone where it is possible to have problems during the simulation is the overlapping zone between overset and background. There is a gap of elements where the solver exchanges the results of the corresponding zones; it is crucial to have a size of mesh very similar in this gap to avoid overflow errors. The blue surface is the symmetry area and the hull in grey is a wall no-slip surface. The step for injection of the air is a wall in no-slip condition (simulation is without air injection) and Inlet mass flow when the simulation involves the air injection.

Boundary Conditions and Mesh
The trimmed cartesian mesh process is used for this simulation. This approach allows defining different blocks for refinement. As shown in Figure 6, the refinement is conducted to keep a low aspect ratio and capture physics phenomena of different zones of  Figure 5 shows the background block with all the boundary conditions defined inside the solver. Each boundary condition is a field function that involves the velocity and pressure component of the fluid. The inlet has a Velocity Field function, it considers the position of the free surface and the velocity of the boat. The free surface is an unknown quantity that is solved each time step. The hidden surface is defined as symmetry in order to simulate only half of the domain. The outlet surface has a pressure boundary condition. It helps to update the free surface also in the outlet.

Boundary Conditions and Mesh
The geometry was defined with the commercial CAD software and all lines of the hull were defined with the higher precision possible (NURBS modeling). The ITTC [35 defines a series of guidelines in order to do a virtual towing tank. All these dimensions are shown in Figure 4.  Figure 5 shows the background block with all the boundary conditions defined inside the solver. Each boundary condition is a field function that involves the velocity and pres sure component of the fluid. The inlet has a Velocity Field function, it considers the posi tion of the free surface and the velocity of the boat. The free surface is an unknown quan tity that is solved each time step. The hidden surface is defined as symmetry in order to simulate only half of the domain. The outlet surface has a pressure boundary condition It helps to update the free surface also in the outlet. In the overset region, the external surfaces are the ones linked with the background region. A zone where it is possible to have problems during the simulation is the overlap ping zone between overset and background. There is a gap of elements where the solve exchanges the results of the corresponding zones; it is crucial to have a size of mesh very similar in this gap to avoid overflow errors. The blue surface is the symmetry area and the hull in grey is a wall no-slip surface. The step for injection of the air is a wall in no-slip condition (simulation is without air injection) and Inlet mass flow when the simulation involves the air injection.

Boundary Conditions and Mesh
The trimmed cartesian mesh process is used for this simulation. This approach allows defining different blocks for refinement. As shown in Figure 6, the refinement is con ducted to keep a low aspect ratio and capture physics phenomena of different zones o In the overset region, the external surfaces are the ones linked with the background region. A zone where it is possible to have problems during the simulation is the overlapping zone between overset and background. There is a gap of elements where the solver exchanges the results of the corresponding zones; it is crucial to have a size of mesh very similar in this gap to avoid overflow errors. The blue surface is the symmetry area and the hull in grey is a wall no-slip surface. The step for injection of the air is a wall in no-slip condition (simulation is without air injection) and Inlet mass flow when the simulation involves the air injection.
The trimmed cartesian mesh process is used for this simulation. This approach allows defining different blocks for refinement. As shown in Figure 6, the refinement is conducted to keep a low aspect ratio and capture physics phenomena of different zones of interest. It is evident that a clear refinement is inside the overset block, in the free-surface zone, and the area around the boat with the Kelvin triangle ( Figure 7).
interest. It is evident that a clear refinement is inside the overset block, in the free-surface zone, and the area around the boat with the Kelvin triangle ( Figure 7).  The hull has meshed with particular attention to the immersed zone ( Figure 8). The rails were captured by the mesh with a local refinement, this increases the number of elements of the mesh. The total number of cells for the overset zone is 4.57 Mln and for the background zone is 0.49 Mln. interest. It is evident that a clear refinement is inside the overset block, in the free-surface zone, and the area around the boat with the Kelvin triangle ( Figure 7).  The hull has meshed with particular attention to the immersed zone ( Figure 8). The rails were captured by the mesh with a local refinement, this increases the number of elements of the mesh. The total number of cells for the overset zone is 4.57 Mln and for the background zone is 0.49 Mln. The hull has meshed with particular attention to the immersed zone ( Figure 8). The rails were captured by the mesh with a local refinement, this increases the number of elements of the mesh. The total number of cells for the overset zone is 4.57 Mln and for the background zone is 0.49 Mln.
J. Mar. Sci. Eng. 2021, 9, x FOR PEER REVIEW interest. It is evident that a clear refinement is inside the overset block, in the fr zone, and the area around the boat with the Kelvin triangle ( Figure 7).  The hull has meshed with particular attention to the immersed zone (Figu rails were captured by the mesh with a local refinement, this increases the num ments of the mesh. The total number of cells for the overset zone is 4.57 Mln a background zone is 0.49 Mln. The wall surface of the hull has meshed with a prism layer of elements that allow capturing the boundary layer appropriately. The parameter that the authors keep unde control is the Wall y+ that for all simulations is always between 30 and 200 ( Figure 9). Finally, thanks to CFD, it is possible to know the area of hull wetted by the water o by the air. It is possible to solve this quantity considering a cut-off method on the volum fraction. The wetted surface by the water is calculated defining an interval of volume frac tion where the hull is considered wet by the water. In this work, the chosen threshold wa a volume fraction comprises between 0.5 and 1. The solver integrates the area where ther are these values of water fraction and so evaluate the amount of hull wetted by the water

Simulations Campaign
The simulations were conducted in two different parts. The first part involves th construction of the resistance curve in function of the velocity without the injection of th air. These results will be compared with experimental tests. During the simulation, th trend of the main physic quantities (drag, trim, and sinkage) was plotted and when thes curves reached a stable condition, the simulation was over.
The velocities under investigation are the same as the experimental tests conducted and reported by Cucinotta et al. [22]. From the first campaign, the recorded pressure and velocity fields of the fluid and the relative position of the yacht in terms of trim and sink age to each Froude number were recorded. These steady states became the initial condi tions of the transient simulation of the second campaign. This methodology helps the con vergence and the speed of the calculus.
During the first campaign of simulation, an uncertainty analysis has been conducted following the methodology proposed by ITTC [35]. The total uncertainty can be divided in three different terms: iterative (UI), grid (UG), and time-step uncertainties (UT). Stern e al. [36] suggested that the value of uncertainty caused by the grid is an order of magnitud greater than other ones and it is recommended for each simulation to evaluate this valu of uncertainty. The method consists in evaluating several n-th physic quantities ( , and ) to three different meshes (1 fine mesh, 2 medium mesh, 3 coarse mesh) with a constant ratio between number of cells ( ) (see Table 3). In this case, the physic quanti ties assessed are Drag, Trim, Sinkage, and Wetted Surface. For each parameter chosen, i can be calculated the difference between medium-fine mesh (2) and between coarse-me dium mesh (3). Thanks to these values, it is possible to evaluate the convergence ratio (4) and the order of accuracy (5). Finally, thanks to CFD, it is possible to know the area of hull wetted by the water or by the air. It is possible to solve this quantity considering a cut-off method on the volume fraction. The wetted surface by the water is calculated defining an interval of volume fraction where the hull is considered wet by the water. In this work, the chosen threshold was a volume fraction comprises between 0.5 and 1. The solver integrates the area where there are these values of water fraction and so evaluate the amount of hull wetted by the water.

Simulations Campaign
The simulations were conducted in two different parts. The first part involves the construction of the resistance curve in function of the velocity without the injection of the air. These results will be compared with experimental tests. During the simulation, the trend of the main physic quantities (drag, trim, and sinkage) was plotted and when these curves reached a stable condition, the simulation was over.
The velocities under investigation are the same as the experimental tests conducted and reported by Cucinotta et al. [22]. From the first campaign, the recorded pressure and velocity fields of the fluid and the relative position of the yacht in terms of trim and sinkage to each Froude number were recorded. These steady states became the initial conditions of the transient simulation of the second campaign. This methodology helps the convergence and the speed of the calculus.
During the first campaign of simulation, an uncertainty analysis has been conducted following the methodology proposed by ITTC [35]. The total uncertainty can be divided in three different terms: iterative (U I ), grid (U G ), and time-step uncertainties (U T ). Stern et al. [36] suggested that the value of uncertainty caused by the grid is an order of magnitude greater than other ones and it is recommended for each simulation to evaluate this value of uncertainty. The method consists in evaluating several n-th physic quantities (S n1 , S n2 , and S n3 ) to three different meshes (1 fine mesh, 2 medium mesh, 3 coarse mesh) with a constant ratio between number of cells (r k ) (see Table 3). In this case, the physic quantities assessed are Drag, Trim, Sinkage, and Wetted Surface. For each parameter chosen, it can be calculated the difference between medium-fine mesh (2) and between coarse-medium mesh (3). Thanks to these values, it is possible to evaluate the convergence ratio R G (4) and the order of accuracy P G (5).

Value of Convergence Ratio
Type of Convergence Oscillatory divergence Finally, with these two values it is possible to obtain the uncertainty with the Richardson formulation (6).
The main objective of this part is to define the reduction of resistance with the injection of the air. The flow rates analyzed were generally ranging from 5500 L/min (liter per minute) to 10,500 L/min, with a step of 1000 L/min. The Froude number studied were between 0.64 and 0.89 and from 5500 to 10,500 L/min for the Froude number between 1.02 and 1.36. Only for low Froude numbers the same high flowrate was not tested, since the difference between low and high flowrate became negligible.
The total number of conditions analyzed is 28. Another important objective of this campaign of simulation is the distribution of the air under the hull during the the air injection. The conditions that were simulated are briefly summarized in Table 4.

Hardware
All the operational phases described were conducted on a workstation with an Intel Xeon 2 GHz, 2 CPU with 16 core, 36 GB of memory, and a Nvidia Quadro M5000 8 GB.

Results of the First Part of the Simulation
The first part of the simulation concerns the model without air injection. The curves of resistance, trim, and sinkage were compared with the ones of experimental tests. In adding to this information, the other purpose of this block of simulations is to have a starting point for the next campaign of simulation. Figure 10 shows the comparison between experimental and numerical simulation. The maximum difference between the two curves is to the maximum velocity. tilation effect. This effect is caused by two main reason: the VoF method, at high Fn and with an overset mesh [37], and the planing hull characteristics (the acute angle caused by the intersection between the hull and the free-surface) [38]. Even if this phenomenon is known and studied, there are very few specific studies on a problem as complex as that relating to a cavity of air injected at high speed under a surface equipped with dynamic motions and with a mesh overset. This effect produces localized areas under the hull with the presence of the air which in reality does not exist. This presence of the air could lead to a reduction of the drag with respect to experimental tests. In this case the effect is localized and produces a difference lower than 7%. Another two crucial parameters compared are the trim and the sinkage. Figure 11 shows the behavior of the trim of the model during the different velocities. The trend of the curves is similar, with peaks of the CFD model more marked than the experimental test. Trim prediction is always the most challenging issue because it depends not only on the pressure but also on its distribution and, in particular, in this case, it is strictly related to the wetted length. However, even if the trend tends to change about Fn 0.8, the maximum difference is lower than 10%. In Figure 12 the sinkage comparison is reported, the two curves have the same trend. Usually, the increase of velocity could lead to a higher probability of numerical ventilation effect. This effect is caused by two main reason: the VoF method, at high Fn and with an overset mesh [37], and the planing hull characteristics (the acute angle caused by the intersection between the hull and the free-surface) [38]. Even if this phenomenon is known and studied, there are very few specific studies on a problem as complex as that relating to a cavity of air injected at high speed under a surface equipped with dynamic motions and with a mesh overset. This effect produces localized areas under the hull with the presence of the air which in reality does not exist. This presence of the air could lead to a reduction of the drag with respect to experimental tests. In this case the effect is localized and produces a difference lower than 7%.
Another two crucial parameters compared are the trim and the sinkage. Figure 11 shows the behavior of the trim of the model during the different velocities. The trend of the curves is similar, with peaks of the CFD model more marked than the experimental test. Trim prediction is always the most challenging issue because it depends not only on the pressure but also on its distribution and, in particular, in this case, it is strictly related to the wetted length. However, even if the trend tends to change about Fn 0.8, the maximum difference is lower than 10%. In Figure 12 the sinkage comparison is reported, the two curves have the same trend.
The CFD results of this campaign confirm that the settings used for the Virtual Towing Tank allow having results very similar to experimental tests. The comparison highlights the same trend for the three different quantities reported. Another important physic quantity is the wetted surface of the boat. This quantity is not a result of the experimental tests, but it can be evaluated with CFD simulations. Figure 13 shows the ratio between the wetted surface area in static condition and during the corresponding speed. The graph highlights a drastic decrease of the wetted surface between 0.7 and 1 of the Froude number. In this condition, the wetted surface is about 65% of the static one. At the Froude number of 0.7, the hull starts to planing so there is a decreasing of the wetted surface. The complete planing is reached at a Froude number of about 1.0.  The CFD results of this campaign confirm that the settings used for the Virtual Towing Tank allow having results very similar to experimental tests. The comparison highlights the same trend for the three different quantities reported. Another important physic quantity is the wetted surface of the boat. This quantity is not a result of the experimental tests, but it can be evaluated with CFD simulations. Figure 13 shows the ratio between the wetted surface area in static condition and during the corresponding speed. The graph highlights a drastic decrease of the wetted surface between 0.7 and 1 of the Froude number. In this condition, the wetted surface is about 65% of the static one. At the Froude number of 0.7, the hull starts to planing so there is a decreasing of the wetted surface. The complete planing is reached at a Froude number of about 1.0.  The CFD results of this campaign confirm that the settings used for the Virtual Towing Tank allow having results very similar to experimental tests. The comparison highlights the same trend for the three different quantities reported. Another important physic quantity is the wetted surface of the boat. This quantity is not a result of the experimental tests, but it can be evaluated with CFD simulations. Figure 13 shows the ratio between the wetted surface area in static condition and during the corresponding speed. The graph highlights a drastic decrease of the wetted surface between 0.7 and 1 of the Froude number. In this condition, the wetted surface is about 65% of the static one. At the Froude number of 0.7, the hull starts to planing so there is a decreasing of the wetted surface. The complete planing is reached at a Froude number of about 1.0. At the last value of this campaign of velocity, as mentioned above, an uncertainty analysis was conducted in order to evaluate the influence of the grid to the final result. At the value of Fn of 1.36, three different meshes were investigated with the number of cells in ratio of 1.4 (as suggested by ITTC [35]). Table 5 shows the results of the physical quan- At the last value of this campaign of velocity, as mentioned above, an uncertainty analysis was conducted in order to evaluate the influence of the grid to the final result. At the value of Fn of 1.36, three different meshes were investigated with the number of cells in ratio of 1.4 (as suggested by ITTC [35]). Table 5 shows the results of the physical quantities investigated to the three different meshes. For all the parameters under investigation, the method used suggests a monotonic convergence (Table 6). Thanks to the monotonic convergence, the safety factor of the Richardson expression (6) as suggested by Stern [36] can be equal to 1. The same table shows that the maximum uncertainty is about 5.1% for the mesh used in this campaign of simulations. The maximum uncertainty is for the value of drag.

The Second Campaign of Simulation
In this part of the simulation, for each number of Froude simulated in the first campaign, the virtual tests were performed adding the air injection. Also in this case, the numerical simulations were compared with experimental tests and in adding, thanks to the Virtual Towing tank, the wetted surface was reported for each condition of air flow. Figure 14 shows the first two flow rates with the relative curves of experimental tests. In both cases, the CFD results of air injection have a trend lower than experimental tests. In general, the air injection to these flow rates produces a decrease of resistance in the magnitude of 1% on average for the experimental tests. For CFD results, this decrease of resistance, on average, is in the magnitude of 5%. The difference between experimental tests and numerical simulations is lower than 15% in all the points. The difference between experimental tests and numerical simulations could be caused by the great mixture of air and water to these values of flow rate. The VoF scheme, with the HRIC model tends to underestimate the mix of air and water in the rails and underestimate the wetted surface, introducing a more positive effect on drag than experimental tests. Thanks to CFD it is possible to evaluate the differences in terms of reduction of the wetted surface with the air injection. The reduction of the wetted surface allows decreasing the frictional resistance of the model. On average, the reduction of wetted surface for 5500 L/min and 6500 L/min is about 60%. Figure 15 shows the trend to different velocities for both the flow rates. The first two flow rates investigated suggest that no great difference there is between 5500 L/min and 6500 L/min. The air under the hull probably is not entirely developed and the effect in terms of drag is very low in both flow rates. magnitude of 1% on average for the experimental tests. For CFD results, this decrease of resistance, on average, is in the magnitude of 5%. The difference between experimental tests and numerical simulations is lower than 15% in all the points. The difference between experimental tests and numerical simulations could be caused by the great mixture of air and water to these values of flow rate. The VoF scheme, with the HRIC model tends to underestimate the mix of air and water in the rails and underestimate the wetted surface, introducing a more positive effect on drag than experimental tests. Figure 14. Experimental and CFD results to air injection 5500 l/min and 6500 L/min. Thanks to CFD it is possible to evaluate the differences in terms of reduction of the wetted surface with the air injection. The reduction of the wetted surface allows decreasing the frictional resistance of the model. On average, the reduction of wetted surface for 5500 L/minand 6500 L/minis about 60%. Figure 15 shows the trend to different velocities for both the flow rates. The first two flow rates investigated suggest that no great difference there is between 5500 L/minand 6500 L/min. The air under the hull probably is not entirely developed and the effect in terms of drag is very low in both flow rates.  Figure 16 shows the drag curves for 7500 L/minand 8500 L/minflow rates and the relative ones without air injection. In both cases, the CFD results of air injection have a trend lower than experimental tests. In general, the air injection to flow rate of 7500 L/min produces a decrease of resistance in the magnitude of 2% on average for the experimental tests. For CFD results, this decrease is about 7%. In the case of 8500 L/min, the reduction is about 8% for experimental tests and 11% for CFD results. The difference between experimental tests and numerical simulations is lower than 15% in all the points.  Thanks to CFD it is possible to evaluate the differences in terms of reduction of the wetted surface with the air injection. The reduction of the wetted surface allows decreasing the frictional resistance of the model. On average, the reduction of wetted surface for 5500 L/minand 6500 L/minis about 60%. Figure 15 shows the trend to different velocities for both the flow rates. The first two flow rates investigated suggest that no great difference there is between 5500 L/minand 6500 L/min. The air under the hull probably is not entirely developed and the effect in terms of drag is very low in both flow rates.  Figure 16 shows the drag curves for 7500 L/minand 8500 L/minflow rates and the relative ones without air injection. In both cases, the CFD results of air injection have a trend lower than experimental tests. In general, the air injection to flow rate of 7500 L/min produces a decrease of resistance in the magnitude of 2% on average for the experimental tests. For CFD results, this decrease is about 7%. In the case of 8500 L/min, the reduction is about 8% for experimental tests and 11% for CFD results. The difference between experimental tests and numerical simulations is lower than 15% in all the points.  Figure 16 shows the drag curves for 7500 L/min and 8500 L/min flow rates and the relative ones without air injection. In both cases, the CFD results of air injection have a trend lower than experimental tests. In general, the air injection to flow rate of 7500 L/min produces a decrease of resistance in the magnitude of 2% on average for the experimental tests. For CFD results, this decrease is about 7%. In the case of 8500 L/min, the reduction is about 8% for experimental tests and 11% for CFD results. The difference between experimental tests and numerical simulations is lower than 15% in all the points.  Figure 17 shows the wetted surface at different Fn for flow rates of 7500 and 8500 L/min. The decrease of the wetted surface area is respectively 62% and 65%. To the flow rate of 8500 L/minthe decrease of drag is appreciable in experimental tests and in CFD simulations. The flow of air under the hull starts to cover a higher percentage of bottom of the hull. This value of flow rate seems to be the turning point, lower than this value the effects are negligible.  Figure 17 shows the wetted surface at different Fn for flow rates of 7500 and 8500 L/min. The decrease of the wetted surface area is respectively 62% and 65%. To the flow rate of 8500 L/min the decrease of drag is appreciable in experimental tests and in CFD simulations. The flow of air under the hull starts to cover a higher percentage of bottom of the hull. This value of flow rate seems to be the turning point, lower than this value the effects are negligible. Figure 17 shows the wetted surface at different Fn for flow rates of 7500 and 8500 L/min. The decrease of the wetted surface area is respectively 62% and 65%. To the flow rate of 8500 L/minthe decrease of drag is appreciable in experimental tests and in CFD simulations. The flow of air under the hull starts to cover a higher percentage of bottom of the hull. This value of flow rate seems to be the turning point, lower than this value the effects are negligible.  Figure 18 shows the drag curves for 9500 L/minand 10500 L/minflow rates and the relative ones without air injection. In this case, the CFD curves have a trend lower than the experimental ones, but the difference is reduced to lower than 10%. In general, the air injection to flow rate of 9500 L/minproduces a decrease of resistance in the magnitude of 11% on average for the experimental tests. For CFD results, this decrease is about 15%. For a flow rate of 10500 L/minthe decrease is 13% for the experimental test and 17% for CFD results.   Figure 18 shows the drag curves for 9500 L/min and 10,500 L/min flow rates and the relative ones without air injection. In this case, the CFD curves have a trend lower than the experimental ones, but the difference is reduced to lower than 10%. In general, the air injection to flow rate of 9500 L/min produces a decrease of resistance in the magnitude of 11% on average for the experimental tests. For CFD results, this decrease is about 15%. For a flow rate of 10,500 L/min the decrease is 13% for the experimental test and 17% for CFD results. Figure 17 shows the wetted surface at different Fn for flow rates of 7500 and 8500 L/min. The decrease of the wetted surface area is respectively 62% and 65%. To the flow rate of 8500 L/minthe decrease of drag is appreciable in experimental tests and in CFD simulations. The flow of air under the hull starts to cover a higher percentage of bottom of the hull. This value of flow rate seems to be the turning point, lower than this value the effects are negligible.  Figure 18 shows the drag curves for 9500 L/minand 10500 L/minflow rates and the relative ones without air injection. In this case, the CFD curves have a trend lower than the experimental ones, but the difference is reduced to lower than 10%. In general, the air injection to flow rate of 9500 L/minproduces a decrease of resistance in the magnitude of 11% on average for the experimental tests. For CFD results, this decrease is about 15%. For a flow rate of 10500 L/minthe decrease is 13% for the experimental test and 17% for CFD results.  Comparing to similar experiences carried out on the same hull but without the presence of the rails (i.e., [31,32]), in this study, the agreement between experimental and numerical data is lower for all the ventilated cases. This result probably depends on the difficulty of the interface models used to predict the real ventilated hull surface. In particular, the HRIC scheme tends to underestimate the air-water mixing in correspondence of hull discontinuity, and consequently to overestimate the ventilated surface, leading a reduction in the drag prediction. This numerical phenomenon must be properly evaluated as it was not a problem for flat bottom hulls. Figure 19 shows the trend to different velocities for flow rates to 9500 and 10,500 L/min. The decrease of the wetted surface area is respectively 68% and 69%. The flow rates are the higher ones and the air under the hull cover almost the entire bottom. Figure 20 shows the trend of drag reduction to different flow rates at the same Fn. The curves show that with the increase of the Froude number the curves tend to an high level of reduction of resistance. The curves show that to Fn number of 0.639 and 0.766, the decrease of resistance is almost equal to zero and also with the increase of the flow rates this behaviour remains the same. the HRIC scheme tends to underestimate the air-water mixing in correspondence of hull discontinuity, and consequently to overestimate the ventilated surface, leading a reduction in the drag prediction. This numerical phenomenon must be properly evaluated as it was not a problem for flat bottom hulls. Figure 19 shows the trend to different velocities for flow rates to 9500 and 10,500 L/min. The decrease of the wetted surface area is respectively 68% and 69%. The flow rates are the higher ones and the air under the hull cover almost the entire bottom.      Figure 20 shows the trend of drag reduction to different flow rates at the same Fn. The curves show that with the increase of the Froude number the curves tend to an high level of reduction of resistance. The curves show that to Fn number of 0.639 and 0.766, the decrease of resistance is almost equal to zero and also with the increase of the flow rates this behaviour remains the same.
To the Fn number of 0.894 the starting point of increase of percentage reduction is the flow rate of 6500 L/min. This starting point it is also visible for the Froude number equal to 1.149 but it is translated to the flow rate of 7500 L/min. To the highest values of Froude number, the curves seem flat and the flow rates does not affect the reduction of drag. To these values of Fn, there is not great need to push the flow rate until to 10,500 in order to obtain the desidered effects.  To the Fn number of 0.894 the starting point of increase of percentage reduction is the flow rate of 6500 L/min. This starting point it is also visible for the Froude number equal to 1.149 but it is translated to the flow rate of 7500 L/min. To the highest values of Froude number, the curves seem flat and the flow rates does not affect the reduction of drag. To these values of Fn, there is not great need to push the flow rate until to 10,500 in order to obtain the desidered effects.
The advantage of this hull grows with increasing speed, up to a peak reduction of 24% (at Fn 1.362 and 9500 L/min). For a speed of Fn < 1, the ventilation is not particularly effective and tends not to lead to benefits. For Fn ≥ 1, on the other hand, the benefits grow rapidly and the most advantageous air flow rates are those at 9500 L/min, above which the gain is reduced again.
In Figure 21 are reported all the wetted surfaces, for each Fn, in three different conditions, respectively, the no-air condition and the minimum and maximum tested flow rates. In the pictures, the water is represented by the color red and the air by the color blue. The figure highlights that to all the Froude numbers, the flow rate of 5500 L/min is not enough to wet the bottom of the hull completely. At this flow rate, for little Froude numbers, the air tends to escape from the sides, while for high Froude numbers there is a tendency of the air to be dragged up to the stern, remaining at the center of the boat. At the same time, at high speeds, air and water are less mixed. For this reason, for the Froude number higher than 1.27, the air does not escape from the sides of the boat and on the bilge area there is a channel of a stable water. From the same figure, it is possible to see that the ideal condition seems to be reached with the condition of 10,500 L/min at the Froude number of 1.15. In this case, the bottom of the hull and the area around the bilge are almost wholly wetted by the air without mixing with water.
bilge area there is a channel of a stable water. From the same figure, it is possible to see that the ideal condition seems to be reached with the condition of 10,500 L/minat the Froude number of 1.15. In this case, the bottom of the hull and the area around the bilge are almost wholly wetted by the air without mixing with water.
The longitudinal rails at low flow rates cannot retain the air and consequently are a disadvantage more than an advantage since they increase the wetted surface.    Flow rate 5500 L/min Flow rate 10,500 L/min Fn = 1.36

No-air
Flow rate 5500 L/min Flow rate 10,500 L/min

Conclusions
The experimental and numerical campaigns carried out on an ACS showed the high potentiality of the hull ventilation, also for planing hulls. The natural under pressure generates under the step, as in an airfoil lower surface, help the air to ventilate the hull without escaping too fast.
The hull reach a greatest advantage of about 24% of drag reduction. In the best conditions, it is clearly visible that the chain is wet (a sign that the air does not escape laterally) and that consequently the air channel completely fills the cavity up to the transom.
Event if it can be considered a good result, the rails under the hull have not been shown to give a great benefit compared to the hull without them. The idea for which they were designed, was that of channeling the flow of air towards the stern, limiting the lateral escaping, has proved to be ineffective. The airflow towards the stern in fact depends mainly on the Froude number. As Fn increases, the air tends to follow the travel direction, in the same way as the hull with the flat bottom [31]. On the contrary, the rails considerably increase the wet surface. The overall balance is therefore negative.
Even if the accordance between experimental and CFD results is good, the hull equipped with rails, compared to the flat ones showed in other papers, reveals a more significant error. This is probably due to greater difficulty in estimating the actual wetted surface in the area of the rails.
In this sense, the VoF scheme with the HRIC model, which well describes the wave surface into the air-water interface, tends to underestimate the mix of air and water in the rails and underestimate the wetted surface. Indeed, the real effect showed by experiments leads to the generation of microbubbles, instead of big, well-formed bubbles. The CFD campaign, therefore, showed that on this kind of boat, when the discontinuities under the hull increase, not only the hull has not great benefits in terms of drag reduction, but also that the numerical error increases.
The obtained results indicate interesting possibilities for future development in the study of biphasic phenomena under the hull, particularly in the presence of abrupt discontinuities. Different models than HRIC can be, in future, developed and applied to search for a better matching between experimental and numerical outcomes. The longitudinal rails at low flow rates cannot retain the air and consequently are a disadvantage more than an advantage since they increase the wetted surface.

Conclusions
The experimental and numerical campaigns carried out on an ACS showed the high potentiality of the hull ventilation, also for planing hulls. The natural under pressure generates under the step, as in an airfoil lower surface, help the air to ventilate the hull without escaping too fast.
The hull reach a greatest advantage of about 24% of drag reduction. In the best conditions, it is clearly visible that the chain is wet (a sign that the air does not escape laterally) and that consequently the air channel completely fills the cavity up to the transom.
Event if it can be considered a good result, the rails under the hull have not been shown to give a great benefit compared to the hull without them. The idea for which they were designed, was that of channeling the flow of air towards the stern, limiting the lateral escaping, has proved to be ineffective. The airflow towards the stern in fact depends mainly on the Froude number. As Fn increases, the air tends to follow the travel direction, in the same way as the hull with the flat bottom [31]. On the contrary, the rails considerably increase the wet surface. The overall balance is therefore negative.
Even if the accordance between experimental and CFD results is good, the hull equipped with rails, compared to the flat ones showed in other papers, reveals a more significant error. This is probably due to greater difficulty in estimating the actual wetted surface in the area of the rails.
In this sense, the VoF scheme with the HRIC model, which well describes the wave surface into the air-water interface, tends to underestimate the mix of air and water in the rails and underestimate the wetted surface. Indeed, the real effect showed by experiments leads to the generation of microbubbles, instead of big, well-formed bubbles. The CFD campaign, therefore, showed that on this kind of boat, when the discontinuities under the hull increase, not only the hull has not great benefits in terms of drag reduction, but also that the numerical error increases.
The obtained results indicate interesting possibilities for future development in the study of biphasic phenomena under the hull, particularly in the presence of abrupt discontinuities. Different models than HRIC can be, in future, developed and applied to search for a better matching between experimental and numerical outcomes.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.