• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.2
• /
• pp.108-120
• /
• 1992

An ADI model for linearized shallow water equation is modified using the method of factorization. In order to show its validity. the computational results are compared both with the analytical solution and with those from existing models, for a rectangualr domain with constant and varying amplitudes at the open boundary. It is shown the accuracy of numerical solutions depends on the size of time step. depth and bottom friction. The modified ADI model is shown to be superior to the existing models such as Leendertse (1971). Butler (1980) and Sheng (1983).

• Proceedings of the Korea Contents Association Conference
• /
• 2011.05a
• /
• pp.423-424
• /
• 2011

Lax-Wendroff 기법을 기반의 천수방정식에 대한 유한요소 모형을 낙동강 상류부 단호교 지역에 적용하였다. 모형의 검증을 위하여 이동경계에 대한 Thacker (1981)의 정확해와 침수-노출 기법(wet-dry scheme)을 이용한 본 모형의 수치해를 비교하여 대체로 잘 일치함을 알게 되었다. 또한, U자형 만곡부를 포함한 수로의 흐름에 적용하여, 기존의 수리실험 결과, 수치모의 결과와 비교하고, 모형의 적용 가능성을 확인하였다. 낙동강 상류지역에 위치한 단호교 지역에 본 모형을 적용하고, 이를 기존 유한차분법에 의한 수치해와 비교하였다.

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

일반적으로 유체와 구조물간 상호작용의 수리동역학적 모의에서는 벽경계조건을 통하여 유동에 대한 구조물의 영향이 반영된다. 하지만 도심지에서 발생한 홍수를 예측하려는 경우 이러한 방법으로는 밀집한 구조물들 사이에 형성된 좁은 길들로 인하여 세밀한 격자망을 요하여 큰 계산량을 유발하고 빠른 예측 속도를 기대할 수 없게 한다. 최근 이러한 문제를 극복하기 위해 성긴 격자망에서도 구조물의 유체에 대한 영향을 반영할 수 있도록 하는 방법들이 큰 관심을 받고 있다. 그 중에서도 다공성 천수방정식은 벽경계조건 대신 다공도(posority)의 개념을 이용한 모형으로 도시범람모의에 있어 계산량과 정확도를 가장 적절하게 타협한 모형으로 보고되고 있다. 이러한 흐름에 맞추어 본 연구는 다공도 천수방정식을 해석하는 수치 기법을 개발하였고, 여기에 최근 쌍곡선계 방정식의 수치적 연구들에서 소개된 주요 특징들이 반영되도록 설계하였다. 우선, WENO 기법과 Runge-Kutaa 기법을 통하여 공간과 시간에 대한 고차 정확도를 만족시켰다. 이 때, 재구성 변수와 알고리즘를 새롭게 제시하여 정상흐름조건에 대한 플럭스항과 생성·소멸항간 절단오차에 의한 비물리적인 흐름생성을 억제하였다. 또한, 수치모의 중 음수심의 발생으로 인하여 수치모형이 불안정해지는 현상을 막기 위해, 양-보존성 제한자를 구축하였다. 마지막으로 도심지에서 즐비한 인위적인 구조물에 의해 나타나는 지형적인 불연속의 효과를 적절하게 반영할수 있도록 정상파 재구축의 단계를 구축하여 수치 기법에 반영하였다. 이렇게 구성된 수치기법은 리만문제의 해석해에 기반하여 기존의 주요 연구들의 결과와 비교되었고, 그 결과 본 연구의 방법이 정확성, 수렴성, 안전성의 측면에서 가장 우수함을 수치적으로 증명하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.21 no.1
• /
• pp.39-44
• /
• 2009

To improve the accuracy of numerical simulation of wave trans- formation across the surf zone, nonlinear shoaling effect based on Shuto's empirical formula and breaking mechanism are induced in the elliptic modified mild slope equation. The variations of shoaling coefficient with relative depth and deep water wave steepness are successfully reproduced and show good agreements with Shuto's formula. Breaking experiments show larger wave height distributions than linear model due to nonlinear shoaling but breaking mechanism shows a little bit larger damping in 1/20 beach slope experiment.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1994.02a
• /
• pp.101-106
• /
• 1994

• Journal of Korea Water Resources Association
• /
• v.36 no.5
• /
• pp.777-785
• /
• 2003

