Numerical modelling of foam Couette flows

A numerical computation based on a tensorial visco-elasto-plastic model based on continuous mechanics is compared to experimental measurements on liquid foams for a bidimensional Couette flow between two glass plates, both in stationary and transient cases. The main features of the model are elasticity up to a plastic yield stress, and viscoelasticity above it. The effect of the friction of the plates is taken into account. The numerical modelling is based on a small set of standard material parameters that are fully characterised. Shear localisation as well as acute transient observations are reproduced and agree with experimental measurements. The plasticity appears to be the fundamental mechanism of the localisation of the flow. Finally, the present approach could be extended from liquid foams to similar materials such as emulsions, colloids or wet granular materials, that exhibit localisation.     

An elasto-visco-plastic model for immortal foams or emulsions

A variety of complex fluids consists in soft, round objects (foams, emulsions, assemblies of copolymer micelles or of multilamellar vesicles--also known as onions). Their dense packing induces a slight deviation from their prefered circular or spherical shape. As a frustrated assembly of interacting bodies, such a material evolves from one conformation to another through a succession of discrete, topological events driven by finite external forces. As a result, the material exhibits a finite yield threshold. The individual objects usually evolve spontaneously (colloidal diffusion, object coalescence, molecular diffusion), and the material properties under low or vanishing stress may alter with time, a phenomenon known as aging. We neglect such effects to address the simpler behaviour of (uncommon) immortal fluids: we construct a minimal, fully tensorial, rheological model, equivalent to the (scalar) Bingham model. Importantly, the model consistently describes the ability of such soft materials to deform substantially in the elastic regime (be it compressible or not) before they undergo (incompressible) plastic creep--or viscous flow under even higher stresses.     

Soft Dynamics simulation. 2. Elastic spheres undergoing a T<Subscript>1</Subscript> process in a viscous fluid

Robust empirical constitutive laws for granular materials in air or in a viscous fluid have been expressed in terms of timescales based on the dynamics of a single particle. However, some behaviours such as viscosity bifurcation or shear localization, observed also in foams, emulsions, and block copolymer cubic phases, seem to involve other micro-timescales which may be related to the dynamics of local particle reorganizations. In the present work, we consider a T₁ process as an example of a rearrangement. Using the Soft Dynamics simulation method introduced in the first paper of this series, we describe theoretically and numerically the motion of four elastic spheres in a viscous fluid. Hydrodynamic interactions are described at the level of lubrication (Poiseuille squeezing and Couette shear flow) and the elastic deflection of the particle surface is modeled as Hertzian. The duration of the simulated T₁ process can vary substantially as a consequence of minute changes in the initial separations, consistently with predictions. For the first time, a collective behaviour is thus found to depend on a parameter other than the typical volume fraction of particles.     

On Sound Propagation in Dispersive Media

Sound propagation in monodisperse emulsions with arbitrary volume concentrations is studied theoretically using a cell model. It is assumed that emulsion cells are bounded by thin and imponderable rigid shells allowing realization of the minimum energy dissipation principle with viscous acoustic losses. Solutions covering many particular cases and wide parameter and variable ranges have been obtained. These solutions are suitable for studying acoustic properties of emulsions and suspensions, marine sediments, fogs, and smoke, as well as elastoviscous materials with solid or liquid inclusions, etc. Sound propagation and absorption in emulsions and suspensions are considered in more detail. The experimental data in the literature is compared.     

Steady drainage in emulsions: corrections for surface Plateau borders and a model for high aqueous volume fraction

We compare extensive experimental results for the gravity-driven steady drainage of oil-in-water emulsions with two theoretical predictions, both based on the assumption of Poiseuille flow. The first is from standard foam drainage theory, applicable at low aqueous volume fractions, for which a correction is derived to account for the effects of the confinement of the emulsion. The second arises from considering the permeability of a model porous medium consisting of solid sphere packings, applicable at higher aqueous volume fractions. We find quantitative agreement between experiment and the foam drainage theory at low aqueous volume fractions. At higher aqueous volume fractions, the reduced flow rate calculated from the permeability theory approaches the master curve of the experimental data. Our experimental data demonstrates the analogy between the problem of electrical flow and liquid flow through foams and emulsions.     

