Next Article in Journal
How Does Public Opinion Influence Production Safety within Small and Medium Enterprises in the Sustainability Context?
Next Article in Special Issue
Design and Test for a New Type of Permanent Magnet Synchronous Generator Applied in Tidal Current Energy System
Previous Article in Journal
Wine Quality Control Parameters and Effects of Regional Climate Variation on Sustainable Production
Previous Article in Special Issue
Finite Element Analysis of Different Infill Patterns for 3D Printed Tidal Turbine Blade
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Evaluation of the Immersed-Solid Method for Simulating the Flow Field around Hydrofoil

1
College of Water Resources and Civil Engineering, China Agricultural University, Beijing 100083, China
2
College of Engineering, China Agricultural University, Beijing 100083, China
3
Beijing Engineering Research Center of Safety and Energy Saving Technology for Water Supply Network System, China Agricultural University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(4), 3493; https://doi.org/10.3390/su15043493
Submission received: 9 January 2023 / Revised: 10 February 2023 / Accepted: 13 February 2023 / Published: 14 February 2023
(This article belongs to the Special Issue Ocean and Hydropower)

Abstract

:
The wall boundary is important in computational-fluid-dynamics simulations. If extremely small leakage, changing leakage or a moving body exists in the simulation case, the difficulty in meshing and solving near-wall flow increases. The immersed-solid method, which inserts a rigid, solid body into the entire fluid domain, was a choice to solve the wall-boundary-solution problems mentioned above, without considering mesh deformation. The purpose of this paper is to verify the effectiveness of the immersed-solid method in the simulation of extremely small leakage, changing leakage or a moving body, and to provide a theoretical basis for the use of the submerged-solid method in engineering. In this study, the NACA0015 hydrofoil was used to check the hydrodynamic characteristics in using the immersed-solid method. The comparative study was conducted at the incidence angle of 8 degrees and a Reynolds number of 5.0 × 105, by using the immersed-solid and traditional no-slip-wall boundary. The results show that the flow striking and separation with pressure rise and drop can be correctly captured using an immersed-solid setup with boundary tracking. However, the accuracy of pressure and velocity field using the immersed-solid method was insufficient. The turbulence-kinetic energy was much higher around the immersed-solid foil body. Generally, the immersed-solid method can qualitatively predict the correct hydrodynamic characteristics. Its convergence ability is better, and it can save approximately 20% of CPU time, even if the grid density is 4.39 times of the traditional no-slip wall. Therefore, the immersed-solid method can be a good choice for engineering-flow cases with complex wall problems.

1. Introduction

Computational fluid dynamics (CFD) is widely used in the performance evaluation [1] and optimization [2] of hydraulic machinery. The wall boundary is the most common boundary type in computational-fluid-dynamics (CFD) simulations, especially considering fluid flows around solid bodies [3,4]. Currently, an accurate solution of near-wall flow is no more difficult. The direct solving of the usage of different wall-function estimations and mixed formulations can be applicable to wall function insensitive to y+ [5], pipe flow [6] and high Reynolds-number flow [7]. However, some near-wall puzzles still exist to be solved in modern-engineering CFD simulations. For example, extremely small leakage is difficult to simulate when it is much smaller than the other regions [8]. Putting mesh topology inside extremely small leakage is a challenge, especially without compromising the accuracy with the simulation cost. Another example is variable leakage, such as leakage in piston-pump housing and gate housing. Blanco et al. [9] used dynamic-grid technology to simulate the flow of a piston pump, and Akoz et al. [10] compared eight different grid systems to obtain a grid that can accurately simulate the flow of a sluice gate. In this case, researchers need to change both the grid-topology distribution and the density. This will strongly affect the time cost of CFD’s pre-process and the solution process. Moving a solid body in fluid is also a difficult problem [11]. If the solid body only translates, the mesh-re-building process can be simplified and well solved. If translation and rotation happen together [12,13], the mesh re-building becomes difficult, and even be unable to be made automatic. The immersed-solid method has become another wall-boundary choice in the commercial software ANSYS CFX, putting a solid domain into the entire fluid domain [14]. It provides a new way for the near-wall flow simulations to solve the puzzles above. Yoon et al. [15] simplified the calculation of small-clearance leakage flow in an external gear pump by the immersed-solid method, overcoming numerical difficulties and saving calculation time. The simulation results show that there is little difference between the numerical predictions and the experimental data. Zhang et al. [16,17] applied the immersed-solid method to the simulation of moving-boundary problems which were difficult to solve by other methods, and simulated the flow field inside the double-helical pump. At present, the immersed-solid method has been widely used in the simulation of extremely small leakage, changing leakage, and the moving body, but the gap between its accuracy and the traditional wall-boundary simulation is not clear, and a good mesh-matching scheme on the fluid–solid interface has not been found.
Hydrofoil, which is a typical, simplified solid body, is chosen as the studied object in this case, due to its standardized profile. It has been widely investigated in both numerical simulations and experiments [18,19,20,21]. For typical hydrofoils, the lift/drag forces, flow-striking and separation point, surface-pressure distribution, cavitation performance and structural-response characteristics are already clear. Therefore, the hydrodynamic characteristics can be well studied for comparing the immersed-solid method with the traditional wall-boundary type [22]. In this study, there are three main points to be solved. Firstly, it is necessary to find a way to mesh the fluid domain, solid domain and the fluid–solid interface when using the immersed-solid method. The mesh independency in using the immersed-solid method and using the traditional wall boundary are checked and compared. Secondly, the flow field should be examined in detail, including the pressure field, velocity field and turbulence field. This will be helpful for the application of the immersed-solid method in engineering flow cases. Thirdly, the convergence ability should be comparatively studied in the immersed-solid case and in the traditional wall-boundary case. The time cost, hardware cost and CFD pre-processing cost need a comprehensive comparison to provide guidance for engineering cases.

