• 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 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 Korea Water Resources Association
• /
• v.44 no.2
• /
• pp.145-156
• /
• 2011

Two dimensional numerical model of high-order accuracy is developed to analyze complex flow including transition flow, discontinuous flow, and wave propagation to dry bed emerging at natural river flow. The bed slope term of two dimensional shallow water equation consisting of integral conservation law is treated efficiently by applying quasi-steady wave propagation scheme. In order to apply Finite Volume Method using Fractional Step Method, MUSCL scheme is applied based on HLL Riemann solver, which is second-order accurate in time and space. The TVD method is applied to prevent numerical oscillations in the second-order accurate scheme. The developed model is verified by comparing observed data of two dimenstional levee breach experiment and dam breach experiment containing structure at lower section of channel. Also effect of the source term is verified by applying to dam breach experiment considering the adverse slope channel.

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

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

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.10 no.2
• /
• pp.1-14
• /
• 2006

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma} $ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.6
• /
• pp.569-574
• /
• 2020

Though the discontinuous Galerkin (DG) method has been developed and applied to shallow water equations mainly in explicit schemes, they have been criticized for the limitation in treatment of bottom friction terms and severe CFL conditions. In this study, an implicit scheme is devised and applied to some representative benchmark problems. The linear triangular elements were employed and the Roe numerical fluxes were adopted for convective fluxes. To preserve TVD property, the slope limiter was employed. As the case studies, the model is applied to the flow around the cylinders and the dam-break flow. Then, the results are compared with the experimental and numerical data of previous studies and good agreements were observed.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.29-34
• /
• 1998

An efficient 3-dimensional compressible solver is developed using the second-order upwind TVD scheme and the multigrid diagonalized ADI method. The multigrid method is improved so that the present DADI algorithm obtains better convergence rates. Results are computed on Cray C90 computer for transonic unsaperated flows past ONERA-M6 wing to demonstrate the accuracy and efficiency. The results show good agreement with experimetal data. A reduction of four orders of residual for 3-dimensional transonic flow is obtained about 99 seconds.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.59-64
• /
• 2002

In this paper, the preconditioned multistage time stepping methods which are popular multigrid smoothers is implemented for the compressible Navier-Stokes calculation with full-coarsening multigrid method. The convergence characteristic of the point-Jacobi and Alternating direction line Jacobi(DDADI) preconditioners are studied. The performance of 2nd order upwind numerical fluxes such as 2nd order upwind TVD scheme and MUSCL-type linear reconstruction scheme are compared in the inviscid and viscous turbulent flow caculations.

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