• Title/Summary/Keyword: HLLC Riemann 해법

Search Result 8, Processing Time 0.03 seconds

A Numerical Analysis of the Shallow Water Equations Using the HLLL Approximate Riemann Solver (HLLL 근사 Riemann 해법을 이용한 천수방정식의 수치해석)

  • Hwang, Seung-Yong;Lee, Sam-Hee
    • 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 기법의 결과는 매우 고무적이었다.

  • PDF

Development of Two-dimensional Finite Volume Model Applicable to Mixed Meshes (혼합격자의 적용이 가능한 2차원 유한체적모형의 개발)

  • Kim, Byung-Hyun;Han, Kun-Yeon;Son, Ah-Long
    • Journal of Korea Water Resources Association
    • /
    • v.44 no.2
    • /
    • pp.109-123
    • /
    • 2011
  • In this study, 2D finite volume model, which can apply to the mixed meshes that is effective to treat the complicated topography such as a natural river, is developed. To do so, an algorithm for finding the neighbouring cell of a computational cell is introduced, and fluxes are computed using the HLLC approximate Riemann solver at each interface between a computational cell and it's neighbouring cells. Moreover, in order to numerically treat the bed slope which has important effect on the balance between flux gradients and sourte terms, different formula to compute the bed slope for rectangular and triangular mesh are applied. The developed model is applied to analyze dam-break in an experimental channel with $90^{\circ}$ bend and Malpasset dam-break in France. The two cases consist of mixed meshes and the suggested method is validated for the experimental channel and natural channel by comparison with the experimental data, field data and computed results.

An Application of the HLLL Approximate Riemann Solver to the Shallow Water Equations (천수방정식에 대한 HLLL 근사 Riemann 해법의 적용)

  • Hwang, Seung-Yong;Lee, Sam Hee
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.32 no.1B
    • /
    • pp.21-27
    • /
    • 2012
  • The HLLL scheme, proposed by T. Linde, determines all the wave speeds from the initial states because the middle wave is evaluated by the introduction of a generalized entropy function. The scheme is considered a genuine successor to the original HLL scheme because it is completely separated form the Roe's linearization scheme unlike the HLLE scheme and does not rely on the exact solution unlike the HLLC scheme. In this study, a numerical model was configured by the HLLL scheme with the total energy as a generalized entropy function to solve governing equations, which are the one-dimensional shallow water equations without source terms and with an additional conserved variable relating a concentration. Despite the limitations of the first order solutions, results to three cases with the exact solutions were generally accurate. The HLLL scheme appeared to be superior in comparison with the other HLL-type schemes. In particular, the scheme gave fairly accurate results in capturing the front of wetting and drying. However, it revealed shortcomings of more time-consuming calculations compared to the other schemes.

Numerical Analysis of Dam-Break Flow in an Experimental Channel using Cut-Cell Method (분할격자기법을 이용한 실험수조 댐붕괴파의 수치모의)

  • Kim, Hyung-Jun;Kim, Jung-Min;Cho, Yong-Sik
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.29 no.2B
    • /
    • pp.121-129
    • /
    • 2009
  • In this study, dam-break flows are simulated numerically by using an efficient and accurate Cartesian cut-cell mesh system. In the system, most of the computational domain is discretized by the Cartesian mesh, while peculiar grids are done by a cutcell mesh system. The governing equations are then solved by the finite volume method. An HLLC approximate Riemann solver and TVD-WAF method are employed to calculation of advection flux of the shallow-water equations. To validate the numerical model, the model is applied to some problems such as a steady flow convergence on an ideal bed, a steady flow over an irregular bathymetry, and a rectangular tank problem. The present model is finally applied to a simulation of dam-break flow on an experimental channel. The predicted water surface elevations are compared with available laboratory measurements. A very reasonable agreement is observed.