2. Numerical Methods

2.1. Turbulence Modeling

The numerical study adopts the Reynolds-averaged Navier–Stokes (RANS) equation. The RANS equation is the equation of average turbulent motion obtained by averaging the physical quantity of turbulence with the sum of an average value and a pulsation value. The eddy-viscosity turbulence model is used to close the averaged equations. It follows the Boussinesq eddy-viscosity hypothesis [23], which establishes the relationship between the eddy viscosity, μt, and the Reynolds stress, τij. In this study, the SST turbulence model [24,25,26], which deals well with the strong freestream and improves the accuracy of adverse pressure gradients by combining the standard k-ε model and Wilcox’s k-ω model, was used. The SST model was based on the turbulence-kinetic-energy, k, equation and the specific-dissipation-rate, ω, equation [27]:
ρ k t + ρ u j k x j = P k β ρ ω k + x j μ + σ k μ t k x j
ρ ω t + ρ u j ω x j = α ω k P k β ρ ω 2 + x j μ + σ ω μ t ω x j + 2 1 F 1 ρ σ ω 2 ω k x j ω x j
where the production term Pk can be expressed by:
P k = τ i j u i x j
τ i j = μ t 2 S i j 2 3 u k x k δ i j 2 3 ρ k δ i j
where δij is the Kronecker number and Sij is the shear-stress tensor, so that
S i j = 1 2 u i x j + u j x i
and
μ t = 1 max a 1 ω , S F 2
F1 is the 1st blending function for variants (for example, φ), by
φ = F 1 φ 1 + 1 - F 1 φ 2
where F1 follows
F 1 = tanh min max k β ω d w , 500 ν d w 2 ω , 4 ρ σ ω 2 k C D k ω d w 2 4
where dw is the distance to the nearest wall, ν is the kinematic viscosity, and CD is
C D k ω = max 2 ρ σ ω 2 1 ω k x j ω x j , 10 20
F2 is the 2nd blending function:
F 2 = tanh max 2 k β ω d , 500 ν d w 2 ω 2
In the above equations, ρ is the density. x, t, and u denote coordinate, time and velocity, respectively. The constants are σk1 = 0.85, σω1 = 0.65, β1 = 0.075 for the k-ω setup and σk2 = 1.00, σω2 = 0.856, β2 = 0.0828 for the k-ε setup. The other constants are β* = 0.09, a1 = 0.31 and α = 0.556.

2.2. Wall Treatments

2.2.1. No-Slip Wall

The no-slip wall boundary has been popularly used in computational-fluid-dynamics (CFD) research [22]. It assumes that fluid molecules may not slip from a wall in an ideal situation [28,29]. Generally, a common approximation for fluid slip can be drawn as:
u u w = β s u n
where u is the fluid velocity, uw is the velocity on the wall, βs is the slip length so that βs = Cs·lmfp where lmfp is the mean free path and Cs is the empirical coefficient [30]. For a no-slip wall-boundary condition, βs is 0, which means u = uw.

2.2.2. Automatic-Wall Functions

If the SST turbulence model was applied, the automatic-wall treatments [6] were also adopted, to empirically model the near-wall flow. As shown in Figure 1, the red line in the figure is the change curve of fluid velocity with the distance from the wall. The near-wall-flow region included 3 main parts: the viscous sublayer, buffer layer and logarithmic layer [31,32]. In the viscous sublayer, there is [33]:
u + = y +
where u+ is the dimensionless velocity of t u+ = u/ut; u is the velocity, ut is the friction velocity of ut = τ w / ρ ., τw is the wall shear stress, and y+ is the dimensionless off-wall distance of y+ = y·ut, where y is the off-wall distance. In the logarithmic layer, there is [33]:
u + = 1 κ ln y + + C +
where κ is the von Karman constant and C+ is the velocity-deviation constant. The automatic-wall functions blend the ω solution and the ut solution by [34]
ω = ω v i s 2 + ω l o g 2
u t = u t v i s 4 + u t l o g 4 4
where superscript/subscript vis represent the viscous sublayer and the buffer layer, respectively, and superscript/subscript log represent the logarithmic layer.
ω v i s = 6 ν β y 2 ; ω log = u t C μ 1 / 4 κ y
u t v i s = u y + ; u t l o g = u log E u y +
where Eu is the model constant of 9.81, and Cμ is the eddy-viscosity coefficient. Here we define y as the mesh-element–centroid off-wall distance in the CFD simulation; if the y value is small, 1/y2 dominates, to emphasize the viscous characters. On the contrary, if the y value is large, 1/y dominates, being mainly logarithmic.

