The Long Bubbles Penetration through Viscoelastic Fluids with Shear-Thinning Viscosity in a Curved Tube

The purpose of this research is to investigate the behavior of a long bubbles penetrating through viscoelastic fluids in a curved tube. The injection gas flow is controlled by a mass flow controller (MFC). The results of the experiments show that the bubbles width approaches constant value at the location six-diameter upstream from the bubbles front. A difference variable is introduced at the six-diameter location to show the shifting deviation of the bubbles in the curved tube. It is shown that, with the same fluid viscosity and the curved angle, the difference is higher when the gas flow rate is higher. Also, the difference increases proportionally when the capillary number and the Weissenberg number increase.


Introduction
The principle of gas-assisted injection molding (GAIM) is to use high-pressure inert gas injection into the molten plastic forming vacuum melt section and push forward to achieve injection molding, packing, cooling, and other processes.Since the gas has an efficient transfer of pressure, it allows the internal pressure consistent throughout.Thus it can eliminate internal stress and prevent product deformation, while significantly reducing the pressure of the mold cavity.In addition, GAIM also can reduce the product weight, eliminate sink marks, improve production efficiency, improve product design freedom, and so on.GAIM has been extensively adopted to manufacture the components of the automobiles, home appliances, communication products, computer chasses, and other items.
Early scholars Saffman and Taylor [1] use tablet mode to shoot the shape of the bubble when water flows through oil and oil flow through glycerin.Taylor [2] discusses Newtonian fluid when the hollow ratio is changed in the straight tube, taken at a distance of 1.5 times the diameter of the tube from bubble tip; the bubble diameter will tend to constant value.Taylor [2] proposed three possible streamline patterns in the bubble front manifold.At high capillary number (Ca) value, the fluid will occur completely through flow field, which does not show the recirculation zone on free surface.
At lower Ca another possible two flow fields will appear, which are transitional stream and recirculation.Kolb and Cerro [3] published research in long bubble-driven flow pattern in square pipe, which found that fractional converge () will be close to 0.64 at high Ca values.Wong et al. [4,5] experimentally and theoretically analyze that the long bubble in the square pipe drives viscosity fluid.Huzyak and Koelling [6] experimentally investigate the differences between Newtonian fluid and non-Newtonian fluid driven by long bubble condition and pointed out in the text that the relationship between Ca and m is irrelevant with pipe diameter in terms of Newtonian fluid, but that will be affected by pipe diameter in terms of viscoelastic fluid.Paper explores the relationship between Deborah number (De) and .Gauri and Koelling [7] continued its experimental architecture and used particle tracking velocimetry to measure flow field of bubble front; the flow field results agree with two flow fields of full transitional stream and recirculation proposed in Taylor [2].Bi and Zhao [8] use four different cross-sectional shapes of the tube, which is circle, triangular, square, and rectangular tubes; the tube was filled with water and placed vertically and gas was injected in the bottom of the pipe, shooting bubble continuous growth figure of different shapes of tubes by camera.The experimental results showed that it have faster growth of the bubble when the pipe diameter is larger.Yang and Chou [9] investigate the distribution of residual wall thickness of GAIM tube member at diameter junction and bend.Gas passing through the bend cavity will tend to the inside, resulting in nonuniform wall thickness inside and outside and empting the area eccentric phenomena after the gas is blown.Yamamoto et al. [10] analyzed Newtonian and non-Newtonian fluids using experimental methods and investigate the effect of Ca and Wi to .Sousa et al. [11] use the experimental model; the non-Newtonian fluid polyacrylamide (PAA) is injected into the vertical tube by the CCD camera captured images of bubbles, and then the shape of the bubble with a dimensionless manner is presented.Hsu et al. [12] simulated two flow fields mentioned by Taylor [2] using numerical analysis method.Their result also found that there is a conversion flow field between the two flow fields.Ha-Ngoc and Fabre [13] use numerical methods to simulate the bubble shape under different conditions with dimensionless bubble length and width.Hsu et al. [14] investigate the long bubble penetrating the silicone oil in the tube; when the viscosity of the silicone oil is fixed, the gas flow rate is higher and the Ca and Reynolds number (Re) are larger.Cheng-Hsing Hsu et al. [15] investigate the long bubble penetrating volatile organic compounds; when the Re values for acetone, methyl alcohol, and butyl acetate were between 4381∼4430, 2460∼2490, and 1940∼1980, respectively, the transient response phenomenon occurred.Chang et al. [16] experimental results showed that bubble velocity, capillary number, and Weissenberg number increase as shear viscosity increases under a constant gas flow by mass flow controller.However, as shear viscosity increases, bubble diameter decreases and m increases.
From the past related studies, we found that most experiments on long bubbles are focused on the bubble contour and the use of runner of round straight tube to observe the bubbles change of the Newtonian fluid, compared with no change in the development of flow path of other shapes.Therefore, the present research uses a non-Newtonian fluid, with the use of three kinds of curved tube of 90 degrees, 135 degrees, and 180 degrees, as observing the flow variation of double-flow runner.

