Next Article in Journal
Fractionation Analysis of Iron in Coastal Rivers to Yantai Sishili Bay with a Bismuth Microrods-Based Electrochemical Sensor
Next Article in Special Issue
Effect of Database Generation on Damage Consequences’ Assessment Based on Random Forests
Previous Article in Journal
1D–3D Coupling Algorithm of Gas Flow for the Valve System in a Compression Ignition Engine
Previous Article in Special Issue
Review on Ship Manoeuvrability Criteria and Standards
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mesh Properties for RANS Simulations of Airfoil-Shaped Profiles: A Case Study of Rudder Hydrodynamics

1
School of Naval Architecture, Ocean and Energy Power Engineering, Wuhan University of Technology, Wuhan 430063, China
2
Intelligent Transportation Systems Research Center, Wuhan University of Technology, Wuhan 430063, China
3
National Engineering Research Center for Water Transport Safety, Wuhan 430063, China
4
Maritime and Transport Technology, Delft University of Technology, Mekelweg 2, 2628 CD Delft, The Netherlands
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(10), 1062; https://doi.org/10.3390/jmse9101062
Submission received: 31 July 2021 / Revised: 3 September 2021 / Accepted: 5 September 2021 / Published: 28 September 2021
(This article belongs to the Special Issue Ship Dynamics and Hydrodynamics)

Abstract

:
A good mesh is a prerequisite for achieving reliable results from Computational Fluid Dynamics (CFD) calculations. Mesh properties include mesh types, computational domain sizes, and node distributions. However, in literature, we found no clear consensus about what these properties should be. In this article, we performed a case study on ship rudders to determine what the suitable mesh properties are for airfoil-shaped profiles. A classic NACA 0012 profile is chosen as an example, and commercial packages ANSYS ICEM are applied for meshing with an ANSYS Fluent solver. With a strategy in consideration of relationships among different mesh properties, a comprehensive parametric investigation is conducted to study the impacts of these properties on the accuracy of rudder hydrodynamic coefficients obtained by CFD methods. The step-by-step study outputs recommended Reynolds numbers, domain sizes, and near- and far-field node distributions for mesh types with distinct topology structures, i.e., C-mesh, O-mesh, H-mesh, and Hybrid-mesh. Specifically, the study shows that a critical Reynolds number is needed for the perspective of efficiency, while a domain extending 60 times of the chord length enables the boundary effects to be negligible. As for node distributions, the near-field nodes should be treated carefully, compared with those in the far-field. After that, corresponding mesh properties for different calculation objectives are illustrated in detail based on the characteristics of mesh types mentioned above. With the proposed strategy for mesh refinements, impacts of different mesh properties on rudder hydrodynamics are clarified and recommended settings are applicable for other airfoil-shaped profiles such as wind turbines and marine propellers.

1. Introduction

Stimulated by developments in computer power and mechanics theories, Computational Fluid Dynamic (CFD) methods were widely applied in ocean engineering to study the hydrodynamic performance of marine structures. Generally, computation domains containing target structures need to be discretized into multiple mesh elements to solve differential equations numerically. In ship maneuvering simulations, the presence of CFD-based tools enables researchers to obtain more precise hydrodynamic coefficients compared with that of traditional empirical methods, which is beneficial to establish individualized and accurate mathematical maneuverability models for different objectives. As the most widely used steering devices for ships, rudders play an important role in the maneuvering of ships [1], and the accuracy of lift coefficients has a great impact on motion predictions. Many CFD simulations were conducted to study rudder hydrodynamics (Liu et al. [2], Badoe et al. [3], Van Nguyen and Ikeda [4]), but few concentrate on mesh generations for rudders. In this paper, the study on mesh properties for RANS simulations of rudder hydrodynamics is performed to present detailed insight into meshes around rudders.
A good mesh is fundamental to the convergence and the accuracy of solutions, while a poor one leads to nonconverging results, inaccurate solutions, and an increase in computation time. Due to the complexities of geometries in structural surfaces and the lack of automatic generation methods, it costs plenty of time to generate meshes in calculation processes. On the one hand, mesh quality is influenced by various factors, but on the other hand, it is unrealistic to adjust all of them to the best states in practice, which requires more computation power and computation time.
A good quality mesh should have a sufficiently large number of cells, implying that a further increase in the number of cells will not lead to significant changes in the solutions. The mesh needs to be fine enough in the regions where fluid characteristics change rapidly. Therefore, mesh properties need to be coordinated to make a trade-off between the accuracy and the time scale. In literature, different mesh properties, i.e., mesh types, domain sizes, and node distributions, are applied in RANS simulations of airfoil-shaped profiles. Lutton [5] indicated that the O-mesh is superior to the C-mesh in determining the pressure coefficients in the vicinity of leading and trailing edges, while the C-mesh shows a better resolution in the wake. Basha [6] constructed a structured O-mesh, a structured C- mesh, and a hybrid mesh on a NACA 0012 airfoil, while the hybrid mesh shows better results in drag coefficient predictions due to finer mesh resolutions in the boundary layer domain. Stuck et al. [7] selected a hybrid mesh with quadrilateral unstructured cells rather than triangular ones to achieve smoother transitions from the boundary mesh to the unstructured outer mesh for RANS simulations on a rudder profile. Parametric grid independence studies are conducted by varying chord-wise and layer-wise grids, respectively, and optimum values are chosen based on variations of integral force coefficients. Turnock et al. [8] carried out decoupled studies on the boundary location and mesh node distributions in different directions for a NACA 0012 section, while related parameters are determined by convergence study.
As no clear consensus about mesh properties was found in current studies, systematic parametric investigations on mesh properties are performed in this paper. As naval architects, our primary interest in CFD lies in applications related to the flow around ships and resultant hydrodynamic forces. In this particular case, we focus on the flow around the rudder and the resulting rudder performance. Kim et al. [9] stated that the NACA series are the most widely used economic rudder profiles. Thus, a classic NACA 0012 airfoil-shape rudder profile is utilized to study mesh properties for its simplicity in the shape, which makes mesh properties easier to change and enhances the applicability of obtained results for other complicated structures.
This paper analyzes the impacts of mesh properties on RANS simulations and proposes suitable meshing strategies for rudder hydrodynamics. Following this introduction, Section 2 further discusses the reason why a 2D rudder profile is selected as the research object in this paper. Section 3 explains definitions of mesh properties in this study. Section 4 illustrates applied research methods in detail, and Section 5 discusses the impacts of different mesh properties based on simulated cases. Finally, a recommended meshing strategy is presented to provide guidance for hydrodynamic analysis in Section 7.

2. The Selection of the Research Object

Bare hulls, propellers, rudders, or full-appendages ships were all studied in CFD methods in literature, though not all of them are suitable for mesh properties studies. There are great curvature changes on the surface of hulls and propellers, and that is why mixture mesh topologies and millions of cells are desired to capture their flow features. The complex nature makes it difficult to perform parametric investigations on a single mesh property economically. What’s more, the geometries of hulls or propellers in the marine industry vary from each other significantly, and rules derived from limited cases may not be applicable for those with different configurations.
Although special rudders drew a growing interest, conventional spade rudders, or those modified but still geometrically similar versions like flap rudders, are still mainstream maneuvering devices, however. Besides, the meshing strategy aiming at a NACA profile can be easily applied to its geometrically similar counterparts. Force coefficients of rudders during ship maneuverings can be obtained based on 2D profiles in open water conditions, which also makes the investigation of 2D profiles beneficial to maneuverability studies. For all these reasons, in this paper, a 2D profile, in other words, a rudder with an infinite aspect ratio, is determined as the research object. According to Fujii [10], Fujii and Tsuda [11,12], the normal force coefficient of actual rudders whose aspect ratios are limited can be expressed as:
C N = 6.13 sin α R Λ G Λ G + 2.25
where C N is the rudder normal force coefficient, and 6.13 is an approximate constant to estimate the partial derivative of C N to sin α R of 2D sections. Liu et al. [2] calculated C N sin α R for different rudder profiles by performing 2D CFD simulations, which can modify the original transformation relation in Equation (1) and obtain force coefficients for 3D rudders in open water conditions. Since rudders generally operate in the propeller slipstream, rudder hydrodynamics affected by the propeller can be estimated based on series model tests by Nienhuis [13], and the experimental data are shown in Figure 1. Compared with that of open water conditions, the stall angle for the rudder increased to about 35 with propeller impacts, while C N sin α R was not significantly changed. Hence, based on 2D simulation results, corrected formulations for a rudder with an aspect ratio Λ G considering propeller impacts can be expressed as:
C L Λ G = k p C L 2 D sin α R sin α R + C L 0 2 D Λ G Λ G + 2.25 C D Λ G = k p C D 2 D sin α R sin α R + C D 0 2 D Λ G Λ G + 2.25 , α R 35
where C L 2 D sin α R , C L 0 2 D , C D 2 D sin α R , C D 0 2 D can be obtained from 2D CFD simulations, and k p can be estimated experimentally. The availability of Equation (2) in ship maneuvering studies was validated by Liu et al. [2,14] against two typical inland vessels. Apart from the transformation relation, another factor that makes the 2D study effective is that span-wise mesh generations can be omitted and cells in the profile section can be focused on. Therefore, 2D rudder profiles are set as the research object in following studies.

3. Mesh Properties

In literature, mesh types and computation domain sizes are generally determined before studying grid independence on certain mesh parameters, and some settings are shown in Table 1 and Table 2. Mesh properties in grid independence study are denoted as near-field and far-field parameters according to locations of distributions. All of these mesh properties are studied in this paper. The detailed descriptions of mesh properties above are in Section 3.1, Section 3.2 and Section 3.3.

3.1. Mesh Types

Based on the connectivity of elements, meshes are identified as structured, unstructured, and hybrid. As regular connectivity can be expressed as a two/three-dimensional array, the structured model is highly efficient in the usage of computational resources. However, regular connectivity restricts the element type to quadrilateral in 2D and hexahedral in 3D, which inherently limits applications of the structured mesh for complex geometries. For simple geometries, a structured mesh may have better convergence and higher resolution than an unstructured mesh.
Instead of a single-block structured mesh, a block-structured or multiple-structured mesh is more widely applied. Three commonly used topologies, C-mesh, H-mesh, and O-mesh, are shown in Figure 2. It is superior in computational efficiency to the unstructured mesh and more flexible in handling the complex geometries than the single-block structured mesh. The choice of the mesh topology has a considerable influence on mesh quality [22], which depends on the domain geometry, the structure of the solution, and the topology in the adjacent block. According to Liseikin [22], in H-mesh, the computational domain is square and two singularities exist in the interior boundary. O-mesh and C-mesh are both in solid squares, while the interior boundary for the former one is smooth and the latter has a singularity the same type as the H-mesh. In this particular case, some modifications are made based on basic topologies above to accurately capture the characteristics of the flow near the rudder, as shown in Figure 2. For brevity, modified types are still marked as C-mesh, H-mesh, and O-mesh, respectively.
The irregular connection of unstructured meshes allows for more freedom in element choices, typically triangular in 2D and tetrahedral in 3D. Compared to that of the structured mesh, the unstructured mesh is more suitable for complex configurations, but it can be highly inefficient. Aftosmis et al. [23] reported that unstructured triangular meshes are 50 times more expensive in both memory and time than that of structured quadrilateral meshes with nearly the same accuracy. Unstructured meshes, however, have advantages over structured meshes in mesh refinement and generation time. In practice, hybrid meshes are applied instead of pure unstructured meshes. Hybrid meshes commonly consist of a structured portion near wall surfaces to capture the boundary layer and an unstructured portion that fills the domain. Thus, hybrid meshes inherit the advantages of both structured and unstructured meshes, such as good orthogonality to wall surfaces, suitability for mesh refinement, flexibility for complex geometries, and fast generation. Figure 3 demonstrates a hybrid triangle mesh, which has quadrilateral elements for the boundary layer and triangular elements for the remainder of the domain.

