• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.376-380
• /
• 2004

측벽이 존재하는 개수로 난류흐름에 대한 DNS 자료를 이용하여 난류의 비등방성을 해석하였다. 측벽의 2등분선(sidewall bisector)에서 난류강도의 분포를 통해 바닥과 자유수면에서의 비등방성을 분석하고, AIM을 도입하여 흐름장 전반에 걸쳐 세부적인 비등방성 해석을 수행하였다. 측벽의 2등분선에서의 난류강도의 분포를 통해 바닥과 자유수면 근처에서 난류강도가 강한 비등방성을 갖는 것을 볼 수 있었다. 또한 3개의 다른 영역에서 AIM의 도입을 통해 측벽이나 바닥근처에서는 난류의 비등방성이 구형 관수로 흐름과 유사한 것으로 나타났으나, 개수로 난류흐름의 주된 특성이 velocity-dip phenomena가 존재하는 영역에서는 구형 관수로 흐름과는 달리 천이영역이 존재한다는 것을 알 수 있었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.138-138
• /
• 2011

식생된 개수로에서 식생의 영향을 파악하기 위해 k-$\in$ 난류 모형을 이용하여 수치모의를 하였다. 식생의 영향을 고려하기 위해 항력항을 추가한 지배방정식을 구성하였으며, 지배방정식을 해석하기 위하여 유한체적법을 사용하였다. 수치모의에서 구한 식생된 개수로의 흐름구조를 기존의 수리실험 결과와 비교하여 비교적 잘 일치함을 확인할 수 있다. 난류의 생성과 소멸을 수치모의한 결과, 부분구간 식생된 경우 식생높이 보다 낮은 구간에서는 후류에 의한 난류 생성이 지배적이며, 식생높이보다 높은 구간에서는 주로 마찰에 의한 난류 생성이 지배적임을 보였다. 기존의 연구들은 식생의 영향을 고려하여 개수로의 흐름을 연구한 예는 드물며, 현재까지 진행되어진 국내의 연구는 난류모형을 이용하여 식생된 개수로에서의 흐름 구조를 모의하였다. 따라서 난류흐름을 모의하는데 가장 보편적인 k-$\in$ 난류모형을 이용하여 식생된 개수로에서 수직방향으로의 흐름구조와 식생의 영향을 해석하는 것은 그 자체로도 의미 있는 연구이며, 앞으로의 환경수리 문제를 해결하기 위해 선행되어야 하는 연구이다. 식생된 개수로에서의 난류구조와 부유사 이동에 대한 식생의 영향을 비정상 1차원 수직모형으로 해석하였으며, 폐합문제를 위해 2-방정식인 k-$\in$ 난류모형을 사용하였다. k-$\in$ 난류모형에 식생에 의한 항력항을 더하여 지배방정식을 구성하였다. 수직방향에 대해 흐름방향 유속 u, 난류에너지 k, 그리고 난류에너지 소산율 $\in$의 분포를 구하고, 부유사에 대한 수송방정식을 풀었다. 식생된 개수로와 식생되지 않은 개수로에서의 유속분포, 난류강도, 레이놀즈 응력 분포와 난류의 생성과 소멸을 구하여 식생이 난류흐름에 미치는 영향을 분석하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.438-443
• /
• 2007

개수로 흐름에 대한 연구는 1700년대 중엽 Chezy에 의해 이론적인 기초를 다졌으며, 광범위하고 조직적인 관측연구는 Darcy(1803-1858)에 의해 약 150년 전에 폭 2m, 길이 600m에 이르는 수로에서 실험관측을 수행하고자 시도하였다. Darcy의 동료이자 후계자인 Bazin은 제한된 조건의 현장관측뿐 아니라 다양한 조건의 수로를 제작하여 실험관측을 수행하였으며 그의 실험자료는 Bazin 자신 뿐만 아니라 Manning이나 Ganguillet와 Kutter 등 여러 연구자들의 경험식 개발에도 이용되었다. Nikuradse(1933)의 균일조도 원형관수로 실험결과로부터 관로 흐름은 층류, 천이층류, 완난류, 천이난류, 전난류 등 다섯 종류의 흐름특성을 가지고 있음을 확인하였고 Bazin(1829-1917)과 varwick(1945)의 실험결과로부터 개수로 흐름에서도 관로 흐름과 마찬가지로 층류, 천이층류, 완난류, 천이난류, 전난류 등 다섯가지의 흐름특성이 존재함을 알 수 있다. 본 연구에서는 기존의 지수형 마찰계수 산정식에서 단순히 조고만의 함수였던 ${\alpha}$에 물의 기본적인 성질인 표면장력, 점성력, 밀도와 자연 하천의 경사, 수심, 수면, 폭, 조고의 영향을 고려한 수심에 관한 무차원수 $Y_h$, 수면 폭에 관한 무차원수 $Y_b$, 조고에 관한 무차원수 $Y_k$를 도입하였다. 따라서 개수로 흐름 해석에 있어서 기존의 마찰계수 산정치보다 여러 영향을 반영하여 정확한 값을 산정할 수 있을 것으로 판단된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.127-127
• /
• 2017

