• Journal of Korea Water Resources Association
• /
• v.44 no.10
• /
• pp.819-830
• /
• 2011

The multi-slope MUSCL, proposed by T. Buffard and S. Clain, determines slopes of conserved variables at each edge of a cell in the linear reconstructions of data. In this study, the second order accurate numerical model was developed according to the multi-slope MUSCL to solve the shallow water equations on the unstructured grids. The HLLL scheme of approximate Riemann solvers was used to calculate fluxes. For the review of the applicability of the developed model, the results of the model were compared to the 'isolated building test' and the 'model city flooding experiment' conducted as part of the IMPACT (Investigation of extreMe flood Processes And unCerTainty) project in Europe. There were limitations to predict abrupt rising of water depths by the resistance of model buildings and water depths at the specific locations among the buildings. But they were identified as the same problems also revealed in results of the other models to the same experiment. On the more refined meshes to the 'model city flooding experiment' simulated results showed good agreement with measurements. It was verified that the developed model simulated well the complex phenomena such as a dam-break problem and the urban inundation by flash floods.

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

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.6
• /
• pp.18-28
• /
• 2006

Numerical analysis of high Mach number flow over a blunt-body poses many difficulties and various numerical schemes have been suggested to overcome the problems. However, the new schemes were used in the limited fields of applications because of the lack of field experience compared to more than 20 years old numerical schemes and the intricacies of modifying the existing code for the special application. In this study, some tips to overcome the numerical difficulties in solving the 3D high-Mach number flows by using Roe's scheme, the most widely used for the past 25 years and adopted in many commercial codes, were examined without a correction of the algorithm or a modification of the CFD code. The well-known carbuncle phenomena of Riemann solvers could be remedied even for an extremely high Mach number by applying the entropy fixing function and a unphysical solution could be overcome by applying a simply modified initial condition regardless of the entropy fixing and grid configuration.

• Journal of Korea Water Resources Association
• /
• v.31 no.3
• /
• pp.279-290
• /
• 1998

The height and speed of the shock wave are critical data in flood-control operations or in the design of channel walls and bridges along rivers with high flow velocities. Therefore, a numerical model is needed for simulating flow discontinuity over a wide range of conditions. In this study, a governing equation. As a Riemann solver Roe(1981)'s one is used. The model employs the modified MUSCL for handling the unstructured grids in this research. this model that adopts the explicit tradditional twl dimmensional dam break problems, two hydraulic dam break model is simulations, and a steady state simulation in a curved channel. Conclusions of this research are as follows : 1) the finite volume method can be combined with the Godonov-type method that is useful for modeling shocks. Hence, the finite volume method is suitable for modeling shocks. 2) The finite volume model combined with the modified MUSCL is successful in modeling shock. Therefore, modified MUSCL is proved to be valid.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.35 no.6
• /
• pp.1237-1246
• /
• 2015

A novel technique was proposed to calculate fluxes accurately by separation of flow area into a part of step face which is dominated by flow resistance of it and an upper part which is relatively less affected by the step face in analyzing shallow-water flows over discontinuous topography. This technique gives fairly good agreement with exact solutions, 3D simulations, and experimental results. It has been possible to directly analyze shallow-water flows over discontinuous topography by the technique developed in this study. It is expected to apply the developed technique to accurate evaluation of overflows over weirs or retaining walls (riverside roads) and areas flooded by the inundation in the city covered in discontinuous topography.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.6
• /
• pp.429-445
• /
• 2008

Two-phase compressible flow fields of air-water are investigated numerically in the fixed Eulerian grid framework. The phase interface is captured via volume fractions of each phase. A way to model two phase compressible flows as a single phase one is found based on an equivalent equation of states of Tait's type for a multiphase cell. The equivalent single phase field is discretized using the Roe‘s approximate Riemann solver. Two approaches are tried to suppress the pressure oscillation phenomena at the phase interface, a passive advection of volume fraction and a direct pressure relaxation with the compressible form of volume fraction equation. The direct pressure equalizing method suppresses pressure oscillation successfully and generates sharp discontinuities, transmitting and reflecting acoustic waves naturally at the phase interface. In discretizing the compressible form of volume fraction equation, phase interfaces are geometrically reconstructed to minimize the numerical diffusion of volume fraction and relevant variables. The motion of a projectile in a water-filled tube which is fired by the release of highly pressurized air is simulated presuming the flow field as a two dimensional one, and several design factors affecting the projectile movement are investigated.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.145-150
• /
• 2005

Generally, Compressible Navier-Stokes codes are used to solve high mach number flows. But, Most of high mach number flows embrace low mach number flows. This phenomenon results in low convergence rate and non-physical solution in CFD analysis. So Many researchers developed preconditioning technique to solve these problems. This Study presents how to modify previous compressible N-S computer code with little changes of structure into preconditioned compressible N-S code applying Roe's Approximate Riemann Solver. And this study show developed preconditioning code is very well operated at all mach number flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.187-194
• /
• 2005

For the analysis of compressible multi-phase and real gas flows, characteristic form of Roe's Riemann solver was derived using real gas equation of state. It was extended to multi component reactive system considering variable specific heat. From this study, it is known that some correction should be made for the use of existing numerical algorithm. 1) Sonic speed and characteristic variable should be corrected with real gas effect. 2) Roe's average was applicable only with the assumption of constant properties. 3) Artificial damping term and characteristic variables should be corrected but their influences may not be significant.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.198-201
• /
• 2010

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