Abstract

Air-jet loom is a textile machine designed to drive the long fiber using a combination flow of high-pressure air from a main nozzle and a series of assistant nozzles. To make the suggestion of how to make the fiber fly with high efficiency and stability in the jet flow, in which vortices also have great influence on fiber movement, the large eddy simulation method was employed to obtain the transient flow field of turbulent jet, and a bead-rod chain fiber model was used to predict long flexible fiber motion. The fluctuation and velocity of fibers with different linear densities in jet flow were studied numerically. The results show that the fluctuation amplitude of a fiber with a linear density of 0.5 × 10−5 kg·m−1 is two times larger than that of a fiber with a linear density of 2.0 × 10−5 kg·m−1. The distance of the first assistant nozzle from the main nozzle should be less than 120 mm to avoid collision between the fiber and the loom. The efficient length of the main nozzle to carry the fiber flying steadily forward is about 100–110 mm. For fibers with a linear density of 0.5 × 10−5 kg·m−1, it is suggested that the distance of the first assistant nozzle from the main nozzle is about 110 mm. With the increase of fiber linear density, the distance could be appropriately increased to 140 mm. The simulation results provide an optimization option for the air-jet loom to improve the energy efficiency by reasonably arranging the first assistant nozzle.

1. Introduction

Carried flight of a long, flexible fiber in a jet flow field is a fundamental science problem for air-jet loom [1], which is designed to drive the yarn using a combination flow of high-pressure air from a main nozzle and a series of assistant nozzles. Different from other pipeline pneumatic conveying systems, the jet flow of the main nozzle will rapidly decelerate in a free area and cannot guarantee a sustained high-speed flight of the fiber. Therefore, assistant nozzles must be added to the system. On the other hand, the presence of the reed groove requires that the fluctuation of the fiber must be less than a certain range. How to make the fiber fly forward rapidly and steadily is always a core issue of fiber flight and is jointly determined by two aspects: the fiber model and the flow simulation.

In general, a long fiber can be modeled as a rigid or flexible fiber for different dynamic research purposes. Yamamoto and Matsuoka [2] proposed a method for simulating the dynamic behavior of rigid and flexible fibers in a flow field. The fiber is regarded as made up of spheres that are lined up and bound to each neighbor. Each pair of bonded spheres can stretch, bend, and twist by changing the bond distance, bond angle, and torsion angle between spheres, respectively. For the rigid fiber, the computed period of rotation and distribution of orientation angle agree with those calculated using Jeffery’s equation with an equivalent ellipsoidal aspect ratio. For the flexible fiber, the period of rotation decreases rapidly with the growth of bending deformation of the fiber and rotation orbits deviate from a circular one for the rigid fiber. De Meulemeester et al. [3] tried to study the dynamic properties of fiber movement and developed a one-dimensional mathematical model, in which the behavior of the fiber is described by Newton’s second law. Lindström and Uesaka [4] proposed a model for flexible fibers in viscous fluid flow; namely, the fibers are modeled as chains of fiber segments which interact with the fluid through viscous and dynamic drag forces. Fiber segments, from the same or different fibers, interact with each other through normal, frictional, and lubrication forces. The simulations using the proposed model successfully reproduced the different regimes of motion for thread-like particles that range from rigid fiber motion to complicated orbiting behavior including coiling and self-entanglement.

Vahidkhah and Abdollahi [5] used the lattice Boltzmann method (LBM) to solve the Newtonian flow field and the immersed boundary method (IBM) to simulate the deformation of the flexible fiber interacting with the flow. The variations of the fiber length during the simulation time for different values of stretching constant are studied. Kabanemi and Hétu [6] carried out a direct simulation study to analyze the effect of fiber rigidity on fiber motion in simple shear flow. The fiber is modeled as a series of rigid spheres connected by stiff springs, which is similar to that used by Yamamoto and Matsuoka [2]. The model correctly predicts the orbit period of fiber rotation, as well as the trend of critical flow strength, versus fiber aspect ratio, during which the breakage occurs in simple shear flow. Meirson and Hrymak [7] extracted rotational friction coefficients from Jeffery’s model, created a general case long flexible fiber orientation model, and applied it in a simple shear flow. Nan et al. [8] presented a linear viscoelastic sphere-chain model based on the discrete element method to quantify the material damping of deformed flexible fibers. A correlation is formulated to quantify the relationship between the damping coefficient of the local bond and that of the flexible fiber. Meulemeester et al. [9] developed a three-dimensional mathematical model of the yarn. The three-dimensional model for the weft insertion on air-jet looms has been successfully tested.