침수식생이 식재된 개수로에서 식생밀도에 따라 유동 및 난류의 특성이 변화된다. 이러한 특성은 식생에서의 유사, 영양물질, 용존 산소 등에 영향을 미치며 수중 생물의 서식에 변화를 준다. 따라서 침수식생이 식재된 개수로 흐름을 이해하는 것은 중요하게 여겨지고 있으며 많은 선행연구자들에 의해 실험 및 수치모의를 통해 활발히 연구되고 있다. 하지만 대부분의 RANS(Reynolds-Averaged Navier-Stokes)를 기반으로 한 선행연구에서는 침수식생의 흐름 특성을 반영하지 않은 모형을 이용하여 정확한 모의 결과를 도출하지 못 하였다. 이에 정확한 침수식생 흐름을 수치모의하기 위해서는 침수식생 흐름의 특성을 반영한 지배방정식을 이용해야 한다. 본 연구의 목적은 침수식생 흐름의 특성을 RANS 모형 중의 하나인 SA (Spalart-Allmaras) model에 반영하고, 식생밀도에 따른 유동 및 난류량을 실측치와 비교하는 것이다. RANS 방정식을 이용하여 난류모델링을 하였으며, 난류폐합문제를 해결하기 위해서 modified SA model을 이용하였다. 침수식생에서의 흐름을 해석하기 위해 운동량방정식에 식생항력을 추가하였다. 선행연구자의 식생수로 실험을 바탕으로 모형검증을 하였으며, 식생밀도에 따라 평균유속 및 난류구조를 확인하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.65-65
• /
• 2021

보를 월류하는 흐름에 대하여 하류 수심이 보의 높이보다 낮고, 보 월류 후 수심의 공액수심보다 작거나 크고 같은 경우를 각각 자유도수와 수중도수라고 한다. 수중도수가 발생하는 경우 하류 수심이 보의 경사면을 따라 흐르는 하강류를 덮게 되며 유속을 감쇄시키기 때문에 최대 유속은 자유도수에 비해 매우 작게 발생한다. 그러나 수중도수는 자유도수에 비해 에너지 감쇄 효율이 낮기 때문에 도수구간의 거리가 증가하게 된다. 따라서 인명피해를 발생시키는 도수구간 재순환영역의 길이를 검토하는 것이나 보 하류 바닥보호공 길이 설계를 위해 수중도수에서의 흐름에 대하여 검토하는 것은 중요하다. 본 연구에서는 k-ω SST 난류모형을 이용하여 보 월류 후 발생하는 수중도수를 수치모의하고 평균흐름과 난류량의 종방향 변화에 대하여 검토하였다. 기존 실험수로에 k-ω SST 난류모형을 사용하여 모형의 적용성을 검토하였다. 다양한 하류 수심을 설정하여 평균흐름과 난류량에 대한 침수도의 영향 및 자유도수 계산결과와 벽면 제트 결과를 함께 비교하였다. 검토 결과 수중도수는 평균흐름과 난류량의 변화율이 자유도수보다 작고 벽면 제트보다는 큰 것을 확인하였다. 또한 침수도가 증가되면서 평균흐름과 난류량의 변화율이 작아지는 것을 확인하였다. 이것은 침수도의 변화에 따른 역압력경사의 차이에 의한 것으로 판단된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.136-136
• /
• 2019