2.2.3. Immersed-Solid Method

The immersed-solid method is enabled in the commercial software ANSYS CFX by immersing a rigid-solid domain into a fluid domain [14]. This immersed-solid domain was more like a region in which a momentum source is added to force the fluid to move with the solid. Figure 2 shows a schematic map of the mesh elements and nodes when using the immersed-solid method. Table 1 lists the acronyms defined for the immersed-solid method.
In ANSYS CFX, there are some important settings for the immersed-solid method, as follows:
  • Wall distance: 0 on the nodes inside the immersed solid; a function of the wall scales on the nodes near the inside of the immersed solid; the physical distance of the nodes near the immersed solid was achieved using the boundary-tracking method for wall functions.
  • k and ω in SST turbulence-model cases: 0 on the nodes inside the immersed solid; scalable-wall treatment was used for ω in the logarithmic region by ωlog = u*2/(Cμ1/2κνy*) where u*= Cμ1/2k1/4 and y* = u*Δy/ν. The value of Δy is the distance between the node near the immersed solid and the immersed-solid boundary. The y* value is lowered to 11.06 for dividing the logarithmic region and the linear viscous region.
  • Friction velocity ut: Based on the scalable-wall treatment, ut in the logarithmic region is utlog = ΔU/[1/κln(y*) + C] = ΔU/u+, where ΔU is the fluid tangential velocity relative to the immersed boundary.
According to the setup above, the velocity on the near-wall nodes can be forced as
U N W F t a n = U I B I S t a n + u / u + F L u / u + N W U F L t a n U I B I S t a n
where superscript F and IS denote the fluid and immersed solid, respectively. The momentum source can be calculated based on the tangential-force velocity (UNWF)tan and the wall-normal velocity (UIBIS)nm:
U N W F = U N W F t a n + U I B I S n m
Hence, fluid can be prevented from getting through the immersed-solid domain.

3. Studied Object

3.1. NACA0015 Profile

The NACA0015 [35] profile belongs to the typical NACA four-digit-profile series. It is symmetric along the meanline. If we define the foil-meanline direction as mf and the thickness direction (the origin from the meanline) as tf, the NACA0015 mf -tf profile can be illustrated in Figure 3, where the total length along mf is denoted as c and the maximum thickness is denoted as tfm. The relationship between tf and mf can be expressed as:
± t f c = t m f 0.2 c 0.297 m f c 0.126 m f c 0.352 m f c 2 + 0.284 m f c 3 0.102 m f c 4
In addition, the leading edge of the NACA0015 profile is an arc, whose radius is rLE, as in the following:
r L E c = 1.102 t m f c 2
Thus, the maximum thickness, tfm, approximately locates at mf/c = 0.3. The incidence angle (between the meanline and rotated meanline) can be denoted as α, the unit representing a degree. It represents the foil-rotation status, based on the rotation centroid.

3.2. Case Detail

In this study, the NACA0015 foil was studied at α = 8 degrees inside an L1 × L2 × L3 cuboid domain, as shown in Figure 4. L1, L2 and L3 were respectively equal to 4.35c, 1.04c and 0.70c. Coordinate components X, Y and Z denote the three orthogonal directions, where X is along the flow direction, Y is perpendicular to the flow direction and Z is along the foil span. Because the foil does not change shape and scale along Z, this 3D case can be treated as an X-Y quasi-2D case, which provides a good case simplification. Fluid flows into the domain at the inlet boundary and out of the domain at the outlet boundary. The four surrounding X-Y and X-Z planes are walls.

3.3. Computational-Fluid-Dynamics Setups

In this case, the hydrodynamic characters were mainly studied. ANSYS CFX was used as the CFD solver by using the unsteady-RANS method with the SST turbulence model, as introduced in the above sections. The fluid medium was set as water at 20 °C under the 1 atm pressure environment. Compressibility, heat transfer and mass transfer were not considered. The advection scheme and the turbulence numeric were set as high resolution. The inlet boundary was set as velocity inlet with uniform Vin distribution, and the pressure at the inlet was at zero gradient. The outlet boundary was set as the static-pressure boundary of relatively 0 Pa, and the velocity at the outlet was at zero gradient. The surrounded walls were set as no-slip type. The foil surface was compared, between the no-slip type and the immersed-solid type. In this case, the Reynolds number, Re, was 5.0 × 105, according to the following equation:
R e = ρ c V i n μ
The unsteady simulation was based on the steady simulation of 600 iterations. The unsteady simulation covered a total of 0.1 s and spent 0.0001 s on each time step. The maximum iteration number for each time step was 10, with the convergence criteria of RMS 1.0 × 10−5.

