• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2005.11a
• /
• pp.46-55
• /
• 2005

For the sea-level performance test of rocket motor designed to operate in the upper atmosphere, ejectors with no induced secondary flow are generally used, which serves dual purposes of evacuating the test cell and performing as a supersonic exhaust diffuser (SED). The main concern of this research is to simulate starting transients in order to visualize evolution of internal shock structures in SED with different initial cell (vacuum chamber) pressures. RANS code with low Reynolds $k-\varepsilon$ turbulence model was employed for these computations. Numerical results were compared with the pressure measurements previously performed [Proceedings of 2004 Annual Conference, KIMST], and showed good agreements with pressure-time history of measured data. In the case of low vacuum chamber pressure, abrupt impingement of the under-expanded supersonic jet from the nozzle onto the diffuser wall was observed, whereas initial impingement point was located downstream and moved slowly upstream in the case of non-vacuum chamber pressure. In spite of initially dissimilar evolution of shock structures, iso-mach contour revealed that the steady shock structures had little difference except the location of flow separation and normal shock.

• Journal of Industrial Technology
• /
• v.22 no.B
• /
• pp.27-34
• /
• 2002

3-D numerical studies are performed to investigate the effect of the air dam height and approaching air velocities on the pressure distribution of notchback road vehicle. For this purpose, the models of test vehicle with four different air dam heights are introduced and PHOENICS, a commercial CFD code, is used to simulate the flow phenomena and to estimate the values of pressure coefficients along the surface of vehicle. The standard $k-{\varepsilon}$ model is adopted for the simulation of turbulence. The numerical results show that the height variation of air dam makes almost no influence on the distribution of the value of pressure coefficient along upper and rear surface but makes strong effects on the bottom surface. That is, the value of pressure coefficient becomes smaller as the height is increased along the bottom surface. Approaching air velocity makes no differences on pressure coefficients. Through the analysis of pressure coefficient on the vehicle surface, one tries to assess aerodynamic drag and lift of vehicle. The pressure distribution on the bottom surface affects more on lift than the pressure distribution on the upper surface of the vehicle does. The increase of air dam height makes positive effects on the lift decrease but no effects on drag reduction.

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.5
• /
• pp.18-23
• /
• 2011

A CFD analysis method is developed and applied for investigating the gas flow and the byproduct particle trajectory in Roots type vacuum pump. The internal fluid flow and thermal fields between the rotors and the housing of vacuum pump are analyzed by using the dynamic mesh, the numerical methods for unsteady 2-D Navier-Stokes equation and the standard k-$\varepsilon$ turbulence model of the Fluent code. Coupled with the flow simulation results, the particle trajectory of the byproduct flowing into the pump with gas stream is analyzed by using discrete phase modeling technique. The CFD analysis results show the pressure, the velocity and the temperature distributions in pump change abruptly due to the rotation of rotors, and back flows are produced due to the strong reverse pressure gradients at rotor/rotor and rotor/housing clearances. The predicted byproduct particle trajectory results also show the particles impinge on the clearance surfaces between the housing and the rotor of pump and then may form the deposit layer causing the failure of pump.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.568-572
• /
• 2000

A numerical investigation was performed to determine the effect of airfoil thickness on the optimum Gurney flap height using NACA 00XX series airfoils. Seven airfoils which have 3% chord thickness difference were used. These were NACA 0006, 0009, 0012, 0015, 0018, 0021, and 0024. A Navier-Stokes code, FLUENT, was used to calculate the flow field about airfoil. The fully turbulent results were obtained using the standard $k-{\varepsilon}$ two-equation turbulence model. To provide a check case fur our computational method, numerical studies for NACA 4412 airfoil were made and compared with already existing experimental data for this airfoil by Wadcock. For every NACA 00XX airfoil, Gurney flap heights ranging from 0.5% to 2.0% chord were changed by 0.5% chord interval and their effects were studied. With the numerical solutions, the relationship between $(L/D)_{max}$ and airfoil thickness as a function of flap height and the relationship between $(L/D)_{max}$ and flap height as a function of airfoil thickness were investigated. The same relationship for $(C_l)_{max}$ also were shown. From these results, the optimum flap size for each airfoil thickness can be determined and vice versa.

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

A three-dimensional coupled turbulent fluid flow and solidification process were analyzed in a continuous casting process of a steel slab with Electromagnetic Brake(EMBR). A revised low-Reynolds number $k-{\varepsilon}$ turbulence model was used to consider the turbulent effects. The enthalpy-porosity relation was employed to suppress the velocity within a mushy region. The electromagnetic field was described by Maxwell equations. Tile application of EMBR to the mold region results in the decrease of the transfer of superheat to the narrow face, the increase of temperature in free surface region and most liquid of submold region, and the higher temperature gradient near the solidifying shell. The increasing magnetic flux density effects mainly to the surface temperature of the solidifying shell of narrow face, hardly to the one of wide face. It is seen that in the presence of EMBR a thicker solidifying shell is obtained at the narrow face of the slab.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.8-13
• /
• 2007