원형 실린더를 주변 흐름에 관한 연구는 오랜 기간 유체역학 전 영역에서 모형실험이나 수치모형으로 광범위하게 연구되었다. 이 흐름은 하천의 교각이나, 바다의 시추선과 같은 수공구조물 주변에서 관측된다. 난류와 와류가 공존하는 복잡한 특성 때문에, 이 흐름은 수공학에서 유사이송, 세굴, 오염물 확산 등에 영향을 준다. 본 연구는 실험실 수로에 설치된 원형 실린더(D=9cm) 후방의 근접 와류 구간에서(x/D<5) 유속을 ADV로 측정한 후, 난류 특성을 Power Spectral Estimation(PSE)와 Singular Spectral Analysis(SSA) 방법으로 연구하였다. PSE는 샘플 스펙트럼의 한계를 보완하고자 자료를 분할하고, window 함수를 적용하여 ensemble 평균을 구하는 경험적 방법이다. PSE를 이용하여 스펙트럼을 계산한 결과, 주 흐름 및 횡방향 흐름은 Inertial subrange에서 Kolmogorov의 가정과 일치하는 추세를 보였다. 그러나 수심방향 흐름의 스펙트럼은 -5/3보다 빠르게 감소하는 추세를 보였다. Inertial subrange 스펙트럼에서 난류 에너지 소산율은 원형 실린더에서 멀어짐에 따라 감소하는 추세를 보였고, 주 흐름방향과 횡방향 흐름은 비슷한 크기를 보였다. 난류 에너지 소산율과 동점성계수를 이용하여 Kolmogorov 길이, 유속, 시간 스케일을 계산했다. 난류의 운동에너지를 계산하기 위해 Triple decomposition 방법 중 하나인 SSA를 적용하였다. SSA는 유속행렬을 이용하여 고윳값과 고유벡터를 계산하고, 유속에서 기여도가 큰 부분을 추출하는 방법이다. SSA를 통해 실린더 후방 흐름에서 와류 성분과 난류 성분을 나누었다. 횡방향 흐름은 강한 와류로 큰 기여도를 갖는 고유벡터가 나타났지만, 주 흐름과 수심방향 흐름은 상대적으로 낮은 기여도를 갖는 고유벡터가 나타났다. 와류를 제외한 흐름에서 난류 운동에너지는 실린더와 멀어짐에 따라 감소하고, 흐름 중앙에서(y/D=0) 가장 큰 값을 보였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.24-24
• /
• 2021

자연하천 바닥 경계층 내에서는 복잡한 난류 구조가 형성되며 이들은 하상에 강한 모멘텀을 전달한다. 바닥 부근에 분포하는 유사 입자들은 경계층 내에서 발생한 난류 흐름으로부터 모멘텀을 전달받아 소류사 혹은 부유사 형태로 이송되게 되며, 이러한 유사 이송 과정을 역학적으로 설명하기 위해서는 경계층 내 유체 흐름에 대한 이해가 선행되어야한다. 경계층 내 난류 흐름 특성이 유사 입자의 움직임에 미치는 영향에 대해 분석하기 위해서는 바닥 경계층 내 고해상도 유속 자료와 유사 움직임을 동시에 포착할 수 있는 기술이 요구된다. 하지만 현재까지 수행된 대부분의 선행 연구들은 점 유속을 측정할 수 있는 음파 도플러 유속계 (Acoustic Doppler Velocimetry) 혹은 2차원 입자 영상 유속계를 활용하였으며, 이들은 복잡한 3차원 난류 흐름 특성을 분석하기에는 한계가 있다. 본 연구의 목적은 실험실 실험을 통해 바닥 경계층 내 3차원 난류 흐름이 유사 이송에 미치는 영향에 대해 조사하는 것이다. 본 연구에서는 유사 주변에서의 고해상도 3차원 흐름 유동장 및 순간적인 유사 움직임에 대해서는 합성 개구 (synthetic aperture) 기반의 3차원 입자 영상 유속계 및 입자 추적 유속계를 활용하여 취득하였다. 취득된 흐름 유동장을 기반으로 레이놀즈 전단응력을 산정하였으며 이를 통해 유체가 하상에 미치는 모멘텀의 크기를 파악하였다. 복잡한 난류 흐름 구조에 대해서는 팔분원 분석 (octant analysis)을 통해 구분했으며, 유사가 움직이는 순간의 유속장을 기반으로 유사 이송을 발생시키는 지배적인 난류 흐름 특성에 대해 규명하였다. 본 연구는 바닥 경계층 내 복잡한 3차원 난류 흐름과 유사 입자의 움직임을 동시에 분석함으로써 기존에 수행되어왔던 선행 연구들의 한계점을 극복하고 보다 명확한 유사 이송의 발생 원인에 대해 분석했다는 점에 의의가 있다.

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

사행수로에서의 원심력은 비선형적인 압력분포를 야기하여 이차류, 편수위 등과 같은 불규칙하고 복잡한 흐름을 발생시킨다. 일반적으로 이들 흐름은 난류이고 매우 3차원적이며 자유수면과의 상호작용이 중요한 역할을 할 수도 있다. 환경, 유사이동, 지형 변화와 관련된 환경 수리학적 관점에서 사행수로에서의 흐름을 이해하고 설계하기 위해서는 이러한 복잡한 3차원 흐름을 정확하게 계산하는 것이 매우 중요하다. 이 연구에서는 유한차분법에 근거한 3차원 흐름해석 모형을 이용하여 사행수로에서의 난류 흐름을 모의하고자 한다. 지배방정식은 3차원 비정상 RANS(Reynolds averaged Navier-Stokes) 방정식이며, 난류 해석을 위해서 공학적으로 널리 이용되고 있는 난류 모형 중 k-omega 모형을 이용한다. 수치모형은 시간과 공간에 대해서 2차 정확도의 이산화 기법을 적용한다. 자유수면의 변동은 이상(two-phase) VOF (volume of fluid) 기법을 이용하여 계산한다. 수치모형의 적용 대상은 기존 문헌에서 제시되어 있는 키노시타 사인곡선을 이용하여 만든 폭 60cm의 사행수로에서 후르드수 0.23 그리고 레이놀즈수 41,700의 조건에서 발생시킨 난류 흐름이다. 적용한 난류모형들을 이용하여 해석한 결과들을 유속벡터분포와 수위의 항으로 비교분석하여 사행수로에서 발생되는 이차류와 편수위 변화 재현에 대한 수치모형의 적용성을 평가하고 각 난류모형들의 특성을 제시한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.313-313
• /
• 2023

