• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.697-706
• /
• 2019

A multiphase flow modeling approach equipped with a hybrid turbulence modeling method is applied to compute the gravity currents in a rectangular channel. The present multiphase solver considers the dense fluid, the less-dense ambient fluid and the air above free surface as three phases with separate flow equations for each phase. The turbulent effect is simulated by the IDDES (improved delayed detach eddy simulation), a hybrid RANS/LES, approach which resolves the turbulent flow away from the wall in the LES mode and models the near wall flow in RANS mode on moderately fine computational meshes. The numerical results show that the present model can successfully reproduce the gravity currents in terms of the propagation speed of the current heads and the emergence of large-scale Kelvin-Helmholtz type interfacial billows and their three dimensional break down into smaller turbulent structures, even on the relatively coarse mesh for wall-modeled RANS computation with low-Reynolds number turbulence model. The present solutions reveal that the modeling approach can capture the large-scale three dimensional behaviors of gravity current head accompanied by the lobe-and-cleft instability at affordable computational resources, which is comparable to the LES results obtained on much fine meshes. It demonstrates that the multiphase modeling method using the hybrid turbulence model can be a promising engineering solver for predicting the physical behaviors of gravity currents in natural environmental configurations.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.187-190
• /
• 2006

In the present study, effects of tree-stream turbulence and surface trip wire on the flow past a sphere at $Re\;=\;0.4\;{\times}\;10^5\;{\sim}\;2.8\;{\times}\;10^5$ are investigated through wind tunnel experiments. Various types of grids are installed upstream of the sphere in order to change the tree-stream turbulence intensity. In the case of surface trip wire, 0.5mm and 2mm trip wires are attached from $20^{\circ}\;{\sim}\;90^{\circ}$ at $10^{\circ}$ interval along the streamwise direction. To investigate the flow around a sphere, drag measurement using a load cell, surface-pressure measurement, surface visualization using oil-flow pattern and near-wall velocity measurement using an I-type hot-wire probe are conducted. In the variation of free-stream turbulence, the critical Reynolds number decreases and drag crisis occurs earlier with increasing turbulence intensity. With increasing Reynolds number, the laminar separation point moves downstream, but the reattachment point after laminar separation and the main separation point are fixed, resulting in constant drag coefficient at each free-stream turbulence intensity. At the supercritical regime, as Reynolds number is further increased, the separation bubble is regressed but the reattachment and the main separation points are fixed. In the case of surface trip wire directly disturbing the boundary layer flow, the critical Reynolds number decreases further with trip wire located more downstream. However, the drag coefficient after drag crisis remains constant irrespective of the trip location.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.44-44
• /
• 2015

저수지나 하천 사면에서 발생하는 산사태와 토석류는 저수지와 하천 수체에 충격을 가한다. 이로 인해 발생하는 수면 충격파는 전파되어 반대편 제방으로 파의 처오름 또는 댐 제체위로의 물넘이로 큰 피해를 줄 수 있다. 최근 외국에서는 2차원 충격파 생성 및 전파의 기본 과정을 구명하기 위한 실험적 연구가 이뤄지고 있으며, 이들 연구들은 충격파의 발생과 전파, 사면활동 물질과 수체의 상호작용 그리고 자유 수면과 유속분표의 발달에 대한 자세한 관측 자료를 제시하고 있다. 아울러 충격파에 영향을 주는 지배 매개변수를 제시하고 있다. 하지만, 이러한 실험적 연구의 최근 진보에도 불구하고, 이들 지배 매개변수를 고려한 충격파 지배공식들은 대상 지역의 복잡한 바닥 지형이나, 평면적 지형 변화를 단순한 추정치로만 고려하게 된다. 따라서 복잡한 지형조건에서 토석류와 수체의 상호작용과 수면 충격파의 전파를 합리적으로 해석하는 데는 한계가 있다. 이 경우 수치모델링 기법을 대안으로 적용할 수 있으나, 수치모델링은 수면에서 충격파의 전파와 수중에서 토석류의 전파를 동시에 모의해야 하고, 뉴턴 유체와 비뉴턴 유체의 특성을 동시에 고려해야하므로 수치해석 연구자들에게는 하나의 큰 도전사항이다. 이 연구는 경계면 포착기법을 이용한 계산유체동력학 기법을 이용하여 사면활동과 이로 인한 정지 수역에서의 충격파의 발생 및 전파를 재현하기 위한 수치 모델링 기법을 개발하는 것이 목적이다. 사면활동과 수면의 경계면을 포착하고 위치를 정립하기 위해서 VOF (volume of fluid) 경계면 재구축 기법을 이용한다. 지배 방정식은 비압축성(incompressible) 질량 보존방정식과 나비어-스톡스(Navier-Stokes) 방정식이며, 서로 다른 유체의 상(phase)애 대한 체적분할이송방정식을 이용한다. 큰와 모의 계열의 난류 모델링 기법을 적용하여 충격파의 전파와 붕괴에 대한 난류의 영향을 고려하였다. 토석류는 비뉴턴 흐름저항 관계식을 적용하여 그 흐름특성을 재현하였다. 이들 지배방정식은 2차 정확도의 유한체적법(finite volume method)을 이용하여 해석한다. 외국의 연구자들이 관측하여 제시한 길이 11 m 그리고 폭 0.5 m의 수로에서 발생한 충격파를 수치적으로 재현하여 개발된 모형의 실제 문제에 대한 적용성을 보여준다.

• 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.32 no.8
• /
• pp.20-28
• /
• 2004

