• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.7
• /
• pp.930-938
• /
• 2018

Underwater noise radiated from submarines is directly related to the probability of being detected by the sonar of an enemy vessel. Therefore, minimizing the noise of a submarine is essential for improving survival outcomes. For modern submarines, as the speed and size of a submarine increase and noise reduction technology is developed, interest in flow noise around the hull has been increasing. In this study, a noise analysis technique was developed to predict flow noise generated around a submarine shape considering the free surface effect. When a submarine is operated near a free surface, turbulence-induced noise due to the turbulence of the flow and bubble noise from breaking waves arise. First, to analyze the flow around a submarine, VOF-based incompressible two-phase flow analysis was performed to derive flow field data and the shape of the free surface around the submarine. Turbulence-induced noise was analyzed by applying permeable FW-H, which is an acoustic analogy technique. Bubble noise was derived through a noise model for breaking waves based on the turbulent kinetic energy distribution results obtained from the CFD results. The analysis method developed was verified by comparison with experimental results for a submarine model measured in a Large Cavitation Tunnel (LCT).

• Water for future
• /
• v.29 no.4
• /
• pp.187-198
• /
• 1996

A numerical model for unsteady dispersion of horizontal line source in turbulent shear flow is developed. A fractional step finite difference method is used which splits the unsteady two-dimensional advective diffusion equation into the longitudinal advection and the vertical diffusion equations, and solves them alternately for half time intervals by the Holly-Preissmann scheme and the Crank-Nicholson scheme, respectively. The developed numerical model is verified using a semi-analytic solution for steady dispersion in turbulent shear flow. Dispersion of an instantaneous plane source in turbulent shear flow is analyzed using the model. The degree of mixing at the same dimensionless time is almost the same regardless of the friction factor, and the travel distance required to reach a certain degree of mixing is inversely proportional to the square root of the friction factor.

• 유체기계공업학회:학술대회논문집
• /
• 1998.02a
• /
• pp.179-185
• /
• 1998

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.105-105
• /
• 2019

보나 여수로와 같은 수공구조물의 주변에서 발생하는 흐름 거동은 구조물 모서리에서 발생하는 흐름분리(flow separation)와 이에 따른 전단층(shear layer)과 재순환(recirculation) 흐름 영역의 발달 그리고 분리된 흐름의 재부착(reattachment)이 특징이다. 공학적으로 난류의 해석에 있어서 이러한 흐름 거동들을 정확하게 예측하는 것은 수공구조물 설계에 있어서 중요하다. 이 연구에서는 흐름 분리와 재순환 영역의 발달 그리고 흐름 재부착을 포함하는 후방계단(backward-facing step) 흐름을 155,000의 레이놀즈수 조건에서 하이브리드 RANS/LES 모델을 적용하여 해석결과를 평가한다. 하이브리드 모델로는 벽에 인접한 격자의 해상도에 상대적으로 민감하지 않은 SST(shear-stress transport) 난류 모델을 이용하는 DES(detached-eddy simulation) 기법을 적용하였다. 계단 높이가 h인 계산영역은 흐름방향 길이가 34h, 높이는 계단 상류와 하류에서 각각 1h와 2h 그리고 폭은 $2{\pi}$이다. 계단은 상류단으로부터 10h 하류부 지점에 위치한다. 경계조건으로 상부와 하부 벽면에 대해서는 비활조건을 적용한다. 상류부 수로에서 완전 발달한 흐름을 재현하기 위해서 유입경계조건은 유입부 하류 $2{\pi}h$ 지점에서 계산된 유속과 난류량을 매핑(mapping)기법을 이용하여 반복적으로 적용한다. 총 3.1백만개와 7.3백만개의 셀로 계산영역을 구현한 두 개의 계산격자 그리고 약 3.1백만개의 셀을 이용했지만 벽면 근처에서의 격자 구성을 다른 방식으로 설정한 두 가지 격자를 이용하여 격자 해상도가 DES 수치해석 결과에 미치는 영향을 분석하였다. 수치해석결과는 본 연구에서 상류단 조건으로 적용한 매핑기법이 대상 수로에서 완전 발달한 흐름을 잘 재현함을 보여주며, 합리적인 DES 해석 결과를 얻기 위해서는 벽에 수직한 방향으로 적절한 격자의 해상도와 분포가 필요함을 보여준다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.1
• /
• pp.24-31
• /
• 2006

