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

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

• 대한조선학회논문집
• /
• 제35권4호
• /
• pp.1-10
• /
• 1998

본 논문에서는 비압축성 Newtonian 점성유동에서 초기에 순간 출발하는 2차원 실린더 주위의 유동을 해석하기 위해서, 와도를 기저로 한 수치해석기법을 제안하고 있다. Helmholtz 분리 형태로 표현된 Navier-Stokes방정식에서 유도되는 와도전달방정식과 압력방정식, 그리고 벡터등식에서 유도되는 속도-와도 관계식을 이 문제의 지배방정식으로 택하고, 경계조건으로는 물체표면에서 와도와 압력의 연성관계와 힘의 평형을 고려한 동적와도경계조건과 동적압력조건이 제시된다. 이 지배방정식과 경계조건을 수치적으로 처리하기 위하여, 와도와 압력이 연성되어 있는 경계조건은 Wu등(1994)이 제안한 대로, 연성관계를 유지한 채로 식을 분리하는 방법을 이용하였고, 와도전달 방정식은 유한체적법으로 계산하였으며, 그 식에 포함된 대류항을 처리하는 방법으로 TVD 방법을 이용하였다. 속도는 Biot-Savart적분항이 포함된 벡터등식에서 panel방법으로 구하고, 압력방정식은 형태가 Poisson방정식이므로 역시 panel방법을 이용하였다. 계산에 사용된 격자로 정규격자를 이용하고, 결과를 다른 수치적, 해석적 결과와 비교하여 그 타당성을 검증하였다.

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

• 한국전산유체공학회지
• /
• 제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).

• 한국군사과학기술학회지
• /
• 제17권4호
• /
• pp.551-559
• /
• 2014

The influence of the battle damage hole on the velocity and vorticity flow field have been studied by using particle image velocimetry. Time averaged velocity and vorticity vector fields in the vicinity of jet are presented. The perforated damage hole on a wing created from a hit by anti-air artillery was modeled as a 10% chord size hole which positioned at quarter chord. At low angles of attack, the vorticity in the forward side of the jet is cancelled due to mixing with the wing surface boundary layer. Stretching of vorticity in the backside of the jet generates a semi-cylindrical vortical layer that enclosing a domain with slow moving reverse flow. Conversely, at higher the angles of attack, the jet vorticity advected away from the wing surface and remains mostly confined to the jet. The mean flow behind the jet has a wake-like structure.

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

• 대한기계학회논문집B
• /
• 제28권9호
• /
• pp.1032-1041
• /
• 2004

In this study we decompose the 3-dimensional velocity field of isotropic turbulent flow into the coherent and the incoherent structure using the discrete wavelet. It is shown that the coherent structure, 3% wavelet modes, has 98% energy and 88% enstrophy and its statistical characteristics are almost same as the original turbulence structure. And it is confirmed that the role of the coherent structure is that it produces the turbulent kinetic energy at the inertia range then transfers energy to the dissipation range. The incoherent structure, with residual wavelet modes, is uncorrelated and has the Gaussian probability density function but it dissipates the kinetic energy in dissipation range. On the procedure, we propose a new but easy way to get the threshold by applying the energy partition percentage concept about coherent structure. The vorticity field extracted from the wavelet-decomposed velocity field has the same structure as the result of the precedent studies which decomposed vorticity field directly using wavelet. Therefore it has been shown that velocity and vorticity field are on the interactive condition.

• 한국해양공학회지
• /
• 제22권3호
• /
• pp.24-33
• /
• 2008

In this study, flow and creation and dissipation of vorticity around rectangular floating breakwater is investigated both experimentally and numerically. The PIV system(Particle image velocimetry) is employed to obtain the velocity field in the vorticity of rectangular structure. The numerical model, combined with ${\kappa}-{\varepsilon}$ turbulence model and the VOF method based on RANS equation, is used to analyze the turbulence structure. In the results of this study, the vorticity is found around conner of rectangular structure at all time domain, and creation and dissipation of vorticity are closely related to wave period. Separation points of phase of vortex due to flow separation for longer period waves are faster then for shorter period waves.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1601-1606
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.