• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.2
• /
• pp.142-151
• /
• 1991

In the present investigation, it has been proposed to utilize a direct contact heat exchanger as an evaporator to solve the difficulties such as scaling, corrosion and law thermal efficiencies, associated with the conventional evaporator. Liquified nitrozen was utilized as a working fluid to investigate basic natures of bubble dynamics in the evaporator, and spray nozzles were adopted to inject liquified nitrozen into the spray column with varying flow rates of dispersed phase fluids. Experimentations were carried out in the range of $6.54{\times}10^{-4}kg/s$ - 0.030 kg/s for dispersed phase flow rates with one, three and five nozzle holes. Observing the bubble dynamics for the evaporator the feasibility of utilizing a direct contact heat exchanger as a LNG evaporator has been evaluated. The results show that no eruption phenomena was observed in the present investigation with $LN_2$ and the interface between $N_2$ bubbles and water was fully turbulent. It is believed that the high injection velocity of $LN_2$ through the spray nozzles provide good mixing effects for both heat and mass transfers between water and $N_2$ bubbles. Ice was formed on the surface of the spray nozzle for higher $LN_2$ flow rates. However, even in this case, it is observed that the ice was detached as soon as it was formed. Under the present experimental conditions, the shapes of $LN_2$ bubbles were in the spherical-cap region according to the Clift, Grace and Weber Graphs. The height of foam region caused by the breakup of larger bubbles keeps increasing with high injection velocities until it reaches it's maximum height.

• Journal of Advanced Marine Engineering and Technology
• /
• v.34 no.1
• /
• pp.76-83
• /
• 2010

The Characteristics of the flow field through cylinders with in-line and staggered arrangements were calculated by vortex method. Vortex distributions and velocity profiles around the cylinders with in-line and staggered arrangements were simulated at the pitch ratio of Pt/D=1.25~2.0 and Reynolds number of Re=$4.0{\times}10^1{\sim}4.0{\times}10^4$. As the results the vortices of clockwise at the upper separation point cylinder and the vortices of anticlockwise at the lower separation point of each cylinder were generated at both in-line and staggered arrangements. The generation of the reverse flow in the rear region of the cylinders was caused by the pitch ratio and Reynolds number, the boundary region was at the pitch ratio of Pt/D=1.5 and Reynolds number of Re=$4.0{\times}10^2{\sim}4.0{\times}10^3$ in case of in-line arrangement and was at the pitch ratio of Pt/D=1.4 and Reynolds number of Re=$4.0{\times}10^1{\sim}4.0{\times}10^2$ in case of staggered arrangement.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.1
• /
• pp.83-93
• /
• 1995

A numerical and experimental investigation on the flow characteristics in the rectangular duct of an MHD propulsion system has been carried out. In numerical analysis, three-dimensional, steady-state, viscous, incompressible electrically conducting fluid flow under the influence of uniformly applied magnetic and electric fields was treated using a finite-difference technique. It was found from the numerical study that when the Lorentz force is weak, the typical parabolic velocity profile under a laminar flow condition changes to an M shaped profile near the electrode region and that the pressure increases linearly from the inlet toward the outlet of the MHD duct under constant electro-magnetic field. In experiment, thrust of the MHD propulsion system can be controlled easily by varying electrode current. The measured pressure gradient along the MHD duct is proportional to the Lorentz force, which is in agreement with the numerical results.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.93-97
• /
• 2005

Air interlacing serves to protect the yarn against damage, strengthens inter-filament compactness or cohesion, and ensures fabric consistency. The air interlacing nozzle is used to introduce intermittent nips to a filament yarn so as to improve its performance in textile processing. This study investigates the effect of interlacing nozzle geometry on the interlacing process. The geometries of interlacing nozzles with multiple air inlets located across the width of a yarn channels are investigated. The basic interlacing nozzle is the yarn channel, with a perpendicular single air inlet in the middle. The yarn channel shapes are cross sections with semicircular or rectangular shapes. This paper presents three doubled sub air inlets with main air inlet and one of them is slightly inclined doubled sub air inlets with main air inlet. The compressed air coming out from the inlet hits the opposing wall of the yarn channel, divides into two branches, flows trough the top side of yarn channel, joins with the compressed air coming out from the sub air inlet and then creates two free jets at both ends of the yarn channel. The compressed air moves in the shape of two opposing directional vortices. The CFD-FASTRAN was used to perform steady simulations of impinging jet flow inside of the interlace nozzles. The vortical structure and the flow pattern such as pressure contour, particle traces, velocity vector plots inside of interlace nozzle geometry are discussed in this paper.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.7
• /
• pp.683-693
• /
• 2013

In this study, the mixing characteristics in lid-driven cavity flows were studied numerically by using a hybrid lattice Boltzmann method (HLBM). First, we compared the numerical results from single-relaxation-time (LB-SRT) and multi-relaxation-time (LB-MRT) models to examine their reliability. In most of the cavity flow, the results from both the LB-SRT and the LB-MRT models were in good agreement with those using a Navier-Stokes solver for Re=100-5000. However, the LB-MRT model was superior to the LB-SRT model for the simulation of higher Reynolds number flows having a geometrical singularity with much lesser spatial oscillations. For this reason, the LB-MRT model was selected to study the mass transport in lid-driven cavity flows, and it was demonstrated that mass transport in the fluid was activated by a recirculation zone in the cavity, which is connected from the top to the bottom surfaces through two boundary layers. Various mixing characteristics such as the concentration profiles, mean Sherwood (Sh) numbers, and velocity were computed. Finally, the detailed transport mechanism and solutions for the concentration profile in the cavity were presented.

• Journal of Environmental Impact Assessment
• /
• v.20 no.4
• /
• pp.539-555
• /
• 2011

