Stability of the plane Couette flow of a disperse medium with a finite volume fraction of the particles

The linear hydrodynamic stability of the plane Couette flow of a suspension with a finite volume fraction of the particles is considered. The two-phase medium flow is described within the framework of the model of mutually penetrating continua which allows for the finiteness of the volume occupied by the particles. In the main flow the phase velocities are the same, while gravity is not taken into account. The stability of disperse flows with both uniform and nonuniform particle distributions is studied. The linearized system of the equations of suspension motion with the no-slip boundary conditions imposed on solid walls is reduced to the eigenvalue problem for an ordinary differential fourth-order equation in the stream function. The eigenvalues are sought using the orthogonolization method. The parametric investigation of the stability characteristics of the disperse flow is performed. It is shown that in the case of the uniform spatial distribution of the particles in the main flow, the presence of an admixture in the flow leads to a slight variation in the wave decay rates, while the flow remains stable for any permissible combinations of the dimensionless governing parameters. In the case of nonuniform distribution of inclusions the flow loses stability already for low Reynolds numbers on a wide range of the dimensionless governing parameters.     

About modelling the shape memory alloy behaviour based on the phase transformation surface identification under proportional loading and anisothermal conditions

Based on the experimental identification of the phase transformation surface, a modelling of shape memory alloys anisothermal behaviour is proposed. Within the framework of the thermodynamics of irreversible processes, two internal variables are chosen: the stress-induced martensite volume fraction and the self-accommodating martensite volume fraction. A special attention is paid to take into account the asymmetry between tension and compression behaviours.     

Contamination of a porous bed by a moving front with a finite porosity jump

The contamination of an initially porous bed by a moving low-concentrated suspension is studied. In the framework of a hyperbolic model taking into account large variations of porosity, several exact solutions are given to the problem with a finite porosity jump at the leading contamination front in the case of flows with plane and cylindrical waves. The fact that the particle front lags behind the carrier liquid front is explained. It is shown that the finiteness of the porosity jump at the contamination wave front causes the deceleration in the jump motion, which can be used to determine experimentally the governing parameters of the model. It is also shown that some solutions with several jumps may exist.     

Numerical simulations of particle migration in suspension flows: Frame-invariant formulation of curvature-induced migration

The objective of this work is to investigate what mechanisms should be employed to qualitatively/quantitatively predict particle migration in a suspension flow. Based on the diffusive flux model originally proposed by Phillips et al. [R.J. Phillips, R.C. Armstrong, R.A. Brown, A.L. Graham, A constitutive equation for concentrated suspensions that account for shear-induced particle migration, Phys. Fluids A 4 (1992) 30–40], we survey the accuracy of three models including original Phillips model (Model I), modified Phillips model with curvature-induced migration mechanism (Model II), and finally the modified Model II with volume-fraction-dependent parameters (Model III). The empirical parameters which appear in the three models are determined by fitting to independent concentric Couette experiments. The accuracy of three models in concentric Couette problem is comparable except that Model III shows more improved predictions near the inner cylinder. However, the predictions of the three models are entirely different on a qualitative level for parallel plate problems and the existence and direction of particle migration are severely model-dependent. Models II and III predict no migration or very slight migration at high volume fraction, which is in good agreement with the previous experiments, whereas Model I predicts inward migration. We show that Model III accurately predicts a solid-free region near the center at low volume fraction, which was experimentally observed.In addition to a survey of migration mechanisms, we developed a frame-invariant curvature-induced migration model applicable to multi-dimensional flows. A transient 2D mixed-order finite element method (FEM) code was implemented to compare the predictions of the three models in a 2D problem. In this work, we considered the eccentric Couette problem, which is often used as a benchmarking problem. Though there is not much difference among the three models, Model III predicts that the particle migration is slightly retarded at high shear rate regions.     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

Discontinuous shear thickening in the presence of polymers adsorbed on the surface of calcium carbonate particles

In the presence of dispersant molecules currently used in cement industry and based on polyethylene oxide (PEO), we found a strong discontinuous shear thickening (DST) at high volume fraction in suspensions of calcium carbonate particles. The transition was reversible and the critical shear rate and shear stress for which this instability appears are reported versus the volume fraction of particles. A model of repulsive forces between polymers, taking into account the thickness of the polymer layer and the density of adsorption on the surface of the particles, can explain the differences of critical stresses observed between these three dispersant molecules. In particular, it explains why a small polymer densely adsorbed can be more efficient to repel the transition at higher stress than a larger molecule less densely adsorbed. Above the transition, we find that the suspension presents a special kind of stick-slip instability with even the presence of a negative shear rate under constant applied stress. A model is proposed which well predicts this regime by taking into account both the inertia of the apparatus and the viscoelasticity of the suspension.     

Structure and stability of the interface between magnetic and conventional fluids. Model of a three-component medium