Experimental Setup.
The architecture of experiment consults that of Hsu et al. [14].Figure 1 shows the experimental setup for the gas penetration.The images of the bubbles are captured by a high-speed camera (CASIO EX-F1) with 6 Mega pixel/60 frames per second in order to observe the change of the flow field involved bubbles and flow type.The camera output pixel is 600 million pixels and the highest animation resolution is 1920 × 1080 pixels.The light source is a halogen spotlight source (LA180-Me) which is characterized by 180 W metal halide lamp and the average brightness of 2.5 million lux.Nitrogen with 99.99% purity is selected as the injection gas.The inlet gas pressure is fixed at 10 kg/cm 2 .The volume rate of injection gas was controlled using a MFC (Brooks 5850) where the maximum output flow rate is 1000 ml/min.The gas was injected at three volume rates 1000 ml/min ( = 100%), 600 ml/min ( = 60%), and 200 ml/min ( = 20%).The tube, which is made of high-hardness heat-resistant glass, is shown in Figure 2. The dimensions of the tube are as follows:  = 800 mm,  = 400 mm,  represent measuring position,  = 50 mm, and  = 8 mm.The curve angle  are 90 ∘ , 135 ∘ and 180 ∘ , respectively.The tube was covered by a viewing box dimensions 50 mm × 50 mm.Glycerin has similar refraction index to that of the glass.A ruler was fixed to the frame for convenience of calculation of the velocity of the bubbles.

Test
Fluids.This work uses viscoelastic fluid.The viscoelastic fluids are 1.5 wt% and 0.5 wt% (carboxymethyl (CMC), polyacrylamide (PAA), and polyacryl sodium (PAS)).The suppliers of the three substances are Chung Yuan Chemicals Incorporation in Taiwan R.O.C. CMC chemical formula is CH 2 CO 2 H, which is a powder-like substance, nontoxic, and easy to dissolve in water, but insoluble in organic solvents and dissolved into a viscous liquid; viscosity will rise due to temperature decline.PAA chemical formula is CH 2 CHCONH 2 , which have low concentration of high viscosity characteristics owing to the intermolecular friction; viscosity decreased with the temperature rise.PAS chemical formula is C 3 H 3 NaO 2 , which is dissolved in water, glycerol, propylene glycol, and other media, and has the temperature stability.With the role of fixed metal ions, it can prevent the negative effects of metal ions on the product and is a variety   of special properties of the surface active agent.These three fluids are water-soluble polymer compounds, characterized by the fluid strain rate and shear stress into a nonlinear relationship.
Rheological measurements were carried out with a coneplate type rheometer (Gemini HR nano).The measurements of rheological properties were carried out at 25 ∘ C, which is the same temperature condition for the gas penetration experiments.All the viscoelastic fluids have shear-thinning viscosity and nonzero first normal stress difference.The shear viscosity of every viscoelastic fluid shows strong shearthinning property.The solutions exhibited shear-thinning behavior as can be seen in the viscosity-shear rate plot of Figure 3.The Carreau-Yasuda viscosity model ( 1) was fitted   to the data and the fittings are represented in the figure by the full lines.Their parameters can be found in Table 1.
where  is the measured viscosity, γ is the applied shear rate,  0 is the viscosity limit when γ → 0, and  ∞ is the viscosity at infinite shear rate, which was set as the solvent viscosity.The parameter  has units of time and  and  are dimensionless parameters described in Sousa et al. [11].Figure 4 shows the shear rate and the first normal stress difference of the test fluids.

