• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.653-662
• /
• 2007

Due to the difficulty raised from the coupling of cavitation modeling with turbulent flow, numerical simulation for two phase flow remains one of the challenging issues in the society. This research focuses on the development of numerical code to deal with incompressible two phase flow around 2D hydrofoil by combing the cavitation model suggested by Kunz et al. with $k-{\varepsilon}$ turbulent model. The simulation results are compared to experimental data to verify the validity of the developed code. Also, the comparison of the calculation results is made with LES results to evaluate the capability of $k-{\varepsilon}$ turbulence model. The calculation results show very good agreement with experimental observations even though this code can not grasp the small scaled bubbles in the calculation wheres LES can hold the real physics. This code will be extended to 3D compressible two phase flow for the study on the fluid dynamics in the inducers and impellers.

• Journal of computational fluids engineering
• /
• v.13 no.1
• /
• pp.35-40
• /
• 2008

Generally flight vehicles have many cavities such as wheel wells, bomb bays and windows on their external surfaces and the flow around these cavities makes separation, vortex, shock and expansion waves, reattachment and other complex flow phenomenon. The flow around the cavity makes abnormal and three-dimensional noise and vibration even thought the aspect ratio (L/D) is small. The cavity giving large effects to the flow might make large noise, cause structural damage or breakage, harm the aerodynamic performance and stability, or damage the sensitive devices. In this study, numerical analysis was performed for cavity flows by the unsteady compressible three dimensional Reynolds-Averaged Navier-Stokes (RANS) equations with Wilcox's $\kappa-\omega$ turbulence model. The MPI(Message Passing Interface) parallelized code was used for calculations by PC-cluster. The cavity has the aspect ratios of 2.5, 3.5 and 4.5 with the W/D ratio of 2 for three-dimensional cavities. The Sound Pressure Level (SPL) analysis was done with FFT to check the dominant frequency of the cavity flow. The dominant frequencies were analyzed and compared with the results of Rossiter's formula and Ahuja& Mendoza's experimental datum.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1154-1164
• /
• 2002

Losses on the turbine consist of the mechanical loss, tip clearance loss, secondary flow loss and blade profile loss etc.,. More than 60 % of total losses on the turbine is generated by the two latter loss mechanisms. These losses are directly related with the reduction of turbine efficiency. In order to provide a new design methodology for reducing losses and increasing turbine efficiency, a two-dimensional axial-type turbine blade shape is modified by the optimization process with two-dimensional compressible flow analysis codes, which are validated by the experimental results on the VKI turbine blade. A turbine blade profile is selected at the mean radius of turbine rotor using on a heavy duty gas turbine, and optimized at the operating condition. Shape parameters, which are employed to change the blade shape, are applied as design variables in the optimization process. Aerodynamic, mechanical and geometric constraints are imposed to ensure that the optimized profile meets all engineering restrict conditions. The objective function is the pitchwise area averaged total pressure at the 30% axial chord downstream from the trailing edge. 13 design variables are chosen for blade shape modification. A 10.8 % reduction of total pressure loss on the turbine rotor is achieved by this process, which is same as a more than 1% total-to-total efficiency increase. The computed results are compared with those using 11 design variables, and show that optimized results depend heavily on the accuracy of blade design.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2929-2934
• /
• 2007

Prediction methods for cavitation noise are presented. At first, direct numerical simulation of cavitating flow noise has been performed, and acoustic analogy equation based on the cavitation noise modeling is derived. For the direct numerical simulation, a density based homogenous equilibrium model is employed to simulate cavitating two-phase flow and the governing equations are solved with high-order numerical schemes to resolve cavitation noise. The compressible Navier-Stokes equations for mixture fluids are discretized with a sixth-order central compact scheme, and the steep gradient of flow variables and supersonic regions are treated with the selective spatial filtering technique. The direct simulation of cavitating flow noise is performed for a 2D circular cylinder at cavitation number 0.7 and 1. The far-field noise is also predicted with the derived analogy equation. Noise spectrum predicted with the equation is well compared with the result of direct numerical simulation and also agree well with the theory.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.49-57
• /
• 2014

In the present study, 2-D gas-kinetic flow solver on unstructured meshes was developed for flows from continuum to transitional regimes. The gas-kinetic BGK scheme is based on numerical solutions of the BGK simplification of the Boltzmann transport equation. In the initial reconstruction, the unstructured version of the linear interpolation is applied to compute left and right states along a cell interface. In the gas evolution step, the numerical fluxes are computed from the evaluation of the time-dependent gas distribution function around a cell interface. Two-dimensional compressible flow calculations were performed to verify the accuracy and robustness of the current gas-kinetic approach. Gas-kinetic BGK scheme was successfully applied to two-dimensional steady and unsteady flow simulations with strong contact discontinuities. Exemplary hypersonic viscous simulations have been conducted to analyze the performances of the gas-kinetic scheme. The computed results show fair agreement with other standard particle-based approaches for both continuum part and transitional part.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.187-193
• /
• 2007

The flight vehicles have cavities such as wheel wells and bomb bays. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves. Resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. In the present study, numerical analysis was performed for cavity flows by the unsteady compressible three dimensional Reynolds-Averaged Navier-Stokes (RANS) equations with Wilcox's ${\kappa}\;-\;{\omega}$ turbulence model. The cavity has the aspect ratios of 2.5, 3.5 and 4.5 for two-dimensional case, same aspect ratios with the W/D ratio of 2 for three-dimensional case. The Mach and Reynolds numbers are 0.53 and 1,600,000 respectively. The flow field is observed to oscillate in the "shear layer mode" with a feedback mechanism. Based on the SPL(Sound Pressure Level) analysis of the pressure variation at the cavity trailing edge, the dominant frequency was analyzed and compared with the results of Rossiter's formula. The MPI(Message Passing Interface) parallelized code was used for calculations by PC-cluster.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.5
• /
• pp.735-745
• /
• 1997

In this study, the compressibility of resin was considered in filling analysis to account for the possible packing type flow. A numerical simulation program employing a hybrid finite element/finite difference scheme was developed to solve Hele-Shaw flow of the compressible viscous fluid at non-isothermal conditions. To advance the melt front, a control volume approach was adopted. Thin complex 3-D shapes of cavities, runners, and sprues were discretized by employing triangular, cylindrical and/or rectangular strip elements. Mass conservation was applied to each control volume to solve for the pressure distribution. Directly applying a constant mass flow rate at the inlet removes calculation of the apparent pressure boundary conditions, resulting in better simulation condition. The Cross model was used to model viscosity and the Tait equation was employed to represent density as a function of temperature and pressure. The validity of the developed program was verified through comparisons with available data in the literature and the effect of compressibility on the pressure distribution was discussed. To reduce computation time, 1-D and 2-D elements were used instead of applying triangular elements and the numerical results were compared to each other.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.4
• /
• pp.112-119
• /
• 1994

This study is investigated into the dynamic effect of the manifold configuration during the gas exchange processes using both simulation and experiment, In theoretical study on the flow analysis, the characteristic method is applied to solve the compressible unsteady flow equation, involving the several steady flow boundary conditions. In order to excute the engine experiment efficiently, a data acquisition system is configured by using A/D converter and PC. Good results which coincided experimental data with simulation output were obtained, and it shows that this simulation method can be applied to obtain the optimal design parameters in the intake and exhaust systems.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.168-173
• /
• 1998

A numerical analysis of axisymetric backward facing step main nozzle flow in air jet loom has been accomplished. To obtain basic design data for an optimum main nozzle for an air-jet loom and to predict the transonic/supersonic flow, a characteristic based upwind flux difference splitting compressible Navier-Stokes method has been used. The wall static pressure of the main nozzle and the flow velocity changes in the nozzle tube were analyzed by changing air tank pressures and acceleration tube lengths. The flow inside the nozzle experiences double choking one at the needle tip and the other at the acceleration tube exit at tank pressures over $4kg_f/cm^2$. The tank pressure $P_t$ leading to the critical condition depends on the acceleration tube length; i.e, $P_t$ is higher for longer acceleration tubes. The $P_t$ value required to bring the acceleration tube exit to the critical condition is nearly constant regardless of acceleration tube length. The round needle tip shape might lead to less total pressure loss when compared with step shape.

• Journal of Korean Society of Water and Wastewater
• /
• v.24 no.2
• /
• pp.231-236
• /
• 2010

A compressible media filtration process with synthetic fiber media was studied for combined sewer overflows (CSOs) treatment. Since the operation performance of fiber media filtration was dependent on the pattern of CSOs, the flow rate of CSOs was investigated and it was characterized by a big fluctuation. Thus, in this study, the fiber media filtration process was tested with wide range of filtration velocity. The removal efficiency was proportion to the increase in compressibility. As the filtration velocity was increased, the treatment efficiency was decreased and consequently leveled off when the velocity exceeded 750 $m^3/m^2$/d. An exponential equation was introduced to express the relationship between the removal efficiency and up-flow velocity. At columm test, six repetition of filtration and backwash cycle did not after the filtering velocity under the constant pressure condition.