• Wind and Structures
• /
• v.29 no.1
• /
• pp.43-53
• /
• 2019

Vortex-induced vibration of three circular cylinders (each of diameter D) in an equilateral triangular arrangement is investigated using the immersed boundary method. The cylinders, with one placed upstream and the other two side-by-side downstream, are free to vibrate in the cross-flow direction. The cylinder center-to-center spacing L is adopted as L/D = 2.0. Other parameters include the Reynolds number Re = 100, mass ratio $m^*=2.0$, reduced velocity $U_r=2{\sim}15$ and damping ratio ${\zeta}=0$. Cylinder vibration responses are dependent on $U_r$ and classified into five regimes, i.e. Regime I ($U_r{\leq}3.2$), Regime II ($3.2<U_r{\leq}5.0$), Regime III ($5.0<U_r{\leq}6.4$), Regime IV ($6.4<U_r{\leq}9.2$) and Regime V ($U_r>9.2$). Different facets of vibration amplitude, hydrodynamic forces, wake patterns and displacement spectra are extracted and presented in detail for each regime.

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

In the present study, a least square weighted residual method and Taylor-Galerkin method were formulated and tested for the discretization of the two hyperbolic type equations of level set method; advection and reinitialization equations. The two approaches were compared by solving a time reversed vortex flow and three-dimensional broken dam flow by employing a four-step splitting finite element method for the solution of the incompressible Navier-Stokes equations. From the numerical experiments, it was shown that the least square method is more accurate and conservative than Taylor-Galerkin method and both methods are approximately first order accurate when both advection and reinitialization phase are involved in the evolution of free surface.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.6
• /
• pp.454-460
• /
• 2009

To analyze the mixture formation process of evaporating diesel spray is important for emissions reduction in actual engines. Then the effects of change in density of ambient gas on spray structure in high temperature and pressure field have been investigated in this study. The ambient gas density was changed from ${\rho}_a=5.0kg/m^3$ to ${\rho}_a=12.3kg/m^3$ with CVC(Constant Volume Chamber). Also, simulation study by modified KIVA-II code was conducted and compared with experimental results. The ambient temperature and injection pressure are kept as 700K and 72MPa, respectively. The images of liquid and vapor phase in the evaporating free spray were simultaneously taken by exciplex fluorescence method. As experimental results, with increasing ambient gas density, the tip penetration of the evaporating free spray decreases due to the increase in the drag force from ambient gas. The spatial structure of a diesel spray can be verified as 2-regions consisted of liquid with momentum decrease and vapor with large-scale vortex. The calculated results obtained by modified KIVA-II code show good agreements with experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.11
• /
• pp.709-717
• /
• 2017

In this study, a model of two-stage Weis-Fogh type water turbine model is proposed, the hydrodynamic characteristics of this water turbine model are calculated by the advanced vortex method. The basic conditions and the motion of each wing are the same as that of the single-stage model previously proposed by the same author. The two wings (NACA0010 airfoils) and both channel walls are approximated by source and vortex panels, and free vortices are introduced from the body surfaces. The distance between the front wing axis and the rear wing axis, and the phase difference between the motion of the two wings, which is in phase and out of phase are set as the calculation parameters. For each case, the unsteady flow fields, pressure fields, force coefficients, and efficiency of the two wings are calculated, and the hydrodynamic characteristics of the proposed water turbine model are discussed.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.10 no.2
• /
• pp.107-111
• /
• 2007

A two-dimensional numerical method for inviscid two-fluid flows with evolution of density interface is developed, and an initially stationary two-dimensional fluid sheet surrounded by another fluid is studied. The interface between two fluids is modeled as a vortex sheet, and the flow field with the evolution of interface is solved by using vortex-in-cell/front-tracking method. The edge of the sheet is pulled back into the sheet due to surface tension and a blob is formed at the edge. This blob and fluid sheet are connected by a thin neck. In the inviscid limit, such process of the blob and neck formation is examined in detail and their kinematic characteristics are summarized with dimensionless parameters. The edge recedes at and the capillary wave propagating into the fluid sheet must be considered for better understanding of the edge receding.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.2
• /
• pp.156-162
• /
• 2008

The static aeroelastic analysis of composite hingeless rotor blades in hover was performed using free-wake method. Large deflection beam theory was applied to analyze blade motions as a one-dimension beam. Anisotropic beam theory was applied to perform a cross-sectional analysis for composite rotor blades. Aerodynamic loads were calculated through a three-dimensional aerodynamic model which is based on the unsteady vortex lattice method. The wake geometry in hover was described using a time-marching free-wake method. Numerical results of the steady-state deflections for the composite hingeless rotor blades were presented and compared with those results based on two-dimensional quasi-steady strip theory and prescribed wake method. It was shown that wakes affect the steady-state deflections.

• The KSFM Journal of Fluid Machinery
• /
• v.13 no.4
• /
• pp.31-36
• /
• 2010

Free energy based lattice Boltzmann method (LBM) has been used to simulate the rising bubble flows with large density ratio. LBM with compact discretization is able to reduce the spurious current of the static bubble test and be satisfied with the Laplace law. The terminal rise velocity and shape of the bubbles are dependent on Eotvos number, Morton number and Reynolds number. For single bubble flows, simulations are executed for various Eotvos number, Morton number and Reynolds number, and the results are agreed well with the experiments. For multiple bubbles, the bubble flow characteristics are related by the vortex pattern of the leading bubble. The coalescence of the bubbles are simulated successfully and the subsequent results are presented. The present method is validated for static, dynamic bubble test cases and compared to the numerical, experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.133-140
• /
• 1998

A Navier-Stokes code based on a unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number $k-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for a domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. Unsteadiness inside boundary layer is entrained when a unsteady vortex impinge on the blade surface. It shoves that local peak value inside the boundary layer and also local minimum near the edge of boundary layer as it developes along the blade surface. The unsteadiness inside the boundary layer is almost isolated from the free stream unsteadiness and being convected at local boundary layer speed, less than the free stream value.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.777-784
• /
• 2000

A Navier-Stokes code based on an unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number ${\kappa}-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for the whole experimental domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. The first harmonics of the velocity in the boundary layer shows local peak value inside the boundary layer and also local minimum near the edge of boundary layer. It is intensified as it develops along the blade surface. This is shown to be caused as the unsteadiness inside the boundary layer is being convected at a speed less than the free stream value. It is also shown that there is negligible mixing of the unsteadiness between the boundary layer and the free stream.

• Journal of Navigation and Port Research
• /
• v.37 no.1
• /
• pp.1-7
• /
• 2013

In this paper the developed prediction technique of wave-making resistance performance for a ship attached with a vertical blade had been verified. Numerical analysis program as a prediction technique had been developed using the Rankine source panel method and the vortex lattice method(VLM). The nonlinearity of the free surface conditions was fully taken into account using the iterative method and the trim and the sinkage of the ship were also considered in the numerical analysis program. Panel cutting method was applied to get hull surface panels. Numerical computations were carried out for a 4000TEU container carrier and the vertical blade was attached 6 different locations astern. To investigate the validity of the numerical analysis program the commercial viscous flow field analysis program FLUENT was used to obtain the viscous flow field around the ship and the model test was performed. The model test results were compared with the numerical analysis results.