The constitutive equation of a dilute suspension of spherical microcapsules

This paper studies the rheological behavior of a dilute suspension of spherical microcapsules, i.e. spherical thin elastic membranes filled with an incompressible liquid. Previous results obtained for the motion of such capsules freely suspended in a simple shear flow are first extended to general linear shear flows. Then the stresslet term in the external flow is computed to order ?², where ?, assumed to be small with respect to unity, is the ratio of viscous to elastic forces acting on the particle. The resulting constitutive equation is of the viscoelastic type and is similar to the one obtained for liquid droplets. It predicts that the microcapsule suspension exhibits a shear dependent viscosity and normal stress effects. The exact dependency of these phenomena on the microscopic parameters of the suspension is explicitly provided by the model.     

Modeling and simulation of three‐dimensional extrusion swelling of viscoelastic fluids with PTT,Giesekus and FENE‐P constitutive models

The investigation of the extrusion swelling mechanism of viscoelastic fluids has both scientific and industrial interest. However, it has been traditionally difficult to afford theoretical and experimental researches to this problem. The numerical methodology based on the penalty finite element method with a decoupled algorithm is presented in the study to simulate three‐dimensional extrusion swelling of viscoelastic fluids flowing through out of a circular die. The rheological responses of viscoelastic fluids are described by using three kinds of differential constitutive models including the Phan‐Thien Tanner model, the Giesekus model, and the finite extensible nonlinear elastic dumbbell with a Peterlin closure approximation model. A streamface‐streamline method is introduced to adjust the swelling free surface. The calculation stability is improved by using the discrete elastic‐viscous split stress algorithm with the inconsistent streamline‐upwind scheme. The essential flow characteristics of viscoelastic fluids are predicted by using the proposed numerical method, and the mechanism of swelling phenomenon is further discussed.Copyright © 2013 John Wiley & Sons, Ltd.     

Effective equations describing the flow of a viscous incompressible fluid through a long elastic tubeEquations efficaces décrivant l'écoulement d'un fluide visqueux incompressible dans un tuyau long et élastique

We study the flow of a viscous incompressible fluid through a long and narrow elastic tube whose walls are modeled by the Navier equations for a curved, linearly elastic membrane. The flow is governed by a given small time dependent pressure drop between the inlet and the outlet boundary, giving rise to creeping flow modeled by the Stokes equations. By employing asymptotic analysis in thin, elastic, domains we obtain the reduced equations which correspond to a Biot type viscoelastic equation for the effective pressure and the effective displacement. The approximation is rigorously justified by obtaining the error estimates for the velocity, pressure and displacement. Applications of the model problem include blood flow in small arteries. We recover the well-known Law of Laplace and provide a new, improved model when shear modulus of the vessel wall is not negligible. To cite this article: S. ?ani?, A. Mikeli?, C. R. Mecanique 330 (2002) 661–666.     

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.     

A numerical study of constitutive models endowed with Pom-Pom molecular attributes

The branched polymer melts are modeled respectively in this investigation by the existing XPP and PTT–XPP models, along with the proposed S-MDCPP (Single/Simplified Modified Double Convected Pom-Pom) model developed on the basis of the existing MDCPP model. A pressure stabilized mass equation is formulated with the finite increment calculus (FIC) process to restrain and further eliminate spurious oscillations of pressure field due to the incompressibility of fluids. The discrete elastic viscous stress splitting (DEVSS) technique is employed, in order to retain an elliptic contribution in the weak form of the momentum equation. An inconsistent streamline-upwind (SU) method is applied to spatially discretize the constitutive equations. The mass, momentum conservation and constitutive equations are discretized and solved by the iterative stabilized fractional step algorithm along with the Crank–Nicolson implicit difference scheme. Thus the finite elements with equal low-order interpolation approximations for velocity–pressure–stress variables can be devised to numerically simulate the viscoelastic contraction flows for branched LDPE melts. The influences of the three viscoelastic constitutive models and the branched arms at the end of the Pom-Pom molecule on the rheological behaviors occurring in this complex flow are discussed. The numerical results demonstrate that the proposed S-MDCPP model is capable of reproducing some properties similar to those predicted by the XPP model in high shear flow and, on the other hand, reproducing some properties similar to those predicted by the PTT–XPP model in high elongational flow. Furthermore, the proposed S-MDCPP model is capable of well identifying the macromolecule topological structures of branched polymer melts.     

