• Journal of Mechanical Science and Technology
• /
• 제19권5호
• /
• pp.1169-1181
• /
• 2005

Experimental and numerical studies on the unsteady wake field behind a square cylinder near a wall were conducted to find out how the vortex shedding mechanism is correlated with gap flow. The computations were performed by solving unsteady 2-D Incompressible Reynolds Averaged Navier-Stokes equations with a newly developed ${\epsilon}-SST$ turbulence model for more accurate prediction of large separated flows. Through spectral analysis and the smoke wire flow visualization, it was discovered that velocity profiles in a gap region have strong influences on the formation of vortex shedding behind a square cylinder near a wall. From these results, Strouhal number distributions could be found, where the transition region of the Strouhal number was at $G/D=0.5{\sim}0.7$ above the critical gap height. The primary and minor shedding frequencies measured in this region were affected by the interaction between the upper and the lower separated shear layer, and minor shedding frequency was due to the separation bubble on the wall. It was also observed that the position (y/G) and the magnitude of maximum average velocity $(u/u_{\infty})$ in the gap region affect the regular vortex shedding as the gap height increases.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2002년도 제18회 학술발표대회 논문초록집
• /
• pp.35-40
• /
• 2002

2차원, 3차원 비압축성 Navier-Stokes 방정식을 이용하여 DCA 압축기 익렬의 수치해석을 수행하고, 여러 가지 입사각에 대해 실험치와 비교.검토하였다. SIMPLE 알고리즘을 적용한 2차원, 3차원 코드는 대류항의 이산화에 하이브리드 도식을, 진동해를 방지하기 위해 집중격자 기법(PWIM)을 사용하였다. 캐스케이드 유동을 예측하는데 있어서 가장 중요한 요소 중의 하나가 난류모델링이다. 이는 캐스케이드 내의 유동이 역압력구배에 의한 박리와 재부착 등의 복잡한 양상을 보이기 때문이다. 본 연구에서는 계산시간의 효율을 위해 k-$\varepsilon$ 벽법칙 모델을 사용하였다.

• Journal of Advanced Marine Engineering and Technology
• /
• 제30권2호
• /
• pp.275-283
• /
• 2006

Experiment and three dimensional numerical investigations of incompressible turbulent flow through square channels with one- and two-sided ribbed walls are performed to determine pressure drop and heat transfer. The CFX(version 5.7) software package is used for the computation. The ribbed walls have a $45^{\circ}$ inclined square rib. Uniform heat flux is maintained on whole inner heat transfer channel area. The numerical results coincide with experimental data that obtained for $7,600{\le}Re{\le}24.900$, the pitch-to-rib height ratio (p/e) of 8.0. and the rib height-to-channel hydraulic diameter ratio ($e/D_h$) of 0.0667. The results show that values of local heat transfer coefficient and friction factor in the channel with two-sided ribbed wall are higher than those in the channel with one-sided ribbed walls.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1996년도 춘계 학술대회논문집
• /
• pp.99-104
• /
• 1996

The jet pump is being used in many fields for several purposes because of its simple construction and easy operation. The characteristics of the geometrical variables, pressure gradient and velocity distribution of the jet pump are studied using the CFD technique. The flow calculations through a bended nozzle. a mixing chamber and a venturi are presented and phenomenological aspects are discussed. This study solve 3-D steady incompressible Navier-Stokes equations using the Iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives and 3rd-order accurate QUICK scheme for the convective terms. The Mark-and-cell concept was applied efficiently to solve continuity equation, which is differenced 2nd-order accurate central differenced scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A O-type of grid system is generated inside a nozzle and venturi of the jet pump. It has concluded that the results of present study properly agree with physical flow phenomena.

• 대한조선학회논문집
• /
• 제47권1호
• /
• pp.38-46
• /
• 2010

This study aims at investigating the effect of the baffle height on the liquid sloshing in the three-dimensional (3D) rectangular tank. In order to simulate the 3D incompressible viscous two-phase flow in the 3D tank with partially filled liquid, the present study has adopted the volume of fluid (VOF) method based on the finite-volume method which has been well verified by comparing with the results of the relevant previous researches. The ratio of the baffle height ($h_B$) to filling level (h) has been changed in the range of $0{\leq}h_B/h{\leq}1.2$ to observe the effect on the impact loads on the side wall and free surface behavior. Generally, as baffle height increases, the impact pressure on the wall decreases and the deformation of free surface becomes weaker. However it seemed that a critical ratio of the baffle height existed to reveal the lowest impact pressure on the wall. Consequently, $h_B/h=0.8$ among $h_B/hs$ considered in the study showed the lowest impact pressure.

• International Journal of Ocean System Engineering
• /
• 제3권2호
• /
• pp.68-76
• /
• 2013

In this study, a 3D numerical model was used to predict nonlinear wave forces on a cylindrical pile installed in a shallow water region. The model was based on solving the viscous and incompressible Navier-Stokes equations for a two-phase flow (water and air) model and the volume of fluid method for treating the free surface of water. A new application was developed based on the cut-cell method to allow easy installation of complicated obstacles (e.g., bottom geometry and cylindrical pile) in a computational domain. Free-surface elevation, water particle velocities, and inline wave forces were calculated, and the results show good agreement with experimental data obtained by the Danish Hydraulic Institute. The simulation results revealed that the proposed model can, without the use of empirical formulas (i.e., Morison equation) and additional wave analysis models, reliably predict non-linear wave forces on an offshore wind turbine foundation installed in a shallow water region.

• 한국전산유체공학회지
• /
• 제10권2호
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• 대한조선학회논문집
• /
• 제51권5호
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

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

• Ocean Systems Engineering
• /
• 제6권3호
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.