• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.93-101
• /
• 1986

The behavior of a plane buoyant jet discharged vertically upward into a stagnant uniform environment is analyzed by continuity, momentum transport equation by numerical scheme. The governing equations are solved by finite difference method employing stream function and vorticity transport and Prandtl's turbulent model. Results for centerline velocities and temperatures, temperature distribution and flow pattern in receiving environment due to buoyant jet in the range of discharge densimetric Froude number of 4 to 32 show good agreement with published data. Spreading rate and dispersion ratio, which are required in integral type analysis of whole range of buoyant jet and have not been obtained yet, are derived in terms of discharge densimetric Froude number and vertical distance from source.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.8 no.3
• /
• pp.185-194
• /
• 1984

The finite analytic method is extended to solve the steady two dimensional Navier-Stokes equation of stream functions and vorticity in polar coordinate system. The method is applied to calculate laminar flows in a sector cavity where the motion is induced by the rotation of the outer wall. Numerical solutions are obtained in the range of Reynolds number 0 to 5000 and aspect ratios 0.50, 1.20, 1.60 and 1.92. The finite analytic method is verfied to be accurate and fast convergent at high Reynolds numbers. It is promising as a numerical method of viscous flows and heat transfer. Flows in sector cavities show different flow structures and formation of secondary vortex with aspect ratios and Reynolds numbers in comparison with rectangular cavities.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.3
• /
• pp.595-609
• /
• 2015

The wake characteristics of Contra-Rotating Propeller (CRP) were investigated using numerical simulation and flow measurement. The numerical simulation was carried out with a commercial CFD code based on a Reynolds Averaged Navier-Stokes (RANS) equations solver, and the flow measurement was performed with Stereoscopic Particle Image Velocimetry (SPIV) system. The simulation results were validated through the comparison with the experiment results measured around the leading edge of rudder to investigate the effect of propeller operation under the conditions without propeller, with forward propeller alone, and with both forward and aft propellers. The evolution of CRP wake was analyzed through velocity and vorticity contours on three transverse planes and one longitudinal plane based on CFD results. The trajectories of propeller tip vortex core in the cases with and without aft propeller were also compared, and larger wake contraction with CRP was confirmed.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.1-1
• /
• 2009

In this presentation, I talk about various fluid simulation methods that have been developed for computer graphics special effects since 1996. They are all based on CFD but sacrifice physical reality for visual plausability and time. But as the speed of computer increases rapidly and the capability of GPU (graphics processing unit) improves, methods for more physical realism have been tried. In this talk, I will focus on four aspects of fluid simulation methods for computer graphics: (1) particle level-set methods, (2) particle-based simulation, (3) methods for exact satisfaction of incompressibility constraint, and (4) GPU-based simulation. (1) Particle level-set methods evolve the surface of fluid by means of the zero-level set and a band of massless marker particles on both sides of it. The evolution of the zero-level set captures the surface in an approximate manner and the evolution of marker particles captures the fine details of the surface, and the zero-level set is modified based on the particle positions in each step of evolution. (2) Recently the particle-based Lagrangian approach to fluid simulation gains some popularity, because it automatically respects mass conservation and the difficulty of tracking the surface geometry has been somewhat addressed. (3) Until recently fluid simulation algorithm was dominated by approximate fractional step methods. They split the Navier-Stoke equation into two, so that the first one solves the equation without considering the incompressibility constraint and the second finds the pressure which satisfies the constraint. In this approach, the first step introduces error inevitably, producing numerical diffusion in solution. But recently exact fractional step methods without error have been developed by fluid mechanics scholars), and another method was introduced which satisfies the incompressibility constraint by formulating fluid in terms of vorticity field rather than velocity field (by computer graphics scholars). (4) Finally, I want to mention GPU implementation of fluid simulation, which takes advantage of the fact that discrete fluid equations can be solved in parallel.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.3
• /
• pp.307-316
• /
• 2012

Numerical simulations of 3D aircraft configurations are performed in order to understand the effects of turbulence models on the prediction of aircraft's aerodynamic characteristics. An in-house CFD code that solves 3D RANS equations and two-equation turbulence model equations are used. The code applies Roe's approximated Riemann solver and an AF-ADI scheme. Van Leer's MUSCL extrapolation with van Albada's limiter is also adopted. Various versions of Menter's $k-{\omega}$ SST turbulence models as well as Coakley's $q-{\omega}$ model are incorporated into the CFD code. Menter's $k-{\omega}$ SST models include the standard model, the 2003 model, the model incorporating the vorticity source term, and the model containing controlled decay. Turbulent flows over a wing are simulated in order to validate the turbulence models contained in the CFD code. The results from these simulations are then compared with computational results from the $3^{rd}$ AIAA CFD Drag Prediction Workshop. Numerical simulations of the DLR-F6 wing-body and wing-body-nacelle-pylon configurations are conducted and compared with computational results of the $2^{nd}$ AIAA CFD Drag Prediction Workshop. Aerodynamic characteristics as well as flow features are scrutinized with respect to the turbulence models. The results obtained from each simulation incorporating Menter's $k-{\omega}$ SST turbulence model variations are compared with one another.

