• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.262-270
• /
• 1987

A numerical method is presented which can solve the unsteady momentum and thermal boundary layers, coupled through the agency of buoyancy force, over a heated circular cylinder impulsively started from rest. By linearizing the nonlinear finite difference equations without sacrificing accuracy, numerical solutions are obtained at each time step without iteration. To get rid of the requirement of excessive number of grid points in the region of reversed flow, special form of transformed variables are used, by which the computational boundary layer thickness is maintained almost constant. These numerical properties enable the method to easily handle the region of reversed flow and how the singularity develops in the interior of the boundary layer. In order to investigated the thermal effects on the skin friction, heat flux, displacement thickness and on the separation, we have successfully solved three different cases of the buoyancy parameter .alpha.(Gr/Re$^{2}$).

• Journal of the Korean Society of Propulsion Engineers
• /
• v.2 no.1
• /
• pp.1-12
• /
• 1998

Industrial ejector system is a facility to transport, to compress or to pump out a low pressure secondary flow by using a high pressure primary flow. An advantage of the ejector system is in its geometrical simplicity, not having any moving part, compared with other fluid machinery. Most of the previous works have been performed experimentally and analytically. The obtained data. are too insufficient to improve our current understanding on the detailed flow field inside the ejector. In order to provide more comprehensive data on this ejector flow field, two-dimensional computations using Reynolds-averaged Navier-Stokes equations were performed for a very wide range of operating pressure ratio of the supersonic ejector with a secondary throat. The current results showed that the supersonic ejector system has an optimum pressure ratio for the secondary flow total pressure to be minimized. The numerical results clearly revealed the shock system, shock/boundary layer interaction, and secondary flow entrainment inside the supersonic ejector.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.11a
• /
• pp.299-304
• /
• 2008

Starting characteristics of the axi-symmetric second throat exhaust diffuser (STED) with zero-secondary flows are numerically investigated. Renolds-Average Navier-Stokes equations with a standard ${\kappa}-{\varepsilon}$ turbulence model incorporated with enhanced wall treatment are solved to simulate the diffusing evolutions of the nozzle plume. Minimum (optimum) starting pressure difference of 20$\sim$25% between 1-D theory and the measured data validated from previous results[5] is also applied to predict the range of an effective diffuser expansion ratio (Ad/At) in this system.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.4 no.2
• /
• pp.151-161
• /
• 2012

The present study numerically investigates the aerodynamic characteristics of two-dimensional wings in the vicinity of the ground by solving two-dimensional steady incompressible Navier-Stokes equations with the turbulence closure model of the realizable k-${\varepsilon}$ model. Numerical simulations are performed at a wide range of the normalized ground clearance by the chord length ($0.1{\leq}h/C{\leq}1.25$) for the angles of attack ($0^{\circ}{\leq}{\alpha}{\leq}10^{\circ}$) in the prestall regime at a Reynolds number (Re) of $2{\times}10^6$ based on free stream velocity $U_{\infty}$ and the chord length. As the physical model of this study, a cambered airfoil of NACA 4406 has been selected by a performance test for various airfoils. The maximum lift-to-drag ratio is achieved at ${\alpha}=4^{\circ}$ and h / C = 0.1. Under the conditions of ${\alpha}=4^{\circ}$ and h / C = 0.1, the effect of the Reynolds number on the aerodynamic characteristics of NACA 4406 is investigated in the range of $2{\times}10^5{\leq}Re{\leq}2{\times}10^9$. As Re increases, $C_l$ and $C_d$ augments and decreases, respectively, and the lift-to-drag ratio increases linearly.

• Journal of the Korean Solar Energy Society
• /
• v.35 no.1
• /
• pp.35-43
• /
• 2015

A numerical simulation on the wake flow of a wind turbine which is a scaled version of a multi-megawatt wind turbine has been performed. Two different inlet conditions of averaged wind speed including one below and one above the rated wind speed were used in the simulation. Steady-state pitch angles of the blade associated with the two averaged wind speeds were imposed for the simulation. The steady state analysis based on the Reynolds averaged Navier-Stokes equations with the method of frame motion were used for the simulation to find the torque of the rotor and the wake field behind the wind turbine. The simulation results were compared with the results obtained from the wind tunnel testing. From comparisons, it was found that the simulation results on the turbine power are pretty close to the experimental values. Also, the wake results were relatively close to the experimental results but there existed some discrepancy in the shape of velocity deficit. The reason for the discrepancy is considered due to the steady state solution with the frame motion method used in the simulation. However, the method is considered useful for solutions with much reduced calculation time and reasonably good accuracy compared to the transient analysis.

• Computers and Concrete
• /
• v.21 no.6
• /
• pp.717-726
• /
• 2018

