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

본 연구에서는 일반 좌표계를 사용함으로써 자유수면이 매시간 변화하는 계산격자를 구성해야하는 문제점을 해결하였고, 또한 유한차분법의 단점인 지형변화의 적용성을 보완할 수 있었다. 자유수면과 동수압을 고려하기 위하여 계산단계를 3단계(정수압 계산단계, 동수압 보정단계, 자유수면 보정단계)로 나누어 해석하였고, 개발된 수치모형을 이용하여 실험값이 존재하는 수중방파제를 지나는 비선형파의 전파문제에 대해 수치모의를 실시하였다. 전반적으로 수치모의에 의한 결과는 실험값과 일치하는 경향을 나타내었다.

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

본 연구에서는 자유수면 흐름을 해석할 수 있는 연직방향에 대해 좌표변환된 3차원 동수압 모형을 제시하였다. 제시한 모형은 자유수면파 동수압의 해석을 위하여, 2중 예측-수정(double predictor-corrector)방법을 적용하였다. 본 연구에서는 자유수면과 동수압을 고려하기 위하여 자유수면 보정단계와 동수압 보정단계로 나누어 정확한 동역학적 경계조건을 적용하는 방법을 제시하였고, 제시한 모형을 이용한 수치모의 결과를 검증하기 위하여 타원형 천퇴에 의한 파의 변형에 대한 수치모의를 실시하였다. 전반적으로 수치모의 결과는 실험자료와 일치하였다.

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

본 연구에서는 자유수면의 계산을 위해 동수압을 해석하는 수치모형을 제시하였다. 동수압과 자유수면을 고려하기 위해, 비점성 Navier-Stokes 방정식을 3단계로 나누어 해석하였다. 제1단계에서는 동수압과 자유수면을 전번 시간단계에서 계산된 값으로 대입하여 차분하였으며 이 차분식은 NGMRES(Newton-Generalized Minimal Residual) 방법을 이용해 음해적으로 해석되었다. 이때 계산된 유속장은 연속방정식을 고려하지 않았으므로 각 계산격자에서 질량보존법칙을 만족하지 않을 수도 있다. 제2단계에서 유속과 동수압 보정항으로 이루어진 식을 연속방정식에 대입하여 얻어지는 타원형 방정식인 동수압-포와송 방정식을 해석하므로 여기서 얻어지는 유속은 질량보존법칙을 만족하게 된다. 마지막 3 단계에서는 자유수면과 최종 유속을 계산하였다. 새로 개발된 수치모형을 검증하기위해 정사각형 탱크에서 수면의 자유 진동 문제에 적용한 결과 수치해는 해석해와 잘 일치하였다.

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

트렌치 수로에서 동수압을 고려한 자유수면 흐름을 해석하고 난류 모형의 적용성을 제시하였다. 본 연구에 사용된 지배방정식으로 비정상 상태의 비압축성 유체에 대한 연속방정식과 비점성 Navier-Stokes 방정식을 사용하였다. 난류완결문제를 해결하기 위해서 $\kappa-\varepsilon$방정식을 사용하여서 난류 와점성을 구할 수 있다. 자유수면과 동수압을 고려하기 위하여 3단계로 나누어서 해석하였다. 제 1단계에서는 운동량방정식을 연직방향에 대해 음해적으로 차분하였다. 제2단계에서는 유속과 동수압 보정항으로 이루어진 식을 연속방정식에 대입하여 타원형 방정식인 동수압-포와송 방정식을 해석하여 얻어지는 유속은 질량보존법칙을 만족하게 된다. 마지막으로 자유수면과 최종유속을 보정 및 계산하였다. 본 연구에서 제시한 수치모형을 검증하기 위해서, 트렌치 수로에서 난류의 흐름에 변화를 고려하기 위해 수치모의를 하였다. 전반적으로 수치모의에 의한 결과와 실험 자료가 일치하는 경향을 보였다.

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

