Abstract

Plastic optical fiber was chosen for information delivery media in smart textile. Cladding layer was peeled off by chemical and mechanical methods to find optimal peeling conditions. Both radial side illumination and longitudinal end-tip illumination were measured for visible light of 627 µm wavelength. A half-cone-shaped jig was manufactured using 3D printing to give various curvature conditions to fibers. Also POFs were embedded in plain weave textile structure to measure the light dissipation effect. The waveguide phenomenon was modeled using discrete ray tracing technique and ray-to-interface collision detection algorithm. Results from the proposed modeling technique showed linear relationship with those from experiment.

1. Introduction

Technology fusion has been an issue for the last decade. Among various approaches in 6T fusion, smart fabric or smart textile seems to be the hottest theme in textile area. The idea of smart textile is to assign additional functionality to garments with IT materials. The new functionality includes not only simple sensing of external environment such as stress/strain [1, 2], temperature [3], or humidity [4], but also active response to the environment [5, 6]. The most common materials used so far for smart textiles were basically conductive materials. They included metals [7], conductive polymers like polyacetylene [8], polyaniline [9], polypyrrole [10], PVDF [11], or semiconductors [12]. This approach used “electrons” as information delivery media.

The above-mentioned electron-based conductive materials have some cons, especially for textile application. Flexibility is the key component to consider for wearable computer. There are ongoing efforts to make a flexible conductive material such as organic light-emitting diodes (OLEDs) [13], graphene [14], graphene-silver nanowire (AgNW) [15], and metal- (Ag-) based nanofiber network [16] or to apply conductive coating like metals like aluminum [17] or single-walled carbon nanotube (SWNT) to nonconductive material surface [18]. Nevertheless, there exists a gap between research and application yet.

We changed the point of view and focused on other type of information delivery media, “light.” Light can deliver multiple number of wavelengths simultaneously, which is the reason why glass optical fiber (GOF) is widely used for telephone and high-speed Internet connection [19]. On the contrary, the high bending stiffness of GOF limits it use in textile application. Instead, plastic optical fiber (POF) seems to be better candidate for our new media, as it can provide fabric structure with both flexibility and light delivery capability. Our goal is to design an active light-emitting garment using POF. The end use can be different depending on the wavelength used. If visible light is used, the POF-textile hybrid structure can be a new kind of fashion item. Otherwise, if ultraviolet light [20] is used, it can be self-sterilized cushions or pillows (Figure 1(a)). Furthermore, light harvesting garment can be thought of (Figure 1(b)). By the way, conventional POF has cylindrical shape and delivers light tip-to-tip direction only. We removed some area of the POF cladding layer by mechanical or chemical treatment to facilitate the textile-thickness-direction light delivery. The light delivery efficiency using the treated POF was modeled theoretically and measured experimentally.

2. Literature Review

Prior art documents were reviewed in both modeling and experimental aspects.

The modeling technique of light transfer in optical fiber, generally known as waveguide theory, was already established in the 1970s. The waveguide theory dealt with mostly rectangular or circular cross-sections [21], which sufficed for GOF or optical devices. And the wavelength of interest was 1550 nm, where the material dispersion is minimized in silica glass [19]. Layers of optical fiber, that is, cladding and core, can have constant (step-index) or nonconstant (gradient-index) level of refractive index [21]. Even if GOF had simple geometry of cylindrical shape, light transfer inside the GOF had to be calculated by solving complex differential equation [22]. Meanwhile, POF has different working range of wavelength, that is, 400∼700 nm [23], and the high curvature makes it difficult to apply boundary condition in the differential equation. Ray tracing can be an alternative technique to predict light transfer if volumetric information of light delivered is only concerned, not detailed waveguide mode information. It assumes that light ray is a discrete series of vectors, and its course can be calculated by tracing the movement of each vector iteratively. Durana et al. used ray tracing to calculate bending loss dependency on POF cladding [24]. Also Lomer et al. used ray tracing to predict lateral polishing of bends in POF-based multipoint liquid-level measurement sensor [25]. But the previous papers still dealt with simple geometry only, and there have been few research on modeling of POF in complex structure like textile weave pattern. For general geometry, Barnhart developed an efficient ray tracing Mathematica® tool kit which can deal with even lenses or graded index fibers [26]. But the Mathematica software is based on exact calculation, and it performance is bounded by hardware specification, so the complex geometry of textile weave structure is difficult to apply.