For the prediction on the onset of oblique shock wave-induced vortex breakdown, computational studies on the Oblique Shock wave/Vortex Interaction (OSVI) are conducted and compared with both experimental results and analytic mode1. A Shock-stable numerical scheme, the Roe scheme with Mach number-based function (RoeM), and a two-equation eddy viscosity-transport approach arc used for three-dimensional turbulent flow computations. The computational configuration is identical to available experiment, and we attempt to ascertain the effect of parameters such as a vortex strength, streamwise velocity deficit, and shock strength at a freestream Mach number of 2.49. Numerical simulations using the k-w SST turbulence model and suitably modeled vortex profiles are able to accurately reproduce many fine features through a direct comparison with experimental observations. The present computational approach to determine the criterion on the onset of oblique shock wave-induced vortex breakdown is found to be in good agreement with both the experimental result and the analytic prediction.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.11
• /
• pp.1169-1176
• /
• 2011

Liquid film flows are classified into waveless laminar, wavy laminar, and turbulent flows depending on the Reynolds number or the flow stability. Since the wavy motions of the film flows are so intricate and nonlinear, studies on them have largely been experimental. Most numerical approaches have been limited to the waveless flow regime. The various free surface-tracking schemes adopted for this problem were used to more accurately estimate the average film thickness, rather than to capture the unsteady wavy motion. In this study, the wavy motions in laminar wavy liquid film flows with Reynolds numbers of 200-1000 were simulated with various numerical schemes based on the volume of fluid (VOF) method for interface tracking. The results from each numerical scheme were compared with the experimental results in terms of the average film thickness, the wave velocity, and the wave amplitude.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.1 no.1
• /
• pp.33-45
• /
• 1997

Axisymmetric shear layers around a free jet is forced by co-flowing and counter-flowing secondary jets from/to an annular tube around the jet nozzle. The jet potential core extends far downstream with co-flowing secondary jets due to inhibited vortex developing and pairing. For counter-flowing cases, the axisymmetric shear layer around the jet transits from convective instability to absolute instability for velocity ratios R＝1.3～l.65 for the uniform velocity jets. Consequently, the jet potential core length increases and the turbulence level in the jet core is reduced significantly. The jets are controlled better with extension collars attached to the outer nozzle exit because the annular secondary flow is guided well by the extension collars. For the vectoring of jet, the annular tube around the jet is divided in two parts and the only one part is used for suction. The half suction makes the different shear layer around the jet and vectoring the jet by Coanda effect. The vectoring and turbulent components are varied significantly by the suction ratio. The experiments are carried out to investigate the characteristics of forced free jets using flow visualization, velocity and turbulence measurements.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.15 no.4
• /
• pp.773-779
• /
• 2011

The width of a submarine's turbulent wake, using Shear-free and Ship wake theory, is proportional to $x^n,\;({\frac{1}{5}}{\leq}n<{\frac{1}{2}})$ If we assume submarine's length, width, velocity are 65m, 6.5m, 6kts respectively, and the minimum diffusion of turbulent wake ; ${\infty}\;x^{1/5}$, the width of wake behind the submarine is about 20m at 1.2km, 30m at 15km when there is no breaking waves on the sea surface. However, in the case of breaking waves, it is very limited to identify submarine's wake on the sea surface because wind generated turbulent wake has higher turbulent kinetic energy than that of submarine's wake. As a result, there is a high possibility to detect submarine's wake on the sea surface in the shallow water such as the Yellow-Sea using a proper detection method such as SAR. This means that in anti-submarine operations, non-acoustic sea surface serveillance applied turbulent wake will be very effective way to detect a submarine in near future. To do this we have to develop exact theory of submarine's turbulent wake above all.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.42 no.7
• /
• pp.535-543
• /
• 2014

A computational simulation for a nonslender BWB UCAV configuration with rounded leading edge and span of 1.0m was performed to analyze its aerodynamic characteristics. The freestream is 50m/s over -4 to 26 degree A.o.A.s. Reynolds number based on the mean chord length is $1.25{\times}10^6$. 3D multi block hexahedral grids are used which allow good grid quality and ease to capture boundary layer. ${\gamma}-Re_{\theta}$ model as well as $k-{\omega}$ SST model is employed to assess the effect of transition for flow behavior. Drag and lift of the UCAV were well predicted while $C_M$ is under predicted at high angle of attacks and influenced by the turbulence models strongly. After assessing pressure distribution, skin friction lines and velocity field around the UCAV configuration, it was found that transition effect should be considered to enhance the prediction of aerodynamic behavior by a vortical flowfield.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.3
• /
• pp.575-583
• /
• 2017

Various types of flow patterns around the stepped weir and spillway, such as the skimming flow over such structures and the wave-type flow with a standing undular hydraulic jump and roller downstream of the structures, are developed in open channels. Unsteady three-dimensional numerical simulations are carried out using a hybrid RANS-LES turbulence modeling approach and the volume of fluid method for resolving free surface fluctuations to represent the turbulent flow including the skimming flow and wave-type flow over a stepped weir installed in a rectangular channel. The comparison of numerical results with an existing experimental measurement reveals that the present numerical simulations reasonably well reproduce the turbulent flow passing the stepped weir, in terms of time-averaged velocity profiles at selected locations downstream of the weir, flow topology characterized by the wave-type and skimming flows, the maximum height and length of the standing wave and the length of reattachment of recirculating zone. The numerical result further elucidates the distinct flow behaviors of the wave-type and skimming flow by presenting instantaneous intense variations of free surface and velocity vectors, the distributions of Reynolds shear stress and turbulent kinetic energy and three-dimensional complex features of coherent structures and total pressure distribution.