Instabilités absolue et convective de l'écoulement de Poiseuille en conduite viscoélastiqueAbsolute and convective instabilities of Poiseuille flow in viscoelastic tube

The temporal stability of Poiseuille flow in viscoelastic tube has been studied by the author in a set of papers. In this paper the spatio-temporal stability of the foregoing system is examined. The method of orthnormalization is used in the numerical computation leading to an efficient algorithm to solve for an infinitesimal axisymmetric and nonaxisymmetric disturbances. Then a new method for eigenmodes search technique is founded which eliminates the need for the recurrence formula, developed previously by the author, in order to find the eigenmodes of the system. This new technique leads to more efficient and tractable algorithm to inspect the dispersion relation of the system in the complex wave number and frequency planes. The nature of the instabilities depends on the boundary conditions imposed at the outside tube surface. In the range of physical parameters taken into account in this paper, it is found that for a noncollapsible tube one of the two axisymmetric unstable modes represents an absolute instability. For collapsible tube, the absolute instability is induced by a set of nonaxisymmetric modes. The ratio of the frequency of these absolute unstable modes in their cusp points are in good agreement with the experimental finding. The absolute instability is introduced by solid-based mode. To cite this article: M. Hamadiche, C. R. Mecanique 330 (2002) 769–775.     

Simulations of magnetorheological suspensions in Poiseuille flow

Particle-level simulations are conducted to study magnetorheological fluids in plane Poiseuille flow. The importance of the boundary conditions for the particles at the channel walls is examined by considering two extreme cases: no friction and infinite coefficient of friction. The inclusion of friction produces Bingham fluid behavior, as commonly observed experimentally for MR suspensions. Lamellar structures, similar to those reported for electrorheological fluids in shear flow, are observed in the post-yield region for both particle boundary conditions. The formation of these lamellae is accompanied by an increase in the bulk fluid velocity. The slip boundary condition produces higher fluid velocities and thicker lamellar structures.     

On the non-linear instability of fiber suspensions in a Poiseuille flow

The non-linear instability characteristics of fiber suspensions in a plane Poiseuille flow are investigated. The evolution equation of the perturbation amplitude analogous to Landau equation is formulated and solved numerically for different fiber parameters. It is found that the equilibrium amplitude decreases with the increase of the fiber aspect ratio and volume fraction, i.e. the addition of the fibers reduces the amplitude of the perturbation, and leads to the reduction of the flow instability. This phenomenon becomes significant for large volume concentration and aspect ratio. The mechanism of the reduction of the flow instability is also analyzed by taking into account of the modification of the mean flow and the energy transfer from the mean flow to the perturbation flow.     

Rheo-optics of hydroxypropylcellulose solutions in Poiseuille flow

The behavior of an aqueous solution of hydroxypropylcellulose in the liquid crystalline phase is investigated when it is flowing in a rectangular channel. Rheological characterization shows that the viscosity vs. shear rate curve follows the typical three region pattern, with the intermediate plateau of region II extending over a relative large range of shear rates (more than one decade). Two complementary rheo-optical determinations are performed. Velocity profiles across the channel thickness are measured by a hydrogen bubble visualization technique. Texture evolution is monitored by in situ optical microscopy. Accurate focusing inside the sample thickness allows observation in real time of the texture at various shear rates, as generated in the Poiseuille type of flow in the channel. It is shown that the velocity profiles can be accurately predicted by assuming that the flow in the channel is purely viscous, and using only the viscosity data described above. It is also shown that the morphology of the texture generated inside the flowing system is a function of the local shear rate. In particular, an elongated structure is observed when the shear rate exceeds the critical value corresponding to the onset of region II in the viscosity curve.     

The linear stability of plane Poiseuille flow under unsteady distortion

This paper investigates the linear stability behaviour of plane Poiseuille flow underunsteady distortion by multiscale perturbation method and discusses further the problemproposed by paper[1].The results show that in the initial period of disturbancedevelopment,the distortion profiles presented by paper[1]will make the disturbances growup,thus augmenting the possibility of instability.     

