Skip to main content

Fluid-fiber-interactions in rotational spinning process of glass wool production

Abstract

The optimal design of rotational production processes for glass wool manufacturing poses severe computational challenges to mathematicians, natural scientists and engineers. In this paper we focus exclusively on the spinning regime where thousands of viscous thermal glass jets are formed by fast air streams. Homogeneity and slenderness of the spun fibers are the quality features of the final fabric. Their prediction requires the computation of the fluid-fiber-interactions which involves the solving of a complex three-dimensional multiphase problem with appropriate interface conditions. But this is practically impossible due to the needed high resolution and adaptive grid refinement. Therefore, we propose an asymptotic coupling concept. Treating the glass jets as viscous thermal Cosserat rods, we tackle the multiscale problem by help of momentum (drag) and heat exchange models that are derived on basis of slender-body theory and homogenization. A weak iterative coupling algorithm that is based on the combination of commercial software and self-implemented code for flow and rod solvers, respectively, makes then the simulation of the industrial process possible. For the boundary value problem of the rod we particularly suggest an adapted collocation-continuation method. Consequently, this work establishes a promising basis for future optimization strategies.

1 Introduction

Glass wool manufacturing requires a rigorous understanding of the rotational spinning of viscous thermal jets exposed to aerodynamic forces. Rotational spinning processes consist in general of two regimes: melting and spinning. The plant of our industrial partner, Woltz GmbH in Wertheim, is illustrated in Figures 1 and 2. Glass is heated upto temperatures of 1,050°C in a stove from which the melt is led to a centrifugal disk. The walls of the disk are perforated by 35 rows over height with 770 equidistantly placed small holes per row. Emerging from the rotating disk via continuous extrusion, the liquid jets grow and move due to viscosity, surface tension, gravity and aerodynamic forces. There are in particular two different air flows that interact with the arising glass fiber curtain: a downwards-directed hot burner flow of 1,500°C that keeps the jets near the nozzles warm and thus extremely viscous and shapeable as well as a highly turbulent cross-stream of 30°C that stretches and finally cools them down such that the glass fibers become hardened. Laying down onto a conveyor belt they yield the basic fabric for the final glass wool product. For the quality assessment of the fabrics the properties of the single spun fibers, that is, homogeneity and slenderness, play an important role. A long-term objective in industry is the optimal design of the manufacturing process with respect to desired product specifications and low production costs. Therefore, it is necessary to model, simulate and control the whole process.

Figure 1
figure 1

Rotational spinning process of the company Woltz GmbH, sketch of set-up. Several glass jets forming part of the row-wise arising fiber curtains are shown in the left part of the disc, they are plotted as black curves. The color map visualizes the axial velocity of the air flow. For temperature details see Figure 2.

Figure 2
figure 2

Rotational spinning process of the company Woltz GmbH, sketch of set-up. Several glass jets forming part of the row-wise arising fiber curtains are shown in the left part of the disc, they are plotted as black curves. The color map visualizes the temperature of air flow. For velocity details see Figure 1.

Up to now, the numerical simulation of the whole manufacturing process is impossible because of its enormous complexity. In fact, we do not long for an uniform numerical treatment of the whole process, but have the idea to derive adequate models and methods for the separate regimes and couple them appropriately, for a similar strategy for technical textiles manufacturing see [1]. In this content, the melting regime dealing with the creeping highly viscous melt flow from the stove to the holes of the centrifugal disk might be certainly handled by standard models and methods from the field of fluid dynamics. It yields the information about the melt velocity and temperature distribution at the nozzles which is of main importance for the ongoing spinning regime. However, be aware that for their determination not only the melt behavior in the centrifugal disk but also the effect of the burner flow, that is, aerodynamic heating and heat distortion of disk walls and nozzles, have to be taken into account. In this paper we assume the conditions at the nozzles to be given and focus exclusively on the spinning regime which is the challenging core of the problem. For an overview of the specific temperature, velocity and length values we refer to Table 1. In the spinning regime the liquid viscous glass jets are formed, in particular they are stretched by a factor 10,000. Their geometry is characterized by a typical slenderness ratio δ=d/l 10 4 of jet diameter d and length l. The resulting fiber properties (characteristics) depend essentially on the jets behavior in the surrounding air flow. To predict them, the interactions, that is, momentum and energy exchange, of air flow and fiber curtain consisting of MN single jets (M=35N=770) have to be considered. Their computation requires in principle a coupling of fiber jets and flow with appropriate interface conditions. However, the needed high resolution and adaptive grid refinement make the direct numerical simulation of the three-dimensional multiphase problem for ten thousands of slender glass jets and fast air streams not only extremely costly and complex, but also practically impossible. Therefore, we tackle the multiscale problem by help of drag models that are derived on basis of slender-body theory and homogenization, and a weak iterative coupling algorithm.

Table 1 Typical temperature, velocity and length values in the considered rotational spinning process, cf

The dynamics of curved viscous inertial jets is of interest in many industrial applications (apart from glass wool manufacturing), for example, in nonwoven production [1, 2], pellet manufacturing [3, 4] or jet ink design, and has been subject of numerous theoretical, numerical and experimental investigations, see [5] and references within. In the terminology of [6], there are two classes of asymptotic one-dimensional models for a jet, that is, string and rod models. Whereas the string models consist of balance equations for mass and linear momentum, the more complex rod models contain also an angular momentum balance, [7, 8]. A string model for the jet dynamics was derived in a slender-body asymptotics from the three-dimensional free boundary value problem given by the incompressible Navier-Stokes equations in [5]. Accounting for inner viscous transport, surface tension and placing no restrictions on either the motion or the shape of the jet’s center-line, it generalizes the previously developed string models for straight [911] and curved [1214] center-lines. However, already in the stationary case the applicability of the string model turns out to be restricted to certain parameter ranges [15, 16] because of a non-removable singularity that comes from the deduced boundary conditions. These limitations can be overcome by a modification of the boundary conditions, that is, the release of the condition for the jet tangent at the nozzle in favor of an appropriate interface condition, [1719]. This involves two string models that exclusively differ in the closure conditions. For gravitational spinning scenarios they cover the whole parameter range, but in the presence of rotations there exist small parameter regimes where none of them works. A rod model that allows for stretching, bending and twisting was proposed and analyzed in [20, 21] for the coiling of a viscous jet falling on a rigid substrate. Based on these studies and embedded in the special Cosserat theory a modified incompressible isothermal rod model for rotational spinning was developed and investigated in [16, 19]. It allows for simulations in the whole (Re, Rb, Fr)-range and shows its superiority to the string models. These observations correspond to studies on a fluid-mechanical ‘sewing machine’, [22, 23]. By containing the slenderness parameter δ explicitely in the angular momentum balance, the rod model is no asymptotic model of zeroth order. Since its solutions converge to the respective string solutions in the slenderness limit δ0, it can be considered as δ-regularized model, [19]. In this paper we extend the rod model by incorporating the practically relevant temperature dependencies and aerodynamic forces. Thereby, we use the air drag model F of [24] that combines Oseen and Stokes theory [2527], Taylor heuristic [28] and numerical simulations. Being validated with experimental data [2932], it is applicable for all air flow regimes and incident flow directions. Transferring this strategy, we model a similar aerodynamic heat source for the jet that is based on the Nusselt number Nu [33]. Our coupling between glass jets and air flow follows then the principle that action equals reaction. By inserting the corresponding homogenized source terms induced by F and Nu in the balance equations of the air flow, we make the proper momentum and energy exchange within this slender-body framework possible.

