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

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.8
• /
• pp.1128-1138
• /
• 2000

Pumping action in ejector systems is generally achieved through the mixing of a high-velocity and high-energy stream with a lower-velocity and lower-energy stream within a duct. The design and performance evaluation of the ejector systems has developed as a combination of scale-model experiments, empiricism and theoretical analyses applicable only to very simplified configurations, because of the generic complexity of the flow phenomena. In order to predict the detailed performance characteristics of such systems, the flow phenomena throughout the operating regimes of the ejector system should be fully understood. This paper presents the computational results for the two-dimensional supersonic ejector system with a second throat. The numerical simulations are based on a fully implicit finite volume scheme of the compressible Reynolds-averaged Navier-Stokes equation in a domain that extends from the stagnation chamber to the diffuser exit. For a wide range of the operating pressure ratio the flow field inside the ejector system is investigated in detail. The results show that the supersonic ejector systems have an optimal throat area for the operating pressure ratio to be minimized.

• Journal of Advanced Marine Engineering and Technology
• /
• v.26 no.5
• /
• pp.596-606
• /
• 2002

The three-dimensional unsteady compressible Euler equation solver with ALE, CFD code, PAM-FLOW based on FEM method has been applied to analyze the flow field around the high speed train which is entering into a channel. From the present study, the pressure and flow transients were calculated and analyzed. The generation of compression wave was observed ahead of train and the high pressure in the gap between the train and the tunnel was also found due to the blockage effects. It was found that abrupt fluctuation in pressure exists in the region from train nose to shoulder of train corresponding to 10％ of total length of train during tunnel entry. Computed time history of aerodynamic forces of train during tunnel entry show that drag coefficient rapidly rises and saturates at about non-dimensional time 0.31. The total increase of drag coefficient before and after tunnel entry is about 1.1%. Transient profile of lift force shows similar pattern to drag coefficient except abrupt drop after saturation and lift force in the tunnel increases 0.08% more than that before tunnel entry.

• Journal of Mechanical Science and Technology
• /
• v.18 no.3
• /
• pp.407-418
• /
• 2004

Prandtl-Meyer expansion flow with homogeneous condensation is investigated experimentally and by numerical computations. The steady and unsteady periodic behaviors of the diabatic shock wave due to the latent heat released by condensation are considered with a view of technical application to the condensing flow through steam turbine blade passages. A passive control method using a porous wall and cavity underneath is applied to control the diabatic shock wave. Two-dimensional, compressible Navier-Stokes with the nucleation rate equation are numerically solved using a third-order TVD (Total Variation Diminishing) finite difference scheme. The computational results reproduce the measured static pressure distributions in passive and no passive Prandtl-Meyer expansion flows with condensation. From both the experimental and computational results, it is found that the magnitude of steady diabatic shock wave can be considerably reduced by the present passive control method. For no passive control, it is found that the diabatic shock wave due to the heat released by condensation oscillates periodically with a frequency of 2.40㎑. This unsteady periodic motion of the diabatic shock wave can be completely suppressed using the present passive control method.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.363-366
• /
• 2010

Among the many different types of wind tunnels, Ludwieg Tube(LT) is the most suitable facility for high Reynolds number testing. Depending on the location of diaphragm, there are two types of LTs. In the present study, a computational work has been carried out to compare the operation characteristics of upstream- and downstream-diaphragm LTs. Two-dimensional, axisymmetric, unsteady, compressible Navier-Stokes equations were solved using a fully implicit finite volume scheme. Based on the present results, the flow mechanism of the starting process was discussed in detail using wave diagrams and characteristics of starting time and working time were investigated.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.10
• /
• pp.649-657
• /
• 2017

A rupture in a high-pressure pipe causes the fluid in the pipe to be discharged in the atmosphere at a high speed resulting in a supersonic jet that generates the compressible flow. This supersonic jet may display complicated and unsteady behavior in general. In this study, Computational Fluid Dynamics (CFD) analysis was performed to investigate the compressible flow generated by a supersonic jet ejected from a high-pressure pipe. A Shear Stress Transport (SST) turbulence model was selected to analyze the unsteady nature of the flow, which depends upon the various gases as well as the diameter of the pipe. In the CFD analysis, the basic boundary conditions were assumed to be as follows: pipe of diameter 10 cm, jet pressure ratio of 5, and an inlet gas temperature of 300 K. During the analysis, the behavior of the shockwave generated by a supersonic jet was observed and it was found that the blast wave was generated indirectly. The pressure wave characteristics of hydrogen gas, which possesses the smallest molecular mass, showed the shortest distance to the safety zone. There were no significant difference observed for nitrogen gas, air, and oxygen gas, which have similar molecular mass. In addition, an increase in the diameter of the pipe resulted in the ejected impact caused by the increased flow rate to become larger and the zone of jet influence to extend further.

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.80.2-80.2
• /
• 2005

• Proceedings of the KSME Conference
• /
• 2005.11a
• /
• pp.81.1-81.1
• /
• 2005