0
Research Papers

A Long Wavelength Model for Manufacturing of Continuous Metal Microwires by Thermal Fiber Drawing From a PreformPUBLIC ACCESS

[+] Author and Article Information
Jingzhou Zhao

Department of Mechanical and
Aerospace Engineering,
University of California, Los Angeles,
404 Westwood Plaza,
Los Angeles, CA 90095;
Department of Mechanical Engineering,
Western New England University,
Springfield, MA 01119
e-mails: Jingzhou.zhao@ucla.edu;
Jingzhou.zhao@wne.edu

Xiaochun Li

Department of Mechanical and
Aerospace Engineering,
University of California, Los Angeles,
404 Westwood Plaza,
Los Angeles, CA 90095
e-mail: xcli@seas.ucla.edu

1Corresponding author.

Contributed by the Manufacturing Engineering Division of ASME for publication in the JOURNAL OF MICRO- AND NANO-MANUFACTURING. Manuscript received March 25, 2017; final manuscript received November 1, 2017; published online December 14, 2017. Assoc. Editor: Takashi Matsumura.

J. Micro Nano-Manuf 6(1), 011003 (Dec 14, 2017) (9 pages) Paper No: JMNM-17-1015; doi: 10.1115/1.4038433 History: Received March 25, 2017; Revised November 01, 2017

Abstract

Thermal drawing from a preform recently emerges as a scalable manufacturing method for the high volume production of continuous metal microwires for numerous applications. However, no model can yet satisfactorily provide effective understanding of core diameter and continuity from process parameters and material properties during thermal drawing. In this paper, a long wavelength model is derived to describe the dynamics of a molten metal micro-jet entrained within an immiscible, viscous, nonlinear free surface extensional flow. The model requires numerical data (e.g., drawing force and cladding profile) be measured in real time. Examination of the boundary conditions reveals that the diameter control mechanism is essentially volume conservation. The flow rate of molten metal is controlled upstream while the flow velocity is controlled downstream realized by solidification of the molten metal. The dynamics of the molten metal jet are found to be dominated by interfacial tension, stress in the cladding, and pressure in the molten metal. Taylor's conical fluid interface solution (Taylor, 1966, “Conical Free Surfaces and Fluid Interfaces,” Applied Mechanics, Springer, Berlin, pp. 790–796.) is found to be a special case of this model. A dimensionless capillary number $Ca=2Fa/γA(0)$ is suggested to be used as the indicator for the transition from continuous mode (i.e., viscous stress dominating) to dripping mode (i.e., interfacial tension dominating). Experimental results showed the existence of a critical capillary number $Cacr$, above which continuous metal microwires can be produced, providing the first ever quantitative predictor of the core continuity during preform drawing of metal microwires.

<>

Introduction

Thermal drawing from a preform recently emerges as a scalable manufacturing method for the high volume production of continuous metal microwires, in a way very similar to continuous casting, but at a much smaller length scale. Indefinitely, long metal wires are produced by melting, forming, and then solidifying directly from bulk metal, with the help of a deformable “mold” made of an amorphous cladding material (e.g., usually made of thermoplastic or glass).

During a typical fiber drawing process, as shown in Fig. 1, the metal embedded preform is aligned and fixed under a preform chuck, which is then slowly lowered into a furnace (sometimes an induction coil) which temperature is carefully controlled and stabilized at a designated value. After the heated preform in the furnace necks down under its own weight (or under external pulling forces sometimes), the bottom portion of the preform is cut away and the fiber is fed through a pinch wheel. The produced fiber filament is ultimately winded onto a spool. The diameter of the fiber $Df$ is usually measured by a laser diameter. The draw-down ratio $Dr$ is defined as the ratio between the fiber pulling speed $U(L)$ and preform feeding speed $U(0)$. At steady-state, the law of mass conservation requires $Dr≡U(L)/U(0)=(Dp/Df)2$, where $Dp$ is the preform diameter. The drawing force or fiber tension, which is often monitored by a load cell, is also an important control parameter and measures the tensile stress along the fiber.

Material selection criteria can be found in the literature [1,2] and are summarized here:

1. (1)The viscosity of the most viscous constituent material (i.e., the cladding) should fall between 103.5 and 107 Poise at the drawing temperature for the process to be controllable. Amorphous materials, such as glass and polymers, are typically used as the support (cladding).
2. (2)The melting temperature of the core metal ($Tm$) should be lower than or at least overlap with the drawing temperature ($T$). Low vapor pressure of the metal is desired and its boiling should be avoided.
3. (3)Chemical reactions between the cladding and core materials should be avoided unless intentionally designed, e.g., for in-fiber synthesis purposes [3].
4. (4)It is desired that cladding and core materials exhibit good adhesion/wetting with each other during and after drawing to avoid cracks, bubbles, and fluid instability of the core material(s).
5. (5)The cladding and core materials should have relatively compatible thermal expansion coefficients in the temperature range up to the drawing temperature.

Following the previous criteria, however, does not guarantee the continuity of the metal core. Establishment of process parameters relies mostly on experiences and often involves trial-and-error.

Low melting temperature metals such as tin (Sn), bismuth (Bi), indium (In), and their alloys have been thermally drawn in polymer cladding, such as polyethersulfone (PES), polysulfone, and polyethylenimine with a softening temperature below 300 °C. The resulting metal fibers with rectangular or circular cross sections have critical dimensions ranging from tens to hundreds of micrometers. They are also thermally drawn along with other functional materials (usually semiconductors or conductive polymers) and serve as conductive electrodes in multimaterial functional fibers, which are in turn utilized as one-dimensional photodetectors [4], thermal sensors [5,6], piezoelectric transducers [7,8], chemical sensors [9], and capacitors [10]. The smallest diameter reported for metal wires that can be reliably drawn into infinitely long filaments is about 4 μm and was achieved from Sn0.95Ag0.05 alloy with PES cladding [11]. Beads, discontinuities, and structural deformation would be observed upon further size reduction. Nevertheless, thermally drawn functional fibers embedding indium wires with diameter approaching 1 μm have been demonstrated [12]. Indium arrays were thermally drawn in polymethyl methacrylate cladding to yield a metamaterial fiber that exhibits electric and magnetic responses for terahertz and far-infrared spectrum.

Despite the fact that reliable drawing of indefinitely long amorphous semiconductor and polymer nanowires has been achieved [11], literature survey above suggests that there exists a fundamental size limit to the diameter of thermally drawn crystalline metal wires below which the metal wires become inherently unstable and extremely difficult to control, if not impossible, by current manufacturing techniques. The physical forces that dominate the break-up, however, are still in dispute. Some attribute the break-up to Plateau-Rayleigh instability [13] while others suggest metal cannot withstand shear stress [11]. Resolution of this challenge meets roadblocks both theoretically and experimentally. No model can yet satisfactorily describe the dynamics of the molten metal core during preform fiber drawing, as pointed out by Zhao et al. [14], although numerous works concerning the modeling and simulation of the compound and microstructured fiber drawing processes exist [1518]. Recent attempt made is based on Tomotika's model that was established nearly a century ago [13,19], which considered the break-up of a cylindrical thread in a surrounding flow that is extending uniformly. Quantitative comparisons between simulation and experiments are emerging [20,21], yet far from being applicable for real-time control. Experimental investigation of the core dynamics, on the other hand, is equally challenging, which require high-speed, high-resolution characterization of high aspect-ratio objects embedded in a thick glass/polymer material under high temperature. A simple yet useful model is called for to elucidate the underlying mechanism of diameter control and to predict in real-time the core continuity from measurable process parameters and material properties.