The paper is structured as follows. We start with the general coupling concept for slender bodies and fluid flows. Therefore, we introduce the viscous thermal Cosserat rod system and the compressible Navier-Stokes equations for glass jets and air flow, respectively, and present the models for the momentum and energy exchange: drag F and Nusselt function Nu. The special set-up of the industrial rotational spinning process allows for the simplification of the model framework, that is, transition to stationarity and assumption of rotational invariance as we discuss in detail. It follows the section about the numerical treatment. To realize the fiber-flow interactions we use a weak iterative coupling algorithm, which is adequate for the problem and has the advantage that we can combine commercial software and self-implemented code. Special attention is paid to the collocation and continuation method for solving the boundary value problem of the rod. Convergence of the coupling algorithm and simulation results are shown for a specific spinning adjustment. This illustrates the applicability of our coupling framework as one of the basic tools for the optimal design of the whole manufacturing process. Finally, we conclude with some remarks to the process.

2 General coupling concept for slender bodies and fluid flows

We are interested in the spinning of ten thousands of slender glass jets by fast air streams, MN=26,950. The glass jets form a kind of curtain that interact and crucially affect the surrounding air. The determination of the fluid-fiber-interactions requires in principle the simulation of the three-dimensional multiphase problem with appropriate interface conditions. However, regarding the complexity and enormous computational effort, this is practically impossible. Therefore, we propose a coupling concept for slender bodies and fluid flows that is based on drag force and heat exchange models. In this section we first present the two-way coupling of a single viscous thermal Cosserat rod and the compressible Navier-Stokes equations and then generalize the concept to many rods. Thereby, we choose an invariant formulation in the three-dimensional Euclidian space E 3 .

Note that we mark all quantities associated to the air flow by the subscript throughout the paper. Moreover, to facilitate the readability of the coupling concept, we introduce the abbreviations Ψ and Ψ that represent all quantities of the glass jets and the air flow, respectively.

2.1 Models for glass jets and air flows

2.1.1 Cosserat rod

A glass jet is a slender body, that is, a rod in the context of three-dimensional continuum mechanics. Because of its slender geometry, its dynamics might be reduced to a one-dimensional description by averaging the underlying balance laws over its cross-sections. This procedure is based on the assumption that the displacement field in each cross-section can be expressed in terms of a finite number of vector- and tensor-valued quantities. In the special Cosserat rod theory, there are only two constitutive elements: a curve specifying the position r:Q E 3 and an orthonormal director triad { d 1 , d 2 , d 3 }:Q E 3 characterizing the orientation of the cross-sections, where Q={(s,t) R 2 |sI(t)=[0,l(t)],t>0} with arclength parameter s and time t. For a schematic sketch of a Cosserat rod see Figure 3, for more details on the Cosserat theory we refer to [6]. In the following we use an incompressible viscous Cosserat rod model that was derived on basis of the work [20, 34] on viscous rope coiling and investigated for isothermal curved inertial jets in rotational spinning processes [16, 19]. We extend the model by incorporating temperature effects and aerodynamic forces. The rod system describes the variables of jet curve r, orthonormal triad { d 1 , d 2 , d 3 }, generalized curvature κ, convective speed u, cross-section A, linear velocity v, angular velocity ω, temperature T and normal contact forces n d α α=1,2. It consists of four kinematic and four dynamic equations, that is, balance laws for mass (cross-section), linear and angular momentum and temperature,

t r = v u d 3 , t d i = ( ω u κ ) × d i , s r = d 3 , s d i = κ × d i , t A + s ( u A ) = 0 , ρ ( t ( A v ) + s ( u A v ) ) = s n + ρ A g e g + f air , ρ ( t ( J ω ) + s ( u J ω ) ) = s m + d 3 × n , ρ c p ( t ( A T ) + s ( u A T ) ) = q rad + q air
(1)

supplemented with an incompressible geometrical model of circular cross-sections with diameter d

J = I ( d 1 d 1 + d 2 d 2 + 2 d 3 d 3 ) , I = π 64 d 4 , A = π 4 d 2

as well as viscous material laws for the tangential contact force n d 3 and contact couple m

n d 3 = 3 μ A s u m = 3 μ I ( d 1 d 1 + d 2 d 2 + 2 3 d 3 d 3 ) s ω .

Rod density ρ and heat capacity c p are assumed to be constant. The temperature-dependent dynamic viscosity μ is modeled according to the Vogel-Fulcher-Tamman relation, that is, μ(T)= 10 p 1 + p 2 / ( T p 3 ) Pa s where we use the parameters p 1 =2.56 p 2 =4,289.18 K and p 3 =(150.74+273.15) K, [33]. The external loads rise from gravity ρAg e g with gravitational acceleration g and aerodynamic forces f air . In the temperature equation we neglect inner friction and heat conduction and focus exclusively on radiation q rad and aerodynamic heat sources q air . The radiation effect depends on the geometry of the plant and is incorporated in the system by help of the simple model

q rad = ε R σ π d ( T ref 4 T 4 )

with emissivity ε R , Stefan-Boltzmann constant σ and reference temperature T ref . Appropriate initial and boundary conditions close the rod system.

Figure 3
figure 3

Special Cosserat rod with Kirchhoff constraint s r= d 3 .

2.1.2 Navier-Stokes equations

A compressible air flow in the space-time domain Ω t ={(x,t)|xΩ E 3 ,t>0} is described by density ρ , velocity v , temperature T . Its model equations consist of the balance laws for mass, momentum and energy,

t ρ + ( v ρ ) = 0 , t ( ρ v ) + ( v ρ v ) = S T + ρ g e g + f jets ,
(2)
t ( ρ e ) + ( v ρ e ) = S : v q + q jets

supplemented with the Newtonian stress tensor S , the Fourier law for heat conduction q

S = p I + μ ( v + v T ) + λ v I , q = k T ,

as well as thermal and caloric equations of state of a ideal gas

p = ρ R T , e = 0 T c p ( T ) d T p ρ

with pressure p and inner energy e . The specific gas constant for air is denoted by R . The temperature-dependent viscosities μ λ , heat capacity c p and heat conductivity k can be modeled by standard polynomial laws, see, for example, [33, 35]. External loads rise from gravity ρ g e g and forces due to the immersed fiber jets f jets . These fiber jets also imply a heat source q jets in the energy equation. Appropriate initial and boundary conditions close the system.

2.2 Models for momentum and energy exchange

The coupling of the Cosserat rod and the Navier-Stokes equations is performed by help of drag forces and heat sources. Taking into account the conservation of momentum and energy, f air and f jets as well as q air and q jets satisfy the principle that action equals reaction and obey common underlying relations. Hence, we can handle the delicate fluid-fiber-interactions by help of two surrogate models, so-called exchange functions, that is, a dimensionless drag force F (inducing f air f jets ) and Nusselt number Nu (inducing q air q jets ). For a flow around a slender long cylinder with circular cross-sections there exist plenty of theoretical, numerical and experimental investigations to these relations in literature, for an overview see [24] as well as, for example, [29, 30, 33, 36] and references within. We use this knowledge locally and globalize the models by superposition to apply them to our curved moving Cosserat rod. This strategy follows a Global-from-Local concept [37] that turned out to be very satisfying in the derivation and validation of a stochastic drag force in a one-way coupling of fibers in turbulent flows [24].

