• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.44-44
• /
• 2011

저수지와 같은 갇혀진 수체는 상류에서 유입되는 오염물질 뿐만 아니라 성층현상에 의해서도 오염될 수 있다. 안정된 성층은 혐기성 조건을 제공하여 바닥에서의 과도한 조류성장을 유발하고 연직순환흐름을 방해한다. 갇혀진 수체에서의 연직순환은 환경적인 문제를 감소시킬 수 있도록 중요한 역할을 한다. 갇혀진 수체에서의 연직순환은 이러한 오염을 줄이는데 중요한 역할을 하는데, 연직순환을 일으키는 인자로는 빛의 입사, 바람, 물의 온도 및 열의 확산 등이 있으며, 그중에서도 가장 중요한 것은 바람의 영향이다. 그러므로 성층화된 흐름에서 바람에 의해 발생하는 연직순환에 대한 수치모형을 개발하고 적용하는 것이 필요하다. 본 연구에서는 수온성층흐름에서의 순환흐름을 해석하기 위하여 유속성분을 계산하고 자유수면 변위와 온도, 염도등과 같은 스칼라양의 해석을 위해 다음과 같은 3단계의 방법을 이용하였다. 첫 번째 단계(정수압 계산단계)에서는 운동량 방정식의 경사항을 음해적으로 해석하고, 두 번째 단계(자유수면 보정단계)에서는 자유수면의 변화를 계산하고 수평방향 유속성분을 계산한다. 예측-수정자 방법(predictor-corrector step method)을 이용하여 자유수면변위와 유속을 구하였다. 마지막으로 세 번째 단계(이송-확산단계)에서는 이송-확산 방정식을 이용하여 스칼라양을 계산하였다. 본 연구에서 제시한 모형의 정확도를 검증하기 위하여 정사각형수조에서 진동하는 자유수면의 해석해와 비교하였고, 성층화된 흐름에서 발생하는 연직순환에 대하여 수치모의를 실시하였다. 그 결과, 본 연구에서 개발된 수치모형이 흐름 내부의 현상을 잘 묘사함을 알 수 있었다.

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

본 연구에서는 일본 전력중앙연구소에서 수행된 wave fission 수리모형실험 자료를 토대로 일차원 FUNWAVE 수치모형을 이용하여 wave fission 현상을 재현하는 수치모의를 수행하였다. FUNWAVE 수치모형은 Boussinesq 방정식을 지배방정식으로 사용하고 있으며 파의 분산효과와 비선형 효과를 고려할 수 있는 수치모형이다. 따라서 wave fission의 주된 발생원인인 분산효과와 비선형효과에 대한 고려를 통해 수치모의 결과는 수리모형실험의 관측치와 상당히 잘 일치함을 확인할 수 있었다. 또한, 본 연구에서는 추가적으로 해수의 흐름이 존재하는 경우를 가정하고 수로 내 일정한 유량의 흐름을 추가하여 wave fission 일차원 수치모의를 수행하였다. 수치모의 결과 파의 진행방향과 반대방향으로 흐름이 존재하는 경우 wave fission으로 인한 수면변위의 크기가 상대적으로 증가함을 확인할 수 있었으며 반대로 파의 진행방향과 동일한 방향으로 흐름이 존재하는 경우 wave fission으로 인한 수면변위의 크기가 상대적으로 감소함을 확인할 수 있었다.

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

하천 수위는 합리적인 수자원의 이용 및 관리를 위해 반드시 필요한 수문 자료이다. 우리나라에서는 수위 측정을 위해 유역 내에 관측소를 설치하여 장비 또는 인력을 통해 수위를 측정하고 있다. 하지만, 많은 관측소를 운영하고 관리하기에는 예산과 인력이 소모되는 한계가 있다. 위성 영상을 통한 시계열 분석은 전지구적 모니터링과 관측 분야에 중요한 역할을 수행할 것으로 기대되고 있으며, 특히 위성 영상자료를 활용한 수자원 분야 연구가 활발히 진행되고 있다. 위성 영상을 활용하여 수면적을 감지하고 수위와 유량을 판별하는 많은 연구가 진행되었지만, 하천 하상의 경사와 단면 형태에 따라 수면적이 변하여 정량적인 수위 추정에는 한계가 존재한다. 본 연구에서는 Sentinel-1의 SAR 영상과 InSAR 기법을 통해 낙동강 유역의 홍수 전후의 하천 수위 변화를 분석하였다. Sentinel-1 IW 모드의 Single Look Complex(SLC) 영상 12장과 ESA 영상 처리 툴인 SNAP을 활용하여 VV(Vertical-Vertical) 데이터의 간섭을 통해 센티미터(cm) 단위지표 변화에 따른 수위 변위를 분석하였다. 위성 영상을 통해 추출한 수위 변위와 계측 수위 및 단면 자료의 정합성을 비교한 결과, 제방과 수체 경계면 식생과 하상 세굴로 인한 오차로 정량적이 수위의 정합성에는 한계가 존재하였지만, 수위의 정량적인 변동성을 확인할 수 있었으며, 수위 변화의 반응속도를 판별할 수 있었다.

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