• 한국가시화정보학회:학술대회논문집
• /
• 2003.11a
• /
• pp.39-42
• /
• 2003

Leading edge extension(LEX) in a highly swept shape applied to a delta wing features the modern air-fighters. The LEX vortices generated upon the upper surface of the wing at high angle of attack enhance the lift force of the delta wing by way of increased negative suction pressure over the surfaces. The present 3-D stereo PIV includes the Identification of 2-D cross-correlation equation, stereo matching of 2-D velocity vectors of two cameras, accurate calculation of 3-D velocity vectors by homogeneous coordinate system, removal of error vectors by a statistical method followed by a continuity equation criterion and so on. A delta wing model with or without LEX was immersed in a circulating water channel. Two high-resolution, high-speed digital cameras$(1280pixel\times1024pixel)$ were used to allow the time-resolved animation work. The present dynamic stereo PIV represents the complicated vortex behavior, especially, in terms of time-dependent characteristics of the vortices at given measuring sections. Quantities such as three velocity vector components, vorticity and other flow information can be easily visualized via the 3D time-resolved post-processing to make the easy understanding of the LEX effect or vortex emerging and collapse which are important phenomena occurring in the field of delta wing aerodynamics.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1672-1677
• /
• 2004

Leading edge extension(LEX) in a highly swept shape applied to a delta wing features the modem air-fighters. The LEX vortices generated upon the upper surface of the wing at high angle of attack enhance the lift force of the delta wing by way of increased negative suction pressure over the surfaces. The present 3-D stereo PIV includes the Identification of 2-D cross-correlation equation, stereo matching of 2-D velocity vectors of two cameras, accurate calculation of 3-D velocity vectors by homogeneous coordinate system, removal of error vectors by a statistical method followed by a continuity equation criterion and so on. A delta wing model with or without LEX was immersed in a circulating water channel. Two high-resolution, high-speed digital cameras($1280pixel{\times}1024pixel$) were used to allow the time-resolved animation work. The present dynamic stereo PIV represents the complicated vortex behavior, especially, in terms of time-dependent characteristics of the vortices at given measuring sections. Quantities such as three velocity vector components, vorticity and other flow information can be easily visualized via the 3D time-resolved post-processing to make the easy understanding of the LEX effect or vortex emerging and collapse which are important phenomena occurring in the field of delta wing aerodynamics.

• Journal of Advanced Marine Engineering and Technology
• /
• v.36 no.8
• /
• pp.1091-1096
• /
• 2012

The aerodynamic characteristics of circular cylinder in a channel are studied to make clear the flow feature by solving the Navier-Stokes equation based on the finite volume method with unstructured grids. Reviews are made on with the vorticity, velocity, dynamic pressure, residual and drag, where the Reynolds numbers are 50 and 100. The flows for $Re{\succeq}50$ shows the vortex shedding in the wake, and the result is the same as the case of moving cylinder. The ground effect of flat bottom results in the growth of vortex, being generated in the upper side of the cylinder and elongated in the rear. As the cylinder approaches to wall, for example 0.6, the cylinder plays as a role of blockage to obstruct the flow between the cylinder and wall. The drag coefficients are compared with others' results to confirm the validity of the present numerical simulation.

• The Journal of the Acoustical Society of Korea
• /
• v.21 no.2
• /
• pp.134-141
• /
• 2002

In this paper, a free wake analysis based on the curved vortex element and CVC wake model is used to predict the aerodynamic performance and noise for HAWT. Also for prediction of RPM, a maximum value through a quadratic regression was suggested. And for a noise prediction, the broadband noise prediction method based on experimental equation was used. The curved vortex element uses a BCVE and an SIVE instead of a straight vertex element. In the CVC wake model, the vortex strengths are assumed to be constant along a span and a vortex filament. The free wake structure made by the curved vortex element and CVC was substituted for a vortex lattice, so it has an advantage for the less calculation time and a depiction of accurate wake structure. For the verification of this program, calculated results are compared with Mr. Kim's experiment model and Zond Z-40FS for performance and with WTS-4 and USWP models for noise. Good agreements are obtained between the predicted and the measured data for the performance and far-field noise spectra.