Submit manuscript...
Journal of
eISSN: 2574-8114

Textile Engineering & Fashion Technology

Research Article Volume 1 Issue 5

Structural optimization of the main nozzle in an air-jet loom based on orthogonal test method

Xiangjun Lan, Dingding Liu, Zhihua Feng

Soochow University, China

Correspondence: Zhihua Feng, School of Mechanical and Electric Engineering, Soochow University, Suzhou 215021, China

Received: March 31, 2017 | Published: April 26, 2017

Citation: Lan X, Liu D, Feng. Structural optimization of the main nozzle in an air-jet loom based on orthogonal test method. J Textile Eng Fashion Techno. 2017;1(5):172-178. DOI: 10.15406/jteft.2017.01.00029

Download PDF

Abstract

Due to having more complicated inner structural features and obvious influence on the air-jet performance arising from the acceleration tube, the structural optimization of the main nozzle for an air-jet loom is still a difficult problem. Naturally, traditional design method for the nozzle will still prolong the time of type approval and increase the design cost. In this paper, orthogonal test method and numerical simulation are utilized to optimize the structure of subsonic acceleration zone; the optimal structure parameters matched with different air supply pressure are obtained. Results indicate that the external diameter of the nozzle needle has a negative influence on the airflow velocity whereas the diameter of the acceleration tube presents positive effect. Both of them play an important influence on the flow field of the main nozzle field and the other structural parameters have little influence on that. This research is also benefit for structural optimization of the main nozzle and energy conservation of an air-jet loom.

Keywords: air-jet loom, main nozzle, structural optimization, flow field, orthogonal test method, rectification process, flow field, jet performance, acceleration tube, numerical integration, adiabatic pipe flow, weft yarn, weft insertion, axial length, needle tip

Introduction

In an air-jet loom, the weft yarn is dragged by the high-speed air flow to pass through the shed to carry out the weft insertion process, where the air flow may reach a critical condition at the nozzle throat zone and even exceed sonic state at outside of the acceleration tube after modulation and rectification process. At the same time, the weft yarn has had a full interaction with high-speed air flow and is accelerated to the needed velocity. So the main nozzle is really the heart of an air-jet loom, and Studies on its pneumatic behaviors have been intermittently pursued over the past few decades. For instance, Ishida and Okajima1,2 experimentally measured the static pressure and air velocity distribution under different acceleration tube length and air supply pressure in the flow field of main nozzle. Their results showed that the attenuate law of air velocity along axial direction and radial direction is similar, and the air velocity of main nozzle exit is not relevant to acceleration tube length. Mohamed and Salama3,4 tested the air velocity and characteristic of the flow turbulence with different main nozzle structures, and they also got the relationship between various geometries of the main nozzle and the jet performance.5 Researched the effect of different acceleration tube geometries on flow field. Due to the complicated inner structure of main nozzle, it is difficult to directly place a sensor to measure the air flow.6 Studied the effect of the length of the acceleration tube and the shape of nozzle needle on flow field through two-dimensional numerical simulation.7−9 systematically researched the flow field characteristics of main nozzle by three-dimensional numerical simulation and some structural parameters of the main nozzle were also dealt with.

In this paper, the drag force of weft yarn is obtained by the theory of one dimensional steady adiabatic pipe flow and equation of gas state. The drag force under different types of air supply pressure is computed by numerical integration, which is an important measurement criterion for the jet performance of the main nozzle. Afterwards number of rectifier tank, external diameter of needle tip, distance between needle tip and rectifier tank, α face projected length and diameter of acceleration tube are optimized by the combination of both numerical simulation and orthogonal test method.

Design objective and measured criterion

Design objective

The main nozzle is a key component of weft insertion system in an air-jet loom. The air chamber and rectifier tank are used to collimate and accelerate the air flow parallel to weft insertion direction. At the same time, a conical sleeve is designed to form a subsonic acceleration zone. When the air flow passed through nozzle throat which is an annual slot, it may exceed sonic speed and produce a negative zone in weft injection region. So the weft yarn is inhaled into main nozzle through weft insertion passage. In the acceleration tube, the air flow is accelerated again so as to obtain a higher speed to obtain the needed kinetic energy. The rotation of advanced air-jet loom has reached 800~1000r/min with weft insertion speed of 3000m/min, and the weft yarn should fly across the shed within 20~30ms. So the main nozzle plays an important role in accelerating yarns and achieving sonic velocity at the main nozzle exit is always the design objective of it

Measurement criterion

Air velocity in the main nozzle: The weft yarn is accelerated by abundant friction with high-speed air flow, where the drag force is depended on the relative speed between yarn and jet flow. The corresponding formula is shown as:

F= 1 2 C f ρπd ( VU ) 2 L MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOrai abg2da9maalaaabaGaaGymaaqaaiaaikdaaaGaam4qamaaBaaabaGa amOzaaqabaGaeqyWdiNaeqiWdaNaamizamaabmaabaGaamOvaiabgk HiTiaadwfaaiaawIcacaGLPaaadaahaaqabeaacaaIYaaaaiaadYea aaa@4801@ (1)

where F MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOraa aa@3960@ is the drag force, C f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4qam aaBaaabaGaamOzaaqabaaaaa@3A69@ is the coefficient of drag force, ρ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaeqyWdi haaa@3A55@ is the air density, d is the weft yarn diameter, V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOvaa aa@3970@ is the air velocity, U MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamyvaa aa@396F@ is the yarn velocity, and L is the yarn length subject to air. From Equation (1) we can find the drag force mainly depends on relative speed and air density in main nozzle. Based on one-dimensional pipe flow theory, the velocity of main nozzle exit can reflect the inner flow condition of the main nozzle, which naturally is the most significant measurement criterion.

Negative pressure zone of the needle tip exit: As the result of high-speed jet flow, negative pressure zone is produced at the needle tip. The weft yarn is inhaled to main nozzle through weft guide passage and then it enters into weft ejection zone and acceleration zone. It is thus clear that a negative pressure zone is necessary for weft insertion in Figure 1.

Figure 1 Schematic of the weft insertion.

Drag force: The friction between yarn and airflow is the motive power for leaping and produces the so-called drag force, which will counteract the friction between yarn and the thread guide, the dynamic tension of ballooning and the static tension in weft accumulator. In order to analyze and calculate the force, we assumed that the weft yarn move along with the axis of acceleration tube, and the drag force exerted on infinitesimal segment is

dF= 1 2 C f ρπd ( VU ) 2 dl MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamizai aadAeacqGH9aqpdaWcaaqaaiaaigdaaeaacaaIYaaaaiaadoeadaWg aaqaaiaadAgaaeqaaiabeg8aYjabec8aWjaadsgadaqadaqaaiaadA facqGHsislcaWGvbaacaGLOaGaayzkaaWaaWbaaeqabaGaaGOmaaaa caWGKbGaamiBaaaa@49F3@ (2)

As the flow velocity at needle tip may exceed the speed of sound whereas the diameter of acceleration tube keeps constant, thus the air flow is a typical Fanno one. Then the differential equation of air velocity along the axial direction is

4fK D dl= 1 M 2 M 2 [ 1+ M 2 ( K1 )/2 ] d M 2 M 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfa4aaSaaae aacaaI0aGaamOzaiaadUeaaeaacaWGebaaaiaadsgacaWGSbGaeyyp a0tbaeqabeGaaaqaamaalaaabaGaaGymaiabgkHiTiaad2eadaahaa qabeaacaaIYaaaaaqaaiaad2eadaahaaqabeaacaaIYaaaamaadmaa baGaaGymaiabgUcaRiaad2eadaahaaqabeaacaaIYaaaamaabmaaba Gaam4saiabgkHiTiaaigdaaiaawIcacaGLPaaacaGGVaGaaGOmaaGa ay5waiaaw2faaaaaaeaadaWcaaqaaiaadsgacaWGnbWaaWbaaeqaba GaaGOmaaaaaeaacaWGnbWaaWbaaeqabaGaaGOmaaaaaaaaaaaa@52F4@ (3)

Where M is the Mach number, D is the acceleration tube diameter, f is the frictional coefficient, and K is the ratio of special heat.

The inlet flow ought to be the same as outlet flow based on the continuity relation, where D is constant. We obtain

ρV= ρ 2 V 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaeqyWdi NaamOvaiabg2da9iabeg8aYnaaBaaabaGaaGOmaaqabaGaamOvamaa BaaabaGaaGOmaaqabaaaaa@408B@ (4)

From the law of conservation of energy, the stagnation temperature of the nozzle exit is equal to that of air tank. Thus, we obtain

C P T 0 = C P T+ V 2 2 /2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4qam aaBaaabaGaamiuaaqabaGaamivamaaBaaabaGaaGimaaqabaGaeyyp a0Jaam4qamaaBaaabaGaamiuaaqabaGaamivaiabgUcaRiaadAfada qhaaqaaiaaikdaaeaacaaIYaaaaiaac+cacaaIYaaaaa@446A@ (5)

Where C P MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4qam aaBaaabaGaamiuaaqabaaaaa@3A53@ is specific heat at constant pressure and T 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamivam aaBaaabaGaaGimaaqabaaaaa@3A49@ is stagnation temperature. Temperature of acceleration tube exit T2 is represented as

T 2 = T 0 V 2 2 /2 C P MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamivam aaBaaabaGaaGOmaaqabaGaeyypa0JaamivamaaBaaabaGaaGimaaqa baGaeyOeI0IaamOvamaaDaaabaGaaGOmaaqaaiaaikdaaaGaai4lai aaikdacaWGdbWaaSbaaeaacaWGqbaabeaaaaa@4394@ (6)

