• 한국안전학회지
• /
• 제10권2호
• /
• pp.56-63
• /
• 1995

A simulation of contaminant dispersion in a water reservoir has been done using 2-D finite difference method(FDM). The steady state velocity field of the reservoir was computed using stream function-vorticity formulation of Wavier-Stokes equation and continuity equation. Based on the computed steady state velocity field, the transient convective diffusion equation of the contaminant dispersion was computed. For the 1m$\times$1m reservoir model with inlet and outlet attached, it was shown that the center of circulation located toward right. For the numerical values of v =0.01($\textrm{cm}^2$/s) and D=0.6($\textrm{cm}^2$/s) and the flow of 50($\textrm{cm}^3$/s ), it was determined that the outflow had to be shut down in 18 seconds to prevent from severe pollution. Also the required time was computed to be 6 seconds for the inflow of 100 ($\textrm{cm}^3$/s). The result of this study is considered, hopefully, to be useful for the design of the water reservoir systems that are the subjects to various contamination.

• 대한기계학회논문집B
• /
• 제24권5호
• /
• pp.733-743
• /
• 2000

We solve the integral representation of the Navier-Stokes equations in a lagrangian view by tracking the particles, which have vortex strengths. We simulate the unsteady viscous flow around an impulsively started cylinder using the vortex particle method. Particles are advanced via the Biot-Savart law for a lagrangian evolution of particles. The particle strength is modified based on the scheme of particle strength exchange. The solid boundary satisfies the no-slip boundary condition by the vorticity generation algorithm. We newly modify the diffusion scheme and the boundary condition for simulating an unsteady flow efficiently. To save the computation time, we propose the mixed scheme of particle strength exchange and core expansion. We also use a lot of panels to ignore the curvature of the cylinder, and not to solve the evaluation of the surface density. Results are compared to those from other theoretical and experimental works.

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

The three-dimensional turbulent cascade flows with and without endwall fences are numerically investigated by solving the incompressible Navier-Stokes equations with a high-Reynolds number $k-{\varepsilon}$ turbulence closure model. A projection method based algorithm is used in the finite-volume formulation, with the second order upwind-differencing scheme for the convective terms. First, assessments on accuracy of the present method are made by comparing the static pressure distributions at the mid-span of the cascade with measured data, and also by confirming the experimental observations on the choice of an optimal fence height for the secondary flow control. In understanding the three-dimensional nature of the secondary flow in turbine cascade, the limiting streamline patterns and the static pressure contours at the suction surface of the blade as well as on the cascade endwall are employed to visualize the effectiveness of the endwall fence for the secondary flow control. Analysis on the streamwise vorticity contour maps along the cascade with the three-dimensional representation of their iso-surfaces reveals the strucuture of the complicated vortical flow in the turbine cascade with endwall fence, and also leads to an understanding on formation of the counter-rotating streamwise vortex over the endwall fence, in explaining the mechanisms of controlling the secondary flow and also for the proper selection of an optimal fence height.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제12권1호
• /
• pp.918-927
• /
• 2020

Flows induced by hybrid CRP pod propulsion systems (CRP-POD) are fundamentally characterized by unsteadiness. This work presents a numerical study on the unsteady flow of a CRP-POD at behind-hull condition based on CFD (Computational Fluid Dynamics). Unsteady RANS method is adopted, coupled with SST k-u turbulence model and sliding mesh method. The propeller thrusts and torques obtained by CFD is validated by model tests and acceptable agreements are obtained. The time histories of shingle-blade loads and pressures near the hull surface are recorded for the analysis of unsteady flow features. The cases of forward propeller alone and aft propeller alone are also computed to distinguish the hull-propeller interaction and propeller-propeller interaction. The results show the blade loads of both forward and aft propellers strongly fluctuate with phase angles. For the forward propeller, the blade load fluctuation is mainly governed by the hull-propeller interaction, while the aft blade load is remarkably affected by the propeller-propeller interaction in terms of the load average and fluctuation pattern. The fields of pressure, vorticity and velocity are also analyzed to reveal the unsteady flow features.

• 한국해안해양공학회지
• /
• 제17권4호
• /
• pp.294-306
• /
• 2005

