• Journal of the Korean Society of Propulsion Engineers
• /
• v.4 no.4
• /
• pp.59-69
• /
• 2000

A numerical analysis based on two-dimensional, incompressible and compressible Navier-Stokes equations was carried out for double circular arc compressor cascade and the results are compared with available experimental data. The incompressible code based on SIMPLE algorithm adopts pressure weighted method and hybrid scheme for the convective terms. The compressible code with preconditioning method involves a upwind-biased scheme for the convective terms and LU-SGS scheme for temporal integration. Several turbulence models are evaluated by computing the turbulent viscous flows; Baldwin-Lomax, standard $\kappa$ -$\varepsilon$, $\kappa$ -$\varepsilon$ Lam. Bremhorst, standard $\kappa$-$\omega$, $\kappa$ -$\omega$ SST model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.9
• /
• pp.1106-1112
• /
• 1999

Substantial losses behind axial flow rotor are generated by the wake, various vortices in the hub region and the tip leakage vortex in the tip region. Particularly, the leakage vortex formed near blade tip is one of the main causes of the reduction of performance, generation of noise and aerodynamic vibration in downstream. In this study, the three-dimensional flow fields in an axial flow rotor were calculated with varying tip clearance under various flow rates, and the numerical results were compared with experimental ones. The numerical technique was based on SIMPLE algorithm using standard $k-{\varepsilon}$ model(WFM) and Launder & Sharma's Low Reynolds Number $k-{\varepsilon}$ model(LRN). Through calculations, the effects of tip clearance and attack angle on the 3-dimensional flow fileds behind a rotor and leakage flow/vortex were investigated. The presence of tip leakage vortex, loci of vortex center and its behavior behind the rotor for various tip clearances and attack angles was described well by calculation.

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

Nonlinear relationship between Reynolds stresses and the rate of strain of nonlinear k-$\varepsilon$models is evaluated theoretically by using the boundary layer assumptions against the turbulence-driven secondary flows in noncircular ducts and then their prediction performance is validated numerically through the application to the fully developed turbulent flow in a square duct. Typical predicted quantities such as mean axial and secondary velocities, turbulent kinetic energy and Reynolds stresses are compared with available experimental data. The nonlinear k-$\varepsilon$ model adopted in a commercial code is found to be unable to predict accurately duct flows with the prediction level of secondary flows one order less than that of the experiment.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.224-232
• /
• 2000

Four turbulent $k-{\varepsilon}$ models (i.e., standard model, modified models with streamline curvature modification and/or preferential dissipation modification) are applied in order to analyze the turbulent flow of wall-attaching offset jet. For numerical convergence, this paper develops a method of slowly increasing the convective effect induced by skew-velocity in skew-upwind scheme (hereafter called Partial Skewupwind Scheme). Even though the method was simple, it was efficient in view of convergent speed, computer memory storage, programming, etc. The numerical results of all models show good prediction in first order calculations (i.e., reattachment length, mean velocity, pressure), while they show some deviations in ·second order (i.e., kinetic energy and its dissipation rate). Like the previous results obtained by upwind scheme, the streamline curvature modification results in better prediction, while the preferential dissipation modification does not.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.741-750
• /
• 1995

The turbulent buoyancy-driven flow in 2-dimensional enclosed cavities heated from the vertical side is numerically calculated for both cases of a Rayleigh number of 5*10$^{10}$ for air and 2.5*10$^{10}$ for water. Three different turbulence models are considered : standard k-.epsilon. model of Ozoe and low-Reynolds-number model of Lam and Bremhorst, and another low-Reynolds-number model of Davidson. The results indicate that the use of low-Reynolds number models is recommended for the indoor airflow computation, and the results from Davidson model are reasonably close to the reported experimental data. A sensitivity study shows that the amounts of wall-heat transfer and the velocity profiles with the Lam and Bremhorst model largely depend on the choice of the wall function for .epsilon..

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.8
• /
• pp.1046-1055
• /
• 1997

An anisotropic k-.epsilon. turbulence model for predicting the rotating flows is proposed with the simple inclusion of a new parameter dealing with the extra straining effects in the .epsilon.-equation. This model is employed to compute the effects of Coriolis forces on fully-developed flow in a rotating channel. The predicted results indicate that the present model captures fairly well the striking rotational-induced effects on the Reynolds stresses and the mean flow distributions, including the argumentation of turbulent transport on the unstable side (pressure surface) of the channel and its damping on the stable side (suction surface).

• Journal of the Korean Society of Propulsion Engineers
• /
• v.11 no.5
• /
• pp.37-50
• /
• 2007

This study investigated the characteristics of spray generated by a liquid coaxial swirl injector used in a combustor of the liquid rocket engine. The linear stability analysis considered long and short wave was introduced in liquid sheet breakup. Through the hydrodynamic analysis the initial liquid sheet thickness spray angle and injection velocity were predicted. To evaluate the effect of turbulence model standard $k-{\varepsilon}$ and RNC $k-{\varepsilon}$ model were applied to numerical calculation and it was known that RNC $k-{\varepsilon}$ model was more applicable to predict spray characteristics. On the basis of this evaluation validation of the developed model was performed with swirl injector installed in LPRE and the predicted results of breakup length, spray angle, and SMD agreed well with experiments qualitatively and quantitatively.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1149-1164
• /
• 1997

In this study, the applicability of the RNG k-.epsilon. model to the analysis of the complex flows is studied. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretized by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-.epsilon. model of three complex flows, i.e., the flow over a backward-facing step and a blunt flat plate, the flow around a 2D model car are compared to these from the standard k-.epsilon. model and experimental data. That of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly and the spray characteristics within a chamber of direct injection model engine are compared to these from the standard k-.epsilon. model and experimental data. The results of reattachment length, separated eddy size, average surface pressure distribution using the RNG k-.epsilon. model show more reasonable trends comparing with the experimental data than those using the modified k-.epsilon. model. Although the predicted rms velocity using the modified k-.epsilon. model is lower considerably than the experimental data in incylinder flow with poppet valve, predicted axial and radial velocity distributions at the valve exit and in-cylinder region show good agreements with the experimental data. The spray tip penetration predicted using the RNG k-.epsilon. model is more close to the experimental data than that using the modified k-.epsilon. model. The application of the RNG k-.epsilon. model seems to have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly and the spray characteristics over the modified k-.epsilon. model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.7
• /
• pp.911-918
• /
• 1997

The turbulent viscous wake flows behind a single airfoil, two-dimensional stationary blade row and three-dimensional rotating blade row were calculated, and the numerical results were compared with experimental ones. The numerical technique was based on the SIMPLE algorithm using three turbulent closure models, standard k-.epsilon. model(WFM), low Reynolds number k-.epsilon. model(LRN) and Reynolds stress model (RSM). In the case of a single airfoil, WFM, LRN and RSM presented fairly good velocity distributions in the wake compared with experimental data. In the case of the stationary blade row, LRN and RSM presented better results than WFM for wake velocity distribution, and especially LRN showed best results among these three turbulent models. In the case of the rotating blade row, WFM and LRN showed fairly good agreement with experimental data of the three-dimensional velocity component distributions in the range from hub to mid span region. LRN was also superior to WFM in accuracy of prediction for the wake velocity distribution as same with the cases of a airfoil and the stationary blade row.

• Polymer(Korea)
• /
• v.35 no.4
• /
• pp.325-331
• /
• 2011

Molecular simulations via the molecular dynamics method have been carried out to an equation of state of penetrable-sphere model fluids over a wide range of packing fraction ${\phi}$ and finite repulsive energy ${\varepsilon}^*$. The resulting simulation data are compared to theoretical predictions from the two limiting cases of high- and low-penetrability approximations available in the literature. A good agreement between theoretical and simulation results is observed ill the case of ${\varepsilon}^*$ <3.0. However, for the highly repulsive energy systems of ${\varepsilon}^*{\geqq}3.0$, where the potential energy barrier is more than two times higher than the particle kinetic energy, a poor agreement is found due to the clustering formation and the non-continuum size effects in the dense systems of ${\phi}{\geqq}0.7$ and ${\varepsilon}^*$=6.0.