In addition to its implication in manufacturing of continuous metal microwires, the model also leads us to believe that metal core fiber drawing from a preform may provide a new way for study and control of fluid interfaces in addition to the existing flow fields such as entrainment problems [22], flow focusing [23], and microfluidic devices [24]. With the various studies pointing out the possibility of obtaining nanoscale droplets by hydrodynamic forces [14], further advancement of knowledge is hindered by the lack of methods for experimental observation and in-line metrology. The fact that the shape of the metal wires is formed at liquid state while can be observed at solid state provides an unexpected yet new way for the study of the complex nonlinear dynamics of fluid flows at micro/nanoscales.

Problem Formulation

During a typical thermal drawing process, as the preform is slowly fed into the furnace at a constant speed and uniform fibers being pulled from below as shown in Fig. 1, a steady-state in Eulerian frame is soon reached such that the temperature and velocity of the cladding vary along the axis of symmetry but do not change with time, if viewed from a fixed laboratory position. We focus on the dynamics of the molten metal core under the condition that the cladding flow is steady, requiring that the diameter of cladding is much larger than the core $A(z)≫πR2(z,t)$ so that the molten metal core is assumed to have negligible effect on the cladding flow except for the region in the vicinity of the centerline (Fig. 2).

We consider preform drawing of Sn microwire in Pyrex cladding. For the exterior cladding flow, inertia, gravity, and surface tension are assumed to be negligible as compared to the viscous term in the cladding flow following Yarin et al. [25]. For the interior metal flow, the corresponding values of the Reynolds, Weber, and Bond numbers of the molten core under the parameters given in Table 1 are calculated as follows: $Re=(2ρU¯R¯/μ0)=0.4$, $We=(2ρU¯2R¯/γ)=2×10−5$, $Bo=(4Δρga2/γ)=9×10−4$, all of which are much smaller than 1, allowing us to neglect the core inertia and gravity effects. The length of the thread length is the length of the neck-down region measured on the drawn preform from the point where its diameter starts to shrink to the point where the shrinkage ceases.

Effect of the thermal gradient requires further consideration. The viscosity of the cladding is very sensitive to temperature and undergoes orders of magnitude variations along the thread line, and thus cannot be neglected. Viscosity change of the molten metal core due to thermal gradient is neglected since $μ0≪μ(z)$. Thermal gradient also introduces interfacial tension gradient along the thread line, which may result in Marangoni flow that can be approximated by Display Formula

(1)$−dγdTdTdz=μ∂uz∂r$

where $uz$ is the surface velocity due to Marangoni flow in the axial direction, assuming a cylindrical interface. A dimensionless number can therefore be defined to gauge the strength of the Marangoni flow $Ma=−(dγ/dT)(4ΔTR¯/μU¯L)=1.64×10−81.64×10−8$ compared to the cladding flow, which in this case is small, allowing us to neglect the Marangoni flow induced by thermal gradient.

The steady-state cladding flow without core is modeled as a quasi-one-dimensional free surface extensional flow with velocity $(−(r/2)(dU(z)/dz),0,U(z))$. We let molten metal of viscosity $μ0$ be entrained from a nozzle or melt front with radius $a$ by the extensional cladding flow with position dependent viscosity $μ(z)$, and solidifies downstream at $z=L$ with radius $R(L,t)$. Focusing on the limit $μ0≪μ(z)$ following Taylor, Acrivos, Sherwood, and others' works [3032], the coupling between the interior metal flow and the exterior cladding flow simplifies to a balance between the interior pressure and the exterior stress, and the entrainment dynamics can be accurately described by a long-wavelength model provided that the slope of the spout is everywhere small.

Extruded polyethersulfone (PES, BASF Ultrason E3010) rods (extruded and distributed by Port Plastics, Inc.) with an outer diameter of 19.05 mm and length of 10 cm were first dehydrated at 150 °C for 5 days under vacuum (2 Torr) to remove moisture. A through hole was then drilled to allow insertion of a Sn metal wire (SRA Soldering Products, Foxborough, MA). The PES preform was then slowly fed into a vertical tube furnace at 250 °C at a speed of 50 μm/s under a vacuum of 40 mTorr to consolidate the preform and ensure seamless contact between the metal core and the cladding. The consolidated preform was used either for thermal drawing or interfacial energy measurement.

Unsteady volume conservation inside the spout gives Display Formula

(2)$π∂R2(z,t)∂t+∂Q(z,t)∂z=0$

where $Q(z,t)$ is the flow rate of core.

The velocity inside the molten metal core is nearly unidirectional and is composed of a plug flow induced by the cladding flow and a pressure driven flow, which has the form Display Formula

(3)$ucore=u(z,r,t)ez=[U(z)−14μ0∂P0(z,t)∂z(R2(z,t)−r2)]ez$

where $P0(z,t)$ is the interior pressure, and $R(z,t)$ is the radius of the metal core.

At $z=0$, the molten metal comes out from the melt front or a nozzle connected to a large reservoir, so the pressure in the metal is determined by the local stress at the core/cladding interface. We thus obtain an expression for $P0(z,t)$Display Formula

(4)$P0(z,t)=2μ(z)(dU(z)dz+1R(z,t)∂R(z,t)∂t+U(z)R(z,t)∂R(z,t)∂z)+Pclad(r,z)+γκ$
which contains contributions from surface tension due to the full curvature $κ=1/R(z,t)(1+(∂R(z,t)/∂z)2)12−∂2R(z,t)/∂z2/(1+(∂R(z,t)/∂z)2)32≈1/R(z,t)−∂2R(z,t)∂z2$ pressure and viscous stress in the cladding.

Flow rate of core is therefore Display Formula

(5)$Q(z,t)=∫su(z,r,t)dS=2π[U(z)R2(z,t)2−R4(z)16μ0ddzP0(z,t)]=πU(z)R2(z,t)−πR4(z,t)8μ0{∂∂z[3μ(z)(dU(z)dz+23R(z,t)∂R(z,t)∂t+2U(z)3R(z,t)∂R(z,t)∂z)+γκ]+dμ(z)dzdU(z)dz}$

where $u(z,r,t)$ is the core velocity, $S(z,t)$ the cross-sectional area of the core.

The unknown quasi-one-dimensional velocity field can be obtained by volume conservation to give $U(z)=A(0)U(0)/A(z)$, where the cross-sectional area of the cladding $A(z)$ can be easily measured. The unknown position-dependent viscosity can be approximated by $μ(z)=−(A(z)F/3A(0)U(0))(dA(z)/dz)−1$[33], knowing draw force $F=3A(z)μ(z)(dU(z)/dz)$.

