• 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 Korea Water Resources Association
• /
• v.38 no.9 s.158
• /
• pp.727-736
• /
• 2005

In this study, the upwind McCormack scheme is introduced to combine the advantage of McCormack scheme, the second order accuracy and simplicity, and the advantage of the upwind scheme, to be applied to the discontinuous flows. This scheme also has another advantage of treating the source terms effectively. This model is approved through applying to the discontinuous flow case with the analytical solution, and the natural river with very strong source terms. Applications of the upwind McCormack scheme developed in this paper show good agreements with the analytical solution without numerical oscillation in existing McCormack scheme. Futhermore, applications to the natural river, the lower Han river with strong variation of bed and width, also show good results in case of both steady flow and unsteady flow. The upwind McCormack scheme in this study will be used for the analysis of flow in natural rivers effectively.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2002.05a
• /
• pp.633-639
• /
• 2002

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.129-137
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.378-386
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.572-583
• /
• 2019

This paper aims to assess the applicability of the Runge Kutta Discontinuous Galerkin-Direct Ghost Fluid Method to the internal explosion inside a water-filled tube, which previously was studied by many researchers in separate works. Once the explosive charge located at the inner center of the water-filled tube explodes, the tube wall is subjected to an extremely high intensity fluid loading and deformed. The deformation causes a modification of the field of fluid flow in the region near the water-structure interface so that has substantial influence on the response of the structure. To connect the structure and the fluid, valid data exchanges along the interface are essential. Classical fluid structure interaction simulations usually employ a matched meshing scheme which discretizes the fluid and structure domains using a single mesh density. The computational cost of fluid structure interaction simulations is usually governed by the structure because the size of time step may be determined by the density of structure mesh. The finer mesh density, the better solution, but more expensive computational cost. To reduce such computational cost, a non-matched meshing scheme which allows for different mesh densities is employed. The coupled numerical approach of this paper has fewer difficulties in the implementation and computation, compared to gas dynamics based approach which requires complicated analytical manipulations. It can also be applied to wider compressible, inviscid fluid flow analyses often found in underwater explosion events.

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

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.52-61
• /
• 2014

The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

• Journal of Korea Water Resources Association
• /
• v.36 no.4
• /
• pp.597-608
• /
• 2003

A transcritical flow occurs when the width and slope of a channel are varying abruptly. In this study, the transcritical flow in a two-dimensional open channel is analyzed by using the shallow-water equations. A weighted average flux scheme that has flux limiter with a total variation diminishing condition is introduced for a second-order accuracy in time and space, and non- spurious oscillations at discontinuous points. A HLLC method with three wane speeds is employed to calculate the Riemann problem. To overcome difficulties resulting from variation of channel sections in a two-dimensional analysis of transcritical flow, the numerical model is developed based on a generalized grid system.

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

Conventional high order interpolation schemes are limitative in several aspects mainly because they need data of neighboring cells at the reconstruction step. However, discontinuous Galerkin method and spectral volume method, two high order flux schemes which will be analyzed and compared in this paper, have an important benefit that they are not necessary to determine the flow gradients from data of neighboring cells or elements. These two schemes construct polynomial of variables within a cell so that even near wall or discontinuity, the high order does not deteriorate.