改进的非定常Bingham岩石蠕变模型及参数辨识

传统Bingham模型适用于描述蠕变的前两个阶段,为了更加准确地描述岩体蠕变的三个阶段,本文提出了一种改进的非定常Bingham蠕变模型,并推导出其蠕变本构方程．该模型在传统Bingham模型的基础上串联一个非定常黏壶元件,并考虑岩体弹性模量及原有黏壶元件蠕变参数的非定常化．基于该模型对泥岩的蠕变试验曲线进行拟合及参数辨识,结果表明该模型具有较高的拟合精度,可以准确地描述泥岩蠕变的三个阶段,并实现了在数值软件中的应用,为实际工程提供了一种简单有效的蠕变本构模型．     

Nonlinear Fractional Order Viscoelasticity at Large Strains

In this paper, we formulate a fractional order viscoelastic model for large deformations and develop an algorithm for the integration of the constitutive response. The model is based on the multiplicative split of the deformation gradient into elastic and viscous parts. Further, the stress response is considered to be composed of a nonequilibrium part and an equilibrium part. The viscous part of the deformation gradient (here regarded as an internal variable) is governed by a nonlinear rate equation of fractional order. To solve the rate equation the finite element method in time is used in combination with Newton iterations. The method can handle nonuniform time meshes and uses sparse quadrature for the calculations of the fractional order integral. Moreover, the proposed model is compared to another large deformation viscoelastic model with a linear rate equation of fractional order. This is done by computing constitutive responses as well as structural dynamic responses of fictitious rubber materials.     

Nonlinear Fractional Order Viscoelasticity at Large Strains

In this paper, we formulate a fractional order viscoelastic model for large deformations and develop an algorithm for the integration of the constitutive response. The model is based on the multiplicative split of the deformation gradient into elastic and viscous parts. Further, the stress response is considered to be composed of a nonequilibrium part and an equilibrium part. The viscous part of the deformation gradient (here regarded as an internal variable) is governed by a nonlinear rate equation of fractional order. To solve the rate equation the finite element method in time is used in combination with Newton iterations. The method can handle nonuniform time meshes and uses sparse quadrature for the calculations of the fractional order integral. Moreover, the proposed model is compared to another large deformation viscoelastic model with a linear rate equation of fractional order. This is done by computing constitutive responses as well as structural dynamic responses of fictitious rubber materials.     

Steady viscoelastic film flow over 2D topography: I. The effect of viscoelastic properties under creeping flow

Two-dimensional, steady flow of a viscoelastic film over a periodic topography under the action of a body force is studied. The exponential Phan-Thien and Tanner (ePTT) constitutive model is used. The conservation equations are solved via the usual mixed finite element method combined with a quasi-elliptic grid generation scheme in order to capture the large deformations of the free surface. The constitutive equation is weighted using the SUPG method and solved via the polymeric stress splitting EVSS-G technique. First, the code is validated by verifying that in isolated topographies the periodicity conditions result in fully developed viscoelastic film flow at the inflow/outflow boundaries and that its predictions for Newtonian fluids over 2D topography under creeping flow conditions coincide with those of previous works. Since the lubrication approximation is not invoked here, the topographical features can have wall segments that form any angle with the main flow, but only slight smoothing of the convex corners assists in reducing the stress singularity there. Thus, steady-state solutions are computed accurately up to high Deborah numbers, resulting in large deformations of the free surface. The magnitude of the capillary ridge in the film before the entrance to a step down of the substrate and of the capillary depression before a step up is increased as De increases up to ～0.7 due to increased fluid elasticity. Above this value they decrease, because increasing De increases also the shear and elongational thinning, which eventually affect them more. Increasing the ratio of solvent to polymer viscosities, β, the elongational parameter, ? and the molecular slip parameter, ξ, monotonically increases their magnitudes and especially that of the capillary ridge, but the mechanisms leading to these changes are different as explained in the text.     

