• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 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.

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

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

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

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

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

• Wind and Structures
• /
• 제14권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.

• Korea-Australia Rheology Journal
• /
• 제19권4호
• /
• pp.211-219
• /
• 2007

Here we demonstrate complex transient behavior of viscoelastic liquid described numerically with the Leonov model in straight and contraction channel flow domains. Finite element and implicit Euler time integration methods are employed for spatial discretization and time marching. In order to stabilize the computational procedure, the tensor-logarithmic formulation of the constitutive equation with SUPG and DEVSS algorithms is implemented. For completeness of numerical formulation, the so called traction boundaries are assigned for flow inlet and outlet boundaries. At the inlet, finite traction force in the flow direction with stress free condition is allocated whereas the traction free boundary is assigned at the outlet. The numerical result has illustrated severe forward-backward fluctuations of overall flow rate in inertial straight channel flow ultimately followed by steady state of forward flow. When the flow reversal occurs, the flow patterns exhibit quite complicated time variation of streamlines. In the inertialess flow, it takes much more time to reach the steady state in the contraction flow than in the straight pipe flow. Even in the inertialess case during startup contraction flow, quite distinctly altering flow patterns with the lapse of time have been observed such as appearing and vanishing of lip vortices, coexistence of multiple vortices at the contraction comer and their merging into one.

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

• Journal of Mechanical Science and Technology
• /
• 제20권10호
• /
• pp.1741-1752
• /
• 2006

A combined Streamline Upwind Petrov-Galerkin method (SUPG) and segregated finite element algorithm for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow is presented. The Streamline Upwind Petrov-Galerkin method is used for the analysis of viscous thermal flow in the fluid region, while the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the presented method is to consistently couple heat transfer along the fluid-solid interface. Four test cases, which are the conjugate Couette flow problem in parallel plate channel, the counter-flow in heat exchanger, the conjugate natural convection in a square cavity with a conducting wall, and the conjugate natural convection and conduction from heated cylinder in square cavity, are selected to evaluate efficiency of the presented method.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2005년도 춘계학술대회 논문집
• /
• pp.43-46
• /
• 2005

Micro molding technology is a promising mass production technology for polymer based microstructures. Mass production technologies such as the micro injection/compression molding, hot embossing, and micro reaction molding are already in use. In the present study, we have developed a numerical analysis system to simulate three-dimensional non-isothermal cavity filling for hot embossing, with a special emphasis on the free surface capturing. Precise free surface capturing has been successfully accomplished with the level set method, which is solved by means of the Runge-Kutta discontinuous Galerkin (RKDG) method. The RKDG method turns out to be excellent from the viewpoint of both numerical stability and accuracy of volume conservation. The Stokes equations are solved by the stabilized finite element method using the equal order tri-linear interpolation function. To prevent possible numerical oscillation in temperature Held we employ the streamline upwind Petrov-Galerkin (SUPG) method. With the developed code we investigated the detailed change of free surface shape in time during the mold filling. In the filling simulation of a simple rectangular cavity with repeating protruded parts, we find out that filling patterns are significantly influenced by the geometric characteristics such as the thickness of base plate and the aspect ratio and pitch of repeating microstructures. The numerical analysis system enables us to understand the basic flow and material deformation taking place during the cavity filling stage in microstructure fabrications.