• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.4 no.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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.10a
• /
• 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.11a
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.575-581
• /
• 1988

To integrate the two-dimensional Navier-Stokes equations numerically in multiply-connected flow regions, the vorticity-stream function formulation is used. The steady stream function value at the surface of the multibody, initially unknown, has been determined interactively by introducing a line integral which requires the single-valuedness of pressure at each interaction step. This procedure is relatively simpler and more efficient than the primitive variable formulation which requires much more computing time and shows poor convergence. Three doubly-connected flow problems are defined and numerically analyzed by the present method. The results have been compared either with earlier existing ones or with the experimental interferograms to demon-strate the validity of the presented method.

• Tribology and Lubricants
• /
• v.18 no.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
• /
• v.17 no.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.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2000.06a
• /
• 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.

• Journal of computational fluids engineering
• /
• v.7 no.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.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• 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
• /
• v.16 no.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.