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

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

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.

Responsiveness of human male volunteers to immunization with ovine follicle stimulating hormone vaccine: Results of a pilot study

A study of 140 days duration was performed to examine if human male volunteers (n = 5) respond to ovine follicle stimulating hormone (oFSH) immunization (administered adsorbed on Alugel on days 1, 20, 40 and 70) by producing antibodies capable of both binding and neutralizing bioactivity of human FSH. The kinetics of antibody production for both the immunogen (oFSH) and the cross-reactive antigen (hFSH) were essentially similar, The volunteers responded only to the first two immunizations, The boosters given on days 40 and 70 were ineffective, probably because of the presence of substantial amounts of circulating antibody to oFSH. Of the antibodies generated to oFSH, 25-45% bound hFSH with a mean binding affinity of 0.65 x 10(9) +/- 0.53 M(-1). The binding capacities at the time of high (30-80 days of immunization) and low (>110 days) titres were 346 +/- 185 and 10.5 +/- 5.8 ng hFSH/ml respectively, During the period of high titre, free serum FSH (value in normal males 1-5 ng/ml) was not monitorable, A 50 mu l aliquot of the antiserum obtained from different volunteers between days 30 and 80 and on day 140 blocked binding of I-125-labelled hFSH to its receptor by 82 +/- 9.7 and 53 +/- 12.2% respectively, The antibody produced was specific for FSH, and no significant change in the values of related glycoprotein hormones (luteinizing hormone/testosterone and thyroid stimulating hormone/thyroxine) were recorded, Seminal plasma transferrin, a marker of Sertoli cell as well as of seminiferous tubular function, showed marked reduction (30-90%) following immunization with oFSH. Considering that endogenous FSH remained neutralized for approximately one sperm cycle only (65 days), the reduction in sperm counts (30-74%) exhibited by some volunteers is encouraging, Immunization with oFSH did not result in any significant changes in haematology, serum biochemistry or hormonal profiles, There was no production of antibodies capable of interacting with non-specific tissues, It is concluded that it should be possible to obtain a sustained long-term blockade of endogenous FSH action in men by using oFSH as an immunogen, This is a prerequisite for obtaining significant reduction in the quality and quantity of spermatozoa produced, thus leading to infertility.

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.

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.