Elastic encapsulation in bicomponent stratified flow of viscoelastic fluids
Chunfeng Zhou2008, Journal of Rheology
Sign up for access to the world's latest research
checkGet notified about relevant papers
checkSave papers to use in your research
checkJoin the discussion with peers
checkTrack your impact
Abstract
The second normal stress difference N 2 experienced by non-Newtonian fluids flowing in a pipe may give rise to secondary flows in the transverse direction. As a result, one component tends to encapsulate the other in stratified flows. In multilayer coextrusion, such secondary flows tend to distort the interface and affect layer uniformity. This paper presents numerical simulations of the elastically driven encapsulation in two-component stratified viscoelastic fluids. The simulations are based on a phase-field theoretical model and use finite elements with adaptive meshing to resolve the moving interfaces. The results suggest two mechanisms for elastic encapsulation: One due to the mismatch of N 2 between the components and the other due to noncircular geometry of the cross section. In circular pipes, the more elastic fluid tends to protrude into the other component in the center of the pipe and become encapsulated. This produces the curtate cycloid interface shape commonly seen in experiments. If the cross section is noncircular, both the geometric effect and the elastic stratification are at work, and the interfacial motion is determined by the competition of these two mechanisms. This understanding provides an explanation for the anomalous encapsulation of the less elastic component by the more elastic one observed in multilayer coextrusion.
Related papers
Phase-field simulations of interfacial dynamics in viscoelastic fluids using finite elements with adaptive meshing
Journal of Computational Physics, 2006
This paper describes a novel numerical algorithm for simulating interfacial dynamics of non-Newtonian fluids. The interface between two immiscible fluids is treated as a thin mixing layer across which physical properties vary steeply but continuously. The property and evolution of the interfacial layer is governed by a phase-field variable / that obeys a Cahn-Hilliard type of convection-diffusion equation. This circumvents the task of directly tracking the interface, and produces the correct interfacial tension from the free energy stored in the mixing layer. Viscoelasticity and other types of constitutive equations can be incorporated easily into the variational phase-field framework. The greatest challenge of this approach is in resolving the sharp gradients at the interface. This is achieved by using an efficient adaptive meshing scheme governed by the phase-field variable. The finite-element scheme easily accommodates complex flow geometry and the adaptive meshing makes it possible to simulate large-scale two-phase systems of complex fluids. In two-dimensional and axisymmetric three-dimensional implementations, the numerical toolkit is applied here to drop deformation in shear and elongational flows, rise of drops and retraction of drops and torii. Some of these solutions serve as validation of the method and illustrate its key features, while others explore novel physics of viscoelasticity in the deformation and evolution of interfaces.
Finite-element simulation of stratified multiphase flows
AIChE Journal, 1987
The analysis of stratified multiphase flow fields is difficult because the position of the interface is unknown aprioriand there is a discontinuity in the normal viscous stress and/or pressure at the interface. A finite-element technique that uses double nodes along the interface has been implemented. The immiscible liquid displacement in a capillary tube was examined in detail. Fountain flow in the advancing fluid, reverse fountain flow in the receding fluid, and a recirculating flow region in the less viscous fluid were determined. In Newtonian flat film coextrusion the entry region, where the two immiscible fluids form an interface, and the characteristic swelling and bending of the extrudate at the die exit were studied. The computational results compare favorably with available experimental observations.
Nonlinear interfacial instability in two-fluid viscoelastic Couette flow
Journal of Non-Newtonian Fluid Mechanics, 2018
Weakly-nonlinear interfacial instabilities in two-fluid planar Couette flow are investigated for the case where one layer is thin. Taking this thin-layer thickness as a small parameter, asymptotic analysis is used to derive a nonlinear evolution equation for the interface height valid for wavelengths that scale with the channel height. Consequently, the influence of the thick layer is felt through a non-local coupling term which is obtained by solving a system of linear equations which are a simplified viscoelastic analogue to the Orr-Sommerfeld equation. The evolution equation allows for the clear identification of the influence of normal stresses at the interface on both the initial instability and the subsequent nonlinear dynamics. Results from numerical simulations illustrate 1) an array of non-stationary states including traveling waves and chaos, 2) competition between elastic instability and instability due to viscosity stratification, and 3) the accuracy of a simplified 'localized' evolution equation (derived using a long-wave approximation to the coupling term) when either the elasticity of the thick-layer fluid is sufficiently weak or the elasticities of the two fluids are sufficiently well-matched.
Incompressible Multiphase flow and Encapsulation simulations using the moment of fluid method
International Journal for Numerical Methods in Fluids
A moment of fluid method is presented for computing solutions to incompressible multiphase flows in which the number of materials can be greater than two. In this work, the multimaterial moment-of-fluid interface representation technique is applied to simulating surface tension effects at points where three materials meet. The advection terms are solved using a directionally split cell integrated semi-Lagrangian algorithm and the projection method is used to evaluate the pressure gradient force term. The underlying computational grid is a dynamic block structured adaptive grid. The new method is applied to multiphase problems illustrating contact line dynamics, triple junctions, and encapsulation in order to demonstrate its capabilities. Examples are given in 2D, 3D axisymmetric (R-Z), and 3D (X-Y-Z) coordinate systems.
CONSTITUTIVE-RELATED INSTABILITIES IN PIPE FLOW OF A VISCOELASTIC FLUID
2000
The analytical solution for the steady-state flow in a pipe of viscoelastic fluids obeying the complete Phan-Thien--Tanner constitutive equation with a linear stress coefficient is derived. The results include the radial profiles of all relevant stresses, of the axial velocity and of the viscosity. The pipe flow is found to be unstable when the pressure gradient exceeds a critical value
The effect of transient viscoelastic properties on interfacial instabilities in superposed pressure driven channel flows
Journal of Non-Newtonian Fluid Mechanics, 1999
The user has requested enhancement of the downloaded file.
An investigation into the interfacial shear stress contribution in two-phase stratified flow
This paper provides a combined theoretical and experimental investigation into the contribution of interfacial shear stress in certain co and counter-current flows in circular pipes. Based on momentum balance two equations were developed for such flows then two fluid systems of significantly different density ratio were experimentally tested to quantify these equations.
3D phase-field simulations of interfacial dynamics in Newtonian and viscoelastic fluids
Journal of Computational Physics, 2010
This work presents a three-dimensional finite-element algorithm, based on the phase-field model, for computing interfacial flows of Newtonian and complex fluids. A 3D adaptive meshing scheme produces fine grid covering the interface and coarse mesh in the bulk. It is key to accurate resolution of the interface at manageable computational costs. The coupled Navier-Stokes and Cahn-Hilliard equations, plus the constitutive equation for non-Newtonian fluids, are solved using second-order implicit time stepping. Within each time step, Newton iteration is used to handle the nonlinearity, and the linear algebraic system is solved by preconditioned Krylov methods. The phase-field model, with a physically diffuse interface, affords the method several advantages in computing interfacial dynamics. One is the ease in simulating topological changes such as interfacial rupture and coalescence. Another is the capability of computing contact line motion without invoking ad hoc slip conditions. As validation of the 3D numerical scheme, we have computed drop deformation in an elongational flow, relaxation of a deformed drop to the spherical shape, and drop spreading on a partially wetting substrate. The results are compared with numerical and experimental results in the literature as well as our own axisymmetric computations where appropriate. Excellent agreement is achieved provided that the 3D interface is adequately resolved by using a sufficiently thin diffuse interface and refined grid. Since our model involves several coupled partial differential equations and we use a fully implicit scheme, the matrix inversion requires a large memory. This puts a limit on the scale of problems that can be simulated in 3D, especially for viscoelastic fluids.
A general non-isothermal model for one-dimensional multilayer coextrusion flows
Polymer Engineering and Science, 1992
A general computation of multilayer coextrusion flow in a flat die geometry is presented. For any given number of layers of different polymers, characterized by their thermal and rheological behaviors, the model permits computation of velocity and temperature fields along the flow in constant or slightly varying geometries. The influence of different operating parameters (wall regulation temperature, flow rate, initial temperatures) on the interface positions and temperature evolutions is evaluated. Theoretical predictions of interface positions are in agreement with experiments carried out on an industrial multimanifold flat die.
Suppression of purely elastic instabilities in the torsional flow of viscoelastic fluid past a soft solid
Physics of Fluids, 2013
Experiments are performed to explore the role of a soft, deformable solid layer on the purely elastic instability in the torsional flow of polymer solutions between two circular discs. The gel layer is placed on the stationary bottom plate of a rheometer, and the polymer solution is placed between the gel and the rotating top disc. The observed variation of viscosity with shear rate (or shear stress) is correlated with the presence or absence of purely elastic instability in the viscometric flow. Earlier work has shown that with increase in shear rate, the torsional flow of a polymer solution between rigid discs undergoes transition from the simple viscometric flow state to elastic turbulence via a sequence of instability modes. We combine rheological observations and flow visualization to show that the deformable solid has a profound effect on the stability of the torsional flow. In marked contrast to flow between rigid plates (where the fluid shows apparent shear-thickening at the onset of instability), the apparent viscosity continues to decrease up to a much larger value of shear rate with the presence of a soft gel. At a fixed shear rate, for flow past a soft gel, the measured stress does not exhibit marked temporal fluctuations that would otherwise be present without the soft gel. Using flow visualization, we show that secondary flow patterns that form after the instability for a rigid surface disappear for flow on soft gel surfaces. In the case of rigid surfaces, the instability is sub-critical and exhibits hysteresis behavior, which again is absent when the flow occurs past a soft solid layer. Our results show that the role of the soft deformable solid is to suppress the purely elastic instability in torsional flows of polymeric liquids for intermediate shear rates. While it is known that soft deformable solids destabilize the flow of Newtonian liquids in the absence of inertial effects, our study shows that the effect of deformability can be opposite in the torsional flow of viscoelastic liquids.
Related papers
Simulations of interfacial dynamics of complex fluids using diffuse interface method with adaptive meshing
2008
A diffuse-interface finite-element method has been applied to simulate the flow of twocomponent rheologically complex fluids. It treats the interfaces as having a finite thickness with a phase-field parameter varying continuously from one phase to the other. Adaptive meshing is applied to produce fine grid near the interface and coarse mesh in the bulk. It leads to accurate resolution of the interface at modest computational costs. An advantage of this method is that topological changes such as interfacial rupture and coalescence happen naturally under a short-range force resembling the van der Waals force. There is no need for manual intervention as in sharp-interface model to effect such event. Moreover, this energy-based formulation easily incorporates complex rheology as long as the free energy of the microstructures is known. The complex fluids considered in this thesis include viscoelastic fluids and nematic liquid crystals. Viscoelasticity is represented by the Oldroyd-B mode...
A numerical simulation of suspension flow using a constitutive model based on anisotropic interparticle interactions
Rheologica Acta, 2000
We report a Brownian con®guration ®eld implementation of a recent constitutive equation for suspensions, reported by . The numerical method is a hybrid technique, which combines a modi®cation of the Brownian con®guration ®eld method described by and the adaptive viscosity split stress formulation proposed by . The implementation is used to examine the¯ow past a sphere in a tube. The relative viscosity derived from the drag force/ sedimentation velocity agrees well with a well-known empiricism. In addition, the ratio of the pressure force to the drag on the sphere seems to be weakly dependent on the volume fraction, and is somewhat higher than Brenner's results of 1962, which were derived for Newtonian¯uids.
Influence of viscoelasticity on the interfacial dynamics of air displacing fluid flows—a computational study
Journal of Non-Newtonian Fluid Mechanics, 2004
A pseudo-solid domain mapping technique, based on a finite strain model, coupled with a DEVSS finite element formulation is applied to study the effects of viscoelasticity on free surface displacement flows. The flow type analyzed is the displacement of a dilute polymeric solution, modeled by a FENE-CR constitutive equation, in a Hele Shaw Cell. Our study indicates the presence of two distinct flow regimes. In the absence of gravity, a recirculation flow at low Ca (Ca < 1.0) and a bypass flow at high Ca (Ca > 1.0), are observed. In the recirculation flow we observe the formation of elastic stress boundary layers in the capillary transition region (the thickness of which decrease with increasing elasticity), an increase in the hydrodynamic film thickness with increasing elasticity and a meniscus invasion when the stresses in the elastic boundary layer become very large. We have also qualitatively studied the effect of channel wall divergence on the flow dynamics and ascertained that the elastic stress boundary layer formation is a purely local phenomenon and largely independent of geometric considerations. For the bypass flow, in addition to the elastic stress boundary layer in the capillary transition region, an additional elastic stress boundary layer near the bubble tip is observed. As the Wi is increased, film thickening is observed but no meniscus invasion occurs. We have also studied the effect of gravity on both flows. Addition of gravity in recirculation flow results in a film thinning effect and an increase in the normal stresses in the capillary transition regime. For the bypass flow, even a small amount of gravity introduces recirculation. Finally, we discuss how both flows are elastically unstable and plausible mechanisms for the onset of instability are suggested.
Displacement of Viscoelastic Liquids in Capillary Tubes with Oldroyd-B Model
2003
Displacement of a liquid in a capillary tube by gas injection occurs in many situations, like enhanced oil recovery, coating of catalytic converters and gas-assisted injection molding. Generally the liquid being displaced is a polymeric solution or dispersion, that are not Newtonian. Viscoelastic forces alter the force balance in various parts of the flow and consenquently change the amount of liquid left attached to the capillary wall. In order to model the effect of the rheological properties in this important flow, the mechanical behavior of the flowing liquid has to be well described by an appropriate constitutive model. Here, the Oldroyd-B differential constitutive equation that approximate viscoelastic behavior of dilute polymer solutions was used, together with momentum and continuity equation, to model the two-dimensional free surface flow near the gas-liquid interface. The equation system was solved with the Finite Element Method. The resulting non-linear system of algebrai...
A computational study of the coalescence between a drop and an interface in Newtonian and viscoelastic fluids
2006
A drop falling onto a fluid-fluid interface may not merge with it at once but may undergo a so-called partial coalescence cascade. Experimental observations of this phenomenon have revealed fascinating features of the process for Newtonian as well as polymeric fluids. In this paper, we describe numerical simulations of partial coalescence based on a phase-field method. Results show that partial coalescence occurs for an intermediate range of drop sizes, and proceeds in two stages: capillary waves propagating along the drop and transforming it into a fluid column, and neck formation on the column and pinch-off of the secondary drop. In the first stage, interfacial energy turns into kinetic energy following film rupture, while in the second, the kinetic energy overcomes an energy barrier due to the initial increase in interfacial area during neck formation. A parametric study establishes a criterion for partial coalescence in terms of a maximum Ohnesorge number that applies to a wide range of fluid densities and viscosities as long as the Bond number is small. Viscoelasticity in either the drop or the matrix tends to delay the pinch-off of the secondary drop, and may even suppress partial coalescence altogether. The underlying mechanism is large tensile polymer stresses resisting the stretching and thinning of the fluid neck. The numerical results are in qualitative, and in some cases quantitative, agreement with prior experiments.
An experimental study on the effect of polymer viscoelasticity on layer rearrangement in coextruded structures
Polymer Engineering & Science, 1998
Previous work on the layer thickness uniformity of coextruded structures has centered primarily on the effects of differing polymer viscosities in the individual layers. These differences in viscosities result in a phenomenon known as "viscous encapsulation" in which the less viscous layer tends to encapsulate the more viscous layer producing nonuniform layer thicknesses. In this study, the effect of polymer viscoelasticity on the layer thickness uniformity of multilayer coextruded structures was investigated by coextruding multilayer structures through die channels of different cross-sectional shapes and observing the location of the interface. In order to minimize the viscosity effects, the coextrusion experiments were conducted with identical materials in each layer that were pigmented to allow observation of the layer interface. It was shown experimentally that coextruding identical polymers through channels of various geometries can cause the layers to rearrange. This layer rearrangement appears to be caused by secondary flows that occur in different channel geometries due to the viscoelastic characteristics of the polymer. Layer rearrangement of identical polymers implies that even resins with viscosities well matched for coextrusion may experience nonuniform layer thicknesses when coextruded through large dies.
Interfacial stability and deformation of two stratified power law fluids in plane poiseuille flow Part I. Stability analysis
Journal of Non-Newtonian Fluid Mechanics, 1990
The interface stability of two stratified power law fluids in plane Poiseuille flow is considered. Stability diagrams are given in terms of viscosity and depth ratios as a function of power law indices. The work of previous investigators for two Newtonian and second-order fluids has been generalized and presented in a similar fashion. It is found that dependence of viscosity on shear rate has a pronounced effect on the interfacial stability regime in a more drastic manner than that of an effective viscosity change. Furthermore, it is found that for power law fluids the effect of shear thinning viscosity in the less viscous fluid mainly affects the stability region for viscosity ratios less than unity, while for viscosity ratios greater than unity shear thinning in both layers affects the stability regime, particularly at small depth ratios.
Computational modeling of multiphase viscoelastic and elastoviscoplastic flows
International Journal for Numerical Methods in Fluids, 2018
In this paper, a three-dimensional numerical solver is developed for suspensions of rigid and soft particles and droplets in viscoelastic and elastoviscoplastic (EVP) fluids. The presented algorithm is designed to allow for the first time three-dimensional simulations of inertial and turbulent EVP fluids with a large number particles and droplets. This is achieved by combining fast and highly scalable methods such as an FFT-based pressure solver, with the evolution equation for non-Newtonian (including EVP) stresses. In this flexible computational framework, the fluid can be modeled by either Oldroyd-B, neo-Hookean, FENE-P, or Saramito EVP models, and the additional equations for the non-Newtonian stresses are fully coupled with the flow. The rigid particles are discretized on a moving Lagrangian grid, whereas the flow equations are solved on a fixed Eulerian grid. The solid particles are represented by an immersed boundary method with a computationally efficient direct forcing method, allowing simulations of a large numbers of particles. The immersed boundary force is computed at the particle surface and then included in the momentum equations as a body force. The droplets and soft particles on the other hand are simulated in a fully Eulerian framework, the former with a level-set method to capture the moving interface and the latter with an indicator function. The solver is first validated for various benchmark single-phase and two-phase EVP flow problems through comparison with data from the literature. Finally, we present new results on the dynamics of a buoyancy-driven drop in an EVP fluid.
Pore-scale numerical simulation of two phase flow of newtonian and viscoelastic fluids
2011
La récupération assistée d'hydrocarbures par injection de polymères est modéliséeà l'échelle du pore. L'écoulement incompressible de deux fluides immiscibles est simulé dans un réseau de microcanaux. Le mouvement du fluide newtonien est régi par leséquations de Stokes. Le modèle Oldroyd-B est utilisé pour simuler le comportement viscoélastique du fluide non newtonien. L'interface entre les deux fluides est suivieà l'aide de la méthode Level Set. La dynamique du point triple s'appuie sur l'analyse théorique de Cox. Des résultats de simulations réalisées en deux dimensions sont présentés.
Use of Double Nodes in the Numerical Simulation of Multilayered Flows of Immiscible Fluids
The use of double nodes in side-by-side flows of immiscible fluids is introduced. The double nodes are required to satisfy the interface boundary conditions of continuous velocities and stresses. The implementation of the technique in finite element simulations of two viscous immiscible fluids is illustrated through a typical example of co-extrusion.
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.