There have also been numerous works on experimental aspects. Side-glowing optical fiber (SOF), in which intentionally inserted defects make Rayleigh scattering inside the core, was used in clinical [27], or photocatalytic [28] applications. It was widely commercialized, but the con is that the direction and position of defects cannot be controlled by end users. Meanwhile, slanted fiber Bragg grating (sFBG) is more elegant way of fabric-thickness-direction light delivery. In sFBG, gratings with constant intervals tilts incident light out of the fiber longitudinal axis [29]. Zhao et al. [30]and Jáuregui et al. [31] used sFBG for temperature sensor and interrogation unit, respectively. sFBG can react to very specific range of wavelength, and the tilting angle can be controlled by adjusting the grating inclination. But the high price of grating imprinting device and process is bottleneck for commercialization. As textile products have generally low range of customer price compared with that of high-tech IT devices, we need a method with reasonable manufacturing price, yet maintaining efficiency and controllability in light delivery. Therefore, this paper used an alternative approach that claddings of POF were peeled off locally and light was delivered through the damaged region. In this way, the direction of light leakage could be easily controlled, and the amount of light delivered could be modeled efficiently with computational techniques.

3. Experimental

3.1. POF Specimen

Commercially available PMMA-based step-index POF (Toray, Raytela®) was chosen for experiment (Tables 1 and 2). Fiber diameters were chosen considering the efficiency of weaving process afterwards. Light intensity was measured before and after cladding removal process.

3.2. Cladding Removal

To find optimal peeling off condition, both chemical and mechanical methods were applied. For chemical approach, dimethyl ketone (DMK or acetone; from Daejung chemical corp.) and dimethyl formamide (DMF; Daejung chemical corp.) were prepared in 100% solution. 2 cm length of POFs were immersed in the solution for 0.5/1/2/3 hours and then water-washed. For mechanical cladding removal treatment, sandpaper (grit size P80; ISO/FEPA) was used.

To verify the cladding-removal effect, fibers were impregnated in epoxy resin, and then their cross-sections were observed using cryo-ultramicrotomer (RMC/PTPC&CRX) and field emission scanning electron microscope (JSM-6500F, JEOL). Conventional refractive index profiler (IFA-100; Interfiber Analysis, LLC) was also tested, but the device did not operate for fibers of irregular diameter of ours.

3.3. Light Intensity Measurement

Red light LED (627 nm; Ocean Optics, Inc.) was used for light source, considering the fact that our POF specimen shows lowest attenuation under 650 nm. Two sets of USB-type flame spectrometer (FLAME-S-VIS-NIR, Ocean Optics corp.) were prepared to measure both the side illumination and end-tip illumination. POF, the light source, and the spectrometers were interconnected using bare fiber adapter (BFA-KIT; Ocean Optics, Inc.) to minimize connection loss. Note that the intensities in the following graphs would be shown in relative normalized value, because the light strength could only be controlled manually.

3.4. Half-Cone Jig Test

To measure light delivery under exact curvature condition, a half-cone-shaped jig was designed and manufactured via 3D CAD tool (123D®, Autodesk, Inc.; Figure 2(a)) and 3D printer (Cube Trio Pro®, 3DSystemsn Inc.; Figure 2(b)). The cone height and bottom radius were designed as 200 mm and 125 mm, respectively, but manufactured as 198 mm and 124 mm due to 3D printer filament’s thermal shrinkage. Side-illumination, which came from the cladding-removed region, and end-tip illumination, which came from the fiber end, were simultaneously measured using the same light source and two spectrometers (Figure 3). The top surface of the exposed cladding-removed region was covered by additional 3D-printed cover, which was developed in our previous work [32].

3.5. POF-Textile Weaving