3.2. Domain Sizes

The computational domain is the geometrical region that bounds the numerical simulation [17]. To obtain highly accurate solutions, the position of boundaries has to be discussed to demonstrate that the interior flow field is not disturbed. Thus, there is the need to study the domain size, which not only minimizes the influences of the boundaries, but also prevents an excessively large domain. Typical mesh types correspond to related topology structures, which lead to different shapes of calculation domains. The geometric parameters for different mesh types are shown in Figure 4, while subscripts, i.e., i n , o u t , and s i d e , stand for distances from the rudder to the inlet, outlet, and side boundary of corresponding mesh types, and O r is the radius of the domain of O-mesh.

3.3. Node Distributions

For structured meshes, since cell distributions vary with node distributions in different directions of the computation domain, convergence studies need to be done by varying related mesh parameters. In this study, node distributions are divided into two categories, including near- and far-field distributions.
Near-field distributions are located in the boundary layer regions near the wall, which are generally refined to capture the flow around the structure. The chord-wise spacing is the distance between two nodes along with the profile, and the layer-wise growth rate is the distances between adjacent mesh points along a mesh line, which determine the boundary layer mesh along with the first mesh height.
Different mesh types share similar near-wall distributions, while it is difficult to define far-field distributions owing to differences in mesh topologies. In this paper, far-field distributions are defined as node distributions that do not contain near-field distributions.

4. Research Approach

A series of parametric studies are conducted to observe the impacts of mesh properties on the prediction of the accuracy of solutions. The classic validation profile NACA 0012 is tested through commercial CFD code ANSYS Fluent. Strategies for investigating mesh properties mentioned in Section 4.1. Numerical methods utilized to study the hydrodynamic performance of the rudder profile are described in Section 4.2, and the error calculation method is mentioned in Section 4.3.

4.1. Strategies

Different mesh properties contribute to mesh quality in varying degrees, and it is the combination of these properties that determines the accuracy of calculation results. Instead of rough global refinement of cells [24], several main impact factors on mesh fineness and the required number of cells are illustrated with simulations. This parametric mesh independence study is chosen because optimization on individual parameters specifically addresses the crucial impact factor of mesh quality, avoiding the waste of cells [7]. The global refinement is more frequently applied for unstructured meshes than structured meshes. However, due to the lack of specific classification on mesh properties, mesh properties above are not independent of each other strictly. Cases with different mesh types also differ from each other in domain sizes, node distributions, and numbers of mesh elements, etc.
As shown in Figure 5, the typical mesh type corresponds to a related topology structure, which leads to different shapes of calculation domains and connection patterns between the profile and domain boundaries. With domain size varying, node distribution in different directions changes. Further, the number of mesh elements is concerned with all mesh properties above.
The strategy in Figure 6 is applied to investigate how different mesh properties affect the calculation results of rudder hydrodynamics. Firstly, cases of C-mesh type with Reynolds numbers in a wide range are tested to obtain the critical Reynolds number ( R e c r i t i c a l ), above which hydrodynamic coefficients tend to be stable. Secondly, the domain sizes of the C-mesh are varied to obtain the proper boundary sizes that are large enough to eliminate boundary effects for the following simulations. Thirdly, cases are divided into 4 groups classified by 4 different mesh types. Since it is inappropriate to compare mesh qualities of different mesh types when other mesh properties are not determined yet, the impacts of node distribution on each mesh type should be studied independently. In this part, the node distribution is the only factor that affects mesh quality for the typical mesh type, and the effects of node distributions are studied based on R e c r i t i c a l and the domain size obtained in the previous stage. Therefore, grid independence studies are conducted on parameters in the near-field, i.e., chord-wise spacings and layer-wise growth rates, and the far-field, i.e., radial nodes, wake nodes, and element sizes. Finally, four mesh types with optimal mesh properties are compared to give recommended mesh settings for rudder hydrodynamics.

4.2. Numerical Methods

The RANS methods presented in this paper use a time-average Reynolds decomposition, which assumes that all the components of the flow velocity and pressure consist of a mean value and a bounded fluctuation to represent turbulence. The governing equations for incompressible viscous flow are as follows:
u ¯ i x i = 0
ρ u ¯ i t + ρ u ¯ j u ¯ i x j = f i p ¯ x i + x j μ u ¯ i x j ρ u i u j ¯ , ( i = 1 , 2 , 3 )
where u ¯ i is the mean velocity component in the x direction, ρ and μ are the density and viscosity coefficient of the fluid, P ¯ is the mean pressure, f i is the body force, while ρ u i u j ¯ is the Reynolds stress. In RANS simulations, turbulence models are introduced to close equations above. The k- ω SST model was designed to give highly accurate predictions of the onset and the amount of flow separation, which provides a better description of flow around the rudder. Therefore, the k- ω SST model is applied for simulations in the following parts. Numerical settings for ANSYS Fluent is shown in Table 3.

4.3. Error Calculation

Routinely, nondimensional coefficients are used to compare the rudder hydrodynamic performance with that of various design choices, while main hydrodynamic characteristics are the lift coefficient C L , the drag coefficient C D , and the moment coefficient C M ([1]), which are expressed as:
C L = L R / 0.5 ρ V R 2 A R
C D = D R / 0.5 ρ V R 2 A R
C M = M R / 0.5 ρ V R 2 A R c
where L R , D R and M R are the lift force, drag force, and moment, respectively, ρ is water density, α is the angle of attack, V R is the rudder inflow speed, A R is the rudder area, and c is the chord length of the rudder profile. The center of pressure is set as 0.25 c from the rudder’s leading edge, based on engineering practice and experimental data [25].
To evaluate the mesh quality, wind tunnel results of the NACA 0012 profile, which are tested by Ladson [25], are taken as the benchmark case. Moreover, three independent compressible CFD codes [19], i.e., CFL3D (NASA LaRC, USA), FUN3D (NASA LaRC, USA), and NTS (NTS, Russia), are cited to make a crosswise comparison with the Fluent solutions. Results of the benchmark and the CFD codes are achieved in essentially incompressible air with a Mach number of 0.15. Thus, they are comparable with the results obtained in incompressible water in this paper. Presented results are compared to that of the benchmark cases in relative differences under various angles of attack at a critical Reynolds number as follows:
C α = C α F l u e n t C α B e n c h m a r k 1
C ( AVE ) = i = 1 n C α i n , i = 1 , 2 , , n
where C α is the relative differences of the calculated lift, drag, or moment coefficients achieved by Fluent and provided by the benchmark data, respectively, and C ( AVE ) is the average relative differences for a series of α .

4.4. Uncertainty Estimation

The uncertainty estimation is important to evaluate the reliability of data in CFD calculations. According to International Towing Tank Conference [26], simulation numerical uncertainty U S N and simulation modeling uncertainty U S M need to be assessed respectively. In this paper, only the numerical errors are considered since modeling errors are difficult to quantify in practical applications [27]. Generally, numerical errors relate to grid size, time step, and other parameters. This paper performs steady simulations and focuses on mesh properties. Accordingly, only the grid uncertainty U G is calculated for uncertainty estimations. With validation uncertainty U V , which is simplified as U G in this study, obtained, U V and the comparison error | E | can be compared. If | E | < U V , the validation is achieved at U V level.
The Grid Convergence Index (GCI) method [28,29] and the correction factor (CF) based method [30] remain popular in uncertainty estimations, however, it is deficient that both methods fail to give reliable estimations when the order of grid convergence p is too high (i.e., p > 2 ) or too small (i.e., p < 0.5 ). In other words, GCI and CF methods require data in the so-called asymptotic range, which means that data on grids are fine enough to give a single dominant term in a power series expansion of the error [31]. Considering that practical flow phenomenons are generally complex, and it is inevitable that p may be not in the range of [ 0.5 , 2 ] , Eça and Hoekstra [31] proposed another GCI method based on least-square roots (LSR-GCI). Such a method allows for error estimations even if the monotonic convergence is not achieved, and it is proven to be more reliable compared with GCI and CF based methods by De Luca et al. [32]. Therefore, the LSR-GCI method [31] is applied for uncertainty estimations in this paper.
For four solutions ϕ i , ( i = 1 , 2 , 3 , 4 ) obtained from systematically refined cases and characterized by typical cell sizes h i , ( i = 1 , 2 , 3 , 4 ) , the discretization error can be estimated with one of the four following equations:
ε ϕ ϕ i ϕ 0 = α h i p
ε ϕ ϕ i ϕ 0 = α h i
ε ϕ ϕ i ϕ 0 = α h i 2
ε ϕ ϕ i ϕ 0 = α 1 h i + α 2 h i 2
where ϕ 0 is the estimate of the exact solution, α is an undetermined constant, and p is the observed order of grid convergence. Besides, the following functions are proposed to select the error estimator:
S R E ϕ 0 , α , p = i = 1 n g ϕ i ϕ 0 + α h i p 2
S R E w ϕ 0 , α , p = i = 1 n g w i ϕ i ϕ 0 + α h i p 2 , ( w i = 1 h i i = 1 n g 1 h i )
S 1 ϕ 0 , α = i = 1 n g ϕ i ϕ 0 + α h i 2
S 1 w ϕ 0 , α = i = 1 n g w i ϕ i ϕ 0 + α h i 2 , ( w i = 1 h i i = 1 n g 1 h i )
S 2 ϕ 0 , α = i = 1 n g ϕ i ϕ 0 + α h i 2 2
S 2 w ϕ 0 , α = i = 1 n g w i ϕ i ϕ 0 + α h i 2 2 , ( w i = 1 h i i = 1 n g 1 h i )
S 12 ϕ 0 , α 1 , α 2 = i = 1 n g ϕ i ϕ 0 + α 1 h i + α 2 h i 2 2
S 12 w ϕ 0 , α 1 , α 2 = i = 1 n g w i ϕ i ϕ 0 + α 1 h i + α 2 h i 2 2 , ( w i = 1 h i i = 1 n g 1 h i )
Details in selecting the error estimator are shown in Figure 7 according to Eça and Hoekstra [31]. ϕ 0 , R E * , α R E * , p R E * are variables that make S R E ϕ 0 , α , p minimum, which can be denoted as ϕ 0 , R E * , α R E * , p R E * = S R E 1 min S R E ϕ 0 , α , p with the inverse function symbol S R E 1 . Similarly, we have ϕ 0 , R E w * , α R E w * , p R E w * , ϕ 0 , 1 * , α 1 * , p 1 * , ϕ 0 , 1 w * , α 1 w * , p 1 w * , ϕ 0 , 2 * , α 2 * , p 2 * , ϕ 0 , 2 w * , α 2 w * , p 2 w * , ϕ 0 , 12 * , α 12 * , p 12 * , ϕ 0 , 12 w * , α 12 w * , p 12 w * . Corresponding standard deviations can be expressed as:
σ R E = i = 1 n g n w i ϕ i ϕ 0 + α h i p 2 n g 3
σ 1 = i = 1 n g n w i ϕ i ϕ 0 + α h i 2 n g 2
σ 2 = i = 1 n g n w i ϕ i ϕ 0 + α h i 2 2 n g 2
σ 12 = i = 1 n g n w i ϕ i ϕ 0 + α 1 h i + α 2 h i 2 2 n g 3
For no-weighted cases, n w i = 1 . Following criteria are applied in selecting the error estimators:
  • Criterion 1 ( S R E , S R E * in Figure 7)
    Select p * corresponding to min σ R E * , σ R E w *
  • Criterion 2 ( S 2 , S 2 w in Figure 7)
    When p * > 2 , small uncertainty quantities can be obtained due to overestimates of the order of accuracy. As recommended by Eça and Hoekstra [31], when p * > 2 , the error estimator is chosen by comparing the standard deviations resulting from Equations (16)–(19). Then, ε ϕ is calculated with the selected estimator. This procedure substitutes the estimated order of accuracy with 1 or 2 because the exponential power of error estimators to be selected (Equations (11) and (12)) is either 1 or 2. Slightly different from Eça and Hoekstra [31], De Luca et al. [32] substitute p as the theoretical order of accuracy p t h , i.e., p = p t h = 2 , when p * > 2 . In fact, the method by De Luca et al. [32] is more conservative, and it is applied in following calculations. After that, the error estimator can be selected according to min σ 2 * , σ 2 w *
  • Criterion 3 ( S 1 , S 1 w , S 2 , S 2 w , S 12 * , S 12 w * in Figure 7)
    Select the error estimator according to min σ 1 * , σ 1 w * , σ 2 * , σ 2 w * , σ 12 * , σ 12 w *

