• 한국전산유체공학회지
• /
• 제4권3호
• /
• pp.44-51
• /
• 1999

A modified version of SST turbulence model is suggested to simulate unsteady separated flows over oscillating airfoils. The original SST model, which shows good performance in predicting various steady flows, often results in oscillatory behavior of aerodynamic loads in large separated flow regions. It is shown that this oscillatory behavior is due to the adoption of the absolute value of vorticity in generalizing the original model. As a remedy, a modification is made such that the vorticity in the original SST model is replaced by strain rate. The present model is verified for a mild separated airfoil flow at fixed angle of incidence and for unsteady flowfields about oscillating airfoils. The results are compared with BSL model and original SST model. It is illustrated that the present model gives a better agreement with the experimental results than other two models.

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

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 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.

• 대한기계학회논문집B
• /
• 제28권7호
• /
• pp.809-817
• /
• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• 대한기계학회논문집
• /
• 제17권9호
• /
• pp.2329-2338
• /
• 1993

An accurate analysis procedure to solve the flow about a flat plate at various incidences has been developed. The Navier-Stokes equations of stream function and vorticity form are solved in a sufficiently large computational domain, in which the grid lines are mutually orthogonal. The details of the flow near the singularity at the tip of the plate is well captured by the analytic solution which is asymptotically matched to the numerically generated outer solution. The solution for each region is obtained iteratively : the solution of one (inner or outer) region uses that of the other as the boundary condition after each cycle. The resulting procedure is accurate everywhere and also computationally efficient as the singularity has been removed. It is applied to the flat plate for a wide range of Re : the results agree very well with the existing computation and experiment.

• 소음진동
• /
• 제7권2호
• /
• pp.215-222
• /
• 1997

The flow and acoustic fields due to a vortex ring interaction with a rigid sphere are simulated numerically. The flow field is regarded as three-dimensional inviscid and incompressible. The vorticity is assumed to be concentrated inside the finite core of vortex filament. The vortex filament curve, described by parabolic blending curve function, is used to effectively solve the modified Biot-Savart equation. The interaction between a vortex ring and a rigid sphere using the parabolic blending curve is calculated. The trajectory of the vortex ring is obtained with several different initial positions between the ring and the sphere. The force variations acting on the sphere are calculated by using the boundary integral method. Finally, we can also obtain the acoustic signals at the far field observation positions from the force variations acting on the rigid surface. We can find that the dipole axis of the directivity patterns are rotated during the interacting phenomena.

• 대한토목학회논문집
• /
• 제36권6호
• /
• pp.1037-1048
• /
• 2016

본 연구에서는 Navier-Stokes 방정식과 자유수면 추적에 VOF법을 채용하는 CADMAS-SURF를 이용하여 파랑과 흐름의 비선형 상호간섭현상을 연구하였다. 파-흐름 공존장에서의 유체거동 해석을 위해 CADMAS-SURF를 수정 및 확장하였고, 계산치를 실험치와 비교하여 본 연구의 타당성을 확인할 수 있었다. 본 연구의 수치수로를 파랑-흐름 공존장에서 주어지는 복잡한 물리 현상 규명과 파-흐름-잠제의 상호간섭해석에 적용하여, 유속장, 와도장, 자유수면과 와도와의 관계 등을 논의하였다.

• Journal of Mechanical Science and Technology
• /
• 제14권12호
• /
• pp.1386-1395
• /
• 2000

The mear field structure of round turbulent jets with initially asymmetric velocity distributions is investigated experimentally. Experiments are carried out using a constant temperature hot-wire anemometry system to measure streamwise velocity in the jets. The measurements are undertaken across the jet at various streamwise stations in a range starting from the jet exit plane and up to a downstream location of twelve diameters. The experimental results include the distributions of mean and instantaneous velocities, vorticity field, turbulence intensity, and the Reynolds shear stresses. The asymmetry of the jet exit plane was obtained by using circular cross-section pipes with a bend upstream of the exit. There pipes used here include a straight pipe, and 90 and 160 degree-bend pipes. Therefore, at the upstream of the upstream of the pipe exit, secondary flow through the bend mean streamwise velocity distribution could be controlled by changing the curvature of pipes. The jets into the atmosphere have two levels of initial velocity skewness in addition to an axisymmetric jet from a straight pipe. In case of the curved pipe, a six diameter-long straight pipe section follows the bend upstream of the exit. The Reynolds number based on the exit bulk velocity is 13,400. The results indicate that the near field structure is considerably modified by the skewness of an initial mean velocity distribution. As the skewness increases, the decay rate of mean velocity at the centerline also increases.

• 한국전산유체공학회지
• /
• 제9권2호
• /
• pp.23-29
• /
• 2004

A numerical simulation of an incompressible cavity flow is conducted using the hybrid turbulence model. The model adopted is a modified type of DES using k- ε two-equation model. Cavity geometry and flow condition are based on Cattafesta's experiment. Computational results are compared with the results of Cattafesta's experiment. The simulation successfully predicts the oscillatory features and the Strouhal number of the oscillation compares very favorably with that of the dominant mode of experimental data. Vorticity contours obtained from the simulation data are consistent with the smoke visualization of the Cattafesta's experiment. The coherent structures of cavity flow are also investigated using Q criterion.