Development of a Numerical Model of Shallow-Water Flow using Cut-cell System (분할격자체계를 이용한 천수흐름 수치모형의 개발)

  • Kim, Hyung-Jun;Lee, Seung-Oh;Cho, Yong-Sik
    • Journal of the Korean Society of Hazard Mitigation
    • /
    • v.8 no.4
    • /
    • pp.91-100
    • /
    • 2008
  • Numerical implementation with a Cartesian cut-cell method is conducted in this study. A Cartesian cut-cell method is an easy and efficient mesh generation methodology for complex geometries. In this method, a background Cartesian grid is employed for most of computational domain and a cut-cell grid is applied for the peculiar grids where the flow characteristics are changed such as solid boundary to enhance the accuracy, applicability and efficiency. Accurate representation of complex geometries can be obtained by using the cut-cell method. The cut-cell grids are constructed with irregular meshes which have various shape and size. Therefore, the finite volume method is applied to numerical discretization on a irregular domain. The HLLC approximate Riemann solver, a Godunov-type finite volume method, is employed to discretize the advection terms in the governing equations. The weighted average flux method applied on the Cartesian cut cell grid for stabilization of the numerical results. To validate the numerical model using the Cartesian cut-cell grids, the model is applied to the rectangular tank problem of which the exact solutions exist. As a comparison of numerical results with the analytical solutions, the numerical scheme well represents flow characteristics such as free surface elevation and velocities in x-and y-directions in a rectangular tank with the Cartesian and cut-cell grids.

Application of mixed mesh for flexible treatment of Topography (지형의 효율적 처리를 위한 혼합격자 적용 기법)

  • Kim, Byung-Hyun;Son, In-Ho;Han, Kun-Yeun
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2010.05a
    • /
    • pp.198-201
    • /
    • 2010
  • 지형이 불규칙한 자연하천에 대해 2차원 격자를 구성할 경우, 사각형 격자만을 사용한다면 지류와 본류의 합류부분에서 격자의 처리가 어려운 문제가 발생할 수 있으며, 삼각형 격자만을 사용하여 지형을 처리한다면 격자수가 많아져 계산시간이 다소 많이 소요되는 어려움이 존재할 수 있다. 혼합격자의 적용이 가능하다면 이러한 어려움은 어느정도 극복할 수 있다. 본 연구에서는 1차정확도 기법인 HLLC 기법을 적용하고, 지형이 복잡한 자연하천에 대한 격자처리의 유연성을 위해 삼각형 및 사각형 격자 그리고 이 두 격자가 혼용된 혼합격자의 적용이 가능한 2차원 유한체적모형을 개발하였다. 그리고 개발모형을 수리모형 실험을 통해 얻어진 실험자료가 존재하는 실험하도 및 실제 자연하천에서의 댐 붕괴에 대해 적용하여 결과를 비교하였다.

  • PDF

Numerical Study on Compressible Multiphase Flow Using Diffuse Interface Method (Diffuse Interface Method를 이용한 압축성 다상 유동에 관한 수치적 연구)

  • Yoo, Young-Lin;Sung, Hong-Gye
    • Journal of Aerospace System Engineering
    • /
    • v.12 no.2
    • /
    • pp.15-22
    • /
    • 2018
  • A compressible multiphase flow was investigated using a DIM consisting of seven equations, including the fifth-order MLP and a modified HLLC Riemann solver to achieve a precise interface structure of liquid and gas. The numerical methods were verified by comparing the flow structures of the high-pressure water and low-pressure air in the shock tube. A 2D air-helium shock-bubble interaction at the incident shock wave condition (Mach number 1.22) was numerically solved and verified using the experimental results.

Handling Method for Flux and Source Terms using Unsplit Scheme (Unsplit 기법을 적용한 흐름율과 생성항의 처리기법)

  • Kim, Byung-Hyun;Han, Kun-Yeon;Kim, Ji-Sung
    • Journal of Korea Water Resources Association
    • /
    • v.42 no.12
    • /
    • pp.1079-1089
    • /
    • 2009
  • The objective of this study is to develop the accurate, robust and high resolution two-dimensional numerical model that solves the computationally difficult hydraulic problems, including the wave front propagation over dry bed and abrupt change in bathymetry. The developed model in this study solves the conservative form of the two-dimensional shallow water equations using an unsplit finite volume scheme and HLLC approximate Riemann solvers to compute the interface fluxes. Bed-slope term is discretized by the divergence theorem in the framework of FVM for application of unsplit scheme. Accurate and stable SGM, in conjunction with the MUSCL which is second-order-accurate both in space and time, is adopted to balance with fluxes and source terms. The exact C-property is shown to be satisfied for balancing the fluxes and the source terms. Since the spurious oscillations in second-order schemes are inherent, an efficient slope limiting technique is used to supply TVD property. The accuracy, conservation property and application of developed model are verified by comparing numerical solution with analytical solution and experimental data through the simulations of one-dimensional dam break flow without bed slope, steady transcritical flow over a hump and two-dimensional dam break flow with a constriction.