• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.11
• /
• pp.1073-1079
• /
• 2009

The proper orthogonal decomposition (POD) method can identify principal modes that optimally capture the energy content from large multi-dimensional data set. In this study unsteady pressure fields on the rotor blade surface of a helicopter in forward flight are expressed by a reduced order model based on the POD method. Special modes containing high energy are analyzed to investigate the aerodynamic characteristics in more efficient way. The CFD simulation of flowfields around helicopter rotor blade in hovering motion is also conducted to validate its prediction with experimental result. In the process 7 modes containing energy ratio 99% from 240 snapshots information are identified and utilized to construct a reduced order model.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.16 no.11
• /
• pp.1077-1083
• /
• 2004

This study is to systematically analyze the effect of various kinds of design factors on the performance of stratified thermal storage tank. Taguchi method, known as a very reasonable tool in the parametric study, is employed in the present work. Three dimensional unsteady numerical experiment is conducted for 4 design parameters of stratified thermal storage tank: inlet Reynolds number, Froude number, diffuser size d with 3 levels (Re=400, 800, 1200, Fr=0.5, 1.0, 2.0 and d=150 mm, 200mm, 300 mm) and diffuser shape with 2 levels. Orthogonal array $L_{18}(2{\times}3^7)$ is adopted for the analysis of variance. The result gives quantitative estimation of the various design parameters affecting the performance and helps to select the main factors for the optimum design of stratified thermal storage tank. Reynolds number is found to be the most dominant parameter and the diffuser shape plays significant role on the performance of stratified thermal storage tank. Based on this finding, the prior questions on the contribution of the diffuser shape proposed by the authors become clear. The optimum condition for the performance is a set of d=300mm, Re=800, and radial regulated plate diffuser. Conformation test shows the repeatability in the analysis and $1.3\%$ difference between the estimated thermocline thickness and that of numerical result.

• Transactions of the Korean hydrogen and new energy society
• /
• v.27 no.5
• /
• pp.547-553
• /
• 2016

This paper presents the effect of a grid system on the performance of a small Savonius wind turbine installed side-by-side. Turbine performance is compared using three different grid systems; tetrahedral grid having a concentrated circular grid around turbine rotors, the tetrahedral grid having a concentrated rectangular grid around turbine rotors and the symmetric grid having a concentrated tetrahedral grid near the turbine rotor blades and a hexahedral grid. The commercial code, SC/Tetra has been used to solve the three-dimensional unsteady Reynolds-averaged Navier-Stokes analysis in the present study. The Savonius turbine rotor has a rotational diameter of 0.226m and an aspect ratio of 1.0. The distance between neighboring rotor tips keeps the same length of the rotor diameter. The variations of pressure and power coefficient are compared with respect to blade rotational angles and rotating frequencies of the turbine blade. Throughout the comparisons of three grid systems, it is noted that the symmetric grid having a concentrated tetrahedral grid near the turbine rotor blades and a hexahedral grid has a stable performance compared to the other ones.

• Fire Science and Engineering
• /
• v.29 no.5
• /
• pp.42-50
• /
• 2015

Three simulation approaches for turbulence were applied for the computation of propane dispersion in a simplified real-scale urban area with one building:, Large Eddy Simulation (LES), Detached Eddy Simulation (DES), and Unsteady Reynolds Averaged Navier-Stokes (RANS). The computations were performed using FLUENT 14, and the grid system was made with ICEM-CFD. The propane distribution depended on the prediction performance of the three simulation approaches for the eddy structure around the building. LES and DES showed relatively similar results for the eddy structure and propane distribution, while the RANS prediction of the propane distribution was unrealistic. RANS was found to be inappropriate for computation of the gas dispersion process due to poor prediction performance for the unsteady turbulence. Considering the computational results and cost, DES is believed to be the optimal choice for computation of the gas dispersion in a real-scale space.

• Fire Science and Engineering
• /
• v.13 no.1
• /
• pp.37-47
• /
• 1999

The spreading fire of very small floating particles after they are ignited is fast and t therefore dangerous. The research on this area has been limited to experiments and global simulations which treat them as dusts or gaseous fuel with certain concentration well m mixed with air. This research attempted micro-scale analysis of ignition of those particles modeling them as liquid droplets. For the beginning, the in-line array of fuel droplets is modeled by two-dimensional, unsteady conservation equations for mass, momentum, energy and species transport in the gas phase and an unsteady energy equation in the liquid phase. They are solved numerically in a generalized non-orthogonal coordinate. The single step chemical reaction with reaction rate controlled by Arrhenius’ law is assumed to a assess chemical reaction numerically. The calculated results show the variation of temperature and the concentration profile with time during evaporation and ignition process. Surrounding oxygen starts to mix with evaporating fuel vapor from the droplet. When the ignition condition is met, the exothermic reactions of the premixed gas initiate a and burn intensely. The maximum temperature position gradually approaches the droplet surface and maximum temperature increases rapidly following the ignition. The fuel and oxygen concentration distributions have minimum points near the peak temperature position. Therefore the moment of ignition seems to have a premixed-flame aspect. After this very short transient period minimum points are observed in the oxygen and fuel d distributions and the diffusion flame is established. The distance between droplets is an important parameter. Starting from far-away apart, when the distance between droplets decreases, the ignition-delay time decreases meaning faster ignition. When they are close and after the ignition, the maximum temperature moves away from the center line of the in-line array. It means that the oxygen at the center line is consumed rapidly and further supply is blocked by the flame. The study helped the understanding of the ignition of d droplet array and opened the possibility of further research.

