Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Nonsimilar mixed convection flow of a non-Newtonian fluid past a vertical wedge

The flow and heat transfer characteristics of a second-order fluid over a vertical wedge with buoyancy forces have been analysed. The coupled nonlinear partial differential equations governing the nonsimilar mixed convection flow have been solved numerically using Keller box method. The effects of the buoyancy parameter, viscoelastic parameter, mass transfer parameter, pressure gradient parameter, Prandtl number and viscous dissipation parameter on the skin friction and heat transfer have been examined in detail. Particular cases of the present results match exactly with those available in the literature.

Vibration and Stability of Fluid Conveying Pipes with Stochastic Parameters

Flexible cantilever pipes conveying fluids with high velocity are analysed for their dynamic response and stability behaviour. The Young's modulus and mass per unit length of the pipe material have a stochastic distribution. The stochastic fields, that model the fluctuations of Young's modulus and mass density are characterized through their respective means, variances and autocorrelation functions or their equivalent power spectral density functions. The stochastic non self-adjoint partial differential equation is solved for the moments of characteristic values, by treating the point fluctuations to be stochastic perturbations. The second-order statistics of vibration frequencies and mode shapes are obtained. The critical flow velocity is-first evaluated using the averaged eigenvalue equation. Through the eigenvalue equation, the statistics of vibration frequencies are transformed to yield critical flow velocity statistics. Expressions for the bounds of eigenvalues are obtained, which in turn yield the corresponding bounds for critical flow velocities.

A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

The two-step positive temperature coefficient in resistance characteristics of n-(Ba,Pb)TiO3 ceramics created through inclusion of grain boundary modifier additives

Donor-doped n-(Ba,Pb)TiO3 polycrystalline ceramics exhibit distinctly two-step positive temperature coefficient of resistance (PTCR) characteristics when formulated with suitable combinations of B2O3 and Al2O3 as grain boundary modifiers by heterogeneous addition. B2O3 or Al2O3 when added singularly resulted in either steep or broad PTCR jumps respectively across the phase transition. The two-step PTCR is attributed to the activation of the acceptor states, created through B2O3 and Al2O3, for various temperature regimes above the Curie point (T-c). The changing pattern of trap states is evident from the presence of Ti4+-O--Al3+ type hole centres in the grain boundary layer regions, identified in the electron paramagnetic resonance (EPR) spectra. That charge redistribution occurs among the inter-band gap defect states on crossing the Curie temperature is substantiated by the temperature coefficient in the EPR results. Capacitance-voltage results clearly show that there is an increase in the density of trap states with the addition of B2O3 and Al2O3. The spread in energy values of these trap states is evident from the large change in barrier height (phi similar or equal to 0.25-0.6 eV) between 500 and 650 K.

Unsteady three-dimensional stagnation point flow of a viscoelastic fluid

The unsteady three-dimensional stagnation point Bow of a viscoelastic fluid has been studied. Both nodal and saddle point regions of How have been considered. The unsteadiness in the Bow field is caused by the free stream velocity which varies arbitrarily with time. The governing boundary layer equations represented by a system of nonlinear partial differential equations have been solved numerically using a finite-difference scheme along with the quasilinearization technique in the nodal point region and a finite-difference scheme in combination with the parametric differentiation technique in the saddle point region. The skin friction coefficients for the viscoelastic fluid are found to be significantly less than those of the Newtonian fluid. The skin friction and heat transfer increase due to suction and reduce due to injection. The heat transfer at the wall increases with the Prandtl number. There is a flow reversal in the y-component of the velocity in the saddle point region. The absolute value of c (<<<0) for which reversal takes place is less than that of the Newtonian fluid. (C) 1997 Elsevier Science Ltd.

Spectral asymptotics of the Helmholtz model in fluid-solid structures

A model representing the vibrations of a coupled fluid-solid structure is considered. This structure consists of a tube bundle immersed in a slightly compressible fluid. Assuming periodic distribution of tubes, this article describes the asymptotic nature of the vibration frequencies when the number of tubes is large. Our investigation shows that classical homogenization of the problem is not sufficient for this purpose. Indeed, our end result proves that the limit spectrum consists of three parts: the macro-part which comes from homogenization, the micro-part and the boundary layer part. The last two components are new. We describe in detail both macro- and micro-parts using the so-called Bloch wave homogenization method. Copyright (C) 1999 John Wiley & Sons, Ltd.

The onset of oscillatory Marangoni convection in a two-layer system of conducting fluid in the presence of a uniform magnetic field

A combination of numerical and analytical techniques is used to analyse the effect of magnetic field and encapsulated layer on the onset of oscillatory Marangoni instability in a two layer system. Oscillatory Marangoni instability is possible for a deformed free surface only when the system is heated from above. It is observed that the existence of a second layer has a positive effect on Marangoni overstability with magnetic field whereas it has an opposite effect without magnetic field.

