A modification of the ITTC57 correlation line for shallow water

The ITTC57 correlation line, which is derived based on the assumption that the water in which ships advance is infinite deep and wide. However, for ships sailing in the waterway with limited water depth, the frictional resistance will be influenced leading to a decreasing accuracy of the prediction with this correlation line. In this study, a modification of the ITTC57 correlation line is proposed to correct the effects in very shallow water specifically for the flat area of the bottom of the ship. Under some assumptions, this area can be simplified to a 2D flat plate with a parallel wall close to it to study how the shallow water conditions of two interacting boundary conditions are affecting the flat plate friction coefficient. Computational fluid dynamics (CFD) calculations are applied to investigate how a friction line specifically in shallow water deviates from the conventional lines. Such deviations may severely affect the extrapolation of a ship model’s resistance to full scale and, therefore, the accuracy of ship’s performance prediction. Cases at ten Reynolds numbers from 105 to 109 are simulated on the 2D flat plate. Seven different distances between the flat plate and the parallel wall were chosen to generate various shallow water conditions, and consequently, a database including frictional resistance coefficients, Reynolds numbers and the distance between those two walls is built. Results indicate that thinner boundary layers are observed in shallow water conditions, and the scale effects which has a significant impact on resistance extrapolation are also observed. Furthermore, the assumption of the zero pressure gradients (ZPG) which is commonly used in deep water is no longer valid in extremely shallow ones. Finally, a modification for the ITTC57 correlations line considering shallow water effects is proposed, which is willing to improve the prediction of the frictional resistance of those ships with a large area of flat bottom and sail in shallow water.


List of symbols
The velocity of the incoming flow V _B The average input velocity of the flow underneath ship's bottom y + y plus, a non-dimensional wall distance for a wallbounded flow κ Von Kármán constant µ The molecular dynamic viscosity µ t The turbulent viscosity ν Kinematic viscosity ( m 2 /s) τ w Wall shear stress

Introduction
An improved understanding of the characteristics of friction/ correlation lines in various conditions can contribute a better model-ship extrapolation, which is in turn beneficial for ship design and optimization. Such prediction is based on the understanding of the boundary layer theory and friction lines derived from a flat plate. In the past century, researchers have provided several friction lines, e.g. Schoenherr [1], Grigson [2] and Katsui [3] proposed friction lines for flat plates and 1957 ITTC [4] correlation line is frequently used to predict the frictional resistance of ships. These methods were based on the results of experiments and/or numerical simulations with an unrestricted incoming flow. In these cases, the pressure gradient along the flat plate is assumed to be zero. However, practically, the incoming flow is usually restricted. For inland ships which sail in rivers and canals, ship resistance is influenced by waterway limitations. The 1987 ITTC [5] indicated that shallow water effects will be noticeable when the ratio of water depth to ship draft is less than 4.0. Various researchers [6][7][8] provided ways to predict the resistance in shallow water but without specifically and physically discussing the changes of friction, which is the major part of the resistance of most inland vessels. As the effects from the lateral direction (i.e. the effects of the sides of the waterway) are much smaller than the vertical direction (i.e. the waterway bottom) [9], this study will focus on the effects on the flat plate friction in pure shallow water only. Additionally, the inflow conditions into the slit between ship's and fairway's bottom are subject to formand flow pattern effects (towards more 2D patterns in very shallow waters), and these effects are subject to ongoing investigations.
For waterways with unlimited width, shallow water mainly affects the bottom area of the ship, and the effects on other wetted surfaces can be ignored. Since the ITTC57 correlation line cannot resolve shallow water effects, especially for the extremely shallow water [10], this correlation line should be corrected for shallow water effects.
Generally, most inland ships have a long parallel midbody and a large area of the flat bottom. The characteristics of the flow close to the bottom are comparable with the flow passing above a 2D flat plate, as shown in Fig. 1 (this is further illustrated in Sect. 2.1). We assume the velocity of far-field incoming water is V (ship-based coordinate system), and the water underneath the ship is accelerated by ΔV (due to the displacement of the ship and/or the limitation of the waterway). This acceleration even happens in deep water but will be more obvious in shallow water. Therefore, in the simplification, the speed of incoming water for the flat plate should be V + ΔV. For illustrating convenience, we keep using the symbol V to represent the velocity for flat plate and revert to the V + ΔV when calculating ship's friction; similarly, we also keep using the symbol L for the length of the plate.
A 2D flat plate has been seen as a reasonable simplification for investigating the physics of friction [2,3]. In the research of Eça and Hoekstra [11], systematic calculations of frictional resistance using CFD have been done on a flat plate with a plenty number of turbulence models. In their research, the parallel boundary was deliberately set far enough away from the flat plate to avoid shallow water effects. They compared the results when the distances between the flat plate and the parallel boundary were 0.25L and 0.5L, where L is the length of the plate, and concluded that the differences were too small to be considered. However, a distance less than 0.25L is quite common in practice, e.g. for some ports and inland waterways, the under-keel clearance may be even smaller than 0.01L [12].
In this study, extremely small under-keel clearances (up to 0.01L) are included. Shallow water effects on the physical details of boundary thickness and pressure gradients are investigated. CFD techniques are applied and it is found that both the pressure gradients along the plate and the friction were affected by the flow limitation. Finally, a regression analysis is made based on the CFD results to propose a numerical friction line for the flat plate in shallow water. Practically, to improve the prediction of ship's friction, after applying ITTC57 correlation line for all wetted surface, this proposed line can be used to correct the shallow water effects on a ship's flat bottom. The combination of these two lines is considered as a suitable method for predicting the frictional resistance in the pure shallow water.
This paper consists of five sections. Section 2 introduces the computational setup and the selection of boundary conditions. Section 3 provides the verification and validation by analyzing the numerical uncertainty and modeling errors. Section 4 presents the results of CFD simulations and the fitting of the numerical friction line; an example of how to apply the method is given. Conclusions are drawn in Sect. 5.

