Auder Ainser

Delphine Dupuy

Gregory P. Panasenko

Ivan Sirakov

Comptes Rendus Mécanique (CRM)

work,ansion isd for thee tubesmériques,de.ion;rtielle duC. R. Mecanique 331 (2003) 609–615Flow in wavy tube structure:asymptotic analysis and numerical simulationAuder Ainsera, Delphine Dupuyb, Gregory P. Panasenkob,c, Ivan Sirakovaa Laboratoire de rhéologie des matières plastiques, Université de Saint-Etienne, 23, rue Paul Michelon, 42023 Saint-Etienne, Franceb Équipe d’analyse numérique, UPRES EA 3058, Université de Saint-Etienne, 23, rue Paul Michelon, 42023 Saint-Etienne, Francec Laboratoire de modélisation en mécanique-CNRS UMR 7607,Université Pierre et Marie Curie-Paris 6, 8, rue du Capitaine Scott, 75015 Paris, FranceReceived 25 May 2003; accepted 6 June 2003Presented by Évariste Sanchez-PalenciaAbstractThis paper deals with the study of the stationary, incompressible, 2D flow of a fluid in a thin wavy tube. In thiswe consider a domain which is the union of two wavy tubes depending on a small parameter. The asymptotic expconstructed. The method of partial asymptotic decomposition is applied. The numerical implementation of this methoextrusion process is developed. The new physical effects are discussed.To cite this article: A. Ainser et al., C. R. Mecanique331 (2003). 2003 Académie des sciences. Published by Éditions scientifiques et médicales Elsevier SAS. All rights reserved.RésuméÉcoulement dans une structure tubulaire ondulée : analyse asymptotique et résolution numérique. Nous considéronsici le mouvement bi-dimensionnel et stationnaire d’un fluide incompressible à l’intérieur d’un domaine constitué dondulés. La méthode de décomposition asymptotique partielle du domaine est mise en place et des résultats nuobtenus pour la modélisation de procédés d’extrusion seront présentés afin de justifier l’application de cette méthoPourciter cet article : A. Ainser et al., C. R. Mecanique 331 (2003). 2003 Académie des sciences. Published by Éditions scientifiques et médicales Elsevier SAS. All rights reserved.Keywords: Computational fluid mechanics; Asymptotic expansion; Stokes equation; Partial asymptotic domain decomposition; ExtrusNumerical solutionMots-clés : Mécanique des fluides numérique ; Dévéloppement asymptotique ; Équation de Stokes ; Décomposition asymptotique padomaine ; Extrusion ; Solution numériqueE-mail addresses: ainser@univ-st-etienne.fr (A. Ainser), delphine.dupuy@univ-st-etienne.fr (D. Dupuy),Grigory.Panasenko@univ-st-etienne.fr (G.P. Panasenko), yvan.sirakov@univ-st-etienne.fr (I. Sirakov).1631-0721/$ – see front matter 2003 Académie des sciences. Published by Éditions scientifiques et médicales Elsevier SAS. All rightsreserved.doi:10.1016/S1631-0721(03)00131-1610 A. Ainser et al. / C. R. Mecanique 331 (2003) 609–615x tubesn utilisantéparer lejustifierge des’intérieurprocédéIn thisby thinymptoticer sub-APDDw in as as welles of the. As arackingonfirmededVersion française abrégéeCette analyse consiste à étudier l’écoulement d’un fluide à l’intérieur d’un domaine constitué de deuondulés dépendant d’un petit paramètreε > 0.Le problème de Stokes dans une structure tube constituée de cylindres fins a été considéré dans [3] ela méthode de décomposition asymptotique partielle du domaine (MAPDD). Cette méthode consiste à sproblème initial en petits sous-problèmes afin de réduire le coût total de la solution numérique.Tout d’abord, il nous faut construire une solution asymptotique pour le problème afin de décrire et del’application de la MAPDD. Cette analyse confirme la localisation des effets de couche limite au voisinazones de transition ainsi que la convergence de la solution asymptotique vers une solution périodique à ldes tubes. La justification numérique proposée ici, est l’application de cette méthode pour simuler und’extrusion de polymère.1. IntroductionThis paper deals with the study of the stationary, incompressible, 2D flow of a fluid in a thin wavy tube.work, we consider a domainGε which is the union of two wavy tubes depending on a small parameterε > 0.Flows in thin layers have been considered in [1,2]. The Stokes problem in tube structure constitutedcylinders has been studied in [3]. The same method is applied in the present study. The proposed partial asdomain decomposition method (MAPDD), described in [4], splits the initial Stokes problem in some smallproblems and essentially reduces the total numerical cost of the total numerical solution.In the first section, we build and justify an asymptotic solution for the Stokes problem. In Section 2, the Mis described and justified. In Section 3, the MAPDD method is validated numerically by analyzing the flowavy tube. This analysis confirmed the localized character of the transition zone between the periodic flowas the rapid convergence of the asymptotic solution toward the periodic one. However, the main advantagproposed method are clearly visible for big 3D problems, where the full solution is often difficult to obtainfinal application, we use the MAPDD method for analyzing the polymer extrusion process.The detailed description of the flow is presented in terms of velocity and pressure distribution, particle tand residence time analysis. The calculations predict the existence of a backward flow which has to be cexperimentally.2. Stokes equations in two wavy tubes2.1. ProblemConsider three setsYi ⊂ [0,1] × (0,1/2] (0  i  2) defined byYi :={(y1, y2) ∈ R2: 0< y1 < 1 and−h+i (y1) < y2 < h+i (y1)}for 0  i  2whereh±i :R → (0,1/2) are functions of classC2(R) and such that there existsm: (∀)t ∈ R, h±i (t)  m> 0 (for0  i  2). We assume that the functionsh±1 andh±2 are 1-periodic.For the boundary, we introduce the notationΓ ±i := {(y1,±h±i (y1)): 0< y1 < 1} for 0  i  2 and we assumthat the curves given by(Γ ±1 − e1)∪ Γ ±0 ∪ (Γ ±2 + e1) are of classC2.Let ε > 0 be a small parameter withε = 1n, n ∈ N.We define the domainGε0, obtained fromY0 by a homothetic contraction in 1/ε times (with respect to 0) anthe domainsGε1 := (⋃nl=1 ε{Y1 − le1})′ andGε2 := (⋃n−1l=1 ε{Y2 + le1})′, where for any setA, Ā is the closureA. Ainser et al. / C. R. Mecanique 331 (2003) 609–615 611oftory.fatemi-nandA′ is the set of the interior points. In this study, we consider the thin domainGε := Gε1 ∪ Gε0 ∪ Gε2. LetΣεi := Gε ∩ {x1 = i} (i = −1,1) be the lateral sides of the domain andΓ ±ε be the upper and lower parts∂Gε \ {Σε−1 ∪Σε1}.The stationary flow of a viscous incompressible fluid through the channelGε is governed by the followingStokes system:{−νvε + ∇pε = f, div vε = 0 inGεvε = 0 onΓ ±ε , vε = ε2ϕε onΣεi , i = −1,1(1)Let H 1div(Gε) be the space of vector valued functions from[H 1(Gε)]2 with vanishing divergence and equalzero on the boundariesΓ ±ε andV (Gε) the subspace ofH 1div(Gε) of functions vanishing on the whole boundaWe denote byWi the space of divergence free functionsu ∈ C∞per((⋃l∈Z(Yi + le1))′), whereC∞per is the space oC∞ functions which are 1-periodic iny1; the closure ofWi with respect to the usual norm of[H 1(Yi)]2 is denotedby W(Yi).Let the functionϕε be such thatϕε = ϕi (x/ε) onGεi with ϕi ∈W(Yi) (for i = 1,2). Moreover, we assume thh+1 (0)∫−h−1 (0)(ϕ1)1(0, y2)dy2 =h+2 (0)∫−h−2 (0)(ϕ2)1(0, y2)dy2 := κ (2)where(·)1 stands the first component.For f ∈ [H−1(Gε)]2 andϕε satisfying (2), problem (1) has a unique solution(vε,pε) ∈ H 1div(Gε)×L2(Gε)/R(see [5]).2.2. Asymptotic expansionIn the following, we suppose thatf(x) = (f1(x1),0)t with f1 a known function ofL2(−1,1). We look for anasymptotic solution in the formvεa(x)= ε2v1(xε) + vbl−1( xε + 1ε e1) + vbl0( xε ) if x ∈Gε1vbl0(xε)if x ∈ Gε0v2(xε) + vbl0( xε ) + vbl1( xε − 1ε e1) if x ∈Gε2(3)For i = 1,2, vi (y)= wαi (y), the unique solution to the following problem: find(wαi ,παi ) ∈ W(Yi)×L2(Yi)/Rsuch that−νwαi + ∇παi = αe1 in (⋃l∈Z (Yi + le1))′, whereα is such that∫ h+i (0)−h−i (0)(wαi )1(0, y2)dy2 = κ .The boundary layers termsvbl−1 andvbl1 are the exponentially decaying solutions of Stokes systems in sinfinite channelΩ1 := (⋃∞l=1(Y1 + le1))′ and Ω2 := (⋃∞l=1(Y2 − le1))′ respectively, with boundary conditiovbl−1 + v1 = ϕ1 on Ω1 ∩ {y1 = 0} andvbl1 + v2 = ϕ2 on Ω2 ∩ {y1 = 0} (see [6]).Denote byΩ = ([⋃∞l=1(Y1 − le1)] ∪ Y0 ∪ [⋃∞l=1(Y2 + le1)])′. The boundary layervbl0 ∈ H 1div(Ω) is anexponentially decaying function which satisfies the Stokes equations with a right-hand side equal tof1e1 onY0.Functionvεa belongs to[H 1div(Gε)]2 and satisfies the boundary conditionvεa = 0 on Γ ±ε and vεa(±1, x2) =ε2[ϕε + vbl0(±1/ε, x2/ε)] onΣε±1. However, it can be modified by adding a functionDε ∈ [H 1div(Gε)]2 such thatDε + vεa = ε2ϕε on ∂Gε, Dε = 0 if −1 + ε < x1 < 1 − ε and‖Dε‖H1(Gε)  ε2cexp(−λ/ε), wherec andλ arepositive constants which do not depend onε.DenoteVεa = vεa + Dε. FunctionVεa is a solution to the following problem:{Vεa ∈ [H 1div(Gε)]2 such thatVεa − ε2ϕε ∈ V (Gε), andν∫Gε ∇Vεa · ∇w =∫Gε f · w + ν∫Gε ∇Dε · ∇w (∀)w ∈ V (Gε)(4)612 A. Ainser et al. / C. R. Mecanique 331 (2003) 609–615cessStokes.but withlymerd twine realenoughl basedting (orcult.crew isframe ofwhichbe takenglectedThe following estimate holds∥∥vε − Vεa∥∥[H1(Gε)]2  cε2 exp(−λ/ε) (5)wherec andλ are two positive constants independent ofε.2.3. Method of asymptotic partial decomposition of domainIntroduce another parameterδ = Kε[| lnε|], with some finiteK ∈ N∗.Let V δ,εdec(Gε) be the spaceVδ,εdec(Gε) := {u ∈ V (Gε): u(x)= ε2[v1(xε)− ϕε(x)]if − 1+ δ < x1 <−δ,u(x)= ε2[v2(xε)− ϕε(x)]if ε + δ < x1 < 1− δ}(6)We define the partially decomposed problem as a variational problem (4) stated on the restricted spaceV δ,εdec(Gε){Find vεd ∈ H 1div(Gε)such thatvεd − ε2ϕε ∈ V δ,εdec(Gε) andν∫Gε∇vεd · ∇w =∫Gεf1e1 · w (∀)w ∈ V δ,εdec(Gε)(7)As in [4], we prove that for allJ ∈ N there existK ∈ N andc ∈ R+ such that∥∥vε − vεd∥∥[H1(Gε)]2  cεJ (8)wherec is a constant which does not depend onε.This estimate shows that the partially decomposed solutionvεd is exponentially close to the solutionvε of theinitial problem (1).The same result holds true in the three-dimensional case, whenYi are the domains{(y1, y2, y3) ∈ R3: 0< y1 <1, −h−(ϕ, y1) < r < h+(ϕ, y1)}, wherer andϕ are the polar coordinates iny2y3-plane andh± are 2π -periodic inϕ and positive functions ofC2(R × [0,1]). The below numerical experiment is developed for an extrusion proin such three-dimensional wavy tube. Another but similar approach was used in [7] to reduce the linearproblem in an infinite domain with some outlets to the Stokes problem set in a truncated bounded domainEstimate (8) is generalized to the case of a wavy tube structure, analogous to the tube structure of [3]Gε1 type domains instead of the thin cylinders of [3].3. Application to the polymer extrusionThe polymer extrusion is an important part of industrial polymer processing. It is mainly used in poforming, blending and mixing. To perform this operation, two main engineering solutions exist: single anscrew extruder. We will consider a model problem with the geometry of the domain that is a part of thtwin screw geometry (Fig. 1), where the main advantages are: the possibility to process polymers with noadhesion at the barrel wall; improved mixing and residence time distribution; development of new materiaon the reactive extrusion. The twin screw extruder contains two intermeshing (or not intermeshing) co-rotacounter-rotating) screws; this geometry is time dependent and therefore the numerical analysis is very diffiTo simplify the problem, we neglect the intermeshing area and therefore only the flow inside one sconsidered as a model problem. An additional advantage of this configuration is that we can change thereference by considering the rotating cylindrical barrel (instead of the screw) with opposite angular velocitymakes the geometry stationary. In this new frame of reference, the centrifugal and Corriolis forces shouldin account, but for high viscosity polymers these inertial forces are not important and generally could be neA. Ainser et al. / C. R. Mecanique 331 (2003) 609–615 613ibed bytonianf theproposesundaryflow andem inre thosew ratemationsvelocitys is doneFig. 1. Twin screw geometry.Fig. 2. Dimensionless pressure and velocity.(as well as all the inertial effects). We should mention that the rheological behavior of polymers is descrcomplex differential or integral viscoelastic constitutive models but in this study we consider them as a Newfluid with constant viscosity which is normalized to 1 in the calculations.Nowadays, it is difficult to solve practical 3D problems in entire flow domain due to complexity ogeometry and excessive memory requirements. The above asymptotic domain decomposition techniquean alternative approach by splitting the flow domain and solving smaller problems joined by appropriate boconditions. For each of two domains, we suppose that the fully developed periodic flow is reached at the inoutflow sections. The velocity distribution of this periodic flow is calculated by solving a small 3D cell problorder to impose velocity field at the inlet and outlet of each sub-domain. The dimensions of the problem aof the twin screw extruder at the laboratory (Clextral BC21). The operating conditions are as follows: flo(700 mm3/s) and screw angular velocity (31.5 rad/s).The problem in each sub-domain is discretized in terms of classical Galerkin formulation. These approxiare relatively expensive for 3D problems but they provided improved accuracy for the pressure and thefields and therefore they were implemented in our code. The resolution of the linearized system of equation614 A. Ainser et al. / C. R. Mecanique 331 (2003) 609–615Fig. 3. Problem domains, finite element mesh, velocity distribution and typical streamline.Fig. 4. Dimensionless pressure distribution and streamline from the back flow.A. Ainser et al. / C. R. Mecanique 331 (2003) 609–615 615tion ofhe finitelly aionhandednnel ofd screws whichs. Theesentedccordingation.3 01.07.odelsr Stokes.for Laméstems,by using the iterative solver based on the bi-conjugate stabilized algorithm of van der Vorst [8] in combinaLU incomplete factorization based on the drop tolerance.In Fig. 3, we present the geometry of the problem together with three sub-domains of decomposition, telement mesh, the extrusion velocity (Vz) distribution at two cross sections and one cutting plane and finatypical streamline with lifetime of about 9 s. The calculatedVz profile proves the periodic character of the solutfar from the three transition zones (Figs. 3 and 2). Due to the difference of the pressure drop in the rightscrew (positive) and in the left handed screw (negative), the pressure driven flow is negative in the chathe right handed screw (Fig. 2 and Fig. 3 cross-section 1) and positive in the channel of the left hande(Fig. 2 and Fig. 3 cross-section 2). From our graphics, we can estimate the size of the transition zoneis for example about 6 mm from the center of the zone 2 (Fig. 3) and corresponds to 1/4 of screw pitch. Thestreamline shape is in agreement with the well known helical motion inside the channels of the screwdimensionless pressure distribution together with another (more exotic) streamline (44 s lifetime) are prin Fig. 4. These computations have proved the existence of a back flow in the right handed screw which, ato our knowledge, has never been mentioned in the literature and is waiting to find its experimental confirmAcknowledgementsThis work was supported by the Region Rhône-Alpes grants, projects nos. 02 020611 01 and 02 02061References[1] S.A. Nazarov, Asymptotic solution of the Navier–Stokes problem on the flow in thin layer fluid, Siberian Math. J. 31 (1990) 296–3[2] A. Bourgeat, E. Marusić-Paloka, Nonlinear effects for flow in periodically constricted channel caused by high injection rate, Math. MMethods Appl. Sci. 8 (3) (1998) 379–405.[3] F. Blanc, O. Gipouloux, G.P. Panasenko, A.M. Zine, Asymptotic analysis and partial asymptotic decomposition of domain foequation in tube structure, Math. Models Methods Appl. Sci. 9 (9) (1999) 1351–1378.[4] G.P. Panasenko, Asymptotic partial decomposition of variational problems, C. R. Acad. Sci. Paris, Sér. IIb 327 (1999) 1185–1190[5] R. Temam, Navier–Stokes Equations, North-Holland, 1984.[6] G.P. Galdi, An Introduction to the Mathematical Theory of the Navier–Stokes Equations, Springer, 1984.[7] S.A. Nazarov, M. Specovius-Neugebauer, Approximation of unbounded domains by bounded domains. Boundary value problemsoperator, St. Petersburg Math. J. 8 (5) (1997) 879–912.[8] H.A. van der Vorst, Bi-CGSTAB, A fast and smoothly converging variant of Bi-CG for the solution of nonsymmetrical linear sySIAM J. Sci. Statist. Comput. 13 (1992) 631–644.

Flow in wavy tube structure: asymptotic analysis and numerical simulation

https://comptes-rendus.academie-sciences.fr/mecanique/item/10.1016/S1631-0721(03)00131-1.pdf

Flow in wavy tube structure: asymptotic analysis and numerical simulation

Abstract

Similar works

Full text

Available Versions

Comptes Rendus Mécanique (CRM)