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 eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids. (C) 2010 Elsevier B.V. All rights reserved.

An implicit technique for solving 3D low Reynolds number moving free surface flows

This paper describes the development of an implicit finite difference method for solving transient three-dimensional incompressible free surface flows. To reduce the CPU time of explicit low-Reynolds number calculations, we have combined a projection method with an implicit technique for treating the pressure on the free surface. The projection method is employed to uncouple the velocity and the pressure fields, allowing each variable to be solved separately. We employ the normal stress condition on the free surface to derive an implicit technique for calculating the pressure at the free surface. Numerical results demonstrate that this modification is essential for the construction of methods that are more stable than those provided by discretizing the free surface explicitly. In addition, we show that the proposed method can be applied to viscoelastic fluids. Numerical results include the simulation of jet buckling and extrudate swell for Reynolds numbers in the range [0.01, 0.5]. (C) 2008 Elsevier Inc. All rights reserved.

Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved

Power-law creep behavior of a semiflexible chain

Rheological properties of adherent cells are essential for their physiological functions, and microrheological measurements on living cells have shown that their viscoelastic responses follow a weak power law over a wide range of time scales. This power law is also influenced by mechanical prestress borne by the cytoskeleton, suggesting that cytoskeletal prestress determines the cell's viscoelasticity, but the biophysical origins of this behavior are largely unknown. We have recently developed a stochastic two-dimensional model of an elastically joined chain that links the power-law rheology to the prestress. Here we use a similar approach to study the creep response of a prestressed three-dimensional elastically jointed chain as a viscoelastic model of semiflexible polymers that comprise the prestressed cytoskeletal lattice. Using a Monte Carlo based algorithm, we show that numerical simulations of the chain's creep behavior closely correspond to the behavior observed experimentally in living cells. The power-law creep behavior results from a finite-speed propagation of free energy from the chain's end points toward the center of the chain in response to an externally applied stretching force. The property that links the power law to the prestress is the chain's stiffening with increasing prestress, which originates from entropic and enthalpic contributions. These results indicate that the essential features of cellular rheology can be explained by the viscoelastic behaviors of individual semiflexible polymers of the cytoskeleton.

Protein-polysaccharide viscoelastic matrices: synergic effects of amylose on lysozyme physical gelation in aqueous dimethylsulfoxide

A synergic effect of amylose on rheological characteristics of lysozyme physical gels evolved out of dimethylsulfoxide-water was verified and analyzed. The dynamics of the gels were experimentally approached by oscillatory rheology. The synergic effect was characterized by a decrease in the threshold DMSO volume fraction required for lysozyme gelation, and by a significant strengthening of the gel structure at over-critical solvent and protein concentrations. Drastic changes in the relaxation and creep curve patterns for systems in the presence of amylose were verified. Complex viscosity dependence on temperature was found to conform to an Arrhenius-like equation, allowing the determination of an activation energy term (Ea, apparent) for discrimination of gel rigidity. A dilatant effect was found to be induced by temperature on the flow behavior of lysozyme dispersions in DMSO-H(2)O in sub-critical conditions for gelation, which was greatly intensified by the presence of amylose in the samples. That transition was interpreted as reflecting a change from a predominant colloidal flow regime, where globular components are the prevailing structural elements, to a mainly fibrillar, polymeric flow behavior.

Experimental analysis of active-passive vibration control using viscoelastic materials and extension and shear piezoelectric actuators

Hybrid active-passive damping treatments combine the reliability, low cost and robustness of viscoelastic damping treatments and the high-performance, modal selective and adaptive piezoelectric active control. Numerous hybrid damping treatments have been reported in the literature. They differ mainly by the relative positions of viscoelastic treatments, sensors and piezoelectric actuators. In this work we present an experimental analysis of three active-passive damping design configurations applied to a cantilever beam. In particular, two design configurations based on the extension mode of piezoelectric actuators combined with viscoelastic constrained layer damping treatments and one design configuration with shear piezoelectric actuators embedded in a sandwich beam with viscoelastic core are analyzed. For comparison purposes, a purely active design configuration with an extension piezoelectric actuator bonded to an elastic beam is also analyzed. The active-passive damping performance of the four design configurations is compared. Results show that active-passive design configurations provide more reliable and wider-range damping performance than the purely active configuration.