지진해일은 많은 인명피해를 입힐 수 있는 위험한 자연재해이며, 예를 들어 각각 약 25만명과 약 2만명의 사상자가 발생하였던 2004년 수마트라 지진해일과 2011년 동일본 지진해일 등이 있다. 우리나라 동해안 또한 향후 지진 발생 가능성이 큰 지진공백역이 존재하여 안전한 지역으로 볼 수 없다. 지진해일 방재대책 수립과 관련된 연구는 지속적으로 이루어지고 있지만 지진해일의 발생빈도는 적고 완벽히 대응하는 것은 현실적으로 불가능하다. 따라서 본 연구에서는 지진해일 방재대책의 가장 기본적인 자료로 이용될 수 있는 지진해일 침수예상도를 효율적인 방법으로 제작하는 것을 연구했다. 현재 우리나라의 지진해일 최대 침수예상도는 과거 및 향후 발생가능한 지진해일의 경우에 대한 모든 범람구역이 고려된 보수적인 방법으로 제작되고 있다. 지진원의 위치와 각 매개변수의 특성에 따라 범람구역이 다양하게 나타날 수 있기 때문에 보수적인 최대 침수예상도는 과도한 범람구역이 고려될 수 있다. 따라서 본 연구에서는 보수적인 최대 침수예상도와 비교하여 AI기술과 로직트리 기법을 통해 더 정확한 최대 침수예상도를 제작하는 것을 목표로 한다. 연구방법은 1) 고려된 모든 지진해일 시나리오에 대한 수치해석 2) 입력자료인 지진해일 초기수면 변위 이미지 증강 3) CNN모델을 활용한 초기수면 변위 이미지 분류 4) 분류된 결과의 범람 구역으로 최대 침수예상도를 제작하였다. 향후 연구결과는 지진해일 재해정보도 제작 및 지진해일 침수예측 모델 개발에 활용될 수 있을 것이다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.3
• /
• pp.204-211
• /
• 1993

Boundary element method is applied to simulate nonlinear water waves using Green's identity formula in a numerical wave flume. A system of linear equations is formulated from the governing equation and free surface boundary conditions in order to calculate velocity potential and water surface elevation at each nodal point. The velocity square terms are included in the dynamic free surface boundary condition. The free surface is treated as a moving boundary. the vertical variation of velocity potential being considered in calculating the time derivative of the velocity potential at the free surface. The present method is applied to simulate solitary wave and Stokes 2nd order wave, and shows excellent agreements with their theoretical values.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.59-59
• /
• 2011

파랑의 전파와 변형에 대한 연구에는 수심방향으로 적분한 2차원방정식인 완경사방정식과 Boussinesq 방정식을 기반으로 한 수치모형을 이용한 연구가 최근까지 가장 활발하게 진행되어 오고 있다. 그러나 실제 구조물의 설계에는 2차원 수치모형에서 고려할 수 없는 수심방향 유속에 기인한 정확도의 문제로 인해 구조물의 형상과 재원을 설계하기 위한 정교한 수치모형실험이 어려워 주로 수리모형실험에 의존해 왔다. 수리모형실험은 실제 현상을 가장 잘 재현해낼 수 있어 신뢰성이 매우 높지만 다양한 실험을 수행하기가 어렵고 많은 시간과 비용이 소요되는 단점이 있다. 이에 따라 최근 수심방향으로 완전한 운동방정식인 Navier-Stokes 방정식을 푸는 3차원 수치모형에 대한 연구가 활발히 진행되고 있다. 이론적으로 매우 우수한 모형이긴 하나 정확도 높은 결과를 얻기 위해서는 매우 조밀한 격자를 필요로 하기 때문에 아직까지 막대한 계산시간이 필요하다는 단점이 있으나 컴퓨터 기술이 급격한 속도로 발전하고 있어 Navier-Stokes 방정식 모형의 적용 가능성은 계속 높아지고 있다. 파랑변형을 다루는 수치모형실험을 수행할 때 외부조파를 사용할 경우 구조물이나 지형에 의해 반사되어 나온 파랑이 조파지점에 도달할 때 실험영역으로 재 반사되는 문제가 발생한다. 이를 해결하기 위해 내부조파기법의 개발에 대한 연구가 필수적이었으며, 자유수면변위를 변수로 사용하는 모형의 경우 그 연구가 매우 활발하게 진행되어 왔다. 한편, Navier-Stokes 방정식 모형의 경우 자유수면변위를 변수로 사용하는 2차원 모형에 비해 상대적으로 연구가 미흡하였다. 본 연구에서는 기존의 연직 2차원 Navier-Stokes 방정식 모형에 사용된 연속방정식에 질량 원천항을 추가하는 내부조파기법을 도입하여 3차원 수치모형에서 고립파를 내부조파하고, 급경사에서의 고립파의 처오름 및 처내림 현상을 수리모형 실험결과와 비교 및 분석하였다. 수치모형은 Navier-Stokes 방정식을 엇갈림 격자체계에서 계산하는 동수압 모형으로서, Two-step projection 기법을 사용하는 유한차분모형을 사용하였다. 본 수치모형은 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 계산하는 LES(large-eddy-simulation) 기반의 수치모형으로, 난류 모델링을 위해 Smagorinsky LES 모형을 사용한다. 또한, 압력장의 계산을 위해 Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하였으며, 자유수면 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 수치모형실험이 전체적으로 수리모형실험에서 관측한 파랑의 처오름 및 처내림 현상을 잘 재현하고 있는 것으로 나타났으며, 정량적인 비교를 통해 수치모형의 성능을 검증하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.34 no.4
• /
• pp.93-99
• /
• 2022