Purely elastic interfacial instabilities in superposed flow of polymeric fluids

Purely elastic interfacial stability of superposed plane Poiseuille flow of polymeric liquids has been investigated utilizing both asymptotic and numerical techniques. It is shown that these instabilities are caused by an unfavorable jump in the first normal stress difference across the fluid interface. To determine the significance of these instabilities in finite experimental geometries, a comparison between the maximum growth rates of purely elastic instabilities with instabilities driven primarily by a viscosity or a combined viscosity and elasticity difference is made. Based on this comparison, it is shown that purely elastic interfacial instabilities can play a major role in superposed flow of polymeric liquids in finite experimental geometries.     

Convective and absolute instabilities in non-Boussinesq mixed convection

The problem of non-Boussinesq mixed convection in a vertical channel formed by two differentially heated infinite plates is investigated and the complete convective/absolute instability boundary is computed for a wide range of physical parameters. A physical insight into the mechanisms causing instabilities is given. In particular, it is shown that the appearance of absolute instability is always dictated by a flow reversal within a channel; however, existence of the flow reversal does not exclude the possibility of convective instability. It is also shown that fluid’s non-linear transport property variations have a dramatic effect on the structure and complexity of spatio-temporal instabilities of the co-existing buoyancy and shear modes as the temperature difference across the channel increases. The validity of the stability results obtained using the procedure described in Suslov (J Comp Phys 212, 188–217, 2006) is assessed using the method of steepest descent. This work was partially supported by a computing grant from the Australian Partnership for Advanced Computing, 2000–2003.     

On the accuracy of the calculation of transient growth in plane Poiseuille flow

This paper addresses the accuracy of numerical methods to compute the transient energy growth of plane Poiseuille flow. We show that using the Chebyshev collocation method to discretize the linearized governing equations in the wall‐normal direction can introduce numerical problems when computing the energy evolution of the flow. We demonstrate that spurious eigenmodes of the discretized linear operator and numerical integration errors are the possible sources of the numerical problems, and we also show that spurious eigenvalues with negative real parts of large magnitude can affect the calculation of energy growth. These difficulties can be avoided by using a spectral Galerkin method where the basis functions satisfy the boundary conditions. Copyright © 2014 John Wiley & Sons, Ltd.     

Thermal instability of a viscoelastic fluid in a fluid-porous system with a plane Poiseuille flow

The thermal convection of a Jeffreys fluid subjected to a plane Poiseuille flow in a fluid-porous system composed of a fluid layer and a porous layer is studied in the paper. A linear stability analysis and a Chebyshev τ-QZ algorithm are employed to solve the thermal mixed convection. Unlike the case in a single layer, the neutral curves of the two-layer system may be bi-modal in the proper depth ratio of the two layers. We find that the longitudinal rolls(LRs) only depend on the depth ratio. Wi...     

Two-phase flow instabilities in a horizontal single boiling channel

A study of the stability of an electrically heated single channel, forced convection horizontal system was conducted by using Freon-11 as the test fluid. Two major modes of oscillations, namely, density-wave type (high frequency) and pressure-drop type (low frequency) oscillations have been observed. The steady-state operating characteristics and stable and unstable regions are determined as a function of heat flux, exit orifice diameter and mass flow rate. Different modes of oscillations and their characteristics have been investigated. The effect of the exit restriction on the system stability has also been studied.A mathematical model has been developed to predict the transient behavior of boiling two-phase systems. The model is based on homogenous flow assumption and thermodynamic equilibrium between the liquid and vapor phases. The transient characteristics of boiling two-phase flow horizontal system are obtained for various heat inputs, flow rates and exit orifice diameters by perturbing the governing equations around a steady state. Theoretical and experimental results have been compared.     

Effect of random longitudinal vibrations on the Poiseuille flow of a complex liquid

