• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.41 no.2
• /
• pp.87-94
• /
• 2017

Because the cavitation flow driven by an underwater propeller corrodes the materials around it and generates a high level of noise, it has become an important topic in engineering research. In this study, computational fluid dynamics techniques are applied to simulate cavitation flow, and the noise in the flow is predicted by applying the acoustic analogy to the predicted flow. The predicted results are compared with measurement results and other predictions in terms of surface pressure distribution and the temporal variation in liquid volume fraction. The predicted results are found to be in good agreement with the measured results. The source of the noise attributed to the time rate of change in the liquid volume fraction around the hydrofoil is modeled as a monopole source, and the source of the noise due to unsteady pressure perturbations on the hydrofoil surface is modeled as a dipole source. Then the predicted noise results are analyzed in terms of directivity and SPL spectrum. The noise caused by unsteady pressure perturbations was dominant in the entire frequency range considered in the study.

• Proceeding of EDISON Challenge
• /
• 2015.03a
• /
• pp.569-574
• /
• 2015

지난 수십 년간 유체역학적인 관점에서 곤충이나 새의 움직임을 모방하기 위해 진동하는 익형(pitching airfoil)과 동적 실속에 관한 많은 연구가 진행되어 왔다. 그러나 유동박리가 일어나지 않는 범위 내에서 진동하는 익형의 특성에 대한 연구는 보기 드물다. 또한 기존의 유동박리가 일어나지 않는 영역에서 익형의 진동 현상에 대해 수행된 연구는 수중과 같이 낮은 레이놀즈수에서 수렴되었기 때문에, 공기 중과 같이 높은 레이놀즈수에서 유동현상과 다른 특성을 보여주고 있을 수 있다. 따라서 본 연구는 높은 레이놀즈수에서의 다양한 환산 진동수, 받음각진폭, 익형에 따른 공력특성을 분석하였다. 그 결과, 익형의 진동으로 인한 양력계수의 차이는 작음을 알 수 있었다. 그러나 높은 환산 진동수에서 익형의 항력계수가 감소하는 경향이 나타나며, 이로 인해 높은 환산 진동수에서 수치적으로 추력이 발생할 수 있음을 확인하였다.

• The KSFM Journal of Fluid Machinery
• /
• v.11 no.3
• /
• pp.7-13
• /
• 2008

The cavitating flow simulation is of practical importance for many engineering systems, such as pump, turbine, nozzle, Infector, etc. In the present work, a solver for two-phase flows has been developed and applied to simulate the cavitating flows past hydrofoils. The governing equation is the two-phase Navier-Stokes equation, comprised of the continuity equation of liquid and vapor phase. The momentum and energy equation is in the mixture phase. The solver employs an implicit, dual time, preconditioned algorithm using finite difference scheme in curvilinear coordinates. An experimental data and other numerical data were compared with the present results to validate the present solver. It is concluded that the present numerical code has successfully accounted for two-phase Navier-Stokes model of cavitation flow.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.26 no.2
• /
• pp.141-147
• /
• 2016

Underwater cavitation is one of the most important issues because it causes not only vibration and erosion of submerged bodies but also significant flow noise problems. In this paper, flow noise due to cavitation flows around the NACA66 MOD hydrofoil is numerically investigated. The cavitation flow simulation is conducted using the Reynolds-Averaged Navier-Stokes equations based on finite difference methods. To capture the cavitation phenomena accurately and effectively, the homogeneous mixture model with the Merkle's cavitation model is applied. The predicted results are compared with available experimental data in terms of pressure coefficients and volume fraction, which confirms the validity of numerical results. Based on flow field analysis results, hydro-acoustic noise field due to the cavitation flow is predicted using the Ffowcs-Williams and Hawkings equation derived from the Lighthill's acoustic analogy. The typical lift dipole propagation patterns are identified.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.17-23
• /
• 2014

Cavitation causes a great deal of noise, damage to components, vibrations, and a loss of efficiency in devices, such as propellers, pump impellers, nozzles, injectors, torpedoes, etc. Thus, the cavitating flow simulation is of practical importance for many engineering systems. In the present work, a two-phase flow solver based on the homogeneous mixture model has been developed. The solver employs an implicit preconditioning, dual time stepping algorithm in curvilinear coordinates. The flow characteristics around Clark-Y hydrofoil were calculated and then validated by comparing with the experimental data. The lift and drag coefficients with changes of angle of attack and cavitation number were obtained. The results show that cavity length and lift, drag coefficient increase with increasing angle of attack.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.20-24
• /
• 2011

In this paper, numerical simulation of cavitation flow for modified NACA66 hydrofoil was made by using the multi-phase RANS equation based on pseudo-compressibility. The Homogeneous mixture model comprised of the mixture continuity, mixture momentum and liquid volume fraction equations was utilized. A vertex-centered finite-volume method was used in conjunction 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing The Spalart-Allmaras one equation model was employed for the closure of turbulence. Reasonable agreements were obtained between the calculation results and the experiment for pressure coefficients on the hydrofoil surface.

• The Journal of the Acoustical Society of Korea
• /
• v.37 no.5
• /
• pp.263-270
• /
• 2018

In most industry fields concerning external flow noise problems, the hybrid computational aeroacoustic techniques based on the FW-H (Ffowcs Williams and Hawkings) equation are widely used for its numerical efficiency. However, when the surface integral form of FW-H equation is used without volume quadrupole sources, it is known to generate significant non-physical noise in a certain case. Especially, in the case of a flow in which the tip vortex cavitation is formed in the distant downstream direction such as flow driven by an underwater propeller, the accuracy in noise prediction becomes poor unless it is not properly modelled. Therefore, in this study, the nonphysical acoustic waves caused by the surface integral form of FW-H equation is reduced by adding the quadrupole correction term. First, to verify the accuracy of the in-house code of FW-H equation, the noise by an axial fan used in the outdoor unit of air conditioner was calculated and compared with the results of ANSYS Fluent. In order to verify the effects of the quadrupole correction term, the noise prediction for isentropic vortex convection is performed and it is confirmed that the error is reduced by the quadrupole correction term. Finally, the noise prediction is performed for the flow field generated by the Clark-Y hydrofoil in underwater. It is confirmed that the error caused by the cavitation passing through the integral surface can be reduced by the quadrupole correction term.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.133-140
• /
• 1998

A Navier-Stokes code based on a unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number $k-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for a domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. Unsteadiness inside boundary layer is entrained when a unsteady vortex impinge on the blade surface. It shoves that local peak value inside the boundary layer and also local minimum near the edge of boundary layer as it developes along the blade surface. The unsteadiness inside the boundary layer is almost isolated from the free stream unsteadiness and being convected at local boundary layer speed, less than the free stream value.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.777-784
• /
• 2000

A Navier-Stokes code based on an unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number ${\kappa}-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for the whole experimental domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. The first harmonics of the velocity in the boundary layer shows local peak value inside the boundary layer and also local minimum near the edge of boundary layer. It is intensified as it develops along the blade surface. This is shown to be caused as the unsteadiness inside the boundary layer is being convected at a speed less than the free stream value. It is also shown that there is negligible mixing of the unsteadiness between the boundary layer and the free stream.