Newly constructed, high-rise dense building areas by urban development can cause changes in local wind fields. Wind fields were analyzed to assess the impact on the local meteorology due to the land use changes during the urban redevelopment called "Eunpyeong new town" in north-western Seoul using CFD_NIMR_SNU (Computational Fluid Dynamics, National Institute of Meteorological Research, Seoul National University) model. Initial value of wind speed and direction use analysis value of AWS (Automatic Weather Station) data during 5 years. In the case of the pre-construction with low rise built-up area, it was simulated that the spatial distribution of horizontal wind fields depends on the topography and wind direction of initial inflow. But, in the case of the post-construction with high rise built-up area, it was analyzed that the wind field was affected by high rise buildings as well as terrain. High-rise buildings can generate new circulations among buildings. In addition, small size vortexes were newly generated by terrain and high rise buildings after the construction. As high-rise buildings act as a barrier, we found that the horizontal wind flow was separated and wind speed was reduced behind the buildings. CFD_NIMR_SNU was able to analyze the impact of high-rise buildings during the urban development. With the support of high power computing, it will be more common to utilize sophisticated numerical analysis models such as CFD_NIMR_SNU in evaluating the impact of urban development on wind flow or channel.

• 유체기계공업학회:학술대회논문집
• /
• 2004.12a
• /
• pp.396-401
• /
• 2004

In this paper, a 1MW HAWT(FIL-1000) rotor blade has been designed by BEMT(Blade Element Momentum Theory) with Prandtl's tip loss. Also, a 3-D flow and performance analysis on the FIL-1000 rotor blade has been carried out by using the 3-D Navier-Stokes commercial solver (CFX-5.7) to provide more efficient design techniques to the large-scale HAWT engineers. The rated power and itsapproaching wind velocity at design point (TSR=7.5) are 1MW and 9.99m/s respectively. The rotor diameter is 54.5m and the rotating speed is 26.28rpm. Airfoils such as FFA W-301, DU91-W-250, DU93-W-210, NACA 63418, NACA 63415 consist of the rotor blade from hub to tip. Recent CFX version, 5.7 was adopted to simulate 3-D flow field and to analyze the performance characteristics of the rotor blade. Entire mesh node number is about 730,000 and it is generated by ICEM-CFD to achieve better mesh quality The predicted maximum power occurringat the design tip speed ratio is 931.45kW. Approaching to the root, the inflow angle becomes large, which causesthe blade to be stalled in the region. Therefore, k-$\omega$ SST turbulence model was used to predict the quantitative flow information more accurately. Application of commercial CFD code to optimum blade design and performance analysis was proved to be more effective environment to HAWT blade designers.

• Wind and Structures
• /
• v.13 no.6
• /
• pp.557-580
• /
• 2010

A new inflow turbulence generation method and a combined dynamic SGS model recently developed by the authors were applied to evaluate the wind effects on 508 m high Taipei 101 Tower. Unlike the majority of the past studies on large eddy simulation (LES) of wind effects on tall buildings, the present numerical simulations were conducted for the full-scale tall building with Reynolds number greater than $10^8$. The inflow turbulent flow field was generated based on the new method called discretizing and synthesizing of random flow generation technique (DSRFG) with a prominent feature that the generated wind velocity fluctuations satisfy any target spectrum and target profiles of turbulence intensity and turbulence integral length scale. The new dynamic SGS model takes both advantages of one-equation SGS model and a dynamic production term without test-filtering operation, which is particular suitable to relative coarse grid situations and high Reynolds number flows. The results of comparative investigations with and without generation of inflow turbulence show that: (1) proper simulation of an inflow turbulent field is essential in accurate evaluation of dynamic wind loads on a tall building and the prescribed inflow turbulence characteristics can be adequately imposed on the inflow boundary by the DSRFG method; (2) the DSRFG can generate a large number of random vortex-like patterns in oncoming flow, leading to good agreements of both mean and dynamic forces with wind tunnel test results; (3) The dynamic mechanism of the adopted SGS model behaves adequately in the present LES and its integration with the DSRFG technique can provide satisfactory predictions of the wind effects on the super-tall building.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.11
• /
• pp.1169-1176
• /
• 2011

Liquid film flows are classified into waveless laminar, wavy laminar, and turbulent flows depending on the Reynolds number or the flow stability. Since the wavy motions of the film flows are so intricate and nonlinear, studies on them have largely been experimental. Most numerical approaches have been limited to the waveless flow regime. The various free surface-tracking schemes adopted for this problem were used to more accurately estimate the average film thickness, rather than to capture the unsteady wavy motion. In this study, the wavy motions in laminar wavy liquid film flows with Reynolds numbers of 200-1000 were simulated with various numerical schemes based on the volume of fluid (VOF) method for interface tracking. The results from each numerical scheme were compared with the experimental results in terms of the average film thickness, the wave velocity, and the wave amplitude.

• The Journal of the Acoustical Society of Korea
• /
• v.37 no.2
• /
• pp.92-98
• /
• 2018

This paper is a study on the generation mechanism of propeller singing based on the cavitation tunnel test, underwater impact test, finite element analysis and computational flow analysis for the model propeller. A wire screen mesh, a propeller and a rudder were installed to simulate ship stern flow, and occurrence and disappearance of propeller singing phenomenon were measured by hydrophone and accelerometer. The natural frequencies of propeller blades were predicted through finite element analysis and verified by contact and non-contact impact tests. The flow velocity and effective angle of attack for each section of the propeller blades were calculated using RANS (Reynolds Averaged Navier-Stokes) equation-based computational fluid analysis. Using the high resolution analysis based on detached eddy simulation, the vortex shedding frequency calculation was performed. The numerical predicted vortex shedding frequency was confirmed to be consistent with the singing frequency and blade natural frequency measured by the model test.