Theory of retardation of magnetic domain walls in rhombic magnetic materials

We calculate the retardation of a magnetic soliton describing a magnetic domain wall by using the generalized phenomenological theory of relaxation. We show that in this theory, based on the real dynamical symmetry of magnetic materials, the dissipation function has a different structure for high and low wall velocities. Finally, we calculate the viscous force of the wall in the Walker model and show that certain features, not discussed in the literature, emerge even when the generalized theory is applied to this simple model. In particular, the dependence of the viscous friction force on the wall velocity may be highly nonlinear and regions of unstable motion may appear. Zh. éksp. Teor. Fiz. 111, 158–173 (January 1997)     

Viscosity minimum in bimodal concentrated suspensions under shear

We study a model of concentrated suspensions under shear in two dimensions. Interactions between suspended particles are dominated by direct-contact viscoelastic forces and the particles are neutrally bouyant. The bimodal suspensions consist of a variable proportion between large and small droplets, with a fixed global suspended fraction. Going beyond the assumptions of the classical theory of Farris (R.J. Farris, Trans. Soc. Rheol. 12, 281 (1968)), we discuss a shear viscosity minimum, as a function of the small-to-large-particle ratio, in shear geometries imposed by external body forces and boundaries. Within a linear-response scheme, we find the dependence of the viscosity minimum on the imposed shear and the microscopic drop friction parameters. We also discuss the viscosity minimum under dynamically imposed shear applied by boundaries. We find a reduction of macroscopic viscosity with the increase of the microscopic friction parameters that is understood using a simple two-drop model. Our simulation results are qualitatively consistent with recent experiments in concentrated bimodal emulsions with a highly viscous or rigid suspended component. Received 28 June 2002 RID="a" ID="a"e-mail: ernesto@pion.ivic.ve     

The foam/emulsion analogy in structure and drainage

The often quoted analogy between foams and emulsions is experimentally tested by studying properties after settling and under forced drainage of oil-in-water emulsions of drop size similar as for bubbles generally used in foam experiments. Observations with regard to structure, water fraction and drainage wave properties confirm the expected similarity in the low flow rate range. However, while for foams a convective circulation on the scale of the container sets in for values of water fraction exceeding about 0.2, no such convection is found in emulsions. Here instabilities are only encountered at water fractions of about 0.4, close to the void fraction of random packings of spheres. These take on the form of descending pulses of increased water fraction and lead to the transition from a frozen to a locally agitated structure.Received: 12 December 2003, Published online: 24 August 2004PACS: 82.70.-y Disperse systems; complex fluids - 47.20.-k Hydrodynamic stability - 47.55.Dz Drops and bubbles     

New slip regimes and the shape of dewetting thin liquid films

We compare the flow behavior of liquid polymer films on silicon wafers coated with either octadecyl-(OTS) or dodecyltrichlorosilane (DTS). Our experiments show that dewetting on DTS is significantly faster than on OTS. We argue that this is tied to the difference in the solid/liquid friction. As the film dewets, the profile of the rim advancing into the undisturbed film is monotonically decaying on DTS but has an oscillatory structure on OTS. For the first time we can describe this transition in terms of a lubrication model with a Navier-slip condition for the flow of a viscous Newtonian liquid.     

Investigations of Viscous Dissipation in Stagnation Point Flow Past a Stretchable Riga Wall: Modern Analysis of Heat Transport