A small hand-loom (Loomini®; Luziuum corp.) of 22 cm × 32.5 cm size was used to make a plain weave fabric in which the POF should be laid. Plain weave means a fabric structure where warp and weft yarns are crossed one by one iteratively like checkerboard [33]. POF and wool/nylon combined yarn (70% : 30%; diameter 3 mm) were used as warp thread, while cotton/acrylic combined yarn (60% : 40%; diameter 2 mm) was used as weft thread. Figure 4 shows the plain weave textile with three kinds of POFs (PGR-FB250, PGR-FB750, and PGR-FB1000) inserted.

4. Modeling

4.1. Assumptions

As our concern is only to find the quantity of light delivered through POF in the fabric-thickness-direction, the following assumptions were used for the simplicity’s sake:(i)Ray is considered as a kind of vector. It behaves like a particle in marching or total internal reflection, and it changes its course in refraction.(ii)Each layer of POF is homogeneous material. Ray-molecule interaction like Rayleigh scattering is not considered.(iii)Rays does not intervene with each other. Harmonic mixing like four-wave mixing does not occur.(iv)The incident light is dissipated as either side illumination or end-tip illumination. Any other geometric loss is ignored.(v)In simulation, POF’s one end is positioned at origin and extends to +x axis direction. Rays start from the first center point of POF.

4.2. Algorithm Procedure

