• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.54-59
• /
• 2010

A parametric study has been carried out to elucidate the characteristics of channel flow with a streamwise-periodic array of cylinders. This flow configuration is relevant to heat exchanger applications. The presence of cylinders in channel flow causes the attached wall boundary layer to separate, leading to significant change in flow instabilities. There exist two kinds of instabilities; flow undergoes a primary instability (Hopf bifurcation) at a lower Reynolds number, and the unsteady two-dimensional flow becomes unstable to three-dimensional disturbances at a higher Reynolds number. We report here the dependencies of the primary instability as well as the flow characteristics of the subsequent unsteady flow, including flow-induced forces and Strouhal number of vortex shedding, on the distance between the cylinder and the channel wall.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.8 s.239
• /
• pp.963-970
• /
• 2005

The combustion instability of turbulent flames is the most important problem of the gas turbine combustor. Thus improved understanding of mechanisms of combustion instability is necessary for the design and operation of gas turbine combustors. In this study, the cause of the combustion instability in a rearward-step dump combustor was investigated with respect to the fuel flow modulation; choked fuel flow, unchoked fuel flow and fully premixed mixture flow. We observed various types of combustion instabilities with respect to the change of equivalence ratio, fuel flow conditions and fuel injection location. Particularly in the unchoked fuel flow condition, it was found that the oscillation time of combustion instability is strongly related to the convection time of the fuel and that the pressure fluctuation in a lab-scale combustor is highly related to the vortex and the equivalence ratio fluctuations due to fuel flow modulation and unmixedness of the fuel and air.

• Kyungpook Mathematical Journal
• /
• v.61 no.4
• /
• pp.831-843
• /
• 2021

In this paper we look at the three dimensional stagnation point flow problem over a rough rotating disk. We study the theoretical behaviour of the stagnation point flow, or forced flow, in the presence of a slip factor in which convective instability stationary modes appear. We make a numerical investigation of the effects of slip on the behaviour of the flow components of the stagnation point flow where the disk is rough. We provide, for the first time in the literature, a complete convective instability analysis and an energy analysis. Suitable similarity transformations are used to reduce the Navier-Stokes equations and the continuity equation into a system of highly non-linear coupled ordinary differential equations, and these are solved numerically subject to suitable boundary conditions using the bvp4c function of MATLAB. The convective instability analysis and the energy analysis are performed using the Chebyshev spectral method in order to obtain the neutral curves and the energy bars. We observe that the roughness of the disk has a destabilising effect on both Type-I and Type-II instability modes. The results obtained will be prominently treated as benchmarks for our future studies on stagnation flow.

• Korea-Australia Rheology Journal
• /
• v.20 no.1
• /
• pp.15-26
• /
• 2008

It has been known that there is a viscoelastic instability in the channel flow past a cylinder at high Deborah (De) number. Some of our numerical simulations and a boundary layer analysis indicated that this instability is related to the shear flow in the gap between the cylinder and the channel walls in our previous work. The critical condition for instability initiation may be related to an inflection velocity profile generated by the normal stress near the cylinder surface. At high De, the elastic normal stress coupling with the streamline curvature is responsible for the shear instability, which has been recognized by the community. In this study, an instability criterion for the flow problem is proposed based on the analysis on the pressure gradient and some supporting numerical simulations. The critical De number for various model fluids is given. It increases with the geometrical aspect ratio h/R (half channel width/cylinder radius) and depends on a viscosity ratio ${\beta}$(polymer viscosity/total viscosity) of the model. A shear thinning first normal stress coefficient will delay the instability. An excellent agreement between the predicted critical Deborah number and reported experiments is obtained.

• Nuclear Engineering and Technology
• /
• v.51 no.1
• /
• pp.310-316
• /
• 2019

Fluidelastic instability is a key issue in steam generator tube bundles subjected in cross-flow. With a low flow velocity, a large amplitude vibration of the tube observed by many researchers. However, the mechanism of this vibration is seldom analyzed. In this paper, the mechanism of cross-flow effects on fluidelastic instability of tube bundles was investigated. Analysis reveals that when the system reaches the critical state, there would be two forms, with two critical velocities, and thus two expressions for the critical velocities were obtained. Fluidelastic instability experiment and numerical analysis were conducted to obtain the critical velocity. And, if system damping is small, with increases of the flow velocity, the stability behavior of tube array changes. At a certain flow velocity, the stability of tube array reaches the first critical state, a dynamic bifurcation occurs. The tube array returns to a stable state with continues to increase the flow velocity. At another certain flow velocity, the stability of tube array reaches the second critical state, another dynamic bifurcation occurs. However, if system damping is big, there is only one critical state with increases the flow velocity. Compared the results of experiments to numerical analysis, it shows a good agreement.

• Nuclear Engineering and Technology
• /
• v.53 no.8
• /
• pp.2499-2508
• /
• 2021