4. Numerical-Accuracy Check and Validation

4.1. No-Slip Wall Case

In the no-slip-wall case (the foil surface was set as the no-slip wall), the domain was meshed by structural hexahedral elements, using ICEMCFD. Six different mesh schemes, whose node number increased proportionally and gradually, as listed in Table 2, were comparatively checked for accuracy by comparing the numerically predicted pressure with the experimentally tested pressure [15] on the foil surface along the X direction. In this case, the dimensionless pressure coefficient, Cp, was used:
C p = p p i n 0.5 ρ V i n 2
where p is the pressure, and pin is the pressure at the inlet boundary. The position along the X direction is also in a relative form, as X/c.
Figure 5 shows the pressure on the foil surface under different mesh schemes by treating the foil surface as no-slip wall boundary. In addition, the near-wall region where y+ was 1.45~16.22 was predicted by automatic-wall functions. The Nos. one to six mesh schemes predicted the pressure well from X/c ≈ 0.4 to X/c ≈ 1.0 by comparing it with the experimental data [15]. However, the main difference occurred in the front 20% X/c region, especially on the foil-suction side. The minimum Cp value for the coarse mesh was larger than −3. Then, it decreased to −3~−3.25 in the simulations for coarse, mid-coarse and mid-fine meshes. When the mesh nodes increased to 257488 (fine), the numerically predicted Cp curve overlapped well with the experimental data. This situation occurred when the mesh nodes increased to 352594 (very fine). Usually, we called this “mesh independence”. Thus, the mesh scheme No. five with 257488 modes can be selected in this study, with good accuracy.

4.2. Immersed-Solid Case

The foil surface was also treated as the immersed-solid boundary. In this situation, four different mesh combinations were discussed, as listed in Table 3. The pressure distribution on the foil surface was also analyzed along the X direction. The dimensionless Cp and relative X/c were still used, as shown in Figure 6.
Figure 6a shows the pressure on the foil surface under different mesh-scheme combinations. The pulsations can be observed on the immersed-solid predicted Cp curves. These pulsations could be caused by the data interpolation on the fluid–solid interface. Thus, the smooth procedures [36] were conducted to create Figure 6b. The No. 1 combination had the same density as the selected No. 5 mesh in the no-slip-wall case. However, it produced a completely wrong Cp curve, compared with the no-slip-wall case and the experimental data. Therefore, the mesh density was refined to No. 2, but no improvement could be observed. Based on the two results, the boundary tracking was set to near-wall approximations on the immersed solid on the fluid side. The “Coarse-B” curve still produced wrong data. On the other hand, the “Fine-B” curve showed an obvious improvement, in that both the leading-edge pressure rise and drop were correctly captured in position. The amplitude of the maximum and minimum Cp values became closer to the experimental data. Therefore, the other flow characteristics will be analyzed in the next section.

5. Comparative-Flow Analysis

5.1. Forces and Torques

The lift force, FY, (along the Y direction) and the drag force, FX, (along the X direction) are analyzed in Table 4. In the no-slip-wall case, the lift force was 113.01 N, which was larger than that in the immersed-solid case, of 90.82 N. The drag force was 4.95 N in the no-slip-wall case, which was smaller than that in the immersed-solid case, of 8.67 N. As a result, the immersed-solid method obtained a smaller lift/drag ratio compared with the no-slip-wall setup, with its lift/drag ratio being approximately 40% of the no-slip-wall case. The torques on the foil on the base of the rotation centroid, TZ, were also predicted, as is also shown in Table 4. In the no-slip-wall case, the foil torque was 3.29 N·m. In the immersed-solid case, the foil torque was 2.31 N·m, which was approximately 70% of the no-slip-wall case. Under the two wall conditions, the difference between the force and torque obtained by numerical simulation is very obvious. In order to understand the reasons for the differences, we need to analyze the internal flow under different circumstances.

5.2. Pressure-Coefficient Cp Field

Figure 7 compares the Cp contour on the X-Y plane at middle Z position in the domain. The contours in the no-slip case and immersed-solid case were similar. Generally, the flow striking on the foil leading edge was correctly captured in both the two cases. The flow separation after the leading edge was also simulated. Therefore, the immersion-solid method can better simulate the pressure distribution in the flow. However, some differences can be found. Firstly, the very low Cp region (Cp < −2.5) on the foil upper side was wider in the no-slip-wall case than in the immersed-solid case. Secondly, Cp increased rapidly in the upper flow field in the no-slip-wall case. However, the Cp increase was relatively slow in the immersed-solid case. Thirdly, there was a high-pressure region in the lower flow field in the no-slip wall case. In the immersed-solid case, the pressure there was relatively low. In the case of the no-slip wall, the pressure difference between the upper and lower sides of the foil was large, which resulted in the lift of the foil being greater than in the case of the immersed solid. Fourthly, in the no-slip-wall case, except for the very low Cp region, the pressure coefficient in other regions was higher than that in the immersed-solid case.

