• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.2
• /
• pp.83-92
• /
• 1998

To evaluate the hydraulic resistance behind bodies in a large scale grid numerical model, a drag stress term which is formulated by the drag force is introduced in the depth-integrated Reynolds equations. And also, the applicability and problems of this model are discussed through various numerical experiments where the analytical solutions exist. In the case of a single body, the error range of velocity difference between analytical and numerical solutions is within $\pm$10% and the wake width behind the body shows a good agreement with the analytical solution. When the drag coefficient and the eddy viscosity are precisely decided, the numerical solutions behind a row of bodies will be efficiently used in real situations.

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

Godunov형 모형을 이용한 홍수모델링에서는 일반적으로 구조적 사각격자나 비구조적 삼각격자가 주로 적용된다. 2차원 수치모형을 이용한 홍수모델링에서 연구유역의 정보가 격자의 노드나 중심에 입력되므로 적용격자의 유형과 생성방법에 따라 모형의 입력자료 오차에 영항을 줄 수 있다. 따라서, 연구유역이 지형 변동성이 심한 지역이거나 흐름형상이나 흐름변동이 심한 구간이라면, 고해상도 격자를 통해 모형의 입력자료 오차를 최소화할 할 수 있다. 본 연구에서는 2가지 유형에 대한 연구를 수행하였다, 첫 번째는 홍수해석을 위한 2차원 모형의 격자형상과 해상도에 따른 홍수위 및 홍수범람범위를 비교·분석하는 연구를 수행하였다. 연구유역은 2000년 10월 29일부터 11월 19일까지 홍수가 발생한 영국의 Severn 강 유역이다. 연구유역의 홍수 모델링을 위한 지형자료는 3m 해상도의 LiDAR(Light Detection And Ranging)를 이용하여 구축하였으며, 격자유형 및 해상도에 따른 2차원 홍수위 및 홍수범람범위를 비교·분석하기 위해서 홍수 발생기간 동안 촬영된 4개(2000년 8월 11, 14, 15, 17일)의 ASAR(Advanced Synthetic Aperture Radar) 영상자료를 활용하였다. 즉, ASAR 영상으로 촬용된 최대범람시기 및 홍수류의 배수기를 활용하여 최대범람범위뿐만 아니라 홍수가 증가하는 시기와 하류단 배수로 인해 홍수가 감소하는 시기를 모두 포함하는 홍수범람범위에 대한 격자유형별 2차원 홍수범람모형의 계산 결과에 대해 비교하였다. 두 번째는 아마존 강 중류유역의 2,500K㎡ 면적에 해당하는 대규모 유역에 대해 SRTM(Shuttle Radar Topography Mission) 지형자료를 이용하여 홍수기와 갈수기에 대해 2차원 모델링을 수행하고 그 결과를 위성자료와 비교하였다.

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

최근 4년간 한반도 근해에서는 53회 이상의 해저지진이 발생하고 있다. 지진해일은 천수효과로 인해 대규모 범람을 일으켜 많은 인명 및 재산피해를 초래할 수 있다. 그 중에서도 원자력발전소의 지진해일에 대한 안전성은 필수적으로 검증되어야한다. 본 연구에서는 Bousssinesq 방정성과 천수방정식을 이용하여 지진해일의 거동을 수치해석하며, 격자 간격과 계산 시간 간격의 세분화에 의한 누적오차를 줄이기 위해 다중격자 연결모형을 이용한다.

• Journal of Korea Water Resources Association
• /
• v.56 no.10
• /
• pp.603-617
• /
• 2023

This study aims to provide a better understanding of the turbulent flow characteristics in swash zone. A double dam-break method is employed to generate the swash zone flow. Comparing with the conventional single dam-break method, a delay between two gate opening can be controlled to reproduce various interactions between uprush and backwash. For numerical simulations, overInterDyMFoam based on OpenFOAM is adopted. Using overInterDyMFoam, interface between two immiscible fluids having different densities (i.e., air and water phases) can be tracked in a moving mesh with multiple layers. Two-dimensional Reynolds-Averaged Navier-Stokes equations are solved with a standard 𝜅-𝜖 turbulence model for momentum and continuity. Numerical model results are validated with laboratory experiment data for the time series of water depth and streamwise velocity. Turbulent kinetic energy distribution is further investigated to identify the turbulence evolution for each flow regime (i.e., uprush, backwash, and swash-swash interaction).

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