5. Results and Discussions

With the strategies and numerical methods illustrated in Section 4, systematic simulations are conducted to study the impacts of mesh properties on hydrodynamic performance of the NACA 0012 rudder profile in this section. For each case with different mesh settings, errors with respect to benchmark data [19,25] are displayed in tabular form. Besides, uncertainty estimations are conducted for cases with varying cell sizes. For colored error tables in this section, the color legends are marked based on the values of relative differences. The smallest difference is the darkest green and the largest one is the darkest red for each mesh group, while other colors between the darkest green to red show the transition. The impacts of Reynolds numbers are presented in Section 5.1, and Section 5.2 discusses the selection of domain sizes taking the C-mesh as an example, and impacts of node distributions and mesh types are discussed in Section 5.3 and Section 5.4.

5.1. Impacts of Reynolds Numbers

Before conducting parametric investigations, the Reynolds number should be checked first. The Reynolds number ( R e ) is defined as the ratio of inertial forces to viscous forces, which represents the viscous similarity among flow patterns and can be expressed as:
R e = ρ U L μ
where U is the inflow velocity, L is the characteristic length, and μ is the viscosity coefficient. For rudder hydrodynamics investigations, due to the limited model size and capacity of the test facility, R e of tests in wind tunnels and towing tanks has to be scaled. Analogously, model-scale rudders are preferred by CFD methods since a fine mesh and a large domain are needed for full-scale simulations. In most cases, the similarity law of viscosity and gravity can not be satisfied simultaneously. Generally, experiments or numerical simulations are conducted under a R e which is large enough to guarantee that the turbulence develops sufficiently. A high R e ensures a fully turbulent flow and reduces the effect of the laminar-turbulent transition, achieving an increase of prediction accuracy compared to that of an analysis at low R e with the same mesh.
Cases with larger R e show less calculation efficiency seeing that more cells around boundaries are required. The height of the first row of cells s varies with R e , as shown in Figure 8, while increasing R e leads to decreasing heights. When discretizing domains into multiples mesh elements, the aspect ratio for a single element should be around 1 to ensure mesh quality. When it comes to cells near walls, aspect ratios are connected with s and chord-wise spacings c . As shown in Figure 9, more nodes along the wall are distributed with smaller s to keep aspect ratios compliant with the requirement of solvers. Thus, a decrease in R e leads to a decrease in the number of cells. To sum up, we suggest that using a minimum R e which is achievable in experiments, ensures full-turbulent conditions, and can be solved with fewer mesh elements.
To determine the threshold of the Reynolds number, above which the turbulence develops completely and force coefficients do not change significantly, cases with varying R e are simulated in this section. Based on the flat-plate boundary layer theory and the calculation method in [33], the height of meshes’ first layer s above the rudder surface depends on y + and R e . According to the application scope of the k- ω SST model, y + is set as 1 for all cases, and corresponding heights are shown in Figure 8.
Common benchmark wind tunnel tests are carried out at R e in the range of 1.00 × 10 5 [34] to 1.00 × 10 7 [25]. As low-Reynolds-number RANS analysis is still challenging [35,36] and high-Reynolds-number simulations may be expensive in computation time, the present work tests the NACA 0012 profile at R e in a range of 2.00 × 10 5 1.89 × 10 7 . Among these aerodynamic experimental results, limited cases with R e = 3.94 × 10 6 1.89 × 10 7 and a small Mach number can be taken as validation results of incompressible water simulations, as the compressibility effects of a fluid with a Mach number smaller than 0.20 is small. Table 4 presents the case settings and results of various R e under distinct α series, as specific values of α are listed in [25].
Figure 10 presents calculation results of hydrodynamic coefficients versus the same rudder angles for different R e , which indicates that the lift curve rises with an increase of the R e while the drag curve decreases. Compared to that of the lift coefficients, the drag coefficients are more sensitive to changes in R e . The drag coefficient under 15 at R e = 1.89 × 10 7 is about half of the value at R e = 4.00 × 10 5 , whereas the lift coefficient is 1.20 times larger. The differences in lift and drag forces between low and high R e increase with an increasing angle of attack. Consistent with findings by Ladson [25] and Molland and Turnock [37], R e = 6.00 × 10 6 can be considered as a threshold value for the mesh generation, above which little variation may be found. Moreover, the y + along the chord for four cases under small and large rudder angles with varying R e are shown in Figure 11. y + values for the case with R e > 5.97 × 10 6 are roughly between 0 and 1, indicating that the k- ω SST turbulence model is applicable for the case. Since ship rudders in propeller slipstream usually experience full turbulence with large R e , low R e conditions are not of interest in current studies. In this paper, we set R e of the following simulations as 6.00 × 10 6 .

5.2. Impacts of Domain Sizes

A C-mesh is applied to investigate the effects of the domain size on RANS solutions by conducting simulations with various sizes against varying cell numbers in this section. Test parameters and related relative differences are given in Table 5. C i n and C o u t are doubled for each case compared with that of the previous one, while the narrowest domain extends 7.5 c upstream and 15 c downstream, which is 16 times smaller than the largest ones, i.e., 120 c and 240 c . Since equivalent mesh nodes correspond to different spacing variations in a variety of domains, different nodes distributions are applied for each domain to determine reasonable meshing strategies.
Pressure distributions in the whole domain and around the rudder are presented in Figure 12. From Figure 12a–c, regions with great gradient changes near the rudder can be observed even for the smallest domain. Comparing Figure 12d–f, neither high-pressure regions near the stagnation points nor regions near the trailing edges show significant distinctions. The impact of domain size on pressure distributions around the rudder is relatively small, which corresponds to minor differences of C L mainly derived from pressure integration.
Compared with that of C L , accuracy of C D is more sensitive to domain sizes in Table 5. With an increase of the domain size, the accuracy of the prediction of C D is more significantly improved than C L . Small domain sizes show poor accuracy in C D but obtained C L is acceptable for engineering applications. The relative differences of C D for α = 15 . 20 drop from 55.19% to 34.79%, indicating that boundary effects can evidently change viscosity components around the rudder, which mainly contribute to drag forces. On the other hand, larger domains, like 60 c / 120 c and 120 c / 240 c , tend to show high and stable precision with more cells. However, further expansion of the domain size may dramatically increase the number of cells and computation time.
From the perspective of efficiency, a domain size that eliminates boundary effects and meanwhile requires fewer mesh nodes is desired for following investigations. Table 5 shows that a domain size of C i n = 30 c and C o u t = 60 c with about 2.50 × 10 5 cells achieves a good balance in both accuracy and efficiency, and Table 6 shows that the validation process is achieved for C L obtained from the case ( ϕ 2 ). This domain size is large enough to obtain sufficiently accurate estimates of lift and drag, and the number of cells is still manageable by common desktop computers. Thus, it is applied for the remainder of the simulations in this paper. Table 5 indicates that even more accurate predictions can be obtained with larger domains. Thus, we suggest applying a large domain when computation power and time are available. If the lift is the only purpose, a small domain size of 15 c / 30 c may be more favorable as it requires less computation power and time than a large domain.

5.3. Impacts of Node Distributions

After determining the test R e and the domain size, raw meshes are ready for the mesh independence study for different mesh types. The common procedure is to carry out simulations on an initial mesh with a residual error in a range of 1.00 × 10 4 1.00 × 10 5 . After that, refine the mesh globally to around 2 or 2 times the initial mesh. Next, run simulations with the refined mesh and compare solutions obtained from the coarse mesh and the fine mesh. Repeat the refinement until the results do not significantly change with a finer mesh. The mesh independence study requires at least 3 solutions to evaluate the convergence of certain inputs [26]. Considering the computation time, it is always better to use the smallest number of cells. In this paper, mesh resolutions are changed by conducting local rather than global refinements, while the grid refinement factor R G is still around 2 for target parameters.
Since the flow near boundaries of the rudder profile changes more violently than that in the far-field and has larger effects on rudder hydrodynamics, near-field node distributions are studied by fixing far-field node distributions for four mesh types. After that, the impacts of far-field distributions are studied with proper near-field distributions.
In this section, with a domain size of 30 c / 60 c , impacts of near-field node distributions are investigated for four mesh types, as shown in Figure 13, Figure 14, Figure 15 and Figure 16. Meshes in different regions are categorized into two parts, which are located in the far-field and near the rudder surface. Near-field meshes are defined as those located along the rudder chord and within the boundary layers regions. The domain for C-mesh in Figure 13 consists of a semicircle and a rectangle, and nodes distributed parallel to the radius are defined as in radial directions, and the rudder in O-mesh is surrounded by radial nodes in Figure 15. The nodes extending from the trailing edge to the outlet of the boundary are defined as in wake directions. Patterned after the definitions above, node distributions for H-mesh are shown in Figure 14, and nodes from the leading edge to the inlet are in inlet directions. For Hybrid-mesh in Figure 16, considering quite different mesh generation methods compared with that of structured meshes, meshes in the refinement area near the rudder are treated as near-field meshes while others are in far-field.

5.3.1. Impacts of Near-Field Node Distributions

