Nonlinear convective stability problems of viscoelastic fluids in finite domains

A Chebyshev pseudospectral method is generalized to solve the nonlinear hydrodynamic stability problems of Rayleigh-Bénard convection of viscoelastic fluids in finite domains, which are compatible with the experimental situations, for the range of viscoelastic parameters where the exchange of stabilities is valid. The effects of box aspect ratio, the Deborah number 5 and the dimensionless retardation time ) on the critical Rayleigh number and convection intensity are investigated. The comparison of these results with the experimental data might be used to guide the selection of constitutive equations and to estimate viscoelastic parameter values. The present technique of hydrodynamic stability analysis is quite versatile and can be employed to solve other hydrodynamic stability problems in finite domains.     

Dispersion of linear gravity waves on a viscoelastic fluid in an horizontal canal

A characteristic equation is derived that describes the spatial decay of linear surface gravity waves on Maxwell fluids. Except at small frequencies, the derived dispersion relation is different from the temporal decay dispersion relation which is normally studied within fluid mechanics. The implications for waves on viscous Newtonian fluids using the two different dispersion relations is briefly discussed. The wave number is measured experimentally as function of the frequency in a horizontal canal. Seven Newtonian fluids and four viscoelastic liquids with constant viscosity have been used in the experiments. The spatial decay theory for Newtonian fluids fits well to the experimental data. The model and experiments are used to determine limits for the Maxwell fluid time numbers for the four viscoelastic liquids. As a result of low viscosity it was not possible within this study to obtain these time numbers from oscillatory experiments. Therefore, a comparison of surface gravity wave experiments with theory is applicable as a method to evaluate memory times of low viscosity viscoelastic fluids.     

Shear wave dispersion and attenuation in periodic systems of alternating solid and viscous fluid layers

The attenuation and dispersion of elastic waves in fluid-saturated rocks due to the viscosity of the pore fluid is investigated using an idealized exactly solvable example of a system of alternating solid and viscous fluid layers. Waves in periodic layered systems at low frequencies are studied using an asymptotic analysis of Rytov’s exact dispersion equations. Since the wavelength of shear waves in fluids (viscous skin depth) is much smaller than the wavelength of shear or compressional waves in solids, the presence of viscous fluid layers necessitates the inclusion of higher terms in the long-wavelength asymptotic expansion. This expansion allows for the derivation of explicit analytical expressions for the attenuation and dispersion of shear waves, with the directions of propagation and of particle motion being in the bedding plane. The attenuation (dispersion) is controlled by the parameter which represents the ratio of Biot’s characteristic frequency to the viscoelastic characteristic frequency. If Biot’s characteristic frequency is small compared with the viscoelastic characteristic frequency, the solution is identical to that derived from an anisotropic version of the Frenkel–Biot theory of poroelasticity. In the opposite case when Biot’s characteristic frequency is greater than the viscoelastic characteristic frequency, the attenuation/dispersion is dominated by the classical viscoelastic absorption due to the shear stiffening effect of the viscous fluid layers. The product of these two characteristic frequencies is equal to the squared resonant frequency of the layered system, times a dimensionless proportionality constant of the order 1. This explains why the visco-elastic and poroelastic mechanisms are usually treated separately in the context of macroscopic (effective medium) theories, as these theories imply that frequency is small compared to the resonant (scattering) frequency of individual pores.     

Slow Motion of a Porous Cylindrical Shell in a concentric cylindrical cavity

This paper concerns the Slow Motion of a Porous Cylindrical Shell in a concentric cylindrical cavity using particle-in-cell method. The Brinkman’s equation in the porous region and the Stokes equation for clear fluid in their stream function formulations are used. The hydrodynamic drag force acting on each porous cylindrical particle in a cell and permeability of membrane built up by cylindrical particles with a porous shell are evaluated. Four known boundary conditions on the hypothetical surface are considered and compared: Happel’s, Kuwabara’s, Kvashnin’s and Cunningham’s (Mehta-Morse’s condition). Some previous results for hydrodynamic drag force and dimensionless hydrodynamic permeability have been verified. Variation of the drag coefficient and dimensionless hydrodynamic permeability with permeability parameter σ, particle volume fraction γ has been studied and some new results are reported. The flow patterns through the regions have been analyzed by stream lines. Effect of particle volume fraction γ and permeability parameter σ on flow pattern is also discussed. In our opinion, these results will have significant contributions in studying, Stokes flow through cylindrical swarms.     

