• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.44-48
• /
• 2009

In this study, a vortex visualization method based on the vorticity magnitude is developed. One of the simplest models for a vortex is a vortex filament with the maximum vorticity on its center. The proposed method is based on the observation of this ideal distribution of vorticity magnitude. Laplacian and Hessian matrix of vorticity magnitude are tested for detecting the local maximum of vorticity magnitude. These ideas were applied to wake flow past a sphere. It was found that the Laplacian method is not able to distinguish vortices from the underlying shear layer clearly, while the Hessian matrix method does not suffer from this problem.

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

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

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

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

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• 한국연소학회:학술대회논문집
• /
• 2002.11a
• /
• pp.49-56
• /
• 2002

A two-dimensional direct numerical simulation is performed to investigate the formation characteristics of a single vortex interacting with $CH_4/N_2$-Air counterflow nonpremixed flame. The numerical method was based on a predictor-corrector scheme for a low Mach number flow. The detailed transport properties and a 16-step augmented reduced mechanism are adopted in this calculation. The budgets of the vorticity transport equation arc examined to reveal the mechanisms leading to the formation, evolution and dissipation of a single vortex interacting with counterflow nonpremixed flame. It is found that the stretching term, which depends on the azimuthal component of vorticity, and radial velocity, mainly generates vortieitv in non-reacting and reacting flows. The viscous and baroclinic torque term destroy the vorticity in non-reacting flow. In addition, the baroclinic torque term due to density and pressure gradient generates vorticity, while viscous and the volumetric expansion terms due to density gradient destroy vorticity in reacting flow.

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