The simulation proceeds in the following order:(i)(Step #0) preparation step: ray vector xω,τ with wavelength ω at time τ is generated. The number of xω,τ follows Gaussian distribution. Orientation of xω,τ are given some deviation, whose maximal value lies inside the acceptance angle of the POF. POF mesh model with number of cylinder center Ncs is prepared. Core and clad cylinders are assigned to each center.(ii)(Step #1) marching step: for each ray vector xω,τ, update its position by time step Δt and get new vector xω,τ + 1. Increase number of iteration by 1.(iii)(Step #2) collision detection step: to check collisions between ray and clad-core interface, search the nearest center point ci from xω,τ + 1. Check the collision in cylinder [cici + 1].(iv)(Step #3) update step: Accept xω,τ + 1 if no collision occurred. Otherwise, modify xω,τ + 1 (see Section 4.6). Display the vectors and prepare the wavelength distribution chart. If the number of iteration exceeds some criteria, stop. Otherwise, go back to step #1.

4.3. Fiber Model

Hollow-cylinder-shaped fiber model was proposed to represent POF with arbitrary curvature and thickness. Each fiber was considered to be a collection of cylinder center point ci’s (Figure 5(a)). Each center point pair, that is, ci and ci + 1, comprised a repeat unit. A hollow-shaped cylinder is assigned one by one for each layer of POF, that is, core and cladding, with different radii (Figure 5(b)). These radii, that is, Rcore and Rclad, would be used for later 3D rendering and ray collision detection. Shape or curvature of the fiber could be controlled by moving control points ci. Figure 5(c) shows simulation model for half-cone test. Each center point was adjusted in its position using cantilever beam analogy. Equations (1) and (2) show conditions needed in finding the deflected new data positions ci′ for a given Rc value.

Figure 6 shows the resultant artificial POFs with varying curvature of radii values using the hollow-cylinder structure.

4.4. Cladding-Removal Simulation

Cladding-removal was simulated by adjusting the radii of cylinders in Figure 5(b). When cladding-removal mode was turned on, the user-selected cylinder segment was heighted (Figure 7(a)). If a user inputted the new diameter of core and cladding, the resultant fiber was updated and displayed (Figure 7(b)). Repeating this procedure, user could remove POF cladding as much as needed (Figures 7(c) and 7(d)).

4.5. Window-Based Simulation Software

The proposed simulation technique was implemented with C++ language (C++ Builder XE5; Embarcadero Technologies, Inc.) under MS Window 7 64 bit O/S. OpenGL 2.0 shading language was used for the 3D display of objects. Micrometer was used as default unit for fiber dimension and ray wavelength representation. Figure 8 shows the screenshot of the simulation software running.

4.6. Ray-to-Interface Collision Detection

The ray vectors had been marched step by step, until they met core-cladding or cladding-air interface. Their final positions were modified based on the Snell’s law [21]. Collision point pc can be found from a simple bilinear interpolation parameter tc,

and p1 denotes a vector end point when ncore = nclad. Vector after refraction v1c,rfr can be known simply by

And then vector end point after total inflection p1,tif becomes

Figure 9 shows the scheme of ray-to-interface collision detection. The detailed description on the parameters is listed in Table 3.

The ray vector was assigned a status flag “INITIAL” when they started to march, “STOPPED” when they go out of the fiber, and “REFLECTED” and “REFRACTED” when they make reflection or refraction, respectively. Each number of flags was counted to make a final artificial photometer graph (Figure 10).

Figures 11 and 12 show an example of ray marching simulation for a virtual half-cone jig test (Figure 11) and a fiber within plain weave structure (Figure 12). The plain weave’s crimp curve was generated by applying sinusoidal function to the center point point ci.

5. Results and Discussion

5.1. Cladding Removal

Several famous solvents, such as benzene, cyclohexane, sodium hydroxide powder, ethanol, hydrogen chloride, were tried as cladding removal agent. But they were not suitable because some did not erode the cladding at all, or others reacted with PMMA, leaving very dirty surface. Only DMF and acetone showed satisfactory result.

Figures 13 and 14 show the cladding removal results using chemical (DMF/acetone) agent for PGR-FB250 specimen. Diameters measured after treatment are displayed in Figure 15 and Table 4. DMF was the most advantageous in the removing speed (Figure 15), and also the cross-sectional boundary after treatment was also very clean throughout the various immersion times (Figure 13). The clean boundary means that dissipated light can have equal radial distribution at side-illumination. However, fast removal time can sometimes be a disadvantage in process control aspect. Our specimens have so small cladding thickness (10 µm in the case of PGR-FB250) that excessive immersion with DMF can unnecessarily damage the core layer. Moreover, DMF is toxic (health danger index 2; NFPA 704) [34, 35] and volatile organic agent, which is not so environment friendly. So acetone was chosen alternatively for our cladding removal agent. It shows about half speed of DMF, and the plotted point in the graph deviates more from the regression line compared with that of DMF. However, it is less toxic (health danger index 1; NFPA 704), and its lower abrasion speed (Figure 15) can be more advantageous in controlling the operation speed.

Also, mechanical peeling off agent was tested. The diameter decreased easily with a few number of abrasions only (Figure 15). The treated boundary was clean (Figure 16), but the fiber dimension was not regular with respect to the radial direction (especially Figure 16(d)). Figure 17 shows the radial distribution of side-illumination for both acetone and sandpaper treatment, respectively. Sandpaper-based method shows less evenness in the graph. Generally, textiles have relatively small thickness compared to area, and the side-illuminated light is expected to be dissipated to the out-of-plane direction. In this case, mechanical abrasion can control the direction of the treated surface, and this radial irregularity can be helpful. However, acetone-based chemical treatment was chosen to guarantee the evenness of abrasion quality in this paper.

5.2. Comparison of Side and End-Tip Measurement

The purpose of half-cone test in Figure 3 was to find proper illumination spot for our specially designed POF. Conventional POF is given light from one end and measured the intensity at the other. Our fiber had another exit, which was the cladding-removed region, and the accuracy of side-illumination measurement method was questioned. So we measured both side illumination and end-tip illumination under the same light source. Figure 18 shows the results for the three kinds of POF specimens, respectively. Theoretically, areal sum of side illumination and end-tip illumination should be equal to that of undamaged fiber. But Figure 18 showed that there was 2∼3% error. This may be due to the connection error between a fiber and the cover structure. The error norm was constant level through the tests, but not negligibly small. So it was natural to choose end-tip measurement type in this paper, where the number of clad-removed region is only one for each fiber. However, if the number of clad-removed sites is multiple, it is worthy to consider side-measurement type instead.

5.3. Effect of Curvature to Illumination

Figures 1921 show the experimental results for half-cone jig tests. The radius of curvature ranged from 0 to 43.2 mm. Firstly, Figure 19 is the end-tip illumination experiment result for PGR-FB250 POF after 3-minute acetone immersion treatment. Note that curve from the undamaged fiber of Figure 18(a) was added to the top most position (thick solid line) as reference. The increase of radius of curvature showed linear relationship with the decrease of end-tip illumination. It seems that global bending within the fiber increased the incidence angle θ0 of (1), leading to more number of reflected light vectors outside the fiber core.

Figure 20 is the result at the same condition for PGR-FB750 fiber. The linear tendency is the same, while wavelength above 635 nm shows more loss. The phenomenon becomes more vivid in Figure 21, which is for PGR-FB1000 fiber. Wavelength under 635 nm was also decreased in this case. The major difference between PGR-FB2500, PGR-FB750, and PGR-FB1000 was the radii. The thicker fiber would be given more bending moment under same radius of moment.

The experimental data were compared with those from simulation. Figure 22 illustrates the integral sum of the graph area of Figures 1921, with simulation results added. The simulation showed highly linear relationship between radius of curvature and end-tip illumination. There was slight difference in the inclination between the simulation and the experiment. The error source could be that radius of curvature was not properly controlled in the simulation because our model uses discrete series of cylinders instead of continuous ones.

Table 5 lists the simulation time for virtual half-cone test of Figure 11. The time was linearly dependent on the number of rays and cylinder center points. Existence of POF curvature delayed calculation, but the degree of curvature did not show linear relationship with the calculation time. The reason is thought to be that collision detection is checked for every ray vector in each frame. An efficient skipping algorithm such as temporal coherence would accelerate the calculation. This means that our ray-tracing-based approach is a suitable method for textile-based POF ray delivery modeling quantitatively and qualitatively.

5.4. Illumination within Weave Structure

Finally, the clad-removed fibers of three different radii were impregnated within a plain-weave textile structure after 3-minute acetone treatment. Plain weave structure was chosen because fibers were given more crimp deformation than in twill or satin. Figure 23 shows the result. It reveals that the thicker fiber showed more loss in the end-tip illumination, which means more lights were dissipated in the clad-removed region. Bigger radius means that the more possibility lies that ray vector can travel through the fiber with larger incident angle out of total internal reflection criteria. The simulation showed also similar relationship. It is also our further work to test other weave types, that is, satin, twill, or arbitrary structure, for side illumination effect of clad-removed POF.

6. Conclusions

Wearable computer needs two necessary components: media that contains information and material that transfer the media. In this paper, light was selected as an alternative media for information delivery, instead of electricity. And as material candidate, POF was chosen for our research because it was the most flexible enough to withstand severe curvature of textile. To find the optimal condition for curved POF light delivery, both experimental and theoretical approaches were taken.

For experiment, the fiber cladding was partially removed by chemical or mechanical treatment for better fiber radial direction dissipation or, namely, side illumination. DMF showed the better peeling effect, but acetone was chosen because of DMF’s toxicity. Optimal clad-removing condition was found that about 2-minute immersion in acetone was enough to remove cladding layer of 250 µm diameter POF. Commercially available POF of three different sizes was used as specimen, and a specially designed half-cone-shaped jig was used for reliable curvature condition. The side illumination could be used to measure radial-direction light dissipation directly, but it showed about 2% error. Alternatively, end-tip illumination was measured throughout the paper, and its amount was in linear relationship with the radius of curvature. A simple plain weave textile structure with POF impregnated was generated with a hand-loom, and it showed that thicker fiber was given more bending rigidity, which leaded to more side illumination of POF.

Also, the proposed method was modeled theoretically. Window-based ray-tracing software was developed to deal with complex textile weave structure geometry. Light ray was considered to be a series of discrete vectors, and the collision between ray vector and material interface was calculated. The half-cone test and weave structure were also simulated in the software. There was a little error between simulation and experimental values because material dispersion was not considered in the simulation. However, the simulated result showed the same curvature-illumination linear relationship like the experimental one. It is our further work to extend the modeling to more complex weave structures, such as twill, satin, and arbitrary Jacquard weaves.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (Ministry of Science, ICT and Future Planning) (no. NRF-2014R1A1A1002463). Also, this study was carried out with the support of the Global Expertise Development Project (10048412) of the Ministry of Trade, Industry, and Energy in 2016.