Hopf bifurcations to quasi-periodic solutions for the two-dimensional plane Poiseuille flow

This paper studies various Hopf bifurcations in the two-dimensional plane Poiseuille problem. For several values of the wavenumber α, we obtain the branch of periodic flows which are born at the Hopf bifurcation of the laminar flow. It is known that, taking α ≈ 1, the branch of periodic solutions has several Hopf bifurcations to quasi-periodic orbits. For the first bifurcation, calculations from other authors seem to indicate that the bifurcating quasi-periodic flows are stable and subcritical with respect to the Reynolds number, Re. By improving the precision of previous works we find that the bifurcating flows are unstable and supercritical with respect to Re. We have also analysed the second Hopf bifurcation of periodic orbits for several α, to find again quasi-periodic solutions with increasing Re. In this case the bifurcated solutions are stable to superharmonic disturbances for Re up to another new Hopf bifurcation to a family of stable 3-tori. The proposed numerical scheme is based on a full numerical integration of the Navier-Stokes equations, together with a division by 3 of their total dimension, and the use of a pseudo-Newton method on suitable Poincaré sections. The most intensive part of the computations has been performed in parallel. We believe that this methodology can also be applied to similar problems.     

Forced convective flow and heat transfer past an unconfined blunt headed cylinder at different angles of incidence

In the present paper, a numerical investigation has been carried out to study the forced convective flow and heat transfer characteristics past a blunt-headed cylinder in crossflow. Employing air as an operating fluid, calculations are carried out for a range of Reynolds number (Re) from 40 to 160. The angle of incidence is varied in the range of 0^∘ ≤ α ≤ 180 ^∘. The thermofluid features of flow and heat transport are analysed in detail for different angles of incidence. To analyse the aerodynamic characteristics, several parameters such as drag and lift coefficients, moment coefficient, Strouhal number, recirculation length, and local time-averaged vorticity flux have been calculated. Furthermore, a stability analysis has been undertaken by using the Stuart Landau equation to enumerate the critical Reynolds number at each angle of incidence. Heat transfer characteristics are studied by computing local and time-averaged values of Nusselt numbers. When compared to a rectangular cylinder, a blunt-headed cylinder exhibits an enhanced heat transfer rate. In the end, an entropy generation analysis has been carried out to study the effects of Re and angle of incidence on the efficiency of thermofluid transport characteristics.     

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode’“dual vortex pair mode” downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f _m in addition to the vortex shedding frequency f _s is identified. Frequency coupling between f _m and f _s occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.     

Numerical computation of end plate effect on Taylor vortices between rotating conical cylinders

End plate effect on Taylor vortices between rotating conical cylinders is studied by numerical method in this paper. We suppose that the inner cone rotates together with the end plate at the top and the outer one as well as the end plate at the bottom remains at rest. It is found that the instability sets in at a critical Reynolds number about Re = 80. Increase Re to about Re = 200 the first single clockwise vortex is formed near the top of the flow system. Further increase Re to about Re = 440 another clockwise vortex is formed under the first one. At about Re = 448 the third vortex is formed which rotates in counterclockwise direction between the first two vortices. With increasing of Re the process continues. Finally, a configuration is obtained with an odd number of vortices in the annulus at about Re = 700, which confirms the experimental observation. Moreover, the local extreme values of pressure and velocity are achieved at the adjacent lines between neighboring vortices or at the medium lines of vortices. The effect of gap size on vortices is also considered. It is shown that the number of vortices increases with decreasing of the gap size and the end plates play an important role in the parity of the number of the vortices.     

Numerical study on viscoelastic fluid flow past a rigid body

Viscoelastic non-Newtonian fluids can be achieved by adding a small amount of polymer additives to a Newtonian fluid. In this paper, numerical simulations are used to investigate the influence of such polymer additives on the behavior of flow past a circular cylinder. A numerical method is proposed that discretizes the non-linear viscoelastic system on a uniform Cartesian grid, with a penalization method to model the presence of the cylinder. The drag of the cylinder and the flow behavior under the effect of different Reynolds numbers ($\mathrm{Re}$), Weissenberg numbers (Wi) and polymer viscosity ratios (ε) are studied. Numerical results show that different flow characteristics are exhibited in different parameter zones. The polymer viscosity ratio plays an important role at low Weissenberg and Reynolds numbers, but as the Reynolds and Weissenberg numbers increase, the influence of ε weakens. The drag force of the cylinder is mostly affected by the Reynolds and Weissenberg numbers. At low Reynolds numbers, the drag of the cylinder and the flow fields are only affected by a large value of Wi when the elastic forces are strong. Non-trivial drag reduction occurs only when there is vortex shedding in the wake flow, whereas drag enhancement happens when the vortex shedding is inhibited.     

Buoyancy-assisted flow of yield stress fluids past a cylinder: Effect of shape and channel confinement

The aiding-buoyancy mixed convection heat transfer in Bingham plastic fluids from an isothermal cylinder of elliptical and circular shape in a vertical adiabatic channel is numerically investigated. For a fixed shape of the elliptical cylinder E = 2 (ratio of major to minor axes), the effect of confinement is studied for three values of blockage ratio, B, defined as the ratio of the channel width to the circumference of the cylinder/π, as 6.5, 2.17 and 1.3. In order to delineate the role of cross-section of the cylinder, results are also presented here for a circular cylinder of the same heat transfer area as the elliptical cylinder. The results presented herein span the range of conditions as: Bingham number, 0 ≤ Bn ≤ 100, Reynolds number, 1 ≤ Re ≤ 40, and Prandtl number, 1 ≤ Pr ≤ 100 over the range of Richardson number Ri = 0 (pure forced convection) to Ri = 10. Extensive results on drag coefficient, local and surface averaged values of the Nusselt number and yield surfaces are presented herein to elucidate the combined effects of buoyancy, blockage ratio and fluid yield stress. The morphology of the yield surfaces shows that the unyielded plug regions formed upstream and downstream of the cylinder grow faster at low Reynolds numbers with the increasing yield stress effects under the weak buoyancy forces, i.e., small values of Grashof or Richardson number. The heat transfer enhancement is observed with the increasing channel-confinement due to the sharpening of the temperature gradients near the surface of the cylinder. The average Nusselt number shows a positive dependence on the Reynolds number, Prandtl number and Richardson number irrespective of the shape of the cylinder or the type of fluid. By employing the modified definitions of the dimensionless parameters (based on the two choices of the overall effective fluid velocity), predictive correlations have been established for estimating the value of the average Nusselt number in a new application.     

Absolute values and real parts for functions in the smirnov class

Let N₊ denote the Smirnov class on the open unit disc D. It is easy to see that for any outer function g in N₊, there exists a function G in N₊ such that |g|; ≤ ReG on δ. We describe such a G. In general, G may not be outer. In this paper, a necessary and sufficient condition on g is given for the existence of an outer function G such that |;g|; < ReG. When g belongs to the Hardy space H¹, G is trivially given as the Herglotz integral of |;g|;.     

Thermal analysis of moving induction heating of a hollow cylinder with subsequent spray cooling: Effect of velocity,initial position of coil,and geometry

In this paper, temperature analysis of the complete process of moving induction heat treatment is performed using numerical methods. A non-linear and transient magneto-thermal coupled problem with a moving coil which is considered as moving heat source, is investigated by an efficient finite-element procedure. A vertical hollow circular cylinder is heated by the moving coil at a given velocity along it, and the heated parts then quenched by a moving water–air spray. The effects of natural convection with air on the both inner and outer surfaces of cylinder, and also radiation of outer surface of cylinder with ambient are taken into account. For quenching of work-piece, a specific kind of atomized spray cooling which utilizes a mixture of water and air with different mass fractions is used. This procedure includes moving boundary conditions, temperature-dependent properties, and change in magnetic permeability of specified alloy at the Curie temperature. Obtained numerical results have been verified by comparison with analytical solutions using Green’s function methods. Also, the effect of velocity, initial position of inductor and inner to outer radius ratio on temperature distribution are investigated.     

Predicting vortex-induced vibration from driven oscillation results

Two-dimensional simulations of flow past both an elastically-mounted cylinder and an externally-driven oscillating cylinder were performed at a Reynolds number of Re = 200. The results were compared to determine if the oscillations of the driven-oscillation model were consistent with the oscillations observed in the elastically-mounted system. It was found that while this is the case, there is considerable sensitivity to input forcing. This sensitivity could explain observed discrepancies between experimental results for the two systems.     

A numerical simulation of flow between two rotating coaxial frustum cones

The behavior of flow between two coaxial frustum cones, with the inner one rotating and the outer stationary, is studied in this paper. It is found that the fluid at the outlet does not flow out directly, but flows up till a certain height. This reflux generates a vortex area with a quite large velocity and pressure magnitude. This reflux area, between Z/H = 0.05 and 0.30, has the trend to move up with increasing Reynolds number Re. The velocity magnitude is linear in the radial direction if the Re is small. This linear relation converts to quasi-quadratic function as the Re increasing. If the frustum cone inclination is small, the flow will tend to be unstable with a quite large velocity and pressure magnitude. Finally, a comparison is made with Taylor–Couette flow.     

The effect of conjugate heat transfer on MHD mixed convection about a vertical slender hollow cylinder

The problem of steady laminar magnetohydrodynamic (MHD) mixed convection heat transfer about a vertical slender hollow cylinder is studied numerically, under the effect of wall conduction. A uniform magnetic field is applied perpendicular to the cylinder. The non-similar solutions using the Keller box method are obtained. The wall conduction parameter, the magnetic parameter and the Richardson number are the main parameters. For various values of these parameters the local skin friction and local heat transfer parameters are determined. The validity of the methodology is checked by comparing the results with those available in the open literature and a fairly good agreement is observed. Finally, it is determined that the local skin friction and the local heat transfer coefficients increase with an increase the magnetic parameter Mn and buoyancy parameter Ri and decrease with conjugate heat transfer parameter p.     

Numerical evaluation of the suppression effect of a free-to-rotate triangular fairing on the vortex-induced vibration of a circular cylinder

The suppression of vortex-induced vibration (VIV) of a circular cylinder with a free-to-rotate triangular fairing in the Reynolds number range of Re = 1100–6100 is numerically investigated using computational fluid dynamics. The unsteady Reynolds-averaged Navier–Stokes equations and the shear stress transport k–ω turbulence model coupled with an improved fourth-order Runge–Kutta method are used to solve the wake flow, the structure's vibration, and the fairing's rotation. The computational model is validated with the available experimental results for a cylinder with an attached short-tail fairing. The numerical results indicate that the triangular fairing has a positive role in suppressing vibration when it achieves a stable position deflected from the flow direction. The suppression effect is sensitive to the incoming flow velocity. The fairing shifts from a stable state to an unstable one when the flow velocity varies. Therefore, maintaining the hydrodynamic stability of the fairing is the key to achieving success in vibration suppression, and the stability is dependent on the characteristic length and the rotational friction. Although the strong flapping of the 70° triangular fairing excites a more vigorous vibration, it may be used as an amplifier of VIV for energy harvesting.     

A theoretical consideration of a free convective boundary layer on an isothermal horizontal conic

Approximate analytical solution of simplified Navier–Stokes and Fourier–Kirchhoff equations describing free convective heat transfer from isothermal surface of horizontal conic of the base angle α has been presented. The solution is based on the typical for natural convection assumption that the normal to the surface component of velocity is negligibly small in comparison with the tangential one. The results obtained for boundary cases of conic under considerations are in good agreement with known solutions for a horizontal cylinder α=π/2 and a vertical round plate α=0.     

Parametric study of liquid droplets impinging on porous surfaces

This work presents a numerical simulation of the fluid dynamics of a liquid droplet during impact/absorption onto a porous medium. The main focus of this paper is on a parametric study of the influence of the governing parameters upon the fluid flow characteristics. The problem is described in a non-dimensional form, and the influence of the main governing parameters is investigated, including their variation along the range of physical configurations of interest. This procedure revealed 7 main governing parameters: Reynolds number (Re), Darcy number (Da), porosity (ε), Froude number (Fr), Weber number (We), contact angle (θ) and the ratio between pore and particle diameter size in the porous substrate (α). The results indicate that the values of Da and Re are more related to the amount of momentum dissipation due to the drag of the solid matrix of the substrate, while the values of We, α and θ can be mainly related to capillary pressure.     

Numerical investigation of natural convection in a rectangular enclosure due to partial heating and cooling at vertical walls