On the other side, as we have mentioned, fiber flying in fluid is also affected by turbulent flow characteristics. Kim et al. [10] analyzed the flow in an air-jet loom by using a time-accurate characteristic-based upwind flux-difference splitting compressible Navier–Stokes method. The unsteady pressure and Mach number behavior along the center line of the main nozzle were analyzed. Andrić et al. [11] analyzed the dynamics of individual flexible fibers in a turbulent flow, the direct numerical simulation of the incompressible Navier–Stokes equations is used to describe the fluid flow in a plane channel, and a one-way coupling is considered between the fibers and the fluid phase. They found that the fiber motion is primarily governed by velocity correlations of the flow fluctuations. In addition, they reported that there is a clear tendency of the thread-like fibers to evolve into complex geometrical configurations in a turbulent flow field, and the fiber inertia has a significant impact on reorientation timescales of fibers suspended in a turbulent flow field. Jin et al. [12] simulated the turbulent flow by solving the Reynolds-averaged Navier–Stokes (RANS) equations and conducted the three-dimensional numerical simulation of the movement of the flexible body. The numerical results show that an unconstrained flexible body would turn over forward along the airflow’s diffusion direction, while a constrained flexible body in the flow field will make a periodic rotation motion along the axis of the flexible body, and the bending deformation is more obvious than that of unconstrained flexible body. Yang et al. [13] studied the two-way coupling turbulent model and rheological properties for fiber suspension in the contraction based on the RANS simulation.

Pei and Yu [14] studied the motional characteristics of the flexible fibers in the airflow inside the Murata vortex spinning (MVS) nozzle. A two-dimensional fluid structure interaction (FSI) model combined with the fiber wall contact is introduced to simulate a single fiber moving in the airflow inside the MVS nozzle. The model is solved using a finite element code ADINA. Based on their simulation results, the formation principle and the influence of some nozzle parameters on the tensile property of the MVS fiber were discussed.

More and more investigators put their interest on the air-jet loom [1518] and other fluid machineries [19, 20], but due to the computation cost gap between scientific researches and engineering needs, only a few of them considered the influence of turbulent fluctuation [21]. The jet flow caused by the main nozzle is a typical free shear turbulence, and there are strong vortices which play an important role on the momentum and energy transport of flow field [22, 23]. The vortices also have great influence on the fiber movement. Large eddy simulation (LES) is a believed mathematical model for turbulence vortex simulation, by directly calculating the large-scale flow motion.

In this paper, we simplified the two-phase flow system and were able to employ the LES method to simulate the development of a vortex in jet flow and the Lagrangian bead-rod model to give the time evolution of a long flexible fiber distribution with different linear densities. Then, the fiber fluctuation and the velocity were discussed to make the suggestion of how to make the fiber fly with high efficiency and stability in an air-jet loom.

2. Mathematical Formulation

2.1. Fluid Flow

The 2D jet flow is shown in Figure 1, in which x and y are the streamwise and cross-stream directions, respectively. The width of the nozzle D is 3.5 mm, the fluid velocity at the nozzle U0 is 240 m/s, and the flow Reynolds number .

The LES equations governing the jet flow obtained by filtering the Navier–Stokes equations are as follows [24]:where is the fluid density, is the filtered velocity, is the filtered pressure, index , is taken as 1, 2 and refers to the , , and is the dynamic viscosity.

The SGS (subgrid stress) tensor and Smargorinsky–Lilly model [24] which are based on the mixing length hypothesis are used to calculate the SGS stress.

2.2. Flexible Fiber Model

A single flexible long fiber is modeled as a bead-rod chain, which is similar to that used by Guo et al. [25]. The fiber model is composed of N beads, which are connected by massless rods (Figure 2). Only the beads are affected by forces, and the rods maintain the configuration of the fiber. Using the model, the chain is allowed to be stretched by changing the distance of adjacent rods.

If the distance between adjacent beads is not equal to the equilibrium distance, the stretching restoring force exerted on the bead will bewhere is the bead radius, is the elastic modulus, is the equilibrium distance, and is the distance variation, which is equal to the transient distance of each two adjacent beads subtracted by the equilibrium distance.