The structure of a flat interphase boundary between a magnetic suspension and a conventional immiscible fluid is investigated within the framework of the model of a three-component medium taking the dependence of the free energy of the system on the concentration gradients into account. It is shown that for certain values of the constitutive parameters the bulk magnetic particle concentration increases significantly inside the interfacial layer, i.e., the particles are significantly adsorbed on the interface. The dependence of the surface tension on the magnetic field strength is determined. It is shown that for certain problem parameters this dependence qualitatively corresponds to that obtained experimentally and described in the phenomenological theory developed by Golubyatnikov and Subkhankulov in 1986. In the case of strong particle adsorption the dependence of the surface tension on the magnetic particle concentration on the phase interface is significantly nonlinear. A refined model of the interface as a two-dimensional continuum with surface magnetization is constructed. Constitutive equations, conditions on the interface, and necessary stability conditions are obtained.     

Stability of equilibrium of a double-layer system with a deformable interface and a prescribed heat flux on the external boundaries

The stability of mechanical equilibrium of a system of two horizontal immiscible-liquid layers with similar densities is studied. The problem is solved for a prescribed heat flux on the external boundaries. Within the framework of a generalized Boussinesq approximation, which takes the interface deformation correctly into account, the onset of convection caused by heating the system from above or below is considered. Two long-wave instability modes attributable to the presence of the deformable interface and the given heat flux on the external boundaries are detected. The system response to monotonic and oscillatory disturbances with finite wavelengths is investigated. A complete stability map is constructed.     

Homogenization-based constitutive models for magnetorheological elastomers at finite strain

This paper proposes a new homogenization framework for magnetoelastic composites accounting for the effect of magnetic dipole interactions, as well as finite strains. In addition, it provides an application for magnetorheological elastomers via a “partial decoupling” approximation splitting the magnetoelastic energy into a purely mechanical component, together with a magnetostatic component evaluated in the deformed configuration of the composite, as estimated by means of the purely mechanical solution of the problem. It is argued that the resulting constitutive model for the material, which can account for the initial volume fraction, average shape, orientation and distribution of the magnetically anisotropic, non-spherical particles, should be quite accurate at least for perfectly aligned magnetic and mechanical loadings. The theory predicts the existence of certain “extra” stresses—arising in the composite beyond the purely mechanical and magnetic (Maxwell) stresses—which can be directly linked to deformation-induced changes in the microstructure. For the special case of isotropic distributions of magnetically isotropic, spherical particles, the extra stresses are due to changes in the particle two-point distribution function with the deformation, and are of order volume fraction squared, while the corresponding extra stresses for the case of aligned, ellipsoidal particles can be of order volume fraction, when changes are induced by the deformation in the orientation of the particles. The theory is capable of handling the strongly nonlinear effects associated with finite strains and magnetic saturation of the particles at sufficiently high deformations and magnetic fields, respectively.     

The bulk electroviscous effect

Electroviscous stresses arise as hydrodynamic flows disturb the ionic (Debye) clouds that screen charged surfaces in electrolyte solutions. The contribution thereof to the effective bulk viscosity (also known as the second or volume viscosity) of two-phase suspensions is quantified here. Specifically, the bulk viscosity of two model suspensions is calculated: (1) a dilute dispersion of rigid charged spherical particles immersed in a compressible electrolyte that undergoes uniform dilatation and (2) a dilute suspension of charged gas bubbles expanding uniformly in an incompressible electrolyte. In both cases, it is assumed that the fluid flow only slightly drives the Debye cloud out of equilibrium, which formally requires that the ratio of the ion diffusion to flow time scales—a Péclet number—is small. For a suspension of rigid particles, the electroviscous contribution to the effective bulk viscosity is proportional to the particle volume fraction and decreases monotonically as the ratio of the particle radius to the Debye length increases. Similar behavior is well known for the electroviscous contribution to the effective shear viscosity of a dilute hard-sphere suspension; a quantitative comparison between the bulk and shear viscosities is made. In contrast, the electroviscous contribution to the bulk viscosity of a dilute suspension of bubbles is independent of the bubble volume fraction and attains a finite value in the limit of vanishing Debye length.     

Electroelastic behavior of piezoelectric composites with coated reinforcements: Micromechanical approach and applications

In this work, a micromechanical model for the estimate of the electroelastic behavior of the piezoelectric composites with coated reinforcements is proposed. The piezoelectric coating is considered as a thin layer with active electroelastic properties different from those of the inclusion and the matrix. The micromechanical approach based on the Green’s functions technique and on the interfacial operators is designed for solving the electroelasticity inhomogeneous coated inclusion problem. The effective properties of a piezoelectric composite containing thinly coated inclusions are obtained through the Mori–Tanaka’s model. Numerical investigations into electroelastic moduli responsible for the electromechanical coupling are presented as functions of the volume fraction and characteristics of the coated inclusions. Comparisons with existing analytical and numerical results are presented for cylindrical and elliptic coated inclusions.     

Magneto-electro-elastic coated inclusion problem and its application to magnetic-piezoelectric composite materials