Approach
The feasibility of the 2D flat plate simplification is first discussed in this section. The shallow water conditions were obtained by adjusting the distance between the flat plate and a parallel wall. In this study, calculations were performed in a RANS solver: ANSYS (TM) Fluent (version 16.2). The Reynolds number in the simulations was varied from 10 5 to 10 9 .

Simplification
In practice, reasonable simplification can help to understand the nature of a phenomenon but with much smaller resources. In this study, we assume the flow passing above a 2D flat plate can represent the flow close to the flat bottom of inland vessels.
Physically, due to the presence of the bow and the stern, part of the water around a ship flows laterally leading to the dissatisfaction of the continuity of fluid on the longitudinal center plane. However, if we change the control volume from 1 to 2, as shown in Fig. 2, the continuity equation can be satisfied.
To further illustrate the assumption, an example is done for an inland vessel. As depicted in Fig. 3, part of the water flows away from the longitudinal center plane at the bow. These effects shall be considered by appropriate corrections imposed to the inflow boundary conditions for the 2D investigated flat bottom plate. Through this, when the water goes into the space underneath the bottom of the ship, its direction is no longer changed and the characteristics of the flow are comparable to that above a 2D flat plate.
Therefore, for the control volume 2, it is rational to apply the simplification mentioned at the beginning of this section.

Computational model
The computational domain is shown in Fig. 4. For conventional viewing angle, the domain is rotated by 180° compared with that in Fig. 1.
The flat plate used in this study is two meters long (L), and the domain stretches L in front of the plate and 1.5L behind it. A geometric progression with a factor of 0.5 is implemented for the distances (D) between the parallel wall and the plate, as depicted in Table 1.
Case 1 simulates the deep water condition for comparison. Case 7 is considered to be extremely shallow water. For example, D/L = 0.01 means a ship with 100 meters   The flow is the fresh water with a density of 999.04 kg/ m 3 and a kinematic viscosity of1.13902 × 10 −6 m/s 2 . This study used ten Reynolds numbers in the range from 10 5 to 10 9 , as shown in Table 2.
Four turbulence models were chosen to simulate the characteristics of the flow: Three fully-turbulent models (Spalart-Allmaras (SA), BSL k-ω and SST k-ω) and one transition model (k-kl-ω).