DEVELOPMENT OF STATIC MIXERS FOR MISCIBLE FLUIDS IN LAMINAR FLOW WITH THE USE OF COMPUTATIONAL FLUID DYNAMICS (CFD)

Static mixers with improved performance were developed from CFD simulations in a stepwise approach. The relevant geometric features of simple mixer designs and the corresponding mixing mechanisms-laminar shear, elongational flow, and distributive mixing-were identified first. This information was used to formulate guidelines for the development of new geometries. The solid elements of the static mixer should: (a) provide restrictions to the flow; (b) deflect the flow; (c) be sequentially rotated around the flow direction to provide symmetry; (d) extend from the center of the pipe to the vicinity of the walls to avoid short-circuiting; and (e) distribute and remix the flow. Based on these guidelines, two improved mixer designs were developed: the DS A-I mixer has a good mixing efficiency and an acceptable pressure drop; the Fins 35 degrees mixer is more efficient and compact, but requires a larger pressure drop. Their performance indicates that their use is possible on industrial applications.

Validation of a manual headspace gas chromatography method for determining volatile compounds in biological fluids

Background: We report the validation of a method for the determination of acetaldehyde, acetone, methanol, and ethanol in biological fluids using manual headspace sample introduction and an acetonitrile internal standard. Method: This method uses a capillary column (I = 30 m, I.D. = 0.25 mm, dF = 0.25 mu m) installed in a gas chromatography-flame ionization detector (GC-FID) apparatus with a run time of 7.5 minutes. Results: Analysis of the retention times and the resolution of the analyte peaks demonstrated excellent separation without widening of the peaks. Precision and accuracy were good (interassay precision < 15% and recovery between 85% and 115%) in both blood and urine. Conclusion: The method was linear (r > 0.09) over the analytical measurement range (AMR) of each analyte.

Detection of Leptospira spp. in semen and vaginal fluids of goats and sheep by polymerase chain reaction

Thirteen goat herds and seven sheep flocks in the state of Rio de Janeiro, Brazil were screened for leptospirosis. From the three herds and three flocks with greatest seroreactivity, 19 goats (16 females and three bucks) and 40 sheep (26 ewes and 14 rams) that were seropositive (specific anti-Leptospira titres >= 400, based on a microscopic agglutination test), were selected for more detailed studies. From those animals, samples of vaginal fluids or semen were collected for bacteriological and molecular assays. For both species of animals, the most prevalent reactions were to serovars Hardjo, Shermani, and Grippotyphosa. Although leptospires were detected by darkfield microscopy in three vaginal fluid samples (from two goats and one ewe), pure isolates were not obtained by bacteriological culture of vaginal fluids or semen. However, seven vaginal fluid samples (from four goats and three ewes) and six semen samples (all from rams) were positive on polymerase chain reaction (PCR). Based on these findings, in addition to analogous findings in cattle, we inferred that there is potential for venereal transmission of leptospirosis in small ruminants. (c) 2008 Elsevier Inc. All rights reserved.

An evaluation of three upwinding approximations for numerical modeling the flow in tubular membrane of Newtonian and non-Newtonian fluids

In this paper, the laminar fluid flow of Newtonian and non-Newtonian of aqueous solutions in a tubular membrane is numerically studied. The mathematical formulation, with associated initial and boundary conditions for cylindrical coordinates, comprises the mass conservation, momentum conservation and mass transfer equations. These equations are discretized by using the finite-difference technique on a staggered grid system. Comparisons of the three upwinding schemes for discretization of the non-linear (convective) terms are presented. The effects of several physical parameters on the concentration profile are investigated. The numerical results compare favorably with experimental data and the analytical solutions. (C) 2011 Elsevier Inc. All rights reserved.