A general purpose program NUFLEX for the analysis 3-D thermo/fluid flow and pre/post processor in complex geometry has been developed, which consists of a flow solver based on FVM and GUI based pre/post processor. The solver employs a general non-orthogonal grid system with structured grid and solves laminar and turbulent flows with standard/RNG $k-{\varepsilon}$ turbulence model. In addition, NUFLEX is incorporated with various physical models, such as interfacial tracking, cavitation, MHD, melting/solidification and spray models. For the purpose of evaluation of the program and testing the applicability, many actual problems are solved and compared with the available data. Comparison of the results with that by STAR-CD or FLUENT program has been also made for the same flow configuration and grid structure to test the validity of NUFLEX.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.811-814
• /
• 2002

Using a mathematical theory, we show that the optimality condition of a turbulent diffuser with maximum pressure recovery at the exit is zero shear stress along the wall. The optimal diffuser shape is designed through iterative procedures by using the $k-{\varepsilon}-{\nu}^{2}-f$ turbulence model for flow simulation. The Reynolds number based on the bulk mean velocity and the channel height at the diffuser entrance is 18,000. We also perform large eddy simulation to validate the shape design results and investigate the flow characteristics near the zero-skin friction wall. Results from large eddy simulation show that the skin friction is slightly higher than zero but is still very small as compared to that of the flat plate boundary layer flow Although the time-averaged wall shear stress is slightly above zero along the diffuser wall, instantaneous flow reversals occur intermittently. The streamwise mein velocity shows an asymptotic behavior of the half-power-law near the wall where the skin friction is close to zero.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.934-939
• /
• 2001

In computational study of the flow in piston engines and the flow through moving valves, the use of moving vertices is essential for modelling flows with moving boundaries. The positions of cell vertices in such cases must be allowed to vary with time. To simulate 3-dimensional port-valve and piston-cylinder of HIMSEN 6H21/32 engine, a commercially available code, STAR-CD, was used. Changes in mesh geometry was specified by PROSTAR commands.(i.e. the Change Grid operation in the EVENTS command module.) Control of the intake flow is expected to play an important role as designers seek to obtain better fuel spray characteristics, fuel mixing and mixture preparation, combustion performance, and emissions reductions to meet national standards. As a result of analysis, velocity fields indicate the presence of a structured flow comprised of one pair of counter-rotating vortices under the intake valve during the early induction process. These flow structures remain visible for most of the intake process. As the piston moves towards BDC, these vortices develops into a larger tumbling motion that dominates the flow structure.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.8 no.3
• /
• pp.181-189
• /
• 2000

The aerodynamic characteristics of large transport vehicle has become more and more important in recent vehicle design to improve driving performance in high speed cruising and raise the product valve with regard to a comfortable driving condition. Hence, detailed knowledge of the flow field around truck coner vane is essential to improve fuel efficiency and reduce the dirt contamination on vehicle body surface. In this study, three-dimensional flow characteristics around corner vane attached to truck cabin were computed for the steady, incompressible, and high speed viscous flow, adopting the RNG k-$\varepsilon$ turbulence model. In order to investigate the influence of configuration and structure of corner vane, computations were carried out for four cases at a high Reynolds number, Re=4.1$\times$106 (based on the cabin height). The global flow patterns, drag coefficient and the distributions such as velocity magnitude, turbulent kinetic energy around the corner vane, were examined. As a result of this study, we could identify the flow characteristics around corner vane for the variation of corner vane length and width. Also, suggest the improved structure to reduce the dirt contamination in cabin side.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.155-161
• /
• 1999

The three-dimensional turbulent cascade flows with and without endwall fences are numerically investigated by solving the incompressible Navier-Stokes equations with a high-Reynolds number $k-{\varepsilon}$ turbulence closure model. A projection method based algorithm is used in the finite-volume formulation, with the second order upwind-differencing scheme for the convective terms. First, assessments on accuracy of the present method are made by comparing the static pressure distributions at the mid-span of the cascade with measured data, and also by confirming the experimental observations on the choice of an optimal fence height for the secondary flow control. In understanding the three-dimensional nature of the secondary flow in turbine cascade, the limiting streamline patterns and the static pressure contours at the suction surface of the blade as well as on the cascade endwall are employed to visualize the effectiveness of the endwall fence for the secondary flow control. Analysis on the streamwise vorticity contour maps along the cascade with the three-dimensional representation of their iso-surfaces reveals the strucuture of the complicated vortical flow in the turbine cascade with endwall fence, and also leads to an understanding on formation of the counter-rotating streamwise vortex over the endwall fence, in explaining the mechanisms of controlling the secondary flow and also for the proper selection of an optimal fence height.