• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.104-110
• /
• 2008

Detached Eddy Simulation (DES) is applied to an axisymmetric base flow at supersonic mainstream. DES is a hybrid approach to modeling turbulence that combines the best features of the Reynolds-averaged Navier-Stokes (RANS) and large-eddy simulation (LES) approaches. In the Reynolds-averaged mode, the model is currently based on either the Spalart-Allmaras (S-A) turbulence model. In the large eddy simulation mode, it is based on the Smagorinski subgrid scale model. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology with less computational cost than that of pure LES and monotone integrated large-eddy simulation (MILES) approaches. The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region, small-eddy motions inside the recirculating region. Comparison of the results shows that it is necessary to resolve approaching boundary layers and free shear-layer velocity profiles from the base edge correctly for the accurate prediction of base flows. The consideration of an empirical constant CDES for a compressible flow analysis may suggest that the optimal value of empirical constant CDES may be larger in the flows with strong compressibility than in incompressible flows.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.3
• /
• pp.29-36
• /
• 2002

A numerical analysis based on two-dimensional and three-dimensional incompressible Wavier-Stokes equations has been carried out for double-circular-arc compressor cascades and the results are compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE algorithm adopt pressure weighted interpolation method for non-staggered grid and hybrid scheme for the convective terms. Turbulence modeling is very important for prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. Considering computation times, $\kappa$-$\varepsilon$ turbulence model with wall function is used.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.416-419
• /
• 2006

Numerical analysis is carried out to identify the wall thinning effect inside the feed water valve. The finite volume method is applied to make analysis for the viscous flows. The commercial cock FLUENT is used for the simulation and the GAMBIT for the grid generation. The RNG $\kappa-\varepsilon$ model is used for the turbulence and the tet-hybrid grid is applied for the modeling. The velocity vector, the pressure contour, the change of residual along the iteration number, and the dynamic head are predicted for the hydrodynamic investigation.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.09b
• /
• pp.842-845
• /
• 2005

A method for the numerical simulation of two-dimensional free-surface flow resulting from the propagation of regular gravity waves over topography with arbitrary bottom shape is presented. The method is based on the numerical solution of the Euler equations subject to the fully nonlinear free-surface boundary conditions and the appropriate bottom, inflow and outflow conditions using a hybrid finite-differences and spectral-method scheme. The formulation includes a boundary-fitted transformation, and is suitable for extension to incorporate large-eddy simulation (LES) and large-wave simulation (LWS) terms for turbulence and breaking wave modeling, respectively. Results are presented for the simulation of the free-surface flow over two different bottom topographies, with constant slope values of 1:10 and 1:20, two different inflow wave lengths and two different inflow wave heights. An absorption outflow zone is utilized and the results indicate minimum wave reflection from the outflow boundary. Over the bottom slope, lengths of waves in the linear regime are modified according to linear theory dispersion, while wave heights remain more or less unchanged. For waves in the nonlinear regime, wave lengths are becoming shorter, while the free surface elevation deviates from its initial sinusoidal shape.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.149-149
• /
• 2020

