• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.2
• /
• pp.479-492
• /
• 2014

In this study, was proposed a numerical scheme imposing exact solutions as the internal boundary conditions for the shallow-water flows over a discontinuous transverse structure such as a step. The HLLL approximate Riemann solver with the MUSCL was used for the test of the proposed scheme. Very good agreement was obtained between simulations and exact solutions for various problems of the shallow-water flows over a step. In addition, results by the numerical model showed good agreement with those of dam-break experiments over a step and stepped chute one. Developed model can simulate the shallow-water flows over discontinuous bottom such as a drop structure without additional rating curve or topography smoothing. Given the proper evaluations for the flow resistance by a step and the energy loss by the nappe flow in the future, could be simulated flooding and drying of the shallow-water flows over discontinuous topography such as a weir or the river road with retaining wall.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.1
• /
• pp.27-33
• /
• 2009

Shock tube flow measurement has been often hampered a finite opening time of diaphragm, but there is no systematic work to investigate its effect on the shock tube flow. In the present study, both the experimental and computational works have been performed on the shock tube flows at low pressure ratios. The computational analysis has been performed using the two-dimensional, unsteady, compressible Navier-Stokes equations, based upon a TVD MUSCL finite difference scheme. It is known that the present computational results reproduce the experimental data with good accuracy and simulate successfully the process of diaphragm opening as a function of time. The concept of an imaginary center is introduced to quantify the non-centered expansion wave due to a finite opening time of diaphragm. The results obtained show that the diaphragm opening time is reduced as the initial pressure ratio of shock tube increases, leading to the effect of a finite opening time of diaphragm to be more remarkable at low pressure ratios.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.2
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• 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 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 Korea Water Resources Association
• /
• v.46 no.2
• /
• pp.139-153
• /
• 2013

For analyzing shallow-water flows over the uneven bottom, the HLLL scheme and the divergence form for bed slope source term (DFB) technique, respectively were applied to the flux gradient and the bottom gradient source terms in a finite-volume model for the shallow water equations. And also the model incorporated the volume/free-surface relationship (VFR) to consider the partially submerged cells (PSC). It was identified that a simpler version of the weighted surface-depth gradient method in the MUSCL was equivalent to the original one in the accuracy for 1D steady flows. It was verified that the flux gradient term and the bottom gradient source term were well-balanced exactly by the VFR for the 1D PSC. The VFR for the triangular PSC settled the problem which the governing equations were not well-balanced by the DFB technique for the 2D PSC. There were good agreements in simulations and experiments for 2D dam-break flows over a triangular sill and a round bump. In addition, the partial dam-break flow was successfully simulated for flooding of roughnesses in an irregular bottom as well as a sloping one. Therefore, this model is expected to be applied to the real river with uneven topography.

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

• The Bulletin of The Korean Astronomical Society
• /
• v.38 no.2
• /
• pp.101.2-101.2
• /
• 2013

We present a new computational fluid dynamic (CFD) simulation code. The code employs the moving and polyhedral unstructured mesh scheme, which is known as a superior approach to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for surfaces of every Voronoi cell to update the hydrodynamic states as well as to move former generated meshes. For the second-order accuracy, the MUSCL-Hancock scheme is implemented. To increase efficiency for generating Voronoi tessellation we also develop the incremental expanding method, by which the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We will discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

• Journal of computational fluids engineering
• /
• v.13 no.1
• /
• pp.7-13
• /
• 2008

Several high resolution schemes such as OSHER, MUSCL, SMART, GAMMA, WACEB and CUBISTA are comparatively studied with respect to the accurate capturing of fluid interfaces throughout the application to two typical test cases of a translation test and a collapsing water column problem with a return wave. It is accomplished by implementing the high resolution schemes in the in-house CFD code(PowerCFD) for computing 3-D flow with an unstructured cell-centered method and an interface capturing method, which is based on the finite-volume technique and fully conservative. The calculated results show that SMART scheme gives the best performance with respect to accuracy and robustness.

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