• Nuclear Engineering and Technology
• /
• v.44 no.3
• /
• pp.283-296
• /
• 2012

The natural convection in a horizontal enclosure heated from the bottom wall, cooled at the top wall, and having a square adiabatic body in the center is studied. Three different Prandtl numbers (0.01, 0.7 and 7) are considered for the investigation of the effect of the Prandtl number on natural convection. Adiabatic boundary conditions are employed for the side walls. A two-dimensional solution for unsteady natural convection is obtained, using an accurate and efficient Chebyshev spectral methodology for different Rayleigh numbers varying over the range of $10_3$ to $10_6$. It had been experimentally reported that the heat transfer mode becomes oscillatory when Pr is out of a specific Pr band beyond the critical Ra. In this study, we reproduced this phenomenon numerically. It was found that when Ra=$10_6$, only the case for intermediate Pr (=0.7) reached a non-changing steady state and the low and high Pr number cases (Pr=0.01 and 7) showed a periodically oscillatory fashion hydrodynamically and thermally. The variation of time- and surface-averaged Nusselt numbers on the hot and cold walls for different Rayleigh numbers and Prandtl numbers are presented to show the overall heat transfer characteristics in the system. Further, the isotherms and streamline distributions are presented in detail to compare the physics related to their thermal behavior.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.119-122
• /
• 2005

Ablation phenomena is very complicated because it includes momentum, energy and mass transfer, chemical reactions as well as phase change. In this paper, ablation at the rocket nozzle throat is modeled as unsteady one dimensional axi-symmetric with proper boundary conditions and field equation is solved numerically. Analytical results are compared with measured ablation data from firing experimental liquid rocket engine. Test variables are combustion pressure and mixture ratio. for low combustion pressure and low mixture ratio, the discrepancy between analysis and experiments are large but for the normal rocket operation range, two results show a simliar trend with maximum discrepancy of $100\%$.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.26-32
• /
• 2006

This paper targets a direct and quantitative prediction of characteristics of unstable waves in a combustion chamber, which employs the governing equations derived in terms of amplification factors of flow variables. A freshly formulated nonlinear acoustic equation is obtained and the analysis of unsteady waves in a rocket engine is attempted. In the present formalism, perturbation method decomposes the variables into time-averaged part that can be obtained easily and accurately and time-varying part which is assumed to be harmonic. Excluding the use of conventional spatially sinusoidal eigenfunctions, a direct numerical solution of wave equation replaces the initial spatial distribution of standing waves and forms the nonlinear space-averaged terms. Amplification factor is also calculated independently by the time rate of changes of fluctuating variables, and is no longer an explicit function for compulsory representation. Employing only the numerical computation, major assumptions inevitably inherent, and in erroneous manner, in up to date analytical methods could be avoided. With two definitions of amplification factor, 1-D stable wave and 3-D unstable wave are examined, and clearly demonstrated the potentiality of a suggested theoretical-numerical method of combustion instability.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.32-38
• /
• 2006

The objective of this study is to evaluate the prediction accuracy of development large eddy simulation(LES) program for turbulent flow behind a bluff-body. The LES solver was implemented on parallel computer consisting 16 processors. To verify the capability of LES code, the results were compared with those of Reynolds Averaged Navier-Stokes(RANS) using standard ${\kappa}-{\varepsilon}$ model as well as experimental data. The results showed that the LES and RANS qualitatively well predicted the experimental results, such as mean axial, radial velocities and turbulent kinetic energy. In the quantitative analysis, however, the LES showed a better prediction performance than RANS. Specially, the LES well described characteristics of the recirculation zones, such as air stagnation point and jet stagnation point. Finally, the unsteady phenomena on the Bluff-body, such as the transition of recirculation region and vorticity, was examined with LES methodology.

• Journal of ILASS-Korea
• /
• v.5 no.1
• /
• pp.13-22
• /
• 2000

The present study has numerically analyzed the vaporization characteristics of fuel droplets in the high temperature convective flow field. The axisymmetric governing equations for mass, momentum, energy, and species are solved by an iterative and implicite unstructured finite-volume method. The moving boundary due to vaporization is handled by the deformable unstructured grid technique. The pressure-velocity coupling in the density-variable flows is treated by the SIMPLEC algorithm. In terms of the matrix solver, Bi-CGSTAB is employed for the numerically efficient and stable convergence. The n-decane is used as a liquid fuel and the initial droplet temperature is 300K. Computations are performed for the nonevaporating and evaporating droplets with the relative interphase velocity(25m/s). The unsteady vaporization process has been simulated up to the nondimensional time, 25. Numerical results indicate that the mathematical model developed in this study succesfully simulates the main features of the droplet vaporization process in the convective environment.