지금까지 침매터널에 대한 대부분의 연구들은 건설기술자의 경험에 바탕을 둔 연구들이었다. 공학적인 관점에서 파랑하중에 의한 침매터널의 안정성은 지반기초의 거동에 큰 영향을 받으므로 파${\cdot}$지반${\cdot}$침매터널의 상호작용에 대한 이해는 매우 중요하다. 본 연구에서는 파${\cdot}$지반${\cdot}$침매터널의 상호작용의 메커니즘에 대한 이해를 위하여 부산-거제간 연결도로 민간투자사업 구간 중 침매터널 구간에 직접수치해석기법을 적용하여 침매터널 주변의 지반 및 피복층 내에서 파랑하중에 의한 간극수압, 와도 및 흐름의 특성을 종합적으로 살펴보았다. 또한 다양한 파랑조건 및 피복층의 평균입경변화에 따른 비선형성의 파랑과 침매터널 및 지반의 상호작용이 침매터널의 안정성과 관련하여 논의되었다.

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

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.690-695
• /
• 2001

The three-dimensional spatial structures of impeller flow created by a six bladed Rushton turbine have identified based on the volumetric velocity information from multi-plane stereoscopic PIV measurements. A total of 10 planes with 2 mm space with a 50 mm by 64 mm size of the field of view were targeted. To reduce the depth of focus, we adopted an angle offset configuration which satisfied the Scheimpflug condition. The distortion compensation procedure was utilized during the in situ calibration. Phase-locked instantaneous data were ensemble averaged and interpolated in order to obtain mean 3-D, volumetric velocity fields on a 60 degree sector of a cylindrical ring volume enclosing the turbine blade. Using the equi-vorticity surface rendering, the spatial structure of the trailing vortices was clearly demonstrated. Detail flow characteristics of the radial jet reported in previous studies of mixer flows were easily identified.

• Wind and Structures
• /
• 제14권2호
• /
• pp.167-186
• /
• 2011

This paper presents the results of measurements relating to the aerodynamic forces on flat square plates which were allowed to rotate at different speeds about their horizontal axis, by modifying the velocity of the incoming flow. A 1 m square test-sheet and a 0.3 m square test-sheet were fitted with a number of pressure sensors in order to obtain information relating to the instantaneous pressure distribution acting on the test-sheet; a compact gyroscope to record the angular velocity during the rotational motion was also implemented. Previous work on autorotation has illustrated that the angular velocity varies with respect to the torque induced by the wind, the thickness and aspect ratio of the test-sheet, any frictional effects present at the bearings, and the vorticity generated through the interaction between the plate and the wind flow. The current paper sets out a method based on the solution of the equation of motion of a rotating plate which enables the determination of angular velocities on autorotating elements to be predicted. This approach is then used in conjunction with the experimental data in order to evaluate the damping introduced by the frictional effects at the bearings during steady autorotation.

• 대한기계학회논문집B
• /
• 제25권7호
• /
• pp.915-922
• /
• 2001

A new phase-average technique using instantaneous velocity fields obtained by a PIV method has been developed. The technique tracks vorticity centers and estimates the value of circulation for a chosen domain. The locations of vortex-centers and the magnitudes of circulation are matched together then showing a sine wave feature due to the periodic vortex shedding from the circular cylinder. Ensemble averaged and phase averaged velocity fields are successfully measured for the circular cylinder wake where Reynolds number is 3900 based on free stream velocity and cylinder diameter. The convection velocities of the vortices center and the vortex shedding frequency were measured by a single hot-wire probe.

• Journal of Advanced Marine Engineering and Technology
• /
• 제34권8호
• /
• pp.1094-1099
• /
• 2010

수중에 있는 구 주위의 3차원 층류유동을 수치 해석하였다. 유한체적법을 기반으로 나비에 스톡스 방정식을 비구조격자로 풀어 물체주위의 동역학 현상을 규명하였다. 보티시티, 속도, 동압, 레지듀얼, 항력계수 등의 데이터를 가지고 분석하였다. 레이놀드 수는 425, 300, 250, 100 이다. 구 사이의 거리가 작은 경우가 구 한 개의 후류 유동보다 안정됨을 보였고 구 사이의 간격이 커짐에 따라 유동현상은 불안정하게 나타났다. 후방의 구는 정체압력의 감소를 보였고 이로 인해 저항이 줄어들게 됨을 알 수 있었다. 항력계수의 비교는 본 수치계산이 타당함을 입증하였다.