A Reynolds stress model for turbulent flows of viscoelastic fluids

A second-order closure is developed for predicting turbulent flows of viscoelastic fluids described by a modified generalised Newtonian fluid model incorporating a nonlinear viscosity that depends on a strain-hardening Trouton ratio as a means to handle some of the effects of viscoelasticity upon turbulent flows. Its performance is assessed by comparing its predictions for fully developed turbulent pipe flow with experimental data for four different dilute polymeric solutions and also with two sets of direct numerical simulation data for fluids theoretically described by the finitely extensible nonlinear elastic - Peterlin model. The model is based on a Newtonian Reynolds stress closure to predict Newtonian fluid flows, which incorporates low Reynolds number damping functions to properly deal with wall effects and to provide the capability to handle fluid viscoelasticity more effectively. This new turbulence model was able to capture well the drag reduction of various viscoelastic fluids over a wide range of Reynolds numbers and performed better than previously developed models for the same type of constitutive equation, even if the streamwise and wall-normal turbulence intensities were underpredicted.

Specific shear-dependent viscoelastic third-grade fluid model

A specific modified constitutive equation for a third-grade fluid is proposed so that the model be suitable for applications where shear-thinning or shear-thickening may occur. For that, we use the Cosserat theory approach reducing the exact three-dimensional equations to a system depending only on time and on a single spatial variable. This one-dimensional system is obtained by integrating the linear momentum equation over the cross-section of the tube, taking a velocity field approximation provided by the Cosserat theory. From this reduced system, we obtain the unsteady equations for the wall shear stress and mean pressure gradient depending on the volume flow rate, Womersley number, viscoelastic coefficient and flow index over a finite section of the tube geometry with constant circular cross-section.

A stress relaxation model for the viscoelastic solids based on the steady-state creep equation

Creep and stress relaxation are inherent mechanical behaviors of viscoelastic materials. It is considered that both are different performances of one identical physical phenomenon. The relationship between the decay stress and time during stress relaxation has been derived from the power law equation of the steady-state creep. The model was used to analyse the stress relaxation curves of various different viscoelastic materials (such as pure polycrystalline ice, polymers, foods, bones, metal, animal tissues, etc.). The calculated results using the theoretical model agree with the experimental data very well. Here we show that the new mathematical formula is not only simple but its parameters have the clear physical meanings. It is suitable to materials with a very broad scope and has a strong predictive ability.

The variant rpoS allele of S-enteritidis strain 27655R does not affect virulence in a chick model nor constitutive curliation but does generate a cold-sensitive phenotype

Adherence of pathogenic Escherichia coli and Salmonella spp. to host cells is in part mediated by curli fimbriae which, along with other virulence determinants, are positively regulated by RpoS. Interested in the role and regulation of curli (SEF17) fimbriae of Salmonella enteritidis in poultry infection, we tested the virulence of naturally occurring S. enteritidis PT4 strains 27655R and 27655S which displayed constitutive and null expression of curli (SEF17) fimbriae, respectively, in a chick invasion assay and analysed their rpoS alleles. Both strains were shown to be equally invasive and as invasive as a wild-type phage type 4 strain and an isogenic derivative defective for the elaboration of curli. We showed that the rpoS allele of 27655S was intact even though this strain was non-curliated and we confirmed that a S. enteritidis rpoS::str(r) null mutant was unable to express curli, as anticipated. Strain 27655R, constitutively curliated, possessed a frameshift mutation at position 697 of the rpoS coding sequence which resulted in a truncated product and remained curliated even when transduced to rpoS::str(r). Additionally, rpoS mutants are known to be cold-sensitive, a phenotype confirmed for strain 27655R. Collectively, these data indicated that curliation was not a significant factor for pathogenesis of S. enteritidis in this model and that curliation of strains 27655R and 27655S was independent of RpoS. Significantly, strain 27655R possessed a defective rpoS allele and remained virulent. Here was evidence that supported the concept that different naturally occurring rpoS alleles may generate varying virulence phenotypic traits. (C) 1998 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Numerical solution of the PTT constitutive equation for unsteady three-dimensional free surface flows