In this part, impacts of near-field node distributions are studied by changing chord-wise spacings and layer-wise growth rates for C-mesh, H-mesh, O-mesh, and Hybrid-mesh. The near-chord refinement is performed by changing the number of nodes in the chord-wise direction for structured meshes and chord-wise spacings for Hybrid-mesh with R G around 2 . For structured meshes, a reference setting is coarsened or refined to perform parametric studies. The quantitative relation between a coarsened or refined case and the reference case can be expressed as:
log R G N c i N c 0 = log R G N c i 0 = i , ( i = 2 , 1 , )
where i denotes the level of refinement with respect to the reference case, and N c i is the number of nodes along the chord. Taking C-mesh for example, the rudder is wrapped by denser cells in chord-wise directions, as shown in Figure 13. As for cells in the boundary layer region of the rudder, much attention should be paid to capture violently changing velocity gradients near the wall. The layer-wise growth rate is investigated to provide a detailed view of mesh resolutions in boundary layer regions.
Cells around the rudder with varying chord-wise spacings for C-mesh are shown in Figure 17. As the chord-wise spacing decreases in Table 7, the number of cells with large aspect ratios in the wake increases, which introduces the difficulty of convergence, especially under large angles of attack. The accuracy of C D is slightly improved with denser chord-wise node distributions, while variations of C L are even smaller. Larger deviations are obtained for C M for small α , probably due to measurement errors in experiments induced by small absolute values. For α = 4 . 18 and α = 8 . 22 , better mesh resolution along the rudder lead to lower precision, which are related to convergence problems due to the fact that small flow features are being solved. Since calculation results tend to stabilize without diffusing after log R G N c i 0 = 2 and 5.48 × 10 5 mesh quadrilaterals show relatively high accuracy, the case with log R G N c i 0 = 3 is selected for the following simulations on C-mesh, and corresponded uncertainty estimations are shown in Table 8.
For H-mesh, besides chord-wise nodes, nodes distributed along with radial directions in the refinement region also count in near-field node distributions. Therefore, chord-wise nodes and radial nodes near the rudder are both refined for H-mesh. Due to discrepancies in refinement regions in Table 9, the NOC of H has a more rapid growth rate compared with C-mesh. Generally, variations of the accuracy of hydrodynamic coefficients are consistent with the increase in the number of cells, and relatively greater improvements can be observed for C D and C M taking log R G N c i 0 = 2 as the dividing line. For the NOC less than 3.00 × 10 5 , C-mesh is a better choice for predictions of α < 8 . 22 when C-mesh and H-mesh share the same order of magnitude. Considering that near-field node distributions of H-mesh are able to capture more precise flow details around the rudder with denser near-field node distributions, H-mesh with fine mesh resolutions can be applied with enough computational resources. In this part, to balance the accuracy and efficiency, log R G N c i 0 = 2 is selected for following simulations on H-mesh, and corresponded uncertainty estimations are shown in Table 10.
For O-mesh in Table 11, the errors of C L and C D are larger than C and H-mesh because the mesh type fails to capture the wake after the trailing edge. Further, variations of layer-wise growth rates are incapable of improving the situation. Therefore, from the perspective of efficiency, further investigations are not carried out on O-mesh considering its obvious defects.
Various chord-wise spacings are applied for Hybrid-mesh, and the results are in Table 12. Compared with that of structured meshes, the calculation accuracy of C L is as good as C and H-mesh, but that of C D is worse than C and H-mesh while still better than O-mesh. C L tends to be stable when the NOC is more than 2.5 × 10 5 , and C D results are gradually improved with smaller chord-wise spacings. In the following part, a chord-wise spacing of 2.23 × 10 4 c is selected, and validation for C L is achieved as indicated in Table 13.
No obvious variations can be observed for both C L , C D and C M in Table 14. The insensitivity to layer-wise growth rates for C-mesh is probably caused by good chord-wise mesh resolutions. A rate of 1.15 encounters convergence problems, while results obtained from r 1.10 are stable. Therefore, r = 1.10 is selected.
For H-mesh in Table 15, the accuracy of C L and C D has subtle improvements within 1%, while variability of C M is slightly larger. The trend shows that a slower growth rate of mesh height leads to better results, but, remarkably, too small rates may reduce calculation efficiencies and contribute to mesh quadrilaterals with too large aspect ratios, which are disadvantageous to the convergence of CFD results. In this case, r = 1.05 is well balanced in accuracy and computation time.
For Hybrid-mesh in Table 16, the connection between the structured region and the unstructured region is influenced by layer-wise growth rates. When r = 1.03 , relatively large aspect ratios in boundary layer regions increases the difficulty of convergence. Therefore, r = 1.10 is selected as a value of transition of convergent conditions.

5.3.2. Impacts of Far-Field Node Distributions

With selected near-field node distributions, nodes in far-field regions need to be determined to match variations of cell sizes near wall boundaries to present global meshing strategies. For structured meshes, similar to Section 5.3.1, nodes in different directions are changed according to characters of their topological structures, which can be expressed as:
log R G N r i N r 0 = log R G N r i 0 = i , ( i = 2 , 1 , , 4 )
log R G N w i N w 0 = log R G N w i 0 = i , ( i = 2 , 1 , , 4 )
where N r i and N w i denote the number of nodes in radial and wake directions. For Hybrid-mesh, cell sizes in far-field regions are varied to change mesh properties.
Number of cells against radial nodes vary over a considerable extent from nearly 3.00 × 10 5 to more than 1.00 × 10 6 for C-mesh and H-mesh, as shown in Table 17 and Table 18. However, such variation has few impacts on magnitudes of any hydrodynamic coefficients, indicating that sufficient mesh resolution along the rudder chord has successfully solved the flow in regions with great gradients. In other words, relatively few nodes in radial directions are capable of ensuring the calculation accuracy. Therefore, the reference case already achieved mesh independence for both C-mesh and H-mesh along with radial directions, and the case is still utilized for benchmark in the following studies.
Similar to impacts of nodes in radial directions, variations of C L , C D and C M against changing cell numbers induced by cell sizes in wake directions are very small, as shown in Table 19 and Table 20. Compared with the trend of accuracy variations in Section 5.3.1, impacts of node distributions in the near-field on CFD results are stronger than that in the far-field, which is reasonable considering the decaying disturbance of the object on the flow field with increased distances from the wall boundary. As long as mesh resolutions in the near-field are handled carefully, fewer mesh elements are needed in the far-field. Mesh elements in unstructured regions for Hybrid-mesh are changed as shown in Table 21, and the magnitude of results variations are the same as C-mesh and H-mesh. Further decreasing element sizes of unstructured meshes are memory-consuming, and divergent cases occur when the size is below 0.0625 c . Therefore, the body element size of 0.125 c is selected.
Table 22, Table 23, Table 24, Table 25 and Table 26 show uncertainty estimations for C-mesh, H-mesh and Hybrid-mesh with varying far-field node distributions. Most uncertainty quantities in Table 22, Table 23, Table 24 and Table 25 are relatively small, corresponding to the fact that further increase in radial or wake directions has little impact on the calculation accuracy.

5.4. Impacts of Mesh Types

Simulations of different mesh types with varying mesh properties are conducted in Section 5.1Section 5.3. Generally, the accuracy of C L is sufficient (within 5 % deviations compared with experimental data) for engineering applications except for O-mesh, while discrepancies of C D between CFD and EFD results are relatively large under high angles of attack. Physically, the lift and drag coefficient are nondimensional forms of lift and drag force, respectively, while the former is mainly generated by the pressure difference between the windward and leeward surfaces, the inviscid component, and the latter primarily contributed by the frictional shear stress, the viscous component. The precise prediction of viscous stress relies on a more exact forecasting model in regions of separated flow, along with an advanced understanding of rudder hydrodynamics, the development of the discretization, and numerical methods [38]. In other words, it is inevitable to observe excessive error of C D with the RANS method and limited computational resources. However, selecting the appropriate mesh type can improve the accuracy of both C L and C D to a certain extent in the current framework. Therefore, based on the simulation results obtained above, the impacts of different mesh types on rudder hydrodynamics are summarized in this section. To evaluate the efficiency and accuracy of those mesh types, the calculation accuracy of C D versus changing numbers of cells in the near-field and far-field is shown in Figure 18 and Figure 19, while blue dotted lines connect minimum errors under each α and red lines mark selected mesh settings for each type. Since there is a common trend between impacts of radial and wake distributions on calculation accuracy for C-mesh and H-mesh, only radial node distributions are selected to perform the estimation.
In Figure 18, for the NOC below 3.00 × 10 5 , C-mesh shows higher accuracy than other mesh types. However, H-mesh with better mesh resolution around the rudder achieves better performance under large angles of attack. What’s more, it can be foreseen that a further increase in nodes may lead to more accurate results if computational power is not an issue. According to the tendency of Hybrid-mesh, the accuracy of C D remained stable, and smaller chord-wise spacing cannot make up for defects of unstructured mesh in solving rudder flows.
In Figure 19, even though variations of C D of Hybrid-mesh are slightly larger than C-mesh and H-mesh, results obtained from these mesh types settle in stable levels with fixed chord-wise spacings. Although making further efforts in far-field mesh resolutions can still improve the accuracy, the consequent mesh elements with small size may dramatically increase computation time, as well as bring difficulties with convergence. C-mesh and H-mesh with NOC around 5.00 × 10 5 are selected, which are about 2.00 × 10 5 less than Hybrid-mesh in the number of cells and show higher accuracy especially under high rudder angles, proving that structured meshes have higher efficiency in computations than hybrid meshes.
Based on the principle of balancing efficiency and accuracy, recommended mesh properties for four mesh types investigated in Section 5.1, Section 5.2 and Section 5.3 are presented in Table 27, and performance of these settings are compared with that of other codes (Langley Research Center [19]) in Table 28, which coincide with H-mesh best. The results of these CFD codes are achieved using the same block-structured C-mesh in a domain of 500 c around the profile and 500 c behind it, which is much larger than the domain in this paper. Hence, the appearance of differences in Table 27 is reasonable, and relatively good results are achieved by the current method efficiently. Structured meshes of C-type and H-type can achieve accurate C L predictions, while O-mesh fails to capture flow features in wake directions because of the difficulties of actualizing clustering grid points in the wake while maintaining a reasonable spacing at the trailing edge [5]. The concentration degree of mesh elements in regions with a large change of solution gradients for H-mesh is greater than that of C-mesh, contributing to higher accuracy in C D predictions when the NOC is above 3.00 × 10 5 . Some unexpected divergent cases can be observed for C-mesh and H-mesh, and no general rules can be derived to explain why that happens, though it probably relates to the processing mechanism of the solver. In systematic parameter studies, it is recommended to fine-tune mesh properties to obtain convergent results. Hybrid-mesh can be generated more simply with less manual intervention, and calculation results are acceptable but still inferior to C-mesh and H-mesh. However, the high adaptability of unstructured meshes to complex geometries extends the scope of application of Hybrid-mesh, which means that such a meshing strategy can be easily transformed into sophisticated configurations like high-lift rudders.

6. Further Validations

To further validate the proposed mesh strategy, more 2D rudder profiles are studied in this section, i.e., NACA 0006, NACA 0009, NACA 0015, and NACA 0018, and the calculated force coefficients are compared with corresponding experimental data. The test conditions are shown in Table 29.
By applying H-mesh and Hybrid-mesh with recommended domain size and nodes distributions, obtained C L and C D are compared with EFD results before stall for tested rudder profiles as shown in Figure 20. For profiles with relatively small thickness ratios, such as NACA 0006 and NACA 0009, Hybrid-mesh failed to obtain credible force coefficients. As for H-mesh, in such cases, the mesh type can achieve acceptable C L calculations, although the scope of α which leads to convergent solutions is limited. Nevertheless, from the perspective of structural strength, thin profiles are generally not used for rudder designs, so no further efforts are made to modify current mesh strategies to adapt thin profiles in this paper. For NACA 0015 and NACA 0018, when α < 5 , some divergence may occur for H-mesh, while Hybrid-mesh shows good convergence. After that, CFD results for H-mesh and Hybrid-mesh agree with EFD ones well, while H-mesh shows better accuracy for α around 15 , demonstrating that the proposed meshing strategy can be used for other rudder profiles.

7. Conclusions