This article scrutinizes the features of viscous dissipation in the stagnation point flow past through a linearly stretched Riga wall by implementing Cattaneo-Christov heat flux model. Viscous dissipation is carried out in Cattaneo-Christov diffusion analysis for the first time in this letter. As a result of Cattaneo-Christov model, some extra terms of viscous dissipation are appeared in the energy equation. These extra terms of viscous dissipation are missing in the literature. On the utilization of suitable transformations, the equations governing the problem are reduced under the boundary layer approximation into the non-linear and dimensionless ordinary differential equations. Convergent approach is utilized to solve the dimensionless governing equations. The solution thus acquired is used to highlight the effects of emerging parameters on velocity distribution and fluid's temperature through the graphs. Features of the drag force (or skin friction co-efficient) are graphically interpreted. It is noticed that the presence of modified Hartman number helps to reduce the fluid's temperature but enhances the velocity profile. Further an enlargement in the value of thermal time relaxation parameter helps to decrease the temperature distribution.     

Investigation of Turbulent Transition in Plane Couette Flows Using Energy Gradient Method

The energy gradient method has been proposed with the aim of better understanding the mechanism of flow transition from laminar flow to turbulent flow. In this method, it is demonstrated that the transition to turbulence depends on the relative magnitudes of the transverse gradient of the total mechanical energy which amplifies the disturbance and the energy loss from viscous friction which damps the disturbance, for given imposed disturbance. For a given flow geometry and fluid properties, when the maximum of the function $K$ (a function standing for the ratio of the gradient of total mechanical energy in the transverse direction to the rate of energy loss due to viscous friction in the streamwise direction) in the flow field is larger than a certain critical value, it is expected that instability would occur for some initial disturbances. In this paper, using the energy gradient analysis, the equation for calculating the energy gradient function $K$ for plane Couette flow is derived. The result indicates that $K$ reaches the maximum at the moving walls. Thus, the fluid layer near the moving wall is the most dangerous position to generate initial oscillation at sufficient high $\operatorname{Re}$ for given same level of normalized perturbation in the domain. The critical value of $K$ at turbulent transition, which is observed from experiments, is about 370 for plane Couette flow when two walls move in opposite directions (anti-symmetry). This value is about the same as that for plane Poiseuille flow and pipe Poiseuille flow (385-389). Therefore, it is concluded that the critical value of $K$ at turbulent transition is about 370-389 for wall-bounded parallel shear flows which include both pressure (symmetrical case) and shear driven flows (anti-symmetrical case).     

Sound absorption of cellular metals with semiopen cells

A combined experimental and theoretical study is presented for the feasibility of using aluminum foams with semiopen cells for sound-absorption applications. The foams are processed via negative-pressure infiltration, using a preform consisting of water-soluble spherical particles. An analytical model is developed to quantify the dependence of pore connectivity on processing parameters, including infiltration pressure, particle size, wetting angle, and surface tension of molten alloy. Normal sound-absorption coefficient and static flow resistance are measured for samples having different porosity, pore size, and pore opening. A theory is developed for idealized semiopen metallic foams, with a regular hexagonal hollow prism having one circular aperture on each of its eight surfaces as the unit cell. The theory is built upon the acoustic impedance of the circular apertures (orifices) and cylindrical cavities due to viscous effects, and the principle of electroacoustic analogy. The predicted sound-absorption coefficients are compared with those measured. To help select processing parameters for producing semiopen metallic foams with desirable sound-absorbing properties, emphasis is placed on revealing the correlation between sound absorption and morphological parameters such as pore size, pore opening, and porosity.     

Non-linear oscillatory rheological properties of a generic continuum foam model: Comparison with experiments and shear-banding predictions

