• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.2
• /
• pp.177-184
• /
• 2008

The dynamic stability of elastically restrained cantilever pipe conveying fluid with crack is investigated in this paper. The pipe, which is fixed at one end, is assumed to rest on an intermediate spring support. Based on the Euler-Bernoulli beam theory, the equation of motion is derived by the energy expressions using extended Hamilton's Principle. The crack section is represented by a local flexibility matrix connecting two undamaged pipe segments. The influence of a crack severity and position, mass ratio and the velocity of fluid flow on the stability of a cantilever pipe by the numerical method are studied. Also, the critical flow velocity for the flutter and divergence due to variation in the support location and the stiffness of the spring support is presented. The stability maps of the pipe system are obtained as a function of mass ratios and effect of crack.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.186-191
• /
• 2004

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.397-400
• /
• 2010

Unstructured adaptive grid flow simulation is applied to the calculation of high speed compressible flow of inert and reactive gas mixtures. Computational results are presented for the case of premixed hydrogen-air supersonic flow over a 2-D wedge. In such a configuration, combustion may be triggered behind the oblique shock wave and transition to an oblique detonation wave is eventually obtained. It is shown that the solution adaptive procedure implemented is able to correctly define the important wave front.

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.4
• /
• pp.352-359
• /
• 2010

In recent years, the market has shown increasing interest in pump-turbines. The prompt availability of pumped storage plants and the benefits to the power system achieved by peak lopping, providing reserve capacity, and rapid response in frequency control are providing a growing advantage. In this context, there is a need to develop pumpturbines that can reliably withstand dynamic operation modes, fast changes of discharge rate by adjusting the variable diffuser vanes, as well as fast changes from pumping to turbine operation. In the first part of the present study, various flow patterns linked to operation of a pump-turbine system are discussed. In this context, pump and turbine modes are presented separately and different load cases are shown in each operating mode. In order to create modern, competitive pump-turbine designs, this study further explains what design challenges should be considered in defining the geometry of a pump-turbine impeller. The second part of the paper describes an innovative, staggered approach to impeller development, applied to a low head pump-turbine project. The first level of the process consists of optimization strategies based on evolutionary algorithms together with 3D in-viscid flow analysis. In the next stage, the hydraulic behavior of both pump mode and turbine mode is evaluated by solving the full 3D Navier-Stokes equations in combination with a robust turbulence model. Finally, the progress in hydraulic design is demonstrated by model test results that show a significant improvement in hydraulic performance compared to an existing reference design.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.1
• /
• pp.101-109
• /
• 2008

The stability of a rotating cantilever pipe conveying fluid with a crack and tip mass is investigated by the numerical method. That is, the effects of the rotating angular velocity, mass ratio, crack severity and tip mass on the critical flow velocity for flutter instability of system are studied. The equations of motion of rotating pipe are derived by using the Euler-Bernoulli beam theory and the extended Hamilton's principle. The crack section of pipe is represented by a local flexibility matrix connecting two undamaged pipe segments. Also, the crack is assumed to be in the first mode of fracture and always opened during the vibrations. When the tip mass and crack are constant, the critical flow velocity for flutter is proportional to the rotating angular velocity of pipe. In addition, the stability maps of the rotating pipe system as a rotating angular velocity and mass ratio ${\beta}$ are presented.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.723-728
• /
• 2008

We performed numerical analysis of base drag phenomenon, when a projectile with backward step flies into atmosphere at supersonic speed. We compared with other researchers. From our previous studies that were 2-dimensional simulation, we found out from sophisticated simulations that need dense mesh points to compare base pressure and velocity profile after from base with experimental data. Therefore, we focus on high order spatial disceretization over 3rd order with TVD such as MUSCL TVD 3rd, 5th, and WENO 5th order, and Limiters such as minmod, Triad. Moreover, we enforce to flux averaging schemes such as Roe, RoeM, HLLE, AUSMDV. In present, one dimensional result of Euler tests, there are Sod, Lax, Shu-Osher and interacting blast wave problems. AUSMDV as a flux averaging scheme with MUSCL TVD 5th order as spatial resolution is good agreement with exact solutions than other combinations. We are carrying out the same approaches into 3-dimensional base flow only candidate flux schemes that are Roe, AUSMDV. Additionally, turbulence models are used in 3-dimensional flow, one is Menter s SST DES model and another is Sparlat-Allmaras DES/DDES model in Navier-Stokes equations.

• Proceedings of the SAREK Conference
• /
• 2008.11a
• /
• pp.546-550
• /
• 2008