This paper describes a comprehensive approach to identifying suitable mesh properties for CFD simulations of airfoil-shaped ship rudders. The goal is to build up an efficient mesh with a sufficient resolution. The suggested mesh properties are summarized in Table 30 for different research objectives. Important conclusions drawn from this work include:
  • The mesh type affects the efficiency of mesh generation and the accuracy of the CFD solutions. The structured mesh gives more accurate drag predictions with fewer cells compared with that of the hybrid mesh. However, an unstructured mesh is more suitable for complex rudder profiles; it can also provide comparable accurate lift predictions with the structured mesh.
  • The domain size has to be reasonably large to fully develop the flow. The drag prediction is more influenced by the domain size than the lift prediction. An increase in domain size leads to an improvement in prediction accuracy.
  • Lift coefficients increase while drag coefficients decrease with an increase of Reynolds numbers in a range of 2 × 10 5 1.89 × 10 7 . The drag coefficient is more sensitive to changes in the Reynolds number than the lift coefficient. A Reynolds number of 6.00 × 10 6 is suggested as a benchmark, above which few changes in lift and drag coefficients were found.
  • The boundary conditions applied in this study are velocity-inlet, no-slip wall, and pressure-outlet. RANS simulation results with these conditions are in good agreement with the published experimental data [25,39,40].
  • In all cases, C D prediction accuracy remains poor compared with that of C L .
We realize that the above conclusions are based on a case study of 2D ship rudders. For unconventional rudder types, such as flap rudders and fishtail rudders, the conclusions should not be applied directly. For instance, the gap of the flap rudder and the concave part of the fishtail rudder may require an additional number of cells. What’s more, attention should be paid to mesh elements in span-wise directions for 3D rudders. Considering larger calculation errors with larger angles of attack, extra efforts are also needed in applications of other CFD methods, such as Large-Eddy Simulation (LES) and Direct Numerical Simulation (DNS), which commonly require a much finer mesh than the RANS method. Therefore, further study is required to determine the most suitable mesh for different rudder types and meshing strategies in 3D configurations.

Author Contributions

Conceptualization, S.L. and J.L.; methodology, S.L. and J.L.; software, S.L. and J.L.; validation, S.L.; formal analysis, R.H.; investigation, R.H.; resources, J.L.; data curation, S.L.; writing—original draft preparation, S.L. and J.L.; writing—review and editing, R.H.; visualization, S.L.; supervision, J.L.; project administration, J.L.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (62003250), Research on Intelligent Ship Testing and Verification (2018473), Program of Marine Economy Development Special Fund (Six Marine Industries) under Department of Natural Resources of Guangdong Province (Project No. GDNRC[2021]59).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting reported results can be found in Ladson 1988, Abbott 1959, Jacobs 1937.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AVEAverage error
NOCNumber of cells
GCIGrid convergence index
CFCorrection factor
α Angle of attack of the rudder
s Height of the meshes’ first layer
c Chord-wise spacing around a rudder profile
C L Relative difference of the calculated lift coefficient and benchmark data
C D Relative difference of the calculated drag coefficient and benchmark data
C M Relative difference of the calculated moment coefficient and benchmark data
δ Thickness of a flat plate boundary layer
EComparison error
ϕ i Solution of ith mesh
cChord length of a rudder profile
N c Number of nodes in chord-wise
N r Number of nodes in radial directions
N w Number of nodes in wake directions
pOrder of accuracy
p t h Theoretical order of accuracy
R e Reynolds number
R G Grid refinement factor
rLayer-wise growth rate in boundary layer regions
U G Grid uncertainty
U V Validation uncertainty
y + Dimensionless wall distance

References

  1. Liu, J.; Hekkenberg, R. Sixty years of research on ship rudders: Effects of design choices on rudder performance. Ships Offshore Struct. 2017, 12, 495–512. [Google Scholar] [CrossRef] [Green Version]
  2. Liu, J.; Quadvlieg, F.; Hekkenberg, R. Impacts of the rudder profile on manoeuvring performance of ships. Ocean. Eng. 2016, 124, 226–240. [Google Scholar] [CrossRef] [Green Version]
  3. Badoe, C.E.; Phillips, A.B.; Turnock, S.R. Influence of drift angle on the computation of hull-propeller-rudder interaction. Ocean. Eng. 2015, 103, 64–77. [Google Scholar] [CrossRef] [Green Version]
  4. Van Nguyen, T.; Ikeda, Y. Hydrodynamic characteristic of rudder sections with high lift force, Part 2: The wedge tail shapes. J. Jpn. Soc. Nav. Archit. Ocean. Eng. 2014, 171–174. [Google Scholar]
  5. Lutton, M.J. Comparison of C-and O-grid Generation Methods Using a NACA 0012 Airfoil. Master’s Thesis, Air Force Institute of Technology, Kaduna, Nigeria, 1989. [Google Scholar]
  6. Basha, W. Accurate drag prediction for transitional external flow over airfoils. Master’s Thesis, Concordia University, Montreal, QC, Canada, 2006. [Google Scholar]
  7. Stuck, A.; Turnock, S.; Bressloff, N. An Evaluation of the RANS Method for the Prediction of Steady Ship Rudder Performance Compared to Wind Tunnel Measurements; Technical Report; University of Southampton: Southampton, UK, 2004. [Google Scholar]
  8. Turnock, S.R.; Holroyd, N.J.; Date, J.C. Appendage design for the America’s Cup using CFD. In Proceedings of the European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS–CFD 2001), Swansea/Wales, UK, 4–7 September 2001. [Google Scholar]
  9. Kim, H.J.; Kim, S.H.; Oh, J.K.; Seo, D.W. A proposal on standard rudder device design procedure by investigation of rudder design process at major Korean shipyards. J. Mar. Sci. Technol. 2012, 20, 450–458. [Google Scholar]
  10. Fujii, H. Experimental researches on rudder performance (1). J. Zosen Kiokai 1960, 107, 105–111. (In Japanese) [Google Scholar] [CrossRef]
  11. Fujii, H.; Tsuda, T. Experimental researches on rudder performance (2). J. Zosen Kiokai 1961, 109, 31–42. (In Japanese) [Google Scholar] [CrossRef]
  12. Fujii, H.; Tsuda, T. Experimental researches on rudder performance (3). J. Zosen Kiokai 1962, 111, 51–58. (In Japanese) [Google Scholar] [CrossRef] [Green Version]
  13. Nienhuis, U. Passieve Manoeuvreerhulpmiddelen: Open Water Proeven met Roer; Technical Report; Maritime Research Institute Netherlands (MARIN): Wageningen, The Netherlands, 1987. (In Dutch) [Google Scholar]
  14. Liu, J.; Hekkenberg, R.; Quadvlieg, F.; Hopman, H.; Bingqian, Z. An integrated empirical manoeuvring model for inland vessels. Ocean. Eng. 2017, 137, 287–308. [Google Scholar] [CrossRef]
  15. Wasberg, C.E.; Reif, B.A.P. Hydrodynamical Simulations in Fluent; Technical Report; Norwegian Defence Research Establishment (FFI): Kjeller, Norway, 2010. [Google Scholar]
  16. Tveiterås, V. Numerical Study of the Interaction of Flow over Two Airfoils in Relative Motion. Master’s Thesis, Norwegian University of Science and Technology, Trondheim, Norway, 2011. [Google Scholar]
  17. Narsipur, S.; Pomeroy, B.; Selig, M. CFD analysis of multielement airfoils for wind turbines. In Proceedings of the 30th AIAA Applied Aerodynamics Conference, New Orleans, LA, USA, 25– 28 June 2012; p. 2781. [Google Scholar]
  18. Eleni, D.C.; Athanasios, T.I.; Dionissios, M.P. Evaluation of the turbulence models for the simulation of the flow over a National Advisory Committee for Aeronautics (NACA) 0012 airfoil. J. Mech. Eng. Res. 2012, 4, 100–111. [Google Scholar]
  19. Langley Research Center. Turbulence Modeling Resource. Available online: http://turbmodels.larc.nasa.gov/naca0012_val.html (accessed on 10 September 2014).
  20. Nordanger, K.; Holdahl, R.; Kvarving, A.M.; Rasheed, A.; Kvamsdal, T. Implementation and comparison of three isogeometric Navier-Stokes solvers applied to simulation of flow past a fixed 2D NACA 0012 airfoil at high Reynolds number. Comput. Methods Appl. Mech. Eng. 2015, 284, 664–688. [Google Scholar] [CrossRef] [Green Version]
  21. Siddiqui, M.S.; Rasheed, A.; Kvamsdal, T. Numerical assessment of RANS turbulence models for the development of data driven Reduced Order Models. Ocean. Eng. 2020, 196, 106799. [Google Scholar] [CrossRef]
  22. Liseikin, V.D. Grid Generation Methods, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  23. Aftosmis, M.; Gaitonde, D.; Tavares, T.S. Behavior of linear reconstruction techniques on unstructured meshes. AIAA J. 1995, 33, 1–51. [Google Scholar] [CrossRef]
  24. Stern, F.; Wilson, R.V.; Coleman, H.W.; Paterson, E.G. Verification and Validation of CFD Simulations; Technical Report; Iowa Institute of Hydraulic Research: Iowa City, IA, USA, 1999. [Google Scholar]
  25. Ladson, C.L. Effects of Independent Variation of Mach and Reynolds Numbers on the Low-Speed Aerodynamic Characteristics of the NACA 0012 Airfoil Section; Technical Report; Langley Research Center: Hampton, VA, USA, 1988. [Google Scholar]
  26. International Towing Tank Conference. ITTC Recommended Procedures and Guidelines 7.5-03-01-01: Uncertainty Analysis in CFD Verification and Validation. 2008. Available online: https://ittc.info/media/4184/75-03-01-01.pdf (accessed on 5 September 2021).
  27. Islam, H.; Soares, C.G. Uncertainty analysis in ship resistance prediction using OpenFOAM. Ocean. Eng. 2019, 191, 105805. [Google Scholar] [CrossRef]
  28. Roache, P.J. Verification and Validation in Computational Science and Engineering; Hermosa: Albuquerque, NM, USA, 1998; Volume 895. [Google Scholar]
  29. Celik, I.B.; Ghia, U.; Roache, P.J.; Freitas, C.J.; Coleman, H.; Raad, P.E. Procedure for estimation and reporting of uncertainty due to discretization in CFD applications. J. Fluids Eng. 2008, 130. [Google Scholar]
  30. Stern, F.; Wilson, R.; Shao, J. Quantitative V&V of CFD simulations and certification of CFD codes. Int. J. Numer. Methods Fluids 2006, 50, 1335–1355. [Google Scholar]
  31. Eça, L.; Hoekstra, M. A procedure for the estimation of the numerical uncertainty of CFD calculations based on grid refinement studies. J. Comput. Phys. 2014, 262, 104–130. [Google Scholar] [CrossRef]
  32. De Luca, F.; Mancini, S.; Miranda, S.; Pensa, C. An extended verification and validation study of CFD simulations for planing hulls. J. Ship Res. 2016, 60, 101–118. [Google Scholar] [CrossRef]
  33. Pointwise. Y+ Calculator. 2020. Available online: http://www.pointwise.com/yplus/ (accessed on 5 September 2021).
  34. Thieme, H. Design of Ship Rudders (Zur Formgebung von Shiffsrudern); Technical Report; Shipbuilding Institute, University of Hamburg: Washington, DC, USA, 1965. [Google Scholar]
  35. Catalano, P.; Tognaccini, R. RANS Analysis of the low-Reynolds number flow around the SD7003 airfoil. Aerosp. Sci. Technol. 2011, 15, 615–626. [Google Scholar] [CrossRef]
  36. Tang, L. Reynolds-Averaged Navier-Stokes simulation of low-Reynolds-number airfoil aerodynamics. J. Aircr. 2008, 45, 848–856. [Google Scholar] [CrossRef]
  37. Molland, A.F.; Turnock, S.R. Marine Rudders and Control Surfaces: Principles, Data, Design and Applications, 1st ed.; Elsevier Butterworth-Heinemann: Oxford, UK, 2007. [Google Scholar]
  38. Roberts, T.W. Prospects for drag prediction using computational fluid dynamics. In Proceedings of the Aerospace Technology Conference and Exposition, Long Beach, CA, USA, 1 October 1990. [Google Scholar]
  39. Abbott, I.H.; Doenhoff, A.E.V. Theory of Wing Sections. Including a Summary of airfoil Data; Dover Publications, Inc.: New York, NY, USA, 1959. [Google Scholar]
  40. Jacobs, E.N.; Sherman, A. Airfoil Section Characteristics as Affected by Variations of the Reynolds Number; Technical Report Archive and Image Library; UNT Libraries Government Documents Department: Denton, TX, USA, 1937. [Google Scholar]
