• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.475-483
• /
• 2008

2D finite element model for analysis of transport of accidentally released pollutants in the flow was developed by SUPG method, and the mass balance of this model was checked though two example problems: line source and point source problem in the straight channel and unidirectional 2D flow field, respectively. All the test cases were simulated with both SUPG and conventional Galerkin method to compare the accuraccy of the numerical mass balance. Test results show that the model with SUPG can adequately conserve the released mass though simulation than the model using Galerkin method, so the developed model verified to be appropriate to solve this accidental mass release problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.6 s.237
• /
• pp.703-710
• /
• 2005

A finite element discretization of the advection and redistancing equations of level set method has been studied. It has been shown that Galerkin spatial discretization combined with Crank-Nicolson temporal discretization of the advection equation of level set yields a good result and that consistent streamline upwind Petrov-Galerkin(CSUPG) discretization of the redistancing equation gives satisfactory solutions for two test problems while the solutions of streamline upwind Petrov-Galerkin(SUPG) discretization are dissipated by the numerical diffusion added for the stability of a hyperbolic system. Furthermore, it has been found that the solutions obtained by CSUPG method are comparable to those by second order ENO method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2667-2677
• /
• 1995

A new numerical algorithm using equal order linear finite element and fractional step method has been developed which is capable of analyzing unsteady fluid flow and heat transfer problems. Streamline Upwind Petrov-Galerkin (SUPG) method is used for the weighted residual formulation of the Navier-Stokes equations. It is shown that fractional step method, in which pressure term is splitted from the momentum equation, reduces computer memory and computing time. In addition, since pressure equation is derived without any approximation procedure unlike in the previously developed SIMPLE algorithm based FEM codes, the present numerical algorithm gives more accurate results than them. The present algorithm has been applied preferentially to the well known bench mark problems associated with steady flow and heat transfer, and proves to be more efficient and accurate.

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

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

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.341-341
• /
• 2016

2차류는 주흐름 방향 유속에 비하여 작지만 유속을 재분포 시켜 흐름의 교란을 발생시키는 매우 중요한 인자이다(Woo, 2001). 2차류는 만곡부, 합류부, 분류부에서 흐름분리와 함께 나타나는 흐름특성 중 하나이지만 각 지점별로 2차류의 분포는 서로 상이하다. 하지만 국내에서는 만곡부에서의 2차류에 대한 영향에 대해서만 주로 연구가 이루어졌으며, 합류부와 분류부에 대해서는 연구가 아직 체계적으로 이루어지지 못하고 있는 실정이다. 따라서 본 연구에서는 합류부와 분류부에서 2차류 영향을 고려한 2차원 수치모형을 적용하여 흐름특성 변화를 수치모의 하였다. 2차류 흐름을 고려한 수치모의를 위해 TELEMAC-2D 수치모형을 이용하였으며, 2차류의 영향을 고려하지 않고 계산한 SUPG기법을 적용한 모의결과에 대하여 Bernard and Schneider(1992)가 제안한 경험계수를 적용 후 그 결과를 실험결과와 비교분석하였다. 합류부(confluence)에서는 본류와 지류의 유량비가 증가할수록 흐름분리구역 길이와 최대폭이 증가하지만, 분류부(Bifurcation)에서는 분류수로 유량비가 증가할수록 흐름분리구역 길이와 최대폭이 감소하는 것을 확인할 수 있었다. 특히, 2차류를 제외한 수치모의에서 실내시험에서 제시한 분류유량비와 수심 및 유속분포를 동시에 충족시키기 어렵다. 반면 2차류 영향을 고려시 분기수로내 통수능을 감소시키는 흐름분리구역과 2차류의 상호작용에 의한 흐름정체효과로 분류유량비가 감소하였으며, 2차류 영향을 고려하여 모의한 결과 분류부 수치모의의 안정성과 정확성을 향상되었다. 하지만 본 모형에서 적용한 2차류 흐름관련 매개변수는 흐름특성에 따라 시행착오법에 따라 산출해야 하는 특성 때문에 향후에 추가적인 연구가 필요하다.