Simulation of particle migration in viscoelastic fluids using the
extended finite element method
Citation for published version (APA):
Choi, Y. J., Hulsen, M. A., & Meijer, H. E. H. (2009). Simulation of particle migration in viscoelastic fluids using the extended finite element method. Poster session presented at Mate Poster Award 2009 : 14th Annual Poster Contest.
Document status and date: Published: 01/01/2009 Document Version:
Publisher’s PDF, also known as Version of Record (includes final page, issue and volume numbers) Please check the document version of this publication:
• A submitted manuscript is the version of the article upon submission and before peer-review. There can be important differences between the submitted version and the official published version of record. People interested in the research are advised to contact the author for the final version of the publication, or visit the DOI to the publisher's website.
• The final author version and the galley proof are versions of the publication after peer review.
• The final published version features the final layout of the paper including the volume, issue and page numbers.
Link to publication
General rights
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain
• You may freely distribute the URL identifying the publication in the public portal.
If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:
www.tue.nl/taverne
Take down policy
If you believe that this document breaches copyright please contact us at:
openaccess@tue.nl
providing details and we will investigate your claim.
Polymer Technology
Simulation of particle migration
in viscoelastic fluids using the
extended finite element method
Y. J. Choi, M. A. Hulsen and H. E. H. Meijer
Materials Technology, Eindhoven University of Technology
/department of mechanical engineering
Introduction
In fluid-particle systems, the apparent properties of suspensions depend on the spatial distribution of particles. Hence, we want to know the evolution of multi-particle structures and the mech-anisms that induce particle migration, especially in polymeric fluids. We aim to develop an eXtended Finite Element Method (XFEM) for the simulation of the flow of viscoelastic fluids with suspended particles.
Numerical Methods
We incorporate an extended finite element method (XFEM) that can decouple the fluid and particle domains completely to capture discontinuities at the interface using a regular (not boundary-fitted) mesh for the whole computational domain in-cluding both the fluid and the particle (Fig. 1(c)).
Figure. 1: Comparison of meshes for different numerical methods.
In the XFEM, the approximation space can be written as
V =Xφi(x)H(s)ai, for all x in the domain,
where aiare degrees of freedom, φiare shape functions and H(s)
is a Heaviside function defined by a scalar level set function s. The XFEM enrichment scheme is shown in Fig. 2 [1].
Figure. 2: The XFEM enrichment for a fluid-particle system.
For moving particles, field variables at the previous time can be undefined near the particle boundary. To overcome this problem, we incorporate a temporary ALE scheme as shown in Fig. 3.
at t=tn at t=tn+1
Eulerian mesh for computation
temporary ALE mesh 1 m 2 m m p m ˆ mesh advects using ( ) where 0 on \ ( ) on ( ) 0 on n n n n n t t P t P t += + ∆ ∇ = Ω = = Γ x x u u u u u m 1 1 Combined Euler : ( ) ( , ) n n nn n n t t + + + − ⋅ = ∆ ∇ ∆+ c u u c f Gc c
Figure. 3: Temporary ALE scheme using semi-implicit Euler.
Results
2 1 1 s s p Geometry - / 4.0 - / 0.4Giesekus constitutive model - mobility parameter =0.01
- viscosity ratio /( ) 0.1
- Deborah number De=
R R a R α β η η η λ = = = + = Ω ⊳ ⊳ x y Ω R1 2 R r a
Figure. 4: Schematic description of a rotating Couette flow.
We consider the motion of a freely suspended particle in rotating Couette flow of a viscoelastic fluid (Fig. 4). The radial migration of a particle released at different initial positions for De=1.0 is shown in Fig. 5. The particle migrates to a stabilized radial po-sition near the outer cylinder regardless of its initial popo-sition. As the Deborah number increases, the rate of migration increases, and the stabilized radial position of the particle shifts toward the outer cylinder (Fig. 6).
* t= Ωt * 1 / r r R = 1 / 1.5 r R= 1 / 2.5 r R= 1 / 3.5 r R=
Figure. 5: Radial position of a particle released at different initial positions for De=1.0.
* t = Ωt t*= Ωt * 0 ( ) a r =3.0 * 0 ( ) b r =3.5 De De De
Figure. 6: The effect of Deborah number on particle migration.
Conclusions
• An extended finite element method has been developed for the simulation of viscoelastic particulate flows. • We find a stabilized radial particle position in rotating
Cou-ette flow for a given Deborah number.
• With increasing Deborah number, the rate of migration increases, and the stabilized radial position shifts toward the outer cylinder.
References:
[1] Y.J. Choi et al., International Journal for Numerical Methods in Engineering,submitted (2009)