• Journal of computational fluids engineering
• /
• v.14 no.3
• /
• pp.76-85
• /
• 2009

The shock wave is deformed and the vortex is elongated simultaneously during the shock-vortex interaction. More precisely, the shock wave is deformed to a S-shape, consisting of a leading shock and a lagging shock by which the corresponding local vortex flows are accelerated and decelerated, respectively: the vortex flow swept by the leading shock is locally expanded and the one behind the lagging shock is locally compressed. As the leading shock escapes the vortex in the order of microseconds, the expanded flow region is quickly changed to a compression region due to the implosion effect. An induced shock is developed here and propagated against the vortex flow. This happens for a strong vortex because the tangential flow velocity of the vortex core is high enough to make the induced-shock wave speed supersonic relative to the vortex flow. For a weak shock, the vortex is basically subsonic and the induced shock wave is absent. For a vortex of intermediate strength, an induced shock wave is developed in the supersonic region but dissipated prematurely in the subsonic region. We have expounded these three shock-vortex interaction patterns that depend on the vortex flow regime using a third-order ENO method and numerical shadowgraphs.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.06a
• /
• pp.414-419
• /
• 2007

This paper describes the prediction of aerodynamic performance and wake of HAWT in normal and yawed flow operation using potential based methods. In order to analyze aerodynamic performance of wind turbine WINFAS program is used, which is based on VLM(Vortex Lattice Method) and CVC(Constant vorticity contour) Free wake model. Some problems of CVC vortex filament method are investigated arid to improve these problems vortex ring wake are introduced in behalf of CVC vortex filament. The prediction results using the vortex lattice wake are compared to experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.3
• /
• pp.62-74
• /
• 1993

The vortex shedding from a circular cylinder placed in a steady uniform stream is simulated by the vortex cloud model of the discrete vortex method. The vorticity created at the cylinder surface is discretely represented by a number of nascent vortices at each time step and the motion of these cumulative vortices is monitored to produce the evolution of the vortex distribution pattern. Convection of vortices was traced by the vortex-in-cell technique and the force coefficients were calculated by both Sarpkaya's formulae and Lee's formulae for comparison. Discussions concerning the interrelation between the computational parameters and some principles for choosing the suitable values are included.

• International Journal of Aeronautical and Space Sciences
• /
• v.6 no.2
• /
• pp.84-96
• /
• 2005

Helical tip vortex is known as stable vortex structure, however the specific frequency component of far wake perturbation induces the vortex pairing in hover and axial flight. It is expected that the tip vortex pairing phenomena may happen in transition flight and very low advance ratio flight so that inflow may be most nonuniform in the low advance ratio flight. The objectives of this paper are that a tip-vortex instability during the transition from hover into very low advance ratio forward flight is numerically predicted to understand a physics by using a time-marching free-wake method. To achieve the objectives, numerical method is firstly validated in typical axial and forward flights cases. Present scheme with trim routine can predict airloads and inflow distribution of forward flight with good accuracy. Then, the transition flight condition is calculated. The rotor used in this wake calculation is a small-scale AH-1G model. By using a tip-vortex trajectory tracking method, the tip-vortex pairing process are clearly observed in transient flight($\mu$=0.03) and disappears at a slightly higher advance ratio($\mu$=0.05). According to the steady flight simulation at $\mu$=0.03, it is confirmed the tip-vortex pairing process is continued in the rear part of rotor disk and not occurs in the front part. Time averaged inflow in this case is predicted as smooth distribution.

• Journal of Institute of Convergence Technology
• /
• v.4 no.1
• /
• pp.37-41
• /
• 2014

For the analysis of lifting surface at high angle of attack, a Nonlinear Vortex Lattice Method(NLVLM) was used. The NLVLM is intented to compute the interactions between lifting surfaces and separated vertical flow. The lifting surfaces are represented by a lattice of discrete vortex rings. And wakes are represented by families of non-lintersecting, semi-infinite vortex line segments. The image method also used to analyze the ground effect. It is found that vortex lines separated from lifting surfaces represent the separated flows successfully. Although the present method is applied for the rectangular wing and delta wing, extensions can be possible for the arbitrary lifting surfaces. The Present results show good agreement with experimental data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.5
• /
• pp.1264-1271
• /
• 1990

The vortex lattice method was adopted to predict the aerodynamic performance of a horizontal axis wind turbine. For this simulation. the rotor blade was divided into many panels both in chordwise and spanwise direction and then replaced by horseshoe vortices. The wake was divided into two parts of near wake and far wake : the near wake was assumed as helical vortex line elements and the far wake was modeled by semi-infinite circular vortex cylinder. The induced velocity components were calculated by the Biot-Savart law. By this way the power coefficient was obtained and represented as a function of the tip speed ratio. The numerical results obtained were compared with those of the other methods and experimental results and showed good agreement with experimental results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.11
• /
• pp.1039-1048
• /
• 2008

new computational procedure for the Non-Linear Vortex Lattice Method (NLVLM) is suggested in this work. Conventional procedures suggested so far usually involves inner iteration loop to update free vortex shape and an under-relaxation based iteration loop to determine the free vortex shape. In this present work, we suggest a new formula based on quasi-steady concept to fix free vortex shape which eliminates the need for inner iteration loop. Further, the ensemble averaging of the induced velocities for a given free vortex segment evaluated at each iteration significantly improves the convergence property of the algorithm without resorting to the under-relaxation technique. Numerical experiments over several low aspect ratio wings are carried out to obtain optimal empirical parameters such as the length of the free vortex segment, the vortex core radius, and the rolled-up wake length.

• Journal of Aerospace System Engineering
• /
• v.10 no.1
• /
• pp.1-7
• /
• 2016

During flight of the aircraft, the vortex merging phenomenon appears under the certain condition between co-rotating vortices which were generated at the wing tip and lifting-surface. And then these merged vortices at both sides show counter-rotating pattern to affect on the downstream of the aircraft. In this paper, the numerical simulations are conducted assuming this phenomenon in two-dimensional co-rotating or counter-rotating vortices pairs. Two-dimensional incompressible Navier-Stokes equations were converted into Vorticity-Streamfunction form and the Galerkin spectral method was adopted. The third order Runge-Kutta method was used for time integration. The effects on the vortex merger and degree of vortex merger were investigated according to time, Reynolds number, and changes in the distance between two vortices.

• Journal of the Society of Naval Architects of Korea
• /
• v.46 no.3
• /
• pp.249-258
• /
• 2009

DVM (Discrete Vortex Method) is a numerical scheme that handles discrete vortex particles to express continuous vorticity field. This scheme is proper to VIV (Vortex Induced Vibration) analysis because there is no need to generate field grids and VIV is caused by separated vorticity from the body. When DVM is applied to VIV analysis, there are some applicable schemes such as using vortex blobs, integral method for computing induced velocity, etc. In this study, the influences of these schemes are investigated and the practical scheme that is appropriate for VIV analysis is proposed.

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