• 대한토목학회논문집
• /
• 제28권5B호
• /
• pp.475-483
• /
• 2008

사고 등으로 인해 하천에 순간적으로 주입되는 오염물질의 2차원 이송-확산을 모의하기 위하여 SUPG법에 의한 2차원 유한요소 모형을 개발하고 주입된 질량과 모의 질량간의 오차를 산정하였다. 지배방정식으로부터 질량 오차 산정식을 유도하고 이를 간단한 예제에 적용하였다. 적용 예제로는 단순 직선수로에서 선원을 통한 질량주입과 일방향의 등질 유속장을 가지는 2차원 공간에서 점원을 통한 질량주입 문제를 가정하였으며, 기존의 Galerkin법과의 비교를 위하여 두가지 수치기법에 대한 모의 검증을 실시하였다. 검증결과 기존의 Galerkin법에 의한 모의 결과에 비하여 SUPG법을 통한 수치모형의 구성이 모의의 정확성뿐만 아니라 주입 질량의 보존에 있어서도 보다 정확한 방법임을 확인하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• 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.

• 대한기계학회논문집B
• /
• 제26권12호
• /
• pp.1663-1673
• /
• 2002

In the present paper, flow and heat transfer characteristics of confined impinging slot jets have been numerically investigated using a SIMPLE-based segregated SUPG finite element method. For laminar jets, it is shown that the skin friction coefficient obtained from the present SUPG formulation approaches the grid-independent Galerkin solution inducing negligible false diffusion in the flow field when a moderate number of grid points are used. For turbulent jets, the k-$\omega$turbulence model is adopted. The streamwise mean velocity and the heat transfer coefficient respectively agree very well with existing experimental data within limited ranges of parameters.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2004년도 학술대회지
• /
• pp.216-221
• /
• 2004

In two dimensional incompressible flaws, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint and hence achieve an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flaw analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flaw showed the present procedure to be stable, very efficient and useful in flaw analyses in conjunction with hierarchical elements.

• 한국해양공학회지
• /
• 제17권5호
• /
• pp.11-18
• /
• 2003

In two-dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure (HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint hence, achieving an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared (CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper and efficient iterative procedure for higher order finite element formulation has not been available, thus far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient, and useful in flow analyses, in conjunction with hierarchical elements.

• Water Engineering Research
• /
• 제4권1호
• /
• pp.45-58
• /
• 2003

Two-dimensional depth-averaged advection-dispersion equation was simulated using FEM. In the straight rectangular channel, the advection-dispersion processes are simulated so that these results can be compared with analyti-cal solutions for the transverse line injection and the point injection. In the straight domain the standard Galerkin method with the linear basis function is found to be inadequate to the advection-dispersion analysis compared to the upwind finite element scheme. The experimental data in the S-curved channel were compared with the result by the numerical model using SUPG(Streamline upwind Petrov-Galerkin) method.

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

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• 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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집B
• /
• pp.603-608
• /
• 2000

The characteristics of laminar and turbulent slot impinging jet flows are examined using segregated FEM with SUPG. Turbulent flows are modeled using $Wilcox^{(1)}$ $k-\;{\omega}$ turbulence model. The results are validated by comparing with velocity field of the existing experimental data. The distance of the target plate from the nozzle varies between 2, 4 and 5 times the slot jet width. Present study shows that the $k-\;{\omega}$ model gives results which agree well with the existing experimental data. In turbulence flows, the velocity profile of present calculation is more accurate than the existing numerical calculations. In laminar flows, We found tertiary vortex which was not found in the previous numerical study by M. $chen^{(6)}$ et al due to the numerical difference.

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