Plugging Eq. (5) into Eq. (2) yields a fourth-order nonlinear differential equation to be solved given $A(z)$ and $F$, both being measurable. While previous works focus on the zero flow rate limit $Q(0,t)→0$ in a linear straining flow [22] or pressure driven flow [34,35], this work focuses on $Q(0,t)=πa2U(0)$ in a nonlinear free surface extensional flow with position-dependent viscosity.

Boundary Conditions and Diameter Control Mechanism

At the nozzle opening or melt front under steady-state drawing, the flow rate of molten metal is controlled at $Q(0,t)=Qcore$. Downstream, due to the fact that the solidified metal core has same velocity as the cladding assuming no slip between solidified core and cladding, the radius of the solidified metal $R(L,t)$ is therefore determined by the drawing speed $U(L)$ of the cladding and the flow rate $Q(L,t)$ of molten metal at the solidification front to give $R(L,t)=Q(L,t)/πU(L)$, neglecting the density change of the metal due to solidification. These constitute the boundary conditions $R(0,t)=a$, $Q(0,t)=Qcore,$ and $R(L,t)=Q(L,t)/πU(L)$.

Under a steady-state drawing, metal wire produced being continuous and uniform requires the flow rate of metal core at the solidification front to be independent of time; thus, $Q(L,t)=Q(L)$. Volume conservation inside the metal jet further requires that the amount of molten metal enters at $z=0$ equals to the amount exits at $z=L$, yielding $Q(L)=Q(0)=Qcore$. The third boundary condition $R(L,t)=Q(L,t)/πU(L)$ is thus turned into $R(L)=Qcore/πU(L)$ under steady-state drawing. Recognizing that $Qcore=πa2U(0)$ for preform drawing, and plugging in the third boundary condition to give $R(L)=aU(0)/U(L)=a1/Dr$, we have shown that the diameter control mechanism for continuous metal microwire production by preform drawing is essentially volume conservation and the metal core diameter does follow the draw-down ratio, if a steady-state can be reached.

To validate the proposed diameter control mechanism, steady-state drawings were done under the parameters listed in Table 2 and cross sections of the resultant metal core fibers are shown in Fig. 3. A linear fit of $R(L)$ with $Qcore/πU(L)$ is plotted in Fig. 4. The slope is $0.9945$ with the coefficient of determination equals to $0.9996$, suggesting $R(L)=Qcore/πU(L)$. We can thus conclude that the diameter of the metal produced indeed follow draw down ratio $Dr=(U(L)/U(0))$, provided that the flow rate is controlled at the nozzle or melt front so that the feeding speed of the metal is the same as the feeding speed of the cladding, a diameter control mechanism that potentially holds valid down to nanoscale in principle but never explicitly stated in the literature. It should be noted that the variation in the value of the vertical axis in Fig. 4 can be attributed to the growth of capillary instability on the metal core during drawing. A cylindrical molten metal thread is inherently unstable due to its high surface tension. A sinusoidal wave therefore grows spontaneously giving rise to the variation of the diameter of the metal core along the thread line. Since the magnitude of the growth of the capillary instability is proportional to the base diameter, the variation of the core diameter appears larger for larger cores. In some cases, the growth of capillary instability would eventually disrupt this assumption and result in slopes not equal to one, which lead us to use this as a criteria to test core continuity. In other words, the slope shall always be equal to one as long as the metal core is continuous.

• (i)Isothermal drawing with constant cladding velocity and zero core flow rate ($U(z)=U(0)$, $μ(z)=μ,Q(z,t)=0$)In the simplest case, when $U(z)=U(0)$, $μ(z)=μ$, and $Q(z,t)=0$, Eq. (5) is reduced to Display Formula
(6)$Q(z,t)=πU(0)R2(z,t)−πR4(z)8μ0{ddz[2μ(1R(z,t)∂R(z,t)∂t+U(0)R(z,t)∂R(z,t)∂z)+γR(z,t)−γd2R(z,t)dz2]}$
Rescaling all radial length scales by $a$, axial length scales by $aμ/4μ0,$ pressure by $(γ/a)$, and solve for steady-state solution, we obtain Display Formula
(7)$R*(z*)=−1U*±(1U*)2−42z*+1$
where $U*=(4μμ0U0/γ)$. We have thus shown that our model includes the conical fluid interface solution first derived by Taylor [36] as a special case.
• (ii)Isothermal drawing with exponential cladding velocity
$(U(z)=U(0)(U(L)U(0))zL,μ(z)=μ)$

The cladding profile is available in an exponential form for isothermal drawing of Newtonian fluid [37] with $U(z)=U(0)(U(L)/U(0))z/L$. In Fig. 5, the analytical cladding profile is plotted against the actual cladding profile drawn under the conditions listed in Table 1.

Plugging the analytical profile into Eq. (5) assuming both cladding and core are at steady-state, and neglecting axial pressure gradient in the cladding we obtain the expression for the flow rate Display Formula

(8)$Q=2π{U0DrzLR2(z)2−R4(z)16μ0ddz[2μU0DrzL(lnDrL+1R(z)dR(z)dz)+γ(1R(z)−d2R(z)dz2)]}$

If we scale $z$ by $L$, $R$ by $a$, and $Q$ by $πa2U(0)$, we obtain Eq. (8) in the dimensionless form Display Formula

(9)$Q*U*=U*Drz*R*2(z*)−R*4(z*)ddz*[CaDrz*(1+1lnDr1R*(z*)dR*(z*)dz*)+1R*(z)−1AR2d2R*(z)dz*2]$

where $Ca=2aμU0lnDr/γL$ is the dimensionless strain rate at the nozzle opening; $U*=4λAR2Ca/lnDr$, where $λ=(μ0/μ)$, $AR=(L/a)$;$Q*$ is the dimensionless flow rate of metal with $Q*=1$ meaning metal is flowing at the same average velocity as the cladding at the melt front. This scaling allows the domain of $z*$ be fixed as $[0,1]$ with $R*(0)=1$ and $R*(1)=Q*/Dr$.

Numerical solution of Eq. (9) is obtained using bvp4c function in Matlab with parameters listed in Table 1, under a third boundary condition

$dR*(z*)dz*|z*=1=0$

The shape of the entrance and exit areas matches the theoretical calculations very well. However, the shape of the calculated metal core during steady-state drawing appears similar to the “die swelling” phenomenon due to the viscoelasticity of the polymer melt typically encountered in melt spinning of polymeric fibers. Here, our theory predicts that “die swelling” may also occur during preform drawing of a Newtonian molten metal within a Newtonian cladding as shown in Fig. 6. Yet such swelling behavior is not observed in the measured core profile, possibly because the calculation was based on isothermal cladding profile while the actual profile is obtained under nonisothermal drawing. Multiple tasks are planned as future works to reduce the discrepancy between experimental and theoretical data. The isothermal cladding profile would first be replaced by the actual cladding profile for the calculation of the core profile to take into account the effect of temperature gradient in the axial direction. Other sources of errors would be screened as well. The validity of the quasi-one-dimensional assumption for the cladding flow would be tested by numerical simulations of the two-dimensional cladding flow taking into account both momentum and energy transport. The measurement method would be improved by installing a high-speed camera to the drawing tower to capture the cladding profile in real time to eliminate the errors introduced by solidification of the cladding and core. Measurements and calculations would be conducted on other material pairs such as glass/metal combinations to eliminate the non-Newtonian effect of the polymer cladding. The ultimate goal is to obtain a valid model that can be used for the feedback control of process parameters during the continuous and reliable production of metal core fibers.