In the present study, a strongly coupled analysis code is developed for transonic flutter analysis. For aerodynamic analysis, two dimensional Reynolds-Averaged Navier-Stokes equation was used for governing equation, and ε-SST for turbulence model, DP-SGS(Data Parallel Symmetric Gauss Seidel) Algorithm for parallelization algorithm. 2 degree-of-freedom pitch and plunge model was used for structural analysis. To obtain flutter response in the time domain, dual time stepping method was applied to both flow and structure solver. Strongly coupled method was implemented by successive iteration of fluid-structure interaction in pseudo time step. Computed results show flutter speed boundaries and limit cycle oscillation phenomena in addition to typical flutter responses - damped, divergent and neutral responses. It is also found that the accuracy of transonic flutter analysis is strongly dependent on the methodology of fluid-structure interaction as well as on the choice of turbulence model.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.2
• /
• pp.97-104
• /
• 2005

A simple but efficient grid generation technique by using the modified compressible form of stream function has been formulated. Transformation of a physical plane to a streamline plane, the Von Mises Transformation, has been widely used to solve the differential equations governing flow phenomena, however, limitation arises in low velocity region of boundary layer, mixing layer and wake region where the relatively large grid spacing is inevitable. Modified Von Mises Transformation with simple mathematical adjustment for the stream function is suggested and applied to solve the confined coaxial turbulent jet mixing with simple $\kappa-\epsilon$ turbulence model. Comparison with several experimental data of axial mean velocity, turbulent kinetic energy, and Reynolds shear stress distribution shows quite good agreement in the mixing layer except in the centerline where the turbulent kinetic energy distributions were somewhat under estimated. This formulation is strongly suggested to be utilized specially for free turbulent mixing layers in axisymmetric flow conditions such as the investigation of mixing behavior, jet noise production and reduction for Turbofan engines.

• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.1
• /
• pp.23-33
• /
• 2020

In this study, the effect of surface roughness distribution and its influence on the inflight icing code was investigated. Previous numerical studies focused on the magnitude of surface roughness, and the effects were only addressed in terms of changes in thermal boundary layers with fully turbulent assumption. In addition, the empirical formula was used to take account the turbulent transition due to surface roughness, which was regarded as reducing the accuracy of ice shape prediction. Therefore, in this study, the turbulent transition model based on the two-equation turbulence model was applied to consider the effects of surface roughness. In order to consider the effect of surface roughness, the transport equation for roughness amplification parameter was applied, and the surface roughness distribution model was implemented to consider the physical properties. For validation, the surface roughness, convective heat transfer coefficient, and ice shape were compared with experimental results and other numerical methodology. As a result, it was confirmed that the excessive prediction of the heat transfer coefficient at the leading edge and the ice horn shape at the bottom of the airfoil were improved accordingly.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.225-234
• /
• 1995

터빈익렬내부의 유동해석을 위해 비압축성 점성유동해석을 이용한 수치 해석 프로그램을 개발하였다. 지배방정식으로는 2차원의 비정상 비압축성 Navier-Stokes 방정식을 일반화된 곡선좌표계로 전환하여 암시적으로(implicitly) 반복적인 시간진행방법을 이용하여 유동해석을 하였다. 지배방정식의 각항들은 시간에 대해 1차의 정확도 그리고 영역에 대해서는 2차의 정확도, 대류항에 대해서는 3차의 정확도를 가지는 Upwind기법을 적용하였다. 특히, 실험적 접근이 매우 어려운 터빈의 정익과 회전하고 있는 동익과의 상호운동을 멀티블럭기법과 데이터 interface를 통해 보다 쉽게 해석할 수 있었다. 본 연구결과는 정익만을 계산한 타 연구자의 결과와의 비교시 매우 일치하였으며 물리적인 유동을 잘 파악할 수 있었다. 난류유동 해석을 위해서 Baldwin-Lomax 모델을 적용하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.134-134
• /
• 2011

3차원 난류 부력젵의 혼합을 큰 와 모의(large-eddy simulation) 기법을 이용하여 수치모의 한다. 개발된 수치모형은 3차원 열동수역학 모형을 이용하여 부력젵의 퍼짐, 자기 보존 그리고 주변류의 연행 등을 포함하는 난류젵의 동적 특성을 분석할 수 있다. 수치해석에서 하부격자규모 (subgrid scale, SGS) 난류 응력은 부력항을 수정한 Smagorinsky 모형을 이용한다. 여과된 엔탈피 수송방정식에서 하부격자규모의 스칼라 플럭스는 상수의 SGS Prandtl 수를 가지는 단순 경사수송 가설에 근거하여 모의한다. 계산된 결과를 실험결과와 비교하며, 결과는 양호하게 일치함을 보여준다. 계산결과에 따르면 부력항의 수정이나 SGS 난류 Prandtl 수는 결과에 큰 영향을 미치지 않지만 SGS 모형 상수인 Cs 값은 부력젵 확산 예측에 중요한 영향을 미치는 것으로 나타났다.