• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.72-78
• /
• 2014

A hybrid particle-mesh method as the combination between the Vortex-In-Cell (VIC) method and penalization method has been achieved in recent years. The VIC method, which is based on the vorticity-velocity formulation, offers particle-mesh algorithms to numerically simulate flows past a solid body. The penalization method is used to enforce boundary conditions at a body surface with a decoupling between body boundaries and computational grids. The main advantage of the hybrid particle-mesh method is an efficient implementation for solid boundaries of arbitrary complexity on Cartesian grids. However, a numerical simulation of flows in large domains is still not too easy. In this study, a multi-domain approach is thus proposed to further reduce computation cost and easily implement it. We validate the implementation by numerical simulations of an incompressible viscous flow around an impulsively started circular cylinder.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.1
• /
• pp.23-30
• /
• 2004

In the Lagrangian vortex particle method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations, a numerical scheme for calculating pressure fields is presented. Implementation of the numerical method is directly connected with the well-established surface panel methods, just by dealing with the dynamic coupling among vorticity field. Assuming the vorticity and the velocity fields are to be calculated in time domain analysis, the pressure calculation for a complete set of solution at present time step is performed in a similar way to the one used in the Eulerian description. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsive started circular cylinder for Reynolds number 550. The comparative study with the Eulerian finite Volume method provides an extensive understanding and application of the mesh-free Lagrangian vortex methods for numerical simulation of viscous flows around arbitrary bodies of general shape.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.262-270
• /
• 1987

A numerical method is presented which can solve the unsteady momentum and thermal boundary layers, coupled through the agency of buoyancy force, over a heated circular cylinder impulsively started from rest. By linearizing the nonlinear finite difference equations without sacrificing accuracy, numerical solutions are obtained at each time step without iteration. To get rid of the requirement of excessive number of grid points in the region of reversed flow, special form of transformed variables are used, by which the computational boundary layer thickness is maintained almost constant. These numerical properties enable the method to easily handle the region of reversed flow and how the singularity develops in the interior of the boundary layer. In order to investigated the thermal effects on the skin friction, heat flux, displacement thickness and on the separation, we have successfully solved three different cases of the buoyancy parameter .alpha.(Gr/Re$^{2}$).

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.518-523
• /
• 2001

Vortex methods were originally conceived as a tool to model the evolution of unsteady, incompressible, high Reynolds number flows of engineering interest. Recently various methods have been proposed for simulating the diffusion in vortex methods for two-dimensional incompressible flows. We test the diffusion schemes of vortex methods. In this paper we directly compare the particle strength exchange scheme with the vorticity redistribution scheme in tenus of their accuracy and computational efficiency. Comparisons between both viscous models described are presented for short-time runs of impulsively started flows past a circular cylinder for Reynolds number of 60. The particle strength exchange scheme has been shown more accurate and efficient than the vorticity redistribution scheme.

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.3
• /
• pp.12-18
• /
• 2012

The Lattice Boltzmann Method has developed for solving the Boltzmann equation in Cartesian domains containing immersed boundaries of arbitrary geometrical complexity moving with prescribed kinematics. When a immersed boundaries are sweeping the fixed fluid node, refilling the node information in a vicinity of fluid nodes is one of the important issues in Lattice Boltzmann Method. In this study, we propose a simple refill algorithm for the particle distribution function based on a proper velocity, density and strain rate to enhance accuracy and stability of the method. The refill scheme based on a asymptotic analysis of LBGK model has improved accuracy than interpolation schemes. The proposed scheme in this study is validated by the simulations of an impulsively started rotating circular cylinder to investigate adaptability for fluid-structure interaction (FSI) problem. This refill scheme has improved stability and accuracy especially at high Reynolds number region.