Measurement of Interfacial Energy

In the previous calculation, the data of interfacial energy between PES and Sn were needed. However, interfacial energy data for liquid–liquid interfaces between metals and polymers are almost nonexistence in the literature. There are increasing needs [13] for such data for the manufacturing of metal micro/nanowires as well as functional multimaterial fibers. Here, the interfacial energy $γ$ between the molten Sn and PES was measured by sessile drop method using the following equation derived by Birdi et al. [38]: Display Formula

(10)$γ=ΔρgH22(1−cos θ)$

where $Δρ$ is the density difference between the drop and its surrounding, $g$ is the gravitational force of acceleration, $H$ is the limiting drop height, and $θ$ is the contact angle.

To minimize oxidation of the molten metal at high temperature, a dripping method was designed and employed as shown in Fig. 7 to avoid exposure of molten metal to ambient air. A consolidated preform was subjected to furnace heating at 350 °C with an alumina crucible placed below. The lower portion of the preform softened and dripped down into the crucible. A single drop of molten Sn surrounded by PES was thus obtained.

Eight droplets with increasing volume were made, whose drop heights are plotted against their volume to make sure that the height of the largest droplets does reach the limiting height $H$ as shown in Fig. 8. Droplets with heights within 5% of the limiting height are used for the measurement of the interfacial energy.

Measurement of contact angles and drop heights was done on the images of the solidified droplet, which is routinely applied for the estimation of surface tension of molten metals [28]. Figure 9 shows optical images of the Sn droplet after removal of the PES by dissolving in Dichloromethane for contact angle measurement.

Images for three droplets were taken using APPR® B/W digital camera and First Ten Angstroms video software. Drop heights were measured by a caliper, and contact angles were measured by ImageJ [39] with the drop shape analysis plugin developed by Stalder et al. [40].

Measurement results are summarized in Table 3. The interfacial energy between Sn and PES at around the melting point of Sn was determined to be $0.319±0.010$ J/m2, which was listed in Table 1 earlier.

Mode Transition and Critical Capillary Number

We suggest a dimensionless capillary number $Ca=2Fa/γA(0)$ to be used as the indicator of the mode transition from continuous entrainment (i.e., viscous stress dominating) to capillary break-up (interfacial tension dominating). The capillary number is defined as the ratio between the stress in the cladding at the melt front $F/A(0)$ and the Laplace pressure due to the interfacial tension $(γ/2a)$, where $F$ is the drawing force, $A(0)$ is the diameter of the preform, $γ$ is the interfacial tension, and $a$ is the initial radius of the metal core. Interfacial tension leads to capillary break-up while viscous stress in the cladding acts against it. $Ca$ therefore measures the likelihood of the core breaking due to interfacial tension. Large $Ca$ ($Ca>1$) suggests viscous stress is dominating resulting in the growth of capillary instability being suppressed. Small $Ca$ ($Ca<1$) suggests interfacial tension dominating; thus, the metal core has higher chance to disintegrate. A continuous core may start to break as capillary number reduces below a critical value. And a discontinuous core may start to become continuous as the capillary number increases. It can be expected that a critical capillary number $Cacr$ exists, which marks the mode transition. Experiments were conducted to measure the critical capillary number as shown in Fig. 10, and the results were plotted against the aspect ratio, which is the ratio between the cladding diameter and the furnace length. For each data point, a metal core preform is continuously drawn with drawing force measured and recorded as the furnace temperature is slowly reduced. Great care was taken to ensure that the mode transition occurred during each run. The exact location where the mode transition occurred is found under a microscope by examining the drawn fibers. The critical time $tcr$ that the mode transition occurred is determined by dividing the length of the drawn fibers by the drawing speed. The value of the drawing force at $tcr$ can then be determined from the recorded data. The critical capillary number appears to be on the order of $0.5$ in the range of aspect ratio studied.

Optical micrographs of the core continuity under continuous mode drawing and dripping mode drawing are shown in Fig. 11. Below $Cacr$, the metal core breaks due to capillary instability and forms a train of stretched droplets that translate downstream as they are deformed in the cladding flow when $Ca=0.20$ as shown in Fig. 11(b), while Fig. 11(a) shows the optical micrograph of a solidified metal core that is continuous and drawn under a capillary number $Ca=0.78$. The critical capillary number, in this case, falls between 0.2 and 0.78.

Conclusions

We have extended previous authors' theoretical work on viscously entrained jet dynamics and derived a long wavelength model that may be used to solve for the dynamics of the molten metal core during the production of continuous metal microwires by thermal fiber drawing from a preform. Examination of the boundary conditions revealed that the diameter control mechanism for continuous microwire production by preform drawing is essentially volume conservation. The flow rate of molten metal is controlled upstream while the flow velocity is controlled downstream realized by solidification of the molten metal. This mechanism works for the diameter control of nanowire in principle as well. Scaling analysis revealed that the dominant physical forces that govern the dynamics of the metal jet are interfacial tension, stress in the cladding, and pressure in the metal. The accuracy of the model can be improved by replacing the cladding flow with a more accurate profile, either obtained analytically [33], by simulation [41] or measurements [42]. Although derived for metal drawing, semiconductor core is applicable as well, as long as it is molten during drawing. Steady-state solution of the model is compared with experiments and discrepancies attributed to the fact that the process is nonisothermal. A general and simple method to measure the liquid–liquid interfacial energy between molten metal and viscous claddings was designed and implemented for Sn/PES interface. A dimensionless capillary number $Ca=2Fa/γA(0)$ is suggested to be used as the indicator for the transition from continuous mode (i.e., viscous stress dominating) to dripping mode (i.e., interfacial tension dominating). Experiments suggest the existence of a critical capillary number above which continuous metal microwires can be produced, providing the first ever quantitative measure of the core continuity during preform drawing of metal microwires based on process parameters and material properties.

Acknowledgements

The authors are grateful to Yu-sheng Kuo, Shenzhen Xu, Zhengxian Qu, Elaheh Alizadehbirjandi, Maher Lagha, Pirouz Kavehpour, and Gan-ce Dai for insightful suggestions.

Funding Data

• National Science Foundation (Grant No. 1449395).

Nomenclature

• $A(z)$ =

preform diameter at z

• Ca =

capillary number

• Cacr =

critical capillary number

• Df =

fiber diameter

• $Dp$ =

preform diameter

• $Dr$ =

draw down ratio

• F =

drawing force

• H =

limiting droplet height

• L =

• Q =

flow rate

• $R(z,t)$ =

core radius at position z at time t

• $Tg$ =

• $Tm$ =

melting point of metal core

• $Tzz$ =

axial stress

• $U(0)$ =

feeding speed

• $U(L)$ =

pulling speed

• $γ$ =

interfacial energy between core and cladding

• $θ$ =

contact angle

• $μ$ =

• $μ0$ =

viscosity of core

• $ρ$ =

density

Appendices

Appendix: Detailed Derivation