2.2.1 Drag forces - f air vs f jets

Let Ψ and Ψ represent all glass jet and air flow quantities, respectively. Thereby, Ψ is the spatially averaged solution of (2). This delocation is necessary to avoid singularities in the two-way coupling. Then, the drag forces are given by

f air ( s , t ) = F ( Ψ ( s , t ) , Ψ ( r ( s , t ) , t ) ) , f jets ( x , t ) = I ( t ) δ ( x r ( s , t ) ) F ( Ψ ( s , t ) , Ψ ( x , t ) ) d s , F ( Ψ , Ψ ) = μ 2 d ρ F ( d 3 , d ρ μ ( v v ) ) ,

where δ is the Dirac distribution. By construction, they fulfill the principle that action equals reaction and hence the momentum is conserved, that is,

I V ( t ) f air ( s , t ) d s = V f jets ( x , t ) d x

for an arbitrary domain V and I V (t)={sI(t)|r(s,t)V}. The (line) force F acting on a slender body is caused by friction and inertia. It depends on material and geometrical properties as well as on the specific inflow situation. The number of dependencies can be reduced to two by help of non-dimensionalizing which yields the dimensionless drag force F in dependence on the jet orientation (tangent) and the dimensionless relative velocity between air flow and glass jet. Due to the rotational invariance of the force, the function

F : S 2 × E 3 E 3

can be associated with its component tuple F for every representation in an orthonormal basis, that is,

F : S R 3 2 × R 3 R 3 , F = ( F 1 , F 2 , F 3 ) with i = 1 3 F i ( τ , w ) e i = F ( i = 1 3 τ i e i , i = 1 3 w i e i )

for every orthonormal basis { e i }.

For F we use the drag model [24] that was developed on top of Oseen and Stokes theory [2527], Taylor heuristic [28] and numerical simulations and validated with measurements [2932]. It shows to be applicable for all air flow regimes and incident flow directions. Let {n,b,τ} be the orthonormal basis induced by the specific inflow situation (τ,w) with orientation τ and velocity w, assuming wτ

n = w w τ τ w n , b = τ × n , w τ = w τ , w n = w 2 w τ 2 .

Then, the force is given by

F ( τ , w ) = F n ( w n ) n + F τ ( w n , w τ ) τ , F n ( w n ) = w n 2 c n ( w n ) = w n r n ( w n ) ,
(3)
F τ ( w n , w τ ) = w τ w n c τ ( w n ) = w τ r τ ( w n )

according to the Independence Principle [38]. The differentiable normal and tangential drag functions c n c τ are