Experimental Method.
In this experiment, the gas is injected into tube filled with various viscoelastic fluids in an isothermal environment.The volume rate of the injection gas was controlled by MFC in this experiment.The profile of the bubble is captured by a high-speed camera.Figure 5 shows a photo before the image processing and the ultimate graphics generated by the image processing.The width of the bubble was estimated by these images.
In this work use the image processing software Pho-toCap as the basis of the image after treatment.In this work, we research the impact of the various flow rates and various viscoelastic fluids on the bubble velocity () and the dimensionless parameters, such as the capillary number (Ca), Weissenberg number (Wi), and dimensionless parameter .
In this work, runner channel shape is curved tube and labeled scale in the box.The total lengths of the flow path are different.The length of 90 ∘ curved tube is 230 cm and the shooting range is 0∼220 cm.The length of 135 ∘ curved tube is 253 cm and the shooting range is 0∼240 cm.The length of 180 ∘ curved tube is 293 cm and the shooting range is 0∼ 270 cm.Shooting the whole image of the bubbles has its degree of difficulty, so we shoot each segment taking 10 cm and repeating three times in this research.Passing time can be obtained from every 1/60 of a second difference between the two images.The average of three speed values get an average of instantaneous velocity in that bubble zone, which is defined as velocity of each segment.
The dimensionless parameters are defined as [10] Ca =  ( γ )   , where  is the average velocity of the bubble, ( γ ) is the shear viscosity of the viscoelastic fluids,  is the surface tension of the viscoelastic fluids, γ is shear rate,  is the density of the viscoelastic fluids,  is the Maxwell type relaxation time,  1 is the first normal stress difference, and  is cure tube representative length.We need to determine a representative shear rate γ = / to evaluate ( γ ) and .

Results and Discussion
3.1.Bubble Velocity.In this work, there are long bubble velocity variations within the different curved tube.Three   kinds of curved tube can be divided into four sections for analysis individually.In 90 ∘ curved tube, the first section is 0 cm to 80 cm where the velocity is progressively increased, the second section is 80 cm to 120 cm where the bubbles begin to enter the tube, and the experiment found just from straight tube into the curved tube (80 cm to 90 cm) rate decreased slightly; afterwards the rate began to increase at the location 110 cm to 120 cm about at 45 degree; the third section is 120 cm to 145 cm in which the rate decreased slightly (120 cm to 130 cm) and then the speed increased again until the end of the curved tube (145 cm); the fourth section is 145 cm to 190 cm in which the rate decreased slightly (145 cm to 150 cm) and then progressively increases until export.The variation of bubble velocity in 90 ∘ curved tube is shown in Figures 6-8.
In 135 ∘ curved tube, the first section is 0 cm to 80 cm where the velocity is progressively increased, the second   section is 80 cm to 130 cm where the bubbles begin to enter the tube and then the velocity is progressively increased about at 120 cm to 130 cm section; the third section is 130 cm to 177 cm in which the rate decreased slightly (120 cm to 130 cm) at the bubble through 67.5 degrees and then the speed increased again until the end of the curved tube (177 cm); the fourth section is 177 cm to 220 cm in which the rate decreased slightly (177 cm to 180 cm) as the bubble flow from curved tube end into straight tube end and then progressively increases until export.The variation of bubble velocity in 135 ∘ curved tube is shown in Figures 9-11    In 180 ∘ curved tube, the first section is 0 cm to 80 cm where the velocity is progressively increased; the second section is 80 cm to 150 cm in which the rate decreased slightly (177 cm to 180 cm) as the bubble from straight tube end into curved tube end and then progressively increase until 140 cm to 150 cm where the bubbles begin to enter the tube and then the velocity is progressively increased about 120 cm to 130 cm stage, the third section is 150 cm to 210 cm in which the rate decreased slightly (120 cm∼130 cm) as the bubble through 90 degrees curved tube and then the speed increased again until the end of the curved tube (210 cm), and the fourth section is 210 cm to 270 cm where the rate decreased slightly   in its junction (210 cm∼220 cm) as the bubble flow from curved tube end into straight tube end and then progressively increase until export.The variation of bubble velocity in 180 ∘ curved tube is shown in Figures 12-14.
The change degree of bubble velocities varies by fluid concentrations.At the same curved tube angle and gas flow rate, the higher concentration has the faster velocity in the same fluid; the results and that of Cheng-Hsing Hsu et al. [15] are the same.The reason is that the higher fluid concentration resulting in bubble development is limited such that the leading edge of the bubble is squeezed, so the bubble diameter getting smaller made bubble velocity along faster.
According to theory, however, 80 to 90 cm section of straight tube (curved tube 0 ∘ ) with the speed should be the same with the end of the three curved tube (90 ∘ , 135 ∘ , and 180 ∘ ).However, the present result is not so; this inquiry is because the exit end of the experiment does not have a damping mechanism as a buffer but is directly connected to the atmospheric pressure and the gas in the tube-driven fluid which required less and less resistance decreases case under the cause of bubbles through the apex of curved tube, after the speed has not decreased but increased.Comparison with the velocity change of the long bubble in curve tube, straight tube [15], and suddenly contracting and expanding tube [16]; these results show that when the gas flow rate is fixed, the fluid viscosity is bigger and the bubbles moves faster in the tube.

