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

A Bubble size distribution model has been developed for the numerical simulation of cryogenic high-speed cavitating flow of the turbo-pumps in the liquid fuel rocket engine. The new model is based on the previous one proposed by the authors, in which the bubble number density was solved as a function of bubble size at each grid point of the calculation domain by means of Eulerian framework with respect to the bubble size coordinate. In the previous model, the growth/decay of bubbles due to pressure difference between bubble and liquid was solved exactly based on Rayleigh-Plesset equation. However, the unsteady heat transfer between liquid and bubble, which controls the evaporation/condensation rate, was approximated by a theoretical solution of unsteady heat conduction under a constant temperature difference. In the present study, the unsteady temperature field in the liquid around a bubble is also solved exactly in order to establish an accurate and efficient numerical simulation code for cavitating flows. The growth/decay of a single bubble and growth of bubbles with nucleation were successfully simulated by the proposed model.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.3
• /
• pp.215-222
• /
• 2009

Jets issuing through small holes in a wall into a freestream has proven effective in the control of flow separation. This technique is known as the vortex generator jet (VGJs) method. If a precursor signal of separation is found, the separation control system using VGJs can be operated just before the onset of separation and the flow field with no separation is always attained. In this study, we measured the flow field and the wall static pressure in a two-dimensional diffuser to find a precursor signal of flow separation. The streamwise velocity measurements were carried out in the separated shear layer and spectral analysis was applied to the velocity fluctuations at some angles with respect to the diffuser. The pattern of peaks in the spectral analysis changes as the divergence angle increases over the angle of which the whole separation occurs. This change in the spectral pattern is related to the enhancement of the growth of shear layer vortices and appears just before the onset of separation. Therefore, the growth of shear layer vortices can be regarded as a precursor signal to flow separation.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.17-22
• /
• 2002

Non-uniform and unsteady inflow into a Horizontal Axis Wind Turbine (HAWT) brings about an asymmetric flow field on the rotor plane and an unsteady aerodynamic load on the blades. In the present paper effects of yawed inflow and wind shear are analyzed by an inviscid aerodynamic model based on the asymptotic acceleration potential method. In the analysis the rotor blades are represented by spanwise and chordwise pressure distribution composed of analytical first-order asymptotic solutions for the Laplace equation. As the actual wind field experienced by a HAWT is turbulent, the effects of the turbulence are also examined using the Veers' model.

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

The 2-dimensional unsteady flows on and around the cambered airfoils were computed by applying LES with the deductive dynamic SGS model. The unsteady flow field were used as inputs to compute the far-field sounds and directivity patterns from rotating blades by a hybrid approach that exploits Farassat's formula. The BEM (Boundary Element Method) was applied to predict the frequency characteristics from the rotating blades for the cases of even- and uneven-pitched fans. The BEM results suggested that the unevenly pitched fan have less pronounced discrete peaks at BEF frequencies, which was confirmed by the experiment.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1340-1346
• /
• 1999

The turbulent shear flow around a surface-mounted vertical fence was investigated using the two-frame PTV system. The Reynolds number based on the fence height(H) was 2950. From this study, it is revealed that at least 400 instantaneous velocity field data are required for ensemble average to get reliable turbulence statistics, but only 100 field data are sufficient for the time-averaged mean velocity information. Various turbulence statistics such as turbulent intensities, turbulence kinetic energy and Reynolds shear stress were calculated from 700 instantaneous velocity vector fields. The fence flow has an unsteady recirculation region behind the fence, followed by a slow relaxation to the flat-plate boundary layer flow. The time-averaged reattachment length estimated from the streamline distribution is about 11.2H. There exists a region of negative Reynolds shear stress near the fence top due to the highly convex (stabilizing) streamline-curvature of the upstream flow. The large eddy structure in the separated shear layer seems to have significant influence on the development of the separated shear layer and the reattachment process.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.5
• /
• pp.454-461
• /
• 2009

