• No results found

A model of rotary spinning process

N/A
N/A
Protected

Academic year: 2021

Share "A model of rotary spinning process"

Copied!
9
0
0

Bezig met laden.... (Bekijk nu de volledige tekst)

Hele tekst

(1)

A model of rotary spinning process

Citation for published version (APA):

Hlod, A., Ven, van de, A. A. F., & Peletier, M. A. (2010). A model of rotary spinning process. (CASA-report; Vol. 1065). Technische Universiteit Eindhoven.

Document status and date: Published: 01/01/2010 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.

(2)

EINDHOVEN UNIVERSITY OF TECHNOLOGY

Department of Mathematics and Computer Science

CASA-Report 10-65

October 2010

A model of rotary spinning process

by

A. Hlod, A.A.F. van de Ven, M.A. Peletier

Centre for Analysis, Scientific computing and Applications

Department of Mathematics and Computer Science

Eindhoven University of Technology

P.O. Box 513

5600 MB Eindhoven, The Netherlands

ISSN: 0926-4507

(3)

A model of rotary spinning process

A. Hlod, A.A.F. van de Ven, and M.A. Peletier

Abstract A rotary spinning process is used to produce aramide fibers. In this process

thin jets of polymer solution emerge from the nozzles of the rotating rotor and flow towards the cylindrical coagulator. At the coagulator the jets hit the water curtain in which they solidify forming fibers. The rotary spinning is described by a steady jet of viscous Newtonian fluid between the rotor and the coagulator. The jet model includes the effects of inertia, longitudinal viscosity, and centrifugal and Coriolis forces. For the jet model the specific type of the boundary conditions depends on the balance between the inertia and viscosity in the momentum transfer through the jet cross-section. Based on that we find two possible flow regimes in rotary spinning: 1) viscous-inertial, where viscosity dominates at the rotor and inertia at the coagulator, 2) inertial, where inertia dominates everywhere in the jet. Moreover, there are two situations where spinning is not possible, either due to lack of a steady-jet solution or because the jet wraps around the rotor. Finally, we characterize the parameter space.

1 Rotary spinning process

A rotary spinning device consists of a rotor and a coagulator [2, 4]. Both the coag-ulator and the rotor have the form of a vertical cylinder. A water curtain falls along

A. Hlod

Dept. of Mathematics and Computer Science, Eindhoven University of Technology, PO Box 513 5600 MB Eindhoven The Netherlands, e-mail: avhlod@gmail.com

A.A.F. van de Ven

Dept. of Mathematics and Computer Science, Eindhoven University of Technology, PO Box 513 5600 MB Eindhoven The Netherlands, e-mail: A.A.F.v.d.Ven@tue.nl

M.A. Peletier

Dept. of Mathematics and Computer Science, Eindhoven University of Technology, PO Box 513 5600 MB Eindhoven The Netherlands, e-mail: M.A.Peletier@TUE.nl

(4)

2 A. Hlod, A.A.F. van de Ven, and M.A. Peletier

the coagulator’s inner wall; see Figure 1. Inside the coagulator the rotor is placed so

Polymer solution Rotor

Coagulator

Water curtain

Washing & drying

Fig. 1 Rotary spinning process

that the symmetry axes of the rotor and the coagulator coincide. The rotor rotates counterclockwise and has small nozzles in its lateral surface. Hot polymer solution is pumped through the rotor’s nozzles, flows to the coagulator under the influence of Coriolis and centrifugal forces, and hits the water curtain at the coagulator wall. The resulting fiber is transported away by the water, then it is washed and cut into small pieces to get pulp.

The rotary spinning problem was presented at the 48th European Study Group Mathematics with Industry in Delft (2004) [2]. To describe the jet, the string model was used. At that time, the equations were not solved due to the assumption that the jet always leaves the nozzle radially. The second study of the rotary spinning process was done in [4, 5]. It has been shown there that the jet orientation at the nozzle is determined by the jet itself. However, understanding why the jet orientation can not be prescribed at the nozzle was missing.

