On the unsteady rotational flow of a generalized Maxwell fluid through a circular cylinder

In this paper, the velocity field and the associated tangential stress corresponding to the rotational flow of a generalized Maxwell fluid within an infinite circular cylinder are determined by means of the Laplace and finite Hankel transforms. Initially, the fluid is at rest, and the motion is produced by the rotation of the cylinder about its axis with a unsteady angular velocity. The solutions that have been obtained are presented under series form in terms of the generalized G _a,b,c (, t)-functions. The similar solutions for the ordinary Maxwell and Newtonian fluids, performing the same motion, are obtained as special cases, when β → 1, respectively β → 1 and λ → 0, from general solutions. Finally, the solutions that have been obtained are compared by graphical illustrations, and the influence of the pertinent parameters on the fluid motion is also underlined by graphical illustrations.     

Exact solutions for the flow of second grade fluid in annulus between torsionally oscillating cylinders

The velocity field and the associated shear stress corresponding to the torsional oscillatory flow of a second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. At time t = 0, the fluid and both the cylinders are at rest and at t = 0 + , cylinders suddenly begin to oscillate around their common axis in a simple harmonic way having angular frequencies ω 1 and ω 2 . The obtained solutions satisfy the governing differential equation and all imposed initial and boundary conditions. The solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for Newtonian fluid are also obtained as limiting cases of our general solutions.     

Exact analytical solutions for axial flow of a fractional second grade fluid between two coaxial cylinders

The velocity field and the adequate shear stress corresponding to the longitudinal flow of a fractional second grade fluid, between two infinite coaxial circular cylinders, are determined by applying the Laplace and finite Hankel transforms. Initially the fluid is at rest, and at time t = 0⁺, the inner cylinder suddenly begins to translate along the common axis with constant acceleration. The solutions that have been obtained are presented in terms of generalized G functions. Moreover, these solutions satisfy both the governing differential equations and all imposed initial and boundary conditions. The corresponding solutions for ordinary second grade and Newtonian fluids are obtained as limiting cases of the general solutions. Finally, some characteristics of the motion, as well as the influences of the material and fractional parameters on the fluid motion and a comparison between models, are underlined by graphical illustrations.     

On the unsteady rotational flow of fractional Oldroyd-B fluid in cylindrical domains

This paper concerned with the unsteady rotational flow of fractional Oldroyd-B fluid, between two infinite coaxial circular cylinders. To solve the problem we used the finite Hankel and Laplace transforms. The motion is produced by the inner cylinder that, at time t=0⁺, is subject to a time-dependent rotational shear. The solutions that have been obtained, presented under series form in terms of the generalized G functions, satisfy all imposed initial and boundary conditions. The corresponding solutions for ordinary Oldroyd-B, fractional and ordinary Maxwell, fractional and ordinary second grade, and Newtonian fluids, performing the same motion, are obtained as limiting cases of general solutions.     

Some exact solutions for rotating flows of a generalized Burgers’ fluid in cylindrical domains

The velocity field and the adequate shear stress corresponding to the flow of a generalized Burgers’ fluid model, between two infinite co-axial cylinders, are determined by means of Laplace and finite Hankel transforms. The motion is due to the inner cylinder that applies a time dependent torsional shear to the fluid. The solutions that have been obtained, presented in series form in terms of usual Bessel functions J₁( ? ), J₂( ? ), Y₁( ? ) and Y₂( ? ), satisfy all imposed initial and boundary conditions. Moreover, the corresponding solutions for Burgers’, Oldroyd-B, Maxwell, second grade, Newtonian fluids and large-time transient solutions for generalized Burgers’ fluid are also obtained as special cases of the present general solutions. The effect of various parameters on large-time and transient solutions of generalized Burgers’ fluid is also discussed. Furthermore, for small values of the material parameters, λ₂ and λ₄ or λ₁, λ₂, λ₃ and λ₄, the general solutions corresponding to generalized Burgers’ fluids are going to those for Oldroyd-B and Newtonian fluids, respectively. Finally, the influence of the pertinent parameters on the fluid motion, as well as a comparison between models, is shown by graphical illustrations.     

Exact solutions for the unsteady rotational flow of an Oldroyd-B fluid with fractional derivatives induced by a circular cylinder