The steady-state cladding flow without core is modeled as a quasi-one-dimensional free surface extensional flow with velocity $(−(r/2)(dU(z)/dz),0,U(z))$, where $U(z)=(A(0)U(0)/A(z))$ by mass conservation and $A(z)$ is the cross-sectional area of the cladding. We let molten metal of viscosity $μ0$ be entrained from a nozzle or melt front with radius $a$ by the extensional cladding flow with position-dependent viscosity $μ(z)$, and solidifies downstream at $z=L$ with radius $R(L,t)$. Focusing on the limit $μ0≪μ(z)$ following Taylor, Acrivos, Sherwood, and others' work [3032], the coupling between the interior flow and the exterior flow simplifies to a balance between the interior pressure and the exterior stress, and the entrainment dynamics can be accurately described by a long-wavelength model provided that the slope of the spout is everywhere small.

Since $μ0≪μ$ and the molten metal core are long and slender, we follow Taylor [31], Zhang [22], and Gordillo et al. [35] in modeling the effect of molten metal core on the cladding flow as a radial flow, $(q(z,t)/2πr)$, corresponding to a line of point sources of unknown intensities situated along $z$-axis. The velocity field in the vicinity of the molten metal core has the form Display Formula

(A1)$(−r2dU(z)dz,0,U(z))+(q(z,t)2πr,0,0)$

where $q(z,t)$ is the strength of the volume flux and has to be determined by the kinematic boundary conditions on the core/cladding interface. The velocity inside the molten metal core is nearly unidirectional and is composed of a plug flow induced by the exterior flow and a pressure-driven flow, which has the form Display Formula

(A2)$ucore=u(z,r,t)ez=[U(z)−14μ0∂P0(z,t)∂z(R2(z,t)−r2)]ez$

where $P0(z,t)$ is the interior pressure, and $R(z,t)$ is the radius of the metal core.

The kinematic boundary condition on the core/cladding interface is $(∂R(z,t)/∂t)+uclad⋅n|r=R(z,t)=0$. Making use of the fact that the core is long and slender, we have at leading order $n|r=R(z,t)=(−1,0,∂R(z,t)/∂z)$. This kinematic boundary condition together with Eq. (A1) yields an expression for the line flux Display Formula

(A3)$q(z,t)2πr|r=R(z,t)=∂R(z,t)∂t+U(z)∂R(z,t)∂z+12R(z,t)dU(z)dz$

The normal stress balance across the spout interface has the form Display Formula

(A4)$σrrext−γκ=σrrint$

where $κ$ is the curvature of the interface, $γ$ is the interfacial energy of the interface, $σrrext$ is the normal stress contribution from flow outside the core, and $σrrint$ is the normal stress contribution from flow inside the core. The exterior normal stress contribution evaluated at the interface is Display Formula

(A5)$σrrext|r=R(z)=−Pclad(r,z)+2μ∂ur(r,z)∂r|r=R(z)=−Pclad(r,z)−2μ(z)(dU(z)dz+1R(z,t)∂R(z,t)∂t+U(z)R(z,t)∂R(z,t)∂z)$

where we have taken the pressure at infinity to be 0. The pressure gradient in the axial direction due to the cladding flow without the line source is

$∂Pclad(r,z)|r=R(z)∂z=μ(z)d2U(z)dz2+2∂μ(z)∂z∂U(z)∂z$
from momentum conservation. The pressure gradient due to the line source is neglected following Gordillo et al. [35] since $∇Psource=0$ in the zero-Reynolds-number limit.

In contrast, the interior normal stress contribution is dominated by the interior pressure $P0$, which can be seen by combining Eq. (A2) with the continuity equation and solve for the deviatoric stress contribution to $σrrint$Display Formula

(A6)$2μ0∂ur(r,z)∂r=2μ0dU(z)dz+12∂2P0(z)∂z2(R(z)2−r2)+R(z)dR(z)dz∂P0(z)∂z$

in which all contributions are small relative to $P0$.

At $z=0$, the molten metal comes out from the melt front or a nozzle connected to a large reservoir, so the pressure in the metal is determined by the local stress at the core/cladding interface. We thus obtain an expression for $P0(z,t)$Display Formula

(A7)$P0(z,t)=2μ(z)(dU(z)dz+1R(z,t)∂R(z,t)∂t+U(z)R(z,t)∂R(z,t)∂z)+Pclad(r,z)+γκ$

which contains contributions from surface tension due to the full curvature

$κ=1R(z,t)(1+(∂R(z,t)∂z)2)12−∂2R(z,t)∂z2(1+(∂R(z,t)∂z)2)32≈1R(z,t)−∂2R(z,t)∂z2$

pressure and viscous stress in the cladding. Unsteady volume conservation inside the spout gives Display Formula

(A8)$π∂R2(z,t)∂t+∂Q(z,t)∂z=0$

where $Q(z,t)$ is the flow rate of core Display Formula

(A9)$Q(z,t)=∫su(z,r,t)dS=2π[U(z)R2(z,t)2−R4(z)16μ0ddzP0(z,t)]=πU(z)R2(z,t)−πR4(z,t)8μ0{∂∂z[3μ(z)(dU(z)dz+23R(z,t)∂R(z,t)∂t+2U(z)3R(z,t)∂R(z,t)∂z)+γκ]+dμ(z)dzdU(z)dz}$

where $u(z,r,t)$ is the core velocity, $S(z,t)$ is the cross-sectional area of the core

$U(z)=A(0)U(0)A(z), μ(z)=−A(z)F3A(0)U(0)(dA(z)dz)−1$

and $F$ the drawing force derived from the tensile stress

$FA(z)=3μ(z)dU(z)dz$

References