In this paper we present the analysis of the rotary spinning model [2, 4, 5] with the boundary conditions derived in [3, Chapter 2], and characterize the complete parameter space. Similar configurations without coagulator are studied in [1, 6–14].

2 Model

We model the rotary spinning process by considering the steady jet (see Figure 2) in two dimensions (neglecting the vertical motion due to gravity). The jet moves in a fixed horizontal frame from the nozzle of the rotor to the contact point with

(5)

A model of rotary spinning process 3

the coagulator. In a fixed reference frame the rotor rotates counter-clockwise with

se Rrot Rcoags R(s) vnozzleRcoag s= 0 φ(s) β(s) FC Fc ex ey

Fig. 2 A schematic picture of the rotary spinning process.

angular velocityΩ. The radii of the rotor and the coagulator are Rrot and Rcoag, respectively. The jet is parameterized by its length s with s= 0 at the nozzle and s= sendat the contact with the coagulator. Note that the jet length sendis unknown in advance. The flow velocity at the nozzle is vnozzle, and at the contact with the coagulator the jet sticks to it having the flow velocityΩRcoag. To describe the jet position in the rotating reference frame of the rotor, we use two sets of coordinates: either the polar coordinates R andβ, with the origin at the center of the rotor, or the arc length of the jet s, and the angleφ the tangent to the jet makes with the radial direction; see Figure 2. In the rotating reference frame, two (inertial) body forces act on the jet, i.e the centrifugal, Fc, and Coriolis, FC, force.

The system describing the stationary jet follows from the conservations of mass and momentum. We scale the flow velocity v with respect toRcoag, both R and

s with respect to Rcoag. Then the system is fully described by three dimensionless

numbers B= 3ν/(ωR2

coag), R0= Rrot/Rcoag, Dr= vnozzle/(ΩRcoag), whereνis the

kinematic viscosity of the fluid. The resulting system is

ξ′(s) = cos(φ(s))R(s)/v(s), (1)

ξ(s)φ′(s) = −R(s)sin(φ(s))/v(s) − sin(φ(s))ξ(s)/R(s) + 2, (2)

v(s) = (v(s)2

−ξ(s)v(s))/B, (3)

R(s) = cos(φ(s)),β′(s) = − sin(φ(s))/R(s), (4)

together with the boundary conditions

v(0) = Dr, v(send) = 1, R(0) = R0,β(0) = 0, R(send) = 1,

(6)

4 A. Hlod, A.A.F. van de Ven, and M.A. Peletier

Here, v(s) is the flow velocity in the jet, ξ(s) represents the momentum transfer

through the jet cross-section. The boundary conditions for the jet orientationφ(s)

are determined by the sign ofξ(s); see [3, Chapter 2].

From (1) it follows thatξ(s) is a strictly increasing function implying that there

are three possibilities for the sign ofξ(s):

• The first situation isξ(s) > 0 everywhere in the jet, and then inertia dominates

everywhere in the momentum transfer through the jet cross-section. We call this jet flow regime inertial. In this case the jet must be aligned with the radial nozzle direction,φ(0) = 0.

• In the second situationξ(s) changes sign from negative to positive, and viscosity

dominates near the nozzle and inertia near the coagulator in the momentum trans-fer through the jet cross-section. We call this jet flow regime viscous-inertial. In this case we can not prescribe any boundary condition for the jet orientation. However, from (2), it follows that at the point s0whereξ(s0) = 0 the jet should be aligned with the direction of the resulting external force acting on the jet, yielding

φ(s0) = arcsin(2v(s0)/R(s0)). (5)

• In the third situation ξ(s) < 0 everywhere in the jet, and viscosity dominates

everywhere in the momentum transfer through the jet cross-section. We call this jet flow regime viscous. For viscous jet we require tangency with the coagulator

φ(send) =π/2. However, the viscous jet situation is not possible in the current

setup because of the following argument.

The border between the parameter regions for the viscous and viscous-inertia jets should satisfy the conditionξ(send) = 0. Then from (5) and the boundary conditions it follows that sin(φ(send)) = 2 leading to a contradiction.

The following solution strategy is suggested. For the equations (1)-(2) we find a first integral

sin(φ(s))ξ(s) =R(s)