This work deals with the development of a numerical technique for simulating three-dimensional viscoelastic free surface flows using the PTT (Phan-Thien-Tanner) nonlinear constitutive equation. In particular, we are interested in flows possessing moving free surfaces. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are considered. The PTT equation is solved by a high order method, which requires the calculation of the extra-stress tensor on the mesh contours. To validate the numerical technique developed in this work flow predictions for fully developed pipe flow are compared with an analytic solution from the literature. Then, results of complex free surface flows using the FIT equation such as the transient extrudate swell problem and a jet flowing onto a rigid plate are presented. An investigation of the effects of the parameters epsilon and xi on the extrudate swell and jet buckling problems is reported. (C) 2010 Elsevier B.V. All rights reserved.

Numerical Solution of the Upper-Convected Maxwell Model for Three-Dimensional Free Surface Flows

This work presents a finite difference technique for simulating three-dimensional free surface flows governed by the Upper-Convected Maxwell (UCM) constitutive equation. A Marker-and-Cell approach is employed to represent the fluid free surface and formulations for calculating the non-Newtonian stress tensor on solid boundaries are developed. The complete free surface stress conditions are employed. The momentum equation is solved by an implicit technique while the UCM constitutive equation is integrated by the explicit Euler method. The resulting equations are solved by the finite difference method on a 3D-staggered grid. By using an exact solution for fully developed flow inside a pipe, validation and convergence results are provided. Numerical results include the simulation of the transient extrudate swell and the comparison between jet buckling of UCM and Newtonian fluids.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Application of the PTT Model to Axisymmetric Free Surface Flows

This work is concerned with numerical simulation of axisymmetric viscoelastic free surface flows using the Phan-Thien-Tanner (PTT) constitutive equation. A finite difference technique for solving the governing equations for unsteady incompressible flows written in Cylindrical coordinates on a staggered grid is described. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are applied. The numerical method is verified by comparing numerical predictions of fully developed flow in a pipe with the corresponding analytic solutions. To demonstrate that the numerical method can simulate axisymmetric free surface flows governed by the PTT model, numerical results of the flow evolution of a drop impacting on a rigid dry plate are presented. In these simulations, the rheological effects of the parameters epsilon and xi are investigated.

NUMERICAL AND ANALYTICAL VERIFICATION OF A MULTISCALE COMPUTATIONAL MODEL FOR IMPACT PROBLEMS IN HETEROGENEOUS VISCOELASTIC MATERIALS WITH EVOLVING DAMAGE

Composites are engineered materials that take advantage of the particular properties of each of its two or more constituents. They are designed to be stronger, lighter and to last longer which can lead to the creation of safer protection gear, more fuel efficient transportation methods and more affordable materials, among other examples. This thesis proposes a numerical and analytical verification of an in-house developed multiscale model for predicting the mechanical behavior of composite materials with various configurations subjected to impact loading. This verification is done by comparing the results obtained with analytical and numerical solutions with the results found when using the model. The model takes into account the heterogeneity of the materials that can only be noticed at smaller length scales, based on the fundamental structural properties of each of the composite’s constituents. This model can potentially reduce or eliminate the need of costly and time consuming experiments that are necessary for material characterization since it relies strictly upon the fundamental structural properties of each of the composite’s constituents. The results from simulations using the multiscale model were compared against results from direct simulations using over-killed meshes, which considered all heterogeneities explicitly in the global scale, indicating that the model is an accurate and fast tool to model composites under impact loads. Advisor: David H. Allen

Numerical simulation of drop impact and jet buckling problems using the eXtended Pom-Pom model

This work presents numerical simulations of two fluid flow problems involving moving free surfaces: the impacting drop and fluid jet buckling. The viscoelastic model used in these simulations is the eXtended Pom-Pom (XPP) model. To validate the code, numerical predictions of the drop impact problem for Newtonian and Oldroyd-B fluids are presented and compared with other methods. In particular, a benchmark on numerical simulations for a XPP drop impacting on a rigid plate is performed for a wide range of the relevant parameters. Finally, to provide an additional application of free surface flows of XPP fluids, the viscous jet buckling problem is simulated and discussed. (C) 2011 Elsevier B.V. All rights reserved.

