• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 1999년도 춘계학술대회논문집
• /
• pp.253-256
• /
• 1999

The Streamline Upwind Petrov-Galerkin(SUPG) finite element method is used to solve the two-dimensional laminar and turbulent flow. The flow is simulated by averaged Navier-Stokes equations with a penalty function approach and the lograithmic(k-$\varepsilon$) turbulent model is employed to take into account its turbulent behavior. The near-wall viscous sub-layer model is employed to approach the dominant viscous effects in the near wall zones. To find a good-enough initial guess of the Newton-Raphson iteration solving Nonlinear Matrix the Incremental method is considered for momentum and the Incomplete logarithmic turbu-lent equations for Turbulence. The validation of our method is investigated in comparision with published experimental data.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2010년도 학술발표회
• /
• pp.749-753
• /
• 2010

본 논문에서는 하천에 유입되는 오염물질 중 대부분을 차지하는 비보존성 오염물질의 확산거동을 분석하기 위해 2차원 수심 평균된 이송분산방정식에 유한요소법을 적용하였다. 수치모형 구성을 위해 SUPG(Streamline-upwind Petrov-Galerkin)법을 이용한 가중잔차법을 사용하였다. 모의대상 수질인자는 BOD와 DO이며, BOD 농도 결과가 DO 농도 계산에서의 입력 자료로 이용되도록 상호 연계를 형성하였다. 모형의 검증을 위하여 직사각형 수로에 선원으로 연속주입하여 얻은 수치해와 해석해를 비교하였다. 비교결과 수치해와 해석해의 결과가 서로 일치하는 것을 볼 수 있었다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• pp.414-418
• /
• 2005

A finite element model is developed for the numerical simulation of bed elevation change in a curved channel. The SU/PG (Streamline-Upwind/Petrov-Galerkin) method is used to solve 2D shallow water equations and the BG (Bubnov-Galerkin) method is used for the Exner equation. For the time derivative terms, the Crank-Nicolson scheme is used. The developed model is a decoupled model in a sense that the bed elevation does not change simultaneously with the flow during the computational time step. The total load formula with is used for the sediment transport model. The slip conditions are described along the lateral boundaries. The effects of gravity force due to geometry change and the secondary flows in a curved channel are considered in the model. For the verification, the model is applied to two laboratory experiments. The first is $140^{\circ}$ bended channel data at Delft Hydraulics Laboratory and the second is $140^{\circ}$ bended channel data at Laboratory of Fluid Mechanics of the Delft University of Technology. The finite element grid is constructed with linear quadrilateral elements. It is found that the computed results are in good agreement with measured data, showing a point bar at the inner bank and a pool at the outer bank.

• 대한기계학회논문집B
• /
• 제32권10호
• /
• pp.754-760
• /
• 2008

Computation of moving interface by the level set method typically requires the reinitialization of level set function. An inaccurate estimation of level set function $\phi$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, an accurate and robust reinitialization process is essential to the simulation of free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates level set function directly using a normal vector on the interface without solving there-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1 splitting/SUPG (Streamline Upwind Petrov-Galerkin) FEM are adopted to discretize advection equation of the level set function and the incompressible Navier-Stokes equation, respectively. Advection equation and re-initialization process of free surface capturing are validated with benchmark problems, i.e., a broken dam flow and timereversed single vortex flow. The simulation results are in good agreement with the existing results.

• 한국수자원학회논문집
• /
• 제42권12호
• /
• pp.1053-1067
• /
• 2009

본 연구의 목적은 하도의 형상이 불규칙한 자연하천에서 2차원 흐름 특성을 해석하고 예측하기 위해 2차 요소를 이용한 정확하고 효과적인 상향가중 유한요소모형의 개발에 있다. 모형의 개발을 위해 선형 삼각형 요소, 선형 사각형 요소와 혼합요소를 적용하였고 2차 삼각형, 사각형 요소와 혼합요소를 적용하여 모형을 개발하였으며, 지배방정식의 수치적분식으로 Gauss Quadrature 방법을 사용하였다. 개발된 모형의 적용성 검증을 위해 하상융기가 있는 수로, U자형 수로 등에 모의를 실시하여 해석해 및 실측치와 비교 검토하였다. 모의 결과 2차 요소가 선형 요소에 비해 보다 정확한 해를 제공하는 것으로 판단되었으며 2차요소를 적용한 상용모형인 RMA-2 모형과 비교한 결과 본 연구 개발 모형이 보다 정확한 해를 나타내는 것을 확인할 수 있었다. 개발된 모형을 향후 자연하천에 적용할 경우 기존의 모형에 비해 향상된 결과를 얻을 수 있을 것으로 판단된다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2006년도 학술발표회 논문집
• /
• pp.425-429
• /
• 2006

The understanding and prediction of the behavior of flow in open channels are important to the solution of a wide variety of practical flow problems in water resources engineering. Recently, frequent drought has increased the necessity of an effective water resources control and management of river flows for reserving instream flow. The objective of this study is to develop an efficient and accurate finite element model based on Streamline Upwind/Petrov-Galerkin(SU/PG) scheme for analyzing and predicting two dimensional flow features in complex natural rivers. Several tests were performed in developed all elements(4-Node, 6-Node, 8-Node elements) for the purpose of validation and verification of the developed model. The U-shaped channel of flow and natural river of flow were performed for tests. The results were compared with these of laboratory experiments and RMA-2 model. Such results showed that solutions of high order elements were better accurate and improved than those of linear elements. Also, the suggested model displayed reasonable velocity distribution compare to RMA-2 model in meandering domain for application of natural river flow. Accordingly, the developed finite element model is feasible and produces reliable results for simulation of two dimensional natural river flow. Also, One contribution of this study is to present that results can lead to significant gain in analyzing the accurate flow behavior associated with hydraulic structure such as weir and water intake station and flow of chute and pool.

• Korea-Australia Rheology Journal
• /
• 제16권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.

• 대한기계학회논문집B
• /
• 제23권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.

• 대한기계학회논문집B
• /
• 제24권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.

• 대한토목학회논문집
• /
• 제31권6B호
• /
• pp.507-514
• /
• 2011

본 논문에서는 하천에 유입되는 열 오염물질의 혼합거동을 분석하기 위해 2차원 수심 평균된 이송-분산방정식에 유한요소법을 적용하여 수치 모형을 개발하였다. 유한요소법의 여러 수치기법 중 SUPG법을 적용하였으며, 복잡한 하천경계를 보다 정확히 재현할 수 있도록 삼각 및 사각 요소망의 혼용이 가능하도록 하였다. 열 오염물질의 거동을 표현하기 위해서 열 교환을 묘사하는 반응항을 평형온도와 수온과의 차이에 비례하는 식으로 나타내고, 열교환 계수 및 평형온도에 따라 수온의 변화가 적용되도록 방정식을 구성하였다. 모형의 검증을 위하여 직사각형 수로에 선원으로 연속주입하여 얻은 수치해와 1차원 정상상태의 해석해를 비교하였다. 비교결과 수치해와 해석해의 결과가 서로 일치하는 것으로 밝혀졌다. 모형의 현장적용을 위해 상수원 보호구역인 팔당댐 하류부터 잠실수중보까지 22.5 km 구간을 대상영역으로 하였다. 구리하수처리장 방류수에 의한 수온 변화를 모의한 결과 수질측정망 측정자료와 비교적 비슷한 경향이 나타났다. 본 연구에서 개발한 수치모형이 열 오염원 유입으로 인한 수온 변화를 잘 표현하는 것을 알 수 있었다.