오픈 소스 유체역학 소프트웨어인 OpenFOAM은 다양한 유체 흐름에 적용 가능한 프로그램들로 구성되어 있다. 이 중 interFoam은 밀도가 다른 두 유체(i.e., 물, 공기) 간의 경계를 추적하는 기법을 기반으로 한 프로그램으로, 파랑의 거동 모의에 주로 쓰이고 있다. 파생형 프로그램으로는 동격자(dynamic mesh) 및 중첩 격자 기법(overset grid method)을 interFoam에 추가한 overInterDyMFoam이 있다. 두 기법을 사용하면 각각 여러 영역에서 유체흐름과 다중 물체 간의 상호작용을 효율적으로 모의할 수 있다. 본 연구에서는 overInterDyMFoam을 사용하여 두 개 수문의 개방 움직임을 구현하고 생성된 파랑이 포말대(swash zone)에 접근하였을 때의 흐름 특성을 조사하였다. 수치모형실험 결과 수문 개방 속도가 댐 붕괴 파랑 흐름 전파속도에 영향을 미치는 사실을 발견하였다. 또한, 처오름과 처내림의 상호작용에 의한 난류 운동 특성을 조사하기 위해 수문 개방시간 간격을 0초~3초로 설정하였다. 수치모형실험 결과는 수리모형실험의 수면 변동 시계열과 속도 시계열 결과와 비교하여 모형의 정확성이 검증되었다.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1995.10a
• /
• pp.122-128
• /
• 1995

전회의 연구(최와 고, 1994)에서 1／6$^{\circ}$ 경.위도 해상도를 갖는 동아시아 해역의 8개 분조의 조석체계를 수치 시뮬레이션으로 산정하였다. 본 연구에서는 DBDBS database의 5분 격자 해상 수심을 그대로 적용하는 5분 격자 체재의 Gulf of Siam(Thailand), 남중국해, 자바해(Java Sea)를 포함하는 모형을 수립하였다. 서측경계는 Malacca해협의 서측안이며 북측 경계는 대만 해협, 동측경계는 대만남측-루손섬북측, 만다나오남측과 Celebes(Sulawesi)의 북동단의 Manado, Celebes섬의 Makasar부근과 Flores섬의 북단을 연결하는 5개소의 개방 경계를 택하였다. (중략)

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

본 연구에서는 3차원 CFD 모형인 FLOW-3D 모형을 이용하여, 군남 홍수조절지 접근수로의 안정성과 웨어의 월류능력을 검토하였다. 이를 위해서 현지측량자료와 1:5,000도의 수치지도를 이용하여 지형과 기본계획 설계(안)의 여수로를 구성하였다. 지형은 모델링 범위를 충분히 포함하도록 홍수조절지의 상 하류에 걸쳐 작성하였으며, 연직범위는 가능 최대 홍수량(PMF)시 저수지의 최고수위를 초과하는 범위까지 재현하였다. 모델링 범위는 웨어와 같은 구조물 주변에서는 조밀한 격자를 사용하여 정밀한 흐름현상을 재현토록 하였으며, 저수지 부분과 하류하천 부분에 대하여는 느슨한 격자망을 적용하여 최종적으로 2가지 격자망으로 모델링 전체범위에 대하여 격자를 구성하였다. 상류 경계조건으로는 저수지 수위를 경계조건으로 설정하였고, 하류 경계조건으로는 군남 홍수조절지가 일반댐과 달리 상 하류의 표고차가 크지 않아 방류시 하류하천 수위 상승으로 인하여 통수능에 영향을 줄 것으로 판단되므로, 방류량 대 하류하천 형성수위 관계식을 산정하여 방류량에 따른 하류하천의 수위를 하류경계조건으로 설정하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.93-93
• /
• 2016

최근 증가하고 있는 기후변화에 의해 설계빈도를 상회하는 강우의 발생빈도가 증가하고 있으며, 이로 인한 도시유역의 내수범람 피해가 증가하고 있다. 도시유역에서 발생하는 침수 피해의 경우 인적 물적 자원이 집중되어 있는 도시의 특성으로 인해 침수로 인한 직접적 피해 규모가 상당할 뿐만 아니라 침수 발생 후 세균 및 박테리아에 의해 발생하는 수인성 전염병의 유행 등과 같은 2차적 피해 또한 심각한 사회적 비용을 초래할 수 있어 도시유역의 침수 피해를 저감시키기 위한 대책이 절실히 요구되어지고 있다. 도시유역의 침수를 예방하기 위한 대책은 구조적 비구조적 대책으로 구분되어 질 수 있으며 구조적 대책의 경우 침수 피해 예방에 직접적인 효과를 낼 수 있다는 장점이 있으나 대규모 사업예산 및 사업 기간으로 인해 직접적 효과를 보기까지 상대적으로 긴 시간이 필요할 뿐만 아니라 사업 진행 중 대상지역 거주민들의 민원으로 인한 갈등 조정 등으로 인해 사업실행에 어려움을 겪고 있다. 이러한 측면에서 비구조적 대책의 일환인 수치해석을 통한 침수피해 재현 및 침수원인 파악을 통한 구조개선 제안은 구조적 대책의 단점을 보완할 수있는 좋은 대안이 될 수 있다. 도시유역의 경우 비도시유역과 대조적인 차이점으로는 높은 비율의 불투수층, 복잡한 지형, 다수의 인공 구조물 및 배수관망 시스템 등을 들 수 있으며, 침수해석 모형의 정확도를 높이기 위해서는 복잡한 지형의 효율적인 처리가 무엇보다 중요하다. 일반적으로 이용되는 2차원 침수해석 모형들은 직교구조 격자 또는 비구조 격자를 이용하여 지형을 묘사하고 있으며 DEM 자료를 직접 사용하는 직교구조 격자의 경우 지형 데이터 생성이 상대적으로 쉽다는 장점이 있으나 복잡한 지형을 표현하기 위해서는 불필요한 지역까지 높은 해상도를 이용해야 하며 이로 인하여 모의시간이 지나치게 길어지는 문제점이 발생한다. 비구조 격자의 경우 상대적으로 복잡한 도시 유역을 잘 묘사할 수 있다는 장점이 있으나 격자망 생성에 필요한 데이터가 많고 격자망 생성에 지나치게 많은 시간과 노력이 소요된다는 단점이 있다. 따라서 본 연구에서는 위에서 언급한 두 가지 방법의 장점만을 취할 수 있도록 메쉬 세분화 기법을 이용한 2차원 침수해석 모형을 개발 하여 복잡한 지형은 고해상도 메쉬를 이용하여 보다 자세히 묘사하고 상대적으로 복잡하지 않은 지형은 저해상도 메쉬를 이용하여 계산시간을 단축시킬 수 있도록 하였다. 수치해석 기법으로는 엇갈림 격자를 이용하는 Leap-Frog 기법과 유한차분 (Finite difference Method)기법을 이용하였다.

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