Linear gravity waves on Maxwell fluids of finite depth

Linear surface gravity waves on Maxwell viscoelastic fluids with finite depth are studied in this paper. A dispersion equation describing the spatial decay of the gravity wave in finite depth is derived. A dimensionless memory (time) number 0 is introduced. The dispersion equation for the pure viscous fluid will be a specific case of the dispersion equation for the viscoelastic fluid as θ=0. The complex dispersion equation is numerically solved to investigate the dispersion relation. The influences of θ and water depth on the dispersion characteristics and wave decay are discussed. It is found that the role of elasticity for the Maxwell fluid is to make the surface gravity wave on the Maxwell fluid behave more like the surface gravity wave on the inviscid fluid.     

A generalized equation for rheology of emulsions and suspensions of deformable particles subjected to simple shear at low Reynolds number

We present analyses to provide a generalized rheological equation for suspensions and emulsions of non-Brownian particles. These multiparticle systems are subjected to a steady straining flow at low Reynolds number. We first consider the effect of a single deformable fluid particle on the ambient velocity and stress fields to constrain the rheological behavior of dilute mixtures. In the homogenization process, we introduce a first volume correction by considering a finite domain for the incompressible matrix. We then extend the solution for the rheology of concentrated system using an incremental differential method operating in a fixed and finite volume, where we account for the effective volume of particles through a crowding factor. This approach provides a self-consistent method to approximate hydrodynamic interactions between bubbles, droplets, or solid particles in concentrated systems. The resultant non-linear model predicts the relative viscosity over particle volume fractions ranging from dilute to the the random close packing in the limit of small deformation (capillary or Weissenberg numbers) for any viscosity ratio between the dispersed and continuous phases. The predictions from our model are tested against published datasets and other constitutive equations over different ranges of viscosity ratio, volume fraction, and shear rate. These comparisons show that our model, is in excellent agreement with published datasets. Moreover, comparisons with experimental data show that the model performs very well when extrapolated to high capillary numbers (C a?1). We also predict the existence of two dimensionless numbers; a critical viscosity ratio and critical capillary numbers that characterize transitions in the macroscopic rheological behavior of emulsions. Finally, we present a regime diagram in terms of the viscosity ratio and capillary number that constrains conditions where emulsions behave like Newtonian or Non-Newtonian fluids.     

Suspensions of monodisperse spheres in polymer melts: particle size effects in extensional flow

Suspensions in polymeric, viscoelastic liquids have been studied in uniaxial extensional flow. The fibre wind-up technique has been used for this purpose. The effects of particle size and particle volume fraction have been investigated, using monodisperse, spherical particles. The results have been compared with shear flow data on the same materials. The values of the relative extensional viscosities at low stretching rates are in agreement with the relative shear viscosities and relative moduli. This indicates that hydrodynamic forces are stronger than the particle interaction forces. At larger strain rates strain hardening occurs; it is suppressed when particles are added. Small aggregating particles reduce the strain hardening more strongly than larger particles; strain hardening can even be totally eliminated. When further increasing the stretching rate, hydrodynamic effects dominate again and the effect of particle size effect on strain hardening disappears.     

Spectrum of density fluctuations in a particle-fluid system— i. monodisperse spheres

