• Journal of the Korean Society of Safety
• /
• v.10 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.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.05a
• /
• 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
• /
• v.12 no.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.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.4
• /
• pp.294-306
• /
• 2005

Most immersed tunnels investigated have been investigated based on the engineer's experience with design and construction. From engineering point of view, it is very important to understand the wave interaction with the seabed and immersed tunnel, since the stability of an immersed tunnel depends largely on the behavior of the seabed foundation. In this study, for the first stage research to find out the mechanism of the wave interaction with the seabed and immersed tunnel, the benchmarking method called as direct numerical simulation (DNS) was employed to analyze comprehensively the wave-induced pore water pressures, vorticity and flows in seabed or inside rubble stone around the immersed tunnel. The immersed tunnel is modeled based on Busan-Geoje fixed link project in Korea, which is now on the stage of planning. Moreover, the nonlinear water wave interaction with an immersed tunnel/its seabed foundation was thoroughly examined with regard to the stabilities of the immersed tunnel subjected to various water wave conditions, median grain size and so forth.

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

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• 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
• /
• v.14 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.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
• /
• v.34 no.8
• /
• pp.1094-1099
• /
• 2010

Three dimensional laminar flows are numerically simulated around the submerged spheres. The finite volume based Navier-Stokes equations with unstructured grids are solved to make clear the hydrodynamic phenomena near and far away from the body. Reviews are made on with the vorticity, velocity, dynamic pressure, residuals, drags, etc. The Reynolds numbers under consideration are 425, 300, 250 and 100. In case of small spacing between spheres, the flow is more stable than that past a single sphere. According to the analysis, the flow past two spheres is found to be unstable as the spacing becomes larger. The rear sphere shows the deduction of stagnation pressure, which results in the decrease of the resistance. The predicted drag coefficients verify that the present numerical calculation is reasonable.