5.3. Velocity-Coefficient Cv Field

Figure 8 compares the velocity field on the X-Y plane at middle Z position in the domain. The velocity was relative to the domain-inlet velocity as a velocity coefficient Cv = V/Vin. The contours in the no-slip case and immersed-solid case were qualitatively similar but quantitatively different. The leading-edge high-velocity region was captured in the two cases. However, the very high Cv region (Cv > 1.5) was wide in the no-slip-wall case but small in the immersed-solid case, and its distribution was basically consistent with the distribution of the very low Cp region. The low velocity after the trailing edge (the wake) was narrow in the no-slip-wall case, but very wide in the immersed-solid case. A high Cv region was in the lower flow field, as indicated in Figure 8b. This shows that the immersed-method underestimates the flow velocity at the leading and trailing edges of the foil, and overestimates the flow velocity of the upper and lower flow fields of the foil.

5.4. Turbulence-Kinetic-Energy, k, Field

Figure 9 shows the turbulence kinetic-energy, k, contour on the X-Y plane at middle Z position. The legend ranged from 0 to 3.9 m2/s2. It can be seen from the figure that the turbulence-kinetic-energy, k, contours in the two cases are quite different. In the no-slip-wall case, there is a relatively high k region from the foil trailing edge to downstream. It was narrow and relatively weaker than in the immersed-solid case. On the contrary, in the immersed-solid case, the high k region was very wide, covering the entire foil and extending to the downstream far field. The highest k (up to 3.7~3.9) occurred at two main locations, the upper surface near the leading edge and the lower surface, near the trailing edge. Compared with the no-slip-wall case, in the immersed-solid case, the turbulent-kinetic energy of the fluid was higher, resulting in a greater flow loss of the fluid and greater drag coefficient on the foil.

6. Convergence-Ability Analysis

The convergence ability was analyzed based on the No. 5 mesh scheme in the no-slip-wall case and the No. 4 mesh combination in the immersed-solid case. The two computations were conducted on the same terminal, using 48-core 2.6 0GHz CPU and 32.0 GB memory. The solver settings were also the same for the two cases. Table 5 compares the total mesh nodes, element number and timestep number. When the cases were converged to RMS 1.0 × 10−5, the clock time until convergence, the clock time per timestep and the clock time per mesh node were also compared.
According to the mesh check in Section 4, the immersed-solid method needed more mesh nodes than using the no-slip wall to obtain equal accuracy. However, the clock time per mesh node in the immersed-solid case was shorter than that in the no-slip-wall case. The immersed-solid case needed a lower timestep number to be converged than the no-slip-wall case. However, the clock time per timestep was longer in the immersed-solid case than that in the no-slip-wall case. Generally, the immersed-solid case converged to RMS 1.0 × 10−5 in 1440.3 s, which was shorter than the 1782.6 s in the no-slip-wall case. The time consumption of the immersed-solid setup in this study took 80.7% of the clock time of the traditional no-slip-wall setup. Because the memory consumption was proportional to the mesh number, the immersed-solid setup in this study consumed 4.39 times the memory of the traditional no-slip-wall setup.

7. Conclusions

As an available option in the commercial software ANSYS CFX, the immersed-solid method can flexibly add wall boundaries in the fluid domain. In this study, under the typical hydrodynamic condition of the flow around the hydrofoil, without considering the hydrofoil movement, the immersed-solid method is compared with the non-slip-wall boundary, and the effectiveness of the immersed-solid method to simulate extremely small leakage, changing leakage, and randomly moving/rotating solid and intermittent attaching/detaching solids is verified, and the following conclusions are drawn. (A) The results of the immersed-solid method were qualitatively correct but quantitatively inaccurate. (B) “Boundary tracking” was crucial for improving the near-wall solution. (C) More mesh nodes were required for the qualitatively correct result when using the immersed-solid method than when using the no-slip wall. (D) The curve-smoothing process was recommended to correct the “pulsation” around the fluid–solid interface, due to data interpolation.
In detail, the flow-striking and separation positions on the solid can be correctly captured by both the no-slip-wall and immersed-solid method. However, the pressure fields and velocity fields performed differently, especially at the downstream location. A much wider wake could be found using the immersed-solid method than in the no-slip-wall case. Because of the added momentum source, the immersed-solid method prevented fluid passing through the hydrofoil. At this time, high turbulence-kinetic energy was triggered around the foil surface, especially at the flow-striking, separation and wake positions.
Generally, the immersed-solid method can accurately simulate the flow field by setting reasonable mesh nodes, especially on the fluid–solid interface. The convergence ability of the immersed-solid method is better than the no-slip-wall setup. If we apply the immersed-solid method to special cases such as extremely small leakage, changing leakage, randomly moving/rotating solid and intermittent attaching/detaching solids, it will provide us with a fast hydrodynamic solution, by overcoming the problems caused by excessive mesh density, mesh deformation or domain change. The research results of this paper have great engineering value, especially in the engineering-flow cases.

Author Contributions

