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

지진해일은 진행속도가 빠르고 파장이 길며 파형의 변화 없이 먼 거리를 진행 할 수 있어 주변지역은 물론 멀리 떨어진 지역에도 심한 범람피해를 야기시킨다. 지진해일의 일반적인 특징으로 장파와 단파가 합성되어 있고 먼 거리를 전파할 경우 분산효과의 역할이 중요하게 된다. 특히 우리나라의 동해안에 영향을 주는 지진해일은 단주기파 성분이 강하고 파장에 비해 먼 거리를 전파하기에 분산을 고려하는 선형 Boussinesq 방정식을 지배방정식으로 사용하는 것이 바람직하다. 하지만 지금까지의 지진해일 전파모의를 위한 모형은 선형 Boussinesq 방정식의 복잡한 계산과 계산시간이 길다는 단점 때문에 선형 천수방정식을 지배방정식으로 사용하고 분산효과는 수치분산을 이용하여 고려해왔다. 지진해일 해석 시 일반적으로 사용되어 오던 기존의 leap-frog 유한차분 모형(Imamura et al., 1988; 조용식, 1996)은 지배방정식으로 선형 천수방정식을 사용하고 파의 분산효과는 수치분산을 이용하여 고려하므로 정해진 시간 간격에 대해 수심에 따라 격자 간격을 적절히 선택해야 하는데 수심이 복잡하게 변하는 경우 격자간격 조정이 불가능하여 분산효과를 정도 높게 고려할 수 없다. 이 문제점을 해결하기 위하여 윤성범 등(2004)은 파동방정식의 인위적인 분산항을 이용하여 Boussinesq 방정식의 분산효과를 고려할 수 있는 능동적인 분산보정기법을 제안하였고 Cho et al.(2007)는 일정한 수심에서 수치적인 분산오차가 Boussinesq 방정식의 물리적인 분산항을 대체하도록 수심, 격자 간격 및 계산 시간 간격 사이의 관계식을 유도하고 Boussinesq 방정식의 분산항과 일치하는 수치분산을 이용하여 실용적인 분산보정기법을 개발하였다. 이에 Ahn(2010)은 현재 컴퓨터의 계산 능력이 향상되어 선형 Boussinesq 방정식을 직접 차분하여 계산하는데 무리가 없다고 판단하여 선형 Boussinesq 방정식을 직접 차분한 모형을 개발하였다. 본 연구에서는 기존의 원해 지진해일 전파모의에 이용되어왔던 선형 천수방정식에 수치분산을 고려한 모형 대신 선형 Boussinesq 방정식의 유한차분 모형을 제안하였으며 기존의 선형 Boussinesq 방정식 모형의 격자와 수심간의 제약을 없애기 위해 차분 기법을 달리 한 2차 정확도의 유한차분 모형을 제안하였다. 검증을 위하여 선형 Boussinesq 방정식의 해석해(Carrier, 1991)와 비교하였다.

• 한국방재학회:학술대회논문집
• /
• 2010.02a
• /
• pp.107.2-107.2
• /
• 2010

본 연구에서는 지진해일의 전파 과정을 모의함에 있어 선형 천수방정식의 수치분산을 이용하는 기법이 아닌 선형 Boussinesq 방정식을 직접 차분하는 유한차분기법을 제안하였다. 기법의 정확성을 검증하기 위하여 Gauss 분포의 초기 자유수면변위를 갖는 문제에 착용하여 선형 Boussinesq 방정식의 해석해와 비교하였다. 그 결과 기존의 선형 천수방정식을 차분화한 수치모형에 비하여 정확한 결과를 제공하였고 분산보정기법을 이용한 수치모형과 동일한 정확도를 보였으나 본 수치모형을 이용했을 때 계산 효율이 개선되었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.424-428
• /
• 2005

본 연구에서는 임의로 변화하는 지형상에 적용시에 보존 특성이 성립하는 쌍곡선형 천수방정식 해석 기법을 개발하였다. 일반적으로 쌍곡선형의 천수방정식은 상류와 사류를 쉽고 정확하게 해석할 수 있고, 또한 Euler 방정식 해석기법을 이용한 다양한 해석기법이 개발되어 있다는 장점을 지니고 있다. 그러나 바닥지형이 변화하는 경우, 생성항과 플럭스항 사이에 수치적 해석기법 차이에서 발생하는 수치적 불균형이 발생하여 수치모형의 적용성이 현저하게 저하된다. 따라서 본 연구에서는 이와 같은 현상을 개선하기 위하여, 기존의 표면경사법을 개선한 기법을 제시하였다. MUSCL-Hancock 기법과 HLLC 근사 Riemann 기법을 이용하였으며, 플럭스항과 수치적 균형을 이루기 위한 이산화기법을 제안하였다. 모형의 검증을 위하여 정상류 상태의 상류와 사류 해석을 수행하였고, 마른바닥에서의 댐붕괴파와 수직한 지형 변화를 갖는 수로상의 서지의 진행 등과 같은 부정류에 대하여 적용하였다. 적용결과, 매우 정확하고 수치적으로 안정된 계산결과를 얻었다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.11 no.4
• /
• pp.197-200
• /
• 1999

An extension of the modified leap-frog scheme is made to solve the nonlinear shallow-water equations. In the extended model. the physical dispersion of the Boussinesq equations is replaced by the numerical dispersion resulted from the leap-frog finite difference scheme. The model is used to simulate propagations of a solitary wave over a constant water depth and a linearly varying water depth. Obtained numerical results are compared with available analytical and other numerical solutions. A reasonable agreement is observed.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 2003.08a
• /
• pp.240-244
• /
• 2003

파랑의 변형 가운데 천수, 굴절, 회절, 반사를 예측하는 수학적 모형은 크게 두 가지 유형으로 나눌 수 있는데, 첫 번째로 파형경사인 ha(k：파수. $\alpha$：진폭)를 비선형의 매개변수로 하는 Stokes 파랑식이 있고, 두 번째로 상대파고인 $\alpha$/h를 비선형의 매개변수로 하고 상대수심인 kh를 분산성의 매개변수로 하는 천수방정식(Shallow water equation)이 있다. 파랑의 변형 가운데 천수, 굴절만을 예측하고 회절, 반사를 예측하지 못하는 수학적 모형으로는 에너지 이송방정식이 있다. (중략)

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.52-56
• /
• 2010

본 연구에서는 지진해일의 전파 과정을 모의함에 있어 선형 천수방정식의 수치분산을 이용하는 기법이 아닌 선형 Boussinesq 방정식을 직접 차분하는 유한차분기법을 제안하였다. 지배방적식과 차분식의 일치성을 해석하기 위해 이산화 오차를 확인하고, 수치해의 안정적 수렴여부를 판단하기 위해 Von neumann 안정성 해석을 수행하였다. 또한 기법의 정확성을 검증하기 위하여 Gauss 분포의 초기 자유수면변위를 갖는 문제에 적용하여 선형 Boussinesq 방정식의 해석해와 비교하였다. 그 결과 기존의 선형 천수방정식을 차분화한 수치모형에 비하여 정확한 결과를 제공하였고 분산보정기법을 이용한 수치모형과 동일한 정확도를 보였으나 본 수치모형을 이용했을 때 비교적 넓은 범위의 조건에서 정확도 높은 결과를 제공하였다.

• 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.

• 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).

• 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)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.