Mesh generation
The number of grid points in the x-direction is identical for most cases (mesh was refined for the shallowest cases and part of the high Re cases), and the distribution of grid points close to the flat plate is similar for all cases. The number of grid points in the y-direction decreases with the decrease of the distance (D). Additionally, the mesh close to the flat plate and that close to the parallel wall was refined from case 2 to case 7, because shear stress was observed on these boundaries.
Denser grids were generated near the flat plate as well as the area close to plate's leading and trailing edge, as shown in Fig. 5 (The number of grid cells in the picture has been reduced for illustration). Grids were also refined close to the moving wall for all cases except for the cases with D/L = 1.00 because the shear stress on the moving wall when D/L = 1.00 is too small to be considered. The "BiGeometric" bunching law [13] was applied to the nodes distribution, indicating that the space expansion between the nodes in each direction is linear.

Inlet boundary
An incompressible, undisturbed flow enters the domain from the inlet with a velocity of U ∞ . The inlet boundary applies Dirichlet conditions for both velocity and turbulence quantities: In Eq. 1, u x and u y are the velocities of upstream flow in x and y direction, respectively. I is the turbulence intensity, u ′ the root-mean-square of the turbulent velocity fluctuations, U ave the Reynolds averaged velocity, t the turbulent viscosity, and the molecular dynamic viscosity. C 1 and C 2 are both constants, which will be chosen in this section.
The static pressure at the inlet boundary is set to zero. The total pressure is not a constant value but will rise to whatever value is necessary to provide the prescribed velocity distribution [14].
In the research of Walters [15], where the three-equation transition model (k-kl-ω) was proposed, three different inlet boundary conditions were applied. Similar sets were chosen in this study (see Table 3).   The friction on the flat plate was recorded and shown in percent compared to a randomly chosen set (Set5: I = 3%, t ∕ = 12), which is listed in Table 4. For the SA model in Fluent, the alteration of turbulence intensity is not available when different is set. Thus, the results of friction were compared to t ∕ = 12 without setting a specific intensity. Comparisons are made for two Reynolds numbers: one is relatively low (Re = 2.51 × 10 6 ) for which the transition from laminar flow to turbulent flow is obvious; another one is high (Re = 1.58 × 10 9 ) for which the turbulent flow is dominant.
Based on Table 4, we found • For the high Reynolds number, the friction is insensitive to different I and t ∕ ; even the maximum difference is less than 0.1%; • For the relatively low Reynolds number (Re = 2.51 × 10 6 ), the calculations with the BSL k-ω and the SST k-ω model also show insensitivity to inlet boundary conditions (maximum 1.2%). For those with the SA model, the deviation is larger but still less than 5%. However, for the k-kl-ω model, the results show big discrepancies especially when t ∕ = 100, because different inlet boundary conditions significantly affect the position of "transition point" from laminar flow to turbulent flow, and the k-kl-ω model can more successively catch these changes; • For the calculations with fully turbulence models, higher turbulent viscosity ratio ( t ∕ ) leads to larger friction.
Consequently, for the BSL k-ω and SST k-ω models, the turbulence intensity and turbulent viscosity ratio at the inlet boundary have minor impacts on results. For the SA model, the influence of boundary conditions is noticeable, and for the k-kl-ω model, the results significantly depend on the boundary conditions. Which model is the most suitable one is decided after verification and validation (see Sect. 3).

Other boundary conditions
The "pressure outlet" was set at the outlet boundary. All derivatives of the flow quantities at x direction were set to zero. The position of outlet boundary should be far enough from the flat plate to avoid influencing the gradient of fluid variables in the domain. Figure 6 shows the pressure gradient at x direction from the trailing edge to the outlet. Three distances (0.5L, 1.0L and 1.5L, L = 2 m) are tested.
For the outlet boundary of 0.5L, 1.0L and 1.5L behind the flat plate, the pressure gradients to the x-direction near the outlet boundary are 1.149, 0.186 and 0.0074, respectively. A value close to zero means the outlet boundary is far enough and will not affect the flow. Therefore, the outlet boundary which is 1.5L behind the plate is enough for this study.
In addition, for the "backflow" (if any) at outlet boundary, Dirichlet conditions were set to the turbulence intensity and turbulence-viscosity ratio with the same values as at the inlet boundary.
The parallel wall was set as a non-slip "moving wall" with the same speed as the free stream. All derivatives of the flow quantities in the y-direction were set to zero, and the speed relative to the parallel wall in the x-direction was set to zero. The flat plate was set as a still, non-slip wall. Dirichlet conditions were set to the velocities, i.e. velocities at x and y directions were zero.
Symmetry conditions were set in front and behind the flat plate. The velocity and the derivatives of all flow quantities at the y-direction were set to zero.

