• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.643-651
• /
• 2008

In the present paper, turbulent open-channel flows over longitudinal bedforms are numerically simulated. The Reynolds- averaged Navier-Stokes equations in curvilinear coordinates are solved with the non-linear $k-{\varepsilon}$ model by Speziale( 1987). First, the developed model is applied to rectangular open channel flows for purposes of model validation and parameter sensitivity studies. It is found that the parameters $C_D$ and $C_E$ are important to the intensity of secondary currents and the level of turbulent anisotropy, respectively. It is found that the non-linear $k-{\varepsilon}$ model can hardly reproduce the turbulence anisotropy near the free surface. However, the overall pattern of the secondary currents by the present model is seen to coincide with measured data. Then, numerical simulations of turbulent flows over longitudinal bedforms are performed, and the simulated results are compared with the experimental data in the literature. The simulated secondary currents clearly show upflows and downflows over the ridges and troughs, respectively. The numerical results of secondary currents, streamwise mean velocity, and turbulence structures compare favorably with the measured data. However, it is observed that the secondary currents towards the troughs were significantly weak compared with the measured data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.281-281
• /
• 2022

바닥에서 생성되는 난류는 순간적으로 강한 모멘텀을 바닥에 전달함과 동시에 바닥에 있는 입자를 움직이게 한다. 경계층 내 난류 운동에 대한 분석은 다양한 유사 이송 문제를 이해하기 위해 필수적이며 이에 따라 많은 선행 연구들은 실험실 실험을 통해 해당 연구를 수행하였다. 본 연구에서는 선행 연구에서 사용하지 못했던 진보된 실험 방법을 활용하여 바닥 경계층 내의 난류 운동에 대해 확인하고 해당 운동에 의해 관성 입자의 움직임이 어떻게 발생하는지에 대하여 물리적으로 설명하고자 한다. 다양한 흐름 조건에서 3가지의 입경 크기를 가지는 모래 입자를 가지고 실험을 수행하였으며, 실험 조건별 고해상도 유속장 및 관성 입자의 움직임은 3차원 입자 영상 유속계 (Particle Image Velocimetry; 이하 PIV)와 입자 추적 유속계 (Particle Tracking Velocimetry; 이하 PTV)를 동시에 적용하여 파악하였다. 취득된 3차원 유속장과 입자 궤적을 기반으로 실험 조건별 흐름 및 입자 거동 특성에 대해 분석하였으며, 관성 입자의 움직임을 발생시키는 3차원 난류 운동은 측정된 유속장에서 산정한 Q-criterion 값을 기반으로 도식화하였다. 측정값 내에는 난류 운동에 대한 정보와 더불어 잡음이 포함되어 있으므로 이를 제거하고자 적합 직교 분해 (Proper Orthogonal Decomposition; 이하 POD) 방법을 적용하였다. 그리고 POD로 추출한 유속장을 통해 바닥면 부근에 존재하는 헤어핀 와류 운동 혹은 와류 묶음과 같은 난류 고유 구조를 파악하였다. 해당 와류 운동들의 3차원 난류 특성을 확인하고자 비등방성 불변 지도(anisotropy invariant map)를 활용하였으며 경계층 내부에서 난류의 형태가 흐름 방향으로 늘어진 럭비공 형태임을 확인하였다. 마지막으로, 입자의 움직임을 발생시키는 난류 이벤트를 결정하고자 사방구 분석 (Quadrant analysis) 기법을 적용하였으며 흐름 조건별로 입자를 움직이게 하는 난류 이벤트는 달라짐을 확인하였다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.1063-1071
• /
• 1999

Fine grid calculations are reported for the developing turbulent flow in a curved duct of square cross-section with a radius of curvature to hydraulic diameter ratio ${\delta}=Rc/D_H=3.357 $ and a bend angle of 720 deg. A sequence of modeling refinements is introduced; the replacement of wall function by a fine mesh across the sublayer and a low Reynolds number algebraic second moment closure up to the near wall sublayer in which the non-linear return to isotropy model and the cubic-quasi-isotropy model for the pressure strain are adopted; and the introduction of a multiple source model for the exact dissipation rate equation. Each refinement is shown to lead to an appreciable improvement in the agreement between measurement and computation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.3
• /
• pp.280-293
• /
• 1998

Numerical study on turbulent and mean structures of a turbulent boundary layer with longitudinal and spanwise pressure gradient is carried out by using Reynolds-stress-model (RSM). The existence of pressure gradient in a turbulent boundary layer causes the skewing or divergence of rates of strain, which contributes to production of turbulent kinetic energy. Also, this augmentation of production due to extra rates of strain can increase the turbulent mixing and cause the anisotropy of turbulent intensities in the outer layer. This paper uses the Reynolds Stress Model to capture anisotropy of turbulent structures effectively and is devoted to compare the results computed by using RSM and the standard k-.epsilon. model with experimental data. It is concluded that the RSM can produce the more accurate predictions for capturing the anisotropy of turbulent structure than the standard k-.epsilon. model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.3
• /
• pp.833-840
• /
• 2014

