• Proceedings of the Korea Water Resources Association Conference
• /
• 1998.05a
• /
• pp.259-264
• /
• 1998

• Proceedings of the Korea Water Resources Association Conference
• /
• 1999.05a
• /
• pp.547-552
• /
• 1999

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.289-294
• /
• 2010

Numerical stud of an oscillating body in incompressible fluid is performed. Stabilized finite element method comprising of Streamline-Upwind/Petrov-Galerkin (SUPG) and Pressure-Stabilizing/Petrov-Galerkin (PSPG) formulations of linear triangular elements was employed to solve 2D incompressible Navier-Stokes equations whereas the motion of the body was considered by incorporating the arbitrary Langrangian-Eulerian(ALE) formulation. An algebraic moving mesh strategy is utilized for obtaining body conforming mesh deformation at each time step. Two tests cases, namely motion of a circular cylinder and of an airfoil in incompressible flow were analyzed. The model is first validated against the stationary cases and then the capability to handle moving boundaries is demonstrated.

• Journal of Mechanical Science and Technology
• /
• v.17 no.7
• /
• pp.986-998
• /
• 2003

In this paper, an adaptive numerical integration scheme, which does not need non-overlapping and contiguous integration meshes, is proposed for the MLPG (Meshless Local Petrov-Galerkin) method. In the proposed algorithm, the integration points are located between the neighboring nodes to properly consider the irregular nodal distribution, and the nodal points are also included as integration points. For numerical integration without well-defined meshes, the Shepard shape function is adopted to approximate the integrand in the local symmetric weak form, by the values of the integrand at the integration points. This procedure makes it possible to integrate the local symmetric weak form without any integration meshes (non-overlapping and contiguous integration domains). The convergence tests are performed, to investigate the present scheme and several numerical examples are analyzed by using the proposed scheme.

• Journal of Korea Water Resources Association
• /
• v.34 no.2
• /
• pp.107-118
• /
• 2001

This study was to develop a tow-dimensional model system for the hydrodynamic analysis and to apply the system on the downstream of Han River. it is performed to design a GIS-based hydrodynamic system for the scientific shallow water profile analysis, and to compare hydrodraulic modeling is the Petrov-Galerkin's finite element method for flow prediction model. This study was to construct a GIS-based river flow system, and it is useful for supporting user's decision making for the on-line status through various analysis. We expect that the results from this study can be used as one of the guidelines for river analysis and management system in order rivers, reservoirs, and estuaries.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.10
• /
• pp.55-65
• /
• 2014

Recent researches on debris flow is focused on understanding its movement mechanism and building a numerical simulator to predict its behavior. However, previous simulators emulating fluid-like debris flow have limitations in numerical stability, geometric modeling and application of various boundary conditions. In this study, depth integration is applied to continuity equation and force equilibrium for debris flow. Thickness of sediment, and average velocities in x and y flow direction are chosen for main variables in the analysis, which improve numerical stability in the area with zero thickness. Petrov-Galerkin formulation uses a discontinuous test function of the weighted matrix from DG scheme. Presented mechanical constitutive model combines fluid and granular behaviors for debris flow. Effects on slope angle, inducing debris height, and bottom friction resistance are investigated for a simple slope. Numerical results also show the effect of embankment at the bottom of the slope. Developed numerical simulator can assess various risk factors for the expected area of debris flow, and facilitate embankment design in order to minimize damage.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.3B
• /
• pp.175-183
• /
• 2012

In this study, Galerkin scheme and SU/PG scheme of Petrov-Galerkin family were applied to the shallow water equations and a finite element model for shallow water flow was developed. Numerical simulations were conducted in several flumes with convection-dominated flow condition. Flow simulation of channel with slender structure in the water course revealed that Galerkin and SU/PG schemes showed similar results under very low Fr number and Re number condition. However, when the Fr number increased up to 1.58, Galerkin scheme did not converge while SU/PG scheme produced stable solutions after 5 iterations by Newton-Raphson method. For the transcritical flow simulation in diverging channel, the present model predicted the hydraulic jump accurately in terms of the jump location, the depth slope, and the flow depth after jump, and the numerical results showed good agreements with the hydraulic experiments carried out by Khalifa(1980). In the oblique hydraulic jump simulation, in which convection-dominated supercritical flow (Fr=2.74) evolves, Galerkin scheme blew up just after the first iteration of the initial time step. However, SU/PG scheme captured the boundary of oblique hydraulic jump accurately without numerical oscillation. The maximum errors quantified with exact solutions were less than 0.2% in water depth and velocity calculations, and thereby SU/PG scheme predicted the oblique hydraulic jump phenomena more accurately compared with the previous studies (Levin et al., 2006; Ricchiuto et al., 2007).

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.237-246
• /
• 1999

A finite element model is studied to simulate unsteady free surface flow based on dynamic wave equation and collocation method. The collocation method is used in conjunction with Hermite polynomials, and resulting matrix equations are solved by skyline method. The model is verified by applying to hydraulic jump, nonlinear disturbance propagation and dam-break flow in a horizontal frictionless channel. The computed results are compared with those by Bubnov-Galerkin and Petrov-Galerkin methods. It is also applied to the North Han River to simulate the floodwave propagation. The computed results have good agreements with those of DWOPER model in terms of discharge hydrographs. The suggested model has proven to be one of the promising scheme for simulating the gradually and rapidly varied unsteady flow in open channels.

• Wind and Structures
• /
• v.14 no.5
• /
• pp.465-480
• /
• 2011

In this paper, a stabilized large eddy simulation technique is developed to predict turbulent flow with high Reynolds number. Streamline Upwind Petrov-Galerkin (SUPG) stabilized method and three-step technique are both implemented for the finite element formulation of Smagorinsky sub-grid scale (SGS) model. Temporal discretization is performed using three-step technique with viscous term treated implicitly. And the pressure is computed from Poisson equation derived from the incompressible condition. Then two numerical examples of turbulent flow with high Reynolds number are discussed. One is lid driven flow at Re = $10^5$ in a triangular cavity, the other is turbulent flow past a square cylinder at Re = 22000. Results show that the present technique can effectively suppress the instabilities of turbulent flow caused by traditional FEM and well predict the unsteady flow even with coarse mesh.