댐붕괴파 (dam-break flow)나 지진해일에 의해 발생하는 undular bore와 충격파 (shock) 현상을 동수압 및 분산효과를 고려하여 수치모의를 수행하였다. 완전비선형 Boussinesq-type equations 모형을 이용하여, 동수압 및 분산 효과를 고려하였다. 방정식은 4차 정확도의 유한체적법을 이용하여 해석하였고, 시간적으로도 4차정확도의 기법을 이용하여 고차미분항에 대한 수치분산을 억제하였다. 다양한 경우의 1차원과 2차원 공간에서의 수치모의를 수행하고 검증을 수행하였다. 그 결과, 완전비선형 Boussinesq-type equations 모형은 천수방정식 (shallow water equations) 기반의 모형에서 재현이 불가능한 undular bore 등을 재현 하는 등, 전반적으로 천수방정식 기반의 모형 보다 물리적으로도 타당하고 정량적으로도 실험결과와 잘 일치하는 경향을 보였다. 즉, 댐붕괴파나 지진해일 등에 의한 범람 모의에 있어 동수압과 분산 효과의 중요성이 공학적으로도 매우 중요한 고려사항 임이 나타났다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.691-696
• /
• 2008

In this paper, a three-dimensional non-hydrostatic pressure model for free surface flows using a normalized vertical coordinate system is presented. To strongly couple the free surface and non-hydrostatic pressure in the momentum equations, a double predictor-corrector method is employed. This research is especially focused on implementing the dynamic boundary condition (a zero pressure condition) at the free surface. This boundary condition can be specified accurately with a small modification to existing models. Numerical results with and without this modification clearly show that a precise implementation of the dynamic boundary condition is paramountly important.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.10 no.1
• /
• pp.103-109
• /
• 2010

In this study, a three-dimensional non-hydrostatic pressure model based on a normalized vertical coordinate system for free surface flows is presented. To strongly couple the free surface and non-hydrostatic pressure with the momentum equations, a double predictor-corrector method is employed. The study is especially focused on implementing the dynamic boundary condition (a zero pressure condition) at the free surface with ignoring of the atmospheric pressure. It is shown that the boundary condition can be specified easily with a slight modification to existing models.

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

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.323-331
• /
• 2019

In this study, a computational model approach for the modeling of hydrodynamic pressures acting on radial gates during strong earthquakes is proposed. The use of the dynamic layering method with the Arbitrary Lagrangian Eulerian (ALE) algorithm and the SIMPLE method for simulating free reservoir surface flow in addition to moving boundary interfaces between the fluid domain and a structure due to earthquake excitation are suggested. The verification and validation of the proposed approach are realized by comparisons performed using the renowned formulation derived by the experimental results for vertical and inclined dam surfaces subjected to earthquake excitation. A parameter study for the truncated lengths of the two-dimensional fluid domain demonstrates that twice the water level leads to efficient and converged computational results. Finally, numerical simulations for large radial gates with different curvatures subjected to two strong earthquakes are successfully performed using the suggested computational model.

• Journal of Wetlands Research
• /
• v.16 no.3
• /
• pp.327-336
• /
• 2014

Submerged thin walls are extreme case of submerged rectangular blocks, and could be used for many purposes in rivers or coastal zones, e.g. to tsunami. To understand flow characteristics including flow and pressure fields around a specific submerged thin wall a numerical model was applied which includes computation of hydrodynamic pressure on ${\sigma}$-coordinate. ${\sigma}$-coordinate has strong merits for simulation of subcritical flow over mild-sloped beds. On the other hand ${\sigma}$-coordinate is quite poor to treat sharp structures on the bed. There have been a few trials to incorporate dynamic pressure in ${\sigma}$-coordinate by some researchers. One of the previous approaches includes process of sloving the Poisson equation. However, the above method includes many high-order terms, and requires long cpu for simulation. Another method SOLA was developed by Hirt et al. for computation of dynamic pressure, but it was valid for straight grid system only. Previous SOLA was modified for ${\sigma}$-coordinate for the present purpose and was adopted in a model system, CST3D. Computed flow field shows reasonable behaviour including vorticity is much stronger than the upstream and downstream of the structure. The model was verified to laboratory experiments at a 2DV flume. Time-average flow vectors were measured by using one-dimensional electro-magnetic velocimeter. Computed flow field agrees well with the measured flow field within 10 % error from the speed point of view at 5 profiles. It is thought that the modified SOLA scheme is useful for ${\sigma}$-coordinate system.