• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.05a
• /
• pp.57-60
• /
• 2006

The instability analysis of submerged gas flow into liquid is studied, which assumes gas and liquid as viscous and irrotational. At low mass flow rate of gas, injected gas plume is collection of bubbles, and increase of gas flow rate makes plume as a jet. It is well known that the transition from bubbling to jetting occurs in the transonic region. But previous works neglect viscous effect of gas flow into liquid. This paper concerns about an application of viscous potential flow theory in cylindrical gas flow into liquid. The growth rate versus wave number and mach number is compared with various condition including inviscid and viscous flow.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.1-6
• /
• 2013

Missile and fighter aircraft have been challenged by low restoring nose-down pitching moment at high angle of attach. The consequence of weak nose-down pitching moment can be resulting in a deep stall condition. Especially, the pressure oscillation has a huge effect on noise generation, structure damage, aerodynamic performance and safety, because the flow has strong unsteadiness at high angle of attack. In this paper, the unsteady aerodynamics coefficients were analyzed at high angle of attack up to 50 degrees around two dimensional NACA0012 airfoil. The two dimensional unsteady compressible Navier-Stokes equation with a LES turbulent model was calculated by OHOC (Optimized High-Order Compact) scheme. The flow conditions are Mach number of 0.3 and Reynolds number of $10^5$. The lift, drag, pressure, entropy distribution, etc. are analyzed according to the angle of attack. The results of average lift coefficients are compared with other results according to the angle of attack. From a certain high angle of attack, the strong vortex formed by the leading edge are flowing downstream as like Karman vortex around a circular cylinder. The primary and secondary oscillating frequencies are analyzed by the effects of these unsteady aerodynamic characteristics.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.3
• /
• pp.201-209
• /
• 2013

Much numerical and experimental research has been done for the flow around an oscillating airfoil. The main research topics are vortex shedding, dynamic stall phenomenon, MAV's lift and thrust generation. Until now, researches mainly have been concentrated on analyzing the wake flow for the variation of frequency and amplitude at a low angle of attack. In this study, wake structures and acoustic wave propagation characteristics were studied for a plunging airfoil at high angle of attack. The governing equations are the Navier-Stokes equation with LES turbulence model. OHOC (Optimized High-Order Compact) scheme and 4th order Runge-Kutta method were used. The Mach number is 0.3, the Reynolds number is, and the angle of attack is from $20^{\circ}$ to $50^{\circ}$. The plunging frequency and the amplitude are from 0.05 to 0.15, and from 0.1 to 0.2, respectively. Due to the high resolution numerical method, wake vortex shedding and pressure wave propagation process, as well as the propagation characteristics of acoustic waves can be simulated. The results of frequency analysis show that the flow has the mixed characteristics of the forced plunging frequency and the vortex shedding frequency at high angle of attack.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.4
• /
• pp.71-78
• /
• 2005

The vaporization characteristics of a liquid heptane droplet in a supercritical nitrogen flow is numerically analyzed. The present model can account for real gas effects, liquid-phase internal circulation, variable thermodynamic properties and high-pressure effects. Time marching method with preconditioning scheme is employed to handle the low Mach number flows in dense heptane droplet region. Computations are made for the wide range of convective velocity and ambient pressure. Numerical results indicate that the droplet deformation becomes stronger by increasing the Reynolds number and it becomes relatively weak by increasing the pressure.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.8-15
• /
• 2005

Aerodynamic sounds generated by a uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method. Thethird-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives, and the second-order-accurate Runge-Kutta scheme is applied for the time marching. We have succeed in capturing very small pressure fluctuations with the same frequency of the Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. It is also apparent that the amplitude of sound pressure is proportional to r /sup -1/2/,r being the distance from the center of the circular cylinder. To investigate the effect of the lattice dependence, furthermore, 2D computations of the tone noises radiated by a square cylinder and NACA0012 with a blunt trailing edge at high incidence and low Reynolds number are also investigate.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2145-2150
• /
• 2003

Ejector system is one of fluid machinery which can entrain the fluid in low pressure part and transport it to the higher pressure part. The ejector system has been widely used for the purpose of obtaining high-vacuum state, fluid transport, thrust augmentation, etc. It can transport a large capacity of fluid with relatively small device of no any moving parts, and thus seldom causes mechanical troubles. However, the conventional ejector system has been pointed out that its overall efficiency is quite low compared with other fluid machinery since it is derived by only the pure shear stresses. In the present study, 4, 6, and 8 lobed petal nozzles with a design Mach number of 1.7 are adopted as a primary nozzle to improve the ejector performance, and are compared with a conventional circular nozzle. The static pressures along the diffuser wall are measured to qualify the flow field inside the supersonic petal ejector system.

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

In this paper a conventional approach for design and analysis of subsonic air vehicle is used. First of all subsonic aerodynamic coefficients are calculated using Computational Fluid Dynamics(CFD) tools and then wind-tunnel model was developed that integrates vehicle components including control surfaces and initial data is validated as well as refined to enhance aerodynamic efficiency of control surfaces. Experimental data and limited computational fluid dynamics solutions were obtained over a Mach number range of 0.5 to 0.8. The experimental data show the component build-up effects and the aerodynamic characteristics of the fully integrated configurations, including control surface effectiveness. The aerodynamic performance of the fully integrated configurations is comparable to previously tested subsonic vehicle models. Mathematical model of the dynamic equations in 6-Degree of Freedom(DOF) is then simulated using MATLAB/SIMULINK to simulate trajectory of vehicle. Effect of altitude on range, Mach no and stability is also shown. The approach presented here is suitable enough for preliminary conceptual design. The trajectory evaluation method devised accurately predicted the performance for the air vehicle studied. Formulas for the aerodynamic coefficients for this model are constructed to include the effects of several different aspects contributing to the aerodynamic performance of the vehicle. Characteristic parameter values of the model are compared with those found in a different set of similar air vehicle simulations. We execute a set of example problems which solve the dynamic equations to find the aircraft trajectory given specified control inputs.

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

• Aerospace Engineering and Technology
• /
• v.4 no.2
• /
• pp.27-35
• /
• 2005

In this study, a performance analysis code was developed for the off-design performance analysis of air turbo ramjet(ATR) engine, and the analyses were performed for the pre-designed ATR engine at several operating points in the envelope. Variable intake and thrust nozzle were assumed to cover the wide envelope. Mathematical models for each components were developed to calculate their off-design performance. Simple design formulas were introduced for some components to explore the performance variation versus the design parameters. As a result, the pre-defined engine couldn't cover the entire mission profile. And it was also found that the effect of the pre-cooler was not very great, especially in the region of low Mach number.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.