Analytical solutions for regular waves generated by bottom wave makers in a 3-dimensional wave basin were derived in this study. Bottom wave makers which have triangular, rectangular and combination of two shapes were adopted. The 3-dimensional velocity potential was derived based on the linear wave theory with the bottom moving boundary condition, kinematic and dynamic free surface boundary conditions in a wave basin. Then, analytical solutions of 3-dimensional particle velocities and free surface displacement were derived from the velocity potential. The solutions showed physically valid results for regular waves generated by bottom wave makers in a wave basin. The analytical solution for obliquely propagating wave generation from bottom wave maker which works like a snake was also derived. Numerical results of the solution agree well with theoretically predicted results.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.36 no.2
• /
• pp.70-79
• /
• 2024

In this study, we introduce a linear spectral method capable of simulating wave generation and transformation caused by a moving bottom in a 3-dimensional space. The governing equations are linear dynamic free-surface boundary conditions and linear kinematic free-surface boundary conditions, which are solved in Fourier space. Solved velocity potential and free-surface displacement should satisfy continuity equation and kinematic bottom boundary condition. For numerical analysis, a 4th order Runge-Kutta method was utilized to analyze the time integral. The results obtained in Fourier space can be converted into velocity potential and free-surface displacement in a real space using inverse Fourier transform. Regular waves generated by various types of moving bottoms were simulated with the linear spectral method. Additionally, obliquely generated regular waves using specified bottom movements were simulated. The results obtained from the spectral method were compared to analytical solutions, showing good agreement between the two.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.6
• /
• pp.368-378
• /
• 2019

Logic trees for probabilistic tsunami hazard assessment include numerous variables to take various uncertainty on earthquake generation into consideration. Results from the hazard assessment vary in different way as more variables are considered in the logic tree. This study is conducted to estimate the effects of various scaling laws and fault parameters on tsunami hazard at the nearshore of Busan. Active fault parameters, such as strike angle, dip angle and asperity, are adjusted in the modelling of tsunami propagation, and the numerical results are used in the sensitivity analysis. The influence of strike angle to tsunami hazard is not as much significant as it is expected, instead, dip angle and asperity show a considerable impact to tsunami hazard assessment. It is shown that the dip angle and the asperity which determine the initial wave form are more important than the strike angle for the assessment of tsunami hazard in the East Sea.

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

본 연구에서는 쇄파대에서 정현파의 쇄파에 대해 수리모형실험과 수치모형실험을 수행하였다. 수치해석 모형에서는 Reynolds 방정식을 지배방정식으로 사용하고 난류해석을 위해 $k-\varepsilon$모델을 적용하였으며, 자유수면변위를 추적하기 위해 VOF기법을 사용하져다. 사면 및 평탄지형상에서 발생하는 쇄파양상을 서로 다르게 설정하기 위해 수심과 입사파의 주기와 파고를 변화시킨다. 발생된 정현파의 파형은 해석해와 잘 일치하였으며, 입사파와 파고계가 설치된 위치에서 측정된 파고비 $H/H_0$는 관측값과 비교해 본 길과 놀은 정확도를 나타내었다.