This study presents numerical simulation of turbulent flows in a rectangular open-channel that has a width-to-depth ratio of 2 using the source code provided by OpenFOAM. Large eddy simulations are carried out by solving the filtered continuity and momentum equations numerically. For the non-isotropic residual stress term, Smagorinsky's (1963) model is used. The flow in the open-channel whose width-to-depth ratio is 2, from experiment of Tominaga et al. (1989), is simulated numerically. Simulation results are compared with measured data by Tominga et al. (1989) and Nezu and Rodi (1985) and with LES data by Shi et al. (1999). Comparisons revealed that the model simulates the mean flow and turbulence statistics well. Specifically, the model reproduced the inner secondary currents located at the corner of sidewall and free surface successfully. In addition, the vortical component of turbulence intensity shows bulged contours towards the bottom edge.

• Nuclear Engineering and Technology
• /
• v.27 no.4
• /
• pp.491-502
• /
• 1995

In accurate prediction of cross-flow based on detailed knowledge of the velocity field in subchannels of a nuclear fuel assembly is of importance in nuclear fuel performance analysis. In this study, the low-Reynolds number k-$\varepsilon$ turbulence model has been adopted in too adjacent subchannels with cross-flow. The secondary flow is accurately estimated by the anisotropic algebraic Reynolds stress model. This model was numerically calculated by the finite element method and has been verified successfully through comparison with existing experimental data. Finally, with the numerical analysis of the velocity Held in such subchannel domain, an analytical correlation of the lateral loss coefficient is obtained to predict the cross-flow rate in subchannel analysis codes. The correlation is expressed as a function of the ratio of the lateral How velocity to the donor subchannel axial velocity, recipient channel Reynolds number and pitch-to-diameter.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.655-659
• /
• 2008

본 연구에서는 3차원 수치모의를 통하여 횡방향 하상형상 및 격자형 이차흐름 구조의 생성 메커니즘을 분석하였다. 이를 위해 곡선좌표계에서의 지배방정식을 구성하고, 난류 폐합을 위해 Speziale(1987)가 제안한 비등방성 k-$\varepsilon$모형을 이용하였다. 또한 Exner 방정식을 이용하여 시간에 따른 하상변동을 예측하였다. 그 결과 바닥과 측벽 사이에서 발생되는 바닥 이차흐름의 하향류에 의해 측벽 부근부터 하상이 침식되고, 침식된 유사량은 이차흐름의 횡방향 유속에 의해 이동되어 퇴적되어, 결국 횡방향으로 연속적인 언덕 및 저면과 같인 하상형상이 생성되는 것으로 나타났다. 또한 시간에 따른 이차흐름 및 바닥 전단력, 하상고의 변화에 대해 살펴보았다.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.1-8
• /
• 1997

Two versions of anisotropic k-ε turbulence model are incorporated in the modified k-ε model of Sohn et al. to avoid the need for the experimental normal stress value in the model and applied to convergent and divergent flows with strong and adverse pressure gradients in the plane of symmetry of a body of revolution. The models are the nonlinear k-ε model of Speziale and the anisotropic model of Nisizima & Yoshizawa. All of the models yield satisfactory results for relatively complex flow on a plane-of-symmetry boundary layer. The results of the models are compared with those results of experimental normal stress value.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.473-478
• /
• 2000

In this study, the characteristics of the three-dimensional turbulent flow in a rotating square sectioned $90^{\circ}$ bend were investigated by numerical simulation and experiment. In the experimental study, the characteristics of a developing turbulent flow are measured using hot-wire anemometer to seize the rotational effects on the flow characteristics and to compare the results of computational simulation with Reynolds stress model. Each refinement is shown to lead to an appreciable improvement in the agreement between measurement and computation.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.771-776
• /
• 2000

Fine grid calculations are reported for the developing turbulent flow in a straight duct and a curved duct of square cross-section with a radius of curvature to hydraulic diameter ratio ${\delta}=R_c/H_H=3.357$ and a bend angle of 180 deg. A sequence of modeling refinements is introduced; the replacement of wall function by a fine mesh across the sublayer and a low Reynolds number second moment closure up to the near wall sublayer in which the non-linear return to isotropy model and the cubic-quasi-isotropy model for the pressure strain are adopted; and the introduction of a multiple source model for the exact dissipation rate equation. Each refinement is shown to lead to an appreciable improvement in the agreement between measurement and computation.