The competing roles of extensional viscosity and normal stress differences in complex flows of elastic liquids

In various attempts to relate the behaviour of highly-elastic liquids in complex flows to their rheometrical behaviour, obvious candidates for study have been the variation of shear viscosity with shear rate, the two normal stress differences N(1) and N(2) especially N(1), and the extensional viscosity eta(E). In this paper, we shall be mainly interested in `constant-viscosity` Boger fluids, and, accordingly, we shall limit attention to N(1) and eta(E). We shall concentrate on two important flows - axisymmetric contraction flow and ""splashing"" (particularly that which arises when a liquid drop falls onto the free Surface of the same liquid). Modem numerical techniques are employed to provide the theoretical predictions. It is shown that the two obvious manifestations of viscoelastic rheometrical behaviour can sometimes be opposing influences in determining flow characteristics. Specifically, in an axisymmetric contraction flow, high eta(E) , can retard the flow, whereas high N(1) can have the opposite effect. In the splashing experiment, high eta(E) can certainly reduce the height of the so-called Worthington jet, thus confirming some early suggestions, but, again, other rheometrical influences can also have a role to play and the overall picture may not be as clear as it was once envisaged.

Genetic Algorithm for the Determination of Linear Viscoelastic Relaxation Spectrum from Experimental Data

Conventional procedures employed in the modeling of viscoelastic properties of polymer rely on the determination of the polymer`s discrete relaxation spectrum from experimentally obtained data. In the past decades, several analytical regression techniques have been proposed to determine an explicit equation which describes the measured spectra. With a diverse approach, the procedure herein introduced constitutes a simulation-based computational optimization technique based on non-deterministic search method arisen from the field of evolutionary computation. Instead of comparing numerical results, this purpose of this paper is to highlight some Subtle differences between both strategies and focus on what properties of the exploited technique emerge as new possibilities for the field, In oder to illustrate this, essayed cases show how the employed technique can outperform conventional approaches in terms of fitting quality. Moreover, in some instances, it produces equivalent results With much fewer fitting parameters, which is convenient for computational simulation applications. I-lie problem formulation and the rationale of the highlighted method are herein discussed and constitute the main intended contribution. (C) 2009 Wiley Periodicals, Inc. J Appl Polym Sci 113: 122-135, 2009

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.

Magnetic Fluids Based on gamma-Fe(2)O(3) and CoFe(2)O(4) Nanoparticles Dispersed in Ionic Liquids

The disclosure of magnetic ionic liquids (MILs) as stable dispersions of surface modified gamma-Fe(2)O(3) or CoFe(2)O(4) nanoparticles (NPs) in the 1-n-butyl-3-methylimidazolium tetrafluoroborate (BMIBF(4)) ionic liquid is reported. The magnetic NPs were characterized by X-ray powder diffraction, transmission electron microscopy, and Raman spectroscopy. The surface modified NPs have proved to form stable dispersions in BMIBF(4) in the absence of water and behave like a magnetic ionic liquid. The MILs have been characterized by Raman spectroscopy, magnetic measurements, and DSC. The stability of the magnetic NPs in BMIBF(4) is consistently explained by assuming the formation of a semiorganized protective layer composed of supramolecular aggregates in the form of [(BMI)(2)(BF(4))(3)](-). A superparamagnetic behavior and saturation magnetization of ca. 18 emu/g for a sample containing 30% w/w maghemite NPs/BMIBF(4) have been inferred from static and dynamic magnetic measurements. DSC results have shown that the MIL composed of 30% w/w CoFe(2)O(4) NPs/BMIBF(4) remains a liquid phase down to -84 degrees C.