• Journal of Advanced Research in Ocean Engineering
• /
• v.3 no.4
• /
• pp.165-173
• /
• 2017

This study aims at validating simulations of the forced and freely vibrating cylinders at Reynolds number of approximately 500 in order to identify the capability of the CFD code, and to establish the analysis process of the vortex-induced vibration (VIV). The direct numerical and large eddy simulations were employed to resolve the various length scales of the vortices, and the morphing technique was used to consider a motion of the circular cylinder. For the forced vibration case, both in- and anti-phase VIV processes were observed regarding the frequency ratio. Namely, when the frequency ratio approaches to unity, the synchronization/lock-in process occurs, leading to substantial increases in drag and lift coefficients. This is strongly linked with the switch in timing of the vortex formation, and this physical tendency is consistent with that of Blackburn and Henderson (J. Fluid Mech., 1999, 385, 255-286) as well as force coefficients. For the free oscillation case, the mass and damping ratio of 50.8 and 0.0024 were considered based on the study of Blackburn et al. (J. Fluid Struct., 2000, 15, 481-488) to allow the direct comparison of simulation results. The simulation results for a peak amplitude of the cylinder and a shedding mode are reasonably comparable to that of Blackburn et al. (2000). Consequently, based on aforementioned results, it can be concluded that numerical methods were successfully validated and the calculation procedure was well established for VIV analysis with reasonable results.

• Ocean Systems Engineering
• /
• v.4 no.3
• /
• pp.185-200
• /
• 2014

The objective of the present simulations is to evaluate the applicability of the standard $k-{\varepsilon}$ turbulence model in engineering practice in the subcritical to supercritical flow regimes. Two-dimensional numerical simulations of flow around a circular cylinder at $Re=1{\times}10^5$, $5{\times}10^5$ and $1{\times}10^6$, had been performed using Unsteady Reynolds-Averaged Navier Stokes (URANS) equations with the standard $k-{\varepsilon}$ turbulence model. Solution verification had been studied by evaluating grid and time step size convergence. For each Reynolds number, several meshes with different grid and time step size resolutions were chosen to calculate the hydrodynamic quantities such as the time-averaged drag coefficient, root-mean square value of lift coefficient, Strouhal number, the coefficient of pressure on the downstream point of the cylinder, the separation angle. By comparing the values of these quantities of adjacent grid or time step size resolutions, convergence study has been performed. Solution validation is obtained by comparing the converged results with published numerical and experimental data. The deviations of the values of present simulated quantities from those corresponding experimental data become smaller as Reynolds numbers increases from $1{\times}10^5$ to $1{\times}10^6$. This may show that the standard $k-{\varepsilon}$ model with enhanced wall treatment appears to be applicable for higher Reynolds number turbulence flow.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.

• 한국가시화정보학회:학술대회논문집
• /
• 2001.12a
• /
• pp.86-99
• /
• 2001

Vortex lock-on or resonance in the flow behind a circular cylinder is investigated from a time-resolved PIV when a single frequency oscillation is superimposed on the mean incident velocity. Measurements are made of the $K\acute{a}rm\acute{a}n$ and streamwise vortices in the wake-transition regime at the Reynolds number 360. Streamwise vortices at the lock-on and natural shedding states are observed, as well as the changes in the wake region with the change of the shedding frequency of lock-on state. When lock-on occurs, the vortex shedding frequency is found to be half the oscillation frequency as expected from previous experiments. At the lock-on state, the $K\acute{a}rm\acute{a}n$ vortices are observed to be more disordered by the increased strength and spanwise wavelength of the streamwise vortices, which leads to a strong three-dimensional motion. Recirculation and vortex formation region at the lock-on state is reduced as the oscillating frequency is increased. By comparing the Reynolds stresses at the lock-on and natural shedding states, $\bar{u'u'}\;and \;\bar{u'u'}$ at the lock-on state are concentrated on the shear layer around the cylinder. The $\bar{u'u'}\;at\;f_o/f_n=2.0$ has a large value near the centerline, compared with that of other cases. Considering the traces of maximum of u', in the wake region near the cylinder, wake width at the lock-on state is wider than that at the natural shedding state.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1417-1426
• /
• 2001