Numerical simulation of some viscoelastic fluids

Numerical simulation of the Oldroyd-B type viscoelastic fluids is a very challenging problem. rnThe well-known High Weissenberg Number Problem" has haunted the mathematicians, computer scientists, and rnengineers for more than 40 years. rnWhen the Weissenberg number, which represents the ratio of elasticity to viscosity, rnexceeds some limits, simulations done by standard methods break down exponentially fast in time. rnHowever, some approaches, such as the logarithm transformation technique can significantly improve rnthe limits of the Weissenberg number until which the simulations stay stable. rnrnWe should point out that the global existence of weak solutions for the Oldroyd-B model is still open. rnLet us note that in the evolution equation of the elastic stress tensor the terms describing diffusive rneffects are typically neglected in the modelling due to their smallness. However, when keeping rnthese diffusive terms in the constitutive law the global existence of weak solutions in two-space dimension rncan been shown. rnrnThis main part of the thesis is devoted to the stability study of the Oldroyd-B viscoelastic model. rnFirstly, we show that the free energy of the diffusive Oldroyd-B model as well as its rnlogarithm transformation are dissipative in time. rnFurther, we have developed free energy dissipative schemes based on the characteristic finite element and finite difference framework. rnIn addition, the global linear stability analysis of the diffusive Oldroyd-B model has also be discussed. rnThe next part of the thesis deals with the error estimates of the combined finite element rnand finite volume discretization of a special Oldroyd-B model which covers the limiting rncase of Weissenberg number going to infinity. Theoretical results are confirmed by a series of numerical rnexperiments, which are presented in the thesis, too.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.

Dynamics of gas bubbles encapsulated by a viscoelastic fluid shell under acoustic fields

The dynamics of a gas-filled microbubble encapsulated by a viscoelastic fluid shell immersed in a Newtonian liquid and subject to an external pressure field is theoretically studied. The problem is formulated by considering a nonlinear Oldroyd type constitutive equation to model the rheological behavior of the fluid shell. Heat and mass transfer across the surface bubble have been neglected but radiation losses due to the compressibility of the surrounding liquid have been taken into account. Bubble collapse under sudden increase of the external pressure as well as nonlinear radial oscillations under ultrasound fields are investigated. The numerical results obtained show that the elasticity of the fluid coating intensifies oscillatory collapse and produces a strong increase of the amplitudes of radial oscillations which may become chaotic even for moderate driving pressure amplitudes. The role played by the elongational viscosity has also been analyzed and its influence on both, bubble collapse and radial oscillations, has been recognized. According to the theoretical predictions provided in the present work, a microbubble coated by a viscoelastic fluid shell is an oscillating system that, under acoustic driving, may experience volume oscillations of large amplitude, being, however, more stable than a free bubble. Thus, it could be expected that such a system may have a suitable behavior as an echogenic agent.

Nonuniformity of the constitutive law parameters for shear rupture and quasistatic nucleation to dynamic rupture: a physical model of earthquake generation processes.

Based on the recent high-resolution laboratory experiments on propagating shear rupture, the constitutive law that governs shear rupture processes is discussed in view of the physical principles and constraints, and a specific constitutive law is proposed for shear rupture. It is demonstrated that nonuniform distributions of the constitutive law parameters on the fault are necessary for creating the nucleation process, which consists of two phases: (i) a stable, quasistatic phase, and (ii) the subsequent accelerating phase. Physical models of the breakdown zone and the nucleation zone are presented for shear rupture in the brittle regime. The constitutive law for shear rupture explicitly includes a scaling parameter Dc that enables one to give a common interpretation to both small scale rupture in the laboratory and large scale rupture as earthquake source in the Earth. Both the breakdown zone size Xc and the nucleation zone size L are prescribed and scaled by Dc, which in turn is prescribed by a characteristic length lambda c representing geometrical irregularities of the fault. The models presented here make it possible to understand the earthquake generation process from nucleation to unstable, dynamic rupture propagation in terms of physics. Since the nucleation process itself is an immediate earthquake precursor, deep understanding of the nucleation process in terms of physics is crucial for the short-term (or immediate) earthquake prediction.