In this work, the rectilinear Poiseuille flow of a complex liquid flowing in a vibrating pipe is analyzed. The pipe wall performs oscillations of small amplitude that can be adequately represented by a weakly stochastic process, for which a quasi-static perturbation solution scheme is suggested. The flow is analyzed using the Bautista–Manero–Puig constitutive equation, consisting on the upper-convected Maxwell equation coupled to a kinetic equation to account for the breakdown and reformation of the fluid structure. A drastic enhancement of the volumetric flow is predicted in the region where the fluid experiences pronounced shear-thinning. Finally, flow enhancement is predicted using experimental data reported elsewhere for wormlike micellar solutions of cetyl trimethyl ammonium tosilate.     

Inlet instabilities in the capillary flow of polyethylene melts

Inlet instabilities in the capillary flow of polyethylene melts were studied in this work. Extrudate distortions in branched polyethylenes, produced by unstable upstream flow, were found to be accompanied by pressure oscillations that do not have their origin in the slip phenomenon, but on polymer compressibility. The absence of slip was clearly evidenced in the experiments, and the differences between pressure oscillations occurring in linear and branched polymers are shown.Pressure oscillations in the capillary flow of branched polyethylenes were found to be made up of at least two components of different frequency and amplitude. These two components were identified with different bulk defects appearing in the extrudates. Information about the dynamics of vortices upstream of the contraction and extrudate distortions is obtained from the analysis of pressure oscillations.The influence of capillary entrance angle on flow curves was also investigated. From the results, it is concluded that the extensional component of the flow in the contraction is the main factor responsible for the slope change usually found in the log-log flow curves of both linear and branched polyethylenes.     

Stokes eigenfunctions and Galerkin projection of the disturbance equations in plane Poiseuille flow: a systematic analytical approach

In the present paper, the L ²-normalized Stokes eigenfunctions for plane Poiseuille flow, which form an orthonormal functional basis for the space of disturbances, are written in a general exponential form. Then, the evolution equations for the disturbances are Galerkin-projected on the considered basis functions, and all the terms of the resulting dynamical system are expressed systematically in analytical form. Finally, a numerical example is given in which the proposed basis functions are used for the simulation of the time evolution of the critical disturbance predicted by the energetic stability theory.     

On the use of genetic algorithm for finding the neutral instability curve in plane Poiseuille flow

A novel method based on genetic algorithm (GA) is proposed, to the best of our knowledge for the first time, for finding the neutral instability curve of the Orr-Sommerfeld equation in (nearly) parallel flows. New concepts such as “proximity of parents” and “gender discrimination” are added to the conventional GA in order for this algorithm to find the neutral instability curve. Certain GA operators such as “crossover” and “mutation” will also be modified in such a way that this algorithm can meet this purpose. To check the applicability of the modified genetic algorithm (MGA) developed in this work in finding the neutral instability curve, the case of plane Poiseuille flow will be used as a benchmark. It will be shown that the modified genetic algorithm developed in this work is well capable of determining the neutral instability curve for this particular flow geometry.     

The effect of internal surface modification on flow instabilities in forced convection boiling in a horizontal tube

This paper presents the experimental result of a study on the effects of heat transfer enhancement on two-phase flow instabilities in a horizontal in-tube flow boiling system. Five different heat transfer surface configurations and five different inlet temperatures are used to observe the effect of heat transfer enhancement and inlet subcooling. All experiments are carried out at constant heat input, system pressure and exit restriction. Dynamic instabilities, namely pressure-drop type, density-wave type and thermal oscillations are found to occur for all the investigated temperatures and enhancement configurations, and the boundaries for the appearance of these oscillations are found. The effect of the enhancement configurations on the characteristics of the boiling flow dynamic instabilities is studied in detail. The comparison between the bare tube and the enhanced tube configurations are made on the basis of boiling flow instabilities. Differences among the enhanced configurations are also determined to observe which of them is the most stable and unstable one. The amplitudes and periods of pressure-drop type oscillations and density-wave type oscillations for tubes with enhanced surfaces are found to be higher than those of the bare tube. The bare tube is found to be the most stable configuration, while tube with internal springs having bigger pitch is found to be the most unstable one among the tested tubes. It is found that system stability increases with decreasing equivalent diameter for the same type heater tube configurations; however, on the basis of effective diameter there is no single result such as stability increase/decrease with increasing/decreasing effective diameter.     

