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

The present paper deals with the continuous work of extending multi-dimensional limiting process (MLP), which has been quite successfully proposed on two- and three-dimensional structured grids, onto the unstructured grids. The basic idea of the present limiting strategy is to control the distribution of both cell-centered and cell-vertex physical properties to mimic a multi-dimensional nature of flow physics, which can be formulated as so called the MLP condition. The MLP condition can guarantee a high-order spatial accuracy without yielding spurious oscillations. Recently, MLP slope limiter was proposed based on the MUSCL-type reconstruction in two-dimensional case and it can be readily extended to three-dimensional case. Through various numerical analyses and extensive computations, it is observed that the proposed limiters are quite effective in controlling numerical oscillations and very accurate in capturing both discontinuous and continuous multi-dimensional flow features on 3-D tetrahedral grids.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5B
• /
• pp.429-439
• /
• 2009

When Catastrophic extreme flood occurs due to dam break, the response time for flood warning is much shorter than for natural floods. Numerical models can be powerful tools to predict behaviors in flood wave propagation and to provide the information about the flooded area, wave front arrival time and water depth and so on. But flood wave propagation due to dam break can be a process of difficult mathematical characterization since the flood wave includes discontinuous flow and dry bed propagation. Nevertheless, a lot of numerical models using finite volume method have been recently developed to simulate flood inundation due to dam break. As Finite volume methods are based on the integral form of the conservation equations, finite volume model can easily capture discontinuous flows and shock wave. In this study the numerical model using Riemann approximate solvers and finite volume method applied to the conservative form for two-dimensional shallow water equation was developed. The MUSCL scheme with surface gradient method for reconstruction of conservation variables in continuity and momentum equations is used in the predictor-corrector procedure and the scheme is second order accurate both in space and time. The developed finite volume model is applied to 2D partial dam break flows and dam break flows with triangular bump and validated by comparing numerical solution with laboratory measurements data and other researcher's data.

• Journal of Korea Water Resources Association
• /
• v.47 no.11
• /
• pp.1061-1066
• /
• 2014

A stable second-order finite volume method was proposed to predict sediment transport under rapidly varied flow conditions such as transcritical flow. For the use under unsteady flow conditions, a sediment transport model was coupled with shallow water equations. HLLC approximate Riemann solver based on a monotone upstream-centered schemes for conservation laws (MUSCL) reconstruction was used for the computation of the flux terms. From the comparisons of dam break flow experiments on erodible beds in one- and two-dimensional channels, good agreements were obtained when proper parameters were provided. Lastly, dam surface erosion problem by overtopped water was simulated. Overall, the numerical solutions showed reasonable results, which demonstrated that the proposed numerical scheme could provide stable and physical results in the cases of subcritical and supercritical flow conditions.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.109-109
• /
• 2019

본 연구에서는 강우-유출 과정 모의를 위한 GPU 기반 확산파 모형을 개발하였다. 확산파 방정식을 풀기위한 수치기법으로는 유한체적법을 이용하였으며, van Leer TVD limiter를 적용한 MUSCL 기법을 이용하여 각 셀의 인터페이스의 물리적 성질을 재구성하여 구하였다. 또한, 침투를 고려하기 위하여 Horton 침투 모형을 이용하였다. 개발된 모형을 이용하여 1D single overland plane과 2D V-shaped overland에서 강우-유출 과정을 모의실험을 하였으며, 각각 해석해와 dynamic wave model을 이용하여 계산된 수치 결과와 비교하여 본 모형의 정확성을 검증하였다. 또한, 1D와 2D의 기복이 심한 지형에 적용하여 강우-유출과정이 본 모형을 통하여 물리적으로 타당한 해석이 가능함을 검증하였다. 마지막으로 복잡한 실제 지형에 적용하였으며, 측정값과의 비교를 통하여 실제 유역에서의 확산파 모형의 적정성을 검증하였다. 또한, 본 연구에서는 NVIDIA사의 GPU인 Geforce GTX 1050과 GPU의 병렬 연산 처리 능력을 활용할 수 있는 NVIDIA사의 CUDA-Fortran을 이용하여 GPU 기반 확산파 모형을 개발하였다. PC windows에서 CPU(Intel i7, 4.70 GHz) 기반 모형 대비 GPU 기반 모형의 계산속도 성능을 비교한 결과, 격자 간격이 증가할수록 CPU 기반 모형 대비 GPU 기반 모형의 연산 효율이 증가하였으며, 격자 간격이 $3200{\times}3200$일 때, CPU 기반 모형 대비 GPU 기반 모형의 연산 효율이 최대 약 150배 증가하였다.

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

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.69-78
• /
• 2008

