• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.267-277
• /
• 1993

A finite element method(FEM) is presented and applied to the two-dimensional bottom turbulent boundary layer. The time-dependent incompressible motion of a viscous fluid is formulated by using the well-known Navier-Stokes equations and vorticity equation in terms of the velocity and pressure fields. The general numerical formulation is based on Galerkin method and solved by introducing the mixing length theory of Prandtl for eddy kinematic viscosity of a turbulent flow field. Numerical computations of the transport of sediment on an arbitrary sea-bed due to wave motion in the turbulent boundary layer are carried out. The results obtained by the FEM made clear the difference in characteristic features between the boundary layer due to oscillatory flow and the boundary layer due to wave motion.

• 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) 기법을 적용하였으며 흐름 조건별로 입자를 움직이게 하는 난류 이벤트는 달라짐을 확인하였다.

• Proceeding of EDISON Challenge
• /
• 2012.04a
• /
• pp.93-96
• /
• 2012

경계층이란 유체와 물체 표면의 마찰로 인해 생성되는 층을 말한다. 경계층은 두께에 따라 층류 경계층, 천이 경계층, 난류 경계층으로 나누어진다. 레이놀즈 수 크기에 따라 경계층은 몇 가지의 층으로 구분되어 진다. 이 계산에서는 경사진 평판 위에서 유동의 현상들이 어떻게 일어나는지 확인하였다. 또한, 경사가 없는 평판위에서 velocity profile과 Blasius solution을 비교하였고, 평판의 뒤쪽에 격자의 간격이 넓음으로 큰 오차가 발생하게 됨을 알 수 있었다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.6
• /
• pp.477-483
• /
• 2016

Velocity profiles of laminar, transition and turbulent boundary layers were investigated by using Particle Image Velocimetry(PIV) measurements on the flat plate at Mach 2.96. The Schlieren visualization and PIV measurements are also used to confirm whether the oblique shock wave generated from the leading edge affects the flow field over the flat plate. The laminar velocity profile measured from the experiment was well matched with the compressible Blasius solution. The velocity profile of the transition boundary layer was well correlated with the theoretical turbulent velocity profile from near the wall and the transition began from Re = $1.41{\times}106$. For the turbulent boundary layer, considering compressibility effects, the Van Driest-transformed velocity satisfies the incompressible log-law. It is found that the log region is extended farther in the wall-normal direction compared to the log region in incompressible boundary layer.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.3
• /
• pp.269-276
• /
• 1985

본 논문에서는 이러한 두꺼운 3차원 난류 경계층의 체계적인 연구를 위하여 적합한 수학선형을 개발하고 저속풍동에서 경계층 실험을 수행하였다. 이 수학선형 은 가능한 실제 선체주위의 유동특성이 잘 나타나도록 설계되었다. 실험을 통하여 전체적 유동을 파악하고 표면에서 평균 속도분포, 마찰저항계수 및 각종 적분변수들을 측정하였다.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.19-26
• /
• 1991

The complete, fully-elliptic Reynolds-averaged Navier-Stokes equations have been solved using a two-layer model, in the $\kappa-\varepsilon$ turbulence model, for the axisymmetric body. Numerically generated boundary-fitted coordinate system and the finite analytic methods are used to solve the governing equations. Calculations are started after the middle body with given inlet conditions. The velocities and the turbulent quantities at the inlet section are specified by solving the boundary layer equations or by standard flat-plate boundary profiles. The effects of the inlet conditions on the solution are investigated.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.1
• /
• pp.69-84
• /
• 2020

Natural shoreline repeats its re-treatment and advance in response to the endlessly varying sea-conditions, and once severely eroded under stormy weather conditions, natural beaches are gradually recovered via a boundary layer streaming when swells are prevailing after storms cease. Our understanding of the boundary layer streaming over surf-zone often falls short despite its great engineering value, and here it should be noted that the most sediments available along the shore are supplied over the surf-zone. In this rationale, numerical simulation was implemented to investigate the hydraulic characteristics of boundary layer streaming over the surf zone in this study. In doing so, comprehensive numerical models made of Spatially filtered Navier-Stokes Eq., LES (Large Eddy Simulation), Dynamic Smagorinsky turbulence closure were used, and the effects of turbulence closure such as Dynamic Smagorinsky in LES and k-ε on the numerically simulated flow field were also investigated. Numerical results show that due to the intrinsic limits of k-ε turbulence model, numerically simulated flow velocity near the bottom based on k-ε model and wall function are over-predicted than the one using Dynamic Smagorinsky in LES. It is also shown that flow velocities near the bottom are faster than the one above the bottom which are relatively free from the presence of the bottom, complying the typical boundary layer streaming by Longuet-Higgins (1957), the spatial scope where boundary layer streaming are occurring is extended well into the surf zone as incoming waves are getting longer. These tendencies are plausible considering that it is the bottom friction that triggers a boundary layer streaming, and longer waves start to feel the bottom much faster than shorter waves.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.56-65
• /
• 2013

Two compression ramp problems and an impinging shock problem are computed to investigate influence of turbulence models and eddy viscosity on the shock-wave / boundary layer interaction. A Navier-Stokes boundary layer generation code was applied to the generation of inflow boundary conditions. Computational results are validated well with the experimental data and effects of turbulence models are investigated. It is shown that the behavior of turbulence (eddy) viscosity directly affects both the extent of the separation and shock-wave positions over the separation.

• Journal of Korea Water Resources Association
• /
• v.44 no.3
• /
• pp.189-198
• /
• 2011

The aim of this study is to compare k-$\varepsilon$ and k-$\omega$ closures under the condition of oscillatory layer flow at high Reynolds number. A one dimensional vertical model incorporated with flow momentum equations and turbulence models (k-$\varepsilon$ and k-$\omega$) is applied to the laboratory measurements in the turbulent oscillatory boundary layer. The numerical simulation reveals that both turbulence models calculate similar velocity profiles and turbulent kinetic energy (TKE). In addition, both deliver high accuracy under the condition of negligible spanwise pressure gradient. Therefore, it is recommended in this study to use k-$\varepsilon$ closure, of which numerical coefficients have been calibrated from many studies, for the cases of straight channel, estuary, and coastal environment where the spanwise pressure gradient is not significant.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.18 no.6
• /
• pp.1049-1056
• /
• 2023

The height of the atmospheric boundary layer indicates the peak developed when turbulence is generated by mixing heat and water vapor, and is generally determined through thermodynamic methods. Wind profilers produce atmospheric information from the scattering of signals sent into the atmosphere. A method for making the spectrum of turbulent components, turbulent kinetic energy dissipation rate, and refractive index structure coefficient was presented to determine the atmospheric boundary layer depth. Compared with the vertical distribution characteristics of potential temperature and specific humidity based on radiosonde data, the determination method of the atmospheric boundary layer height from wind profiler output was evaluated as very useful.