A comprehensive numerical investigation on the natural convection in a rectangular enclosure is presented. The flow is induced due to the constant partial heating at lower half of the left vertical wall and partial cooling at upper half of the right vertical wall along with rest walls are adiabatic. In this investigation the Special attention is given to understand the effect of aspect ratio and heat source intensity i.e. Rayleigh number, Ra, on the fluid flow configuration as well as on the local and average heat transfer rates. The range of Rayleigh (Ra) and aspect ratio (A) is taken [10³, 10⁶] and [0.5, 4] respectively. The results are presented in terms of stream function (ψ), temperature (θ) and heat transfer rates (local Nusselt numbers Nu_L, and average Nusselt numbers Nu). The numerical experiments show that increasing of Ra implies the enhancement of thermal buoyancy force, which in turn increases the thermal convection in the cavity. As a result, the local as well as average heat transfer rate is expected to increase. The local transfer rate (Nu_L) is increases in the small region near the left vertical wall of the left wall of the cavity and after that it is decreases in the middle portion of heated region. And, it start to increase near to the middle point of left wall. It is also observed that the local heat transfer is increases as increases the aspect ratio. The average heat transfer rate (Nu) is increases as the aspect ratio A increases from 0.5 to 1 and beyond that it is decreases smoothly. It is also found that the heat transfer rate attains its maximum value at aspect ratio one.     

On the bifurcation analysis of a food web of four species

This paper is concerned with bifurcations of equilibria and the chaotic dynamics of a food web containing a bottom prey X, two competing predators Y and Z on X, and a super-predator W only on Y. Conditions for the existence of all equilibria and the stability properties of most equilibria are derived. A two-dimensional bifurcation diagram with the aid of a numerical method for identifying bifurcation curves is constructed to show the bifurcations of equilibria. We prove that the dynamical system possesses a line segment of degenerate steady states for the parameter values on a bifurcation line in the bifurcation diagram. Numerical simulations show that these degenerate steady states can help to switch the stabilities between two far away equilibria when the system crosses this bifurcation line. Some observations concerned with chaotic dynamics are also made via numerical simulations. Different routes to chaos are found in the system. Relevant calculations of Lyapunov exponents and power spectra are included to support the chaotic properties.     

Symmetry and bifurcation near families of solutions

Let X, Z and Λ be Banach spaces, M: X × Λ → Z a C¹-function, and assume that the equation M(x, λ) = 0 has a family of solutions for λ = 0. In this paper we consider the bifurcation of solutions from this family, for ¦λ¦ small, under the condition that both the unperturbed (λ = 0) and the perturbed (λ ≠ 0) equations have certain symmetry properties. The problem is reduced by the Liapunov-Schmidt method, and the bifurcation equations are solved by a straightforward use of the symmetry. As an application we obtain existence of certain periodic solutions for the undamped Duffing equation, a result recently obtained by Schmitt and Mazzanti using different methods.     

Pressure and flow in two dimensional numerical bifurcation models

The finite element method has been used to solve the Navier-Strokes equations for steady flow conditions in bifurcations. The results are presented as pressure, velocity and streamline plots at different Reynolds number. The three bifurcations considered have rigid walls and bifurcation angles of 0°, 20° and 180°. For the bifurcation with branch angles 0° and 20° there is flow separation along the inner wall of the outlet branches and large spatial pressure variations, these phenomena being more pronounced at the higher Reynolds numbers. For the bifurcation with a branch angle of 180° the high pressure gradients occured at the outer corner and for the high Reynolds number a vortex formation developed downstream of this corner.     

Streamwise oscillations of a cylinder beneath a free surface: Free surface effects on vortex formation modes

A computational study of a viscous incompressible two-fluid model with an oscillating cylinder is investigated at a Reynolds number of 200 and at a dimensionless displacement amplitude of A=0.13 and for the dimensionless forcing cylinder oscillation frequency-to-natural vortex shedding frequency ratios, f/f₀=1.5,2.5,3.5. Specifically, two-dimensional flow past a circular cylinder subject to forced in-line oscillations beneath a free surface is considered. The method is based on a finite volume discretization of the two-dimensional continuity and unsteady Navier-Stokes equations (when a solid body is present) on a fixed Cartesian grid. Two-fluid model based on improved volume-of-fluid method is used to discretize the free surface interface. The study focuses on the laminar asymmetric flow structure in the near wake region and lock-on phenomena at a Froude number of 0.2 and for the dimensionless cylinder submergence depths, h=0.25, 0.5 and 0.75. The equivorticity patterns and pressure distribution contours are used for the numerical flow visualization. The code validations in special cases show good comparisons with previous numerical results.