When immersed in the unsteady flow field which is calculated in the last section, the fiber is subjected to hydrodynamic forces, which are also changed with time. In this paper, only drag force is considered in order to decrease the computational cost. Other hydrodynamic forces such as Basset history term, additional mass, slip-rotational lift force, and fluid inertia are negligible. For bead , the drag force is contributed by fiber sections and . It can be calculated as follows:where and are the drag forces acting on bead , which are devoted to the fiber section and , respectively. The drag acting on the fiber section can be expressed as follows:where and are the fluid and fiber velocities at the mass center of the fiber section , and is the drag coefficient associated with Reynolds number, which can be represented as follows:where is the Reynolds number of equivalent spherical and is a constant between 0.4 and 0.45. According to Newton’s second law, the equations of motion for the bead that constitute the fiber are as follows:where is the mass of the bead and is the position vector of the bead .

3. Results and Discussion

3.1. Computation Conditions

The computational domain covers the nozzle width D = 3.5 mm, fully developed boundary conditions at the outlet, and static surrounding environment conditions are assumed. The local mesh is shown in Figure 3. Fiber motion is solved in one-way coupling between the fiber and flow; fiber-wall and fiber-fiber interactions are neglected.

3.2. Method Verification
3.2.1. Jet Flow Field Verification

According to the experimental results of [26, 27], the streamwise velocity satisfies the self-preservation distribution at and the profile is expressed as a Gaussian curve. The self-preservation profiles of the mean velocity are shown in Figure 4, where the data are normalized by the centerline mean velocity and the half-width , which is the distance between the position where the streamwise velocity and the center line of the jet. The figure shows that the air flow which is mentioned in Section 2.1 has been successfully simulated.

3.2.2. Flexible Fiber Model Verification

As the limitation of the visible area and the resolution of the high-speed camera, an experiment about the motion of microscale fixed flexible fiber was carried out to verify the long flexible fiber model. The fixed fiber length Lf  = 100 mm, and the linear density (mass per unit length) . As shown in Figure 5, the noncontact test bed consists of an air jet loom, high-speed camera, and professional image analysis software Image Pro-Plus. The fiber is fixed at the nozzle so its motion state can be recorded by the high-speed camera after each parameter adjustment of the loom.

The transient fiber shape comparison of numerical simulation results and that recorded by the high-speed camera are shown in Figure 6. It can be seen that the It can be seen that the free end of the fiber fluctuates significantly all the time, but the fluctuation of the fiber segment close to the main nozzle is not obvious. fluctuates significantly all the time, but the fluctuation of the fiber segment close to the main nozzle is not obvious. This is due to the turbulence characteristics of the jet flow, fiber flexibility, and the shape characteristics of the fiber.

There are some less significant differences between the simulation models used in this paper and the actual fiber properties: the bending modulus of the fiber, which is much lower than the elasticity modulus, is not taken into account; the hydrodynamic forces such as Basset force and slip-rotational lift force, which are relative small compared to the drag force, are neglected as well. Despite all the simplifications, the experimental and simulation results are still in relatively good agreement (Figure 6). Therefore, the following simulation of the flight fiber is reasonable and credible.

3.3. Discussions and Analysis

When a micro- or nanoscale free fiber is being carried by the high-velocity airflow after being injected into the jet flow field, as described above, the flexible fiber model assumes that the fiber mass is concentrated at a series of beads, and the intervals between every two beads can vary as a result of the force-displacement balance on the beads. Subsequently, the transient equivalent radius of the fiber varies too; they can even be reduced from microscale to nanoscale. Fiber linear density is another important characteristic and has a significant effect on calculating the force acted on the beads which would determine the fiber motion in the jet flow. So the fiber motion with different linear densities of , 1.0 × 10−5, 1.5 × 10−5, and 2.0 × 10−5 kg·m−1 are separately studied.

The initial velocity of the fiber is . The initial length of static free fiber is Lf = 40 mm. The initial bead radius is given as R = 100 . The two ends of the fibers are marked as end A and B, respectively. At the initial stage, end A is at the position x = 40 mm and end B is at x = 0 mm (Figure 7). The motional characteristic of the fiber is stable at the initial stage. The farther from the nozzle, the lower the air velocity. So, the velocity of end B will be faster than that of end A after a period of time. As a result, fiber bending deformation is observed.