Verification and validation
Numerical simulations shall be verified and validated before application. This is done in Sect. 3.1 and 3.2 respectively.

Verification
The verification gives the numerical error and uncertainty of a simulation. Typically, a numerical error is the discrepancy between a numerical result and the exact solution. According to Roache [16], the numerical error contains three components: round-off error, iterative error, and discretization error.
Round-off error results from the finite precision of the computers, but the double precision format can usually keep this error negligible [17]. The iterative error is a consequence of the non-linearity of the mathematical equations. Using double-precision scheme and a sufficient number of iterations can normally reduce this error to the level of the round-off error. In this study, the convergence criteria of all residuals are set to 10 − 7 . This does not mean the iterative error is 10 −7 , and actually, this usually keeps the iterative error at the level of 10 − 9 . The discretization error is generated when transforming the partial differential equations into algebraic equations and mostly dominates the numerical error. Therefore, in this article, only the discretization error is considered.
A grid refinement study is commonly used to estimate the discretization error [16]. At least four grids are recommended [18] to justify whether the results are in the "asymptotic region". In this verification, four geometrically similar grids (G1-G4) with a refinement ratio of 1.6 were generated for the Case1 (D/L = 1.0, the deep case). Both the number of nodes (Table 5) and the distance between the nodes (e.g. the y + , see Table 6) had a ratio as close as possible to 1.6. No wall functions were used, and even the coarsest grid had a y + less than 1.
In this grid refinement study, a turbulence intensity of 3% and a turbulent viscosity ratio of 12 were used. Two Reynolds numbers (Re = 2.51 × 10 6 and Re = 1.58 × 10 9 ) were chosen for the verification with three turbulence models, as shown in Fig. 7.
The scale of different grids can be represented by a symbol h i , which is defined as follows: where h i is the symbol representing the scale of Grid i; N xi and N yi are the numbers of Grid i at x and y directions, respectively. The value of h 1 is set to 1 by default. The value of h i is useful for presenting the results of grid refinement study (i.e. the x-axis of Fig. 7).
Although this study applied the three-equation k-kl-ω model for Re = 2.51 × 10 6 , the results showed oscillatory convergence ( Table 7).
The symbol R is defined as where C fi is the frictional resistance coefficient of Grid i (i ≥ 3). For R < 0 and |R| < 1 , C f has an oscillatory convergence [18]. Following the procedure of numerical uncertainty analysis proposed by Eça and Hoekstra [18], the uncertainty of C f for the finest grid (G1) and the corresponding observed order of accuracy p is shown in Table 8. Since only the monotonously convergence was considered, the uncertainty analysis of the calculations with the k-kl-ω model was excluded in this table.
Based on Fig. 5 and Table 8, it can be observed  • In the grid refinement study, the value of C f depends on Reynolds number and the selected turbulence models. C f increases with the refinement of the mesh. Calculations with SA model generate the highest C f and SST k-ω leads to the lowest C f for both Reynolds numbers; • The chosen turbulence models affect the C f more than the density of grids. Those changes of C f become smaller with the mesh refinement; • Calculations at Re = 1.58 × 10 9 show smaller uncertainty than that at Re = 2.51 × 10 6 , but both of their uncertainties are smaller than 1%.
Consequently, the finest grid (Grid 1), which has been verified in this section, was chosen for the further systematic simulations.

Validation
After estimating the numerical uncertainty, the deviations of simulations from experimental data (modeling error) are estimated.
In this section, the results of global friction (C f ) on the flat plate are validated using the friction lines proposed by Prandtl--Schlichting [19] and Katsui et al. [3]. Different inlet boundary conditions listed in Table 3 are compared and the best suitable one is chosen in Table 9. Finally, the local skin friction and the mean flow velocity profile, are validated with the experimental data of Nagib et al. [20].

