• Korea-Australia Rheology Journal
• /
• 제16권3호
• /
• pp.117-128
• /
• 2004

The effect of molecular parameters on the steady vortex behaviors in the inertial viscoelastic flow past a cylinder has been investigated. FENE-CR model was considered as a constitutive equation. A recently developed iterative solution method (Kim et al., (in press)) was found to be successfully applicable to the computation of inertial viscoelastic flows. The high-resolution computations were carried out to understand the detailed flow behaviors based on the efficient iterative solution method armed with ILU(0) type pre-conditioner and BiCGSTAB method. The discrete elastic viscous split stress-G/streamline upwind Petrov Galerkin (DEVSS-G/SUPG) formulation was adopted as a stabilization method. The vortex size decreased as elasticity increases. However, the vortex enhancement was also observed in the case of large extensibility, which means that the vortex behavior is strongly dependent upon the material parameters. The longitudinal gradient of normal stress was found to retard the formation of vortex, whereas the extensional viscosity played a role in the vortex enhancement. The present results are expected to be helpful for understanding the inertial vortex dynamics of viscoelastic fluids in the flow past a confined cylinder.

• 대한기계학회논문집B
• /
• 제26권2호
• /
• pp.261-268
• /
• 2002

Turbulent flow past a square cylinder confined in a channel is numerically investigated by Large Eddy Simulation(LES). The main objectives of this study are to verify the experimental results of Nakagawa et al.[Exp. in Fluids, Vol. 27, 3, pp. 284∼294, 1999] by LES and to obtain related flow information in detail. The LES results obtained are in excellent agreement with the experiment both qualitatively and quantitatively. The passive paticles numerically released into the flow field clearly show the barman vortex street. However, the vortices shed from the cylinder are significantly affected by the presence of the plates. Futhermore, periodic and alternating vortex-rollups are observed in the vicinity of the plates. The rolled-up vortex is convected downstream together with the corresponding Karman vortex forming a counter-rotating vortex pair. It is also revealed that the cylinder greatly enhances mixing process of the flow.

• Korea-Australia Rheology Journal
• /
• 제20권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.