• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.6
• /
• pp.326-334
• /
• 2017

We obtain height of waves overtopping on a porous breakwater using both the one-layer and two-layer Boussinesq equations. The one-layer Boussinesq equations of Lee et al. (2014) are used and the two-layer Boussinesq equations are derived following Cruz et al. (1997). For solitary waves overtopping on a porous breakwater, we find through numerical experiments that the height of waves overtopping on a low-crested breakwater (obtained by the Navier-Stokes equations) are smaller than the height of waves passing through a high-crest breakwater (obtained by the one-layer Boussinesq equations) and larger than the height of waves passing through a submerged breakwater (obtained by the two-layer Boussinesq equations). As the wave nonlinearity becomes smaller or the porous breakwater width becomes narrower, the heights of transmitting waves obtained by the one-layer and two-layer Boussinesq equations become closer to the height of overtopping waves obtained by the Navier-Stokes equations.

• 한국방재학회:학술대회논문집
• /
• 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)와 비교하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.2
• /
• pp.41-49
• /
• 2019

We approximately obtain heights of cnoidal waves overtopping on a porous breakwater using both the one-layer Boussinesq equations (Vu et al., 2018) and the two-layer Boussinesq equations (Huynh et al., 2017). For cnoidal waves overtopping on a porous breakwater, we find through numerical experiments that the heights of cnoidal waves overtopping on a low-crested breakwater (obtained by the Navier-Stokes equations) are smaller than the heights of waves passing through a high-crested breakwater (obtained by the one-layer Boussinesq equations) and larger than the heights of waves passing through a submerged breakwater (obtained by the two-layer Boussinesq equations). As the cnoidal wave nonlinearity becomes smaller or the porous breakwater width becomes narrower, the heights of transmitting waves obtained by the one-layer and two-layer Boussinesq equations become closer to the height of overtopping waves obtained by the Navier-Stokes equations.

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

본 연구에서는 GPU 가속화 기반의 Boussinesq 모형인 Celeris Advecnt에 수심 적분된 2차원 이송-확산방정식을 추가하여 인터렉티브 시스템 기반의 추적자 이동 모형을 개발하였다. Celeris Advent는 최초로 개발된 인터렉티브 시스템을 갖춘 Boussinesq 모형으로, 시뮬레이션 중에 사용자가 모형의 파라미터뿐 아니라 모델 도메인 내 수위 및 수심을 바꿀 수 있다. 이를 통해 사용자는 모의가 진행되는 도중에 모델의 안정성 및 효율성을 위해 시간 간격을 조정할 수 있을 뿐 아니라 방파제 설치 등과 같은 지형 변화를 고려하기 위해 도메인 내 격자별 수심을 조정할 수 있다. 본 연구에서는 연안에서의 추적자 이동 모의를 위해 Boussinesq 방정식과 더불어 이송-확산방정식을 풀이하는 추적자 이동 모형을 개발하였다. 추적자의 확산항의 경우 분자 자체의 확산과 더불어 쇄파에 따른 난류 확산을 고려하였다. 난류 확산계수는 슈미트 수를 1로 두어 와동점성계수와 동일하게 두었으며, 와동점성은 단순화된 형태의 쇄파모형을 고려하여 계산하였다. 쇄파모형의 고려로 인해 이송-확산방정식과 더불어 운동량 방정식에서도 쇄파에 따른 운동량 소산이 고려되었다. 마지막으로, 추적자 농도에 대한 인터렉티브 시스템을 추가하여, 모델 구동 중에도 사용자가 수심적분된 추적자 농도를 조정할 수 있도록 하였다. 기수행된 2개의 수리실험 조건과 관측값을 이용하여 벤치마크 테스트를 수행하였으며, 관측값과 대체로 일치하는 것을 확인하였다.

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

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.446-456
• /
• 2007

Numerical analyses of stem waves, the interaction between incident and reflected waves of obliquely incident regular waves along a vertical wall in a constant water depth, are presented. For the numerical model of the analysis, the two-layer Boussinesq equations developed by Lynett and Liu(2004a,b) are employed. Numerical results are compared with both laboratory measurements and those obtained using parabolic approximation model. The overall comparisons between the results from the two numerical models and the experiments are good. However, the two-layer Boussinesq model is more accurate than the parabolic approximation model as the angle of incident waves increases. In particular, the higher harmonic generation due to the wave nonlinearity is captured only in the Boussinesq model.

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

Madsen et al. (2002)이 제안한 일차원 고차 Boussinesq 방정식에 대하여 불연속갤러킨 유한요소법(Discontinuous Galerkin Finite Element Method)을 적용하였다. 연속적인 Boussinesq 방정식에서 각 요소경계에 불연속을 허용할 수 있도록 공간차분하고, 시간방향으로 4차 Runge-Kutta 시간적분법, 각 요소사이에는 Lax-Friedrichs 수치흐름률을 사용하였다. 계산영역의 양쪽에 불필요한 파랑의 반사를 억제하도록 흡수층을 설치하였으며, 영역 내부에서 조파할 수 있도록 하였다. Luth et al.(1994)의 수중잠제 실험에 적용하여 관측값과 잘 일치함을 확인하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.4
• /
• pp.274-281
• /
• 2002

Numerical analysis for the Bragg reflection due to a sinusoidally and a doubly-sinusoidally varying seabeds was performed by using a couple of ordinary differential equations derived from the Boussinesq equations. Incident waves are a train of cnoidal waves. The effects of the dispersion and shape of seabed were investigated. It is shown that the reflection of a sinusoidally varying seabed is enhanced by increasing the dispersion and the amplitude of a seabed. The reflection of waves over a doubly-sinusoidally varying seabed can also be enhanced by increasing the amplitude of seabed decreasing the difference of wave numbers of seabed components.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.1
• /
• pp.21-31
• /
• 2005

In this study, derivation is made of a one-grid source function for the extended Boussinesq equations of Nwogu (1993) in order to generate nonlinear waves internally. The energy velocity approach used in the line source method is verified analytically by the fractional step splitting method. The source function method is verified by generating accurately nonlinear waves as well as linear waves for horizontally one-dimensional cases. It is found that numerical solutions by the source function method are the same as those by the line source method.

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

불규칙파를 사용하여 설계 자료로 이용하기 위해서는 설계해역에서 불규칙파의 파랑변형을 예측할 수 있는 수치모형의 개발이 선행되어야 한다. 비선형 불규칙파의 거동을 해석할 수 있는 Boussinesq 방정식은 상대파고인 $\alpha$/h($\alpha$는 수면의 진폭, h는 수심임)를 비선형의 매개변수로 하고 상대수심인 kh(k는 파수임)를 분산성의 매개변수로 하여 섭동법을 사용하여 유도된다. Boussinesq 식은 수심이 일정한 경우에 Boussinesq(1872)가 비선형 항을 O($\alpha$/h,(kh)$^2$)까지 포함하여 처음으로 개발하였고 수심의 변화가 완만한 경우에 Peregrine(1967)이 개발하였다. (중략)