The Analysis of the Bubble in Equal Diameter Position.
According to past experiments, the bubbles want to investigate the rate of change in the hollow straight tube when the tube is taken at 1.5 times the diameter of the bubbles from the tip, because it rapidly approaches parallel on both sides of the outer edge of the bubble [2].This experiment is to investigate the case of the bubble in the curved tube; the 1.5 times value of the experience in straight tube is not applicable.It must be found within a bubble developing in the curved tube where its width tends to have a constant value.First, the coordinate system of the software PhotoCap is set to Pixel units and then loading image to software, in which case the screen can display the image at any point in the coordinate position.
Image processing sequences of this experiment are as follows.
Step 1. Load original image into software as shown in Figure 15.
Step 2. Detect boundary as shown in Figure 16.Step 3. Adjust contrast as shown in Figure 17.
Step 4. Enlarge image as shown in Figure 18.
It must first confirm the same thickness of up and down the wall to avoid errors when shooting in Step 4. To use the inner diameter converse to get the relationship between the pixel unit and mm, 1 pixel is approximately equal to 0.038 mm.In this work, the analysis of the image assumed that each adjacent pixel of curved tube is parallel.Fixing the -axis values can get the bubbles width for required analysis.Spacing 100 pixel measured the bubble width once.If the bubble in three successive intervals gap within 1 Pixel, they use approximation looking back close to the constant starting value.The distance from starting point to bubble head divided by the diameter of tube is the desire to find multiple.In grayscale judgments, if the boundary grayscale value is less than 100 it is regarded as 0 (black) and more than 100 substituted 255 (white).However, judging boundary grayscale value needs experience to judge.The end of the experiment takes 100 cm-110 cm as grayscale judgments range.Required multiples after processing for each situation are shown in Table 3.Multiple changes from 1.53 times to 5.95 times, so this experiment fixes pipe diameters, which is six multiples from bubble head to analysis difference.
With different curved tube and solution concentration, multiple will change.Most of the cases has the same trends as follows: (1) Under the same solution and concentration, the larger angle of curved tube has the greater multiples.
(2) Under the same angle of curved tube, the larger solution concentration has the greater multiples.
PhotoCap software samples nine kinds of conditions at random, and the percentage error of this experimental method compares with that of MATLAB process (a flow rate of 100%, e.g.) within about 7% (the maximum relative error of 6.86%; the minimum relative error of 0.7%; on average 3.24% relative error).

Difference. The difference formula is seen in Yang and
Chou [9] as follows: where  is difference;  is outer wall thickness;  is inner wall thickness; and  is average thickness.It gets a higher gas pressure which can reduce the extent of gas eccentricity and difference.Equation ( 3) can analyze the offset conditions of gas hollowed zone, and the greater difference has the greater eccentricity.When increasing the gas emptied eccentric resistance has smaller difference.The inner and outer wall thickness in this experiment are defined as shown in Figure 19.In this work, computational fluid difference takes PAS 0.5%, PAS 1.0%, PAS 1.5%, and PAA 1.5% to analysis.The research measures, respectively, the values from the bubble wall inside and outside at the location six-diameter upstream from the bubbles front.To take 15  significant relationship with the fluid flow and viscosity.In the same curved tube and viscosity of the fluid, higher flow has greater difference, and on the contrary lower flow has smaller difference.Investigating the causes, greater flow has faster velocity and greater kinetic energy of gas outward emptied fluid similar to the effect of centrifugal force.Therefore, larger flow has larger difference relatively.Lower velocity when flow is lower has smaller kinetic energy of gas outward emptied fluid.The bubbles difference will be smaller.Cheng-Hsing Hsu et al. [15] discussed fractional ratio of long bubble using silicone oil of different viscosity.Larger viscosity of silicone oil has higher fractional ratio where fractional convergence equals radius of circular tube minus radius of bubble all divided by radius of circular tube.
In this work, relatively the same fluid, but not at the same viscosity when the flow rate is fixed, has the greater eccentricity when the viscosity is larger.The result is the same as Cheng-Hsing Hsu et al. [15].Investigating the causes, greater fluid viscosity limits the development of long bubbles, so the bubble front is smaller which make the gap larger inside and outside and the eccentricity also increases.It was found that maximum difference which is 45 ∘ , 67.5 ∘ , and 90 ∘ is in the middle of angle in the 90 ∘ , 135 ∘ , and 180 ∘ curved tube.The reason is that curved tube curvature becomes large in bubble forming process when flowing from straight tube into curved tube, and bubble front develop into elongated bubble.Outer wall thickness becomes larger, and thus also followed eccentricity becomes larger.In the difference error aspect, the relative mean errors of three experiments of 90 ∘ , 135 ∘ , and 180 ∘ curved tube are 0.38%, 0.53%, and 0.51%, respectively, Angle (degree) and the maximum relative errors are 0.89%, 1.09%, and 1.10%, respectively.

