• 한국전산유체공학회:학술대회논문집
• /
• 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.

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

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

The unsteady supersonic flow over two- and three-Dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k - w turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. Numerical method is upwind TVD scheme based on the flux vector split with the Van Leer limiters, and time accuracy is used explicit 4th stage Runge-Kutta scheme. Cavity flows are Comparison of two- and three-dimensional. The cavity has a L/D ratio of 3 for two-dimensional case. and same L/D and W/D ratio is 1 for three-dimensional case. The Mach and Reynolds numbers are held constant at 1.5 and 450000 respectively. For the three-dimensional case, the flow field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follow Rossiter's formula. On the other hand, the self-sustained oscillating flow transitions to a 'wake mode' for the two-dimensional simulation, with more violent fluctuations inside the cavity.

• Journal of Mechanical Science and Technology
• /
• v.15 no.3
• /
• pp.366-373
• /
• 2001

A numerical analysis of shock wave/boundary layer interaction in transonic/supersonic axial flow compressor cascade has been performed by using a characteristics upwind Navier-Stokes method with various turbulence models. Two equation turbulence models were applied to transonic/supersonic flows over a NACA 0012 airfoil. The results are superion to those from an algebraic turbulence model. High order TVD schemes predicted shock wave/boundary layer interactions reasonably well. However, the prediction of SWBLI depends more on turbulence models than high order schemes. In a supersonic axial flow cascade at M=1.59 and exit/inlet static pressure ratio of 2.21, k-$\omega$ and Shear Stress Transport (SST) models were numerically stables. However, the k-$\omega$ model predicted thicker shock waves in the flow passage. Losses due to shock/shock and shock/boundary layer interactions in transonic/supersonic compressor flowfields can be higher losses than viscous losses due to flow separation and viscous dissipation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1383-1393
• /
• 2014

Recently, with rapid improvement in computer hardware and theoretical development in the field of computational fluid dynamics, high-order accurate schemes also have been applied in the realm of computational hydraulics. In this study, numerical solutions of 1D shallow water equations are presented with TVD Runge-Kutta discontinuous Galerkin (RKDG) finite element method. The transcritical flows such as dam-break flows due to instant dam failure and transcritical flow with bottom elevation change were studied. As a formulation of approximate Riemann solver, the local Lax-Friedrichs (LLF), Roe, HLL flux schemes were employed and MUSCL slope limiter was used to eliminate unnecessary numerical oscillations. The developed model was applied to 1D dam break and transcritical flow. The results were compared to the exact solutions and experimental data.

• 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 for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.21-29
• /
• 2002

In the present study, a passive control method, using a porous wall and cavity system, is applied to the shock wave/boundary layer interactions in transonic moist air flow. The two-dimensional, unsteady, compressible, Navier-Stokes equations, which are fully coupled with a droplet growth equation, are solved by the third-order MUSCL type TVD finite difference scheme. Baldwin-Lomax model is employed to close the governing equations. In order to investigate the effectiveness of the present control method, the total pressure loss of the flow and the time-dependent behaviour of shock motions are analyzed in detail. The computed results show that the present passive control method considerably reduces the total pressure losses due to the shock wave/boundary layer interaction in transonic moist air flow and suppresses the unsteady shock wave motions over the airfoil as well. It is also found that the location of the porous ventilation significantly affects the control effectiveness.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.45-50
• /
• 2000

We study a conceptual numerical model on shock-vortex interaction setting an impulsive shock in a compressible vertex. Navier-Stokes equations are solved for the investigation of interactive structure and acoustic wave propagation. The rotationally symmetric vortex enforces two compression-expansion pairs resultantly forming a quadrupolar shape. These compressive and expansive waves cylindrically propagate to the far field and turn to acoustic waves. Using a fine uniform Cartesian grid system and a TVD-high resolution method, the flow data irl: precisely obtained to extend our interest to the sound source.