• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제4권2호
• /
• pp.111-133
• /
• 2000

A Legendre spectral tau approximation scheme for solving the two-dimensional stationary incompressible Stokes equations is considered. Based on the vorticity-stream function formulation and variational forms, boundary value and normal derivative of vorticity are computed. A factorization technique for matrix stems based on the Schur decomposition is derived. Several numerical experiments are performed.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2003년도 가을 학술발표회 논문집
• /
• pp.437-443
• /
• 2003

The object of this paper is the treatment of how to make the vorticity boundary condition instead of pressure in the primitive variable case. An improved algorithm for solving the vorticity-stream function equation is presented. The linear finite element approximation for the solution of Wavier-Stokes and Stokes flows is constructed. Not only regular domain but also complicate domain can be analyze d, using this formulation.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.49-55
• /
• 2009

This paper evaluates performances of a recently developed divergence-free finite element method based on Hermite interpolated stream functions. Velocity bases are derived from Hermite interpolated stream functions to form divergence-free basis functions. These velocity basis functions constitute a solenoidal function space, and the simple gradient of the Hermite functions constitute an irrotational function space. The incompressible Navier-Stokes equation is orthogonally decomposed into a solenoidal and an irrotational parts, and the decoupled Navier-Stokes equations are projected onto their corresponding spaces to form proper variational formulations. To access accuracy and convergence of the present algorithm, three test problems are selected. They are lid-driven cavity flow, flow over a backward-facing step and buoyancy-driven flow within a square enclosure. Hermite interpolation functions from cubic to quintic are chosen to run the test problems. Numerical results are shown. In all cases it has shown that the present method has performed well in accuracies and convergences. Moreover, the present method does not require an upwinding or a stabilized term.

• 대한기계학회논문집
• /
• 제12권3호
• /
• pp.575-581
• /
• 1988

본 연구에서는 앞서의 '압력의 단가성' 개념을 재도입하여 유동함수의 경계치 설정에 수치적으로 유용하게 쓸 수 있음을 보였고, 이를 다중연결된 물체에 관한 2차 원 응용문제들에서 수렴성이 나쁘고 많은 계산시간을 요구하는 종래의 원시변수들을 이용한 수식화 과정을 벗어나, 이제는 간편하고 명확한 유동함수-와도의 방법으로 비 교적 적은 컴퓨터 시간으로 유동장의 계산이 가능하게 되었음을 보였다.

• Tribology and Lubricants
• /
• 제18권6호
• /
• pp.411-417
• /
• 2002

Two dimensional slow viscous flow between a rotating cylinder and a translating plate is investigated using Stokes' approximation. An exact formal expression of the stream function is obtained by using the bipolar cylinder coordinates and Fourier series expansion. From the stream function obtained, the streamline patterns are shown and the pressure distribution in the flow field is determined. By integrating the stress distributions on the cylinder, the farce and the moment exerted on the cylinder are calculated. The flow rate through the gap between the cylinder and the plate is also determined as a function of the distance between the cylinder and the plate. Special attention is directed to the case of very small distance between the cylinder and the plate concerned with the lubrication theory and the minimum pressure is calculated to explain a possible cavitation.

• Tribology and Lubricants
• /
• 제17권5호
• /
• pp.391-398
• /
• 2001

Two dimensional slow viscous flow around two counter-rotating equal cylinders is investigated based on Stokes'approximation. An exact formal expression of the stream function is obtained by using the bipolar cylinder coordinates and Fourier series expansion. From the stream function obtained, the streamline patterns around the cylinders are shown and the pressure distribution in the flow field is determined. By integrating the stress distributions on the cylinder, the force and the moment exerted on the cylinder are calculated. The flow rate through the gap between the two cylinders is also determined as the distance between two cylinders varies. Special attention is directed to the case of very small distance between two cylinders concerned with the lubrication theory and the minimum pressure is calculated to explain a possible cavitation.

• 한국윤활학회:학술대회논문집
• /
• 한국윤활학회 2000년도 제31회 춘계학술대회
• /
• pp.142-149
• /
• 2000

Two dimensional slow viscous flow around two counter-rotating equal cylinders is Investigated based on Stokes' approximation. An exact formal expression of the stream function is obtained by using the bipolar cylinder coordinates and Fourier series expansion. From the stream function obtained, the streamline patterns around the cylinders are shown and the pressure distribution In the flow field is determined. By Integrating the stress distribution on the cylinder, the force and the moment exerted on the cylinder are calculated. The flow rate through the gap between the two cylinders is determined as the distance between two cylinders vary. It Is also revealed that the velocity at the far field has finite non-zero value. Special attention is directed to the case of very small distances between two cylinders by way of the lubrication theory.

• 한국전산유체공학회지
• /
• 제7권3호
• /
• pp.53-59
• /
• 2002

A two-dimensional laminar flow through a channel with a pair of symmetric vertical fins is investigated. At far up- and down-stream from the fins, the plane Poiseuille flow exists in the channel. The Stokes flow for this channel is first investigated analytically and then the other laminar flows by numerical method. For analytic method, the method of eigen function expansion and collocation method are employed. In numerical solution for laminar flows, finite difference method(FDM) is used to obtain vorticity and stream function. From the results, the streamline patterns are shown and the additional pressure drop due to the attached fins and the force exerted on the fin are calculated. It is clear that the force depends on the length of fins and Reynolds number. When the Reynolds number exceeds a critical value, the flow becomes asymmetric. This critical Reynolds number Re/sub c/ depends on the length of the fins.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2009년도 학술발표회 초록집
• /
• pp.516-523
• /
• 2009

The artificial compressibility (AC) method for the incompressible Navier-Stokes equations in the generalized curvilinear coordinates using the primitive form is implemented. The main advantage of the AC approach is that the resulting system of equations resembles the system of compressible N-S equations and can thus be integrated in time using standard, well-established time-marching methods. The errors, which are the odd-even oscillation, for pressure field in using the artificial compressibility can be eliminated by using the $4^{th}$ order artificial dissipation term which is explicitly included. Even though this paper focuses exclusively on 2D laminar flows to validate and assess the performance of this solver, this numerical method is general enough so that it can be readily extended to carry out 3D URANS simulation of engineering flows. This algorithm yields practically identical velocity profiles and primary vortex and secondary vortices that are in excellent overall agreement with the results of the vorticity-stream function formulation (Ghia et al., 1982). However, the grid resolution have to be required to be large enough to express the various vortices.

• Journal of applied mathematics & informatics
• /
• 제16권1_2호
• /
• pp.383-405
• /
• 2004

We study the incomprssible Navier Stokes equations for the flow inside contraction geometry. The governing equations are expressed in the vorticity-stream function formulations. A rectangular computational domain is arised by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of acurvilinear coordinate system by transforming the governing equation into computational plane. The transformed equations are approximated using central differences and solved simultaneously by successive over relaxation iteration. The time dependent of the vorticity equation solved by using explicit marching procedure. We will apply the technique on several irregular-shapes.