In this research article, the unsteady rotational flow of an Oldroyd-B fluid with fractional derivative model through an infinite circular cylinder is studied by means of the finite Hankel and Laplace transforms. The motion is produced by the cylinder, that after time t=0⁺, begins to rotate about its axis with an angular velocity Ωt ^p . The solutions that have been obtained, presented under series form in terms of the generalized G-functions, satisfy all imposed initial and boundary conditions. The corresponding solutions that have been obtained can be easily particularized to give the similar solutions for Maxwell and Second grade fluids with fractional derivatives and for ordinary fluids (Oldroyd-B, Maxwell, Second grade and Newtonian fluids) performing the same motion, are obtained as limiting cases of general solutions. The most important things regarding this paper to mention are that (1) we extracted the expressions for the velocity field and the shear stress corresponding to the motion of Second grade fluid with fractional derivatives as a limiting case of our general solutions corresponding to the Oldroyd-B fluid with fractional derivatives, this is not previously done in the literature to the best of our knowledge, and (2) the expressions for the velocity field and the shear stress are in the most simplified form, and the point worth mentioning is that these expressions are free from convolution product and the integral of the product of the generalized G-functions. Finally, the influence of the pertinent parameters on the fluid motion, as well as a comparison between models, is shown by graphical illustrations.     

RETRACTED ARTICLE: Flow of fractional Maxwell fluid between coaxial cylinders

This paper deals with the study of unsteady flow of a Maxwell fluid with fractional derivative model, between two infinite coaxial circular cylinders, using Laplace and finite Hankel transforms. The motion of the fluid is produced by the inner cylinder that, at time t = 0⁺, is subject to a time-dependent longitudinal shear stress. Velocity field and the adequate shear stress are presented under series form in terms of the generalized G and R functions. The solutions that have been obtained satisfy all imposed initial and boundary conditions. The corresponding solutions for ordinary Maxwell and Newtonian fluids are obtained as limiting cases of general solutions. Finally, the influence of the pertinent parameters on the fluid motion as well as a comparison between the three models is underlined by graphical illustrations.     

Steady mixed convection boundary-layer flow over a vertical flat surface in a porous medium filled with water at 4°C: variable surface heat flux

The steady mixed convection boundary-layer flow over a vertical impermeable surface in a porous medium saturated with water at 4°C (maximum density) when the surface heat flux varies as x ^m and the velocity outside the boundary layer varies as x ^(1+2m)/2, where x measures the distance from the leading edge, is discussed. Assisting and opposing flows are considered with numerical solutions of the governing equations being obtained for general values of the flow parameters. For opposing flows, there are dual solutions when the mixed convection parameter λ is greater than some critical value λ _c (dependent on the power-law index m). For assisting flows, solutions are possible for all values of λ. A lower bound on m is found, m > −1 being required for solutions. The nature of the critical point λ _c is considered as well as various limiting forms; the forced convection limit (λ = 0), the free convection limit (λ → ∞) and the limits as m → ∞ and as m → −1.     

Asymptotic approximation for the solution of a boundary-value problem with varying type of boundary conditions in a thick two-level junction

We consider a mixed boundary-value problem for a Poisson equation in a plane two-level junction Ω_ε that is the union of a domain Ω₀ and a large number 3N of thin rods with thickness of order . The thin rods are divided into two levels depending on their length. In addition, the thin rods from each level are ε-periodically alternated. The homogeneous Dirichlet conditions and inhomogeneous Neumann conditions are given on the sides of the thin rods from the first level and the second level, respectively. Using the method of matched asymptotic expansions and special junction-layer solutions, we construct an asymptotic approximation for the solution and prove the corresponding estimates in the Sobolev space H ¹(Ω_ε) as ε → 0 (N → +∞). Published in Neliniini Kolyvannya, Vol. 9, No. 3, pp. 336–355, July–September, 2006.     

Exact solutions for the flow of a generalized Oldroyd-B fluid induced by a constantly accelerating plate between two side walls perpendicular to the plate