A method for calculating the density autocorrelation ′(x)′(x + r) for a homogeneous particle-fluid system in both physical and Fourier transform space has been developed. The density autocorrelation was related to two quantities, the Overlap function which is defined as the volume of intersection of two spheres as a function of the separation distance and the radial distribution function (RDF) of the particles. In dimensionless co-ordinates, the parameter that characterizes the density autocorrelation is the volume fraction of particles, ₁, , or equivalently the dimensionless mean separation distance (normalized by the particle diameter), . For an isotropic randomly distributed system of particles, the density autocorrelation was observed to oscillate with the correlation distance r, with a wavelength that was proportional to λ. The Fourier transform of the autocorrelation likewise oscillated with the wavenumber k, however the effect of changes in the particle volume fraction was limited to the first peak only. Subsequent peaks were more closely associated with the Overlap function.

The results for the density autocorrelation were extended to a particle-fluid system which experienced an asymptotically large pressure gradient. This initially produced a uniform relative motion between the two fields. In this limit, other higher-order moments such as the Reynolds stress can be related to the density autocorrelation in a straightforward manner. Moreover the spectral shapes of all moments collapse onto the density autocorrelation spectrum in this limit. It was pointed out that the uniform relative motion will eventually become unstable because of hydrodynamic forces on the particles induced by the relative motion. This effect was estimated by introducing a mildly attractive force into the RDF. The results demonstrated that the induced hydrodynamic force promoted a shift in the density spectrum toward small k (large scale) indicating an alternative mechanism for growth in the integral length scale.     

Euler-Euler CFD modeling of fluidized bed＂ Influence of specularity coefficient on hydrodynamic behavior

Euler-Euler two-fluid model is used to simulate the hydrodynamics of gas-solid flow in a bubbling flu- idized bed with Geldert B particles where the solid property is calculated by applying the kinetic theory of granular flow （KTGF）. Johnson and Jackson wall boundary condition is used for the particle phase, and different amount of slip between particle and wall is given by varying the specularity coefficient （φ） from 0 to 1. The simulated particle velocity, granular temperature and particle volume fraction are compared to investigate the effect of different wall boundary conditions on the hydrodynamic behavior, Some of the results are also compared with the available experimental data from the literature. It was found that the model predictions are sensitive to the specularity coefficient. The hydrodynamic behavior deviated sig- nificantly for φ = 0 and φ = 0.01 with maximum deviation found at φ = 0 i.e. free-slip condition. However, the overall bed height predicted by all the conditions is similar.     

钛酸钙系电流变液在直流电场下损耗模量的研究

研究了直流电场下钛酸钙系电流变液在动态振荡剪切模式下损耗模量的变化。用流变仪测试了不同颗粒体积分数材料在不同的温度和电场下损耗模量随应变的关系曲线。讨论了颗粒体积分数、温度、电场强度以及频率对材料损耗模量的影响。以理论推导和试验数据拟合的方法给出了损耗模量与颗粒体积分数、振荡频率、温度以及电场强度等参数的半经验数学关系式。理论值对比试验结果表明,本文损耗模量表达式与试验结果符合较好,可以用于预测直流电场下钛酸钙系电流变液在动态振荡剪切模式下颗粒体积分数、温度、电场强度、频率和剪切应变对损耗模量的影响。     

Laminar flow in three-dimensional square–square expansions

In this work we investigate the three-dimensional laminar flow of Newtonian and viscoelastic fluids through square–square expansions. The experimental results obtained in this simple geometry provide useful data for benchmarking purposes in complex three-dimensional flows. Visualizations of the flow patterns were performed using streak photography, the velocity field of the flow was measured in detail using particle image velocimetry and additionally, pressure drop measurements were carried out. The Newtonian fluid flow was investigated for the expansion ratios of 1:2.4, 1:4 and 1:8 and the experimental results were compared with numerical predictions. For all expansion ratios studied, a corner vortex is observed downstream of the expansion and an increase of the flow inertia leads to an enhancement of the vortex size. Good agreement is found between experimental and numerical results. The flow of the two non-Newtonian fluids was investigated experimentally for expansion ratios of 1:2.4, 1:4, 1:8 and 1:12, and compared with numerical simulations using the Oldroyd-B, FENE-MCR and sPTT constitutive equations. For both the Boger and shear-thinning viscoelastic fluids, a corner vortex appears downstream of the expansion, which decreases in size and strength when the elasticity of the flow is increased. For all fluids and expansion ratios studied, the recirculations that are formed downstream of the square–square expansion exhibit a three-dimensional structure evidenced by a helical flow, which is also predicted in the numerical simulations.     

