• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.537-537
• /
• 2012

천수방정식의 수치모형은 하천의 유량예측, 홍수범람해석, 해일 모의 등에 널리 사용되고 있고, 그 결과는 수자원 관리, 재난 대책 등 정책적인 의사결정에 있어 유용한 자료로 활용되고 있다. 이처럼 천수방정식의 수치모형은 연구목적뿐만 아니라 실생활에 있어서도 큰 영향을 미치고 있으며, 이에 보다 정확하고 효율적인 수치모형의 구축이 수리/수자원/방재분야에서 중요한 영역이 되고 있다. 본 연구는 정확하고 안정적인 수치모의를 위해 천수방정식의 수치모형에 MLP(Multi dimensional Limiting Process)기법을 적용한 후 다차원 모의 시 MLP의 수치 진동 제어 성능을 검증하고자 하였다. MLP기법은 다차원에서 수치진동을 억제할 수 있도록 개발된 기법으로, 기존 TVD 제어자(limiter)과 MLP의 차이점은 기존 제어자들이 흐름이 발생하는 셀 경계면에서 재구성된 값이 Maximum Principle을 만족시키도록 제어자를 유도하는데 반해, MLP는 셀 절점에서 Maximum Principle을 만족시키도록 제어자를 유도한다는데 있다. MLP기법은 압축성 유체를 표현하는 2, 3차원 오일러 방정식에 적용되어 기존의 제어자들에 비해 안정적이며 정확한 수치모의를 가능하게 하는 것이 검증되었다. 하지만 천수방정식에 적용된 예는 없으며, 이에 본 연구는 천수방정식에 MLP를 적용하고 천수방정식 수치모형 검증에 주로 사용되는 수치모의를 통해 MLP의 수치 진동 제어 성능을 검증하였다. 모의 결과, MLP는 2차원 천수방정식에 있어서도 기존의 제어자들과 비교하여 수치진동을 보다 잘 제어하는 것으로 판단된다. MLP의 사용으로 인해 불연속면 근처에서 정확도가 향상되었고 수치진동이 발생하지 않아 보다 안정적인 모의가 가능하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.2B
• /
• pp.123-130
• /
• 2012

MLP (Multi dimensional Limiting Process) is implemented to simulate shallow water flows, and its performance over conventional TVD limiters in multidimensional flows is verified through several numerical simulations. MLP was developed to control oscillations for multi-dimensional compressible flows and proved to improve accuracy, efficiency and robustness in compressible flows. In this study, we applies MLP to modeling shallow water equations(SWEs) given that the SWEs are amenable to be solved using the large range of numerical methods developed to deal with compressible flows and MLP has been yet used for SWEs. Simulation results through the benchmark tests show that MLP has favorable features such as numerical oscillation control and convergence behaviors comparable to the conventional limiters. Both numerical accuracy and stability are improved in multi-dimensional discontinuous flows.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.9
• /
• pp.12-26
• /
• 2004

A new treatment of cell-interface flux in AUSM-type methods is introduced to reduce the numerical dissipation. Through analysis of TVD limiters, a criterion for the more accurate prediction of cell-interface state is found out and M-AUSMPW+ is developed by determining the transferred property newly and appropriately within the criterion. The superiority of M-AUSMPW+ is clearly revealed in multi-dimensional flow problems. It can eliminate numerical dissipation effectively in a non-flow aligned grid system. As a result, M-AUSMPW+ is shown to be much more accurate and effective than other previous schemes in multi-dimensional problems. Through a stationary contact discontinuity, a vortex flow, a shock wave/boundary layer interactions and viscous shock tube problems, it is verified that accuracy of M-AUSMPW+ is improved.

• Journal of Korea Water Resources Association
• /
• v.36 no.4
• /
• pp.597-608
• /
• 2003

A transcritical flow occurs when the width and slope of a channel are varying abruptly. In this study, the transcritical flow in a two-dimensional open channel is analyzed by using the shallow-water equations. A weighted average flux scheme that has flux limiter with a total variation diminishing condition is introduced for a second-order accuracy in time and space, and non- spurious oscillations at discontinuous points. A HLLC method with three wane speeds is employed to calculate the Riemann problem. To overcome difficulties resulting from variation of channel sections in a two-dimensional analysis of transcritical flow, the numerical model is developed based on a generalized grid system.

• Journal of Korea Water Resources Association
• /
• v.36 no.5
• /
• pp.777-785
• /
• 2003

A numerical model for the solution of two-dimensional free surface flow is developed on unstructured grid. By using fractional step method, the two-dimensional shallow water equations (SWE) are treated as two one-dimensional problems. Thus, it is possible to simulate computational hydraulic problems with higher computational efficiency. The one-dimensional problems are solved using upwind TVD version of second-order Weighted Averaged Flux (WAF) scheme with HLLC approximate Riemann solver. The numerical oscillations which are common with second-order numerical scheme are controlled by exploiting WAF flux limiter, Some idealized test problems are solved using this model and very accurate and stable solutions are obtained. It can be concluded as an efficient implement for the computation of SWE including dam break problems that concerning discontinuities, subcritical and supercritical flows and complex domain.

• Journal of Korea Water Resources Association
• /
• v.43 no.1
• /
• pp.105-117
• /
• 2010

This paper presents a new and simple technique to perform MUSCL reconstruction for solving 2D shallow water equations. The modified MUSCL scheme uses weighted area ratio to apply non-uniform grid in stead of the previous method that equally distributed the difference of conservation variables to each interface. The suggested method can physically reconstruct conservation variables in case of uniform grid as well as non-uniform grid. In this study, Unsplit scheme applicable to unstructured grid is used and efficient slope limiter of TVD scheme is used to control numerical oscillation which can be occurred in modified MUSCL scheme. For accurate and efficient treatment of bed slope term, the modified MUSCL scheme is coupled with the surface gradient method. The finite volume model applied to suggested scheme is verified through a comparison between numerical solution and laboratory measurements data such as the simulations of isolated building test case and Bellos's dam break test case.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.3
• /
• pp.69-78
• /
• 2007

In this study, a time-dependent aspect of an embankment failure is considered to simulate a flood inundation map and calculate overtopping discharge induced by an embankment failure. A numerical model has been developed by solving the two dimensional nonlinear shallow water equations with a finite volume method on unstructured grids. To analyze a Riemann problem, the HLLC approximate Riemann solver and the Weighted Averaged Flux method are employed by using a TVD limiter and the source term treatment is also employed by using the operator splitting method. Firstly, the numerical model is applied to a dam break problem and a sloping seawall. Obtained numerical results show good agreements with experimental data. Secondly, the model is applied to a flow induced by an embankment failure by assuming that the width and elevation of embankment are varied with time-dependent functions. As a result of the comparison with each numerical overtopping discharge, established flood inundation discharges in the previous studies are overestimated than the result of the present numerical model.