Unsteady axisymmetric stagnation-point flow of a viscous fluid on a cylinder

The unsteady viscous flow in the vicinity of an axisymmetric stagnation point of an infinite circular cylinder is investigated when both the free stream velocity and the velocity of the cylinder vary arbitrarily with time. The cylinder moves either in the same direction as that of the free stream or in the opposite direction. The flow is initially (t = 0) steady and then at t > 0 it becomes unsteady. The semi-similar solution of the unsteady Navier-Stokes equations has been obtained numerically using an implicit finite-difference scheme. Also the self-similar solution of the Navier-Stokes equations is obtained when the velocity of the cylinder and the free stream velocity vary inversely as a linear function of time. For small Reynolds number, a closed form solution is obtained. When the Reynolds number tends to infinity, the Navier-Stokes equations reduce to those of the two-dimensional stagnation-point flow. The shear stresses corresponding to stationary and the moving cylinder increase with the Reynolds number. The shear stresses increase with time for the accelerating flow but decrease with increasing time for the decelerating flow. For the decelerating case flow reversal occurs in the velocity profiles after a certain instant of time. (C) 1999 Elsevier Science Ltd. All rights reserved.

MHD flow over a moving plate in a rotating fluid with magnetic field, Hall currents and free stream velocity

The non-similar boundary layer flow of a viscous incompressible electrically conducting fluid over a moving surface in a rotating fluid, in the presence of a magnetic field, Hall currents and the free stream velocity has been studied. The parabolic partial differential equations governing the flow are solved numerically using an implicit finite-difference scheme. The Coriolis force induces overshoot in the velocity profile of the primary flow and the magnetic field reduces/removes the velocity overshoot. The local skin friction coefficient for the primary flow increases with the magnetic field, but the skin friction coefficient for the secondary flow reduces it. Also the local skin friction coefficients for the primary and secondary flows are reduced due to the Hall currents. The effects of the magnetic field, Hall currents and the wall velocity, on the skin friction coefficients for the primary and secondary flows increase with the Coriolis force. The wall velocity strongly affects the flow field. When the wall velocity is equal to the free stream velocity, the skin friction coefficients for the primary and secondary flows vanish, but this does not imply separation. (C) 2002 Published by Elsevier Science Ltd.

Stability of fluid flow past a membrane

The stability of fluid flow past a membrane of infinitesimal thickness is analysed in the limit of zero Reynolds number using linear and weakly nonlinear analyses. The system consists of two Newtonian fluids of thickness R* and H R*, separated by an infinitesimally thick membrane, which is flat in the unperturbed state. The dynamics of the membrane is described by its normal displacement from the flat state, as well as a surface displacement field which provides the displacement of material points from their steady-state positions due to the tangential stress exerted by the fluid flow. The surface stress in the membrane (force per unit length) contains an elastic component proportional to the strain along the surface of the membrane, and a viscous component proportional to the strain rate. The linear analysis reveals that the fluctuations become unstable in the long-wave (alpha --> 0) limit when the non-dimensional strain rate in the fluid exceeds a critical value Lambda(t), and this critical value increases proportional to alpha(2) in this limit. Here, alpha is the dimensionless wavenumber of the perturbations scaled by the inverse of the fluid thickness R*(-1), and the dimensionless strain rate is given by Lambda(t) = ((gamma) over dot* R*eta*/Gamma*), where eta* is the fluid viscosity, Gamma* is the tension of the membrane and (gamma) over dot* is the strain rate in the fluid. The weakly nonlinear stability analysis shows that perturbations are supercritically stable in the alpha --> 0 limit.

Unsteady mhd flow and heat transfer on a rotating disk in an ambient fluid

An unsteady flow and heat transfer of a viscous incompressible electrically conducting fluid over a rotating infinite disk in an otherwise ambient fluid are studied. The unsteadiness in the flow field is caused by the angular velocity of the disk which varies with time. The magnetic field is applied normal to the disk surface. The new self-similar solution of the Navier-Stokes and energy equations is obtained numerically. The solution obtained here is not only the solution of the Navier-Stokes equations, but also of the boundary layer equations. Also, for a simple scaling factor, it represents the solution of the flow and heat transfer in the forward stagnation-point region of a rotating sphere or over a rotating cone. The asymptotic behaviour of the solution for a large magnetic field or for a large independent variable is also examined. The surface shear stresses in the radial and tangential directions and the surface heat transfer increase as the acceleration parameter increases. Also the surface shear stress in the radial direction and the surface heat transfer decrease with increasing magnetic field, but the surface shear stress in the tangential direction increases. (C) 2002 Editions scientifiques et medicales Elsevier SAS. All rights reserved.