Rheological characterization of bimodal colloidal dispersions

In order to investigate the effect of the particle size distribution on the rheological properties of concentrated colloidal dispersions both steady-state shear and oscillatory measurements have been performed on well-characterized bimodal dispersions of sterically stabilized PMMA particles. Replacing a minor amount of large particles by small ones in a concentrated dispersion, keeping the total effective volume fraction constant, decreases the viscosity quite drastically. On the other hand, replacing a small amount of small particles by big ones hardly effects the viscosity at all. This behavior can be attributed to the deformability of the stabilizing polymer layer. A procedure is proposed to calculate the limiting viscosities in a bimodal colloidal dispersion starting from the characteristics of the monodisperse systems. A good agreement has been obtained between the calculated values and the experimental results. The linear viscoelastic properties of the concentrated dispersions have been investigated by means of oscillatory measurements. The plateau values of the storage modulus for the bimodal dispersions decrease with an increasing fraction of the coarse particles. By substituting the bimodal dispersion by an equivalent monodisperse system the storage modulus can be superimposed on the values for the monodisperse suspensions when plotted as a function of the mean interparticle distance.     

Shear-augmented solute dispersion during drug delivery for three-layer flow through microvessel under stress jump and momentum slip-Darcy model

Nanoparticle(drug particle) dispersion is an important phenomenon during nanodrug delivery in the bloodstream by using multifunctional carrier particles. The aim of this study is to understand the dispersion of drug particle(nanoparticle) transport during steady blood flow through a microvessel. A two-phase fluid model is considered to define blood flow through a microvessel. Plug and intermediate regions are defined by a non-Newtonian Herschel-Bulkley fluid model where the plug region appears due to the aggregation of red blood cells at the axis in the vessel. The peripheral(porous in nature)region is defined by the Newtonian fluids. The wall of the microvessel is considered to be permeable and characterized by the Darcy model. Stress-jump and velocity slip conditions are incorporated respectively at the interface of the intermediate and peripheral regions and at the inner surface of the microvessel. The effects of the rheological parameter, the pressure constant, the particle volume fraction, the stress jump constant, the slip constant,and the yield stress on the dispersion are analyzed and discussed. It is observed that the non-dimensional pressure gradient and the yield stress enhance the dispersion rate of the nanoparticle, while the opposite trends are observed for the velocity slip constant, the nanoparticle volume fraction, the rheological parameter, and the stress-jump constant.     

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.     

A dynamic investigation of fiber-reinforced viscoelastic materials

In recent years, several mathematical models have been proposed to describe the quasi-static response of fiber-reinforced materials, consisting of continuous, elastic fibers embedded in a linear viscoelastic matrix. By assuming that geometric dispersion (dispersion resulting from the internal geometry of the material) is small in comparison to viscoelastic dispersion (dispersion resulting from the viscoelastic nature of the material), these proposed constitutive equations can be extended from a quasi-static regime to a dynamic regime. Here, we examine how the extension to the dynamic regime may be accomplished, compare the results with a theoretical model that includes geometric dispersion, and use the results of an experimental program to evaluate the models. In general, the quasi-static constitutive equations predict phase velocities that are larger than that predicted by the model which contains geometric dispersion and attenuation coefficients that are lower; and, the experimental results agree with the theoretical predictions, provided the fibers were spread more or less uniformly over the cross section.     

Wave propagation analysis of porous functionally graded piezoelectric nanoplates with a visco-Pasternak foundation

This study investigates the size-dependent wave propagation behaviors under the thermoelectric loads of porous functionally graded piezoelectric(FGP) nanoplates deposited in a viscoelastic foundation. It is assumed that(i) the material parameters of the nanoplates obey a power-law variation in thickness and(ii) the uniform porosity exists in the nanoplates. The combined effects of viscoelasticity and shear deformation are considered by using the Kelvin-Voigt viscoelastic model and the refined hi...     