In this work, a micromechanical model for the estimate of the magneto-electro-elastic behavior of the magnetic-piezoelectric composites with coated reinforcements is proposed. The coating is considered as a thin layer with properties different from those of the inclusion and the matrix. The micromechanical approach based on the Green’s functions techniques and on the interfacial operators is designed for solving the magneto-electro-elastic inhomogeneous coated inclusion problem. The effective magneto-electro-elastic properties of the composite containing thinly coated inclusions are obtained through the Mori–Tanaka’s model. Numerical investigations into magneto-electro-elastic moduli responsible for the magneto-electric coupling are presented as functions of the volume fraction and characteristics of the coated inclusions. Comparisons with existing models are presented for various shape and orientation of the coated inclusions.     

A model for the shear viscosity of non-colloidal suspensions with Newtonian matrix fluids

We present a model for the shear viscosity of non-colloidal suspensions with Newtonian matrix fluids. The model is based on the original idea first presented by Brinkman (Applied Sci Research A1:27-34. 1947) for the viscous force exerted by a flowing fluid on a dense swarm of spherical particles. In particular, we consider an inertialess suspension in which the mean flow is driven by a pressure difference, and simultaneously, the suspension is subject to simple shear. Assuming steady state, incompressibility and taking into account a resistance force which is generated due to the presence of the particles in the flow, the three-dimensional governing equations for the mean flow around a single spherical particle are solved analytically. Self-consistency of the model provides a relationship between the resistance parameter and the volume fraction of the solid phase. A volume, or an ensemble, averaging of the total stress gives the bulk properties and an expression for the relative (bulk) viscosity of the suspension. The viscosity expression reduces to the Einstein limit for dilute suspensions and agrees well with empirical formulas from the literature in the semi-dilute and concentrated regimes. Since the model is based on a single particle and its average interaction with the other particles is isotropic, no normal stress differences can be predicted. A possible method of addressing this problem is provided in the paper.     

A physically-based constitutive model for anisotropic damage in rubber-toughened glassy polymers during finite deformation

The present work focuses on the development of a physically-based model for large deformation stress-strain response and anisotropic damage in rubber-toughened glassy polymers. The main features leading to a microstructural evolution (regarding cavitation, void aspect ratio, matrix plastic anisotropy and rubbery phase deformation) in rubber-toughened glassy polymers are introduced in the proposed constitutive model. The constitutive response of the glassy polymer matrix is modelled using the hyperelastic-viscoplastic model of [Boyce et al., 1988] and [Boyce et al., 2000]. The deformation mechanisms of the matrix material are accounted for by two resistances: an elastic-viscoplastic isotropic intermolecular resistance acting in parallel with a visco-hyperelastic anisotropic network resistance, each resistance being modified to account for damage effects by void growth with a variation of the void aspect ratio. The effective contribution of the hyperelastic particles to the overall composite behaviour is taken into account by treating the overall system in a composite scheme framework. The capabilities of the proposed constitutive model are checked by comparing experimental data with numerical simulations. The deformation behaviour of rubber-toughened poly(methyl methacrylate) was investigated experimentally in tension at a temperature of 80 °C and for different constant true strain rates monitored by a video-controlled technique. The reinforcing phase is of the soft core-hard shell type and its diameter is of the order of one hundred nanometers. The particle volume fraction was adjusted from 15% to 45% by increments of 5%. The stress-strain response and the inelastic volumetric strain are found to depend markedly on particle volume fraction. For a wide range of rubber volume fractions, the model simulations are in good agreement with the experimental results. Finally, a parametric analysis demonstrates the importance of accounting for void shape, matrix plastic anisotropy and rubber content.     

Effect of structural transformation and deformation nonlinearity on the stability of a shape memory alloy rod

Within the framework of a model of nonlinear deformations of shape memory alloys (SMA) under phase and structural transformations and for different statements of the problem, an analytical solution of the problem of stability of an SMA rod undergoing a direct martensitic phase transformation under the action of a compressive load is obtained. It is shown that taking account of the nonlinearity of the deformation process and structural transformation in the transition into the adjacent form of equilibrium significantly changes the solution for sufficiently flexible rods. At the same time, taking into account the strains developed in a phase transition is topical for thick-walled SMA elements.     

Simulation of the motion of thin rotating jets of a viscous incompressible fluid

A regular procedure is proposed for deriving approximate equations of motion of straight thin rotating jets of a viscous incompressible fluid, and similarity parameters of such flows are established. The problem of the stability of free steady motion of a finite jet is considered in the framework of a model that takes into account in the zeroth approximation the effects of viscosity, the rotation of the jet, and capillarity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 51–59, March–April, 1984.     

Batchelor’s theory extended to elongated cylindrical or ellipsoidal particles

Extensions to Batchelor’s theory have been derived to take into account different shaped particles while relating extensional viscosity enhancement to three parameters — shape, volume fraction and particle aspect ratio. The extended theory now allows calculation of the extensional viscosity enhancement, at a given volume fraction of particles, for either ellipsoidal or cylindrical particles. The formula improves the predictive capability of Batchelor’s theory when compared with measurements found in the literature for different rod-like polymer solutions.