A numerical model for the solution of two-dimensional free surface flow is developed on unstructured grid. By using fractional step method, the two-dimensional shallow water equations (SWE) are treated as two one-dimensional problems. Thus, it is possible to simulate computational hydraulic problems with higher computational efficiency. The one-dimensional problems are solved using upwind TVD version of second-order Weighted Averaged Flux (WAF) scheme with HLLC approximate Riemann solver. The numerical oscillations which are common with second-order numerical scheme are controlled by exploiting WAF flux limiter, Some idealized test problems are solved using this model and very accurate and stable solutions are obtained. It can be concluded as an efficient implement for the computation of SWE including dam break problems that concerning discontinuities, subcritical and supercritical flows and complex domain.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.4
• /
• pp.229-239
• /
• 2019

A model of landslides is developed using the shallow water equations to simulate time-dependent performance of landslides. The shallow water equations are derived using the (b, s) coordinate system which can be applied in both river and ocean. The finite volume scheme employing the HLL approximate Riemann solver and the total variation diminishing (TVD) limiter is applied to deal with the numerical discontinuities occurring in landslides. For dam-break water flow and debris flow, numerical results are compared with analytical solutions and experimental data and good agreements are observed. The developed landslide model is successfully applied to predict subaerial and submarine landslides. It is found that the subaerial landslide propagates faster than the submarine landslide and the speed of propagation becomes faster with steeper bottom slope and less bottom roughness.

• Journal of Korea Water Resources Association
• /
• v.38 no.6 s.155
• /
• pp.461-469
• /
• 2005

A second-order upwind scheme is used to investigate the run-up heights of tsunamis in the East Sea and the predicted results are compared with the field data and results of a first-order upwind scheme. In the numerical model, the governing equations solved by the finite difference scheme are the linear shallow-water equations in deep water and nonlinear shallow-water equations in shallow water. The target events are 1983 Central East Sea Tsunami and 1993 Hokkaido Nansei Oki Tsunami. The predicted results represent reasonably well the run-up heights of tsunamis in the East Sea. And, the results of simulation are used for the design of inundation map.

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

본 연구에서는 2차원 비선형 천수모형과 수심평균된 스칼라 이송모형을 해석하는 수치모형에 대해 기술하였다. 수치모형의 정확성을 보장함과 동시에 안정성을 높이기 위해 유한체적법, 플럭스 재구성 및 minmod 제한자를 사용하였다. 비선형 천수방정식의 이송항과 바닥 경사항은 계산된 수심의 양수 보존과 흐름의 정상 상태를 보장하기 위한 second order well-balanced positivity preserving central-upwind method를 이용하여 수치적으로 이산화되었다. 마찬가지로, 이송-확산 방정식 내 이송항은 동일한 2차 풍상차분법을 통해 수치적으로 풀이하였다. 격자점 경계면에서의 불연속으로 인한 수치진동을 방지하기 위해 이송항의 계산에 포함된 보존항의 차이로 인해 발생하는 스칼라의 수치확산을 최소화하기 위해 무차원의 비소산함수를 도입하였다. 또한, 확산항은 유한차분법을 이용하여 이산화하였다. 제안된 수치모형은 시간미분항의 계산을 위해 오일러 기법을 적용하여 계산된 수심 및 스칼라의 양수 보존여부와 함께 정지된 흐름의 정상 상태의 보존여부를 확인하였다. 제안된 수치모형의 해석 정확성을 평가하기 위해 1, 2차원 공간 내 다양한 흐름 조건에서의 해석해를 이용한 3개의 벤치마크 테스트를 수행하였다. 평균 제곱근 오차(Root Mean Squared Error, RMSE)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 2002.08a
• /
• pp.287-293
• /
• 2002

조석을 수치모의하기 위해 North carolina 대학의 R.A. Luettich와 Notre Dame 대학의 J.J. Westerink가 개발한 ADvanced CIRCulation model for oceanic, coastal and estuarine waters(ADCIRC) 유한요소모형 (Luettich, et al., 1992)은 수심 적분된 2 차원 모형(2DDI)과 3차원 모형(3DL)으로 구성되어 있는데, 그 중 2차원 ADCIRC 유한요소모형은 천수방정식에서 연속방정식과 운동방정식을 합성하여 수면변위에 대해 하나의 식으로 표현한 파동방정식(wave equation)을 지배방정식으로 사용하고 있다. (중략)