The occurrence of shear bands in a complex fluid is generally understood as resulting from a structural evolution of the material under shear, which leads (from a theoretical perspective) to a non-monotonic stationary flow curve related to the coexistence of different states of the material under shear. In this paper we present a scenario for shear-banding in a particular class of complex fluids, namely foams and concentrated emulsions, which differs from other scenarios in two important ways. First, the appearance of shear bands is shown to be possible both without any intrinsic physical evolution of the material (e.g. via a parameter coupled to the flow such as concentration or entanglements) and without any finite critical shear rate below which the flow does not remain stationary and homogeneous. Secondly, the appearance of shear bands depends on the initial conditions, i.e. the preparation of the material. In other words, it is history dependent. This behaviour relies on the tensorial character of the underlying model (2D or 3D) and is triggered by an initially inhomogeneous strain distribution in the material. The shear rate displays a discontinuity at the band boundary whose amplitude is history dependent and thus depends on the sample preparation.     

Negative differential viscosity in magnetic suspensions

Recent experiments have shown that the dependence of the macroscopic viscous stress on the mean velocity gradient during the Couette flow of concentrated magnetic suspensions in an external magnetic field is N-shaped. As the field strength is decreased, the amplitude of the N-shaped curve decreases and in the absence of the field, the stress monotonically increases with the shear velocity. A model is proposed to explain the shape of the rheological curve. The model assumes that the magnetic field initiates the formation of dense aggregates in a suspension, which connect the opposite walls of a measurement cell. In the Couette flow, the friction of aggregates on the cell walls causes their deviation from the applied magnetic field through an angle determined by the velocity of the relative motion of the walls. For large enough velocities, the aggregates are detached from the wall and are destroyed by viscous forces. It is shown that the friction of aggregates on cell walls results in the initial increasing and decreasing part of the N-shaped rheogram, while the flow after the detachment of aggregates corresponds to its right increasing part.     

计入固液界面作用的润滑热力学模型与分析

摩擦与润滑过程是典型的能量耗散过程, 在机理上与非平衡热力学中的熵增、耗散结构等理论颇有相似之处. 通过热力学分析可以对一些典型的摩擦磨损过程做出合理的机理揭示与推测. 本文利用热力学理论对典型的润滑过程进行了建模分析. 采用分离压模型表征和计入了微尺度下的固液界面作用, 揭示分析了润滑热力学模型与润滑状态Stribeck曲线的联系. 从分析计算结果来看, 润滑Stribeck曲线的摩擦系数最低点与系统热力学上的熵增率最低点具有相当好的对应关系, 而润滑状态从弹流润滑向薄膜润滑的转变过程, 可以用耗散结构理论加以机理解释. 文中的热力学模型和方法能够有效地体现出润滑过程中多物理要素跨尺度非线性耦合的作用, 对实际工程与实验有着重要的指导作用.     

Shear viscosity of glass-forming melts in the liquid-glass transition region

A new approach to interpreting the hole-activation model of a viscous flow of glass-forming liquids is proposed. This model underlies the development of the concept on the exponential temperature dependence of the free energy of activation of a flow within the range of the liquid-glass transition in complete agreement with available experimental data. The “formation of a fluctuation hole” in high-heat glass-forming melts is considered as a small-scale low-activation local deformation of a structural network, i.e., the quasi-lattice necessary for the switching of the valence bond, which is the main elementary event of viscous flow of glasses and their melts. In this sense, the hole formation is a conditioned process. A drastic increase in the activation free energy of viscous flow in the liquid-glass transition region is explained by a structural transformation that is reduced to a limiting local elastic deformation of the structural network, which, in turn, originates from the excitation (critical displacement) of a bridging atom like the oxygen atom in the Si-O-Si bridge. At elevated temperatures, as a rule, a necessary amount of excited bridging atoms (locally deformed regions of the structural network) always exists, and the activation free energy of viscous flow is almost independent of temperature. The hole-activation model is closely connected with a number of well-known models describing the viscous flow of glass-forming liquids (the Avramov-Milchev, Nemilov, Ojovan, and other models).     

da Vinci fluids,catch-up dynamics and dense granular flow