The friction proposed by Katsui et al. [3]:
For both Re = 2.51 × 10 6 and Re = 1.58 × 10 9 , the difference of C f in unrestricted condition (D/L = 1.0) is compared with these two friction lines, which is shown in Table 9.
From Table 9, it can be derived that where ∶ a = 0.042612 lg Re + 0.56725.  For full-scale ships, turbulent boundary layer starts immediately after the flow reaches the bow and dominates the whole boundary layer. Therefore, a similar proportion of turbulent boundary layer is expected for the geometrically similar model ships. To achieve this, turbulence stimulations (such as sand, tripwire, and pins [21]) are applied. Consequently, using fully turbulence models is more comparable with the physics of ship model tests than using a transition model. Additionally, according to Table 9, the BSL k-ω model with I = 6% and t ∕ = 100 has a good agreement with the fully turbulent method at both low and high Reynolds numbers. This set of turbulence model and boundary conditions are, therefore, chosen for further systematic simulations.
Nagib et al. [20] evaluated many previously proposed empirical local friction formulas based on the more recent experimental datasets (Nagib et al. [22] and Österlund [23]). Through proper modifications, they concluded that those formulas can accurately describe the local friction. The mean velocity profiles in the turbulent boundary layers were also shown in a graph against a large range of momentum thickness Reynolds number (Re θ ), which provides validating data for this research.
Based on the fitting of Nagib, the local Reynolds number (Re x ) can be converted to the momentum thickness Reynolds number (Re θ ) using The CFD calculations of local friction coefficients (C fx ) are validated in Fig. 8. Compared the fitting results of Nagib, the CFD calculations have up to 6% higher values for C fx . If the classical local friction line, i.e. Prandtl-Schiliching, is used, the CFD results generally underestimate the C fx by about 4%. The data of Nagib is based on more recent experiments, which are assumed to be more reliable. Nevertheless, an error less than 6% is practically acceptable for CFD calculations unless an extremely rigorous requirement is posed.
Additionally, the mean velocity profiles are validated in Fig. 9. The CFD results agree well with the experiments, especially for the viscous sublayer (y + < 8, roughly) and loglaw region (50 < y + < 500, roughly). In other regions, CFD calculations tend to underestimate the velocity, which means a thinner boundary layer and thus a higher velocity gradient, which can also explain the discrepancies of the C fx in Fig. 8. Based on the validation, CFD techniques applied in this study have acceptable errors and can be trusted for further calculations.

Results and analysis
The systematic simulations involve ten Reynolds numbers and seven shallow water conditions, as shown in Table 10. The number of grid cells in the x-direction (Nx) is unchanged from Re = 3.9811 × 10 5 to Re = 2.5119 × 10 8 , and to stabilize the calculation, Nx was increased to 827 (finer mesh) at higher Reynolds numbers. The number of grid cells in the y-direction (Ny) was varied according to the distance between the flat plate and the parallel wall.

The physical effects of shallow water on friction
Physically, the flow will be accelerated due to the displacement of the boundary layer, leading to a thinner boundary layer especially in extremely shallow water conditions. To clarify the changes of boundary layers, a square area close to the flat plate was chosen, as shown in Fig. 10. The boundary layer thicknesses at the trailing edge of the plate at two Reynolds numbers are shown in Fig. 11. The space in the y-direction is amplified 5 times to clearly illustrate the physical details in the boundary layer.
From Fig. 11, some remarks can be derived: Fig. 9 Mean velocity profiles and log-law diagnostic function against y + (U + is the dimensionless velocity defined by U + = u∕ √ w ∕ , where u is the flow velocity parallel to the plate and τ w is the wall shear stress; in the logarithmic law, κ = 0.394 according to Nagib) • For both Re = 2.51 × 10 6 and Re = 1.58 × 10 9 , the maximum flow speed increases and boundary thickness decreases with the decrease of D/L; • When D/L = 0.01, where the distance between the flat plate and the parallel wall has the same order of magnitude as the boundary layer thickness, the boundary layer thickness decreases more significantly than that for D/L ≥ 0.02.
As a consequence, the velocity gradient normal to the flat plate increases leading to a rise of the local friction (C fx ). Figure 12 shows the results of C fx at Re = 2.51 × 10 6 and Re = 1.58 × 10 9 in different shallow conditions. Some remarks can be drawn from Fig. 12: • For both Reynolds numbers, limited water depth has an increasing influence on C fx from the leading edge to the trailing edge; • The C fx increases rapidly with the decrease of D/L; • Influence of shallow water on local friction is more significant at the low Reynolds number. For example, at the trailing edge, when D/L = 0.01, the increase of C fx can be more than 50% at Re = 2.51 × 10 6 but is only 20% at Re = 1.58 × 10 9 . This discrepancy shows important scale effects when extrapolating resistance in shallow water from model scale to full scale.
Consequently, the total frictional resistance, which is the sum of the local frictions, is also increased on shallow water conditions and depends on D/L as well.