On the modeling of fluidity loss phenomena in Couette and Poiseuille flows of elastic liquids

Some effects of the possible relaxation transition from viscoelastic liquid state to highly elastic solid state were theoretically and numerically investigated in the shear situations, within the approach proposed in papers [1, 2, 5, 16]. It was found that for a single Maxwellian model the constitutive equations developed in [1, 2, 5] are not valid at elevated shear stresses. Some new aspects of the possible rheological behavior of elastic liquids in subcritical (before transition) and supercritical (after transition) regimes were demonstrated. The mechanism of fluidity loss studied in this paper could serve as a possible trigger mechanism for the melt flow instabilities.     

Poiseuille flow of Leslie–Ericksen discotic liquid crystals: solution multiplicity, multistability, and non-Newtonian rheology

Computational modeling of the steady capillary Poiseuille flow of flow-aligning discotic nematic liquid crystals (DNLCs) using the Leslie–Ericksen (LE) equations predicts solution multiplicity and multistability. The phenomena are independent of boundary conditions. The steady state solutions are classified into: (a) primary, (b) secondary, and (c) hybrid. Primary solutions exist for all orientation boundary conditions and all flow rates, and are characterized by a flow-alignment angle that is closest to the anchoring angle at the bounding surface. Secondary solutions exist for all orientation boundary conditions and flow rates above a certain critical value. The secondary solutions are characterized by a flow-alignment angle which can be either the nearest neighbor below the primary solution or any multiple of π above. Hybrid solutions interpolate between the primary and the nearest secondary solutions, and hence exhibit two alignment angles. All solutions are stable to planar finite amplitude perturbations. Hybrid solutions are unstable to front propagation and lead to primary or secondary solutions. The non-Newtonian rheology of the primary and secondary solutions is characterized by non-classical shear thinning and thickening apparent viscosity behavior. Well-aligned monodomains can lead to shear thickening, thinning, or a sequence of both. The degree of rheological uncertainty is present for planar and homeotropic anchoring conditions. The non-Newtonian rheology of non-aligned samples leads to shear thinning and lack the uncertainty of well-aligned samples, since the apparent viscosity becomes insensitive to orientation.     

The Convective Boundary-Layer Flow on a Reacting Surface in a Porous Medium

The convective boundary-layer flow on an impermeable vertical surface in a fluid-saturated porous medium is considered where the flow results from the heat released by an exothermic catalytic reaction on the surface converting a reactive component within the convective fluid to an inert product. The reaction is modelled by first-order kinetics with an Arrhenius temperature dependence. Numerical solutions of the governing equations are obtained for a range of parameter values. These show, for large activation energies, that localized rapid changes in wall temperature and localized high reaction rates occur a little way from the leading edge. Asymptotic expansions, valid at large distances from the leading edge, are derived, the form that these expansions take is qualitatively different depending on whether or not reactant consumption is included in the model.     

The hydrodynamic force and torque on a bounded sphere in Poiseuille flow

The lattice Boltzmann (LB) method is used to study the hydrodynamic force and torque acting on a sphere held stationary between parallel plates in pressure‐driven flow. This and associated flow configurations are explored in this paper. LB results are in excellent agreement with existing theory and numerical results for simple pressure‐driven flow between parallel plates, for flow through a periodic medium of spheres [Zick AA, Homsy GM. Stokes flow through periodic arrays of spheres. Journal of Fluid Mechanics 1982; 115: 13], and for the force and torque acting on a sphere held fixed at the quarter vertical position in a pressure‐driven flow between parallel plates. In the latter case, LB calculations reveal a screening effect caused by neighboring periodic images of the test sphere. It is shown that the test sphere is hydrodynamically decoupled from its periodic images when separated by approximately 30 sphere radii. LB results for force and torque as a function of sphere height and flow cell height are also reported. Copyright © 2001 John Wiley & Sons, Ltd.