One Equation Model for Turbulent Channel Flow with Second Order Viscoelastic Corrections

A modified second order viscoelastic constitutive equation is used to derive a k–l type turbulence closure to qualitatively assess the effects of elastic stresses on fully-developed channel flow. Specifically, the second order correction to the Newtonian constitutive equation gives rise to a new term in the momentum equation involving the time-averaged elastic shear stress and in the turbulent kinetic energy transport equation quantifying the interaction between the fluctuating elastic stress and rate of strain tensors, denoted by P _w , for which a closure is developed and tested. This closure is based on arguments of isotropic turbulence and equilibrium in boundary layer flows and a priori P _w could be either positive or negative. When P _w is positive, it acts to reduce the production of turbulent kinetic energy and the turbulence model predictions qualitatively agree with direct numerical simulation (DNS) results obtained for more realistic viscoelastic fluid models with memory which exhibit drag reduction. In contrast, P _w  < 0 leads to a drag increase and numerical breakdown of the model occurs at very low values of the Deborah number, which signifies the ratio of elastic to viscous stresses. Limitations of the turbulence model primarily stem from the inadequacy of the k–l formulation rather than from the closure for P _w . An alternative closure for P _w , mimicking the viscoelastic stress work predicted by DNS using the Finitely Extensible Nonlinear Elastic-Peterlin fluid model, which is mostly characterized by P _w  > 0 but has also a small region of negative P _w in the buffer layer, was also successfully tested. This second model for P _w leads to predictions of drag reduction, in spite of the enhancement of turbulence production very close to the wall, but the equilibrium conditions in the inertial sub-layer were not strictly maintained.     

An investigation of non-newtonian flow past a sphere

Any experimental work on the flow of a polymer solution or any theoretical analysis on the basis of a visoelastic constitutive equation does not always bring out viscoelastic effects but may be showing a non-Newtonian viscosity effect. Therefore, in order to obtain a clear understanding about viscoelastic effects, it is desirable to have a sufficient knowledge of the non-Newtonian viscosity effect. To facilitate this, finite-difference numerical solutions of non-Newtonian flow were carried out using a non-Newtonian viscous model for the Reynolds numbers of 0.1, 1.0, 20 and 60.Drag force measurements and flow visualization experiments were also performed over a wide range of experimental conditions using polymer solutions. The present work appears to support the following idea: When compared with the Newtonian case on the basis of DV_∞P/η₀, where η₀ is the zero shear viscosity, it is on account of the non-Newtonian viscosity that the friction and pressure drags decrease, that the separating vortices behind the sphere become larger, and that no shift occurs in the streamlines. On the other hand, it is due to viscoelasticity that the normal force drag increases, that the separating vortices behind the sphere become smaller, and that an upstream shift occurs in the streamlines.     

Three‐dimensional simulation of planar contraction viscoelastic flow by penalty finite element method

The planar contraction flow is a benchmark problem for the numerical investigation of viscoelastic flow. The mathematical model of three‐dimensional viscoelastic fluids flow is established and the numerical simulation of its planar contraction flow is conducted by using the penalty finite element method with a differential Phan‐Thien–Tanner constitutive model. The discrete elastic viscous split stress formulation in cooperating with the inconsistent streamline upwind scheme is employed to improve the computation stability. The distributions of velocity and stress obtained by simulation are compared with that of Quinzani's experimental results detected by laser–doppler velocimetry and flow‐induced birefringence technologies. It shows that the numerical results agree well with the experimental results. The numerical methods proposed in the study can be well used to predict complex flow patterns of viscoelastic fluids. Copyright © 2009 John Wiley & Sons, Ltd.     

