• Journal of computational fluids engineering
• /
• v.4 no.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.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.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.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.05a
• /
• 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).

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.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.

• Journal of KSNVE
• /
• v.7 no.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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.36 no.6
• /
• pp.1037-1048
• /
• 2016

In this study, nonlinear wave interaction in the presence of a uniform current is studied using numerical model, named CADMAS-SURF which is based on the Navier-Stokes equations coupled with Volume of Fluid for tracking free surface deformation. The original CADMAS-SURF developed for interaction of wave with structure is modified/extended to simulate nonlinear fluid dynamic motions within wave-current coexisting field. The capability of Numerical Wave-Current Tank (NWCT) in this study is validated by comparing with available existing laboratory experiments for both wave-following and wave-opposing current. The numerical results for interaction between wave and current are shown to be in good agreement with experimental data. Then, this study focused on the dynamic motions of the water velocity, surface elevation and vorticity within combined wave-current field in demonstrating complex nonlinear physical phenomena due to interaction between wave and current. In addition, NWCT is applied to simulate a more complex wave-current-structure field for wave propagating over a submerged breakwater associated with current. Detailed discussion including characteristics of velocity and vorticity fields and the relation between free surface and vorticity are given.

• Journal of Mechanical Science and Technology
• /
• v.14 no.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.

• Journal of computational fluids engineering
• /
• v.9 no.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.