In this study, 3-dimensional Computational Fluid Dynamics (CFD) analysis was induced to simulate air flow and particle motion in the axial flow cyclone separator. The commercialized CFD code FLUENT was used to visualize pressure drop and particle collection efficiency inside the cyclone. We simulated 4 cyclone models with different shape of vane, such as turning angle or shape of cross section. For the air flow simulation, we calculated the flow field using standard ${\kappa}-{\varepsilon}$ turbulence viscous model. Each model was simulated with different inlet or outlet boundary conditions. Our major concern for the flow filed simulation was pressure drop across the cyclone. For the particle trajectory simulation, we adopted Euler-Lagrangian approach to track particle motion from inlet to outlet of the cyclone. Particle collection efficiencies of various conditions are calculated by number based collection efficiency. The result showed that the rotation angle of the vane plays major roll to the pressure drop. But the smaller rotation angle of vane causes particle collection efficiency difference with different inlet position.

• Journal of the Society of Naval Architects of Korea
• /
• v.59 no.3
• /
• pp.183-191
• /
• 2022

A design method for a propeller type rim-driven axial-flow turbine for a micro-hydropower system is presented. The turbine consists of pre-stator, impeller and post-stator, where the pre-stator plays a role as a guide vane to provide circumferential velocity to the on-coming flow, and the impeller as a rotational power generator by absorbing angular momentum of the flow. BEM(Blade Element Method), which is based on the turbine Euler equation, is employed to design the pre-stator and impeller blades. NACA 66 thickness form and a=0.8 mean camber line, which is widely accepted as a marine propeller blade section, is used for the pre-stator and turbine blade section. A CFD method, derived from the discretization of the RANS equations, is applied for the analysis of the designed turbine system. The design conditions of the turbine is confirmed by the CFD calculation. Turbine characteristic curve is calculated by the CFD method, in order to provide the performance characteristics at off-design operation conditions. The proposed procedures for the design of a propeller type rim-driven axial-flow turbine are established and confirmed by the CFD analysis.

• Nuclear Engineering and Technology
• /
• v.54 no.10
• /
• pp.3909-3918
• /
• 2022

In the event of a severe accident, the reactor core may melt due to insufficient cooling. the high-temperature core melt will have a strong interaction (FCI) with the coolant, which may lead to steam explosion. Steam explosion would pose a serious threat to the safety of the reactors. Therefore, the study of steam explosion is of great significance to the assessment of severe accidents in nuclear reactors. This research focuses on the development of a two-dimensional steam explosion integrated analysis code called SEINA. Based on the semi-implicit Euler scheme, the three-phase field was considered in this code. Besides, the influence of evaporation drag of melt and the influence of solidified shell during the process of melt droplet fragmentation were also considered. The code was simulated and validated by FARO L-14 and KROTOS KS-2 experiments. The calculation results of SEINA code are in good agreement with the experimental results, and the results show that if the effects of evaporation drag and melt solidification shell are considered, the FCI process can be described more accurately. Therefore, it is proved that SEINA has the potential to be a powerful and effective tool for the analysis of steam explosions in nuclear reactors.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2292-2314
• /
• 2006

Efficient numerical method is developed for the prediction of aerodynamic noise generation and propagation in low Mach number flows such as aeolian tone noise. The proposed numerical method is based on acoustic/viscous splitting techniques of which acoustic solvers use simplified linearised Euler equations, full linearised Euler equations and nonlinear perturbation equations as acoustic governing equations. All of acoustic equations are forced with immersed surface dipole model which is developed for the efficient computation of aerodynamic noise generation and propagation in low Mach number flows in which dipole source, originating from unsteady pressure fluctuation on a solid surface, is known to be more efficient than quadrupole sources. Multi-scale overset grid technique is also utilized to resolve the complex geometries. Initially, aeolian tone from single cylinder is considered to examine the effects that the immersed surface dipole models combined with the different acoustic governing equations have on the overall accuracy of the method. Then, the current numerical method is applied to the simulation of the aeolian tones from twin cylinders aligned perpendicularly to the mean flow and separated 3 diameters between their centers. In this configuration, symmetric vortices are shed from twin cylinders, which leads to the anti-phase of the lift dipoles and the in-phase of the drag dipoles. Due to these phase differences, the directivity of the fluctuating pressure from the lift dipoles shows the comparable magnitude with that from the drag dipoles at 10 diameters apart from the origin. However, the directivity at 100 diameters shows that the lift-dipole originated noise has larger magnitude than, but still comparable to, that of the drag-dipole one. Comparison of the numerical results with and without mean flow effects on the acoustic wave emphasizes the effects of the sheared background flows around the cylinders on the propagating acoustic waves, which is not generally considered by the classic acoustic analogy methods. Through the comparison of the results using the immersed surface dipole models with those using point sources, it is demonstrated that the current methods can allow for the complex interactions between the acoustic wave and the solid wall and the effects of the mean flow on the acoustic waves.