Concrete pipelines are the most efficient and safe means for gas and oil transportation over a long distance. The use of nano materials and nono-engineering can be considered for enhancing concrete pipelines properties. the tests show that $SiO_2$ nanoparticles can improve the mechanical behavior of concrete. Moreover, severe hazard for pipelines is seismic ground motion. Over the years, scientists have attempted to understand pipe behavior against earthquake most frequently via numerical modeling and simulation. Therefore, in this paper, the dynamic response of underwater nanocomposite submerged pipeline conveying fluid is studied. The structure is subjected to the dynamic loads caused by earthquake and the governing equations of the system are derived using mathematical model via Classic shell theory and Hamilton's principle. Navier-Stokes equation is employed to calculate the force due to the fluid in the pipe. As well, the effect of external fluid is modeled with an external force. Mori-Tanaka approach is used to estimate the equivalent material properties of the nanocomposite. 1978 Tabas earthquake in Iran is considered for modelling seismic load. The dynamic displacement of the structure is extracted using differential quadrature method (DQM) and Newmark method. The effects of different parameters such as $SiO_2$ nanoparticles volume percent, boundary conditions, thickness to radius ratios, length to radius ratios, internal and external fluid pressure and earthquake intensity are discussed on the seismic response of the structure. From results obtained in this paper, it can be found that the dynamic response of the pipe is increased in the presence of internal and external fluid. Furthermore, the use of $SiO_2$ nanoparticles in concrete pipeline reduces the displacement of the structure during an earthquake.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.6
• /
• pp.474-481
• /
• 2018

This paper is to compare by numerical analysis the flow characteristics and propulsion performance of stern with the shape change of K-duct, a pre-swirl duct developed by Korea Research Institute of Ships & Ocean Engineering (KRISO). First, the characteristics of the propeller and the resistance and self-propulsion before and after the attachment of the K-duct to the ship were verified and the validity of the calculation method was confirmed by comparing this result with the model test results. After that, resistance and self-propulsion calculations were performed by the same numerical method when the K-duct was changed into five different shapes. The efficiency of the other five cases was compared using the delivery horsepower in the model scale and the flow characteristics of the stern were analyzed as the velocity and pressure distributions in the area between the duct end and the propeller plane. For the computation, STAR-CCM +, a general-purpose flow analysis program, was used and the Reynolds Averaged Navier-Stokes (RANS) equations were applied. Rigid Body Motion (RBM) method was used for the propeller rotating motion and SST $k-{\omega}$ turbulence model was applied for the turbulence model. As a result, the tangential velocity of the propeller inflow changed according to the position angle change of the stator, and the pressure of the propeller hub and the cap changes. This regulated the propeller hub vortex. It was confirmed that the vortex of the portion where the fixed blade and the duct meet was reduced by blunt change.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.1
• /
• pp.76-84
• /
• 2019

For decades, solitary waves have commonly been used to simulate tsunami conditions in numerical studies. However, the main component of a tsunami waveform acts at completely different spatial and temporal distributions than a solitary waveform. Thus, this study applied a 2-D numerical wave tank that included a non-reflected tsunami generation system based on Navier-Stokes equations (LES-WASS-2D) to directly simulate the run-up of a tsunami-like solitary wave on a slope. First, the waveform and velocity due to the virtual depth factor were applied to the numerical wave tank to generate a tsunami, which made it possible to generate the wide waveform of a tsunami, which was not reproduced with the existing solitary wave approximation theory. Then, to validate the applied numerical model, the validity and effectiveness of the numerical wave tank were verified by comparing the results with the results of a laboratory experiment on a tsunami run-up on a smooth impermeable 1:19.85 slope. Using the numerical results, the run-up characteristics due to a tsunami-like solitary wave on an impermeable slope were also discussed in relation to the volume ratio. The maximum run-up heights increased with the ratio of the tsunami waveform. Therefore, the tsunami run-up is highly likely to be underestimated compared to a real tsunami if the solitary wave of the approximation theory is applied in a tsunami simulation in a coastal region.

• Journal of Convergence for Information Technology
• /
• v.9 no.8
• /
• pp.155-163
• /
• 2019

Laminar mixed convection of a nanofluid consists of water and $Al_2O_3$ in a horizontal circular tube has been studied numerically. Two-phase mixture model has been used to investigate hydrodynamic and thermal behaviors of the nanofluid with variables physical properties. Three dimensional Navier-Stokes, energy and volume fraction equations have been discretized using the finite volume method. The Brownian motions of nanoparticles have been considered to determine the thermal conductivity and dynamic viscosity of $Al_2O_3$-Water nanofluid, which depend on temperature. The calculated results show good agreement with the previous numerical data. Results show that in a given Reynolds number (Re), increasing solid nanoparticles volume fraction and Richardson number (Ri) increases the convective heat transfer coefficient and wall shear stress.

• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.4
• /
• pp.361-372
• /
• 2019

The flow pattern of air layers and skin-friction drag reduction by air injection are investigated to find the suitable multiphase flow model using unstructured finite-volume CFD solver for the Reynolds-averaged Navier-Stokes equations. In the present computations, two different multiphase flow modeling approaches, such as the Volume of Fluid (VOF) and the Eulerian Multi-Phase (EMP), are adopted to investigate their performances in resolving the two-phase flow pattern and in estimating the frictional drag reduction. First of all, the formation pattern of air layers generated by air injection through a circular opening on the bottom of a flat plate are investigated. These results are then compared with those of MMkiharju's experimental results. Subsequently, the quantitative ratios of skin-friction drag reduction including the behavior of air layers, within turbulent boundary layers in large scale and at high Reynolds number conditions, are investigated under the same conditions as the model test that has been conducted in the US Navy's William B. Morgan Large Cavitation Channel (LCC). From these results, it is found that both VOF and EMP models have similar capability and accuracy in capturing the topology of ventilated air cavities so called'air pockets and branches'. However, EMP model is more favorable in predicting quantitatively the percentage of frictional drag reduction by air injection.