Tao, G. , Stolyarov, A. M. , and Abouraddy, A. F. , 2012, “ Multimaterial Fibers,” Int. J. Appl. Glass Sci., 3(4), pp. 349–368.
Ma, Z. , Hong, Y. , Ding, S. , Zhang, M. , Hossain, M. , and Su, M. , 2011, “ Three-Dimensional Micro/Nanomaterials Generated by Fiber-Drawing Nanomanufacturing,” Three-Dimensional Nanoarchitectures, W. Zhou and Z. Wang , eds., Springer, New York, pp. 117–132.
Orf, N. D. , Shapira, O. , Sorin, F. , Danto, S. , Baldo, M. A. , Joannopoulos, J. D. , and Fink, Y. , 2011, “ Fiber Draw Synthesis,” Proc. Natl. Acad. Sci., 108(12), pp. 4743–4747.
Bayindir, M. , Sorin, F. , Abouraddy, A. F. , Viens, J. , Hart, S. D. , Joannopoulos, J. D. , and Fink, Y. , 2004, “ Metal-Insulator-Semiconductor Optoelectronic Fibres,” Nature, 431(7010), pp. 826–829. [PubMed]
Bayindir, M. , Shapira, O. , Saygin-Hinczewski, D. , Viens, J. , Abouraddy, A. F. , Joannopoulos, J. D. , and Fink, Y. , 2005, “ Integrated Fibres for Self-Monitored Optical Transport,” Nat. Mater., 4(11), pp. 820–825.
Bayindir, M. , Abouraddy, A. F. , Arnold, J. , Joannopoulos, J. D. , and Fink, Y. , 2006, “ Thermal-Sensing Fiber Devices by Multimaterial Codrawing,” Adv. Mater., 18(7), pp. 845–849.
Egusa, S. , Wang, Z. , Chocat, N. , Ruff, Z. M. , Stolyarov, A. M. , Shemuly, D. , Sorin, F. , Rakich, P. T. , Joannopoulos, J. D. , and Fink, Y. , 2010, “ Multimaterial Piezoelectric Fibres,” Nat. Mater., 9(8), pp. 643–648. [PubMed]
Chocat, N. , Lestoquoy, G. , Wang, Z. , Rodgers, D. M. , Joannopoulos, J. D. , and Fink, Y. , 2012, “ Piezoelectric Fibers for Conformal Acoustics,” Adv. Mater., 24(39), pp. 5327–5332. [PubMed]
Gumennik, A. , Stolyarov, A. M. , Schell, B. R. , Hou, C. , Lestoquoy, G. , Sorin, F. , McDaniel, W. , Rose, A. , Joannopoulos, J. D. , and Fink, Y. , 2012, “ All-in-Fiber Chemical Sensing,” Adv. Mater., 24(45), pp. 6005–6009. [PubMed]
Lestoquoy, G. , Chocat, N. , Wang, Z. , Joannopoulos, J. D. , and Fink, Y. , 2013, “ Fabrication and Characterization of Thermally Drawn Fiber Capacitors,” Appl. Phys. Lett., 102(15), p. 152908.
Yaman, M. , Khudiyev, T. , Ozgur, E. , Kanik, M. , Aktas, O. , Ozgur, E. O. , Deniz, H. , Korkut, E. , and Bayindir, M. , 2011, “ Arrays of Indefinitely Long Uniform Nanowires and Nanotubes,” Nat. Mater., 10(7), pp. 494–501. [PubMed]
Tuniz, A. , Lwin, R. , Argyros, A. , Fleming, S. C. , and Kuhlmey, B. T. , 2012, “ Fabricating Metamaterials Using the Fiber Drawing Method,” J. Visualized Exp., (68), p. e4299.
Alchalaby, A. , Lwin, R. , Al-Janabi, A. H. , Trimby, P. W. , Fleming, S. C. , Kuhlmey, B. T. , and Argyros, A. , 2016, “ Investigation of Plateau–Rayleigh Instability in Drawn Metal–Polymer Composite Fibers for Metamaterials Fabrication,” J. Lightwave Technol., 34(9), pp. 2198–2205.
Zhao, J. , Javadi, A. , Lin, T.-C. , Hwang, I. , Yang, Y. , Guan, Z. , and Li, X. , 2016, “ Scalable Manufacturing of Metal Nanoparticles by Thermal Fiber Drawing,” ASME J. Micro Nano-Manuf., 4(4), p. 041002.
Suman, B. , and Tandon, P. , 2010, “ Fluid Flow Stability Analysis of Multilayer Fiber Drawing,” Chem. Eng. Sci., 65(20), pp. 5537–5549.
Jasion, G. , Shrimpton, J. , Chen, Y. , Bradley, T. , Richardson, D. , and Poletti, F. , 2015, “ MicroStructure Element Method (MSEM): Viscous Flow Model for the Virtual Draw of Microstructured Optical Fibers,” Opt. Express, 23(1), pp. 312–329. [PubMed]
Pone, E. , Dubois, C. , Gu, N. , Gao, Y. , Dupuis, A. , Boismenu, F. , Lacroix, S. , and Skorobogatiy, M. , 2006, “ Drawing of the Hollow All-Polymer Bragg Fibers,” Opt. Express, 14(13), pp. 5838–5852. [PubMed]
Stokes, Y. M. , Buchak, P. , Crowdy, D. G. , and Ebendorff-Heidepriem, H. , 2014, “ Drawing of Micro-Structured Fibres: Circular and Non-Circular Tubes,” J. Fluid Mech., 755, pp. 176–203.
Tomotika, S. , 1936, “ Breaking Up of a Drop of Viscous Liquid Immersed in Another Viscous Fluid Which Is Extending at a Uniform Rate,” Proc. R. Soc. London. Ser. A, Math. Phys. Sci., 153(879), pp. 302–318.
Xue, S. , Barton, G. W. , Fleming, S. , and Argyros, A. , 2017, “ Analysis of Capillary Instability in Metamaterials Fabrication Using Fibre Drawing Technology,” J. Lightwave Technol., 35(11), pp. 2167–2174.
Xue, S. , Barton, G. , Fleming, S. , and Argyros, A. , “ Heat Transfer Modelling of the Capillary Fibre Drawing Process,” ASME J. Heat Transfer, 139(7), p. 072001.
Zhang, W. W. , 2004, “ Viscous Entrainment From a Nozzle: Singular Liquid Spouts,” Phys. Rev. Lett., 93(18), p. 184502. [PubMed]
Ganan-Calvo, A. M. , Gonzalez-Prieto, R. , Riesco-Chueca, P. , Herrada, M. A. , and Flores-Mosquera, M. , 2007, “ Focusing Capillary Jets Close to the Continuum Limit,” Nat. Phys., 3(10), pp. 737–742.
Christopher, G. F. , and Anna, S. L. , 2007, “ Microfluidic Methods for Generating Continuous Droplet Streams,” J. Phys. D: Appl. Phys., 40(19), p. R319.
Yarin, A. , Gospodinov, P. , Gottlieb, O. , and Graham, M. , 1999, “ Newtonian Glass Fiber Drawing: Chaotic Variation of the Cross-Sectional Radius,” Phys. Fluids, 11(11), pp. 3201–3208.
BASF Corporation, 2017, “ Datasheet From Manufacturer (BASF),” BASF Corporation, Wyandotte, MI, accessed Nov. 20, 2017,
Wang, L. , Wang, Q. , Xian, A. , and Lu, K. , 2003, “ Precise Measurement of the Densities of Liquid Bi, Sn, Pb and Sb,” J. Phys.: Condens. Matter, 15(6), p. 777.
Eustathopoulos, N. , Nicholas, G. , and Drevet, B. , 1999, Wettability at High Temperatures, Elsevier Science, Oxford , UK.
Gancarz, T. , Moser, Z. , Gąsior, W. , Pstruś, J. , and Henein, H. , 2011, “ A Comparison of Surface Tension, Viscosity, and Density of Sn and Sn–Ag Alloys Using Different Measurement Techniques,” Int. J. Thermophys., 32(6), pp. 1210–1233.
Acrivos, A. , and Lo, T. , 1978, “ Deformation and Breakup of a Single Slender Drop in an Extensional Flow,” J. Fluid Mech., 86(4), pp. 641–672.
Taylor, G. , 1934, “ The Formation of Emulsions in Definable Fields of Flow,” Proc. R. Soc. London. Ser. A, 146(858), pp. 501–523.
Sherwood, J. , 1984, “ Tip Streaming From Slender Drops in a Nonlinear Extensional Flow,” J. Fluid Mech., 144, pp. 281–295.
Taroni, M. , Breward, C. , Cummings, L. , and Griffiths, I. , 2013, “ Asymptotic Solutions of Glass Temperature Profiles During Steady Optical Fibre Drawing,” J. Eng. Math., 80(1), pp. 1–20.
Castro-Hernández, E. , Campo-Cortés, F. , and Gordillo, J. M. , 2012, “ Slender-Body Theory for the Generation of Micrometre-Sized Emulsions Through Tip Streaming,” J. Fluid Mech., 698, pp. 423–445.
Gordillo, J. M. , Sevilla, A. , and Campo-Cortés, F. , 2014, “ Global Stability of Stretched Jets: Conditions for the Generation of Monodisperse Micro-Emulsions Using Coflows,” J. Fluid Mech., 738, pp. 335–357.
Taylor, G. , 1966, “ Conical Free Surfaces and Fluid Interfaces,” Applied Mechanics, Springer, Berlin, pp. 790–796.
Middleman, S. , 1977, Fundamentals of Polymer Processing, McGraw-Hill, New York.
Birdi, K. , Vu, D. , and Winter, A. , 1988, “ Interfacial Tension of Liquids From the Height and Contact Angle of a Single Sessile Drop,” Colloid Polym. Sci., 266(9), pp. 849–854.
Schneider, C. A. , Rasband, W. S. , and Eliceiri, K. W. , 2012, “ NIH Image to ImageJ: 25 Years of Image Analysis,” Nat. Methods, 9(7), p. 671. [PubMed]
Stalder, A. , Kulik, G. , Sage, D. , Barbieri, L. , and Hoffmann, P. , 2006, “ A Snake-Based Approach to Accurate Determination of Both Contact Points and Contact Angles,” Colloids Surf. A, 286(1), pp. 92–103.
Reeve, H. , and Mescher, A. , 2003, “ Effect of Unsteady Natural Convection on the Diameter of Drawn Polymer Optical Fiber,” Opt. Express, 11(15), pp. 1770–1779. [PubMed]
Zaporojan, S. , Plotnic, C. , Calmicov, I. , and Larin, V. , 2010, Knowledge-Based Intelligent System Advancements: Systemic and Cybernetic Approaches, IGI Global, Hershey, PA.
View article in PDF format.

