• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.5
• /
• pp.415-423
• /
• 2015

In this study, we present a grid refinement model in the lattice Boltzmann method (LBM) for two-dimensional incompressible fluid flow. That is, the model combines the desirable features of the lattice Boltzmann method and stream function-vorticity formulations. In order to obtain an accurate result, very fine grid (or lattice) is required near the solid boundary. Therefore, the grid refinement model is used in the lattice Boltzmann method for stream function-vorticity formulation. This approach is more efficient in that it can obtain the same accurate solution as that in single-block approach even if few lattices are used for computation. In order to validate the grid refinement approach for the stream function-vorticity formulation, the numerical simulations of lid-driven cavity flows were performed and good results were obtained.

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

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

• Journal of Electrical Engineering and Technology
• /
• v.4 no.2
• /
• pp.249-254
• /
• 2009

This article deals with the analysis of a coupling between stationary Maxwell's equations, the transient state Navier-Stokes and thermal equations. The resolution of these equations is obtained by introducing the magnetic vector potential A, the vorticity ${\xi}$, the stream function ${\psi}$ and the temperature T. The flux density, the electromagnetic thrust, the electric power density, the velocity, the pressure and the temperature are graphically visualized. Also, the influence of the frequency is presented.

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.6
• /
• pp.1163-1170
• /
• 1992

A computer program was developed to analyze the two-dimensional unsteady incompressible viscous flow behind a rectangular bluff body between two parallel plates. The Peaceman-Rachford alternating direction implicit numerical method and Wachspress parameter were adopted to solve the governing equations in vorticity-transport and stream function formulation. The steady state flow and the vortex flow behind a rectangular bluff body in a chemical were investigated for Reynolds numbers of 200 and 500. The vortex shedding was generated by a physical pertubation numerically imposed at the center of the flow field for a short time. It was observed that the perturbed flow became periodic after a transient period.

• The Korean Journal of Rheology
• /
• v.2 no.1
• /
• pp.33-45
• /
• 1990

This study analyzes the planar 4:1 contraction flow of viscoelastic fluids with retardation time using finite volume method. To consider separately the elasticity effect of the viscoelastic fluid without shear thinn-ing effect, Oldroyd B liquid model is adopted for the numerical simulation. Instead of the stream function-vorticity formulation, SIMPLER algorithm with staggered grid system which incorporates primitive variable has been introduced in discretizing the momentum equations. An upwind corrected scheme has been used in discetizing the constitutive equations for the non-Newtonian part of the stress. The size of the corner vortex is shown to be slightly influenced by the Weissenberg number. However as the Weissenberg number is increased the chang-ing of the vortex shape agrees qualitatively well with some experimental studies.

• Journal of the Korean Society of Safety
• /
• v.10 no.2
• /
• pp.56-63
• /
• 1995

A simulation of contaminant dispersion in a water reservoir has been done using 2-D finite difference method(FDM). The steady state velocity field of the reservoir was computed using stream function-vorticity formulation of Wavier-Stokes equation and continuity equation. Based on the computed steady state velocity field, the transient convective diffusion equation of the contaminant dispersion was computed. For the 1m$\times$1m reservoir model with inlet and outlet attached, it was shown that the center of circulation located toward right. For the numerical values of v =0.01($\textrm{cm}^2$/s) and D=0.6($\textrm{cm}^2$/s) and the flow of 50($\textrm{cm}^3$/s ), it was determined that the outflow had to be shut down in 18 seconds to prevent from severe pollution. Also the required time was computed to be 6 seconds for the inflow of 100 ($\textrm{cm}^3$/s). The result of this study is considered, hopefully, to be useful for the design of the water reservoir systems that are the subjects to various contamination.