Understanding thixotropic and antithixotropic behavior of viscoelastic micellar solutions and liquid crystalline dispersions. I. The model

A simple model consisting of the Upper Convected Maxwell constitutive equation and a kinetic equation for destruction and construction of structure, first proposed by Fredrickson in 1970, is used here to reproduce the complex rheological behavior of viscoelastic systems that also exhibit thixotropy and rheopexy under shear flow. The model requires five parameters that have physical significance and that can be estimated from rheological measurements. Several steady and unsteady flow situations were analyzed with the model. The model predicts creep behavior, stress relaxation and the presence of thixotropic loops when the sample is subjected to transient stress cycles. Such behavior has been observed with surfactant-based solutions and dispersions. The role of the characteristic time for structure built up, λ, in the extent and shape of the thixotropic loops is demonstrated.     

Steady laminar flow of viscoelastic fluid in a curved pipe of circular cross-section with varying curvature

The steady laminar flow of viscoelastic fluid through pipes of circular cross-section, whose center-line curvature varies locally, is analyzed theoretically. The flows in three kinds of pipes whose center-lines are specified by

as the examples of once, twice, and periodically curved pipes, respectively, are discussed in comparison with purely viscous flow. The analysis is valid for any other two-dimensionally curved pipes, when the center-line curvature is small. In addition, the reason why the secondary flow of a viscoelastic fluid in a curved pipe of circular cross-section is stronger than that of a purely viscous fluid is explained. In the present paper, the White—Metzner model is employed as the constitutive equation.     

A rheological model for materials which support coexistent shear rates

In this work we study a version of the three constant differential-type Oldroyd constitutive relation which allows distinct objective time derivatives for the extra stress and the stretching. We integrate the constitutive equation and determine an equivalent history integral representation for this model for the general class of viscometric motions. For certain choices of the material parameters and initial conditions, we find that this model allows for the development of shear rate discontinuities in the flow domain as a steady viscometric flow is achieved. Correspondingly, we also give evidence that intense shear rate oscillations may occur during the transient period as an impulsively started viscometric flow in a channel tends to a steady state under a constant critical shear stress. This critical shear stress lies in an interval of values for which the material experiences the phenomenon of “flow yielding”. A qualitative comparison with experimental data is made for certain creams and greases. The material instabilities inherent in this constitutive theory for viscometric motions are suggestive of the instabilities that occur in many viscoelastic fluids such as sharkskin patterns, wavy fracture, and spurt flow.     

Rheological characterization of continuous fiber—reinforced viscous fluid

The present paper describes a micromechanical technique to determine rheological properties of viscous fluid reinforced with unidirectional continuous fibers. Fluid viscosity is described by a shear thinning model and high viscosity is considered for continuous fibers having considerable rigidity compared to net fluid. The microstructure is identified by a representative volume element that is subjected to equivalent macroscopic deformation fields. The energy balance and periodicity conditions are considered to relate deformation and stress in macro and micro-levels. It is shown that response of viscous fluid reinforced with rigid fibers depends on deformation history as well as rate-of-deformation in the transverse intraply shear and transverse squeeze flows. An orthotropic viscous constitutive equation is derived to describe response of such materials. The material viscosities are evaluated for viscous fluid reinforced with different fiber volume fractions during deformation applied in different rates of deformation. The results are used to derive the functions predicting effective anisotropic viscosities of reinforced fluid.     

VISCOELASTIC CONSTITUTIVE MODEL RELATED TO DEFORMATION OF INSECT WING UNDER LOADING IN FLAPPING MOTION

Flexible insect wings deform passively under the periodic loading during napping flight. The wing flexibility is considered as one of the specific mechanisms on improving insect flight performance. The constitutive relation of the insect wing material plays a key role on the wing deformation, but has not been clearly understood yet. A viscoelastic constitutive relation model was established based on the stress relaxation experiment of a dragonfly wing (in vitro). This model was examined by the finite element analysis of the dynamic deformation response for a model insect wing under the action of the periodical inertial force in flapping. It is revealed that the viscoelastic constitutive relation is rational to characterize the biomaterial property of insect wings in contrast to the elastic one. The amplitude and form of the passive viscoelastic deformation of the wing is evidently dependent on the viscous parameters in the constitutive relation.     

