• Journal of The Korean Astronomical Society
• /
• v.53 no.2
• /
• pp.49-57
• /
• 2020

We present numerical simulations of decaying hydrodynamic turbulence initially driven by solenoidal (divergence-free) and compressive (curl-free) drivings. Most previous numerical studies for decaying turbulence assume an isothermal equation of state (EOS). Here we use a polytropic EOS, P ∝ ρ^γ, with polytropic exponent γ ranging from 0.7 to 5/3. We mainly aim at determining the effects of γ and driving schemes on the decay law of turbulence energy, E ∝ t^-α. We additionally study probability density function (PDF) of gas density and skewness of the distribution in polytropic turbulence driven by compressive driving. Our findings are as follows. First of all, we find that even if γ does not strongly change the decay law, the driving schemes weakly change the relation; in our all simulations, turbulence decays with α ≈ 1, but compressive driving yields smaller α than solenoidal driving at the same sonic Mach number. Second, we calculate compressive and solenoidal velocity components separately and compare their decay rates in turbulence initially driven by compressive driving. We find that the former decays much faster so that it ends up having a smaller fraction than the latter. Third, the density PDF of compressively driven turbulence with γ > 1 deviates from log-normal distribution: it has a power-law tail at low density as in the case of solenoidally driven turbulence. However, as it decays, the density PDF becomes approximately log-normal. We discuss why decay rates of compressive and solenoidal velocity components are different in compressively driven turbulence and astrophysical implication of our findings.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.259-267
• /
• 2017

Flow fields around a KARI-11-180 airfoil, SDM and transonic body are numerically simulated by using an unstructured meshes based compressible flow solver developed at KAIST. RANS equations are solved to analyse the flow fields and Roe's FDS method is adopted to evaluate convective fluxes. Turbulence effect of the flow fields is modeled by a SA model, SST model and ${\gamma}-{\widetilde{Re}}_{{\theta}t}$ model. It is found that smaller drag coefficients are predicted for the KARI-11-180 airfoil when a transition phenomenon is considered and small deviations exist between CFD and EFD results. For the SDM, flow separation is observed at a leading edge and calculated aerodynamic properties show similar tendencies to experimental results. A shock wave on main wings of the transonic body is successfully captured by the present flow solver at a Mach number 0.9. Estimated pressure profiles by means of the present CFD method also agree well with those of wind tunnel results.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.66-71
• /
• 2013

In this paper, the supersonic flows around space launch vehicles have been numerically simulated by using a 3-D RANS flow solver. The focus of the study was made for investigating plume-induced flow separation(PIFS). For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras model was employed for the closure of turbulence. The Gauss-Seidel iteration was used for time integration. To validate the flow solver, calculation was made for the 0.04 scale model of the Saturn-5 launch vehicle at the supersonic flow condition without exhaust plume, and the predicted results were compared with the experimental data. Good agreements were obtained between the present results and the experiment for the surface pressure coefficient and the Mach number distribution inside the boundary layer. Additional calculations were made for the real scale of the Saturn-5 configuration with exhaust plume. The flow characteristics were analyzed, and the PIFS distances were validated by comparing with the flight data. The KSLV-1 is also simulated at the several altitude conditions. In case of the KSLV-1, PIFS was not observed at all conditions, and it is expected that PIFS is affected by the nozzle position.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.169-172
• /
• 2002

When a plane shockwave reflects ken a concave wall, it is focused at a certain location, resulting in extremely high local pressure and temperature. This focusing is due to a nonlinear phenomenon of shock wave. The focusing phenomenon has been extensively applied to many diverse folds of engineering and medical treatment as well. In the current study, the focusing of shock wave over a reflector is numerically investigated using a CFD method. The Harten-Yee total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The incident shock wave Mach number $M_{s}\;of\;1.1{\~}l.3$ is applied to the parabolic reflectors with several different depths. Detailed focusing characteristics of the shock wave are investigated in terms of peak pressure, gasdynamic and geometrical foci. The results obtained are compared with the previous experimental results. The results obtained show that the peak pressure of shock wave focusing and its location strongly depend on the magnitude of the incident shock wave and depth of parabolic reflector. It is also found that depending up on the depth of parabolic reflector, the weak shock wave focusing process can classified into three distinct patterns : the reflected shock waves do not intersect each other before and after focusing, the reflected shock waves do not intersect each other before focusing, but intersect after focusing, and the reflected shock waves intersect each other before and after focusing. The predicted Schlieren images represent the measured shock wave focusing with a good accuracy.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.20 no.1
• /
• pp.20-27
• /
• 2016