Conceptualization, R.X. and R.T.; writing—original draft preparation, Z.H. and D.Z.; writing—review and editing, R.X. and D.Z.; supervision, R.T. and D.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the Open Research Subject of Key Laboratory of Fluid and Power Machinery (Xihua University), Ministry of Education (grant number LTDL-2022009).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request from the corresponding author.

Acknowledgments

The authors would like to acknowledge the help from Zhengwei Wang and his research group at Tsinghua University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kan, K.; Zhang, Q.; Zheng, Y.; Xu, H.; Xu, Z.; Zhai, J.; Muhirwa, A. Investigation into Influence of Wall Roughness on the Hydraulic Characteristics of an Axial Flow Pump as Turbine. Sustainability 2022, 14, 8459. [Google Scholar] [CrossRef]
  2. Tong, Z.; Xin, J.; Ling, C. Many-Objective Hybrid Optimization Method for Impeller Profile Design of Low Specific Speed Centrifugal Pump in District Energy Systems. Sustainability 2021, 13, 10537. [Google Scholar] [CrossRef]
  3. Wendt, J.F.; Anderson, J.D. Basic Philosophy of CFD. In Computational Fluid Dynamics: An Introduction; Spinger: Berlin/Heidelberg, Germany, 2009; pp. 5–8. [Google Scholar]
  4. Kajishima, T.; Taira, K. Numerical Simulation of Fluid Flows. In Computational Fluid Dynamics; Springer: Berlin/Heidelberg, Germany, 2017; pp. 1–22. [Google Scholar]
  5. Grotjans, H.; Menter, F.R. Wall Functions for General Application CFD Codes. Comput. Fluid Dyn. 1998, 1, 1112–1117. [Google Scholar]
  6. Kader, B.A. Temperature and Concentration Profiles in Fully Turbulent Boundary Layers. Int. J. Heat Mass Transf. 1981, 24, 1541–1544. [Google Scholar] [CrossRef]
  7. Patel, V.C.; Rodi, W.; Scheuerer, G. Turbulence Models for Near-Wall and Low Reynolds Number Flows—A Review. AIAA J. 1985, 23, 1308–1319. [Google Scholar] [CrossRef]
  8. Ha, T.W.; Choe, B.S. Numerical Simulation of Rotordynamic Coefficients for Eccentric Annular-Type-Plain-Pump Seal Using CFD Analysis. J. Mech. Sci. Technol. 2012, 26, 1043–1048. [Google Scholar] [CrossRef]
  9. Blanco, A.M.; Oro, J.F. Unsteady Numerical Simulation of an Air-Operated Piston Pump for Lubricating Greases Using Dynamic Meshes. Comput. Fluids 2012, 57, 138–150. [Google Scholar] [CrossRef]
  10. Akoz, M.S.; Kirkgoz, M.S.; Oner, A.A. Experimental and Numerical Modeling of a Sluice Gate Flow. J. Hydraul. Res. 2009, 47, 167–176. [Google Scholar] [CrossRef]
  11. Qun-feng, Z.; Pan-pan, Y.A.N.; Jian-ping, S.; Wan-fa, H.E. Numerical Simulation of Inlet and Outlet Pressures Influencing Internal Flow of a Vane Pump. Acta Armamentarii 2014, 35, 1223. [Google Scholar]
  12. Adkins, D.; Yan, Y.Y. CFD Simulation of Fish-like Body Moving in Viscous Liquid. J. Bionic. Eng. 2006, 3, 147–153. [Google Scholar] [CrossRef]
  13. Jin, C.; Xu, K.; Chen, S. A Three Dimensional Gas-Kinetic Scheme with Moving Mesh for Low-Speed Viscous Flow Computations. Adv. Appl. Math. Mech 2010, 2, 746–762. [Google Scholar] [CrossRef]
  14. Kajishima, T.; Takeuchi, S. Simulation of Fluid-Structure Interaction Based on an Immersed-Solid Method. J. Mech. Eng. Sci. 2013, 5, 555–561. [Google Scholar] [CrossRef]
  15. Yoon, Y.; Park, B.-H.; Shim, J.; Han, Y.-O.; Hong, B.-J.; Yun, S.-H. Numerical Simulation of Three-Dimensional External Gear Pump Using Immersed Solid Method. Appl. Therm. Eng. 2017, 118, 539–550. [Google Scholar] [CrossRef]
  16. Zhang, W.; Jiang, Q.; Bois, G.; Li, H.; Liu, X.; Yuan, S.; Heng, Y. Experimental and Numerical Analysis on Flow Characteristics in a Double Helix Screw Pump. Energies 2019, 12, 3420. [Google Scholar] [CrossRef]
  17. Qiu, R.; Zhang, W.; Jiang, Q.; Shen, K.; Gread, B.; Yuan, S. Research on Full Flow Passage Characteristics of Double Screw Pump Based on Immersion Entity Method. In Proceedings of the Journal of Physics: Conference Series; IOP Publishing: Bristol, UK, 2022; Volume 2186, p. 012015. [Google Scholar]
  18. Cervone, A.; Bramanti, C.; Rapposelli, E.; d’Agostino, L. Thermal Cavitation Experiments on a NACA 0015 Hydrofoil. J. Fluids Eng. 2006, 128, 326–331. [Google Scholar] [CrossRef]
  19. Dreyer, M.; Decaix, J.; Münch-Alligné, C.; Farhat, M. Mind the Gap: A New Insight into the Tip Leakage Vortex Using Stereo-PIV. Exp. Fluids 2014, 55, 1849. [Google Scholar] [CrossRef]
  20. Nowruzi, H.; Ghassemi, H.; Ghiasi, M. Performance Predicting of 2D and 3D Submerged Hydrofoils Using CFD and ANNs. J. Mar. Sci. Technol. 2017, 22, 710–733. [Google Scholar] [CrossRef]
  21. Hu, C.; Chen, G.; Yang, L.; Wang, G. Large Eddy Simulation of Turbulent Attached Cavitating Flows around Different Twisted Hydrofoils. Energies 2018, 11, 2768. [Google Scholar] [CrossRef]
  22. Day, M.A. The No-Slip Condition of Fluid Dynamics. Erkenntnis 1990, 33, 285–296. [Google Scholar] [CrossRef]
  23. Schmitt, F.G. About Boussinesq’s Turbulent Viscosity Hypothesis: Historical Remarks and a Direct Evaluation of Its Validity. Comptes Rendus Mécanique 2007, 335, 617–627. [Google Scholar] [CrossRef]
  24. Menter, F.R.; Kuntz, M.; Langtry, R. Ten Years of Industrial Experience with the SST Turbulence Model. Turbul. Heat Mass Transf. 2003, 4, 625–632. [Google Scholar]
  25. Launder, B.E.; Spalding, D.B. The Numerical Computation of Turbulent Flows. In Numerical Prediction of Flow, Heat Transfer, Turbulence and Combustion; Elsevier: Amsterdam, The Netherlands, 1983; pp. 96–116. [Google Scholar]
  26. Wilcox, D.C. Formulation of the K-w Turbulence Model Revisited. AIAA J. 2008, 46, 2823–2838. [Google Scholar] [CrossRef]
  27. Menter, F.R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
  28. Shu, J.-J.; Bin Melvin Teo, J.; Kong Chan, W. Fluid Velocity Slip and Temperature Jump at a Solid Surface. Appl. Mech. Rev. 2017, 69, 020801. [Google Scholar] [CrossRef]
  29. Shu, J.-J.; Teo, J.B.M.; Chan, W.K. A New Model for Fluid Velocity Slip on a Solid Surface. Soft Matter 2016, 12, 8388–8397. [Google Scholar] [CrossRef]
  30. Morris, D.L.; Hannon, L.; Garcia, A.L. Slip Length in a Dilute Gas. Phys. Rev. A 1992, 46, 5279–5281. [Google Scholar] [CrossRef]
  31. Pope, S.B. The Scales of Turbulent Motion. Turbul. Flows 2000, 1, 195–197. [Google Scholar]
  32. Pöschl, T. Mechanische Ähnlichkeit. In Lehrbuch der Technischen Mechanik für Ingenieure und Physiker; Springer: Berlin/Heidelberg, Germany, 1930; pp. 307–311. [Google Scholar]
  33. Yang, X.I.A.; Xia, Z.-H.; Lee, J.; Lv, Y.; Yuan, J. Mean Flow Scaling in a Spanwise Rotating Channel. Phys. Rev. Fluids 2020, 5, 074603. [Google Scholar] [CrossRef]
  34. Menter, F.; Esch, T. Elements of Industrial Heat Transfer Predictions. In Proceedings of the 16th Brazilian Congress of Mechanical Engineering (COBEM), Uberlândia, Brazil, 26–30 November 2001; Volume 109, p. 650. [Google Scholar]
  35. Abbott, I.H.; Von Doenhoff, A.E. Theory of Wing Sections; Dover Publications: Mineola, NY, USA, 1959; pp. 112–115. [Google Scholar]
  36. Arlinghaus, S. Practical Handbook of Curve Fitting; CRC Press: Boca Raton, FL, USA, 1994; ISBN 0-8493-0143-2. [Google Scholar]