Passive containment cooling systems (PCCSs) have been actively studied to improve the inherent safety of nuclear power plants. Hered, we present two concepts, open-loop PCCS (OL-PCCS) and closed-loop PCCS (CL-PCCS), applicable to the PWR with a concrete-type containment. We analyzed the heat-removal performance and flow instability of these PCCS concepts using the GOTHIC code. In both cases, PCCS performance improved when a passive containment cooling heat exchanger (PCCX) was installed in the lower part of the containment building. The OL-PCCS was found to be superior in terms of heat-removal performance. However, in terms of flow instability, the OL-PCCS was more vulnerable than the CL-PCCS. In particular, the possibility of flow instability was higher when the PCCX was installed in the upper part of the containment. Therefore, the installation location of the OL-PCCS should be restricted to minimize flow instability. Conversely, a CL-PCCS can be installed without any positional restriction by adjusting the initial system pressure within the loop, which eliminates flow instability. These results could be used as base data for the thermo-hydraulic evaluation of PCCS in PWR with a large dry concrete-type containment.

• Journal of Mechanical Science and Technology
• /
• v.19 no.9
• /
• pp.1811-1820
• /
• 2005

Thermo-acoustic instability of methane/ air flames in an industrial gas-turbine combustor is numerically investigated adopting CFD analysis. The combustor has 37 EV burners through which methane and air are mixed and then injected into the chamber. First, steady fuel! air mixing and flow characteristics established by the burner are investigated by numerical analysis with single burner. And then, based on information on the flow data, the burners are modeled numerically via equivalent swirlers, which facilitates the numerical analysis with the whole combustion system including the chamber and numerous burners. Finally, reactive flow fields within the chamber are investigated numerically by unsteady analysis and thereby, spontaneous instability is simulated. Based on the numerical results, scaling analysis is conducted to find out the instability mechanism in the combustor and the passive control method to suppress the instability is proposed and verified numerically.

• Nuclear Engineering and Technology
• /
• v.51 no.8
• /
• pp.2026-2033
• /
• 2019

Fluidelastic instability of a tube array is a key factor of the security of a nuclear power plant. An unsteady model of the fluidelastic instability of a tube array subjected to two-phase flow was developed to analyze the fluidelastic instability of tube bundles in two-phase flow. Based on this model, a computational program was written to calculate the eigenvalue and the critical velocity of the fluidelastic instability. The unsteady model and the program were verified by comparing with the experimental results reported previously. The influences of void fraction and the tube's material properties on the critical velocity were investigated. Numerical results showed that, with increasing the void fraction of the two-phase flow, the tube array becomes more stable. The results indicate that the critical velocities of the tube array made of stainless are much higher than those of the other two tube arrays within void fraction ranging from 20% to 80%.

• Journal of Mechanical Science and Technology
• /
• v.17 no.9
• /
• pp.1388-1396
• /
• 2003

The hydrodynamic instability of the three-dimensional boundary layer on a rotating disk introduces a periodic modulation of the mean flow in the form of stationary cross flow vortices. The instability labeled Type II by Faller occurs first at lower Reynolds number than that of well known Type I instability. Detailed numerical values of the amplification rates, neutral curves and other characteristics of the two instabilities have been calculated over a wide range of parameters. Presented are the neutral stability results concerning the two instability modes by solving the appropriate linear stability equations reformulated not only by considering whole convective terms but also by correcting some errors in the previous stability equations. The present stability results agree with the previously known ones within reasonable limit. Consequently, the flow is found to be always stable for a disturbance whose dimensionless wave number is greater than 0.75. Some spatial amplification contours have been computed for the stationary disturbance wave, whose azimuth angle $\varepsilon$= 11.29$^{\circ}$ to 15$^{\circ}$ and for the moving disturbance wave, whose azimuth angle $\varepsilon$ = 12.5$^{\circ}$ to 15$^{\circ}$. Also, some temporal amplification contours have been computed for the stationary disturbance wave, whose azimuth angle $\varepsilon$= 11.29$^{\circ}$ to 15$^{\circ}$ and for the moving disturbance wave, whose azimuth angle $\varepsilon$= 12$^{\circ}$ to 15$^{\circ}$. The flow instability was observed by using a white titanium tetrachloride gas over rotating disk system. When the numerical results are compared to the present experimental data, the numerical results agree quantitatively, indicating the existence of the selective frequency mechanism.

• Wind and Structures
• /
• v.26 no.6
• /
• pp.355-367
• /
• 2018

This paper studies the aerodynamic stability of a tensioned, geometrically nonlinear orthotropic membrane structure with hyperbolic paraboloid in sag direction. Considering flow separation, the wind field around membrane structure is simulated as the superposition of a uniform flow and a continuous vortex layer. By the potential flow theory in fluid mechanics and the thin airfoil theory in aerodynamics, aerodynamic pressure acting on membrane surface can be determined. And based on the large amplitude theory of membrane and D'Alembert's principle, interaction governing equations of wind-structure are established. Then, under the circumstance of single-mode response, the Bubnov-Galerkin approximate method is applied to transform the complicated interaction governing equations into a system of second-order nonlinear differential equation with constant coefficients. Through judging the frequency characteristic of the system characteristic equation, the critical velocity of divergence instability is determined. Different parameter analysis shows that the orthotropy, geometrical nonlinearity and scantling of structure is significant for preventing destructive aerodynamic instability in membrane structures. Compared to the model without considering flow separation, it's basically consistent about the divergence instability regularities in the flow separation model.