Internal aeroacoustics of an axial fan used for circulating cold air in refrigerators are computed by using the hybrid method where CFD, acoustic analogy and BEM techniques are utilized. The unsteady flow field around the axial fan is predicted by solving the incompressible RANS equations with the conventional CFD techniques. Then, main noise sources are extracted from this unsteady flow field predictions using Acoustic Analogy. Lastly, BPF noise generated from an axial fan are predicted using these modeled sources combined with the tailed Green function techniques, which are numerically solved by the BEM technique. This hybrid model is validated by comparing the prediction with the experiment. Then, parameter studies are carried out, which suggest a capability of the current method as a design tool for the low-noise of the current axial fan system in a refrigerator.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.3
• /
• pp.257-268
• /
• 2014

The objective of the present study is to investigate thermal diffusion and radiation effects on unsteady MHD flow past a linearly accelerated vertical porous plate with variable temperature and also with variable mass diffusion in presence of heat source or sink under the influence of applied transverse magnetic field. The fluid considered here is a gray, absorbing/emitting radiation but a non-scattering medium. At time t > 0, the plate is linearly accelerated with a velocity $u=u_0t$ in its own plane. And at the same time, plate temperature and concentration levels near the plate raised linearly with time t. The dimensionless governing equations involved in the present analysis are solved using the closed analytical method. The velocity, temperature, concentration, skin-friction, the rate or heat transfer and the rate of mass transfer are studied through graphs in terms of different physical parameters like magnetic field parameter (M), radiation parameter (R), Schmidt parameter (Sc), Soret number (So), Heat source parameter (S), Prandtl number (Pr), thermal Grashof number (Gr), mass Grashof number (Gm) and time (t).

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.41-46
• /
• 2006

When new weapons are introduced, the target points estimation is one of the important objectives in the flight test as well as the safe separation. The prediction methods help to design the flight test schedule. However, the incremental aerodynamic coefficients in the aircraft flow field so-called BSE are difficult to predict. Generally, the semiempirical methods such as the grid methods, IFM and Flow TGP using database are used for estimation of BSE. However, these methods are quasi-steady methods using static aerodynamic loads. Nowadays the time-accurate CFD method is often used to predict the store separation event. In the current process, the incremental aerodynamic coefficients in BSE regime are calculated directly, and the elimination of delta coefficients is checked simultaneously. This stage can be used for the initial condition of Flow TGP with freestream database. Two dimensional supersonic and subsonic store separation problems have been simulated and incremental coefficients are calculated. The results show the time when the store gets out of BSE region.

• Journal of the Korean Society of Visualization
• /
• v.8 no.1
• /
• pp.13-18
• /
• 2010

An experimental investigation of the flow around a vertical fence was carried out using a PIV velocity field measurement technique. The vertical fence was embedded in a turbulent boundary layer. The instantaneous velocity fields measured at cross-sectional planes reveal complex longitudinal vortices that vary in size and strength, developing from the upstream location. In the instantaneous vorticity and velocity field data, the shear flow separated from the fence top is highly turbulent and shows unsteady flow characteristics. The topography of the ensemble averaged velocity fields, especially the separation bubble formed behind the fence, shows that the spatial distributions of streamwise velocity (U) and vertical (V) are symmetric, the spanwise velocity (W) is skew-symmetric with respect to the central xy-plane(z=0).

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2012.05a
• /
• pp.309-317
• /
• 2012

Recent years have witnessed the use of micro shock tube in various engineering applications like micro combustion, micro propulsion, particle delivery systems etc. The flow characteristics occurring in the micro shock tube shows a considerable deviation from that of well established conventional macro shock tube due to very low Reynolds number and high Knudsen number effects. Also the diaphragm rupture process, which is considered to be instantaneous process in many of the conventional shock tubes, will be crucial for micro shock tubes in determining the near diaphragm flow field and shock formation. In the present study, an axi-symmetric CFD method has been applied to simulate the micro shock tube, with Maxwell's slip velocity and temperature jump boundary conditions. The effects of finite diaphragm rupture process on the flow field and the shock formation was investigated, in detail. The results show that the shock strength attenuates rapidly as it propagates through micro shock tubes.