2

− C1

R(s) ,

where C1= R0for the inertial jet, C1= R(s0) for the viscous-inertial jet. Then sys-tem is solved using the shooting method.

The parameter space is described by the three dimensionless numbers R0, Dr, and B. In Figure 3 we present the partitioning of the parameter space in the R0, Dr-plane for different B. In this Dr-plane we observe the regions of inertial jet, viscous-inertial jet, and two regions of nonexistence of a jet solution. For the parameter regions ”no solution 1” the cause of nonexistence is that the jet does not reach the coagulator and wraps round the rotor. For the parameter regions ”no solution 2” the jet reaches the coagulator, but the flow velocity at the coagulator can not be matched, indicating the unsteady jet. In the region ”no solution” both causes of nonexistence are possible. The borders between the viscous-inertial and inertial regimes are obtain from the conditionξ(0) = 0 and the borders of the nonexisting

regions are obtained using monotonicity properties of the jet solution together with the conditionφ(send) =π/2.

(7)

A model of rotary spinning process 5 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 D r R0 B= 0.01 inertial jet viscous-inertial jet no solution 1 no solution 2 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 D r R0 B= 0.15 inertial jet viscous-inertial jet no solution 1 no solution 2 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 D r R0 B= 0.2617 inertial jet no solution 1 no solution 2 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 D r R0 B= 0.4 inertial jet no solution

Fig. 3 Parameter regions in the R0, Dr-plane for different B.

For B close to zero (see the plot for B= 0.01) the region of the inertial jet

occupies almost the whole plane except a narrow area near Dr= 0. The

viscous-inertial jet region forms long and narrow strip along the line Dr= 0. The

nonex-istence regions ”no solution 1” and ”no solution 2” are small areas concentrated near the points Dr= 0, R0= 0, and Dr = 0, R0= 1, respectively. With increasing B (see the plot for B= 0.15) the border of the inertial jet region rises, the

viscous-inertial jet region becomes higher and shorter, and the nonexistence regions grow. For B= 0.2617 the viscous-inertial jet regions is just disappeared (shrinked to one

curve). For larger B> 0.2617 (see the plot for B = 0.4) the viscous-inertial jet is

absent, the inertial jet region becomes smaller and the nonexistence region ”no so-lution” expands.

3 Conclusions

The steady jet model in rotary spinning configuration is studied. The jet is described by a system of ODE’s on an unknown domain. The scaled system describing the jet is characterized by three dimensionless parameters. For the jet in rotary spinning we distinguish three situations: the inertial jet, the viscous-inertial jet, and a steady jet solution does not exist. There are two causes of the nonexistence of the jet between the rotor and the coagulator 1) the jet wraps around the rotor, and 2) the flow velocity at the coagulator can not be reached, indicating unsteady jet.

(8)

6 A. Hlod, A.A.F. van de Ven, and M.A. Peletier

Acknowledgements The authors would like to acknowledge Teijin Aramid, a part of the Teijin

group of companies.

References

1. Decent, S.P., King, A.C., Wallwork, I.M.: Free jets spun from a prilling tower. Journal of Engineering Mathematics 42(3), 265 – 282 (2002)

2. den Decker, P., Knoester, H., Meerman, H., Dekker K. van Horssen, W., Vuik, C., Wesseling, P., Prokert, G., van ’t Hof, B., van Beckum, F.: The rotor spinning process for fibre production. In: Proceedings of the 48th European Study Group Mathematics with Industry (Delft, 15-19 March 2004), pp. 35–48 (2004)

3. Hlod, A.: Curved jets of viscous fluid: Interactions with a moving wall. Ph.D. thesis, Eind-hoven University of Technology (2009)

4. Kolk, E.: Mathematical models for a rotor spinning process. Interim report, TU Delft (2005) 5. Kolk, E.: Mathematical models for a rotor spinning process. Final report, TU Delft (2005) 6. Marheineke, N., Wegener, R.: Asymptotic model for the dynamics of curved viscous fibres

with surface tension. Journal of Fluid Mechanics 622(-1), 345–369 (2009)

