• Korea-Australia Rheology Journal
• /
• v.16 no.3
• /
• pp.117-128
• /
• 2004

The effect of molecular parameters on the steady vortex behaviors in the inertial viscoelastic flow past a cylinder has been investigated. FENE-CR model was considered as a constitutive equation. A recently developed iterative solution method (Kim et al., (in press)) was found to be successfully applicable to the computation of inertial viscoelastic flows. The high-resolution computations were carried out to understand the detailed flow behaviors based on the efficient iterative solution method armed with ILU(0) type pre-conditioner and BiCGSTAB method. The discrete elastic viscous split stress-G/streamline upwind Petrov Galerkin (DEVSS-G/SUPG) formulation was adopted as a stabilization method. The vortex size decreased as elasticity increases. However, the vortex enhancement was also observed in the case of large extensibility, which means that the vortex behavior is strongly dependent upon the material parameters. The longitudinal gradient of normal stress was found to retard the formation of vortex, whereas the extensional viscosity played a role in the vortex enhancement. The present results are expected to be helpful for understanding the inertial vortex dynamics of viscoelastic fluids in the flow past a confined cylinder.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.43-57
• /
• 1999

The fluid flow and heat transfer in a thin liquid film are investigated numerically. The flow Is assumed to be two-dimensional laminar and surface tension is considered. The most important characteristics of this flow is the existence of a hydraulic jump through which the flow undergoes very sharp and discontinuous change. Arbitrary Lagrangian-Eulerian(ALE) method is used to describe moving free boundary and a modified SIMPLE algorithm based on streamline upwind Petrov-Galerkin(SUPG) finite element method is used for time marching iterative solution. The numerical results obtained by solving unsteady full Navier-Stokes equations are presented for planar and radial flows subject to constant wall temperature or constant wall heat flux, and compared with available experimental data. It Is discussed systematically how the inlet Reynolds and Froude numbers and surface tension affect the formation of a hydraulic jump. In particular, the effect of temperature dependent fluid properties is also discussed.

• Structural Engineering and Mechanics
• /
• v.76 no.4
• /
• pp.505-512
• /
• 2020

This paper presents an error estimation technique for 2-D crack analysis by an enriched natural element (more exactly, enriched Petrov-Galerkin NEM). A bare solution was approximated by PG-NEM using Laplace interpolation functions. Meanwhile, an accurate quasi-exact solution was obtained by a combined use of enriched PG-NEM and the global patch recovery. The Laplace interpolation functions are enriched with the near-tip singular fields, and the approximate solution obtained by enriched PG-NEM was enhanced by the global patch recovery. The quantitative error amount is measured in terms of the energy norm, and the accuracy (i.e., the effective index) of the proposed method was evaluated using the errors which obtained by FEM using a very fine mesh. The error distribution was investigated by calculating the local element-wise errors, from which it has been found that the relative high errors occurs in the vicinity of crack tip. The differences between the enriched and non-enriched PG-NEMs have been investigated from the effective index, the error distribution, and the convergence rate. From the comparison, it has been justified that the enriched PG-NEM provides much more accurate error information than the non-enriched PG-NEM.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.244-244
• /
• 2022

4차 산업혁명에 따라 유역 및 하천관리 사업, 각종 풍수해 예방사업 분야에 다양한 스마트기술이 도입되고 있으나 건설현장 침수 피해 사고는 지속적으로 발생하고 있다. 굴착공사 현장에서 발생할 수 있는 침수피해를 사전에 예측하기 위해서는 공정별로 변화하는 현장상황을 반영하여 다양한 강우 시나리오를 기반으로 침수 예측 모델링이 선행되어야한다. 따라서 본 연구에서는 2차원 동수역학 모형인 HDM-2D 모형을 기반으로 굴착공사 현장 침수피해 예측 모델을 개발하여 굴착공사현장 침수 예·경보 시스템에 탑재하고자 한다. 침수피해 예측 모델은 천수방정식을 Petrov-Galerkin stabilized scheme 으로 이산화하여 해석하는 수평 2차원 동수역학 흐름해석모델로서 수로 및 지표면 등 다양한 지형 상황에서의 물 흐름을 상세하게 해석할 수 있다. 지형자료 생성 이후 경계조건 부여를 통해 수행되며 침수발생지역의 유속, 수심, 수위를 취득할 수 있다. 배수지나 굴착공사 현장에 2차원 흐름해석을 적용하는 경우 지형의 경사나 배치가 공간에 따라 변화하므로 불연속적인 흐름을 유발하여 모의결과의 계산 오차를 검토해야 한다. 2차원 침수피해 예측 모형의 정확성을 확인하기 위해 지면 돌출부가 있는 흐름 문제와 테스트베드 대상지에 침수해석 모형을 적용하였다. 돌출부 흐름 문제의 경우 돌출부를 지나며 발생하는 유속과 수심 모의 결과를 상용모형과 비교검증 하였으며 테스트베드 대상지역에 침수피해예측모형을 적용했을때 지형 경사에 따른 흐름의 변화와 침수양상을 확인할 수 있었다.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1053-1067
• /
• 2009

The objective of this study is to develop an efficient and accurate quadratic finite element model based on Streamline Upwind/Petrov Galerkin (SU/PG) scheme for analyzing and predicting two dimensional flow features in complex natural rivers. For a development of model, quadratic tin, quadrilateral and mixed elements as well as linear tin, quadrilateral and mixed elements were used in the model. Also, this model was developed through reinforcement of Gauss Quadrature which was necessary to integral of governing equation. Several tests for bottom-rising channel and U-type channel were performed for the purpose of validation and verification of the developed model. Such results showed that solutions of second order elements are better accurate and improved than those of linear elements. Results obtained by the developed model and RMA-2 model are compared, and the results for the developed model were better accurate than those of RMA-2 model. In the future if the developed model is applied in natural rivers, it can provide better accurate results than those of existing model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.6B
• /
• pp.507-514
• /
• 2011

A numerical model has been developed by employing a finite element method to simulate the depth-averaged 2-D dispersion of the heat pollutant, which is an important pollutant material in natural streams. Among the finite element methods, the Streamline Upwind/Petrov Galerkin (SUPG) method was applied. Also both linear and quadratic elements can be applied so that irregular river boundaries can be easily represented. To show the movement of heat pollutants, the reaction term describing heat transfer was represented as an equation in which sink/source term is proportional to the difference between the equilibrium temperature and water surface temperature. The equation was expressed so that the water surface temperature changes according to the temperature transfer coefficient and the equilibrium temperature. For the calibration of the model developed, analytic and numerical results from a case of rectangular channel with full width continuous injection have been compared in a steady state. The comparisons showed that the numerical results were in good agreement with analytical solutions. The application site was selected from the downstream of Paldang dam to Jamsil submerged weir, and overall length of this site is about 22.5 km. The change of water temperature caused by the discharge from the Guri sewage treatment plant has been simulated, and results were similar to the observed data. Overall it is concluded that the developed model can represent the water temperature changes due to heat transport accurately. But the verification using observed data will further enhance the validity of the model.