• Journal of the Korean Society for Marine Environment & Energy
• /
• v.20 no.1
• /
• pp.37-44
• /
• 2017

Complexity of multiphase flows due to existence of more than two interface including free-surface in one system, cannot be simulated easily. Since more than two fluids affect to flows and disturb interface, non-linearities such as instabilities can be appeared. Among several instabilities on multiphase flows, one of representative is Rayleigh-taylor instability. In order to examine in importance of density disparity, several cases with numerous Atwood number are set. Moreover, investigation of influence on initial disturbance were also considered. Moving particle simulation (MPS) method, which was employed in this paper, was not widely used for multiphase problem. In this study, by adding new particle interaction models such as self-buoyance correction, surface tension, and boundary condition at interface models, MPS were developed having more strength of physics and robust. By applying newly developed multiphase MPS, considered cases are performed and compared each other. Additionally, though existence of disagreement of magnitude of rising velocity between theoretical values from linear potential theory and that of numerical simulation, agreement of tendency can be proved of similarity of result. the discordance of magnitude can be explained due to non-linear effects on numerical simulation which was not considered in theoretical result.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.11
• /
• pp.897-905
• /
• 2008

Selective catalytic reduction(SCR) is known as one of promising methods for reducing $NO_x$ emissions in diesel exhaust gases. $NO_x$ emissions react with ammonia in the catalyst surface of SCR system at working temperature of catalyst. In this study, to raise the reacting temperature when the exhaust gas temperature is too low, a heater is located at the bottom of SCR reactor. At an ambient temperature, ammonia is radially injected perpendicular to the exhaust gas flow at inlet pipe and uniformly mixed in the mixing area after being impinged against the wall. To predict the turbulent model inside the mixing area of SCR system, the standard ${\kappa}\;-\;{\varepsilon}$ model is applied. This work investigates numerically the effects of induced heat on the gaseous flow. The results show that the Taylor-$G{\ddot{o}}rtler$ type vortex is generated after the gaseous flow impinges the wall in which these vortices influence the temperature distribution. The addition of heat disturbs the flow structure in bottom area and then stretching flow occurs. Vorticity strand is also formed when heat is continuously increased. Constriction process takes place, however, when a further heat input over a critical temperature is increased and finally forms shed vortex which is disconnected from the vorticity strand. The strong vortex restricts the heat transport in the gaseous flow.

• Korean Chemical Engineering Research
• /
• v.53 no.5
• /
• pp.609-613
• /
• 2015

The onset of instability induced by impulsive spin-down of the rigid-body flow placed in the gap between two coaxial cylinders is analyzed by using the energy method. In the present stability analysis the growth rate of the kinetic energy of the base state and also that of disturbances are taken into consideration. In the present system the primary flow is a transient, laminar one. But for the Reynolds number equal or larger than a certain one, i.e. $Re{\geq}Re_G$ secondary motion sets in, starting at a certain time. For $Re{\geq}Re_G$ the dimensionless critical time to mark the onset of vortex instabilities, ${\tau}_c$, is here presented as a function of the Reynolds number Re and the radius ratio ${\eta}$. For the wide gap case of small ${\eta}$, the transient instability is possible in the range of $Re_G{\leq}Re{\leq}Re_S$. It is found that the predicted ${\tau}_c$-value is much smaller than experimental detection time of first observable secondary motion. It seems evident that small disturbances initiated at ${\tau}_c$ require some growth period until they are detected experimentally.

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

An experimental study is made of turbulent shear flows in a nearly two-dimensional 90.deg. curved duct by using the hot-wire anemometer. The Reynolds normal and shear stresses, triple velocity products, integral length scales, Taylor micro length scales and dissipation length scales are measured and analyzed. For a positive shear at the inlet, the afore-mentioned turbulence quantities are all suppressed. However, when the inlet shear flow is negative, they are augmented, i.e., the convex curvature suppresses the turbulence whereas the concave curvature augments it. It is found that the curvature effects are rather sensitive to the triple velocity products than the Reynolds stresses. The evolution of turbulence under the curvature with the different shear conditions is well described by the modified curvature parameter S' and the non-dimensional development time ${\tau}$.'

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.568-571
• /
• 2008

The computation of moving interface by the level set method typically requires reinitializations of level set function. An inaccurate estimation of level set function ${\phi}$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, accurate and robust reinitialization process is essential to the free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates directly level set function ${\phi}$ using a normal vector in the interface without solving the re-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1splitting FEM are adopted to discretize advection equation of the level set function and the Navier-Stokes equation, respectively. Advection equation of free surface and re-initialization process are validated with benchmark problems, i.e., a broken dam flow and time-reversed single vortex flow. The simulation results are in good agreement with the existing results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.10
• /
• pp.754-760
• /
• 2008