We introduce and study a da Vinci fluid, a fluid whose dissipation is dominated by solid friction. We analyse the flow rheology of a discrete model and then coarse-grain it to the continuum. We find that the model gives rise to behaviour that is characteristic of dense granular fluids. In particular, it leads to plug flow. We analyse the nucleation mechanism of plugs and their development. We find that plug boundaries generically expand and we calculate the growth rate of plug regions. In systems whose internal effective dynamic and static friction coefficients are relatively uniform we find that the linear size of plug regions grows as (time)^1/3 . The suitability of the model to granular materials is discussed.     

A numerical analysis on flow field and heat transfer by interaction between a pair of vortices in rectangular channel flow

This paper is a numerical study of the effect of flow field and heat transfer created by interactions between a pair of vortices generated by a vortex generator in a rectangular channel flow. In order to analyze the vortices produced by the vortex generator, the pseudo-compressibility method is introduced into the Navier–Strokes (NS) equation of a three-dimensional unsteady, incompressible viscous flow. A two-layer k–ε turbulence model is used on the flat plate three-dimensional turbulence boundary to predict the turbulence characteristics of the vortices. The computational results accurately predict the vortex characteristics, which are related to Reynolds stress, turbulent kinetic energy, and flow field. Also, in the prediction of thermal boundary layers, skin friction characteristics, and heat transfers, the present results are reasonably close to the experimental results obtained by other researchers.     

Rapid Plateau border size variations expected in three simple experiments on 2D liquid foams<Superscript>⋆</Superscript>

Up to a global scaling, the geometry of foams squeezed between two solid plates (2D GG foams) essentially depends on two independent parameters: the liquid volume fraction and the degree of squeezing (bubble thickness to diameter ratio). We describe it in two main asymptotic regimes: fully dry floor tiles, where the Plateau border radius is smaller than the distance between the solid plates, and dry pancakes, where it is larger. We predict a rapid variation of the Plateau border radius in one part of the pancake regime, namely when the Plateau border radius is larger than the inter-plate distance but smaller than the geometric mean of that distance and the bubble perimeter. This rapid variation is not related to any topological change in the foam: in all the regimes we consider, the bubbles remain in mutual lateral contact through films located at mid-height between both plates. We provide asymptotic predictions in different types of experiments on such 2D GG foams: when foam is being progressively dried or wetted, when it is being squeezed further or stretched, when it coarsens through film breakage or through inter-bubble gas diffusion. Our analysis is restricted to configurations close to equilibrium, as we do not include stresses resulting from bulk viscous flow or from non-homogeneous surfactant concentrations. We also assume that the inter-plate distance is sufficiently small for gravity to be negligible. The present work does not provide a method for measuring small Plateau border radii experimentally, but it indicates that large (and easily observable) Plateau borders should appear or disappear rather suddenly in some types of experiments with small inter-plate gaps. It also gives expected orders of magnitude that should be helpful for designing experiments on 2D GG foams.     

Hydrodynamic instability of premixed flame propagating in narrow planar channel in the presence of gas flow

The paper analyses the hydrodynamic instability of a flame propagating in the space between two parallel plates in the presence of gas flow. The linear analysis was performed in the framework of a two-dimensional model that describes the averaged gas flow in the space between the plates and the perturbations development of two-dimensional combustion wave. The model includes the parametric dependences of the flame front propagation velocity on its local curvature and on the combustible gas velocity averaged along the height of the channel. It is assumed that the viscous gas flow changes the surface area of the flame front and thereby affects the propagation velocity of the two-dimensional combustion wave. In the absence of the influence of the channel walls on the gas flow, the model transforms into the Darrieus–Landau model of flame hydrodynamic instability. The dependences of the instability growth rate on the wave vector of disturbances, the velocity of the unperturbed gas flow, the viscous friction coefficients and other parameters of the problem are obtained. It is shown that the viscous gas flow in the channel can lead, in some cases, to a significant increase in instability compared with a flame propagating in free space. In particular, the instability increment depends on the direction of the gas flow with respect direction of the flame propagation. In the case when the gas flow moves in the opposite direction to the direction of the flame propagation, the pulsating instability can appear.