An Eulerian formulation for large deformations of anisotropic elastic and viscoelastic solids and viscous fluids

An Eulerian formulation has been developed for the constitutive response of a group of materials that includes anisotropic elastic and viscoelastic solids and viscous fluids. The material is considered to be a composite of an elastic solid and a viscous fluid. Evolution equations are proposed for a triad of vectors m _i that represent the stretches and orientations of material line elements in the solid component. Evolution equations for an orthonormal triad of vectors s _i are also proposed to characterize anisotropy of the fluid component. In particular, for an elastic solid it is shown that the material response is totally characterized by the functional form of the strain energy and by the current values of m _i , which are measurable in the current state of the material. Moreover, it is shown that the proposed Eulerian formulation removes unphysical arbitrariness of the choice of the reference configuration in the standard formulation of constitutive equations for anisotropic elastic solids.     

An adaptive viscoelastic stress splitting scheme and its applications: AVSS/SI and AVSS/SUPG

We report an adaptive viscoelastic stress splitting (AVSS) scheme, which was found to be extremely robust in the numerical simulation of viscoelastic flow involving steep stress boundary layers. The scheme is different from the elastic viscous split stress (EVSS) in that the local Newtonian component is allowed to depend adaptively on the magnitude of the local elastic stress. Two decoupled versions of the scheme were implemented for the Upper Convected Maxwell (UCM) model: the streamline integration (AVSS/SI), and the streamline upwind Petrov-Galerkin (AVSS/SUPG) methods of integrating the stress. The implementations were benchmarked against the known analytic Poiseuille solution, and no upper limiting Weissenberg number was found (the computation was stopped at Weissenberg number of O(10⁴)). The flow past a sphere in a tube was solved next, giving convergent results up to a Weissenberg number of 3.2 with the AVSS/SI and 1.55 with the AVSS/SUPG (both were decoupled schemes; without the adaptive scheme, the limiting Weissenberg number for the decoupled streamline integration method was about 0.3). These results show that the decoupled AVSS is more stable than the corresponding EVSS, and the SI is more robust than SUPG in solving the constitutive equation of hyperbolic type. In addition, we found a very long stress wake behind the sphere, and a weak vortex in the rear stagnation region at a Weissenberg number above W_i of about 1.6, which does not seem to increase in size or strength with increasing W_i.     

A FENE-P k–ε turbulence model for low and intermediate regimes of polymer-induced drag reduction

A new low-Reynolds-number k–ε turbulence model is developed for flows of viscoelastic fluids described by the finitely extensible nonlinear elastic rheological constitutive equation with Peterlin approximation (FENE-P model). The model is validated against direct numerical simulations in the low and intermediate drag reduction (DR) regimes (DR up to 50%). The results obtained represent an improvement over the low DR model of Pinho et al. (2008) [A low Reynolds number k–ε turbulence model for FENE-P viscoelastic fluids, Journal of Non-Newtonian Fluid Mechanics, 154, 89–108]. In extending the range of application to higher values of drag reduction, three main improvements were incorporated: a modified eddy viscosity closure, the inclusion of direct viscoelastic contributions into the transport equations for turbulent kinetic energy (k) and its dissipation rate, and a new closure for the cross-correlations between the fluctuating components of the polymer conformation and rate of strain tensors (NLT_ij). The NLT_ij appears in the Reynolds-averaged evolution equation for the conformation tensor (RACE), which is required to calculate the average polymer stress, and in the viscoelastic stress work in the transport equation of k. It is shown that the predictions of mean velocity, turbulent kinetic energy, its rate of dissipation by the Newtonian solvent, conformation tensor and polymer and Reynolds shear stresses are improved compared to those obtained from the earlier model.