Stability of the flow of a fluid through a flexible tube at high Reynolds number

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 = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (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 do), 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 fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε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 Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) 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 [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 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.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Study of Ion Transport in Lithium Perchlorate-Succinonitrile Plastic Crystalline Electrolyte via Ionic Conductivity and in Situ Cryo-Crystallography

Ion transport mechanism in lithium perchlorate (LiClO4)-succinonitrile (SN), a prototype of plastic crystalline soft matter electrolyte is discussed in the context of solvent configurational isomerism and ion solvation. Contributions of both solvent configurational isomerism and ion solvation are reflected in the activation energy for ion conduction in 0-1 M LiClO4-SN samples. Activation energy due to solvent configurational changes, that is, trans-gauche isomerism is observed to be a function of salt content and decreases in presence of salt (except at high salt concentrations, e.g. 1 M LiClO4-SN). The remnant contribution to activation energy is attributed to ion-association. The X-ray diffraction of single crystals obtained using in situ cryo-crystallography confirms directly the observations of the ionic conductivity measurements. Fourier transform infrared spectroscopy and NMR line width measurements provide additional support to our proposition of ion transport in the prototype plastic crystalline electrolyte.

Nonlinear transport in semiconducting polymers at high carrier densities

Conducting and semiconducting polymers are important materials in the development of printed, flexible, large-area electronics such as flat-panel displays and photovoltaic cells. There has been rapid progress in developing conjugated polymers with high transport mobility required for high-performance field-effect transistors (FETs), beginning(1) with mobilities around 10(-4) cm(2) V-1 s(-1) to a recent report(2) of 1 cm(2) V-1 s(-1) for poly(2,5-bis(3-tetradecylthiophen-2-yl) thieno[3,2-b] thiophene) (PBTTT). Here, the electrical properties of PBTTT are studied at high charge densities both as the semiconductor layer in FETs and in electrochemically doped films to determine the transport mechanism. We show that data obtained using a wide range of parameters (temperature, gate-induced carrier density, source-drain voltage and doping level) scale onto the universal curve predicted for transport in the Luttinger liquid description of the one-dimensional `metal'.

A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Functional programming systems revisited

Functional Programming (FP) systems are modified and extended to form Nondeterministic Functional Programming (NFP) systems in which nondeterministic programs can be specified and both deterministic and nondeterministic programs can be verified essentially within the system. It is shown that the algebra of NFP programs has simpler laws in comparison with the algebra of FP programs. "Regular" forms are introduced to put forward a disciplined way of reasoning about programs. Finally, an alternative definition of "linear" forms is proposed for reasoning about recursively defined programs. This definition, when used to test the linearity of forms, results in simpler verification conditions than those generated by the original definition of linear forms.

Mechanism of glycine uptake by the silk glands of the silkworm Bombyx mori L

The transport of glycine in vitro into the silk glands of the silkworm has been studied. Glycine accumulates inside the tissue to a concentration higher than that present outside, indicating an active transport mechanism. The kinetics of uptake show a biphasic curve and two apparent Km values for accumulation, 0.33 mM and 5.00 mM. The effect of inhibitors on the energy metabolism of glycine transport is inconclusive. Exchange studies indicate the existence of two pools inside the gland, one that is easily removed by exchange and osmotic shock, and the other which is not. The results obtained conform with the carrier model of Britten and McClure concerning the amino-acid pool in E. coli.

Simulation of laminar flow in a three-dimensional lid-driven cavity by lattice Boltzmann method

Purpose - The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid-driven flow in a three-dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two-dimensional (2D) cavity. Design/methodology/approach - The second-order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25-4) and Reynolds numbers (0.01-1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier-Stokes (NS) equations and with available experimental data. Findings - The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side-walls in transverse direction (in 3D), due to which the vorticity in the core-region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry-plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid-motion are observed; these are weak at low. Reynolds numbers, but become quite strong at high Reynolds numbers. Originality/value - The findings will be useful in the study of variety of enclosed fluid flows.

Smoldering combustion of "incense" sticks - Experiments and modeling

This paper is concerned with the experimental and modeling studies on the smoldering rates of incense sticks as a function of ambient oxygen fraction in air, the flow velocity and size. The experimental results are obtained both for forward and reverse smolder conditions. The results are explained on the basis of surface combustion due to diffusion of oxygen to the surface by both free and forced convection supporting the heat transfer into the solid by conduction, into the stream by convection and the radiant heat transfer from the surface. The heat release at the surface is controlled by the convective transport of the oxidizer to the surface. To obtain the diffusion rates particularly for the reverse smolder, CFD calculations of fluid flow with along with a passive scalar are needed; these calculations have been made both for forward and reverse smolder. The interesting aspect of the CFD calculations is that while the Nusselt umber for forward smolder shows a clear root( Re-u) dependence ( Re-u = Flow Reynolds Number), the result for reverse smolder shows a peak in the variation with Reynolds number with the values lower than for forward smolder and unsteadiness in the flow beyond a certain flow rate. The results of flow behavior and Nusselt number are used in a simple model for the heat transfer at the smoldering surface to obtain the dependence of the smoldering rate on the diameter of the incense stick, the flow rate of air and the oxygen fraction. The results are presented in terms of a correlation for the non-dimensional smoldering rate with radiant flux from the surface and heat generation rate at the surface. The correlations appear reasonable for both forward and reverse smolder cases.

The oscillation of a spheroid along its axis in a non-Newtonian fluid

The present paper investigates the nature of the fluid flow when a spheroid is suspended in an infinitely extending elastico-viscous fluid defined by the constitutive equations given by Oldroyd or Rivlin and Ericksen, and is made to perform small amplitude oscillations along its axis. The solution of the vector wave equation is expressed in terms of the solution of the corresponding scalar wave equation, without the use of Heine's function or spheroidal wave functions. Two special cases (i) a sphere and (ii) a spheroid of small ellipticity, are studied in detail.

Non-disruptive and null-deflection mass flow measurement by a pressure compensation technique

Accurate mass flow measurement is very important in various monitoring and control applications. This paper proposes a novel method of fluid flow measurement by compensating the pressure drop across the ends of measuring unit using a compensating pump. The pressure drop due to the flow is balanced by a feedback control loop. This is a null-deflection type of measurement. As the insertion of such a measuring unit does not affect the functioning of the systems, this is also a non-disruptive flow measurement method. The implementation and design of such a unit are discussed. The system is modeled and simulated using the bond graph technique and it is experimentally validated. (C) 2009 Elsevier Ltd. All rights reserved.

The Evolution of Tribomaterial During Sliding: A Brief Introduction

This brief introductory article summarizes key findings from experiments and from computer simulations concerning the dramatic changes that commonly occur adjacent to sliding interfaces. We conclude that a wide range of observed features depends on a few basic processes (plastic deformation, interactions with the environment (including the counterface) and mechanical mixing) and that sliding leads to flow patterns similar to those expected in fluid flow.

Mechanism of glycine uptake by the silk glands of the silkworm Bombyx mori L.

The transport of glycine in vitro into the silk glands of the silkworm has been studied. Glycine accumulates inside the tissue to a concentration higher than that present outside, indicating an active transport mechanism. The kinetics of uptake show a biphasic curve and two apparent Km values for accumulation, 0.33 mM and 5.00 mM. The effect of inhibitors on the energy metabolism of glycine transport is inconclusive. Exchange studies indicate the existence of two pools inside the gland, one that is easily removed by exchange and osmotic shock, and the other which is not. The results obtained conform with the carrier model of Britten and McClure concerning the amino-acid pool in E. coli.

Transport and magnetic properties of conducting polyaniline doped with BX3 (X=F, Cl, and Br)

We report transport and magnetic properties of a different class of highly conducting polyaniline, doped with boron trihalides BX3 (X=F, Cl, and Br). In order to understand the transport mechanism we analyze the temperature dependence of resistivity of a large number of samples, made by pelletizing doped polyaniline powder and by doping films of polyaniline. We find that the charge transport in this class of conducting polyaniline is driven by the charging-energy limited transport of charge carriers, in contrast to the quasi-one-dimensional variable range hopping conduction prevalent in conventional proton-doped polyaniline samples. Magnetic susceptibility provides further insight into the unusually high intrinsic conductivity behavior.