Numerical Study of Effects of Initial Flow Field on Taylor Vortex Flow

In concentrically rotating double cylinders consisting of a stationary outer cylinder and a rotating inner cylinder, Taylor vortex flow appears. Taylor vortex flow occurs in journal bearings, various fluid machineries, containers for chemical reaction, and other rotating components. Therefore, the analysis of the flow structure of Taylor vortex flow is highly effective for its control. The main parameters that determine the modes of Taylor vortex flow of a finite length are the aspect ratio Γ, Reynolds number Re. Γ is defined as the ratio of the cylinder length to the gap length between cylinders, and Re is determined on the basis of the angular speed of the inner cylinder. Γ was set to be 3.2, 4.8 and 6.8, and Re to be values in the range from 100 to 1000 at intervals of 100. Thus far, a large number of studies on Taylor vortex flow have been carried out; however, the effects of the differences in initial conditions have not yet been sufficiently clarified. In this study, we changed the initial flow field between the inner and outer cylinders in a numerical analysis, and examined the resulting changes in the mode formation and bifurcation processes. In this study, the initial speed distribution factor α was defined to be a function of the initial flow field and set to be 1.0, 0.999, 0.9 and 0.8 for the calculation. As a result, a difference was observed in the final mode depending on the difference in α for each Γ. From this finding, non-uniqueness, which is a major characteristic of Taylor vortex flow, was confirmed. However, no regularities regarding the difference in mode formation were found and the tendency of the mode formation process was not specified. Moreover, the processes of developing the vortex resulting in different final modes were monitored over time by visual observation. Similar flow behaviors were initially observed after the start of the calculation. Then, a bifurcation point, at which the flow changed to a mode depending on α, was observed, and finally the flow became steady. In addition, there was also a difference in the time taken for the flow to reach the steady state. These findings are based on only visual observation. Accordingly, a more detailed analysis at each lattice point and a comparison of physical quantities, such as kinetic energy and enstrophy, will be our future tasks.