From the ideal gas state equation (i.e. P=ρRT MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamiuai abg2da9iabeg8aYjaadkfacaWGubaaaa@3DE0@ ) and equations (4) and (6), air density of the acceleration tube is

ρ 2 = P 2 /[ R( T 0 V 2 2 / C P ) ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaeqyWdi 3aaSbaaeaacaaIYaaabeaacqGH9aqpcaWGqbWaaSbaaeaacaaIYaaa beaacaGGVaWaamWaaeaacaWGsbWaaeWaaeaacaWGubWaaSbaaeaaca aIWaaabeaacqGHsislcaWGwbWaa0baaeaacaaIYaaabaGaaGOmaaaa caGGVaGaam4qamaaBaaabaGaamiuaaqabaaacaGLOaGaayzkaaaaca GLBbGaayzxaaaaaa@4A76@ (7)

Where P 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamiuam aaBaaabaGaaGOmaaqabaaaaa@3A47@ and V 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOvam aaBaaabaGaaGOmaaqabaaaaa@3A4D@ are static pressure and air velocity of nozzle exit, respectively. Then from equations (2), (3), (4) and (7), we obtain

dF= πdD P 2 V 2 C f 8fKRC( T 0 V 2 2 /2 C P ) ( CMU ) 2 M 1 M 2 M 2 [ 1+ M 2 ( K1 )/2 ] d M 2 M MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamizai aadAeacqGH9aqpdaWcaaqaaiabec8aWjaadsgacaWGebGaamiuamaa BaaabaGaaGOmaaqabaGaamOvamaaBaaabaGaaGOmaaqabaGaam4qam aaBaaabaGaamOzaaqabaaabaGaaGioaiaadAgacaWGlbGaamOuaiaa doeadaqadaqaaiaadsfadaWgaaqaaiaaicdaaeqaaiabgkHiTiaadA fadaqhaaqaaiaaikdaaeaacaaIYaaaaiaac+cacaaIYaGaam4qamaa BaaabaGaamiuaaqabaaacaGLOaGaayzkaaaaaiabgkci3oaalaaaba WaaeWaaeaacaWGdbGaamytaiabgkHiTiaadwfaaiaawIcacaGLPaaa daahaaqabeaacaaIYaaaaaqaaiaad2eaaaGaeyOiGCBbaeqabeGaaa qaamaalaaabaGaaGymaiabgkHiTiaad2eadaahaaqabeaacaaIYaaa aaqaaiaad2eadaahaaqabeaacaaIYaaaamaadmaabaGaaGymaiabgU caRiaad2eadaahaaqabeaacaaIYaaaamaabmaabaGaam4saiabgkHi TiaaigdaaiaawIcacaGLPaaacaGGVaGaaGOmaaGaay5waiaaw2faaa aaaeaadaWcaaqaaiaadsgacaWGnbWaaWbaaeqabaGaaGOmaaaaaeaa caWGnbWaaWbaaeqabaaaaaaaaaaaaa@6F19@ (8)

Where R MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOuaa aa@396C@ is gas constant and C MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4qaa aa@395D@ is velocity of sound.

As the result of short length of acceleration tube, coefficient Cf can be taken as a constant, thus, we have

F= πdD C f 5.6fRC P 2 V 2 T 0 V 2 2 /2 C P ×[ F( M 2 )F( M 1 ) ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOrai abg2da9maalaaabaGaeqiWdaNaamizaiaadseacaWGdbWaaSbaaeaa caWGMbaabeaaaeaacaaI1aGaaiOlaiaaiAdacaWGMbGaamOuaiaado eaaaGaeyOiGC7aaSaaaeaacaWGqbWaaSbaaeaacaaIYaaabeaacaWG wbWaaSbaaeaacaaIYaaabeaaaeaacaWGubWaaSbaaeaacaaIWaaabe aacqGHsislcaWGwbWaa0baaeaacaaIYaaabaGaaGOmaaaacaGGVaGa aGOmaiaadoeadaWgaaqaaiaadcfaaeqaaaaacqGHxdaTdaWadaqaai aadAeadaqadaqaaiaad2eadaWgaaqaaiaaikdaaeqaaaGaayjkaiaa wMcaaiabgkHiTiaadAeadaqadaqaaiaad2eadaWgaaqaaiaaigdaae qaaaGaayjkaiaawMcaaaGaay5waiaaw2faaaaa@5EB1@ (9)

Wherem

F(M)= U 2 /3 M 3 +CU/ M 2 +(K+1) U 2 /2 C 2 /M+(K+1)UCln(M)+ [( K 2 1) U 2 2(K+1) C 2 ]/ 8(K1) ×arctan( (K1)/2 M) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOrai aacIcacaWGnbGaaiykaiabg2da9iabgkHiTiaadwfadaahaaqabeaa caaIYaaaaiaac+cacaaIZaGaamytamaaCaaabeqaaiaaiodaaaGaey 4kaSIaam4qaiaadwfacaGGVaGaamytamaaCaaabeqaaiaaikdaaaGa ey4kaSIaaiikaiaadUeacqGHRaWkcaaIXaGaaiykaiaadwfadaahaa qabeaacaaIYaaaaiaac+cacaaIYaGaeyOeI0Iaam4qamaaCaaabeqa aiaaikdaaaGaai4laiaad2eacqGHRaWkcaGGOaGaam4saiabgUcaRi aaigdacaGGPaGaamyvaiaadoeaciGGSbGaaiOBaiaacIcacaWGnbGa aiykaiabgUcaRGqaaiaa=bcacaGGBbGaaiikaiaadUeadaahaaqabe aacaaIYaaaaiabgkHiTiaaigdacaGGPaGaamyvamaaCaaabeqaaiaa ikdaaaGaeyOeI0IaaGOmaiaacIcacaWGlbGaey4kaSIaaGymaiaacM cacaWGdbWaaWbaaeqabaGaaGOmaaaacaGGDbGaai4lamaakaaabaGa aGioaiaacIcacaWGlbGaeyOeI0IaaGymaiaacMcaaeqaaiabgEna0k GacggacaGGYbGaai4yaiaacshacaGGHbGaaiOBaiaacIcadaGcaaqa aiaacIcacaWGlbGaeyOeI0IaaGymaiaacMcacaGGVaGaaGOmaaqaba GaamytaiaacMcaaaa@8090@ (10)

Substituting K into the above formula, we have

F(M)= U 2 /3 M 3 +CU/ M 2 +1.2 U 2 C 2 /M+2.4UCln(M)+(0.24 U 2 1.2 C 2 )× 5 arctan(M/ 5 ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOrai aacIcacaWGnbGaaiykaiabg2da9iabgkHiTiaadwfadaahaaqabeaa caaIYaaaaiaac+cacaaIZaGaamytamaaCaaabeqaaiaaiodaaaGaey 4kaSIaam4qaiaadwfacaGGVaGaamytamaaCaaabeqaaiaaikdaaaGa ey4kaSIaaGymaiaac6cacaaIYaGaamyvamaaCaaabeqaaiaaikdaaa GaeyOeI0Iaam4qamaaCaaabeqaaiaaikdaaaGaai4laiaad2eacqGH RaWkcaaIYaGaaiOlaiaaisdacaWGvbGaam4qaiGacYgacaGGUbGaai ikaiaad2eacaGGPaGaey4kaSIaaiikaiaaicdacaGGUaGaaGOmaiaa isdacaWGvbWaaWbaaeqabaGaaGOmaaaacqGHsislcaaIXaGaaiOlai aaikdacaWGdbWaaWbaaeqabaGaaGOmaaaacaGGPaGaey41aq7aaOaa aeaacaaI1aaabeaaciGGHbGaaiOCaiaacogacaGG0bGaaiyyaiaac6 gacaGGOaGaamytaiaac+cadaGcaaqaaiaaiwdaaeqaaiaacMcaaaa@6E62@ (11)

At the beginning of weft insertion, the velocity of weft yarn is equal to zero (i.e. U=0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamyvai abg2da9iaaicdaaaa@3B2F@ ), so the static weft drag force is

F= πdD C f C 5.6fR P 2 V 2 T 0 V 2 2 /2 C P ×{ 1 M 1 1 M 2 +1.2 5 [arctan( M 1 / 5 )arctan( M 2 / 5 )]} MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOrai abg2da9maalaaabaGaaeiWdiaadsgacaWGebGaam4qamaaBaaabaGa amOzaaqabaGaam4qaaqaaiaaiwdacaGGUaGaaGOnaiaadAgacaWGsb aaaiabgkci3oaalaaabaGaamiuamaaBaaabaGaaGOmaaqabaGaamOv amaaBaaabaGaaGOmaaqabaaabaGaamivamaaBaaabaGaaGimaaqaba GaeyOeI0IaamOvamaaDaaabaGaaGOmaaqaaiaaikdaaaGaai4laiaa ikdacaWGdbWaaSbaaeaacaWGqbaabeaaaaGaey41aqRaai4Eamaala aabaGaaGymaaqaaiaad2eadaWgaaqaaiaaigdaaeqaaaaacqGHsisl daWcaaqaaiaaigdaaeaacaWGnbWaaSbaaeaacaaIYaaabeaaaaGaey 4kaSIaaGymaiaac6cacaaIYaWaaOaaaeaacaaI1aaabeaacaGGBbGa ciyyaiaackhacaGGJbGaaiiDaiaacggacaGGUbGaaiikaiaad2eada Wgaaqaaiaaigdaaeqaaiaac+cadaGcaaqaaiaaiwdaaeqaaiaacMca cqGHsislciGGHbGaaiOCaiaacogacaGG0bGaaiyyaiaac6gacaGGOa GaamytamaaBaaabaGaaGOmaaqabaGaai4lamaakaaabaGaaGynaaqa baGaaiykaiaac2facaGG9baaaa@73C5@ (12)

From Equation (11) we can find that the drag force is direct proportional to parameters of D MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamiraa aa@395E@ , C f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4qam aaBaaabaGaamOzaaqabaaaaa@3A69@ , d MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamizaa aa@397E@ and P 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamiuam aaBaaabaGaaGOmaaqabaaaaa@3A47@ . As the acceleration tube length L MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamitaa aa@3966@ has an effect on M 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamytam aaBaaabaGaaGymaaqabaaaaa@3A43@ and M 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamytam aaBaaabaGaaGOmaaqabaaaaa@3A44@ , it is also a significant influencing factor for drag force which is the most important measured criterion of the jet performance of the main nozzle. This paper researched the cotton yarn and the relative parameters are defined as follows: d MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamizaa aa@397E@ =0.3352mm, C f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4qam aaBaaabaGaamOzaaqabaaaaa@3A69@ =0.03, f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOzaa aa@3980@ =0.004, C MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4qaa aa@395D@ =340m/s, R MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOuaa aa@396C@ =287J/Kg K, T 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamivam aaBaaabaGaaGimaaqabaaaaa@3A49@ =296K, C P MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4qam aaBaaabaGaamiuaaqabaaaaa@3A53@ =1004.5J/Kg˙K, and K MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaam4saa aa@3965@ =1.4.

The stress range of weft yarn in main nozzle is divided into four parts as shown in Figure 2. The weft guide passage is the first stress area where the yarn is brought into the main nozzle by negative pressure. The backflow zone is the second section where the drag force is decreased slightly because of the counterblast. The airflow mixing zone is the third section where the yarn is dragged by forward airflow again. The weft acceleration zone is the last section where the yarn is accelerated to the needed velocity with the high-speed airflow. The calculation method of each section is illustrated as follows. The static drag force of weft yarn in previous three sections is obtained by numerical integration with axial length in Equation (2), where the air density ρ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaeqyWdi haaa@3A55@ and velocity V are obtained from numerical simulation. And then in the last section, the drag force is obtained from the Equation (11), where the parameters of both ends of the acceleration tube such as P 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamiuam aaBaaabaGaaGOmaaqabaaaaa@3A47@ , V 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamOvam aaBaaabaGaaGOmaaqabaaaaa@3A4D@ , M 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamytam aaBaaabaGaaGymaaqabaaaaa@3A43@ and M 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGWj0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaqcfaOaamytam aaBaaabaGaaGOmaaqabaaaaa@3A44@ gained by numerical simulation.

From the previous research in consistency between numerical simulation and experiment, the weft drag force calculated by means of numerical simulation is possible and reliable.

Figure 2 Schematic of the weft traction force in stepwise method.

Analysis of structural parameter

In this paper, a combination of numerical simulation and orthogonal test method are used to optimize the key structure of main nozzle such as number of rectifier tank, external diameter of needle tip, distance between needle tip and rectifier tank, α face projected length and diameter of acceleration tube. So it reflected the influence of multiple parameters changed meanwhile on flow field and the most optimal structure is obtained.

Figure 3 shows the main structural parameters of the nozzle for optimization. As the flow passage is formed by nozzle needle and nozzle body, where the airflow is rectified and accelerated to high speed, we select these parameters as research object and make further improvement on jet performance of the main nozzle.

Figure 3 Primary structural parameters of main nozzle.

Orthogonal experiment design

Introduction of orthogonal test method

In order to confirm the effect of main nozzle structure on inner flow field, it is need to carry out thousands of experiment. For example, when there are five factors and four levels on each of them, it is needed to test 1024times.

As the result of the limit of manufacturing and experimental condition, it is impossible to carry out such an experiment; however, orthogonal test method can meet the demand.

Orthogonal test method is one of multiple factors and multiple levels, where only typical experiment is carried out from the full combination method and it has the characteristic of uniform distribution. Each level of factors is distributed more reasonable in experimental program. So we could obtain the most satisfied results by the least consumption. Orthogonal test method has such advantages compared with full combination project: reducing the consumption of manpower and material resource, determining the effect of each factor on measured criterion, and confirming the experiment orientation for future.

The process of orthogonal test method is explained as follows: determining experiment factor and level, selecting appropriate orthogonal table, listing the test scheme and result, analyzing the measured criterion, obtaining the optimal model and verifying result. With the improvement of computational fluid dynamics it is feasible to use software CFD for numerical simulation of flow field and structure optimization. So we adopt Fluent to do numerical experiment of orthogonal program, in order to study the influence of each parameter on nozzle flow field and obtain the optimal structures under different air supply pressure.

Numerical scheme and evaluation criterion

In orthogonal test, the main structure parameters influencing the flow field are selected as research object and others keep the same. In this paper, number of rectifier tank (i.e. ), external diameter of needle tip (i.e. ), distance between needle tip and rectifier tank (i.e. ),  face projected length (i.e. l) and diameter of acceleration tube (i.e. ) are the optimal factors which are in four levels. As shown in Table 1, there are sixteen experiment scheme confirmed by orthogonal table L16 (45). Then the three-dimensional geometry model of the scheme from orthogonal table is established by Pro/E, and the flow condition is illustrated from the result of numerical simulation by software Fluent. In the paper, air velocity of the main nozzle exit and static weft drag force are the evaluation criterion of jet performance and they are obtained by the method above mentioned.

Factor Level

A

B

C

D

E

n

D

d

L

l

1

6

3

2

12

3

2

8

3.5

2.2

14

4

3

10

4

2.4

16

5

4

12

4.5

2.6

18

6

Table 1 Level of each factor in orthogonal test

Numerical result and analysis

The comprehensive balance method is always used to analyze the result of multiple-index orthogonal test. It means that a subtraction between maximum and minimum of measured criterion is done and the range is obtained. From the analyzed range, we have obtained the relationship between each factor and evaluation index. Then the general optimal model is obtained with the comparison between theoretical method and practical experience. In this paper, the structure optimization under the air supply pressure of 0.3MPa is carried out and the analytical method could be also used to other conditions.

Table 2 shows the air velocity of nozzle exit and static weft drag force of the different conditions from numerical simulation. k Express the sum of all experimental indexes under the level i, and is the average. It reflects the affection degree of each factor in evaluation index. The orthogonal test shows the effect of five factors in orders is B>C>A>D>E in Table 1. So it is reasonable to select the air velocity of nozzle exit as an evaluation criterion. At the same time, the static drag force is relative to air velocity and density, so that each factor has the same influence law on the two evaluation index. As shown in Figures 4 & 5, the curves of the nozzle exit velocity and static drag force varies with the factors and levels. As a whole, each factor has the same influence law on the two evaluation index. As the air supply pressure is low, the average of air velocity of the nozzle exit is lower than 300m/s in all levels of each factor. In this paper, the maximum static drag force of cotton yarn is 0.193N, which is close to the experimental value.1 There are some conclusions obtained as follows:

  1. The evaluation index decreases with the increase of rectifier tank number, but the amplitude of variation is smaller than 10%. It is indicated that the factor A has a small influence on flow field of the main nozzle. Considering the process of manufacturing, the air velocity of nozzle exit and weft drag force reach the maximum at the level of 1.
  2. As shown in Table 1, the diameter of acceleration tube has a maximum effect on above two evaluation index. The air velocity and weft drag force is augmented rapidly with the increase of the diameter. When the diameter increases to 4mm, the evaluation index almost becomes invariant, even if the diameter of acceleration tube is infinitely increased; the evaluation is still not invariant except the sharp increase of energy consumption. Considering the energy consumption and weft insertion effect, the factor B is the best choice at the level of 3.
  3. As the throat area decreases with the increase of the external diameter of needle tip, the airflow rate and density is reduced and the air velocity and weft drag force decreases rapidly. So the evaluation index decreases with the increase of factor C which has a great effect on flow field second to factor B. The external diameter of acceleration tube is the best at the level of 2.
  4. With the increase of distance between needle tip and rectifier tank, the air velocity and density decreases firstly and then keeps the same. When  increases, the weft ejection zone moves to exit and the energy consumption becomes larger, so the evaluation index has decreased gradually. To a certain degree, the evaluation index remained stable. It is the best for weft insertion that the factor D is at level of 1.
  5. The projected length range of α face is the least, so it has the least effect on air velocity and weft drag force, where the amplitude of variation is smaller than 10%. As shown in Figure 4, the longer the l with narrow acceleration zone is, the greater the velocity acceleration is. When l is 6mm, the evaluation index is the highest.

The optimal structure combination of the main nozzle is obtained from the range analysis (i.e. A1B3C1D1E4), where the air velocity and drag force reach the maximum. Considering the factor affecting relationship, the assembly of A3B3C1D2E4 in numerical simulation is the closest to the optimal scheme. It is indicated that the analysis method of structure optimization for the main nozzle is correct and reliable from (Table 2) (Table 3). When the throat area is oversize, there is an air leakage phenomenon in weft guide passage, and the ballooning is produced in weft inlet by counterblast. So the weft drag force is decreased and sometimes even lead failure of weft insertion and weft stoppage. In another word, the air leakage phenomenon is not allowed.

From the Table 2, there is an air leakage phenomenon in optimal structure combination of numerical simulation. Taking into account energy consumption and process, expanding the external diameter of nozzle needle (i.e. decreasing throat area) is taken to prevent the leakage phenomenon at weft inlet. With the combination of numerical simulation and range analysis, there is not an air leakage phenomenon in assembly of A4B3C2D4E1 and the evaluation index is close to the optimal scheme of range analysis, so that it is feasible for expanding the external diameter of nozzle needle to 2mm. In the end, the assembly of A1B3C2D1E4 is selected as the optimal structure under the air supply pressure of 0.3MPa. With the verifying of the numerical simulation, the air velocity of nozzle exit of the model reached 303.7m/s (i.e. 0.9M), which satisfied the weft insertion requirement.

Experiment Number

Effect Factor

Outlet Velocity V(m/s)

Leakage

A

B

C

D

E

1

1

1

1

1

1

273.32

N

2

1

2

2

2

2

283.18

N

3

1

3

3

3

3

288.56

N

4

1

4

4

4

4

285.31

N

5

2

1

2

3

4

243.1

N

6

2

2

1

4

3

296.38

N

7

2

3

4

1

2

277.96

N

8

2

4

3

2

1

294.52

N

9

3

1

3

4

2

204.69

N

10

3

2

4

3

1

236.92

N

11

3

3

1

2

4

311.08

Y

12

3

4

2

1

3

303.5

Y

13

4

1

4

2

3

156.41

N

14

4

2

3

1

4

266.89

N

15

4

3

2

4

1

297.47

N

16

4

4

1

3

2

301.11

Y

K1

1130.38

877.53

1181.91

1121.73

1102.25

Sum of the indicator of every effect factor

K2

1111.97

1083.38

1127.3

1045.21

1066.95

K3

1056.24

1175.08

1054.67

1069.7

1044.89

K4

1021.89

1184.49

956.6

1083.86

1106.39

k1

282.59

219.38

295.47

280.43

275.56

Average of the indicator of every effect factor

k2

277.99

270.84

281.82

261.3

266.73

k3

264.06

293.77

263.66

267.42

261.22

k4

255.47

296.12

239.15

270.96

276.59

R

27.12

76.74

56.32

19.12

15.37

Optimal Case

A1

B3

C2

D1

E4

Table 2 Range analysis method of outlet velocity (p=0.3MPa)

Figure 4 Outlet velocity changed with factor (p=0.3MPa)

Figure 5 Weft traction force changed with factor (p=0.3MPa).

Experiment Number

Effect Factor

Static Weft Drag Force F(N)

A

B

C

D

E

1

1

1

1

1

1

0.142

2

1

2

2

2

2

0.15

3

1

3

3

3

3

0.154

4

1

4

4

4

4

0.142

5

2

1

2

3

4

0.109

6

2

2

1

4

3

0.158

7

2

3

4

1

2

0.145

8

2

4

3

2

1

0.16

9

3

1

3

4

2

0.078

10

3

2

4

3

1

0.108

11

3

3

1

2

4

0.183

12

3

4

2

1

3

0.174

13

4

1

4

2

3

0.051

14

4

2

3

1

4

0.137

15

4

3

2

4

1

0.159

16

4

4

1

3

2

0.167

K1

0.588

0.381

0.65

0.599

0.569

Sum of the indicator of every effect factor

K2

0.572

0.553

0.592

0.544

0.541

K3

0.543

0.641

0.529

0.538

0.536

K4

0.514

0.642

0.446

0.536

0.572

k1

0.147

0.095

0.163

0.15

0.142

Average of the indicator of every effect factor

k2

0.143

0.138

0.148

0.136

0.135

k3

0.136

0.16

0.132

0.134

0.134

k4

0.128

0.161

0.112

0.134

0.143

R

0.019

0.065

0.051

0.016

0.008

Optimal Case

A1

B3

C2

D1

E4

Table 3 Range analysis method of weft traction force (p=0.3MPa)

Comparison and conclusion

Figures 6-8 show the comparison of the air velocity, Mach number and static pressure between the optimal model and the original model,1 where the air supply pressure is 0.3MPa. The diameter of acceleration tube, external diameter of needle tip, throat area, and rectifier tank number for the original model are 3mm, 1.9mm, 4.2mm2, and 8, respectively, and are 4mm, 2.2mm, 8.75mm2 and 6, respectively for the optimal model. Other structural parameters keep constant. From these figures, we can find that the air flow velocity and static pressure have the growing tendency with the increase of acceleration tube diameter and nozzle throat area. As shown in Figure 6 & 7, the axial velocity of these two models has little difference in weft guide passage. After being slightly accelerated in cone-shape area, the air velocity remains the same in the straight round pipe and there is not the air leakage phenomenon, which is suitable for weft insertion. As a result of the mixture between the high-pressure flow and atmosphere, the air velocity changed rapidly in weft ejection zone and the backflow zone is produced in two models, where the maximum reaches 0.6M. The length of backflow zone of optimal model is a bit longer than that of the original model, but there is little influence on weft insertion. In acceleration zone, the air velocity of two models is increased gradually, which is accelerated from 273.8m/s of original model to 303.7m/s of optimal model at the nozzle exit. When the air flow is jet into atmosphere, the air velocity is descend rapidly. As the throat area of optimal model is larger than that of original model, the flux is greater than that of the latter and the decrement of air velocity is lower. From above analysis, the air velocity at weft acceleration zone and free jet zone of optimal model is higher than that of the original model, which is benefit for weft insertion.

Figure 6 Axial velocity of optimized model and original model.

Figure 7 Mach number of optimized model and original model.

As shown in Figure 8, the axial static pressure of the two models is stable at an atmosphere in weft guide passage and free jet zone. As a result of the high-speed airflow in nozzle throat, a backflow zone is produced in that area of the two models. The negative pressure of optimal model is lower than that of original model, which is helpful for weft suction from outside. In weft acceleration zone, the air velocity is increased gradually with the decline of static pressure. But the static pressure in acceleration tube of optimal model is higher than that of original model, so that the weft drag force is greater. As for cotton yarn, the weft drag force is enlarged from 0.142N for the original model to 0.192N for the optimized one, and the increasing range is 30%. From above analysis, we can obtain some conclusions on the structure optimization of the main nozzle under the air supply pressure of 0.3MPa:

  1. The axial air velocity and static pressure increases considerably in the assembly of A1B3C2D1E4, and the negative pressure of weft ejection zone are lower, which is benefit for the weft insertion.
  2. The weft drag force is greater in the assembly of A1B3C2D1E4, which is helpful for the increase of weft velocity.

Taking the air velocity, static pressure and weft drag force into account, the exit velocity of the main nozzle in the assembly of A1B3C2D1E4 has reached high subsonic, the weft drag force is greatly improved and there is no air leakage phenomenon in nozzle entrance. The optimal main nozzle is suitable for weft insertion of high-speed air-jet loom and the analytical method could be used for optimization of other types of main nozzles or different air supply pressure.

Figure 8 Static pressure of optimized model and original model.

Acknowledgements

None.

Conflict of interest

Author declares there is no conflict of interest in publishing the article.

References

  1. Ishida M, Okajima A. Flow Characteristics of the Main Nozzle in an Air-Jet Loom Part I: Measuring Flow in the Main Nozzle. Text Res J. 1994;64(1):10‒20.
  2. Ishida M, Okajima A. Flow Characteristics of the Main Nozzle in an Air-Jet Loom Part II: Measuring High Speed Jet Flows from the Main Nozzle and Weft Drag Forces. Text Res J. 1994;64(2):88‒100.
  3. Mohamed MH, Salama M. Mechanics of a Single Nozzle Air-Jet Filling Insertion System Part I: Nozzle Design and Performance. Text Res J. 1986;56(11):683‒690.
  4. Mohamed MH, Salama M. Mechanics of a Single Nozzle Air-Jet Filling Insertion System Part II: Velocity Distribution and Design of the Air Guide System. Text Res J. 1986;56(12):721‒726.
  5. Jeong SY, Kim KH, Choi JH, et al. Design of the main nozzle with different acceleration tube and diameter in an air-jet loom. International Journal of Precision Engineering and Manufacturing. 2005;6(1):23‒30.
  6. Oh TH, Kim SD, Song DJ. A Numerical Analysis of Transonic/Supersonic Flows in the Axisymmetric Main Nozzle of an Air-Jet Loom. Text Res J. 2001;71(9):783‒790.
  7. Dong-xu Y, Yuan Y, Zhi-hua F. Analysis on the Flow Field of the Main Nozzle in an Air-Jet Loom Based on FLUENT. Journal of Soochow University (Engineering Science Edition, in Chinese). 2007;27(2):14‒17.
  8. Guo Jie, Feng Zhihua, Zeng Tingwei. 3D analysis on the inner flow field of the main nozzle in air-jet loom based on FLUENT. Journal of Soochow University (Engineering Science Edition, in Chinese). 2009;29(2):38‒42.
  9. Xiang-Fei L, Zhi-hua F, Zhong-kui S, et al. Calculation and analysis of weft insertion force for air-jet loom’s main nozzle based on Fluent. J Text Res (in Chinese). 2011;32(9):125‒129.
Creative Commons Attribution License

©2017 Lan, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.