• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.43-53
• /
• 2001

본 논문은 충격파 불안정성이 나타나지 않는 충격파 안정적인 수치기법의 개발을 목표로 하고 있다. Roe의 수치기법은 유동의 수치계산에 있어 높은 정확도를 보장하지만 carbuncle 현상과 같은 충격파 불안정성이 나타나는 것으로 알려져 있다. Roe의 수치기법과 HLLE 수치기법의 수치점성을 비교하여 충격파 불안정성의 원인을 살펴보았으며, Roe의 수치기법에 나타나는 반감쇠항에 마하수의 함수인 조절함수 f와 g를 도입하여 충격파 안정성을 획득하였다. 본 논문에서 제안된 수치기법을 다양한 유동문제에 적용하여 수치기법의 충격파 안정성과 정확성을 검증하였다

• 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 기법의 결과는 매우 고무적이었다.

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

• Journal of computational fluids engineering
• /
• v.6 no.2
• /
• pp.40-46
• /
• 2001

Two-dimensional steady flowfields generated by slot injection into supersonic flow are numerically simulated by the integration of Navier-stokes equation with two-equation κ-turbulence model. High-order upwind scheme is used on unstructured adaptive meshes. The numerical results are compared with experimental data in terms of surface static pressure distributions, the length of the upstream separation region, and the height of the Mach surface for steady flowfields with a Mach number of 3.71 and a unit Reynolds number of 5.83×10/sup 6//m.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.99-104
• /
• 2001

This paper deals with the development of shock stable scheme that is free from shock instability. Roe's FDS is known to preserve good accuracy but to suffer from shock instability, i.e. the carbuncle phenomenon. As the first step toward the shock stable scheme, Roe's FDS is compared with HLLE scheme to identify the source of shock instability. Then control function f is introduced into the pressure term in Roe's FDS to cure shock instability. Various numerical tests concerned with shock instability are performed to demonstrate the shock stability of the proposed scheme.