• Ocean Systems Engineering
• /
• 제11권1호
• /
• pp.43-58
• /
• 2021

In the present work, a 3D numerical model is proposed to study local scouring around single vertical piers with different cross-section shapes under steady-current flow. The model solves the flow field and sediment transport processes using a coupled approach. The flow field is obtained by solving the Unsteady Reynolds Averaged Navier-Stokes (URANS) equations in combination with the k-ω SST turbulence closure model and the sediment transport is considered using both bedload and suspended load models. The proposed model is validated against the empirical measurements of local scour around single vertical piers with circular, square, and diamond cross-section shapes obtained from the literature. The measurement of scour depth in equilibrium condition for the simulations reveal the differences of 4.6%, 6.7% and 13.1% from the experimental measurements for the circular, square, and diamond pier cases, respectively. The model displayed a remarkable performance in the prediction of scour around circular and square piers where horseshoe vortices (HSVs) have a leading impact on scour progression. On the other hand, the maximum deviation was found in the case of the diamond pier where HSVs are weak and have minimum impact on the formation of local scour. Overall, the results confirm that the prediction capability of the present model is almost independent of the strength of the formed HSVs and pier cross-section shapes.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of Mechanical Science and Technology
• /
• 제20권8호
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• 한국수자원학회논문집
• /
• 제41권7호
• /
• pp.669-679
• /
• 2008

본 연구는 수리실험으로 얻어진 권파에 의한 월파수괴의 유속을 댐붕괴흐름과 비교하여 거동의 유사성을 검토하였다. 댐붕괴흐름은 해석해가 간략하고 월파 거동과 유사함으로 인해 월파의 유속산정에 이용되어왔다. 월파는 일반적으로 많은 연행기포로 인해 기존의 유속측정기법을 적용하는데 제한을 받게 되므로, 본 실험에서는 기포나 기포조직모양을 이용한 기포영상유속계를 이용하여 월파 유속을 측정하였다. 실험결과로부터 월파의 유속단면을 검토하였고, 단면의 최대유속과 수심평균유속을 댐붕괴흐름의 1차원 해석해인 Ritter의 해와 비교하였다. 해석해와의 비교를 위해 댐붕괴 전 초기수심을 파랑조건과 월파의 측정값으로부터 유도하였다. 파랑조건으로부터 추정된 붕괴 전 초기 수심을 이용한 해석해는 측정된 유속의 분포형태에 있어서는 차이를 보였으나, 월파수괴의 전면속도로부터 산정된 초기수심을 이용한 해석해는 유속크기의 비교에서 좋은 일치를 보였다.

• 대한토목학회논문집
• /
• 제44권3호
• /
• pp.329-336
• /
• 2024

우리나라에 설치된 다기능보의 경우 가동보 구간은 수문의 부분 개방 시 수문 개방 높이와 하류의 관리수위에 의해 자유도수 또는 수중도수 형태의 흐름이 발생한다. 본 연구에서는 수문 아래를 지나 흐르며 발생하는 수중도수를 수치모의하고 평균흐름, 난류량, 그리고 상대수심 등에 대하여 분석하였다. 수치모의를 위하여 unsteady Reynolds-averaged Navier-Stokes 방정식과 volume of fluid 기법, 그리고 k-ω SST 난류모형을 이용하였다. 기존에 수행된 다른 연구자들의 실험 결과를 이용하여 수치모형을 검증하여 수치모형이 도수에서 발생하는 이상흐름을 적절히 모의하는 것을 검토하였다. 또한 내부 평균흐름 및 난류량의 분포에 대하여 모의하여 분포 형태에 대해 분석하고, 자유수면과 재순환영역의 길이 등을 분석하였다.

• Geomechanics and Engineering
• /
• 제8권6호
• /
• pp.811-828
• /
• 2015

Proper modelling of the basal resistance terms is key in simulating the motion of fluidized granular flow. In this paper, standard depth-averaged governing equations of granular flow are used together with the classical Coulomb, Voellmy, and velocity dependent friction models (VDFM). A high-resolution modified TVDLF method is implemented to solve the partial differential equations without numerical oscillations. The effects of basal resistance terms on the motion of granular flows such as geometric shape evolution, travel times and final deposits are analyzed. Based on the numerical results, the predictions of the front and rear end positions and developing length of granular flow with Coulomb friction model show excellent agreements with experiment results reported by Hutter et al. (1995), and illustrate the validity of the numerical approach. For the Voellmy model, the higher value of turbulent coefficient than reality may obtain more reasonable predicted runout for the small-scale avalanche or granular flow. The energy exchange laws indicate that VDFM is different from the Coulomb and Voellmy models, although the flow characteristics of both three models fit the measurements and observations very well.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2006년도 학술발표회 논문집
• /
• pp.1925-1930
• /
• 2006