Numerical simulations of dispersive mixing of viscoelastic suspensions in a four-roll mill

The dispersive mixing of particles suspended in Newtonian and viscoelastic fluids in a four-roll mill is studied by direct numerical simulations. A fictitious domain method is used to handle the particle motion. To quantify the mixing, a proper mixing distribution function is defined. The combined effect of fluid rheology and particle-particle/particle-wall hydrodynamic interactions is addressed. At variance with the Newtonian case where the particle distribution remains uniform, the viscoelasticity-induced migration leads to a significant segregation process. The effect of the Deborah number (the product of the fluid relaxation time and the roll angular velocity), shear-thinning, particle concentration, and size on the microstructure evolution is thoroughly investigated.     

On viscoelastic effects in non-Newtonian steady flows through porous media

This paper presents a mathematical model for describing approximately the viscoelastic effects in non-Newtonian steady flows through a porous medium. The rheological behaviour of power law fluids is considered in the Maxwell model of elastic behaviour of the fluids. The equations governing the steady flow through porous media are derived and an analytical solution of these equations in the case of a simple flow system is obtained. The conditions for which the viscoelastic effects may become observable from the pressure distribution measurements are shown and expressed in terms of some dimensionless groups. These have been found to be relevant in the evaluation of viscoelastic effects in the steady flow through porous media.     

The dynamics and dissolution of gas bubbles in a viscoelastic fluid

This paper reports an experimental study of the motion of dissolving and non-dissolving gas bubbles in a quiescent viscoelastic fluid. The objective of the investigation was to determine the influence of the abrupt transition in bubble velocity, which had been observed at a critical radius of approx. on the rate of mass transfer. Thus, a range of bubble sizes from an equivalent (spherical) radius of 0.2–0.4 cm was employed using CO₂ gas, and five different fluids, including one Newtonion glycerine/water solution and four viscoelastic solutions of Separan AP30 in water (0.1, 0.5, 1% by weight) and in a water/glycerine mixture.The experimental data on bubble velocity shows that the discontinuous increase with bubble volume observed previously for air bubbles in viscoelastic fluids, does not occur for dissolving CO₂ bubbles—presumably due to the continuous decrease in bubble volume. Instead, a very steep but definitely continuous transition is found. Mass transfer rates are found to be significantly enhanced by viscoelasticity, and comparison with available theoretical results shows that the increase is greater than expected for purely viscous, power-law fluids. We conclude that a fully viscoelastic constitutive model would be necessary for a successful analysis of the dissolution of a gas bubble which is translating through a (high molecular weight) polymer solution.     

Hydrodynamic force and torque models for a particle moving near a wall at finite particle Reynolds numbers

This research work is aimed at proposing models for the hydrodynamic force and torque experienced by a spherical particle moving near a solid wall in a viscous fluid at finite particle Reynolds numbers. Conventional lubrication theory was developed based on the theory of Stokes flow around the particle at vanishing particle Reynolds number. In order to account for the effects of finite particle Reynolds number on the models for hydrodynamic force and torque near a wall, we use four types of simple motions at different particle Reynolds numbers. Using the lattice Boltzmann method and considering the moving boundary conditions, we fully resolve the flow field near the particle and obtain the models for hydrodynamic force and torque as functions of particle Reynolds number and the dimensionless gap between the particle and the wall. The resolution is up to 50 grids per particle diameter. After comparing numerical results of the coefficients with conventional results based on Stokes flow, we propose new models for hydrodynamic force and torque at different particle Reynolds numbers. It is shown that the particle Reynolds number has a significant impact on the models for hydrodynamic force and torque. Furthermore, the models are validated against general motions of a particle and available modeling results from literature. The proposed models could be used as sub-grid scale models where the flows between particle and wall can not be fully resolved, or be used in Lagrangian simulations of particle-laden flows when particles are close to a wall instead of the currently used models for an isolated particle.