The high-order numerical method based on the adaptive mesh refinement(AMR) on the quadrilateral unstructured grids has been developed in this paper. This adaptive-grid method, originally developed with MUSCL-TVD scheme, is now extended to the WENO (weighted essentially no-oscillatory) scheme with the Runge-Kutta time integration of fifth order in spatial and temporal accuracy. The multidimensional interpolation was studied in the preliminary research, which allows us to maintain the same order of accuracy for the computation of numerical flux between two adjacent cells of different levels. Some standard benchmark tests are done to validate this method for checking the overall capacity and efficiency of the present adaptive-grid technique.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.7-14
• /
• 2004

It is well known that preconditioning methods are efficient for convergence acceleration at compressible low Mach number flows. In this study, the original Euler equations and three preconditioners nondimensionalized differently are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as flux discretization and time integration respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning one produces Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.128-134
• /
• 2005

In the CFD area, the numerical analysis of high Mach number flow over a blunt-body poses many issues. Various numerical schemes have been developed to cover the issues, but the traditional schemes are still used widely due to the complexities of new schemes and intricacy of modifying the established codes. In the present study, the well-known Roe's FDS based on TVD-MUSCL scheme is used for the solution of very high Mach number three-dimensional flows posing carbuncle and non-physical phenomena in numerical analysis. A parametric study was carried out to account for the effects of the entropy fixing, grid configurations and initial condition. The carbuncle phenomena could be easily overcome by the entropy fixing, and the non-physical solution could be eliminated by the use of the modified initial condition regardless of entropy fixing and grid configurations.

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

지구 온난화 등에 의한 이상기후 현상으로 인해 야기되는 대규모 호우 등의 기상이변 현상은 댐 및 제방 붕괴와 같은 비상상황을 발생시키고 있고, 실제로 최근 10년간 태풍 등의 기상이변 현상으로 인해 낙동강의 여러 지류하천의 제방이 붕괴되는 피해가 발생했으며, 잇따른 피해들은 이 분야에 대한 종합적인 연구의 필요성을 증대시켰다. 본 연구의 목적은 이상홍수 및 국지성 호우에 의해서 하천 제방의 붕괴로 인한 제내지에서의 비상상황 발생에 대비하기 위해, 제내지에서의 범람홍수 해석결과를 통한 홍수위험지도 제작을 위한 tool로 활용함으로써 피해예상지역 내 주민 등의 신속한 대처를 통해 주민의 생명과 재산을 보호하기 위함에 있다. 이러한 하도 및 제내지에서의 범람홍수 양상을 효율적으로 계산하기 위해 실제 제방붕괴사례를 포함하는 사상에 대해 1차원 하천흐름 해석을 실시하였으며, 이를 통해 산정되는 제방붕괴 유출유량을 통해 제내지에 대한 2차원 범람홍수해석을 실시하였다. 1차원 제방붕괴 해석을 위해 FLDWAV 모형을 적용하였으며, 2차원 범람해석 모형으로 흐름의 전파양상을 정확하게 반영할 수 있는 상류이송기법인 Godunov 기법과 수치적인 계산 이전에 인접자료의 값을 이용하여 자료를 재구성하는 MUSCL(Monotone Upstream-centered Schemes for Conservation Laws) 기법을 사용하여 개발된 고정확도 유한체적모형을 적용하였다. 실제 제방붕괴 사상을 적용하기 위해 남강의 제방붕괴 사례를 고려하였으며, 2003년과 2006년에 각각 발생한 태풍 매미와 에위니아 사상에 대해 1차원 하천흐름해석 및 2차원 홍수범람해석을 실시하였다. 1차원 하천흐름해석에 대해서 하천 내에 위치한 수위표에서 관측된 실측수위를 통해 검증을 실시하였으며, 2차원 범람홍수모형에 의해 산정된 홍수범람범위는 침수흔적도를 통해 검증하였다. 본 연구에서 개발되어 적용된 2차원 범람해석 모형을 국가홍수위험지도 제작에 대해 활용할 수 있다면, 정확도 높은 통합홍수방재시스템 구축에 기여할 수 있을 것으로 기대된다.