The validity of zero-pressure gradient (ZPG) assumption
An assumption in the previous research [1-4, 11, 24] about the flow passing a flat plate is that the pressure gradient along the plate is zero. This makes it easier to derive theoretical solutions for friction lines and usually has a good agreement with experiments.
For an incompressible, two-dimensional laminar flow, the control equations within the boundary layer are Blasius [25] (translated by NACA [26]) set the term − −1 × p∕ x to zero in Eq. 2 and derived the Blasius solution.
For a two-dimensional, fully-developed turbulent flow, based on the ZPG assumption and a large amount of experiments results, Schlichting [19] fitted a formula for C f in fully turbulent flow: However, the validity of the ZPG in shallow water should be tested. In this study, two Reynolds numbers were again chosen: Re = 2.51 × 10 6 and Re = 1.58 × 10 9 . Sample points were picked on three straight lines which close to the boundary layer, with the same length of the flat plate and offset the plate by d = 5 mm, 10 mm, and 15 mm, as shown in Fig. 13. Values of pressure gradient and velocity gradient were recorded at these points.
To test the ZPG assumption, the problem becomes to compare the order of magnitude of u u∕ x and − −1 × p∕ x . The ZPG assumption is valid only when they have the same order. Figure 14 presents the values of u u∕ x and − −1 × p∕ x at different conditions.
When D/L is 1.00, − −1 × p∕ xare close to zero, which can be ignored compared to u u∕ xfor both Reynolds numbers. However, when D/L drops to 0.01, the magnitude of − −1 × p∕ xincreases significantly and even reaches the same order as u u∕ x. Therefore, when D/L ratio is small enough, the item − −1 × p∕ xcan no longer be ignored, both at low and high Reynolds numbers. In other words, the ZPG assumption is invalid at extremely shallow water conditions. Thus, those methods based on the zero-pressure gradient assumption are insufficient.

Shallow water effects on the frictional resistance
Results of the total frictional resistance coefficients of the flat plate and their increase compared with deep water condition (in percent) are presented in Fig. 15. Based on Fig. 15  As a result, the influence of shallow water conditions on the C f again has scale effects, and the effect is larger at low Re because the boundary layer is thinner will is less likely to be influenced at high Re.
In practice, the scale effects at shallow water conditions may cause large errors in engineering. For instance, for an inland ship with 100 meters long and sails at 18 km/h (relative to the flow under the ship bottom), its Reynolds number is about 4.4 × 10 8 , and it sails in shallow water with the D/L = 0.01. When using a 1/25 scaled model, which has a Reynolds number about 3.5 × 10 6 , to extrapolate the total resistance, C f of the bottom part at full scale will increase by about 11% compared to deep water, but it will increase by about 28% at model scale according to this study. This discrepancy will lead to an unreliable resistance extrapolation if the IITC57 correlation line is applied. Consequently, a new friction line considering shallow water effect is needed to correct the shallow water effects.