Ca and Wi.
Using PAS 1.5%, PAA 1.5%, and PAS 0.5% fluid and taking 90 ∘ position discuss the relationships of Wi and Ca with difference.Figure 29 shows the case when larger Ca has larger bubble difference.Figure 30 shows the case when larger Wi has larger bubble difference.The relationship between Ca and Wi with difference is proportional.The viscoelastic fluids shear viscosity increases when the bubble speed drops, resulting in increased resistance to gas emptied eccentric; the difference is reduced.The viscoelastic fluids shear viscosity drops when the bubble speed increases, resulting in reduced resistance to gas emptied eccentric; the difference is increased.

Conclusions
In this experiment using viscoelastic fluids and the experimental results obtained at six times the diameter of the bubble of the pipe, the diameter of the bubbles will tend to constant.When the gas flow rate is fixed, the shear viscosity is greater and the relative speed of the bubble is larger.Outlet end directly connected to the atmospheric pressure will push the fluid resistance getting smaller, so that the back end speed of the elbow is not decreased but increased.
In the curved tube angle corresponds to the difference aspect, when the viscosity of the fluid is the same, the flow is higher and its difference will be greater; on the contrary the lower the flow, the smaller the difference.Capillary number   and Weissenberg number are proportional to the difference of bubbles in the tip of the curved tube.

Figure 3 :
Figure 3: The shear viscosity as a function of the shear rate for the studied solutions and Carreau-Yasuda model fitted curves.

Figure 4 :
Figure 4: The first normal stress difference as a function of the shear rate for the studied solutions.

Figure 5 :
Figure 5: Bubbles image before and after processing.

Figure 6 :
Figure 6: The PAS velocity-position diagram of various flow in 90 ∘ curved tube.

Figure 7 :
Figure 7: The PAA velocity-position diagram of various flow in 90 ∘ curved tube.

Figure 8 :
Figure 8: The CMC velocity-position diagram of various flow in 90 ∘ curved tube.

Figure 9 :
Figure 9: The PAS velocity-position diagram of various flow in 135 ∘ curved tube. .

Figure 10 :
Figure 10: The PAA velocity-position diagram of various flow in 135 ∘ curved tube.

Figure 11 :
Figure 11: The CMC velocity-position diagram of various flow in 135 ∘ curved tube.

Figure 12 :
Figure 12: The PAS velocity-position diagram of various flow in 180 ∘ curved tube.

Figure 13 :
Figure 13: The PAA velocity-position diagram of various flow in 180 ∘ curved tube.

Figure 14 :
Figure 14: The CMC velocity-position diagram of various flow in 180 ∘ curved tube.

Figure 15 :
Figure 15: The original image loading software.

Figure 19 :
Figure 19: Difference is defined within the outer gap.

180∘
Figure 21: The PAS 0.5% angle of curved tube: difference diagram of various flows in 135 ∘ curved tube.

Figure 23 :Figure 24 :
Figure 23: The PAS 1.5% angle of curved tube: difference diagram of various flows in 135 ∘ curved tube.

Figure 25 :
Figure 25: The PAS 1.5% angle of curved tube: difference diagram of various flows in 180 ∘ curved tube.
Figure 27: The PAA 1.5% angle of curved tube: difference diagram of various flows in 135 ∘ curved tube.

Figure 29 :
Figure 29: The PAS capillary number: difference diagram of various flows in 180 ∘ curved tube.

Figure 30 :
Figure 30: The PAS Weissenberg number: difference diagram of various flows in 180 ∘ curved tube.

Table 1 :
Parameters of the Carreau-Yasuda viscosity model for the studied solutions.

Table 2 :
Physical properties of viscoelastic fluids.

Table 2
summarizes both surface tensions and the density  of the test fluids at 25 ∘ C. The surface tension was measured with a Du Nouy (DST30) and the density  was measured with (DA-130N).

Table 3 :
The multiples of the diameter from the tip width of the long bubbles approached to a constant value.
Figure 18: Zoom in picture.