7. Panda, S.: The dynamics of viscous fibers. Ph.D. thesis, Technische Universit¨at Kaiserslautern (2006)

8. Panda, S., Marheineke, N., Wegener, R.: Systematic derivation of an asymptotic model for the dynamics of curved viscous fibers. Mathematical Methods in the Applied Sciences 31(10), 1153–1173 (2008)

9. Parau, E.I., Decent, S., Simmons, M., Wong, D.C.Y., King, A.C.: Nonlinear viscous liquid jets from a rotating orifice. Journal of Engineering Mathematics 57(2), 159–179 (2007)

10. Parau, E.I., Decent, S.P., King, A.C., Simmons, M.J.H., Wong, D.C.: Nonlinear travelling waves on a spiralling liquid jet. Wave Motion 43(7), 599–618 (2006)

11. Partridge, L., Wong, D.C.Y., Simmons, M.J.H., Parau, E.I., Decent, S.P.: Experimental and theoretical description of the break-up of curved liquid jets in the prilling process. Chemical Engineering Research and Design 83(A11), 1267–1275 (2005)

12. Uddin, J., Decent, S.P., Simmons, M.J.: The instability of shear thinning and shear thickening spiralling liquid jets: linear theory. Transactions of the ASME Journal of Fluids Engineering

128(5), 968–975 (2006)

13. Wallwork, I.M., Decent, S.P., King, A.C., Schulkes, R.M.S.M.: The trajectory and stability of a spiralling liquid jet. part 1. inviscid theory. Journal of Fluid Mechanics 459, 43–66 (2002) 14. Wong, D.C.Y., Simmons, M.J.H., Decent, S.P., Parau, E.I., King, A.C.: Break-up dynamics

and drop size distributions created from spiralling liquid jets. International Journal of Multi-phase Flow 30(5), 499–520 (2004)

(9)

PREVIOUS PUBLICATIONS IN THIS SERIES:

Number Author(s)

Title

Month

10-61

10-62

10-63

10-64

10-65

V. Savcenco

B. Haut

E.J.W. ter Maten

R.M.M. Mattheij

D.J. Bekers

S.J.L. van Eijndhoven

L.M.J. Florack

F. Schilder

R. Christensen

J. Sylvest

K. Kumar

M. Rudnaya

V. Savcenco

B. Sobiegraj

C. Please

R. Style

A. Hlod

A.A.F. van de Ven

M.A. Peletier

Time domain simulation

of power systems with

different time scales

Spectral analysis of

integral-differential

operators applied in linear

antenna modeling

Scale space

representations locally

adapted to the geometry of

base and target manifold

Design parameters for a

siphon system

A model of rotary spinning

process

Oct. ‘10

Oct. ‘10

Oct. ‘10

Oct. ‘10

Oct. ‘10

Ontwerp: de Tantes, Tobias Baanders, CWI

Referenties

GERELATEERDE DOCUMENTEN

• The final published version features the final layout of the paper including the volume, issue and page numbers.. Link

The major differentiators in the place order on supplier process are cost (affecting finished product costs), quality and reliability (availability and security of supply

Its focus on rituals that affect forgiveness between God and human beings, as well as between human agents is important, since the Gospel of Matthew recounts both the

For the outcomes part, it is considered how eventual data from components can help the organization in the need/wish of further preventive maintenance philosophy working and how

We hebben besloten om alleen die geneesmiddelen te beoordelen die een financieel risico vormen voor de basisverzekering of indien de kwaliteit van zorg een risico vormt?. We zijn

Een megabedrijf is dan de extreme ver - beelding van de landbouw waarvan men dacht dat die niet bestond.’ &gt;&gt; Hoe moeten onderzoekers en andere betrokkenen het een volgende

Duizenden jaren geleden moeten de eerste hoenders in Zuid-Oost Azië, op zoek naar voedselresten, de sprong hebben gewaagd van de rand van het woud naar menselijke nederzettingen.

“Dat zijn eendimensionale gewas- en gelddenkers.” Na zijn studie aan de heao heeft hij in 2007 en 2008 eerst een jaartje als backpacker de wereld rondgereisd, door Australië