• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• 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.13 no.2
• /
• pp.42-53
• /
• 2009

Gas turbine engine simulation program has been developed. In compressor and turbine, 2-D NS implicit code is used with k-$\omega$ SST turbulent model. In combustor, 0-D lumped method chemical equilibrium code is adopted under the limitations, the products are only 10 species of molecular and air-fuel is perfectly mixed state with 100% combustion efficiency at constant pressure. Fluid properties are shared on interfaces between engine components. The outlet conditions of compressor have been used as the inlet condition of combustor. The inlet condition of turbine comes from the compressor The back pressure in compressor outlet is transferred by the inlet pressure of turbine. Unsteady phenomena at rotor-stator in compressor and turbine is covered by mixing-plane method. The state of engine can be determined only by given inlet condition of compressor, outlet condition of turbine, equivalence ratio and rotating speed.

• 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.

• Computers and Concrete
• /
• v.24 no.2
• /
• pp.125-135
• /
• 2019

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 SiO2 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 SiO2 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 SiO2 nanoparticles in concrete pipeline reduces the displacement of the structure during an earthquake.

• Journal of Ocean Engineering and Technology
• /
• v.36 no.3
• /
• pp.143-152
• /
• 2022

To reach a port, a ship must pass through a shallow water zone where seabed effects alter the hydrodynamics acting on the ship. This study examined the maneuvering characteristics of an autonomous surface ship at 3-DOF (Degree of freedom) motion in deep water and shallow water based on the in-port speed of 1.54 m/s. The CFD (Computational fluid dynamics) method was used as a specialized tool in naval hydrodynamics based on the RANS (Reynolds-averaged Navier-Stoke) solver for maneuvering prediction. A virtual captive model test in CFD with various constrained motions, such as static drift, circular motion, and combined circular motion with drift, was performed to determine the hydrodynamic forces and moments of the ship. In addition, a model test was performed in a square tank for a static drift test in deep water to verify the accuracy of the CFD method by comparing the hydrodynamic forces and moments. The results showed changes in hydrodynamic forces and moments in deep and shallow water, with the latter increasing dramatically in very shallow water. The velocity fields demonstrated an increasing change in velocity as water became shallower. The least-squares method was applied to obtain the hydrodynamic coefficients by distinguishing a linear and non-linear model of the hydrodynamic force models. The course stability, maneuverability, and collision avoidance ability were evaluated from the estimated hydrodynamic coefficients. The hydrodynamic characteristics showed that the course stability improved in extremely shallow water. The maneuverability was satisfied with IMO (2002) except for extremely shallow water, and collision avoidance ability was a good performance in deep and shallow water.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3003-3008
• /
• 2007

Critical nozzle has been frequently employed to measure the flow rate of various gases, but hydrogen gas, especially being at high-pressure condition, was not nearly dealt with the critical nozzle due to treatment danger. According to a few experimental data obtained recently, it was reported that the discharge coefficient of hydrogen gas through the critical nozzle exceeds unity in a specific range of Reynolds number. No detailed explanation on such an unreasonable value was made, but it was vaguely inferred as real gas effects. For the purpose of practical use of high-pressure hydrogen gas, systematic research is required to clarify the critical nozzle flow of high-pressure hydrogen gas. In the present study, a computational fluid dynamics(CFD) method has been applied to predict the critical nozzle flow of high-pressure hydrogen gas. Redlich-Kwong equation of state that take account for the forces and volume of molecules of hydrogen gas were incorporated into the axisymmetric, compressible Navier-Stokes equations. A fully implicit finite volume scheme was used to numerically solve the governing equations. The computational results were validated with some experimental data available. The results show that the coefficient of discharge coefficient is mainly influenced by the compressibility factor and the specific heat ratio, which appear more remarkable as the inlet total pressure of hydrogen gas increases.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.14-19
• /
• 2008

Effects of the Reynolds and Knudsen numbers on a micro-viscous pump are studied by using a Navier-Stokes code based on a finite volume method. The micro viscous pump consists of a circular rotor and a two-dimensional channel. The channel walls are treated by using a slip velocity model. The Reynolds number is studied in the range of $0.1{\sim}50$. The Knudsen number varies from 0.01 to 0.1. Numerical solutions show that the pump works efficiently when two counter rotating vortices formed on both sides of the rotor have the same size and intensity. As the Reynolds number increases, the size and intensity of the vortex on the inlet side of the pump decrease. It disappears when the Reynolds number is larger than about Re=20. The characteristics of the performance of the pump is shown to deteriorate, in terms of mean velocity and pressure rise, as the Reynolds number increases. The Knudsen number shows a different effect on the characteristics of the pump. As it increases, the mean velocity and pressure rise decrease but the characteristics of the vortex flow remains unchanged, unlike the effect of Reynolds number.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.8
• /
• pp.764-772
• /
• 2009

In this study, fluid/structure coupled analyses have been conducted for 3-D stator and rotor configuration. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate fluid/structure responses of general stator-rotor configurations. To solve the fluid/structure coupled problems, fluid domains are modeled using the structural grid system with dynamic moving and local deforming techniques. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras(S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3-D turbine blades for fluid-structure interaction(FSI) problems. Detailed fluid/structure analysis responses for stator-rotor interaction flow conditions are presented to show the physical performance and flow characteristics.