References

Tao, G. , Stolyarov, A. M. , and Abouraddy, A. F. , 2012, “ Multimaterial Fibers,” Int. J. Appl. Glass Sci., 3(4), pp. 349–368.
Ma, Z. , Hong, Y. , Ding, S. , Zhang, M. , Hossain, M. , and Su, M. , 2011, “ Three-Dimensional Micro/Nanomaterials Generated by Fiber-Drawing Nanomanufacturing,” Three-Dimensional Nanoarchitectures, W. Zhou and Z. Wang , eds., Springer, New York, pp. 117–132.
Orf, N. D. , Shapira, O. , Sorin, F. , Danto, S. , Baldo, M. A. , Joannopoulos, J. D. , and Fink, Y. , 2011, “ Fiber Draw Synthesis,” Proc. Natl. Acad. Sci., 108(12), pp. 4743–4747.
Bayindir, M. , Sorin, F. , Abouraddy, A. F. , Viens, J. , Hart, S. D. , Joannopoulos, J. D. , and Fink, Y. , 2004, “ Metal-Insulator-Semiconductor Optoelectronic Fibres,” Nature, 431(7010), pp. 826–829. [PubMed]
Bayindir, M. , Shapira, O. , Saygin-Hinczewski, D. , Viens, J. , Abouraddy, A. F. , Joannopoulos, J. D. , and Fink, Y. , 2005, “ Integrated Fibres for Self-Monitored Optical Transport,” Nat. Mater., 4(11), pp. 820–825.
Bayindir, M. , Abouraddy, A. F. , Arnold, J. , Joannopoulos, J. D. , and Fink, Y. , 2006, “ Thermal-Sensing Fiber Devices by Multimaterial Codrawing,” Adv. Mater., 18(7), pp. 845–849.
Egusa, S. , Wang, Z. , Chocat, N. , Ruff, Z. M. , Stolyarov, A. M. , Shemuly, D. , Sorin, F. , Rakich, P. T. , Joannopoulos, J. D. , and Fink, Y. , 2010, “ Multimaterial Piezoelectric Fibres,” Nat. Mater., 9(8), pp. 643–648. [PubMed]
Chocat, N. , Lestoquoy, G. , Wang, Z. , Rodgers, D. M. , Joannopoulos, J. D. , and Fink, Y. , 2012, “ Piezoelectric Fibers for Conformal Acoustics,” Adv. Mater., 24(39), pp. 5327–5332. [PubMed]
Gumennik, A. , Stolyarov, A. M. , Schell, B. R. , Hou, C. , Lestoquoy, G. , Sorin, F. , McDaniel, W. , Rose, A. , Joannopoulos, J. D. , and Fink, Y. , 2012, “ All-in-Fiber Chemical Sensing,” Adv. Mater., 24(45), pp. 6005–6009. [PubMed]
Lestoquoy, G. , Chocat, N. , Wang, Z. , Joannopoulos, J. D. , and Fink, Y. , 2013, “ Fabrication and Characterization of Thermally Drawn Fiber Capacitors,” Appl. Phys. Lett., 102(15), p. 152908.
Yaman, M. , Khudiyev, T. , Ozgur, E. , Kanik, M. , Aktas, O. , Ozgur, E. O. , Deniz, H. , Korkut, E. , and Bayindir, M. , 2011, “ Arrays of Indefinitely Long Uniform Nanowires and Nanotubes,” Nat. Mater., 10(7), pp. 494–501. [PubMed]
Tuniz, A. , Lwin, R. , Argyros, A. , Fleming, S. C. , and Kuhlmey, B. T. , 2012, “ Fabricating Metamaterials Using the Fiber Drawing Method,” J. Visualized Exp., (68), p. e4299.
Alchalaby, A. , Lwin, R. , Al-Janabi, A. H. , Trimby, P. W. , Fleming, S. C. , Kuhlmey, B. T. , and Argyros, A. , 2016, “ Investigation of Plateau–Rayleigh Instability in Drawn Metal–Polymer Composite Fibers for Metamaterials Fabrication,” J. Lightwave Technol., 34(9), pp. 2198–2205.
Zhao, J. , Javadi, A. , Lin, T.-C. , Hwang, I. , Yang, Y. , Guan, Z. , and Li, X. , 2016, “ Scalable Manufacturing of Metal Nanoparticles by Thermal Fiber Drawing,” ASME J. Micro Nano-Manuf., 4(4), p. 041002.
Suman, B. , and Tandon, P. , 2010, “ Fluid Flow Stability Analysis of Multilayer Fiber Drawing,” Chem. Eng. Sci., 65(20), pp. 5537–5549.
Jasion, G. , Shrimpton, J. , Chen, Y. , Bradley, T. , Richardson, D. , and Poletti, F. , 2015, “ MicroStructure Element Method (MSEM): Viscous Flow Model for the Virtual Draw of Microstructured Optical Fibers,” Opt. Express, 23(1), pp. 312–329. [PubMed]
Pone, E. , Dubois, C. , Gu, N. , Gao, Y. , Dupuis, A. , Boismenu, F. , Lacroix, S. , and Skorobogatiy, M. , 2006, “ Drawing of the Hollow All-Polymer Bragg Fibers,” Opt. Express, 14(13), pp. 5838–5852. [PubMed]
Stokes, Y. M. , Buchak, P. , Crowdy, D. G. , and Ebendorff-Heidepriem, H. , 2014, “ Drawing of Micro-Structured Fibres: Circular and Non-Circular Tubes,” J. Fluid Mech., 755, pp. 176–203.
Tomotika, S. , 1936, “ Breaking Up of a Drop of Viscous Liquid Immersed in Another Viscous Fluid Which Is Extending at a Uniform Rate,” Proc. R. Soc. London. Ser. A, Math. Phys. Sci., 153(879), pp. 302–318.
Xue, S. , Barton, G. W. , Fleming, S. , and Argyros, A. , 2017, “ Analysis of Capillary Instability in Metamaterials Fabrication Using Fibre Drawing Technology,” J. Lightwave Technol., 35(11), pp. 2167–2174.
Xue, S. , Barton, G. , Fleming, S. , and Argyros, A. , “ Heat Transfer Modelling of the Capillary Fibre Drawing Process,” ASME J. Heat Transfer, 139(7), p. 072001.
Zhang, W. W. , 2004, “ Viscous Entrainment From a Nozzle: Singular Liquid Spouts,” Phys. Rev. Lett., 93(18), p. 184502. [PubMed]
Ganan-Calvo, A. M. , Gonzalez-Prieto, R. , Riesco-Chueca, P. , Herrada, M. A. , and Flores-Mosquera, M. , 2007, “ Focusing Capillary Jets Close to the Continuum Limit,” Nat. Phys., 3(10), pp. 737–742.
Christopher, G. F. , and Anna, S. L. , 2007, “ Microfluidic Methods for Generating Continuous Droplet Streams,” J. Phys. D: Appl. Phys., 40(19), p. R319.
Yarin, A. , Gospodinov, P. , Gottlieb, O. , and Graham, M. , 1999, “ Newtonian Glass Fiber Drawing: Chaotic Variation of the Cross-Sectional Radius,” Phys. Fluids, 11(11), pp. 3201–3208.
BASF Corporation, 2017, “ Datasheet From Manufacturer (BASF),” BASF Corporation, Wyandotte, MI, accessed Nov. 20, 2017,
Wang, L. , Wang, Q. , Xian, A. , and Lu, K. , 2003, “ Precise Measurement of the Densities of Liquid Bi, Sn, Pb and Sb,” J. Phys.: Condens. Matter, 15(6), p. 777.
Eustathopoulos, N. , Nicholas, G. , and Drevet, B. , 1999, Wettability at High Temperatures, Elsevier Science, Oxford , UK.
Gancarz, T. , Moser, Z. , Gąsior, W. , Pstruś, J. , and Henein, H. , 2011, “ A Comparison of Surface Tension, Viscosity, and Density of Sn and Sn–Ag Alloys Using Different Measurement Techniques,” Int. J. Thermophys., 32(6), pp. 1210–1233.
Acrivos, A. , and Lo, T. , 1978, “ Deformation and Breakup of a Single Slender Drop in an Extensional Flow,” J. Fluid Mech., 86(4), pp. 641–672.
Taylor, G. , 1934, “ The Formation of Emulsions in Definable Fields of Flow,” Proc. R. Soc. London. Ser. A, 146(858), pp. 501–523.
Sherwood, J. , 1984, “ Tip Streaming From Slender Drops in a Nonlinear Extensional Flow,” J. Fluid Mech., 144, pp. 281–295.
Taroni, M. , Breward, C. , Cummings, L. , and Griffiths, I. , 2013, “ Asymptotic Solutions of Glass Temperature Profiles During Steady Optical Fibre Drawing,” J. Eng. Math., 80(1), pp. 1–20.
Castro-Hernández, E. , Campo-Cortés, F. , and Gordillo, J. M. , 2012, “ Slender-Body Theory for the Generation of Micrometre-Sized Emulsions Through Tip Streaming,” J. Fluid Mech., 698, pp. 423–445.
Gordillo, J. M. , Sevilla, A. , and Campo-Cortés, F. , 2014, “ Global Stability of Stretched Jets: Conditions for the Generation of Monodisperse Micro-Emulsions Using Coflows,” J. Fluid Mech., 738, pp. 335–357.
Taylor, G. , 1966, “ Conical Free Surfaces and Fluid Interfaces,” Applied Mechanics, Springer, Berlin, pp. 790–796.
Middleman, S. , 1977, Fundamentals of Polymer Processing, McGraw-Hill, New York.
Birdi, K. , Vu, D. , and Winter, A. , 1988, “ Interfacial Tension of Liquids From the Height and Contact Angle of a Single Sessile Drop,” Colloid Polym. Sci., 266(9), pp. 849–854.
Schneider, C. A. , Rasband, W. S. , and Eliceiri, K. W. , 2012, “ NIH Image to ImageJ: 25 Years of Image Analysis,” Nat. Methods, 9(7), p. 671. [PubMed]
Stalder, A. , Kulik, G. , Sage, D. , Barbieri, L. , and Hoffmann, P. , 2006, “ A Snake-Based Approach to Accurate Determination of Both Contact Points and Contact Angles,” Colloids Surf. A, 286(1), pp. 92–103.
Reeve, H. , and Mescher, A. , 2003, “ Effect of Unsteady Natural Convection on the Diameter of Drawn Polymer Optical Fiber,” Opt. Express, 11(15), pp. 1770–1779. [PubMed]
Zaporojan, S. , Plotnic, C. , Calmicov, I. , and Larin, V. , 2010, Knowledge-Based Intelligent System Advancements: Systemic and Cybernetic Approaches, IGI Global, Hershey, PA.

