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

댐 제방 등의 붕괴로 인하여 발생하는 급격한 유량의 변화와 흐름영역의 변화로 인한 천이류 및 도수의 발생, 불규칙한 하천단면에서 갈수기 저수기의 흐름해석은 기존의 수치해법의 한계로 인하여 수리모형실험 및 경험식 또는 단면의 단순화 등에 의존하고 있는 실정이다. 본 연구에서는 자연하천에서 비선형 흐름율 계산에 불연속초기조건의 해석해인 Riemann 근사해법을 사용하여 수치적으로 안정되고 정확한 1차원 모형을 개발하고자 한다. 이를 위하여 유한체적법을 사용하였고, 수위와 유량의 계산을 위하여 요구되는 유한체적을 유출입하는 흐름율의 계산에 HLL Riemann 해법을 사용하였으며, MUSCL 기법으로 2차 정확도기법으로 확장하였다. Riemann 해법을 통하여 계산된 비선형의 흐름율과 보존 특성을 만족시켜줄 수 있는 하상 및 하폭변화로 인한 생성항을 처리하는 기법을 제안함으로서 새로운 1차원 수치해석모형을 개발하였다. 개발된 모형의 실제하천의 적용성을 확인하기 위하여 하상과 하폭이 변화하는 부정류 흐름에 적용하여 모형의 적용성 및 정확성을 검증하였다.

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

Riemann 문제는 천수방정식과 같은 쌍곡선형 방정식과 단일한 도약에 의해 불연속인 어떤 점의 좌 우에서 상수인 자료로 구성되는 초기치 문제로서 그 해법은 Godunov 방법과 같이 정확해에 의하면 정확 Riemann 해법, 근사 기법에 의하면 근사 Riemann 해법으로 불린다. 지금까지 이용되는 근사 Riemann 해법으로는 1981년에 P. L. Roe가 제안한 Roe의 선형화 기법과 1983년에 A. Harten, P. D. Lax, 그리고 B. van Leer가 제안한 HLL 기법의 수정 기법들이다. 최대 및 최소 파속만 고려하는 것으로 알려진 HLL 기법은 1988년에 B. Einfeldt의 제안에 의해 두 파속의 결정에서 Roe의 선형화 기법에 따른 고유치와 비교하는 것으로 수정되었다(HLLE 기법). 또한, 1994년에 E. F. Toro 등은 접촉파를 고려하기 위해 선형화된 지배방정식의 정확해로부터 중앙 파속을 고려하는 기법을 제안하였고, 이를 HLLC 기법으로 불렀다. 2002년에 T. Linde는 중앙 파속을 평가하기 위해 일반화된(수학적) 엔트로피 함수를 도입하였으며, van Leer는 이를 HLLL 기법으로 불렀다. 이 기법에서는 접촉파의 평가를 위해 보존변수에 대한 일반화된 엔트로피 함수로부터 중앙 파속이 유도되며, 이것과 특성 속도의 비교를 통해 최대 및 최소 파속이 결정된다. 따라서 이 기법에서는 모든 파속이 초기치로부터 결정되므로 HLLE 기법과 달리 Roe의 선형화 기법과 완전히 결별되고 HLLC 기법과 달리 정확해에 의존되지 않는 점에서 HLLL 기법은 모태인 HLL 기법의 온전한 계승으로 볼 수 있다. HLLL 기법은 여러 분야에 적용된 바 있으나, 수공학 분야에 적용된 사례는 알려진 바 없다. 이는 천수방정식에 대한 (물리적) 엔트로피 함수가 명확하지 않기 때문인 것으로 보인다. 이 연구에서는 보존변수로부터 정의되는 총 에너지를 일반화된 엔트로피 함수로 간주하여 모형을 구성하고, 정확해가 알려진 1차원 문제에 대해 적용성을 검토하였다. 정확해가 알려진 경우에 대해 모의한 결과, 1차 정도 수치해의 한계에도 불구하고, HLLL 기법의 결과는 대체로 정확해와 잘 일치하였으며 그 외의 HLL-형 기법의 그것에 비해 우수한 것으로 나타났다. 특히, 물이 빠져 바닥이 드러나는 상태에 대한 접촉 파속의 추정에서 Riemann 불변량을 이용하는 HLLC 기법에 비해 물이 빠지는 전선을 더 정확하게 포착하는 HLLL 기법의 결과는 매우 고무적이었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.322-322
• /
• 2017

천수 흐름에 대한 수치해석에서 매우 작은 수심의 발생은 해가 불안정해지는 주요 원인 중 하나이며, 경사면이 잠기고 드러나는 그 전선에서 그 현상은 더욱 두드러질 수 있다. 특히, 지배 방정식이 보존형으로 기술되는 경우, 흐름률이나 생성항의 계산에서 수심에 의한 나눗셈이 불가피하므로 보존변수를 정확하게 계산하는 것이 해의 안정성을 도모하기 위한 관건이 된다. 이러한 기대에 부응할 수 있는 수치해법으로 흐름률을 정확한 계산할 수 있는 Riemann 해법을 들 수 있다. 또한, 생성항을 정확하게 계산할 수 있도록 계산 격자를 적절하게 구성하고 그 격자가 완전히 잠기지 않을 경우에 대해 물리적으로 타당하게 처리할 필요가 있다. 이 연구에서는 흐름률의 계산에 근사 Riemann 해법을 적용하여 포물면 지형을 지나는 천수 흐름에 대해 모의하였다. 1981년에 W. C. Thacker는 회전 포물면 위의 천수 문제에 대해 천수방정식의 정확해를 처음으로 유도하였다. 이 문제는 지형의 잠김과 드러남이 다수의 계산 격자에서 지속적으로 이루어지기 때문에 천수흐름의 수치 모의에서 극도로 혹독한 조건의 시험으로 알려져 있다. 회전 포물면 위의 천수 문제에 대해 근사 Riemann 해법에 따른 자료의 재구축 방법, 잠김과 드러남의 처리 등에 대해 검토하였다.