The fitting of a numerical friction line
This section proposes a numerical friction line for correcting the shallow water effects on ship's bottom based on the described CFD calculations.
A regression analysis using the method of least squares was applied. During the procedure, choosing a suitable function model is essential to the quality of the fitting. Inspired by the ITTC57 [4] formula, we chose the function model as follows: The a, b and c are constants defined in the fitting. In this study, we followed two steps in the regression analysis: first, fit a numerical friction line in deep water condition using the results of D/L = 1.00; Second, include the parameter D/L based on the numerical results of other values of D/L.
The result for the first fitting with the R-squared value (a statistical measure of how close the data are to the fitted regression line) of 0.9996 is shown as follows: The C f in shallow water conditions depends on Reynolds number and D/L, and additionally, when D/L is close to 1, which approaches the deep water, the numerical friction line should approach the Eq. 3.
Based on the results of the simulations, a numerical friction line in shallow water conditions is proposed using regression analysis: Re is Reynolds number; D is the distance between flat plate and parallel wall; L is the length of the flat plate.
Errors between the simulations and the new proposed friction line are shown in Table 11.
The errors are mostly less than 1%, except when D/L = 0.01(still less than 3%). This error table can be referred when using this new friction line.

Application and a case study
The previous section proposed a modification of the ITTC57 correlation line. However, since a ship also has non-horizontal wetted surfaces, this method cannot be applied directly. In this paper, it is suggested to modify the frictional resistance on the flat bottom area with the proposed friction line and apply the ITTC57 correlation line for the non-horizontal wetted surface. Therefore, we propose the following steps to use this method in the prediction of the frictional resistance of an actual ship: As a case study, the above steps are applied for a 1/30 scaled model ship, whose prototype is an 86-m long inland ship. The free surface is not considered in CFD calculations to eliminate wave effects on friction. Details of the ship parameters can be found in [10]. L here is the length of the flat bottom.
To achieve the average velocity (V _B ) at the leading edge of the wide flat bottom (as indicated in Fig. 3), a systematic CFD calculation has been done for various incoming velocities and water depths. The results of V _B are shown in Fig. 16, and an empirical formula is regressed for this specific case with an uncertainty of 2.5%: For h/T ≤ 4.0, Following the steps shown at the beginning of this section, the frictional resistance of this ship model is predicted at 0.8 and 1.0 m/s in various water depths, as shown in Table 12.
For all the three different water depths, results of C f using the proposed method agrees better with the values derived from CFD than ITTC57. With the errors less than 3%, the improved correlation line is practically considered applicable to successfully predict the friction including shallow water effects on ship's flat bottom. (4)

Conclusions
This study has proposed a numerical friction line for correcting shallow water effects on ship's bottom on shallow water conditions using CFD calculations. A moving plate approaching a flat plate was applied to build the shallow water conditions. Three fully turbulence models (Spalart-Allmaras, BSL k-ω, and SST k-ω) and one transition model (k-kl-ω) were applied to verify and validate the numerical calculations. Practically, the BSL k-ω was the most suitable model, and the turbulence intensity of 6% with the turbulent viscosity ratio of 100 was the best inlet boundary condition in this study.
Based on the calculations in shallow water conditions, we concluded • For low Reynolds numbers, D/L significantly influences the friction on the flat plate. When lg(Re) = 5.6, the increase of C f can reach almost 50% compared to the deep water condition; • The influence of shallow water on friction has scale effects, which diminishes with the increase of Reynolds number. For a relatively high Reynolds number (e.g. lg(Re) = 9.2), the increase of C f is only 10% compared to the deep water condition; • The zero-pressure gradient assumption is invalid in extremely shallow water.
Reynolds number is no longer the only factor influencing friction at shallow water conditions. A numerical friction line considering the shallow water effects and the scale effects is proposed. After applying this proposed line to modify the frictional resistance on ship's bottom can For resistance extrapolation in shallow water from ship models to full-scale ships, the proposed friction line can also correct the scale effect on the friction on ship's bottom and thus, increase the reliability of the model-ship extrapolation of resistance.
In the future, more factors will be considered to improve the prediction of ship's frictional resistance in shallow water. Ship waves, trim and sinkage are all possible factors influencing ship friction. The physical explanation proposed in this paper, which shows the changes in the boundary layer, is applicable if those factors are considered in the future research. The modification to include the parameters of the waterway also shows the way to update the friction or model-ship correlation line in shallow water.