Figures

Fig. 1

Schematic and process parameters of metal core fiber drawing from a preform

Fig. 2

Image of the neck-down region (a) and the schematic (b) corresponding to an axisymmetric free surface extensional flow in the cladding entraining an immiscible molten metal from a nozzle (or melt front) located at z=0. The entrained molten metal core has radius R(z,t). Downstream at z=L, the metal solidifies with diameter R(L,t).

Fig. 7

Schematic of the procedure to obtain a Sn sessile drop surrounded by PES

Fig. 9

A solidified Sn droplet after removing the PES for contact angle measurements

Fig. 6

Measured and calculated core profile

Fig. 5

Fig. 4

Experimental result supporting the proposed diameter control mechanism

Fig. 3

Cross-sectional images of the metal core fibers obtained from process parameters given in Table 2

Fig. 8

Drop height versus drop volume

Fig. 10

Capillary number versus aspect ratio

Fig. 11

Continuous mode (a) and dripping mode (b)

Fig. 12

An axisymmetric free surface extensional flow (−(r/2)(dU(z)/dz),0,U(z)) in the cladding entrains an immiscible molten metal from a nozzle (or melt front) located at z=0. The entrained molten metal core has radius R(z,t). Downstream at z=L, the metal solidifies with diameter R(L,t).

Tables

Table 1 Material properties and process parameters
Table 2 Process parameters for experimental validation of the diameter control mechanism
Table 3 Interfacial energy measurement

Errata

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections