• Journal of Power System Engineering
• /
• v.3 no.4
• /
• pp.16-21
• /
• 1999

Numerical simulation on two-dimensional turbine cascade flow has been performed using compressible Navier-Stokes equations. The flow equations are written in a cartesian coordinate system, then mapped into a generalized body-fitted ones. All direction of viscous terms are incoporated and turbulent effects are modeled using the extended ${\kappa}-{\epsilon}$ model. Equations are discretized using control volume SIMPLE algorithm on the nonstaggered grid sysetem. Applications are made at a VKI turbine cascade flow in atransonic wind-tunnel and compared to experimental data. Present numerical results are shown to be in good agreement with the experimental results and simulate the compressible viscous flow characteristics inside the turbine blade passage.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.184-190
• /
• 2000

A well-known pressure correction method, a SIMPLE algorithm is extended to treat compressible flows. Collocated grids are used and density is linked to pressure via an equation of state. The influence of pressure on density in the case of compressible flows is implicitly incorporated into the extended SIMPLE algorithm. The first-order Upwind and high-order Quick scheme are compared with respect to an accuracy and convergence time at all speeds. The extended method is verified on a number of test cases and the results we compared with other numerical results available in the literature. The calculated results show that the Quick scheme improves accuracy at all speed and also reduces the calculation time at supersonic flows, compared with the Upwind scheme.

• Journal of computational fluids engineering
• /
• v.5 no.2
• /
• pp.1-8
• /
• 2000

A well-known pressure correction method, a SIMPLE algorithm, is extended to treat compressible flows. Collocated grids are used and density is linked to pressure via an equation of state. The influence of pressure on density in the case of compressible flows is implicitly incorporated into the extended SIMPLE algorithm. The first-order Upwind and high-order Quick scheme are compared with respect to an accuracy and convergence time at all speeds. The extended method is verified on a number of test cases and the results are compared with other numerical results available in the literature. The calculated results show that the Quick scheme improves accuracy at all speed and also reduces the calculation time at supersonic flows, compared with the Upwind scheme.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.329-332
• /
• 2002

A theoretical study is made of the steady flow of a compressible fluid in a rapidly rotating finite cylinder. Flow is generated by imposing mechanical and/or thermal disturbances at the rotating endwall disks. Both the Ekman and Rossby numbers are small. A detailed consideration is given to the energy budget for a control volume in the Ekman boundary layer. A combination of physical variables, which is termed the energy contents, consisting of temperature and modified angular momentum, emerges to be relevant. The distinguishing features of a compressible fluid, in contrast to those of an incompressible fluid, are noted. For the Taylor-Proudman column to be sustained, in the interior, it is shown that the net energy transport between the solid disk wall and the interior fluid should vanish. Physical rationalizations are facilitated by resorting to the concept of the afore-stated energy content.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.52-61
• /
• 2014

The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

• Journal of computational fluids engineering
• /
• v.17 no.3
• /
• pp.99-107
• /
• 2012

In this paper, the shock waves in compressible solids and liquids are simulated using a six-equation diffuse interface multiphase flow model that is extended to the Cochran and Chan equation of state. A pressure relaxation method based on a volume fraction function and a pressure-correction equation are newly implemented to the six-equation model. The developed code has been validated by a shock tube problem with liquid nitromethane and an impact problem of a copper plate on a solid explosive. In addition, a new problem, an impact of a copper plate on liquid nitromethane, has been solved. The present code well shows the wave structures in compressible solids and liquids without any numerical oscillations and overshoots. After the impact of a solid copper plate on liquid, two shock waves (one propagates into liquid and the other into solid) are generated and a material interface moves to the impacting direction. The computational results show that the shock velocity inside the liquid linearly increases with the impact velocity.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.1
• /
• pp.132-141
• /
• 1995

Characteristics of pressure-drop oscillations(PDO) in a boiling channel were studied numerically and compared with experimental data. Effects of initial and boundary conditions on PDO were investigated in terms of oscillation period and amplitude. The period and amplitude of PDO increased with increasing of the compressible volume in the surge tank and the heat input. PDO occurred within the specific range of the fluid temperature, at which oscillation period and amplitude diminished rapidly with the increase of the fluid temperature. The increase of the loss coefficient in fluid supply line resulted in slightly longer oscillation period and larger amplitude. Numerical results showed good agreement with the experimental data.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.305-310
• /
• 2001

This paper presents the computational method for analyzing the compressible flow fields in a high voltage gas circuit breaker. There are many difficult problems in analyzing the gas flow in GCB due to complex geometry, moving boundary, shock wave and so on. In particular, the distortion problem of the grid due to the movement of moving parts can be worked out by the fixed grid technique. Numerical simulations are based on a fully implicit finite volume method of the compressible Reynolds-averaged Navier-Stokes equations to obtain the pressure, density, and velocity through the entire interruption process. The presented method is applied to the real circuit breaker model and the pressure in front of the piston is good agreement with the experimental one.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.1
• /
• pp.77-87
• /
• 2003

The residual stress and birefringence in injection-molded plastic parts can be divided into the flow-induced residual stress and birefringence produced in flowing stage, the thermally-induced residual stress and birefringence produced in cooling stage. However, the physics involved in the generation of the thermally-induced residual stress and birefringence still remains to be understood. Because polymer experiences viscoelastic history near the glass-transition temperature it is hard to model the entire process. Volume relaxation phenomenon was included to predict the final thermally-induced residual stress and birefringence in quenched plastic parts more accurately. The present study focused on comparing the predicted values far thermally-induced residual stress and birefringence with and without volume relaxation behavior (or free volume theory) under free and constrained quenching conditions. As a result, tile residual stress remained as a tensile stress at the center and as a compressible stress near the surface for the free quenching cases. In contract the residual stress remained as a compressible stress at the center and as a tensile stress near the surface fur the constrained quenching cases. The residual birefringence remained as minus values at the center and as plus values near the surface for the free quenching cases. Interestingly the residual birefringence showed minus values in entire zone for the constrained quenching cases. In the prediction of birefringence only the case including free volume theory showed the correct result for the distribution of birefringence in thickness direction.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3003-3008
• /
• 2007

Critical nozzle has been frequently employed to measure the flow rate of various gases, but hydrogen gas, especially being at high-pressure condition, was not nearly dealt with the critical nozzle due to treatment danger. According to a few experimental data obtained recently, it was reported that the discharge coefficient of hydrogen gas through the critical nozzle exceeds unity in a specific range of Reynolds number. No detailed explanation on such an unreasonable value was made, but it was vaguely inferred as real gas effects. For the purpose of practical use of high-pressure hydrogen gas, systematic research is required to clarify the critical nozzle flow of high-pressure hydrogen gas. In the present study, a computational fluid dynamics(CFD) method has been applied to predict the critical nozzle flow of high-pressure hydrogen gas. Redlich-Kwong equation of state that take account for the forces and volume of molecules of hydrogen gas were incorporated into the axisymmetric, compressible Navier-Stokes equations. A fully implicit finite volume scheme was used to numerically solve the governing equations. The computational results were validated with some experimental data available. The results show that the coefficient of discharge coefficient is mainly influenced by the compressibility factor and the specific heat ratio, which appear more remarkable as the inlet total pressure of hydrogen gas increases.