하천에서 물 흐름이 보와 댐과 같은 수공구조물을 지날 때 일반적으로 흐름상태에 다양하고 급진적인 변화가 발생한다. 특히 흐름이 구조물을 지나면서 사류(supercritical flow)로 변하고 다시 상류(subcritical flow)로 복원되면서 일어나는 도수(hydraulic jump) 현상은 수위의 급변화, 흐름 에너지 소산, 변동성이 강한 압력 분포 등이 특징이다. 이러한 흐름 특성들은 보나 여수로와 같은 수공구조물 자체의 성능뿐만 아니라 이들 수공구조물의 하류에서 발생하는 국부세굴로 인해 구조물의 안정성에 부정적인 영향을 줄 수 있다. 따라서 수공구조물을 설계할 때는 이들 구조물을 통과하는 흐름의 비정상 난류 흐름 특성을 정확하게 해석하여 반영하여야 한다. 이 연구에서는 k-omega SST 난류 모형과 자유수면의 급격한 변동을 해석하기 위한 하이브리드-VOF(hybrid volume of fluid)기법을 이용하여 도수현상을 수치적으로 재현하고자 한다. 기존 CFD(computational fluid Dynamics) 모델링에서는 자유수면 변동의 영향을 고려하기 위해 VOF 기법을 많이 사용하였다. 하지면 전통적인 VOF 기법은 다상흐름(multiphase flow)을 오직 부피분율(volume fraction)의 함수로만 고려하며 모의함으로써 강한 수면변동뿐만 아니라 공기연행(air entrainment)를 동반하는 난류 흐름을 모의하는데는 한계가 있다. 이 연구에서 이용하는 Eulerian 기법인 하이브리드 VOF 기법은 물과 공기의 각 상에 대하여 흐름 특성들을 개별적으로 계산하기 때문에 공기연행을 포함한 급변류 흐름에서 전통적인 VOF 기법보다 적용성이 우수하다. 이와 같은 난류모형과 자유수면 포착기법을 이용하여 3차원 비정상 난류 흐름 수치모형을 구축하여 수공구조물 주변에서 발생하는 강한 공기연행과 난류 특성를 보이는 급변류를 수치적으로 재현한다. 이 연구는 계산된 수치해석 결과를 기존 수리실험 결과와 비교하여 수치모형의 적용성을 평가하고 도수 현상에서 발생하는 독특한 흐름 특성을 제시한다.

• 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.

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

The flow in the meandering channel is characterized by the spiral motion of secondary currents that typically cause the erosion along the outer bank. Hydraulic structures, such as spur dike and groyne, are commonly installed on the channel bottom near the outer bank to mitigate the strength of secondary currents. This study is to investigate the effects of submerged vanes installed in a $90^{\circ}$ meandering channel on the development of secondary currents through three-dimensional numerical modeling using the hybrid RANS/LES method for turbulence and the volume of fluid method, based on OpenFOAM open source toolbox, for capturing the free surface at the Froude number of 0.43. We employ the second-order-accurate finite volume methods in the space and time for the numerical modeling and compare numerical results with experimental measurements for evaluating the numerical predictions. Numerical results show that the present simulations well reproduce the experimental measurements, in terms of the time-averaged streamwise velocity and secondary velocity vector fields in the bend with submerged vanes. The computed flow fields reveal that the streamwise velocity near the bed along the outer bank at the end section of bend dramatically decrease by one third of mean velocity after the installation of vanes, which support that submerged vanes mitigate the strength of primary secondary flow and are helpful for the channel stability along the outer bank. The flow between the top of vanes and the free surface accelerates and the maximum velocity of free surface flow near the flow impingement along the outer bank increases about 20% due to the installation of submerged vanes. Numerical solutions show the formations of the horseshoe vortices at the front of vanes and the lee wakes behind the vanes, which are responsible for strong local scour around vanes. Additional study on the shapes and arrangement of vanes is required for mitigate the local scour.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.2
• /
• pp.79-93
• /
• 2015

Three-dimensional structures of large ocean rings in the Gulf Stream region are investigated using the HYbrid Coordinate Ocean Model (HYCOM). Numerically simulated flow structures around four selected cyclonic and anticyclonic rings are compared with two different horizontal resolutions: $1/12^{\circ}$ and $1/48^{\circ}$. The vertical distributions of Lagrangian Coherent Structures (LCSs) are analyzed using Finite Size Lyapunov Exponent (FSLE) and Okubo-Weiss parameters (OW). Curtain-shaped FSLE ridges are found in all four rings with extensions of surface ridges throughout the water columns, indicating that horizontal stirring is dominant over vertical motions. Near the high-resolution rings, many small-scale flow structures with size O(1~10) km are observed while these features are rarely found near the low-resolution rings. These small-scale structures affect the flow pattern around the rings as flow particles move more randomly in the high-resolution models. The dispersion rates are also affected by these small-scale structures as the relative horizontal dispersion coefficients are larger for the high-resolution models. The absolute vertical dispersion rates are, however, lower for the high-resolution models, because the particles tend to move along inclined eddy orbits when the resolution is low and this increases the magnitude of absolute vertical dispersion. Since relative vertical dispersion can reduce this effect from the orbital trajectories of particles, it gives a more reasonable magnitude range than absolute dispersion, and so is recommended in estimating vertical dispersion rates.