3.3.1. Effect of Fiber Linear Density on Fiber Fluctuation

Due to different densities, the time of fiber flying across the simulation domain is different. The fiber movement is analyzed until the end A of fiber arrives at the calculating boundary.

The motion and bending deformation of fibers with different linear densities over time are shown in Figure 8. As the jet flow goes on, the vortex keeps rolling-up, transporting and mixing with each other (Figure 1), and the fluctuation of vortex driven on fiber becomes more and more strong. By comparing Figures 8(a)8(d), it can be seen that the fiber fluctuation relating with linear density is more obvious than that of fiber with linear density , which shows that the smaller linear density causes the more unstable fiber motion. This is because the fiber transverse acceleration generates from the velocity difference between transverse velocity of flow field and that of fiber. Moreover, the smaller linear density results in the larger transverse acceleration and the larger transverse acceleration lead to the more unstable fiber motion.

The transverse velocity of end A of fibers with different linear densities is shown in Figure 9. The fiber transverse velocity is 0 m/s at the initial time, and it becomes larger with the increasing of running time. The end transverse velocity of fiber with linear density increases to about 2 m/s after 0.0012 s. At the same time, the end transverse velocity of fiber with linear density increases to about 18 m/s. Thus, the smaller fiber linear density leads to the larger transverse velocity. In other words, the fiber transverse fluctuation becomes more obvious.

In order to learn more about the influence of air velocity on fiber fluctuations, several equidistant points are selected on the fiber (Figure 10). To represent the amplitude of fiber fluctuation, the fluctuation standard deviation of fiber is calculated using (7) when the end A of fiber arrives at 40 mm, 80 mm, 100 mm, 120 mm, and 140 mm (i.e., xA = 40 mm, 80 mm, 100 mm, 120 mm, and 140 mm) from the main nozzle:where is the vertical distance from the point to the center line and is the sum of the selected points in the fiber (n = 41). The fluctuation standard deviation is listed in Table 1.

As shown in Table 1, the fiber is affected by the airflow and the fiber fluctuation standard deviation increases gradually in the motion process with time. When the fiber with different linear densities arrives at the same location in jet flow, the smaller fiber linear density causes the greater fluctuation standard deviation. The fluctuation amplitude of fibers with is two times larger than that of fibers with . That is to say, the fiber with smaller linear density will fluctuate more obviously as it flies across the computation domain. As the cross section of the reed groove in the loom is about 5 mm × 5 mm, so when the fluctuating amplitude of the fiber end is larger than 2.5 mm, that is, after the head of fiber with arrived at xA = 120 mm and the head of fiber with arrived at xA = 140 mm, the fiber will much more likely impact with the groove if there is no another assistant jet flow. So it is strongly recommended that the distance of the first assistant nozzle from the main nozzle is less than 120 mm for the fiber whose linear density is less than 0.5 × 10−5 kg·m−1. With the increase of fiber linear density, the distance could be appropriately increased.

3.3.2. Effect of Fiber Linear Density on Fiber Velocity

The x-velocity () and distance from the nozzle (xA) of end A with different linear densities and time are shown in Figure 11. The corresponding flow axial velocity (at the positions where the end A is) is given in Figure 12. Within the running time of 0–0.0002 s, the fiber velocity increases obviously. Then, the fiber is in a state of slow acceleration during the time of 0.0002–0.0004 s. At the next stage, there is no longer a significant change on fiber velocity. After about 0.0008 s, the fiber is in a state of deceleration. It can be explained as follows. The fiber velocity is relatively low at the initial moment. When the end A velocity is around 70 m/s, the air velocity of flow field is around 160 m/s (Figure 12). There is a large velocity difference between the flow field and fiber. The fiber is in rapid acceleration because of the large drag force as listed previously. It means the fiber gets a lot of energy from the flow in the period of 0–0.0002 s, and the end A flies about 15 mm forward in this ultra-short time. Then, the velocity of flow field, where the fiber is, decreases rapidly, which we can see from Figure 12. During the time of 0.0002–0.0005 s, the velocity of airflow is larger than the fiber, but the velocity difference decreases. So the fiber is in a state of slow acceleration. As the velocity difference between the flow field and fiber keeps decreasing in 0.0005–0.0008 s, the influence of flow field on fiber velocity becomes unapparent and the fiber remains at the same velocity magnitude especially for the fibers with a relatively large linear density. At t = 0.0008 s, the fiber end A arrives around 100 mm away from the main nozzle, which we can see from Figure 11, and the velocity of airflow is lower than that of the fiber. The fiber cannot get energy from the flow any more for all the fibers simulated here. In other words, the efficient length of the main nozzle to carry the fiber flying rapidly forward is about 100–110 mm. So the addition of an assistant nozzle is suggested in this place.