흐름의 연직분포를 가정하여 수면에서 바닥까지 적분하여 평균화한 방정식에 기초한 준2차원 흐름모형을 연안해역에서의 취송류해석을 위해 적용되었다. 바람장의 변동에 따른 취송류장의 변동이 크게 나타나므로 본 연구에서는 육상지형의 영향을 고려할 수 있는 3차원 중규모 지역기상장 모형을 통해 연안해역에서의 시간에 따른 바람장변동을 수치해석 하였다. 육상에서의 지표면 열수지 변화를 위해 식생모델을 적용하였다. 기상장과 취송류 해석을 위한 수치모형실험결과 바람장의 변동에 따른 연안역에서의 흐름특성을 효과적으로 예측하고 있음을 확인하였다.

• Water Engineering Research
• /
• 제1권3호
• /
• pp.187-197
• /
• 2000

This paper describes an adaptive quadtree-based 2DH wave-current interaction model which is able to predict wave breaking, shoaling, refraction, diffraction, wave-current interaction, set-up and set-down, mixing processes (turbulent diffusion), bottom frictional effects, and movement of the land-water interface at the shoreline. The wave period-and depth-averaged governing equations are discretised explictly by means of an Adams-Bashforth second-order finite difference technaique on adaptive hierarchical staggered quadtree grids. Grid adaptation is achieved through seeding points distributed according to flow criteria(e.g. local current gradients). Results are presented for nearshore circulation at a sinusoidal beach. Enrichment permits refined modelling of important localised flow features.

• 한국환경과학회지
• /
• 제23권7호
• /
• pp.1223-1232
• /
• 2014

Density currents have been easily observed in environmental flows, for instance turbidity currents and pollutant plumes in the oceans and rivers. In this study, we explored the propagation dynamics of density currents using the FLOW-3D computational fluid dynamics code. The renormalization group (RNG) $k-{\varepsilon}$ scheme, a turbulence numerical technique, is employed in a Reynold-averaged Navier-Stokes framework (RANS). The numerical simulations focused on two different types of intrusive density flows: (1) propagating into a two-layer ambient fluid; (2) propagating into a linearly stratified fluid. In the study of intrusive density flows into a two-layer ambient fluid, intrusive speeds were compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting propagation speed of the density currents. We also numerically reproduced the effect of the ratio of current depth to the overall depth of fluid. The numerical model provided excellent agreement with the analytical values. It was also clearly demonstrated that RNG $k-{\varepsilon}$ scheme within RANS framework is able to accurately simulate the dynamics of density currents. Simulations intruding into a continuously stratified fluid with the various buoyancy frequencies are carried out. These simulations demonstrate that three different propagation patterns can be developed according to the value of $h_n/H$ : (1) underflows developed with $h_n/H=0$ ; (2) overflows developed when $h_n/H=1$ ; (3) intrusive interflow occurred with the condition of 0 < $h_n/H$ < 1.

• 한국수자원학회논문집
• /
• 제43권9호
• /
• pp.823-833
• /
• 2010

침수식생 개수로 흐름은 식생영역과 상부영역에 서로 다른 흐름구조를 보인다. 즉, 식생영역에서 전단으로 인해 생성되는 난류는 억제되며 비교적 균일한 유속 분포를 보이며 상부영역에서는 일반 개수로 흐름과 유사한 흐름구조를 보인다. 이와 같이 두 상이한흐름구조가 결합된 복잡한 흐름특성으로 인해 침수식생 개수로흐름은 공학적인 관심의 대상이 되어왔다. 본 연구에서는 침수식생 개수로 흐름의 층적분 모형의 비교 분석을 수행하였다. 일반적으로 식생흐름의 층적분 모형은 층의 수에 따라 2층 및 3층모형으로 구분한다. 즉, 전체 수심을 식생영역과 상부영역으로 구분하는 2층모형과 식생영역을 바닥 조도의 영향 유무에 따라 내부 및 외부 식생영역으로 구분하는 3층모형으로 분류된다. 본 연구에서는 2층모형과 3층모형을 비교하였다. 다양한 실험조건에 적용한 결과, 3층모형이 식생영역에서 유속의 변화를 고려할 수 있으나 결과는 레이놀즈응력 분포에 민감하며, 적분된 유속은 2층모형에 의한 예측 결과가 더욱 정확한 것으로 나타났다. 3층모형에서 내부 식생영역의 결과가 전체 흐름구조에 미치는 영향이 무시할 수 있으므로 이 점을 착안하여 식생영역에서 유속 변화가 고려되는 수정 2층모형을 제시하였다. 수정 2층모형에서 가정하는 레이놀즈응력 분포는 상부영역에서는 선형, 식생영역에서는 멱함수 형으로 변화한다. 다양한 조건에 적용한 결과, 수정 2층모형이 대체로 기존의 모형과 비슷한 정도의 예측을 수행하나 식생밀도가 매우 작은 흐름의 경우 예측 결과가 불량한 것으로 나타났다.