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

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

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Solitary waves on a free surface of a heated Maxwell fluid

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Coeficientes de perda de carga em válvulas e acessórios no escoamento de fluídos não-newtonianos puramente viscosos e viscoelásticos

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

Estudos teóricos e numéricos de escoamentos com escorregamento

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

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.

Creep and recovery of magnetorheological fluids: experiments and simulations

A direct comparative study on the creep-recovery behavior of conventional MR fluids is carried out using magnetorheometry and particle-level simulations. Two particle concentrations are investigated (ϕ=0.05 and 0.30) at two different magnetic field strengths (53 kA•m-1 and 173 kA•m-1) in order to match the yield stresses developed in both systems for easier comparison. Simulations are mostly started with random initial structures with some additional tests of using preassembled single chains in the low concentration case. Experimental and simulation data are in good qualitative agreement. The results demonstrate three regions in the creep curves: i) In the initial viscoelastic region, the chain-like (at ϕ=0.05) or percolated three-dimensional network (at ϕ=0.30) structures fill up the gap and the average cluster size remains constant; ii) Above a critical strain of 10 %, in the retardation region, these structures begin to break and rearrange under shear. At large enough imposed stress values, they transform into thin sheet-like or thick lamellar structures, depending on the particle concentration; iii) Finally in the case of larger strain values either the viscosity diverges (at low stress values) or reaches a constant low value (at high stress values), showing a clear bifurcation behavior. For stresses below the bifurcation point the MR fluid is capable to recover the strain by a certain fraction. However, no recovery is observed for large stress values.

Processing and hygrothermal effects on viscoelastic behavior of glass fiber/epoxy composites

Fiber reinforced epoxy composites are used in a wide variety of applications in the aerospace field. These materials have high specific moduli, high specific strength and their properties can be tailored to application requirements. In order to screening optimum materials behavior, the effects of external environments on the mechanical properties during usage must be clearly understood. The environmental action, such as high moisture concentration, high temperatures, corrosive fluids or ultraviolet radiation (UV), can affect the performance of advanced composites during service. These factors can limit the applications of composites by deteriorating the mechanical properties over a period of time. Properties determination is attributed to the chemical and/or physical damages caused in the polymer matrix, loss of adhesion of fiber/resin interface, and/or reduction of fiber strength and stiffness. The dynamic elastic properties are important characteristics of glass fiber reinforced composites (GRFC). They control the damping behavior of composite structures and are also an ideal tool for monitoring the development of GFRC's mechanical properties during their processing or service. One of the most used tests is the vibration damping. In this work, the measurement consisted of recording the vibration decay of a rectangular plate excited by a controlled mechanism to identify the elastic and damping properties of the material under test. The frequency amplitude were measured by accelerometers and calculated by using a digital method. The present studies have been performed to explore relations between the dynamic mechanical properties, damping test and the influence of high moisture concentration of glass fiber reinforced composites (plain weave). The results show that the E' decreased with the increase in the exposed time for glass fiber/epoxy composites specimens exposed at 80 degrees C and 90% RH. The E' values found were: 26.7, 26.7, 25.4, 24.7 and 24.7 GPa for 0, 15, 30, 45 and 60 days of exposure, respectively. (c) 2005 Springer Science + Business Media, Inc.

Influence of viscoelastic characteristics of the colloidal suspensions on the formation of SnO2 supported membranes

SnO2 supported membranes have been prepared by sol-casting on alumina tubular substrate, using aqueous colloidal suspensions prepared by sol-gel route. The viscoelastic behaviour during sol ageing was analyzed by dinamic rheologial measurements. The complex viscosity and the storage and loss moduli have been followed during the sol-gel transition and the results have been correlated with the linear aggregates growth and the scalar percolation models. The scanning electron microscopy has evidenced that the homogeneity and thickness of the membrane depend on the sol ageing time. Crack-free and homogeneous membranes have been obtained for coated layers of 0.5μm thickness. © 1997 Trans Tech Publications.

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.

Numerical study of ultrasound induced non-linear shape and size bubble oscillations in viscoelastic media

The theoretical study of forced bubble oscillations is motivated by the importance of cavitation bubbles and oscillating encapsulated microbubbles (i.e. contrast agents) in medical sciences. In more details,theoretical studies on bubble dynamics addressing the sound-bubble interaction phenomenon provide the basis for understanding the dynamics of contrast agent microbubbles used in medical diagnosis and of non-linearly oscillating cavitation bubbles in the case of high-intensity ultrasound therapy. Moreover, the inclusion of viscoelasticity is of vital importance for an accurate theoretical analysis since most biological tissues and fluids exhibit non-Newtonian behavior.

Capillary break-up rheometry of low-viscosity elastic fluids

We investigate the dynamics of the capillary thinning and break-up process for low viscosity elastic fluids such as dilute polymer solutions. Standard measurements of the evolution of the midpoint diameter of the necking fluid filament are augmented by high speed digital video images of the break up dynamics. We show that the successful operation of a capillary thinning device is governed by three important time scales (which characterize the relative importance of inertial, viscous and elastic processes), and also by two important length scales (which specify the initial sample size and the total stretch imposed on the sample). By optimizing the ranges of these geometric parameters, we are able to measure characteristic time scales for tensile stress growth as small as 1 millisecond for a number of model dilute and semi-dilute solutions of polyethylene oxide (PEO) in water and glycerol. If the final aspect ratio of the sample is too small, or the total axial stretch is too great, measurements are limited, respectively, by inertial oscillations of the liquid bridge or by the development of the well-known beads-on-a-string morphology which disrupt the formation of a uniform necking filament. By considering the magnitudes of the natural time scales associated with viscous flow, elastic stress growth and inertial oscillations it is possible to construct an operability diagram characterizing successful operation of a capillary break-up extensional rheometer. For Newtonian fluids, viscosities greater than approximately 70 mPas are required; however for dilute solutions of high molecular weight polymer, the minimum Viscosity is substantially lower due to the additional elastic stresses arising from molecular extension. For PEO of molecular weight 2.10(6) g/mol, it is possible to measure relaxation times of order 1 ms in dilute polymer solutions with zero-shear-rate viscosities on the order of 2-10 mPas.

Simulations of pendant drop formation of a viscoelastic liquid

A modified Volume-of-Fluid (VOF) numerical method is used to predict the dynamics of a liquid drop of a low viscosity dilute polymer solution, forming in air from a circular nozzle. Viscoelastic effects are represented using an Oldroyd-B model. Predicted drop shapes are compared with experimental observations. The main features, including the timing of the shape evolution and the bead-on-a-string effect, are well reproduced by the simulations. The results confirm published conclusions of the third author, that the deformation is effectively Newtonian until near the time of Newtonian pinch-off and that the elastic stress becomes large in the pinch region due to the higher extensional flow there.