• 한국방재학회:학술대회논문집
• /
• 2009.02b
• /
• pp.152.1-152.1
• /
• 2009

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

댐 제방 붕괴파는 갑작스러운 유량의 증가가 발생하여 불연속적인 흐름특성을 가지는 충격파(shock wave)가 전파되며, 갈수기 저수기에는 중소하천의 상류, 여울과 소에서의 흐름 또는 낙차공이나 보, 댐 여수로 등의 수공구조물에서 부분적인 사류 흐름이 발생된다. 이 때 흐름은 한계수위를 통과하게 되므로 기존 수치해법의 적용에 어려움이 존재한다. 본 연구에서는 실제하천에 적용될 수 있는 1차원 HLL, Roe Riemann 근사해법들을 간단히 소개하고, 시간공간적으로 2차의 고정확도 기법으로 확장하는 방법에 대하여 소개하였다. 각 기법을 정확해가 존재하는 댐붕괴 및 마른하도 전파의 경우에 적용하여 각 기법의 적용성 및 정확성을 비교하였다. 그리고 기존 Lax-Friedrichs 기법과 Lax-Wendroff 기법의 적용결과를 비교하였다. 적용결과 Lax-Friedrichs 기법을 제외한 모든 기법이 정확해와 잘 일치하였으며 특히 HLL 기법을 2차 정확도로 확장한 WAF 기법이 가장 높은 정확도로 계산되었다. 그러나 비선형 생성항이 존재하는 실제하천에 있어서 MUSCL 기법을 이용한 2차 정확도 기법이 합리적일 것으로 판단된다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.3
• /
• pp.183-195
• /
• 2007

Comprehensive mathematical comparison of numerical dissipation vector was made for a compressible and the preconditioned version Roe's Riemann solvers. Choi and Merkle type preconditioning method was selected from the investigation of the convergence characteristics of the various preconditioning methods for the flows over a two-dimensional bump. The investigation suggests a way of migration from a compressible code to a preconditioning code with a minor changes in Eigenvalues while maintaining the same code structure. Von Neumann stability condition and viscous Jacobian were considered additionally to improve the stability and accuracy for the viscous flow analysis. The developed code was validated through the applications to the standard validation problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.365-367
• /
• 2010

In this paper, we presented the exact Riemann solver for the air-water two-phase shock tube problems where the strength of the propagated sock wave is moderately weak. The shock tube has a diaphragm in the middle which separates water medium in the left and air medium in the right. By rupturing the diaphragm, various waves such as rarefaction wave, shock wave and contact discontinuity are propagated into water and air. Both fluids are treated as compressible, with the linearized equations of state. We used the isentropic relations for the air and water assuming a weak shock wave. We solved the shock tube problem considering a high pressure in the water and a low pressure in the air. The numerical results cleary showed a left-traveling rarefaction wave in the water, a right-traveling shock wave in the air, and the right-traveling material interface.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.73-80
• /
• 2010

In this paper, we present the exact Riemann solver for the compressible liquid-gas two-phase shock tube problems. We hereby consider both isentropic and non-isentropic two-phase flows. The shock tube has a diaphragm in the mid-section which separates the liquid medium on the left and the gas medium on the right. By rupturing the diaphragm, various waves are observed on the phasic field variables such as pressure, density, temperature and void fraction in the form of rarefaction wave, shock wave and material interface (contact discontinuity). Both phases are treated as compressible fluids using the linearized equation of state or the stiffened-gas equation of state. We solve several shock tube problems made of a high/low pressure in the liquid and a low/high pressure in the gas. The wave propagations are well resolved by the exact Riemann solutions.

• Journal of Korea Water Resources Association
• /
• v.42 no.4
• /
• pp.271-279
• /
• 2009

The objective of this study is to develop the scheme to apply one-dimensional finite volume method (FVM) to natural river with complex geometry. In the previous study, FVM using the Riemann approximate solver was performed successfully in the various cases of dam-break, flood propagation, etc. with simple and rectangular cross-sections. We introduced the transform the natural into equivalent rectangular cross-sections. As a result of this way, the momentum equation was modified. The accuracy and applicability of newly developed scheme are demonstrated by means of a test example with exact solution, which uses triangular cross-sections. Secondly, this model is applied to natural river with irregular cross-sections and non-uniform lengths between cross-sections. The results shows that the aspect of flood propagation, location and height of hydraulic jump, and numerical solutions of maximum water level are in good agreement with the measured data. Using the developed scheme in this study, existing numerical schemes conducted in simple cross-sections can be directly applied to natural river without complicated numerical treatment.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.761-772
• /
• 2008

The object of this study is to develop the model that solves the numerically difficult problems in hydraulic engineering and to demonstrate the applicability of this model by means of various test examples, such as, verification in the gradually varied unsteady condition, three steady flow problems with the change of bottom slope with exact solution, and frictional bed with analytical solution. The governing equation of this model is the integral form of the Saint-Venant equation satisfying the conservation laws, and finite volume method with the Riemann solver is used. The evaluation of the mass and momentum flux with the HLL Riemann approximate solver is executed. MUSCL-Hancock scheme is used to achieve the second order accuracy in space and time. This study introduce the new and simple technique to discretize the source terms of gravity and hydrostatic pressure force due to longitudinal width variation for the balance of quantity between nonlinear flux and source terms. The results show that the developed model's implementation is accurate, robust and highly stable in various flow conditions with source terms, and this model is reliable for one-dimensional applications in hydraulic engineering.