Figure 1. Law of the near-wall velocity distribution. u: velocity; y: off-wall distance.
Figure 1. Law of the near-wall velocity distribution. u: velocity; y: off-wall distance.
Sustainability 15 03493 g001
Figure 2. Schematic map of the mesh elements and nodes when using the immersed-solid method.
Figure 2. Schematic map of the mesh elements and nodes when using the immersed-solid method.
Sustainability 15 03493 g002
Figure 3. The NACA0015 hydrofoil profile.
Figure 3. The NACA0015 hydrofoil profile.
Sustainability 15 03493 g003
Figure 4. The 3D-flow domain with NACA0015 foil at incidence angle of 8 degrees.
Figure 4. The 3D-flow domain with NACA0015 foil at incidence angle of 8 degrees.
Sustainability 15 03493 g004
Figure 5. Pressure on foil surface under different mesh schemes. (a) Cp distribution on foil surface along X direction; (b) Enlarged view within X/c = 0~0.15 and Cp = −0.5~−3.5.
Figure 5. Pressure on foil surface under different mesh schemes. (a) Cp distribution on foil surface along X direction; (b) Enlarged view within X/c = 0~0.15 and Cp = −0.5~−3.5.
Sustainability 15 03493 g005
Figure 6. Pressure on foil surface under different mesh-scheme combinations. (a) Cp distribution on foil surface along X direction; (b) Smoothed Cp curves.
Figure 6. Pressure on foil surface under different mesh-scheme combinations. (a) Cp distribution on foil surface along X direction; (b) Smoothed Cp curves.
Sustainability 15 03493 g006
Figure 7. The Cp contour on the X–Y plane at middle Z position. (a) No-slip wall case; (b) Immersed-solid case.
Figure 7. The Cp contour on the X–Y plane at middle Z position. (a) No-slip wall case; (b) Immersed-solid case.
Sustainability 15 03493 g007
Figure 8. The Cv contour on the X–Y plane at middle Z position. (a) No-slip-wall case; (b) Immersed-solid case.
Figure 8. The Cv contour on the X–Y plane at middle Z position. (a) No-slip-wall case; (b) Immersed-solid case.
Sustainability 15 03493 g008
Figure 9. The k contour on the X–Y plane at middle Z position. (a) No-slip-wall case; (b) Immersed-solid case. Unit of k: [m2/s2].
Figure 9. The k contour on the X–Y plane at middle Z position. (a) No-slip-wall case; (b) Immersed-solid case. Unit of k: [m2/s2].
Sustainability 15 03493 g009
Table 1. The acronyms defined for the immersed-solid method.
Table 1. The acronyms defined for the immersed-solid method.
AcronymsMeaning
IWIn-wall node
IW*In-wall point
NWNear-wall node
IBPoint on the immersed-solid boundary
FLFluid node
FL*Fluid point
NIBGElements between NW and IW nodes
NIBFElements between NW and FL nodes
Table 2. Mesh schemes checked in the no-slip-wall case.
Table 2. Mesh schemes checked in the no-slip-wall case.
No.Node NumberDescription
130,884Very Coarse
262,496Coarse
3107,652Mid-Coarse
4175,150Mid-Fine
5257,488Fine
6352,594Very Fine
Table 3. Mesh combinations checked in the immersed-solid case.
Table 3. Mesh combinations checked in the immersed-solid case.
No.Node NumberDescription
FluidSolidTotal
1234,24851,720291,368Coarse
2753,984378,8401,132,824Fine
3234,24851,720291,368Coarse with Boundary Tracking (Coarse-B)
4753,984378,8401,132,824Fine with Boundary Tracking (Fine-B)
Table 4. Lift force, drag force and torque on hydrofoil.
Table 4. Lift force, drag force and torque on hydrofoil.
No-Slip-Wall CaseImmersed-Solid Case
Lift Force, FY113.01 [N]90.82 [N]
Drag Force, FX4.95 [N]10.48 [N]
Lift/Drag Ratio, FY/FX22.838.67
Torque, TZ3.29 [N·m]2.31 [N·m]
Table 5. Convergence Ability Based on the Steady-State Simulation.
Table 5. Convergence Ability Based on the Steady-State Simulation.
No-Slip-Wall CaseImmersed-Solid Case
Total Mesh-Node Number, Nmn257,4881,132,824
Total Mesh-Element Number241,5421,077,732
Timesteps Converged to RMS, 1.0 × 10−5 Scc23681
Clock Time until Convergence, tcc1782.6 [s]1440.3 [s]
Clock Time per Timestep, tcc/Scc7.55 [s]17.78 [s]
Clock Time per Mesh Node, tcc/Mmn0.0069 [s]0.0013 [s]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, Z.; Zhu, D.; Xiao, R.; Tao, R. Comparative Evaluation of the Immersed-Solid Method for Simulating the Flow Field around Hydrofoil. Sustainability 2023, 15, 3493. https://doi.org/10.3390/su15043493

AMA Style

Hu Z, Zhu D, Xiao R, Tao R. Comparative Evaluation of the Immersed-Solid Method for Simulating the Flow Field around Hydrofoil. Sustainability. 2023; 15(4):3493. https://doi.org/10.3390/su15043493

Chicago/Turabian Style

Hu, Zilong, Di Zhu, Ruofu Xiao, and Ran Tao. 2023. "Comparative Evaluation of the Immersed-Solid Method for Simulating the Flow Field around Hydrofoil" Sustainability 15, no. 4: 3493. https://doi.org/10.3390/su15043493

APA Style

Hu, Z., Zhu, D., Xiao, R., & Tao, R. (2023). Comparative Evaluation of the Immersed-Solid Method for Simulating the Flow Field around Hydrofoil. Sustainability, 15(4), 3493. https://doi.org/10.3390/su15043493

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