Figure 1. Comparison of the hydrodynamic coefficients of a rudder in open water conditions to those of the rudder in propeller slipstream. Data are adopted from Nienhuis [13].
Figure 1. Comparison of the hydrodynamic coefficients of a rudder in open water conditions to those of the rudder in propeller slipstream. Data are adopted from Nienhuis [13].
Jmse 09 01062 g001
Figure 2. Original and modified mesh types for the block-structured mesh.
Figure 2. Original and modified mesh types for the block-structured mesh.
Jmse 09 01062 g002
Figure 3. The mesh type of Hybrid-mesh.
Figure 3. The mesh type of Hybrid-mesh.
Jmse 09 01062 g003
Figure 4. The geometric parameters for different mesh types.
Figure 4. The geometric parameters for different mesh types.
Jmse 09 01062 g004
Figure 5. Relationship between different mesh properties.
Figure 5. Relationship between different mesh properties.
Jmse 09 01062 g005
Figure 6. Strategy applied to study different mesh properties.
Figure 6. Strategy applied to study different mesh properties.
Jmse 09 01062 g006
Figure 7. Procedure for error and uncertainty estimations by GCI-LSR method recommended by Eça and Hoekstra [31], De Luca et al. [32].
Figure 7. Procedure for error and uncertainty estimations by GCI-LSR method recommended by Eça and Hoekstra [31], De Luca et al. [32].
Jmse 09 01062 g007
Figure 8. Different meshes’ height of first layer with varying R e and y + = 1 .
Figure 8. Different meshes’ height of first layer with varying R e and y + = 1 .
Jmse 09 01062 g008
Figure 9. Impacts of height of first mesh s on aspects ratios.
Figure 9. Impacts of height of first mesh s on aspects ratios.
Jmse 09 01062 g009
Figure 10. Impacts of R e on lift and drag coefficients.
Figure 10. Impacts of R e on lift and drag coefficients.
Jmse 09 01062 g010
Figure 11. y + along rudder chord with varying R e .
Figure 11. y + along rudder chord with varying R e .
Jmse 09 01062 g011
Figure 12. Pressure contours of different domain sizes with α = 15 . 20 .
Figure 12. Pressure contours of different domain sizes with α = 15 . 20 .
Jmse 09 01062 g012
Figure 13. Node distributions for C-mesh.
Figure 13. Node distributions for C-mesh.
Jmse 09 01062 g013
Figure 14. Node distributions for H-mesh.
Figure 14. Node distributions for H-mesh.
Jmse 09 01062 g014
Figure 15. Node distributions for O-mesh.
Figure 15. Node distributions for O-mesh.
Jmse 09 01062 g015
Figure 16. Node distributions for Hybrid-mesh.
Figure 16. Node distributions for Hybrid-mesh.
Jmse 09 01062 g016
Figure 17. Cells around the rudder with varying chord-wise spacings for C-mesh.
Figure 17. Cells around the rudder with varying chord-wise spacings for C-mesh.
Jmse 09 01062 g017
Figure 18. Calculation deviations of C D compared with experimental benchmark data adopted from [25] varying with numbers of cells, while changing chord-wise spacings in the near-field (blue dotted lines connect minimum errors under each α and red lines mark selected mesh settings for each type).
Figure 18. Calculation deviations of C D compared with experimental benchmark data adopted from [25] varying with numbers of cells, while changing chord-wise spacings in the near-field (blue dotted lines connect minimum errors under each α and red lines mark selected mesh settings for each type).
Jmse 09 01062 g018
Figure 19. Calculation deviations of C D compared with experimental benchmark data adopted from [25] varying with numbers of cells, while changing radial node distributions in near-field (blue dotted lines connect minimum errors under each α and red lines mark selected mesh settings for each type).
Figure 19. Calculation deviations of C D compared with experimental benchmark data adopted from [25] varying with numbers of cells, while changing radial node distributions in near-field (blue dotted lines connect minimum errors under each α and red lines mark selected mesh settings for each type).
Jmse 09 01062 g019
Figure 20. Comparison of results of H and Hybrid mesh with experimental benchmark data adopted from [39,40] for NACA 0006, NACA 0009, NACA 0015 and NACA 0018.
Figure 20. Comparison of results of H and Hybrid mesh with experimental benchmark data adopted from [39,40] for NACA 0006, NACA 0009, NACA 0015 and NACA 0018.
Jmse 09 01062 g020
Table 1. Mesh properties studies in literature: objects, Reynolds numbers, mesh types, and domain sizes.
Table 1. Mesh properties studies in literature: objects, Reynolds numbers, mesh types, and domain sizes.
No.SettingsObject Re Mesh TypeDomain Size
LiteratureTestedSelected
1Lutton [5]NACA 0012 2.00 × 10 6 C, O
2Turnock et al. [8]NACA 0012 6.00 × 10 6 Outlet 15 c
3Stuck et al. [7]NACA 0012 6.00 × 10 6 Hybrid C,
Hybrid triangle and quadrilateral
Hybrid quadInlet 15 c
Outlet 15 c
4Basha [6]NACA 0012 2.00 × 10 6 O, C, HybridHybrid
5Wasberg and Reif [15]NACA 0009,
NACA 65209
7.04 × 10 6 HHInlet 29 c
Outlet 30 c
6Tveiterås [16]NACA 0012 2.00 × 10 6 ,
3.00 × 10 6
Hybrid, CHybridOutlet 19 c
7Narsipur et al. [17]Multi-element
airfoils
2.51 × 10 6 HybridHybrid
8Eleni et al. [18]NACA 0012 3.00 × 10 6 CCSide 20 c
Outlet 30 c
9Van Nguyen and Ikeda [4]NACA 0024 6.00 × 10 6 HybridHybrid
10Langley Research Center [19]NACA 0012 6.00 × 10 6 CCOutlet 500 c
11Nordanger et al. [20]NACA 0012 3.00 × 10 6 Multiple patchMultiple patchInlet 15 c
Outlet 20 c
12Siddiqui et al. [21]NACA 0015 1.96 × 10 6 CCInlet 8 c
Outlet 20 c
Table 2. Mesh properties studies in literature: node distributions and turbulence models.
Table 2. Mesh properties studies in literature: node distributions and turbulence models.
No.SettingsNodes DistributionTurbulence Model
LiteratureChordBoundary LayerFar Field
1Lutton [5]
2Turnock et al. [8]244 nodes y + = 30 49 radial cells
99 wake cells
k- ϵ standard and RNG
3Stuck et al. [7]200 nodes t 3 δ k- ϵ RNG
4Basha [6]60 nodes y + < 5 SA
5Wasberg and Reif [15] y + 1
s = 1.50 2.00 × 10 5
SA
6Tveiterås [16]340 nodes y + = 1.7 3
s = 3.00 × 10 5 , r = 1.14
Transition SST
7Narsipur et al. [17] y + 1 k- ω SST
8Eleni et al. [18] y + 0.2
s = 1.00 × 10 5
k- ω SST
9Van Nguyen and Ikeda [4] y + 1 k- ω SST
10Langley Research Center [19]
11Nordanger et al. [20]127 nodes y + = 30 SA
12Siddiqui et al. [21]k- ω SST
Table 3. Numerical settings for RANS simulations with ANSYS.
Table 3. Numerical settings for RANS simulations with ANSYS.
PropertySetting
Turbulence modelk- ω SST
Pressure-velocity couplingCoupled
GradientLeast squares cell based
PressureSecond order
MomentumSecond order upwind
Turbulent kinetic energySecond order upwind
Specific dissipation rateSecond order upwind
Table 4. Comparison of results of hydrodynamic coefficients under different R e with that of experimental benchmark data adopted from [25] in relative differences.
Table 4. Comparison of results of hydrodynamic coefficients under different R e with that of experimental benchmark data adopted from [25] in relative differences.
Re % C D ( α 0 . 00 ) % C L ( α 4 . 18 ) % C D ( α 4 . 18 ) % C L ( α 8 . 22 ) % C D ( α 8 . 22 ) % C L ( α 13 . 10 ) % C D ( α 13 . 10 ) % C L ( α 15 . 20 ) % C D ( α 15 . 20 )
3.94 × 10 6 32.263.3349.481.6924.03−0.6741.61−1.9939.59
5.97 × 10 6 25.671.9717.531.2944.17−1.8450.87−3.2844.76
6.00 × 10 6 25.572.0117.361.3243.61−1.6449.30−2.9142.49
8.86 × 10 6 22.101.7320.070.0926.34−2.3251.85−3.3349.86
1.19 × 10 7 38.37−0.4717.540.0025.81−1.5741.30−3.2432.23
1.89 × 10 7 47.951.135.67−0.889.18−1.5125.84−1.9919.82
Table 5. Comparison of results of C-mesh based on various domain sizes with that of experimental benchmark data adopted from [25] in relative differences.
Table 5. Comparison of results of C-mesh based on various domain sizes with that of experimental benchmark data adopted from [25] in relative differences.
C in / C out NOC % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C L ( AVE ) % C D ( AVE )
7.5 c / 15 c 160,293
203,144−1.7568.46−2.5755.192.16 61.82
251,065−1.4665.74−21.04
423,225−2.1071.32−3.2859.882.69 65.60
524,605−1.6567.15−2.3853.522.01 60.33
15 c / 30 c 160,293
203,144−2.6043.80
251,065−1.3150.88
423,225−1.9656.69
524,605−1.5052.33−2.4042.501.95 47.41
30 c / 60 c 160,293
203,144−1.4547.99−2.4939.871.97 43.93
251,065−1.1944.95−1.8734.191.53 39.57
423,225−1.8350.77−3.2744.6642.55 47.71
524,605−1.3746.32−2.3337.851.85 42.08
60 c / 120 c 160,293−0.8341.54−1.2830.841.06 36.19
203,144−1.3545.53−2.4037.901.87 41.72
251,065
423,225−1.7548.12−3.2242.642.49 45.38
524,605−1.2943.63−2.2835.761.78 39.70
120 c / 240 c 160,293−0.9242.50
203,144−1.3344.90−2.3737.431.85 41.16
251,065−1.0641.43−1.8532.641.45 37.03
423,225−1.7146.89
524,605−1.2442.38−2.2534.791.75 38.59
Table 6. Uncertainty estimations for 4 cases ( 30 c / 60 c , NOC = 203,144, 251,065, 423,225, 524,605, α = 15 . 2 ) based on mesh dependency.
Table 6. Uncertainty estimations for 4 cases ( 30 c / 60 c , NOC = 203,144, 251,065, 423,225, 524,605, α = 15 . 2 ) based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 2.0010.022.492.000.2139.872.000.122.13
ϕ 2 16.181.870.5434.190.120.81
ϕ 3 17.733.270.6244.660.216.39
ϕ 4 11.572.330.3737.850.091.91
Table 7. Comparison of results of C-mesh based on various chord-wise spacings with that of experimental benchmark data adopted from [25] in relative differences with 670 nodes in chord-wise direction for reference setting.
Table 7. Comparison of results of C-mesh based on various chord-wise spacings with that of experimental benchmark data adopted from [25] in relative differences with 670 nodes in chord-wise direction for reference setting.
log R G N c i 0 NOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
−2169,8601.8518.10−60.521.0247.98−22.60
−1203,5341.9817.39−62.911.3743.79−30.17−1.2846.156.94−2.1036.680.521.6836.0025.13
0251,0651.9917.36−62.481.4242.89−30.23
1318,4131.9817.47−61.991.2844.17−26.44−1.8150.8913.49
2413,4752.0117.36−62.581.3243.61−27.09−1.6449.3011.47−2.9142.494.261.9738.1926.35
3548,0222.0317.15−63.171.4442.09−29.86−1.0543.734.52−1.7634.33−1.741.5734.3224.82
4738,2952.0317.15−62.901.4242.16−28.96−1.1244.265.42−1.8835.22−1.121.6134.7024.60
51,007,0912.0217.15−62.631.3942.23−28.06−1.1844.796.31−1.5634.66−0.841.5434.7124.46
Table 8. Uncertainty estimations for C-mesh ( log R G N c i 0 = 2 , 3 , 4 , 5 , α = 15 . 2 ) with varying chord-wise spacings based on mesh dependency.
Table 8. Uncertainty estimations for C-mesh ( log R G N c i 0 = 2 , 3 , 4 , 5 , α = 15 . 2 ) with varying chord-wise spacings based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 2.002.022.911.041.2842.491.230.704.26
ϕ 2 0.461.761.1134.320.591.74
ϕ 3 0.061.880.7235.220.371.11
ϕ 4 0.041.560.6734.650.340.84
Table 9. Comparison of results of H-mesh based on various chord-wise spacings with that of experimental benchmark data adopted from [25] in relative differences with 386 nodes in chord-wise direction for reference setting.
Table 9. Comparison of results of H-mesh based on various chord-wise spacings with that of experimental benchmark data adopted from [25] in relative differences with 386 nodes in chord-wise direction for reference setting.
log R G N c i 0 NOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
−2178,2092.5318.68−83.921.0861.15−34.30−3.4784.9923.30−6.2881.6113.113.3461.6138.66
−1208,2452.6117.19−83.641.4850.20−39.06−2.5564.6619.18−4.6060.3310.912.8148.0938.20
0255,8522.7317.41−85.961.8245.53−45.29−1.4551.548.66
1333,3062.8417.90−88.341.9743.93−48.35−1.0547.094.81−2.0937.972.961.9936.7236.11
2462,2882.8417.75−87.832.0542.18−49.32−0.6741.651.45−1.3030.54−0.371.7233.0334.74
3683,8622.8717.24−87.772.0841.06−49.42−0.6440.211.63−1.2829.36−0.091.7231.9734.73
41,079,7602.8917.06−88.122.1340.32−50.54−0.4237.89−1.11−0.8526.08−2.601.5730.3435.59
Table 10. Uncertainty estimations for H-mesh ( log R G N c i 0 = 1 , 2 , 3 , 4 , α = 15 . 2 ) with varying chord-wise spacings based on mesh dependency.
Table 10. Uncertainty estimations for H-mesh ( log R G N c i 0 = 1 , 2 , 3 , 4 , α = 15 . 2 ) with varying chord-wise spacings based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 2.002.202.081.120.4237.960.950.212.96
ϕ 2 0.801.300.3130.540.160.37
ϕ 3 0.071.280.2129.350.310.09
ϕ 4 0.08 0.850.1426.080.102.60
Table 11. Comparison of results of O-mesh based on various chord-wise spacings to that of experimental benchmark data adopted from [25] in relative differences with 296 nodes in chord-wise direction for reference setting.
Table 11. Comparison of results of O-mesh based on various chord-wise spacings to that of experimental benchmark data adopted from [25] in relative differences with 296 nodes in chord-wise direction for reference setting.
log R G N c i 0 NOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
−2116,0127.66−0.45−87.3025.0773.3816.88−54.321016.56−430.26−76.12597.97−111.5640.79422.09161.50
−1163,8523.4815.45−77.2922.23100.942.36−2.7870.684.52−48.31540.63−67.6819.20181.9337.96
0193,752−0.3928.12−82.6120.05132.553.36−3.3584.925.53−3.9966.75−3.116.9578.0823.65
1331,292−1.2434.48−96.7719.16142.69−8.57−3.80106.31−6.21−5.0096.24−9.147.3094.9330.17
2472,800−4.0044.66−87.6417.65163.98−11.24−3.90110.78−11.54−4.5290.11−11.437.52102.3830.46
3666,77017.96156.92−6.51−4.74124.83−16.92−52.91842.17−653.75
4944,840−7.2248.92−111.48−8.0391.7317.80
Table 12. Comparison of results of Hybrid-mesh based on various chord-wise spacings with that of experimental benchmark data adopted from [25] in relative differences.
Table 12. Comparison of results of Hybrid-mesh based on various chord-wise spacings with that of experimental benchmark data adopted from [25] in relative differences.
c NOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
2.23 × 10 3 165,8013.8926.16−139.772.9358.92−102.81−0.4566.54−27.97−1.8856.56−17.502.2952.0572.01
1.12 × 10 3 220,2302.8124.61−98.771.9756.27−65.66−0.2857.21−19.72−0.3943.43−20.391.3645.3851.13
7.88 × 10 4 265,5463.1523.68−111.242.3353.73−77.53−0.2655.59−20.22−1.0342.59−12.241.6943.9055.31
5.88 × 10 4 314,0703.1123.18−108.702.2352.39−71.74−0.6055.97−13.90−0.8842.25−13.511.7143.4551.96
4.46 × 10 4 374,8532.9323.23−102.732.0953.57−66.59−0.6455.73−12.36−1.0241.83−12.421.6743.5948.53
3.15 × 10 4 476,6622.8922.13−100.242.1452.34−66.90−0.6054.14−11.16−1.3542.18−8.661.7442.7046.74
2.23 × 10 4 620,9952.7322.04−94.711.8451.30−57.25−0.9854.76−4.10−1.3841.61−6.551.7342.4240.66
Table 13. Uncertainty estimations for Hybrid-mesh ( c = 5.88 × 10 4 c , 4.46 × 10 4 c , 3.15 × 10 4 c , 2.23 × 10 4 c , α = 15 . 2 ) with varying chord-wise spacings based on mesh dependency.
Table 13. Uncertainty estimations for Hybrid-mesh ( c = 5.88 × 10 4 c , 4.46 × 10 4 c , 3.15 × 10 4 c , 2.23 × 10 4 c , α = 15 . 2 ) with varying chord-wise spacings based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 2.008.080.881.210.1342.242.000.1513.50
ϕ 2 6.581.020.1241.830.0612.42
ϕ 3 5.581.350.1242.180.038.66
ϕ 4 5.961.380.0841.600.0046.54
Table 14. Comparison of results of C-mesh based on various layer-wise growth rates with that of experimental benchmark data adopted from [25] in relative differences.
Table 14. Comparison of results of C-mesh based on various layer-wise growth rates with that of experimental benchmark data adopted from [25] in relative differences.
rNOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
1.20548,0221.9617.06−63.801.3143.17−28.81−1.1746.094.41−1.7935.84−2.751.5635.5424.94
1.15548,022
1.10548,0222.0317.17−63.061.4242.31−29.35−1.1744.875.98−2.0135.40−0.401.6634.9424.70
1.05548,0222.0317.16−62.981.4242.31−29.26−1.1844.886.15−2.0235.42−0.301.6634.9424.67
1.03548,0222.0317.16−63.191.4342.29−29.67−1.1544.805.60−1.9535.21−0.811.6434.8624.82
Table 15. Comparison of results of H-mesh based on various layer-wise growth rates that of experimental benchmark data adopted from [25] in relative differences.
Table 15. Comparison of results of H-mesh based on various layer-wise growth rates that of experimental benchmark data adopted from [25] in relative differences.
rNOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
1.2046,2882.8616.54−87.452.0541.22−48.90−0.7241.522.02−1.3830.81−0.101.7532.5234.62
1.1546,2882.8717.12−87.992.0841.55−49.74−0.6141.140.50−1.1729.96−1.481.6832.4434.93
1.1046,2882.8417.75−87.822.0542.20−49.30−0.6741.661.46−1.3030.55−0.371.7233.0434.74
1.0546,2882.8617.52−87.962.0741.74−49.57−0.6341.040.94−1.2329.92−0.781.7032.5534.81
1.0346,2882.8717.40−88.052.0841.60−49.77−0.6040.870.51−1.1829.74−1.181.6832.4034.88
Table 16. Comparison of results of Hybrid-mesh based on various layer-wise growth rates with that of experimental benchmark data adopted from [25] in relative differences.
Table 16. Comparison of results of Hybrid-mesh based on various layer-wise growth rates with that of experimental benchmark data adopted from [25] in relative differences.
rNOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
1.20607,9872.3021.80−81.671.2350.70−36.72−1.6454.727.37−2.1444.18−2.791.8342.8532.14
1.15611,4912.2422.13−80.521.3651.11−41.85−1.1654.54−1.17−1.4341.82−6.331.5542.4032.47
1.10620,9952.7322.03−94.781.8451.32−57.42−0.9854.76−4.09−1.8743.50−4.281.8642.9040.14
1.05652,2332.8717.74−98.811.7848.03−54.62−1.2454.52−0.37−1.6941.02−5.831.8940.3339.91
1.03660,2813.376.31−109.962.1638.29−63.27−1.0050.00−3.25
Table 17. Comparison of results of C-mesh based on various radial nodes with that of experimental benchmark data adopted from [25] in relative differences with 150 nodes in radial directions for reference setting.
Table 17. Comparison of results of C-mesh based on various radial nodes with that of experimental benchmark data adopted from [25] in relative differences with 150 nodes in radial directions for reference setting.
log R G N r i 0 NOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
−2272,1721.8917.08−63.181.1646.20−19.87−1.2747.335.08−1.9037.84−2.641.5537.1122.69
−1386,1902.0217.17−64.081.4142.79−30.40−1.0945.273.84−1.7535.82−2.591.5735.2625.23
0548,0222.0317.17−63.061.4242.31−29.35−1.1744.865.98−2.0136.06−0.401.6635.1024.70
1776,0582.0317.17−63.061.4242.31−29.35−1.1744.865.98−2.0136.06−0.401.6635.1024.70
21,099,722
31,555,7942.0417.12−63.041.4342.15−29.31−1.1744.646.04−2.0035.88−0.361.6634.9524.69
42,203,1222.0417.10−63.021.4442.08−29.29−1.1644.546.06−2.0035.80−0.341.6634.8824.68
Table 18. Comparison of results of H-mesh based on various radial nodes to experimental benchmark data adopted from [25] in relative differences with 150 nodes in radial directions for reference setting.
Table 18. Comparison of results of H-mesh based on various radial nodes to experimental benchmark data adopted from [25] in relative differences with 150 nodes in radial directions for reference setting.
log R G N r i 0 NOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
−2314,0882.8417.74−87.762.0542.20−49.31−0.6741.661.45−1.3030.52−0.361.7233.0334.72
−1375,3442.8417.75−87.822.0542.20−49.30−0.6741.651.45−1.3030.55−0.371.7233.0434.74
0462,2882.8417.75−87.832.0542.20−49.30−0.6741.651.45−1.3030.55−0.381.7233.0434.74
1584,8002.8417.75−87.832.0542.20−49.31−0.6741.651.45−1.3030.55−0.381.7233.0434.74
2758,6882.0542.19−49.30−0.6741.651.45−1.3030.54−0.38
31,002,7122.0542.19−49.31−0.6741.641.45−1.3030.53−0.37
41,351,488−0.6741.641.45−1.3030.54−0.37
Table 19. Comparison of results of C-mesh based on various wake nodes with that of experimental benchmark data adopted from [25] in relative differences with 300 nodes in wake directions for reference setting.
Table 19. Comparison of results of C-mesh based on various wake nodes with that of experimental benchmark data adopted from [25] in relative differences with 300 nodes in wake directions for reference setting.
log R G N r i 0 NOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
−2503,3223.5432.94−41.901.4242.31−29.37−1.1744.865.97−2.0136.07−0.412.0339.0419.41
−1521,7982.0317.17−63.071.4242.31−29.37−1.1744.865.97−2.0136.07−0.411.6635.1024.70
0548,0222.0317.17−63.071.4242.31−29.36−1.1744.865.97−2.0136.06−0.411.6635.1024.70
1584,9742.0317.17−63.070.3941.54−4.57−1.1744.865.98−2.0136.06−0.401.4034.9118.51
2637,4222.0317.16−63.071.4242.30−29.35−1.1744.845.98
3711,6242.0317.16−63.071.4242.29−29.35−1.1844.836.01−2.0136.04−0.381.6635.0824.70
4816,222
Table 20. Comparison of the results of H-mesh based on various wake nodes with that of experimental benchmark data adopted from [25] in relative differences with 300 nodes in wake directions for reference setting.
Table 20. Comparison of the results of H-mesh based on various wake nodes with that of experimental benchmark data adopted from [25] in relative differences with 300 nodes in wake directions for reference setting.
log R G N r i 0 NOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
−2405,7382.8417.75−87.832.0741.74−49.57−0.6741.651.45−1.3030.55−0.381.7232.9234.81
−1429,1122.8417.75−87.782.0741.74−49.57−0.6741.651.45−1.3130.57−0.351.7232.9334.79
0462,2882.8417.75−87.802.0741.74−49.57−0.6741.651.45−1.3130.57−0.351.7232.9334.79
1509,0362.8417.74−87.792.0741.74−49.57−0.6741.651.45−1.3130.56−0.351.7232.9234.79
2575,3882.8417.75−87.802.0741.74−49.57−0.6741.651.45−1.3130.57−0.351.7232.9334.79
3669,261−2296.892.0741.74−49.57−0.6741.631.46−1.3130.54−0.33
4801,588−8910.582.0741.74−49.57−0.6741.631.47−1.3130.54−0.33
Table 21. Comparison of the results of Hybrid-mesh based on various element sizes with that of experimental benchmark data adopted from [25] in relative differences.
Table 21. Comparison of the results of Hybrid-mesh based on various element sizes with that of experimental benchmark data adopted from [25] in relative differences.
Element SizeNOC % C L ( 4 . 18 ) % C D ( 4 . 18 ) % C M ( 4 . 18 ) % C L ( 8 . 22 ) % C D ( 8 . 22 ) % C M ( 8 . 22 ) % C L ( 13 . 10 ) % C D ( 13 . 10 ) % C M ( 13 . 10 ) % C L ( 15 . 20 ) % C D ( 15 . 20 ) % C M ( 15 . 20 ) % C L ( AVE ) % C D ( AVE ) % C M ( AVE )
0.17678c621,1172.7122.02−94.431.8451.32−57.41−0.9954.35−3.76−1.4141.77−6.381.7442.3740.49
0.12500c719,0432.5921.88−90.641.8051.40−55.08−0.7253.35−6.48−1.1040.96−8.891.5541.9040.27
0.08839c919,2552.7721.42−96.732.0251.73−61.37−0.7654.01−6.93−1.2440.11−7.051.7041.8243.02
0.06250c1,314,7632.7021.71−93.171.8851.95−57.74−0.8654.26−5.91−1.1642.14−9.811.6542.5141.66
0.04419c2,106,0772.7020.96−91.301.9250.73−57.52−0.5852.55−9.67−1.4141.24−7.211.6541.3741.43
Table 22. Uncertainty estimations for C-mesh with varying radial node distributions ( log R G N c i 0 = 2 , 1 , 0 , 1 ) based on mesh dependency.
Table 22. Uncertainty estimations for C-mesh with varying radial node distributions ( log R G N c i 0 = 2 , 1 , 0 , 1 ) based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 0.232.871.902.000.0637.840.560.332.63
ϕ 2 2.691.750.0435.810.232.58
ϕ 3 1.772.000.0236.060.180.40
ϕ 4 1.272.000.0236.060.160.40
Table 23. Uncertainty estimations for H-mesh with varying radial node distributions ( log R G N c i 0 = 2 , 1 , 0 , 1 ) based on mesh dependency.
Table 23. Uncertainty estimations for H-mesh with varying radial node distributions ( log R G N c i 0 = 2 , 1 , 0 , 1 ) based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 2.000.00021.302.000.000330.562.000.00010.38
ϕ 2 0.00061.300.000230.550.00010.37
ϕ 3 0.00041.300.000130.550.00010.37
ϕ 4 0.00071.300.000130.550.00010.37
Table 24. Uncertainty estimations for C-mesh with varying wake node distributions ( log R G N c i 0 = 2 , 1 , 0 , 1 ) based on mesh dependency.
Table 24. Uncertainty estimations for C-mesh with varying wake node distributions ( log R G N c i 0 = 2 , 1 , 0 , 1 ) based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 2.000.0032.002.000.00136.062.000.0010.42
ϕ 2 0.0022.000.00136.060.0010.42
ϕ 3 0.0012.000.00136.060.0010.42
ϕ 4 0.0012.000.00136.060.0010.42
Table 25. Uncertainty estimations for H-mesh with varying wake node distributions ( log R G N c i 0 = 2 , 1 , 0 , 1 ) based on mesh dependency.
Table 25. Uncertainty estimations for H-mesh with varying wake node distributions ( log R G N c i 0 = 2 , 1 , 0 , 1 ) based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 2.000.01021.302.000.000230.552.000.00060.38
ϕ 2 0.00571.310.000630.570.00020.35
ϕ 3 0.00321.310.000630.570.00010.35
ϕ 4 0.00171.310.000730.560.00010.35
Table 26. Uncertainty estimations for Hybrid-mesh (Element size = 0.17678 c , 0.12500 c , 0.08839 c , 0.06250 c , α = 15 . 2 ) with varying far-field node distributions based on mesh dependency.
Table 26. Uncertainty estimations for Hybrid-mesh (Element size = 0.17678 c , 0.12500 c , 0.08839 c , 0.06250 c , α = 15 . 2 ) with varying far-field node distributions based on mesh dependency.
Solution C L C D C M
p % U G % | E | p % U G % | E | p % U G % | E |
ϕ 1 2.000.141.412.000.1641.770.860.636.38
ϕ 2 0.171.100.1540.960.618.88
ϕ 3 0.141.240.2440.110.577.04
ϕ 4 0.121.160.2042.140.439.81
Table 27. Recommended mesh properties for different mesh types.
Table 27. Recommended mesh properties for different mesh types.
Mesh TypesC-MeshH-MeshO-MeshHybrid-Mesh
Mesh Properties
R e 6.00 × 10 6 6.00 × 10 6 6.00 × 10 6 6.00 × 10 6
Domain size 30 c / 60 c 30 c / 60 c 60 c / 60 c 30 c / 60 c
Chord-wise spacing1895 nodes692 nodes658 nodes 1.78 × 10 3 c
Layer-wise growth rate1.101.051.10
Far-field node distributionRadial-150 nodesRadial-150 nodesRadial-∖Element size- 0.125 c
Wake-300 nodesWake-300 nodesWake-∖
Number of cells 5.48 × 10 5 4.62 × 10 5 7.19 × 10 5
Table 28. Comparison of mesh types with recommended settings to that of other CFD codes adopted from [19] in relative differences.
Table 28. Comparison of mesh types with recommended settings to that of other CFD codes adopted from [19] in relative differences.
Codes α ( ) C L C D
C-MeshH-MeshHybrid-MeshC-MeshH-MeshHybrid-Mesh
CFL3D0\\\0.91%1.44%3.47%
10−0.57%−0.01%0.05%7.88%7.18%16.02%
15−1.37%−0.72%−0.63%11.05%6.80%15.94%
FUN3D0\\\1.03%1.57%3.60%
10−1.13%−0.58%−0.53%6.42%5.72%14.45%
15−1.64%−0.99%−0.90%8.31%4.17%13.09%
NTS0\\\0.91%1.44%3.47%
10−0.45%0.11%0.17%6.59%5.89%14.63%
15−1.58%−0.93%−0.84%12.67%8.36%17.64%
Table 29. More profiles studied.
Table 29. More profiles studied.
Profile Re Reference
NACA 0006 6.00 × 10 6 Abbott and Doenhoff [39]
NACA 0009 6.00 × 10 6 Abbott and Doenhoff [39]
NACA 0015 3.26 × 10 6 Jacobs and Sherman [40]
NACA 0018 3.26 × 10 6 Jacobs and Sherman [40]
Table 30. Recommended mesh properties for different objective.
Table 30. Recommended mesh properties for different objective.
ObjectivePropertiesSettings
Accurate prediction for both C L and C D Mesh typeH-mesh
R e Actual Reynolds number whenever possible
Domain sizeAs large as possible or at least 30 c around and 60 c after the profile
Chord-wise spacingMore than 692 nodes for c = 1 m
Layer-wise growth rate1.05
Number of radial nodes150 nodes
Number of wake nodes300 nodes
Number of mesh cellsNo less than 4.62 × 10 5
Accurate C L only with less mesh numberMesh typeC-mesh
R e 6 × 10 6
Domain size 30 c around and 60 c after the profile
Chord-wise spacing670 nodes for c = 1 m
Layer-wise growth rate1.10
Number of radial nodes150 nodes
Number of wake nodes300 nodes
Number of mesh cellsAbout 2.51 × 10 5
Acceptable C L and C D with simple meshing strategyMesh typeHybrid-mesh
R e 6 × 10 6
Domain size 30 c around and 60 c after the profile
Chord-wise spacing 2.23 × 10 4 c
Layer-wise growth rate1.10
Element size 0.125 c
Number of mesh cellsAbout 7.19 × 10 5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, S.; Liu, J.; Hekkenberg, R. Mesh Properties for RANS Simulations of Airfoil-Shaped Profiles: A Case Study of Rudder Hydrodynamics. J. Mar. Sci. Eng. 2021, 9, 1062. https://doi.org/10.3390/jmse9101062

AMA Style

Lu S, Liu J, Hekkenberg R. Mesh Properties for RANS Simulations of Airfoil-Shaped Profiles: A Case Study of Rudder Hydrodynamics. Journal of Marine Science and Engineering. 2021; 9(10):1062. https://doi.org/10.3390/jmse9101062

Chicago/Turabian Style

Lu, Suli, Jialun Liu, and Robert Hekkenberg. 2021. "Mesh Properties for RANS Simulations of Airfoil-Shaped Profiles: A Case Study of Rudder Hydrodynamics" Journal of Marine Science and Engineering 9, no. 10: 1062. https://doi.org/10.3390/jmse9101062

APA Style

Lu, S., Liu, J., & Hekkenberg, R. (2021). Mesh Properties for RANS Simulations of Airfoil-Shaped Profiles: A Case Study of Rudder Hydrodynamics. Journal of Marine Science and Engineering, 9(10), 1062. https://doi.org/10.3390/jmse9101062

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop