• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.2
• /
• pp.54-63
• /
• 2009

Large eddy simulation was conducted for flow development in a chamber with wall injection which simulates the cold flow in an idealized hybrid rocket motor. It was found that a peculiar timescale, roughly corresponding to St~0.5, resides in the flowfield resulting from the interaction between the main oxidizer and wall injected flows. However, the fact that this time characteristics is absent in the temperature field in the vicinity of the wall indicates that even a small regression rate renders the passive scalar, such as temperature, dissimilar to the velocity field. This implies that a classical approach, which assumes that constant turbulent Prandtl number, should be replaced by a more sophisticated turbulence models to accurately predict the temperature field in the hybrid motor.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.1
• /
• pp.64-72
• /
• 2010

Due to the interaction between main oxidizer flow and the wall injected flow resulting from the regression process, a specific time characteristics identified in the frequency spectrum of streamwise velocity is generated in the hybrid rocket motor. In order to understand the response of the turbulent flow to two different types of external momentum forcing, LES analysis was conducted without considering the combustion. It turns out that both concentrated and distributed forcings do not lead to the disastrous resonance phenomenon. Energy contents are enhanced due to the added momentum but the peak frequency was not modified in the turbulent flow near the end of the rocket motor. Natural frequency of the flow system should be taken into account to further pursue the instability issue by using external forcing.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.8
• /
• pp.839-847
• /
• 2012

Numerical simulations for channel flows with $Re_{\tau}$ = 180, 395 and 590 have been performed to investigate the hairpin packet formation process in wall-bounded turbulent flows. Using direct numerical simulation databases, the initial flow fields are given by the conditionally averaged flow field with the second quadrant event specified at the buffer layer. By tracking the initial vortex development, the change in the initial vortex to an ${\Omega}$-shaped vortex and th generation of a secondary hairpin vortex were found to occur with time scales based on the wall units. In addition, at the time when the initial vortex has grown to the channel center, the inclination angle of the hairpin vortex packet is approximately $12{\sim}14^{\circ}$, which is similar for all three Reynolds numbers. Finally, numerical simulations of the evolution of two adjacent hairpin vortices have been performed to examine the interaction between the adjacent vortex packets.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.3
• /
• pp.345-352
• /
• 1985

본 연구에서는 내류(internal flow)에서 유로가 비대칭으로 급확대될 경우의 박리현상과 유동현상을 고찰하였다. 비대칭 급확대 채널에서의 층류영역에서 난류영 역까지의 유동현상을 B.F. Armaly, C.E. Thomas는 실험적 해석과 유한요소법을 사용하 여 이론적 해석을 하였고 Donald. M. Kuehn, Denham & Patrick, Kwon, Patrick J. Ro- ache, Anand Kumar등은 같은 모델에 대해서 실험적 해석과 유한차분법을 사용하여 이 론적 해석을 하였으며 지금까지의 유한해석법에 의한 연구는 입구와 출구조건이 같은 경우 및 밀폐 공간 혹은 한면의 속도가 주어지는 밀폐공간등에 대해 수행되어 왔으나 본 연구에서는 입구와 출구조건이 같지 않은 2차원, 비대칭 급확대 채널에서의 유동현 상을 유한해석법으로 해석하여 실험치와 비교하였다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.9
• /
• pp.1302-1310
• /
• 2002

A direct numerical simulation in turbulent curved channel flow is performed. The drifting Taylor-Gortler vortices are identified by applying a conditional averaging. A new algorithm is proposed based on the wavelet transform of the wall information. A continuous wavelet transform with Marr wavelets is employed to decompose the flow signals at a chosen length scale. An active cancellation is applied to attenuate the Taylor-Gortler vortices and to reduce the wall skin friction.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.118-124
• /
• 2008

The characteristics of heat transfer in a two-dimensional channel obstructed by a square rod is investigated by a turbulence model. The computation is made for the six cases of different rod positions between channel walls. To analyze the wall heat transfer, the heat flux of channel walls is set as a constant value and the $k-{\epsilon}-f_{\mu}$ model is employed. Downstream the square rod, the flow recirculation region appear and they are varied by the rod position. The enhancement of the turbulent heat transfer to the wall is induced by the flow instability using a square rod. The averaged heat transfer rate is maximized at a specific rod position. Finally, the effects of square rod on unsteady flows are scrutinized with the frequency analysis.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.859-866
• /
• 2003

The Lagrangian dispserion of fluid particles in inhomogeneous turbulence is investigated by a direct numerical simulation of turbulent channel flow. Fluid particle velocity and acceleration along a particle trajectory are computed by employing several interpolation schemes such as linear interpolation, high-order Lagrange polynomial interpolation and the Hermite interpolation schemes. The performances of the schemes are evaluated through comparison of errors in computed particle positions, velocities and accelerations against spectral interpolation. Adopting the four-point Hermite interpolation in the homogeneous directions and Chebyshev polynomials in the wall-normal direction appears to produce most reliable Lagrangian statistics including acceleration correlations with a reasonable amount of computational overhead.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.867-876
• /
• 2003

The Lagrangian dispersion of fluid particles in inhomogeneous turbulence is investigated by a direct numerical simulation of turbulent channel flow. Four points Hermite interpolation in the homogeneous direction and Chebyshev polynomials in the inhomogeneous direction is adopted to simulate the fluid particle dispersion. An inhomogeneity of Lagrangian statistics in turbulent boundary layer is investigated by releasing many particles at several different wall-normal locations and tracking those particles. The fluid particle dispersions and Lagrangian structure functions of velocity are scaled by the Kolmogorov similarity. The auto-correlations of velocity and acceleration are shown at the different releasing locations. Effect of initial particle location on the dispersion is analyzed by the probability density function at the several downstreams and time instants.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.8
• /
• pp.1060-1067
• /
• 2001

A large-scale bursting event has been analyzed in a turbulent channel flow using a data obtained from a direct numerical simulation (DNS). Large-scale, plume-like structures have been frequently observed in many experimental results, but the origin of those structures is far from being fully understood. It is believed that those large scale events occur occasionally but contribute significantly to the generation of Reynolds shear stress in the outer layer. This paper attempts to give detailed examples of those large-scale motions observed in a turbulent channel flow at relatively low Reynolds number.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.9
• /
• pp.1227-1236
• /
• 2001

A systematic analysis is made of suboptimal control for drag reduction. The influence of the amplitude of actuation (A) and the time scale of actuation ($\Delta$t(sub)a(sup)+) is evaluated. Two wall sensing variables are employed (∂w/∂y│(sub)w and ∂p/∂z│(sub)w) with two wall actuations (${\Phi}$$_2$and ${\Phi}$$_3$). To test the suboptimal control, direct numerical simulations of turbulent channel flow at Re(sub)$\tau$=100 are performed in a spectral domain. It is found that the effect of A and $\Delta$t(sub)a(sup)+∼1. The near-wall behaviors of flow structure are analyzed to characterize the drag reduction. The size effect of the sensor/actuator is examined.