이 연구에서는 침수식생 조건에서 식생 형태 별 frontal area, solid volume fraction이 유속 분포에 미치는 영향을 분석하고, 흐름측정결과로부터 식생 형태에 따른 난류흐름특성을 분석하기 위하여 수행 되었다. 식생흐름 구현을 위하여 5 cm의 간격으로 총 257개의 모형식생을 전체 영역에 배치했다. 유속측정위치는 수위측정결과에 따라 흐름이 안정화되는 구간에서 연직방향으로 17개 지점에서 측정한 후 앙상블 평균하여 분석했다. Branch의 유무에 따라 Type I과 II로 구분하여 각 식생에 대해 유속의 연직분포를 측정한 결과, Branch가 없는 Type I에서는 유속이 지속적으로 감소하는 반면, Type 2에서는 Frontal area가 급격히 증가하는 Branch 구간에서 유속이 급격히 감소한 후 Trunk 구간에서 유속이 다시 증가하는 변화를 보였다. Velocity Spectrum 분석 결과, 모든 지점에 대해 평균한 결과 고주파수 영역에서 -5/3 law를 따르는 것으로 나타나 전체 영역에서 isotropic & homogeneous 난류흐름이 발생함을 확인했다. 난류흐름특성 계산결과, Turbulent kinetic energy(k)를 mean kinetic energy(K)로 무차원화하여 연직분포를 비교했을 때 -k/K는 모두 식생에 근접하며 증가했다. Shear production(P_s)의 계산결과로부터 전단흐름에 의한 난류운동에너지 생성영향분석결과, Type I과 II가 식생경계의 mixing interface 부근에서 급격히 증가하는 분포를 보였으며, 이는 시간평균유속분포에서 분석한 결과와 일치한다. Wake production(P_w)의 연직분포계산결과, P_s와 유사하게 식생경계 부근에서 상승하는 결과가 나타났으며, 이는 식생경계에서 발생하는 Large scale eddy로 인해 발생함을 알 수 있다. 마지막으로 x-방향 난류확산계수로부터 scale factor(α_x)의 연직분포를 계산한 결과, 식생경계부근의 mixing interface에서 증가한 후 식생영역 내에서 감소하는 분포를 나타냈다. z-방향 난류확산 계수의 scale factor(α_z)는 α_x에 비해 작게 계산되었다. 이러한 결과는 오염물질의 연직확산이 식생경계에서 증가한 후 식생 내부에서 감소하여 오염물질, 부유사 등의 축적이 이뤄질 것으로 예상된다. 이는 가지로 인해 식생저항이 증가할 경우 용존성 물질의 혼합이 감소하여 식생의 저장대 효과가 증가함을 의미한다.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1059-1069
• /
• 2023

Dam-break flow occurs when an elevated dam suddenly collapses, resulting in the catastrophic release of rapid and uncontrolled impounded water. This study compares laminar and turbulent closure models for simulating three-dimensional dam-break flows using OpenFOAM. The Reynolds-Averaged Navier-Stokes (RANS) model, specifically the k-ε model, is employed to capture turbulent dissipation. Two scenarios are evaluated based on a laboratory experiment and a modified multi-layered block obstacle scenario. Both models effectively represent dam-break flows, with the turbulent closure model reducing oscillations. However, excessive dissipation in turbulent models can underestimate water surface profiles. Improving numerical schemes and grid resolution enhances flow recreation, particularly near structures and during turbulence. Model stability is more significantly influenced by numerical schemes and grid refinement than the use of turbulence closure. The k-ε model's reliance on time-averaging processes poses challenges in representing dam-break profiles with pronounced discontinuities and unsteadiness. While simulating turbulence models requires extensive computational efforts, the performance improvement compared to laminar models is marginal. To achieve better representation, more advanced turbulence models like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are recommended, necessitating small spatial and time scales. This research provides insights into the applicability of different modeling approaches for simulating dam-break flows, emphasizing the importance of accurate representation near structures and during turbulence.