Turbulent flow characteristics in the near wake of a square cylinder have been studied experimentally by using a Digital PIV method. Experiments are performed at the Reynolds numbers of 1600 and 3900 based on the free-stream velocity and the square height. The ensemble averaged turbulence statistics are acquired from 2030 realizations of instantaneous fluctuating velocity field after the conventional Reynolds decomposition. The differences in turbulent intensity and Reynolds shear stress profiles fur both oases indicate that the effect of Reynolds number seems to be descernible mainly due to the occurrence of transition in the separated shear layer. Because of the periodic nature of vortex shedding process, transverse velocity fluctuations contribute dominantly , to turbulent kinetic energy distribution. A comparison with previous LDV data obtained at much higher Reynolds number shows a fairly good agreement each other. It turns out that the effect of Reynolds number diminishes as increasing Reynolds number, which is a well-known feature of a sharp-edged bluff body wake. The streamwise variation of turbulence intensities are compared with those from a circular cylinder along the centerline at the same Reynolds number. The overall magnitudes and the decay rates of turbulence intensities are quite similar, but some differences are noticeble especially in the transverse intensity variation.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.2
• /
• pp.7-14
• /
• 2011

Numerical simulations were performed for the evaluation of wave and current loads on a fixed cylindrical substructure model for an ocean wind turbine using the ANSYS-CFX package. The numerical wave tank was actualized by specifying the velocity at the inlet and applying momentum loss as a wave damper at the end of the wave tank. The Volume-Of-Fluid (VOF) scheme was adopted to capture the air-water interface. An accuracy validation of the numerical wave tank with a truncated vertical circular cylinder was accomplished by comparing the CFD results with Morison's formula, experimental results, and potential flow solutions using the higher-order boundary element method (HOBEM). A parametric study was carried out by alternately varying the length and amplitude of the wave. As a meaningful engineering application, in the present study, three kinds of conditions were considered, i.e., cases with current, waves, and a combination of current and progressive waves, passing through a cylindrical substructure model. It was found that the CFD results showed reasonable agreement with the results of the HOBEM and Morison's formula when only progressive waves were considered. However, when a current was included, CFD gave a smaller load than Morison's formula.

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

• 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.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.4
• /
• pp.25-34
• /
• 1986

A two-dimensional linear method has been developed for the motion and the second-order steady force arising from the hydrodynamic coupling between waves and currents in the presence of a body of arbitrary shape. Interaction between the incident wave and current in the absence of the body lies in the realm beyond our interest. A Fredholm integral equation of the second kind is employed in association with the Haskind's potential for a steadily moving source of pulsating strength located in or below the free surface. The numerical calculations at the preliminary stage showed a significant fluctuation of the hydrodynamic forces on the surface-piercing body. The problem is approximately solved by using the asymptotic Green function for $U^2{\rightarrow}0$. The original Green function, however, is applied for the fully submerged body. Numerical calculations are made for a submerged and for a half-immersed circular cylinder and extensively for the mid-ship section of a Lewis-form. Some of the results are compared with other analytical results without any available experimental data. The current has strong influence on roll motion near resonance. When the current opposes the waves, the roll response are generally negligible in the low frequency region. The current has strong influence on roll motion near resonance. When the current opposes the wave, the roll response decreases. When the current and wave come from the same direction, the roll response increases significantly, as the current speed increases. The mean drift forces and moment on the submerged body are more affected by current than those on the semi-immersed circular cylinder or on the ship-like section in the encounter frequency domain.

• Wind and Structures
• /
• v.22 no.2
• /
• pp.253-272
• /
• 2016

The paper concerns with the method and results of wind tunnel investigations of the Strouhal number (St) of a stationary iced cable model of cable-supported bridges with respect to different angles of wind attack. The investigations were conducted in the Climatic Wind Tunnel Laboratory of the Czech Academy of Sciences in $Tel{\check{c}}$. The methodology leading to the experimental icing of the inclined cable model was prepared in a climatic section of the laboratory. The shape of the ice on the cable was registered by a photogrammetry method. A section of an iced cable model with a smaller scale was reproduced with a 3D printing procedure for subsequent aerodynamic investigations. The St values were determined within the range of the Reynolds number (Re) between $2.4{\cdot}10^4$ and $16.5{\cdot}10^4$, based on the dominant vortex shedding frequencies measured in the wake of the model. The model was oriented at three principal angles of wind attack for each of selected Re values. The flow regimes were distinguished for each model configuration. In order to recognize the tunnel blockage effect the St of a circular smooth cylinder was also tested. Good agreement with the reported values in the subcritical Re range of a circular cylinder was obtained. The knowledge of the flow regimes of the airflow around an iced cable and the associated St values could constitute a basis to formulate a mathematical description of the vortex-induced force acting on the iced cable of a cable-supported bridge and could allow predicting the cable response due to the vortex excitation phenomenon.