The unsteady flow of an incompressible generalized Oldroyd-B fluid induced by a constantly accelerating plate between two side walls perpendicular to the plate has been studied using Fourier sine and Laplace transforms. The obtained solutions for the velocity field and shear stresses, written in terms of the generalized G and R functions, are presented as sum of the similar Newtonian solutions and the corresponding non-Newtonian contributions. For α = β = 1 and λ_r → λ these solutions are going to the corresponding Newtonian solutions. Furthermore, the solutions for generalized Maxwell fluids as well as those for ordinary Oldroyd-B and Maxwell fluids, performing the same motion, are also obtained as limiting cases of our general solutions. In the absence of the side walls, namely when the distance between the two walls tends to infinity, the solutions corresponding to the motion over an infinite constantly accelerating plate are recovered. For λ_r → 0 and β → 1, these solutions reduce to the known solutions from the literature. Finally, the effect of the material parameters on the velocity profile is spotlighted by means of the graphical illustrations.     

Motion of a sphere in an electrically conducting rotating fluid

The combined effect of rotation and magnetic field is investigated for the axisymmetric flow due to the motion of a sphere in an inviscid, incompressible electrically conducting fluid having uniform rotation far upstream. The steady-state linearized equations contain a single parameter α=1/2βR _m, β being the magnetic pressure number and R _m the magnetic Reynolds number. The complete solution for the flow field and magnetic field is obtained and the distribution of vorticity and current density is found. The induced vorticity is O(α⁴) and the current density is O(R _m) on the sphere.     

Combined heat and mass transfer along a vertical moving cylinder with a free stream