The prediction of combustion instability is important to avoid an obvious threat to the structural safety and the motor performance because it affects the apparent response function of the propellant, the burning rate, and a mean flow Mach number at the local surface. The combustion instability occurs in case acoustic waves were coupled with the combustion/flow dynamic frequency. In this paper, an acoustic instability model is derived from the nonlinear wave equation for analysing acoustic dynamics in solid rocket motors. The chamber pressure and burning rate effects on combustion instability have been investigated.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.11a
• /
• pp.93-98
• /
• 2007

The air velocity flowing in inner combustion chamber of SCramjet is supersonic and the time of its stay is very short as a few milliseconds. Within this short time, fuel injection， air-fuel mixing, and combustion process should be accomplished. Several methods are suggested for mixing enhancement. Among these, cavity is selected to study for mixing characteristics. The numerical simulation is performed in the case of freestream Mach number of 2.5 and cavity located in front of fuel jet injection. 3 different sized cavities of the same length-height ratio were used in order to recognize the effect about cavity size. Also, the case without cavity was analyzed to find the effect of cavity. Used code compared with the result of experiment under identical conditions and it was verified. Through this comparison and verification, mixing enhancement by cavity size could be confirmed.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.11a
• /
• pp.350-353
• /
• 2007

In this study, the genetic algorithm has been used for the real-time defect diagnosis on the operation of the aircraft gas-turbine engine. The component elements of the gas-turbine engine for consideriation of the performance deterioration is consist of the compressor, the gas generation turbine and the power turbine, repectively. Compared to the on-design point on the sea-level condition, the learning data has been increased 200 times in case of the off-design conditions for the altitude, the flight mach number and the fuel consumption. Therefore, enormous learning time has been required for the satisfied convergence. The optimum division has been proposed to decrease learning time as well as to obtain high accuracy. As results, the RMS errors of the defect diagnosis using the genetic algorithm have been estimated under 5 %.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.1
• /
• pp.9-17
• /
• 2013

A considerable deal of work has been carried out to get an insight into the gas-solid suspension flows and to specify the particle motion and its influence on the gas flow field. In this paper an attempt is made to develop an analytical model to study the effect of nozzle inlet/exit pressure ratio, particle/gas loading and the particle diameter effect on gas-solid suspension flow. The effect of the particle/gas loading on the mass flow, Mach number, thrust coefficient and static pressure variation through the nozzle is analyzed. The results obtained show that the presence of particles seems to reduce the strength of the shock wave. It is also found that smaller the particle diameter is, bigger will be the velocity as bigger particle will have larger slip velocity. The suspension flow of smaller diameter particles has almost same trend as that of single phase flow with ideal gas as working fluid. Depending on the ambient pressure, the thrust coefficient is found to be higher for larger particle/gas loading or back pressure ratio.

• Transactions of the Korean hydrogen and new energy society
• /
• v.20 no.1
• /
• pp.31-37
• /
• 2009

An ejector is a machine utilized for mixing fluid, maintaining a vacuum, and transporting fluid. The Ejector enhances system efficiency, are easily operated, have a mechnically simple structure, and do not require a power supply. Because of these advantages, the ejector has been applied to a variety of industrial fields such as refrigerators, power plants and oil plants. In this work, an ejector was used to safely recycle anode tail gas in a 5 kW Molten Carbonate Fuel Cell system at KEPRI(Korea Electric Power Research Institute). In this system, the ejector is placed at mixing point between the anode tail gas and the cathode tail gas or the fresh air. Commercial ejectors are not designed for the actual operating conditions for our fuel cell system. A new ejector was therefore designed for use beyond conventional operating limits. In this study, the entrainment ratio is measured according to the diametrical ratio of nozzle to throat in the designed ejector. This helps to define important criteria of ejectors for MCFC recycling.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.42 no.12
• /
• pp.1028-1036
• /
• 2014

A design procedure is established for hypersonic nozzles by using MOC(Method of Characteristics) and CFD. The inviscid nozzle contour is designed by MOC, then BLC(Boundary Layer Correction) is made by evaluating the boundary layer thickness from viscous CFD analysis. By comparing various definitions of the boundary layer thicknesses, it seems that the boundary layer thickness of 95% speed of the maximum value at the cross section satisfies best the design Mach number. Design procedure is as follow; MOC design, grid generation, inviscid analysis, viscous analysis, BLC and viscous analysis for confirmation and post-processing. All procedures are made automatically by using the batch processing.