• 대한조선학회논문집
• /
• 제51권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.

• 대한조선학회논문집
• /
• 제41권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년도 추계 학술대회논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제3권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년도 춘계 학술대회논문집
• /
• 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년도 춘계 학술대회논문집
• /
• 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.

• 대한조선학회논문집
• /
• 제40권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
• /
• 제4권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.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2009년도 학술발표회 초록집
• /
• 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.

• 대한토목학회논문집
• /
• 제13권2호
• /
• pp.267-277
• /
• 1993

본 논문에서는 유한요소법을 2차원 난류경계층에 적용하였으며, 점성유체의 시간의존 비압축성 운동을 시간과 압력장(場)에서 Navier-Stokes방정식과 vorticity방정식을 이용하여 정식 화하였다. 수치계산방법은 Galerkin방법에 기초하였으며, 난류 경계층의 eddy kinematic viscosity에 대해서는 Prandtl의 혼합거리이론을 도입하였다. 난류 경계층에서 파동에 의한 임의 저면에서 저질의 이동을 수치계산하였다. 유한 요소법에 의해 얻어진 결과는 진동흐름에 의한 경계층과 파동에 의한 경계층에서의 특성의 차이를 분명히 하였다.