c n ( w n ) = { 4 π S w n [ 1 w n 2 S 2 S / 2 + 5 / 16 32 S ] , w n < w 1 , exp ( j = 0 3 p n , j ln j w n ) , w 1 w n w 2 , 2 w n + 0.5 , w 2 < w n , c τ ( w n ) = { 4 π ( 2 S 1 ) w n [ 1 w n 2 2 S 2 2 S + 1 16 ( 2 S 1 ) ] , w n < w 1 , exp ( j = 0 3 p τ , j ln j w n ) , w 1 w n w 2 , γ w n , w 2 < w n ,

with S( w n )=2.0022ln w n , transition points w 1 =0.1 w 2 =100, amplitude γ=2. The regularity involves the parameters p n , 0 =1.6911 p n , 1 =6.7222 10 1 p n , 2 =3.3287 10 2 p n , 3 =3.5015 10 3 and p τ , 0 =1.1552 p τ , 1 =6.8479 10 1 p τ , 2 =1.4884 10 2 p τ , 3 =7.4966 10 4 . To be also applicable in the special case of a transversal incident flow wτ and to ensure a realistic smooth force F, the drag is modified for w n 0. A regularization based on the slenderness parameter δ matches the associated resistance functions r n r τ (3) to Stokes resistance coefficients of higher order for w n 1, for details see [24].

2.2.2 Heat sources - q air vs q jets

Analogously to the drag forces, the heat sources are given by

q air ( s , t ) = Q ( Ψ ( s , t ) , Ψ ( r ( s , t ) , t ) ) , q jets ( x , t ) = I ( t ) δ ( x r ( s , t ) ) Q ( Ψ ( s , t ) , Ψ ( x , t ) ) d s , Q ( Ψ , Ψ ) = 2 k ( T T ) Nu ( v v v v d 3 , π 2 d ρ μ v v , μ c p k ) .

The (line) heat source Q acting on a slender body also depends on several material and geometrical properties as well as on the specific inflow situation. The number of dependencies can be reduced to three by help of non-dimensionalizing which yields the dimensionless Nusselt number Nu in dependence of the cosine of the angle of attack, Reynolds and Prandtl numbers. The Reynolds number corresponds to the relative velocity between air flow and glass jet, the typical length is the half jet circumference.

For Nu we use a heuristic model. It originates in the studies of a perpendicular flow around a cylinder [33] and is modified for different inflow directions (angles of attack) with regard to experimental data. A regularization ensures the smooth limit for a transversal incident flow in analogon to the drag model for F in (3). We apply

(4)
Nu t u ( Re , Pr ) = 0.037 Re 0.9 Pr Re 0.1 + 2.443 ( Pr 2 / 3 1 ) , h ( c , Re ) = { c Re / δ h , Re < δ h , c , Re δ h .

2.3 Generalization to many rods

In case of k slender bodies in the air flow, we have Ψ i , i=1,,k, representing the quantities of each Cosserat rod, here k=MN. Assuming no contact between neighboring fiber jets, every single jet can be described by the stated rod system (1). Their multiple effect on the air flow is reflected in f jets and q jets . The source terms in the momentum and energy equations of the air flow (2) become

f jets ( x , t ) = i = 1 k I i ( t ) δ ( x r i ( s , t ) ) F ( Ψ i ( s , t ) , Ψ ( x , t ) ) d s , q jets ( x , t ) = i = 1 k I i ( t ) δ ( x r i ( s , t ) ) Q ( Ψ i ( s , t ) , Ψ ( x , t ) ) d s .

3 Models for special set-up of rotational spinning process

In the rotational spinning process under consideration the centrifugal disk is perforated by M rows of N equidistantly placed holes each (M=35, N=770). The spinning conditions (hole size, velocities, temperatures) are thereby identical for each row, see Figures 1 and 2. The special set-up allows for the simplification of the general model framework. We introduce the rotating outer orthonormal basis { a 1 (t), a 2 (t), a 3 (t)} satisfying t a i =Ω× a i , i=1,2,3, where Ω is the angular frequency of the centrifugal disk. In particular, Ω=Ω a 1 and e g = a 1 (gravity direction) hold. Then, glass jets and air flow become stationary, presupposing that we consider spun fiber jets of certain length. In particular, we assume the stresses to be vanished at this length. Moreover, the glass jets emerging from the rotating device form dense curtains for every spinning row. As a result of homogenization, we can treat the air flow as rotationally invariant and each curtain can be represented by one jet. This yields an enormous complexity reduction of the problem. The homogenization together with the slender-body theory makes the numerical simulation possible.

3.1 Transition to stationarity

3.1.1 Representative spun jet of certain length

For the viscous Cosserat rods (1), the mass flux Q is constant in the stationarity, that is, uA=Q/ρ=const. We deal with Ω-adapted linear and angular velocities, v Ω =vΩ×r and ω Ω =ωΩ, which fulfill the explicit stationarity relations

v Ω = u d 3 , ω Ω = u κ

resulting from the first two equations of (1). Moreover, fictitious Coriolis and centrifugal forces and associated couples enter the linear and angular momentum equations. Using the material laws we can formulate the stationary rod model in terms of a boundary value problem of first order differential equations. Thereby, we present it in the director basis { d 1 , d 2 , d 3 } for convenience (see (5) and compare to [19] except for the temperature equation). Note that to an arbitrary vector field z= i = 1 3 z ˘ i a i = i = 1 3 z i d i E 3 , we indicate the component tuples corresponding to the rotating outer basis and the director basis by z ˘ =( z ˘ 1 , z ˘ 2 , z ˘ 3 ) R 3 and z=( z 1 , z 2 , z 3 ) R 3 , respectively. The director basis can be transformed into the rotating outer basis by the tensor-valued rotation R, that is, R= a i d i = R i j a i a j E 3 E 3 with associated orthogonal matrix R=( R i j )=( d i a j )SO(3). Its transpose and inverse matrix is denoted by R T . For the components, z=R z ˘ holds. The cross-product z×R is defined as mapping (z×R): R 3 R 3 yz×(Ry). Moreover, canonical basis vectors in R 3 are denoted by e i i=1,2,3, for example, e 1 =(1,0,0). Then, the stationary Cosserat rod model stated in the director basis for a spun glass jet reads

s r ˘ = R T e 3 , s R = κ × R , s κ = ρ 3 Q κ n 3 μ + 4 π ρ 2 3 Q 2 u μ P 3 / 2 m , s u = ρ 3 Q u n 3 μ , s n = κ × n + Q u κ × e 3 + ρ 3 u n 3 μ e 3 + 2 Q Ω ( R e 1 ) × e 3 + Q Ω 2 1 u R ( e 1 × ( e 1 × r ˘ ) ) + Q g 1 u R e 1 R f ˘ air , s m = κ × m + n × e 3 + ρ 3 u μ P 3 m Q 12 π n 3 μ P 2 κ Q Ω 12 π n 3 μ u P 2 ( R e 1 ) Q 2 Ω 4 π ρ 1 u P 2 ( κ × R e 1 ) Q 2 4 π ρ 1 u 2 P 2 ( u κ + Ω R e 1 ) × ( u κ + Ω R e 1 ) , s T = 1 c p Q ( q rad + q air )
(5)

with q rad =2 π ε R σ Q / ρ ( T ref 4 T 4 )/ u and diagonal matrix P k =diag(1,1,k)kR. For a spun jet emerging from the centrifugal disk at s=0 with stress-free end at s=L, the equations are supplemented with

(cf. Table 1). Considering the jet as representative of one spinning row, we choose the nozzle position to be (H,R,0) with respective height H R is here the disk radius. The initialization R(0) prescribes the jet direction at the nozzle as ( d 1 , d 2 , d 3 )(0)=( a 1 , a 3 , a 2 ).

Remark 1 The rotationsRSO(3)can be parameterized for example in Euler angles or unit quaternions[39]. The last variant offers a very elegant way of rewriting the second equation of (5). Define

R ( q ) = ( q 1 2 q 2 2 q 3 2 + q 0 2 2 ( q 1 q 2 q 0 q 3 ) 2 ( q 1 q 3 + q 0 q 2 ) 2 ( q 1 q 2 + q 0 q 3 ) q 1 2 + q 2 2 q 3 2 + q 0 2 2 ( q 2 q 3 q 0 q 1 ) 2 ( q 1 q 3 q 0 q 2 ) 2 ( q 2 q 3 + q 0 q 1 ) q 1 2 q 2 2 + q 3 2 + q 0 2 ) ,
q=( q 0 , q 1 , q 2 , q 3 )

withq=1then we have s q=A(κ)qwith skew-symmetric matrix

A ( κ ) = 1 2 ( 0 κ 1 κ 2 κ 3 κ 1 0 κ 3 κ 2 κ 2 κ 3 0 κ 1 κ 3 κ 2 κ 1 0 ) .

3.1.2 Rotationally invariant air flow

Due to the spinning set-up the jets emerging from the rotating device form row-wise dense curtains. As a consequence of a row-wise homogenization, the air flow (2) can be treated as stationary not only in the rotating outer basis { a 1 (t), a 2 (t), a 3 (t)}, but also in a fixed outer one. Because of the symmetry with respect to the rotation axis, it is convenient to introduce cylindrical coordinates (x,r,ϕ)R× R + ×[0,2π) for the space and to attach a cylindrical basis { e x , e r , e ϕ } with e x = a 1 to each space point. The components to an arbitrary vector field z E 3 are indicated by z ˆ =( z x , z r , z ϕ ) R 3 . Then, taking advantage of the rotational invariance, the stationary Navier-Stokes equations in (x,r) simplify to

x ( ρ v x ) + 1 r r ( r ρ v r ) = 0 , x ( ρ v x 2 ) + 1 r r ( r ρ v r v x ) = x p + x ( 2 μ x v x + λ v ˆ ) + 1 r r ( r μ ( x v r + r v x ) ) ρ g + ( f x ) jets , x ( ρ v x v r ) + 1 r r ( r ρ v r 2 ) 1 r ρ v ϕ 2 = r p + x ( μ ( x v r + r v x ) ) + 2 r r ( r μ r v r ) + r ( λ v ˆ ) 2 r 2 μ v r + ( f r ) jets , x ( ρ v x v ϕ ) + 1 r r ( r ρ v r v ϕ ) + 1 r ρ v r v ϕ = x ( μ x v ϕ ) + 1 r 2 r ( r 3 μ r ( 1 r v ϕ ) ) + ( f ϕ ) jets ,
(6)

with v ˆ = x v x +( r (r v r ))/r and equipped with appropriate inflow, outflow and wall boundary conditions, cf. Figures 1 and 2.

3.2 Exchange functions

To perform the coupling between (5) and (6), we have to compute the exchange functions in the appropriate coordinates. These calculations are simplified by the rotational invariance of the problem. As introduced, we use the subscripts ˘ and ˆ to indicate the component tuples corresponding to the rotating outer basis { a 1 (t), a 2 (t), a 3 (t)} and the cylindrical basis { e x , e r , e ϕ }, respectively. Essentially for the coupling are the jet tangent and the relative velocity between air flow and glass jet, they are

τ ˘ = R T e 3 , τ ˆ = ( τ ˘ 1 , r ˘ 2 τ ˘ 2 + r ˘ 3 τ ˘ 3 r ˘ 2 2 + r ˘ 3 2 , r ˘ 2 τ ˘ 3 r ˘ 3 τ ˘ 2 r ˘ 2 2 + r ˘ 3 2 )

and

v ˘ rel = ( v x , r ˘ 2 v r r ˘ 3 ( v ϕ Ω r ) r ˘ 2 2 + r ˘ 3 2 , r ˘ 3 v r + r ˘ 2 ( v ϕ Ω r ) r ˘ 2 2 + r ˘ 3 2 ) u τ ˘ , v ˆ rel = ( v x , v r , v ϕ Ω r ) u τ ˆ .

Then, the drag forces are

f ˘ air ( s ) = F ˘ ( Ψ ( s ) , Ψ ( r ˘ 1 ( s ) , r ˘ 2 2 ( s ) + r ˘ 3 2 ( s ) ) ) , f ˆ jets ( x , r ) = N 2 π I 1 r δ ( x r ˘ 1 ( s ) ) δ ( r r ˘ 2 2 ( s ) + r ˘ 3 2 ( s ) ) F ˆ ( Ψ ( s ) , Ψ ( x , r ) ) d s , F ˘ ( Ψ , Ψ ) = 2 Q π ρ μ 2 ρ 1 u F ( τ ˘ , 2 Q π ρ ρ μ 1 u v ˘ rel ) , F ˆ ( Ψ , Ψ ) = 2 Q π ρ μ 2 ρ 1 u F ( τ ˆ , 2 Q π ρ ρ μ 1 u v ˆ rel )

and the heat sources

q air ( s ) = Q ( Ψ ( s ) , Ψ ( r ˘ 1 ( s ) , r ˘ 2 2 ( s ) + r ˘ 3 2 ( s ) ) ) , q jets ( x , r ) = N 2 π I 1 r δ ( x r ˘ 1 ( s ) ) δ ( r r ˘ 2 2 ( s ) + r ˘ 3 2 ( s ) ) Q ( Ψ ( s ) , Ψ ( x , r ) ) d s , Q ( Ψ , Ψ ) = 2 k ( T T ) Nu ( v ˘ rel v ˘ rel τ ˘ , π Q ρ ρ μ v ˘ rel u , μ c p k ) .

Here, f ˆ jets and q jets represent the homogenized effect of the N glass jets emerging from the equidistantly placed holes in an arbitrary spinning row. Correspondingly, system (5) with f ˘ air and q air describes one representative glass jet for this row. To simulate the full problem with all MN glass jets in the air, jet representatives Ψ i , i=1,,M for all M spinning rows with the respective boundary and air flow conditions have to be determined. Their common effect on the air flow is

4 Numerical treatment

The numerical simulation of the glass jets dynamic in the air flow is performed by an algorithm that weakly couples glass jet calculation and air flow computation via iterations. This procedure is adequate for the problem and has the advantage that we can combine commercial software and self-implemented code. We use FLUENT, a commercial finite volume-based software by ANSYS, that contains the broad physical modeling capabilities needed to describe air flow, turbulence and heat transfer for the industrial glass wool manufacturing process. In particular, a pressure-based solver is applied in the computation of (6). To restrict the computational effort in grid refinement needed for the resolution of the turbulent air streams we consider alternatively a stochastic k-ω turbulence model. (For details on the commercial software FLUENT, its models and solvers we refer to http://www.fluent.com.) Note that the modification of the model equations has no effect on our coupling framework, where the exchange functions are incorporated by UDFs (user defined functions). For the boundary value problem of the stationary Cosserat rod (5), systems of nonlinear equations are set up via a Runge-Kutta collocation method and solved by a Newton method in MATLAB 7.4. The convergence of the Newton method depends thereby crucially on the initial guess. To improve the computational performance we adapt the initial guess iteratively by solving a sequence of boundary value problems with slightly changed parameters. The developed continuation method is presented in the following. Moreover, to get a balanced numerics we use the dimensionless rod system that is scaled with the respective conditions at the nozzle. The M glass jet representative are computed in parallel. The exchange of flow and fiber data between the solvers is based on interpolation and averaging, as we explain in the weak iterative coupling algorithm.

4.1 Collocation-continuation method for dimensionless rod boundary value problem

The computing of the glass jets is based on a dimensionless rod system. For this purpose, we scale the dimensional equations (5) with the spinning conditions of the respective row. Apart from the air flow data, (5) contains thirteen physical parameters, that is, jet density ρ, heat capacity c p , emissivity ε R , typical length L, velocity U and temperature θ at the spinning hole as well as hole diameter D and height H, centrifugal disk radius R, rotational frequency Ω, reference temperature for radiation T ref and gravitational acceleration g. The typical jet viscosity is chosen to be μ 0 =μ(θ). These induce various dimensionless numbers characterizing the fiber spinning, that is, Reynolds number Re as ratio between inertia and viscosity, Rossby number Rb as ratio between inertia and rotation, Froude number Fr as ratio between inertia and gravity and Ra as ratio between radiation and heat advection as well as h and ϵ as length ratios between jet length, hole height, diameter and disk radius, respectively,

Re = ρ U R μ 0 , Rb = U Ω R , Fr = U g R , Ra = 4 ε R σ θ 3 R ρ c p U D , = L R , h = H R , ϵ = D R .

In addition, we introduce dimensionless quantities that also depend on local air flow data, similarly to the Nusselt number in (4)

A 1 = 4 μ 2 R π ρ ρ U 2 D 3 , A 2 = ρ U D μ , A 3 = 8 k R π ρ c p θ D 2 , A 4 = μ c p k .

Here, A 4 is the Prandtl number of the air flow. To make (5) dimensionless we use the following reference values:

s 0 = L , r 0 = R , κ 0 = R 1 , u 0 = U , T 0 = θ , μ 0 = μ ( T 0 ) , n 0 = π μ 0 U D 2 / ( 4 R ) , m 0 = π μ 0 U D 4 / ( 16 R 2 ) .

We choose the disk radius R as macroscopic length scale in the scalings, since it is well known by the set-up. As for L, we consider jet lengths where the stresses are supposed to be vanished. In general, R and L are of same order such that the parameter ϵ can be identified with the slenderness ratio δ of the jets, cf. Introduction. The last two scalings for n 0 and m 0 are motivated by the material laws and the fact that the mass flux is Q=πρU D 2 /4. Then, the dimensionless system for the stationary viscous rod has the form

1 s r ˘ = R T e 3 = τ ˘ , 1 s R = κ × R , 1 s κ = 1 3 μ κ n 3 + 4 3 μ u P 3 / 2 m , 1 s u = 1 3 μ u n 3 , 1 s n = κ × n + Re u ( κ × e 3 + 1 3 μ n 3 e 3 ) + 2 Re Rb ( R e 1 ) × e 3 + Re Rb 2 1 u R ( e 1 × ( e 1 × r ˘ ) ) + R e F r 2 1 u R e 1 Re A 1 u R F ( τ ˘ , A 2 1 u v ˘ rel ) , 1 s m = κ × m + 4 ϵ 2 n × e 3 + Re 3 μ ( u P 3 m 1 4 n 3 P 2 κ ) Re 4 Rb 1 u P 2 ( 1 3 μ R e 1 n 3 + κ × R e 1 ) Re 4 ( 1 u 2 P 2 ( u κ + 1 Rb R e 1 ) ) × ( u κ + 1 Rb R e 1 ) , 1 s T = Ra 1 u ( T ref 4 T 4 ) + A 3 ( T T ) Nu ( v ˘ rel v ˘ rel τ ˘ , π A 2 2 v ˘ rel u , A 4 ) ,
(7)

with

r ˘ ( 0 ) = ( h , 1 , 0 ) , R ( 0 ) = e 1 e 1 e 2 e 3 + e 3 e 2 , κ ( 0 ) = 0 , u ( 0 ) = 1 , T ( 0 ) = 1 , n ( 1 ) = 0 , m ( 1 ) = 0 .

Here, T ref and the air flow associated T and v ˘ rel are scaled with θ and U, respectively. System (7) contains the slenderness parameter ϵ (ϵ1) explicity in the equation for the couple m and is hence no asymptotic model of zeroth order. In the slenderness limit ϵ0, the rod model reduces to a string system and their solutions for ( r ˘ , τ ˘ ,u,N= n 3 ,T) coincide. Only these jet quantities are relevant for the two-way coupling, as they enter in the exchange functions. However, the simpler string system is not well-posed for all parameter ranges, [15, 16]. Thus, it makes sense to consider (7) as ϵ-regularized string system, [19]. We treat ϵ as moderate fixed regularization parameter in the following to stabilize the numerics, in particular we set ϵ=0.1.

For the numerical treatment of (7), systems of non-linear equations are set up via a Runge-Kutta collocation method and solved by a Newton method. The Runge-Kutta collocation method is an integration scheme of fourth order for boundary value problems, that is, s z=f(s,z)f:[a,b]× R n R n with g(z(a),z(b))=0. It is a standard routine in MATLAB 7.4 with adaptive grid refinement (solver bvp4c.m). Let a= s 0 < s 1 << s N =b be the collocation points in [a,b] with h i = s i s i 1 and denote z i =z( s i ). Then, the nonlinear system of (N+1) equations, S( z h )=0, for the discrete solution z h = ( z i ) i = 0 , , N is set up via

S 0 ( z h ) = g ( z 0 , z N ) = 0 , S i + 1 ( z h ) = z i + 1 z i h i + 1 6 ( f ( s i , z i ) + 4 f ( s i + 1 / 2 , z i + 1 / 2 ) + f ( s i + 1 , z i + 1 ) ) = 0 , z i + 1 / 2 = 1 2 ( z i + 1 + z i ) h i + 1 8 ( f ( s i + 1 , z i + 1 ) f ( s i , z i ) )

for i=0,,N1. The convergence and hence the computational performance of the Newton method depends crucially on the initial guess. Thus, we adapt the initial guess iteratively by help of a continuation strategy. We scale the drag function F with the factor C F 2 and the right-hand side of the temperature equation with C T and treat Re, Rb, Fr, C F and C T as continuation parameters. We start from the solution for (Re,Rb,Fr,, C F , C T )=(1,1,1,0.15,,0) which corresponds to an isothermal rod without aerodynamic forces that has been intensively numerically investigated in [19]. Its determination is straight forward using the related string model as initial guess. Note that we choose so small to ensure that the glass jet lies in the air flow domain. The actual continuation is then divided into three parts. First, (Re,Rb,Fr,CF) are adjusted, then C T and finally . In the continuation we use an adaptive step size control. Thereby, we always compute the interim solutions by help of one step and two half steps and decide with regard to certain quality criteria whether the step size should be increased or decreased.

4.2 Weak iterative coupling algorithm

The numerical difficulty of the coupling of glass jet and air flow computations, S jets and S air , results from the different underlying discretizations. Let I h denote the rod grid used in the continuation method and I Δ be an equidistant grid of step size Δs with respective jet data ΨΔ for data exchange. Moreover, let Ω h denote the finite volume mesh with the flow data Ψ , V for the cell V, (so the chosen mesh realizes the necessary averaging). For the air associated exchange functions, the flow data is linearly interpolated on I h . Precisely, the linear interpolation L with respect to r ˘ ( s j ), s j I h is performed over all VN( s j ), where N( s j ) is the set of the cell containing r ˘ ( s j ) and its direct neighbor cells,

f ˘ air ( s j ) F ˘ ( Ψ ( s j ) , L N ( s j ) [ Ψ , V ] ) , q air ( s j ) Q ( Ψ ( s j ) , L N ( s j ) [ Ψ , V ] ) .

For the jet associated exchange functions entering the finite volume scheme, we need the averaged jet information for every cell V Ω h . We introduce I Δ , V ={ s j I Δ | r ˘ ( s j )V} and | I V |=Δs| I Δ , V |, then the averaging E with respect to V is performed over the I Δ , V -associated data,

2 π | V | V r f ˆ jets ( x , r ) d x d r N | I V | | V | F ˆ ( E V [ Ψ Δ ] , Ψ , V ) , 2 π | V | V r q jets ( x , r ) d x d r N | I V | | V | Q ( E V [ Ψ Δ ] , Ψ , V ) .

The ratio | I V |/|V| can be considered as the jet length density for the cell V. In case of M jet representatives, we deal with I Δ , V , i and | I V , i | for i=1,,M. Consequently, we have I Δ , V = i = 1 M I Δ , V , i and | I V |= i = 1 M | I V , i |. Note, that the interpolation and averaging approximation strategies have the disadvantage that they are qualitatively different. Thus, momentum and energy conservation are only ensured for very fine resolutions.

Summing up, the algorithm that we use to couple glass jet S jets and air flow S air computations has the form:

Algorithm 1 Generate flow mesh Ω h

Perform flow simulation S air without jets to obtain Ψ ( 0 )

Initialize k=0

Do

  • Compute: Ψ i ( k ) = S jets ( Ψ ( k ) ) for i=1,,M where flow data is linearly interpolated on I h

  • Interpolate jet data on equidistant grid I Δ

  • Find for every cell V in Ω h the relevant rod points I Δ , V and average the respective data

  • Compute: Ψ ( k + 1 ) = S air ( Ψ ( k ) )

  • Update: k=k+1

while Ψ ( k ) Ψ ( k 1 ) >tol

Remark 2 From the technical point of view, the efficient management of the simulation and coupling routines is quite demanding. In a preprocessing step we generate the finite volume mesh Ω h via the software Gambit and save it in a file that is available for FLUENT and MATLAB. The program of Algorithm 1 is then realized with FLUENT as master tool. After the air flow simulation FLUENT starts MATLAB. MATLAB governs the parallelization of the jets computation via MATLAB executables. Collecting the jets information, it provides the averaged jets data on Ω h in a file. FLUENT reads in this data and performs a new air flow simulation with immersed jets.

5 Results

In this section we illustrate the applicability of our asymptotic coupling framework to the given rotational spinning process. We show the convergence of the weak iterative coupling algorithm and discuss the effects of the fluid-fiber-interactions.

For all air flow simulations we use the same finite volume mesh Ω h whose refinement levels are initially chosen according to the unperturbed flow structure, independently of the glass jets. This implies a very fine resolution at the injector of the turbulent cross flow which is coarsen towards the centrifugal disk. For mesh details see Figure 4. The turbulent intensity is visualized in Figure 5. As expected it is high at the injector and moderate in the remaining flow domain. In particular, it is less than 2% in the region near the centrifugal disk where the glass jets will be presumably located. Thus, we neglect turbulence effects on the jets dynamics in the following. However, note that such effects can be easily incorporated by help of stochastic drag models [24, 37, 40] that are based on RANS turbulence descriptions (for example, k-ϵ model or k-ω model). For the jet computations the grid I h is automatically generated and adapted by the continuation method in every iteration. To ensure that sufficient jet points lie in each flow cell and a proper data exchange is given we use an equidistant grid I Δ with appropriate step size Δs (at minimum 2 jet points per interacting flow cell).

Figure 4
figure 4

Finite volume mesh Ω h for air flow computations (mesh detail).

Figure 5
figure 5

Turbulent intensity of the air flow.

The weak iterative coupling algorithm is fully automated. Each iteration starts with the same initialization. There is no parameter adjustment. The algorithm turns out to be very robust and reliable in spite of coarse flow meshes. For our set-up an air flow simulation takes around 30 minutes CPU-time, and the computation of a single jet takes approximatively just as long. The algorithm converges within 12-14 iterations. Figure 6 shows the relative L 2 -error of all jet curve components over the number of iterations k, that is,

i = 1 M r ˘ j , i ( k ) r ˘ j , i f L 2 ( I i ) r ˘ j , i f L 2 ( I i ) , with r ˘ j , i f final solution , j = 1 , 2 , 3 .
Figure 6
figure 6

Convergence of the weak iterative coupling algorithm. Relative error of all M curve coordinates in L 2 (I)-norm over number of iterations, plotted in logarithmic scale.

The effects of the fluid-fiber interactions and the necessity of the two-way coupling procedure for the rotational spinning process can be concluded from the following results. Figure 7 shows the swirl velocity of the air flow and the location of the immersed glass jets over the iterations. In the unperturbed flow without the glass jets there is no swirl velocity. In fact, the presence of the jets cause the swirl velocity, since the jets pull the flow with them. Moreover, the jets deflect the downwards directed burner flow, as seen in Figures 8 and 9. The jets behavior looks very reasonable. Trajectories and positions are as expected. Furthermore, their properties, that is, velocity u and temperature T, correspond to the axial flow velocity and flow temperature, which implies a proper momentum and heat exchange. For jet details we refer to Figures 10, 11, 12 and 13. They show the influence of the spinning rows. The jet representative of the highest spinning row is warmer than the one of the lowest row which implies better stretching capabilities. It is also faster and hence thinner (A= u 1 ). This certainly comes from the fact that the highest jet is longer affected by the fast hot burner flow. However, in view of quality assessment, slenderness and homogeneity of the spun fiber jets play an important role. This requires the optimal design of the spinning conditions, for example, different nozzle diameters or various distances between spinning rows. But for this purpose, also the melting regime has to be taken into account in modeling and simulation which is left to future research.

Figure 7
figure 7

Illustration of iterative coupling procedure. Iteration results for air swirl velocity and immersed glass jets (plotted as white curves).

Figure 8
figure 8

Final simulation result. Glass jets and air flow in given rotational spinning process. The color map visualizes the axial velocity of the air flow. In addition, the immersed M=35 glass jet representatives are colored with respect to their corresponding quantity u. For temperature information see Figure 9. Moreover, the dynamics and properties of the highest and lowest jets are shown in detail in Figures 10, 11, 12 and 13.

Figure 9
figure 9

Final simulation result. Glass jets and air flow in given rotational spinning process. The color maps visualize axial velocity and temperature of the air flow, respectively. In addition, the immersed M=35 glass jet representatives are colored with respect to their corresponding quantity T. For velocity information see Figure 8. Moreover, the dynamics and properties of the highest and lowest jets are shown in detail in Figures 10, 11, 12 and 13.

Figure 10
figure 10

Dynamics of the jets emerging from the highest and lowest spinning rows of the centrifugal disk. Side view of the plant. See associated Figures 11, 12 and 13.

Figure 11
figure 11

Dynamics of the jets emerging from the highest and lowest spinning rows of the centrifugal disk. Top view of the plant. See associated Figures 10, 12 and 13.

Figure 12
figure 12

Dynamics of the jets emerging from the highest and lowest spinning rows of the centrifugal disk. Jets velocity u(s), s[0,L]. See associated Figures 10, 11 and 13.

Figure 13
figure 13

Dynamics of the jets emerging from the highest and lowest spinning rows of the centrifugal disk. Temperature T(s), s[0,L]. See associated Figures 10, 11 and 12.

6 Conclusion

The optimal design of rotational spinning processes for glass wool manufacturing involves the simulation of ten thousands of slender viscous thermal glass jets in fast air streams. This is a computational challenge where direct numerical methods fail. In this paper we have established an asymptotic modeling concept for the fluid-fiber interactions. Based on slender-body theory and homogenization it reduces the complexity of the problem enormously and makes numerical simulations possible. Adequate to problem and model we have proposed an algorithm that weakly couples air flow and glass jets computations via iterations. It turns out to be very robust and converges to reasonable results within few iterations. Moreover, the possibility of combining commercial software and self-implemented code yields satisfying efficiency off-the-shelf. The performance might certainly be improved even more by help of future studies. Summing up, our developed asymptotic coupling framework provides a very promising basis for future optimization strategies.

In view of the design of the whole production process the melting regime must be taken into account in modeling and simulation. Melting and spinning regimes influence each other. On one hand the conditions at the spinning rows are crucially affected by the melt distribution in the centrifugal disk and the burner air flow, regarding, for example, cooling by mixing inside, aerodynamic heating outside. On the other hand the burner flow and the arising heat distortion of the disk are affected by the spun jet curtains. This obviously demands a further coupling procedure.

References

  1. Klar A, Marheineke N, Wegener R: Hierarchy of mathematical models for production processes of technical textiles. Z. Angew. Math. Mech. 2009, 89: 941–961. 10.1002/zamm.200900282

    Article  MATH  MathSciNet  Google Scholar 

  2. Pearson JRA: Mechanics of Polymer Processing. Elsevier, London; 1985.

    Google Scholar 

  3. Decent SP, King AC, Simmons MJH, Parau EI, Wallwork IM, Gurney CJ, Uddin J: The trajectory and stability of a spiralling liquid jet: Viscous theory. Appl. Math. Model. 2009,33(12):4283–4302. 10.1016/j.apm.2009.03.011

    Article  MATH  MathSciNet  Google Scholar 

  4. Parau EI, Decent SP, King AC, Simmons MJH, Wong D: Nonlinear viscous liquid jets from a rotating orifice. J. Eng. Math. 2006, 57: 159–179.

    Article  MathSciNet  Google Scholar 

  5. Marheineke N, Wegener R: Asymptotic model for the dynamics of curved viscous fibers with surface tension. J. Fluid Mech. 2009, 622: 345–369.

    Article  MATH  MathSciNet  Google Scholar 

  6. Antman SS: Nonlinear Problems of Elasticity. Springer, New York; 2006.

    Google Scholar 

  7. Entov VM, Yarin AL: The dynamics of thin liquid jets in air. J. Fluid Mech. 1984, 140: 91–111. 10.1017/S0022112084000525

    Article  MATH  Google Scholar 

  8. Yarin AL: Free Liquid Jets and Films: Hydrodynamics and Rheology. Longman, New York; 1993.

    MATH  Google Scholar 

  9. Cummings LJ, Howell PD: On the evolution of non-axissymetric viscous fibres with surface tension inertia and gravity. J. Fluid Mech. 1999, 389: 361–389.

    Article  MATH  MathSciNet  Google Scholar 

  10. Dewynne JN, Howell PD, Wilmott P: Slender viscous fibers with inertia and gravity. Q. J. Mech. Appl. Math. 1994, 47: 541–555. 10.1093/qjmam/47.4.541

    Article  MATH  MathSciNet  Google Scholar 

  11. Dewynne JN, Ockendon JR, Wilmott P: A systematic derivation of the leading-order equations for extensional flows in slender geometries. J. Fluid Mech. 1992, 244: 323–338. 10.1017/S0022112092003094

    Article  MATH  MathSciNet  Google Scholar 

  12. Decent SP, King AC, Wallwork IM: Free jets spun from a prilling tower. J. Eng. Math. 2002, 42: 265–282. 10.1023/A:1016127207992

    Article  MATH  MathSciNet  Google Scholar 

  13. Panda S, Marheineke N, Wegener R: Systematic derivation of an asymptotic model for the dynamics of curved viscous fibers. Math. Methods Appl. Sci. 2008, 31: 1153–1173. 10.1002/mma.962

    Article  MATH  MathSciNet  Google Scholar 

  14. Wallwork IM, Decent SP, King AC, Schulkes RMSM: The trajectory and stability of a spiralling liquid jet. Part 1. Inviscid theory. J. Fluid Mech. 2002, 459: 43–65.

    Article  MATH  MathSciNet  Google Scholar 

  15. Götz T, Klar A, Unterreiter A, Wegener R: Numerical evidence for the non-existence of solutions to the equations describing rotational fiber spinning. Math. Models Methods Appl. Sci. 2008,18(10):1829–1844. 10.1142/S0218202508003200

    Article  MATH  MathSciNet  Google Scholar 

  16. Arne W, Marheineke N, Meister A, Wegener R: Numerical analysis of Cosserat rod and string models for viscous jets in rotational spinning processes. Math. Models Methods Appl. Sci. 2010,20(10):1941–1965. 10.1142/S0218202510004738

    Article  MATH  MathSciNet  Google Scholar 

  17. Hlod A, Aarts ACT, van de Ven AAF, Peletier MA: Mathematical model of falling of a viscous jet onto a moving surface. Eur. J. Appl. Math. 2007, 18: 659–677.

    Article  MATH  MathSciNet  Google Scholar 

  18. Hlod, A., Aarts, A.C.T., van de Ven, A.A.F., Peletier, M.A.: Three flow regimes of viscous jet falling onto a moving surface. ArXiv:0811.2574 (2008)

    Google Scholar 

  19. Arne, W., Marheineke, N., Wegener, R.: Asymptotic transition of Cosserat rod to string models for curved viscous inertial jets. Math. Mod. Meth. Appl. Sci. 21(10) (2011, in press)

    Google Scholar 

  20. Ribe NM: Coiling of viscous jets. Proc. R. Soc., Math. Phys. Eng. Sci. 2004, 2051: 3223–3239.

    Article  MathSciNet  Google Scholar 

  21. Ribe NM, Habibi M, Bonn D: Stability of liquid rope coiling. Phys. Fluids 2006., 18:

    Google Scholar 

  22. Chiu-Webster S, Lister JR: The fall of a viscous thread onto a moving surface: A ‘fluid-mechanical sewing machine’. J. Fluid Mech. 2006, 569: 89–111.

    Article  MATH  MathSciNet  Google Scholar 

  23. Ribe NM, Lister JR, Chiu-Webster S: Stability of a dragged viscous thread: Onset of ‘stitching’ in a fluid-mechanical ‘sewing machine’. Phys. Fluids 2006., 18:

    Google Scholar 

  24. Marheineke N, Wegener R: Modeling and application of a stochastic drag for fiber dynamics in turbulent flows. Int. J. Multiph. Flow 2011, 37: 136–148. 10.1016/j.ijmultiphaseflow.2010.10.001

    Article  Google Scholar 

  25. Tomotika S, Aoi T, Yosinobu H: On the forces acting on a circular cylinder set obliquely in a uniform stream at low values of Reynolds number. Proc. R. Soc., Math. Phys. Eng. Sci. 1953,219(1137):233–244. 10.1098/rspa.1953.0142

    Article  MATH  MathSciNet  Google Scholar 

  26. Batchelor GK: Slender-body theory for particles of arbitrary cross-section in Stokes flow. J. Fluid Mech. 1970,44(3):419–440. 10.1017/S002211207000191X

    Article  MATH  MathSciNet  Google Scholar 

  27. Cox RG: The motion of long slender bodies in a viscous fluid. Part 1. General theory. J. Fluid Mech. 1970,44(4):791–810. 10.1017/S002211207000215X

    Article  MATH  Google Scholar 

  28. Taylor GI: Analysis of the swimming of long and narrow animals. Proc. R. Soc., Math. Phys. Eng. Sci. 1952, 214: 158–183. 10.1098/rspa.1952.0159

    Article  MATH  Google Scholar 

  29. Sumer BM, Fredsoe J: Hydrodynamics Around Cylindrical Structures. World Scientific Publishing, London; 2006.

    MATH  Google Scholar 

  30. Schlichting H: Grenzschicht-Theorie. Verlag G. Braun, Karlsruhe; 1982.

    MATH  Google Scholar 

  31. Schewe G: On the force fluctuations acting on a circular cylinder in cross-flow from subcritical up to transcritical Reynolds numbers. J. Fluid Mech. 1983, 133: 265–285. 10.1017/S0022112083001913

    Article  Google Scholar 

  32. Tritton DJ: Experiments on the flow past circular cylinder at low Reynolds number. J. Fluid Mech. 1959, 6: 547–567. 10.1017/S0022112059000829

    Article  MATH  Google Scholar 

  33. VDI-Wärmeatlas. Springer, Düsseldorf; 2006.

  34. Ribe NM, Huppert HE, Hallworth MA, Habibi M, Bonn D: Multiple coexisting states of liquid rope coiling. J. Fluid Mech. 2006, 555: 275–297.

    Article  MATH  MathSciNet  Google Scholar 

  35. Ferziger JH, Perić M: Computational Methods for Fluid Dynamics. Springer, Berlin; 2002.

    Book  MATH  Google Scholar 

  36. Zdravkovich MM: Flow Around Circular Cylinders, Vol 1: Fundamentals. Oxford University Press, New York; 1997.

    Google Scholar 

  37. Marheineke N, Wegener R: Fiber dynamics in turbulent flows: General modeling framework. SIAM J. Appl. Math. 2006,66(5):1703–1726. 10.1137/050637182

    Article  MATH  MathSciNet  Google Scholar 

  38. Hoerner, S.F.: Fluid-dynamic Drag. Practical Information on Aerodynamic Drag and Hydrodynamic Resistance. Published by the author. Obtainable from ISVA (1965)

    Google Scholar 

  39. Mahadevan L, Keller JB: Coiling of flexible ropes. Proc. R. Soc., Math. Phys. Eng. Sci. 1996, 452: 1679–1694. 10.1098/rspa.1996.0089

    Article  MATH  MathSciNet  Google Scholar 

  40. Marheineke N, Wegener R: Fiber dynamics in turbulent flows: specific Taylor drag. SIAM J. Appl. Math. 2007, 68: 1–23. 10.1137/06065489X

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors would like to acknowledge their industrial partner, the company Woltz GmbH in Wertheim, for the interesting and challenging problem. This work has been supported by German Bundesministerium für Bildung und Forschung, Schwerpunkt ‘Mathematik für Innovationen in Industrie und Dienstleistungen’, Projekt 03MS606 and by German Bundesministerium für Wirtschaft und Technologie, Förderprogramm ZIM, Projekt AUROFA 114626.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Nicole Marheineke.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The success of this study is due to the strong and fruitful collaboration of all authors. Even in details it is a joint work. However, special merits go to WA for the numerical analysis of Cosserat rods; to NM for modeling, investigating the asymptotic coupling concept and drafting the manuscript; to JS for conceptualizing and implementing the weak coupling software, performing the simulations and designing the visualizations; and to RW for developing the model framework, investigating the asymptotic coupling concept and implementing the continuation method for the jets. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Arne, W., Marheineke, N., Schnebele, J. et al. Fluid-fiber-interactions in rotational spinning process of glass wool production. J.Math.Industry 1, 2 (2011). https://doi.org/10.1186/2190-5983-1-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2190-5983-1-2

Keywords