혼합층(Mixed layer)은 온도가 일정한 수심층으로, 해수표면에 작용하는 바람의 영향으로 인하여 해수가 위아래로 섞여 형성된다. 이러한 혼합층은 영양염의 순환과 산소의 공급 등과 함께 일차생산량을 결정하는 중요한 요인이 될 수 있으며 혼합층 두께의 변동은 양식 산업에 영향을 미칠 수 있다. 최근에는 기후변화로 인한 해수면 상승 및 해수온 상승 등이 지속되고 있으며, 이러한 현상은 해양생태계의 변화를 초래하여 수산업의 피해를 유발할 수 있다(강원연구원, 2017). 이에 국립수산과학원, 기상청, 국립해양조사원 등 유관기관에서는 정선해양 수온 관측 및 해수순환모델을 이용하여 혼합층의 분석을 수행하고 있으나 격자 구축 및 초기·경계장 설정의 한계가 존재하여 정밀하고 정확한 혼합층 분석에는 어려움이 있다. 이에 본 연구에서는 비정형격자를 사용하여 격자 구축에 제약이 없는 SCHISM (Semi-implicit Cross-scale Hydroscience Integrated System Model)을 이용하여 우리나라 연안해역의 계절변화 및 기후변동성에 따른 혼합층 두께의 변화를 검토하고자 한다. 연구대상지는 서해·동해·남해를 포함한 우리나라 전체 연안 해역(위도: 32°N ~ 39°N, 경도: 124°E ~ 132°E)으로 선정하였으며, 격자크기 100 ~ 3,000 m인 삼각격자로 격자를 구축하였다. 혼합층을 분석하기 위하여 수직격자 층은 50층으로 SZ(Sigma Z coordinate system)좌표계를 사용하였다. 초기·경계장은 FES(Finite Element Solution)2014, HYCOM(Hybrid Coordinate Ocean Model) 및 대기모델 결과를 이용하여 설정하였다. 수치모형 검증을 위하여 수온관측소에서 수심별 측정한 수온 값과 SCHISM 결과 값을 비교하였고, 상대오차가 약 10% 이내로 나타나 모형의 정확도를 확인하였다. 최종적으로 해수면 상승 및 해수온 상승 시나리오를 고려하여 계절별 연안해역의 혼합층 두께의 변화 양상에 대하여 검토하였다. 향후에는 보다 정밀한 대기모델과의 혼합모형 구축 및 다양한 수심 별 관측자료를 활용한다면 실무에서 적용 가능한 혼합층 분석 및 수산업 피해 발생 지역에 대한 피해저감 대책 수립이 가능할 것으로 판단된다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.3
• /
• pp.253-265
• /
• 2007

This study applied the previous results of the NAO model, a tidal correction model, to the open boundary condition for the behavior of tidal elevation and tidal currents in the South Sea. This study used the EFDC model considering the wetting and drying problem and using the $\sigma-coordinate$ as a vertical coordinate and generated two mesh cases of the constant grid size of 2.0 km and the variable grid size of $0.5\sim2.0km$. The numerical results for the tides showed that the predicted results were in quite good agreements with the observational data acquired from the tidal stations of the NORI. The predicted tides were observed to propagate from the east area to the west area in the South Sea. The verification results reveal that the numerical results are more correlated with the measured tidal data as the grid size decreases. The grid size of 2 km results in proper simulation of tidal currents in wide waterway and offshore area whereas the numerical results from the grid size of 0.5 km tend to somewhat underestimate the tidal currents affected by narrow waterway and topography in inner-bay.