Computation of moving interface by the level set method typically requires the reinitialization of level set function. An inaccurate estimation of level set function $\phi$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, an accurate and robust reinitialization process is essential to the simulation of free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates level set function directly using a normal vector on the interface without solving there-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1 splitting/SUPG (Streamline Upwind Petrov-Galerkin) FEM are adopted to discretize advection equation of the level set function and the incompressible Navier-Stokes equation, respectively. Advection equation and re-initialization process of free surface capturing are validated with benchmark problems, i.e., a broken dam flow and timereversed single vortex flow. The simulation results are in good agreement with the existing results.

• Journal of ILASS-Korea
• /
• v.17 no.1
• /
• pp.35-44
• /
• 2012

To improve the mixing and atomizing performance at the center region of the conventional coaxial shear injector spray, the concept of a coaxial porous injector was invented. This novel injection concept for liquid rocket engines utilizes the Taylor-Culick flow in the cylindrical porous tube. The 2-dimensional injector, which can be converted in three injection configurations, was fabricated, and several cold flow tests using water-air simulant propellant was performed. The hydraulic characteristics and the effects of a gas flow condition on the spray pattern and the Sauter mean diameter (SMD) was analyzed for each configuration. The atomizing mechanism of coaxial porous injector was different with the coaxial shear injector, and it was explained by the momentum of the gas jet, which is injected normally against the center liquid column, and by the secondary disintegration at the wavy interface of liquid jet, which was generated at the recessed region. The SMD of 2D coaxial porous injector, which has higher gas momentum, was measured and it shows better atomizing performance at the center and outer side of spray than the 2D coaxial shear injector.

• Tunnel and Underground Space
• /
• v.12 no.4
• /
• pp.277-283
• /
• 2002

This experimental study concerns the characteristics of a helical flow in a concentric annulus with a diameter ratio of 0.52, whose outer cylinder is stationary and inner one is rotating. The pressure losses and skin friction coefficients have been measured for the fully developed flow of Non-Newtonian fluid, aqueous solution of sodium carbomethyl cellulose (CMC) and bentonite with inner cylinder rotational speed of 0~400 prm. Also, the visualization of helical flows has been performed to observe the unstable waves. The results of present study reveal the relation of the Reynolds number Re and Rossby number Ro with respect to the skin friction coefficients. In somehow, they show the existence of flow instability mechanism. The pressure losses increase as the rotational speed increases, but the gradient of pressure losses decreases as the Reynolds number increases in the regime of transition and turbulence. And the increase of flow disturbance by Taylor vortex in a concentric annulus with rotating inner cylinder results in the decrease of the critical Reynolds number with the increase of skin friction coefficient.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.4
• /
• pp.227-236
• /
• 2016

In the present experimental study, the effect of gas- and liquid-injected methods on the formation of bubble and liquid slug at the merging micro T-junction of a square microchannel with dimensions $600{\mu}m{\times}600{\mu}m$ was investigated. Nitrogen and water were used as test fluids. The superficial velocities of the liquid and gas were in the range of 0.05 - 1 m/s, and 0.1 - 1 m/s, respectively, where the Taylor flow was observed. The bubble length, liquid slug length, bubble velocity, and bubble generation frequency were measured by analyzing the images captured using a high-speed camera. Under similar inlet superficial velocity conditions, in the case of gas injection to the main channel at the merging T-junction (T_gas-liquid), the lengths of the bubble and liquid slug were longer, and the bubble generation frequency was lower than in the case of liquid injection to the main channel at the merging T-junction (T_liquid-gas). On the other hand, in both cases, the bubble velocity was almost the same. The previous correlation proposed using experimental data for T_liquid-gas had predicted the present experimental data of bubble length, bubble velocity, liquid slug length, and bubble generation frequency for T_gas-liquid to be ~24%, ~9%, ~39%, ~55%, respectively.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.12
• /
• pp.48-58
• /
• 2006

The spray and combustion characteristics of liquid jet in cross flow with variation of injection angle are numerically studied. Numerical analysis was carried out using KIVA code, which may be used to generate numerical solutions to spray and chemical reactive fluid problem in three space dimensions and modified to be suitable for simulating liquid jet ejected into the cross flow. Wave model and Kelvin- Helmholtz(KH) /Rayleigh-Taylor(RT) hybrid model were used for the purpose of analyzing liquid column, ligament, and the breakup of droplet. Penetration length increases as flow velocity decreases and injection velocity increases. Numerical error increases as inflow velocity increases. The results of flame propagation contour in combustion chamber and local temperature distribution, combustion emissions were obtained.