• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.328-333
• /
• 2014

A hybrid particle-mesh method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations is a combination of the Vortex-In-Cell(VIC) method for convection and the penalization method for diffusion. The key feature of the numerical methods is to determine velocity and vorticity fields around a solid body on a temporary grid, and then the time evolution of the flow is computed by tracing the convection of each vortex element using the Lagrangian approach. Assuming that the vorticity and velocity fields are to be computed in time domain analysis, pressure fields are estimated through a complete set of solutions at present time step. It is possible to obtain vorticity and velocity fields prior to any pressure calculation since the pressure term is eliminated in the vorticity-velocity formulation. Therefore, pressure field is explicitly treated by solving a suitable Poisson equation. In this paper, we propose a simple way to numerically implement the vorticity-velocity-pressure formulation including a penalty term. For validation of the proposed numerical scheme, we illustrate the early development of viscous flows around an impulsive started circular cylinder for Reynolds number of 9500.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.47-52
• /
• 2002

We present an improved Lagrangian vortex method in 2-D incompressible unsteady viscous flows, which is based on a mesh-free integral approach of the velocity-vorticity formulation. Vorticity fields are represented by discrete vortex blobs that are updated by the Lagrangian vorticity transport with the particle strength exchange scheme. Velocity fields are expressed in a form of the Helmholtz decomposition, which are calculated by a fast algorithm of the Biot-Savart integration with a smoothed kernel and by a well-established panel method. No-slip condition is enforced through viscous diffusion of vorticity from a solid body into field. The vorticity flux is determined in such a way that spurious slip velocity vanishes. Through the comparison with the existing finite volume scheme for the transient vortical flows around an impulsively started cylinder at Reynolds number Re=550, we would obtain a more accurate scheme for vortex methods in complicated flows.

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

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.37-42
• /
• 2003

A vorticity-velocity integro-differential formulation of incompressible Wavier-Stokes equations is described, focusing on a scheme for calculating pressure fields in application of the Lagrangian vortex method in connection with panel methods. It deals with the dynamic coupling among velocity, vorticity and pressure, and the Helmholtz decomposition of the velocity field, through a comparative study with the Eulerian finite volume method, we provide an extensive understanding of the Lagrangian vortex methods for numerical simulations of viscous flows around arbitrary bodies.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme 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. The 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 present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

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

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.161-165
• /
• 2004

In the vortex particle method based on the vorticity-velocity formulation for solving the Wavier-Stokes equations, the unsteady, incompressible, viscous laminar flow over a NACA 0012 foil is simulated. By applying an operator-splitting method, the 'convection' and 'diffusion' equations are solved sequentially at each time step. The convection equation is solved using the vortex particle method, and the diffusion equation using the particle strength exchange(PSE) scheme which is modified to avoid a spurious vorticity flux. The scheme is improved for variety body shape using one image layer scheme. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsively started NACA 0012 foil for Reynolds number 550.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.2
• /
• pp.21-27
• /
• 2003

The phenomena of vortex shedding around a cylinder oscillating harmonically in a fluid at rest are investigated by a two-dimensional numerical simulation of the Navier-Stokes equations. The simulation is based on a vorticity-velocity integro-differential formulation dealing with vorticity, velocity and pressure variables. Three combinations of Reynolds number(Re) and Keulegan-Carpenter number(KC) were taken to investigate the associated vortex development around the cylinder in the different flow regimes. Drag and lift forces are computed to describe their dominant frequency modulation which is related to the vortex shedding and to the harmonic motion of the cylinder.

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

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