Besides, the difference of velocity distribution is obvious with different fiber linear densities. According to Newton's second law, the lower linear density leads to greater fiber acceleration. So when the velocity of airflow is greater than that of the fiber, the greater fiber acceleration causes the greater velocity. But, when the velocity of airflow is lower than that of the fiber, the velocity of fiber with a lower density falls quickly. It means that the fiber with lower linear density is more sensitive to the flow, and the flow velocity change along the fiber will easily cause a nonuniform velocity distribution at different segments of the fiber.

The x-velocity of ends A and B of the fiber with different linear densities is shown in Figure 13. The velocity difference between end A and B of the fiber with liner density is about 35 m/s when the end A arrives at the calculating boundary, and the fiber with linear density is only about 10 m/s. The nonuniform velocity distribution leads to obvious bending deformation of the fiber (Figure 8(a)), which is undesirable for pneumatic conveying especially jet loom. In order to guarantee the stability of movement of the fiber with lower density, a more stable velocity distribution of flow is needed. On the other hand, the decrease of the distance between first the assistant nozzle and the main nozzle can reduce the axial velocity change and improve the stability of flow [28]. So, in order to make the flight more stable, for the fiber with lower linear density, the distance of first assistant nozzle from the main nozzle should be appropriately smaller. But, it is obviously worse for saving energy. Actually, in industrial production, people will try their best to add the spacing to make full use of the high-speed air jet flow.

The discussions above are summarized as follows. First, as we have discussed before, for fiber with linear density , if there is an assistant nozzle within 120 mm from the main nozzle, it will be highly likely for the fiber to impact with the groove. Second, the efficient length of the main nozzle to carry the fiber forward is less than 110 mm. Third, to make full use of the high-speed air jet flow, the distance between the main nozzle and the first assistant nozzle cannot be very small, although the smaller nozzle spacing could lead to a more stable fiber flight. Considering all the factors, there is only a narrow range of suggestion distance to set the first assistant nozzle. For fiber with a linear density , it is 110 mm, and when the fiber linear density increases, the distance could be appropriately increased to 140 mm.

4. Conclusions

To make the suggestion of how to make the fiber fly with high efficiency and stability in a jet flow, we employed the LES method to simulate the development of vortices in the jet flow and Lagrangian bead-rod model to give the time evolution of long flexible fiber distribution with different linear densities. The fluctuation and velocity of fiber in jet flow were then studied numerically, and the simulation results can provide an optimization option for the air-jet loom to improve the energy efficiency by reasonably arranging the first assistant nozzle. The results are as follows.(1)As the primary vortex rolls up, transports, and mixes with each other, the fiber fluctuation becomes stronger and the motion becomes more unstable as the linear density decreases.(2)The fluctuation amplitude of a fiber with is two times larger than that of a fiber with . It is as large as 3.34 mm when the fiber with arrives at xA = 120 mm.(3)The distance of the first assistant nozzle from the main nozzle should be less than 120 mm to avoid the collision between the fiber and the loom. With the increasing fiber linear density, the distance could be appropriately increased.(4)The efficient length of the main nozzle to carry the fiber flying steadily forward is about 100–110 mm. So an assistant nozzle should be added in this place.(5)To save energy, according to (2) and (3), the suggested distance between the main nozzle and the first assistant nozzle is 110 mm for the fiber with . When the fiber linear density increases, the distance could be appropriately increased to 140 mm.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 51576180 and 51676173), Fluid Engineering Innovation Team of Zhejiang Sci-Tech University (ZSTU) of China (Grant no. 11132932611309), National Scholarship Foundation of China, and 521 Talents Fostering Program Funding of Zhejiang Sci-Tech University of China.