The mixed convection flow over a continuous moving vertical slender cylinder under the combined buoyancy effect of thermal and mass diffusion has been studied. Both uniform wall temperature (concentration) and uniform heat (mass) flux cases are included in the analysis. The problem is formulated in such a manner that when the ratio λ(= u _w/(u _w + u _∞), where u _w and u _∞ are the wall and free stream velocities, is zero, the problem reduces to the flow over a stationary cylinder, and when λ = 1 it reduces to the flow over a moving cylinder in an ambient fluid. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. We have also obtained the solution using a perturbation technique with Shanks transformation. This transformation has been used to increase the range of the validity of the solution. For some particular cases closed form solutions are obtained. The surface skin friction, heat transfer and mass transfer increase with the buoyancy forces. The buoyancy forces cause considerable overshoot in the velocity profiles. The Prandtl number and the Schmidt number strongly affect the surface heat transfer and the mass transfer, respectively. The surface skin friction decreases as the relative velocity between the surface and free stream decreases. Received on 17 May 1999     

Stokes＇ first problem for a viscoelastic fluid with the generalized Oldroyd-B model

The flow near a wall suddenly set in motion for a viscoelastic fluid with the generalized Oldroyd-B model is studied. The fractional calculus approach is used in the constitutive relationship of fluid model. Exact analytical solutions of velocity and stress are obtained by using the discrete Laplace transform of the sequential fractional derivative and the Fox H-function. The obtained results indicate that some well known solutions for the Newtonian fluid, the generalized second grade fluid as well as the ordinary Oldroyd-B fluid, as limiting cases, are included in our solutions. The project supported by the National Natural Science Foundation of China (10272067), the Doctoral Program Foundation of the Education Ministry of China (20030422046), the Natural Science Foundation of Shandong Province, China (Y2006A14) and the Research Foundation of Shandong University at Weihai. The English text was polished by Keren Wang.     

Stability of a fluid-saturated porous medium contained in a vertical cylinder heated from below by forced convection

A linear stability analysis determining the critical Rayleigh number R _c for onset of convection in a bounded vertical cylinder containing a fluid-saturated porous medium is performed for insulated sidewalls, isothermal top surface, and bottom surface heated by forced convection. This Newtonian heating of the bottom surface involves a Biot number Bi that allows consideration of the continuum of boundary conditions ranging from constant heat flux, with global minimum R _min=27.096 found as Bi→0, to isothermal, with global minimum R _min=4π² found as Bi→ ∞. In both cases and for most cylinder aspect ratios, incipient convection sets in as an asymmetric mode, though islands of aspect ratio exist where the onset mode is symmetric. Sample three-dimensional renderings of disturbance temperature distributions showing preferred modes at onset of convection for fixed Bi are provided and an analytical fit to R _min as a function of Bi is given.     

Mixed convection flow of second grade fluid along a vertical stretching flat surface with variable surface temperature

In the present study we have explored the effects of thermal buoyancy on flow of a viscoelastic second grade fluid past a vertical, continuous stretching sheet of which the velocity and temperature distributions are assumed to vary according to a power-law form. The governing differential equations are transformed into dimensionless form using appropriate transformations and then solved numerically. The methods here employed are (1) the perturbation method together with the Shanks transformation, (2) the local non-similarity method with second level of truncation and (3) the implicit finite difference method for values of ξ ( = Gr _x /Re _x ², defined as local mixed convection parameter) ranging in [0, 10]. The comparison between the solutions obtained by the aforementioned methods found in excellent agreement. Effects of the elasticity parameter λ on the skin-friction and heat transfer coefficients have been shown graphically for the fluids having the values of the Prandtl number equal to 0.72, 7.03 and 15.0. Effects of the viscoelastic parameter and the mixed convection parameter, ξ, on the temperature and velocity fields have also been studied. We notice that with the increase in visco-elastic parameter λ, velocity decreases whereas temperature increases and that velocity gradient is higher than that of temperature. On leave of absence from the Department of Mathematics, University of Dhaka, Bangladesh.     

Thermal convection in a porous layer: Effects of anisotropy and surface boundary conditions

The theory describing the onset of convection in a homogeneous porous layer bounded above and below by isothermal surfaces is extended to consider an upper boundary which is partly permeable. The general boundary condition p + λ ∂p/∂n = constant is applied at the top surface and the flow is investigated for various λ in the range 0 ⩽ λ < ∞. Estimates of the magnitude and horizontal distribution of the vertical mass and heat fluxes at the surface, the horizontally-averaged heat flux (Nusselt number) and the fraction of the fluid which recirculates within the layer are found for slightly supercritical conditions. Comparisons are made with the two limiting cases λ → ∞, where the surface is completely impermeable, and λ = 0, where the surface is at constant pressure. Also studied are the effects of anisotropy in permeability, ξ = K _H /K _V , and anisotropy is thermal conductivity, η = k _H /k _V , both parameters being ratios of horizontal to vertical quantities. Quantitative results are given for a wide variety of the parameters λ, ξ and η. In the limit ξ/η → 0 there is no recirculation, all fluid being converted out of the top surface, while in the limit ξ/η → ∞ there is full recirculation.     

Mixed Convection Boundary Layer Flow from a Horizontal Circular Cylinder Embedded in a Porous Medium Filled with a Nanofluid

Steady mixed convection boundary layer flow from an isothermal horizontal circular cylinder embedded in a porous medium filled with a nanofluid has been studied for both cases of a heated and cooled cylinder. The resulting system of nonlinear partial differential equations is solved numerically using an implicit finite-difference scheme. The solutions for the flow and heat transfer characteristics are evaluated numerically for various values of the governing parameters, namely the nanoparticle volume fraction φ and the mixed convection parameter λ. Three different types of nanoparticles are considered, namely Cu, Al₂O₃ and TiO₂. It is found that for each particular nanoparticle, as the nanoparticle volume fraction φ increases, the magnitude of the skin friction coefficient decreases, and this leads to an increase in the value of the mixed convection parameter λ which first produces no separation. On the other hand, it is also found that of all the three types of nanoparticles considered, for any fixed values of φ and λ, the nanoparticle Cu gives the largest values of the skin friction coefficient followed by TiO₂ and Al₂O₃. Finally, it is worth mentioning that heating the cylinder (λ > 0) delays separation of the boundary layer and if the cylinder is hot enough (large values of λ > 0), then it is suppressed completely. On the other hand, cooling the cylinder (λ < 0) brings the boundary layer separation point nearer to the lower stagnation point and for a sufficiently cold cylinder (large values of λ < 0) there will not be a boundary layer on the cylinder.     

Mixed convection from a cylinder oscillating vertically in a quiescent fluid

 The problem of heat convection from a vertically oscillating cylinder in a quiescent fluid is investigated. The governing equations of motion and energy are solved numerically in a non-inertial frame of references to determine the flow field and heat transfer characteristics under different conditions. The main dominating parameters are Keulegan–Carpenter number, KC, frequency parameter, β, Grashof number, Gr and Prandtl number, Pr. The ranges considered for these parameters are KC ≤ 10, β≤40 and Gr ≤ 10⁵ while Prandtl number is kept constant. The study revealed that the effect of amplitude and frequency of oscillation on heat transfer is strongly influenced by the Grashof number range. In the forced convection regime (Gr = 0), the increase of KC creates extensive vortex motion at all cylinder positions that leads to a significant increase in heat transfer. A similar trend, but with a lesser extent, is also observed for the increase of β. However, at high Grashof numbers, the effect of oscillation on heat convection is only significant at large values of KC. Received on 5 June 2000 / Published online: 29 November 2001     

Bifurcations of Heteroclinic Orbits

The search for traveling wave solutions of a semilinear diffusion partial differential equation can be reduced to the search for heteroclinic solutions of the ordinary differential equation ü − cu̇ + f(u) = 0, where c is a positive constant and f is a nonlinear function. A heteroclinic orbit is a solution u(t) such that u(t) → γ ₁ as t → −∞ and u(t) → γ ₂ as t → ∞ where γ ₁, γ ₂ are zeros of f. We study the existence of heteroclinic orbits under various assumptions on the nonlinear function f and their bifurcations as c is varied. Our arguments are geometric in nature and so we make only minimal smoothness assumptions. We only assume that f is continuous and that the equation has a unique solution to the initial value problem. Under these weaker smoothness conditions we reprove the classical result that for large c there is a unique positive heteroclinic orbit from 0 to 1 when f(0) = f(1) = 0 and f(u) > 0 for 0 < u < 1. When there are more zeros of f, there is the possibility of bifurcations of the heteroclinic orbit as c varies. We give a detailed analysis of the bifurcation of the heteroclinic orbits when f is zero at the five points −1 < −θ < 0 < θ < 1 and f is odd. The heteroclinic orbit that tends to 1 as t → ∞ starts at one of the three zeros, −θ, 0, θ as t → −∞. It hops back and forth among these three zeros an infinite number of times in a predictable sequence as c is varied.     

Fully developed flow of a viscoelastic film down a vertical cylindrical or planar wall

The one-dimensional, gravity-driven film flow of a linear (l) or exponential (e) Phan-Thien and Tanner (PTT) liquid, flowing either on the outer or on the inner surface of a vertical cylinder or over a planar wall, is analyzed. Numerical solution of the governing equations is generally possible. Analytical solutions are derived only for: (1) l-PTT model in cylindrical and planar geometries in the absence of solvent, b o [(h)\tilde]_s/([(h)\tilde]_s +[(h)\tilde]_p)=0\beta\equiv {\tilde{\eta}_s}/\left({\tilde{\eta}_s +\tilde{\eta}_p}\right)=0, where [(h)\tilde]_p\widetilde{\eta}_p and [(h)\tilde]_s\widetilde{\eta}_s are the zero-shear polymer and solvent viscosities, respectively, and the affinity parameter set at ξ = 0; (2) l-PTT or e-PTT model in a planar geometry when β = 0 and x 1 0\xi \ne 0; (3) e-PTT model in planar geometry when β = 0 and ξ = 0. The effect of fluid properties, cylinder radius, [(R)\tilde]\tilde{R}, and flow rate on the velocity profile, the stress components, and the film thickness, [(H)\tilde]\tilde{H}, is determined. On the other hand, the relevant dimensionless numbers, which are the Deborah, De=[(l)\tilde][(U)\tilde]/[(H)\tilde]De={\tilde{\lambda}\tilde{U}}/{\tilde{H}}, and Stokes, St=[(r)\tilde][(g)\tilde][(H)\tilde]²/([(h)\tilde]_p +[(h)\tilde]_s )[(U)\tilde]St=\tilde{\rho}\tilde{g}\tilde{\rm H}^{2}/\left({\tilde{\eta}_p +\tilde{\eta}_s} \right)\tilde{U}, numbers, depend on [(H)\tilde]\tilde{H} and the average film velocity, [(U)\tilde]\widetilde{U}. This makes necessary a trial and error procedure to obtain [(H)\tilde]\tilde{H} a posteriori. We find that increasing De, ξ, or the extensibility parameter ε increases shear thinning resulting in a smaller St. The Stokes number decreases as [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to zero for a film on the outer cylindrical surface, while it asymptotes to very large values when [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to unity for a film on the inner surface. When x 1 0\xi \ne 0, an upper limit in De exists above which a solution cannot be computed. This critical value increases with ε and decreases with ξ.