Microsoft Word - Vol 12_no1_2011_June.doc IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 53 DESIGN AND FABRICATION OF MICRONOZZLES K.H. CHEAH AND J.K. CHIN Department of Chemical and Environmental Engineering, University of Notitngham - Malaysia Campus, Jalan Broga, 43500 Semenyih, Selangor Darul Ehsan, Malaysia. E-mails: Jit-Kai.Chin@nottingham.edu.my, Cheah.Keanhow@nottingham.edu.my ABSTRACT: Micronozzle, a key component in micropropulsion system, has been designed and fabricated. Quasi 1D inviscid theory was used in designing a series of conical micronozzles of different expander half-angles (10°-50°). Aerospike micronozzle, a promising candidate to achieve high performance propulsion system, was designed with Angelino method (or Approximate method). Both micronozzles were fabricated using soft lithography, an inexpensive and relatively simple technique comparing to well-established deep reactive ion etching (DRIE) technique, with polydimethylsiloxane (PDMS) as structural material. Micronozzles with two different nozzle throat width, 53.5 µm and 107 µm, were fabricated for comparison. Microscopic inspections reveal 107µm is the more producible nozzle throat width with current equipments. The PDMS-based micronozzle can be used as cold gas microthruster system for micro- and nanosatellites. ABSTRAK: Mikronozzle, satu komponen kunci dalam sistem pendorongan-mikro telah direkabentuk and dibina. Teori quasi 1D inviscid telah digunakan dalam perekaan beberapa siri conical micronozzles dengan expander sudut-separa yang berbeza (10°- 50°). Aerospike micronozzle, calon yang menjanjikan pencapaian prestasi tinggi untuk sistem dorongan telah direkabentuk dengan menggunakan Angelino method (atau Approximate method). Kedua-dua micronozzles tersebut telah dibuat dengan menggunakan soft lithography, suatu teknik yang relatif murah and senang, dengan polydimethylsiloxane (PDMS) sebagai bahan struktur. Micronozzles dengan dua dimensi nozzle throat yang berbeza, 53.5µm and 107µm, telah dibuat untuk perbandingan. Pemeriksaan mikroskopik mendedahkan bahawa 107µm adalah dimensi yang boleh dibuat dengan peralataan yang sedia ada. KEY WORDS: Micronozzles, Micropropulsion, Microfabrication, Soft lithography 1. INTRODUCTION Microsatellites (<20 kg) have provided the space community a low cost platform in performing space education and demonstration of new technologies in the last decade. The launching cost reduces significantly due to their small size and volume. Evolution of these satellites to perform actual scientific missions requires more pronouns manoeuvrability and fast response attitude control system for precise pointing and formation flying [1]. In order to fulfil the requirements of future missions, micropropulsion system with low thrust (µN to mN), system mass (<1 kg), volume (<1 cm3) and power (<1 W) is required [2]. IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 54 Previous numerical studies show that performance of micronozzle, a key component in micropropulsion system, is strongly affected by viscous effect [3, 4]. A deep micronozzle is preferred to mitigate viscous loss due to merging of viscous boundary layers at end walls. Aerospike micronozzle with virtual wall boundaries could further reduce viscous loss [5]. Deep reactive ion etching (DRIE), a proven but complicated technique in fabrication of micronozzles, would cause photomask breakdown over time and reduced wall straightness as the depth increases [6]. Moreover, the technique requires expensive equipments and tooling. In contrast, soft lithography, an inexpensive and relatively simple microstructures fabrication technique, enables fabrication of ultra-thick structure as thick as 1.5 mm with good vertical wall profile [7, 8]. In the present study, we first design two types of micronozzle, conical and aerospike, with the aim of improved efficiency through deep micronozzle approach. Micronozzles are then fabricated using soft lithography with PDMS as structural material. Surface contour and geometry accuracy of fabricated micronozzles are examined under inverted light microscope and scanning electron microscope (SEM). 2. DESIGN OF MICRONOZZLES Designing micronozzles requires some unique considerations. Flow characteristics within micronozzles are different from their macro-scale counterparts. The flow is characterized by low Reynolds number, Re, where  VL Re or `Re A Lm    (1) For a typical micronozzle operates at mass flow rate in the order of µ kg/s and characteristic length in µm scale, Re is within the range of 101 to 103 suggesting the significance of viscous effect and laminar flow behaviour [9]. Macro-nozzles experience small loss (less than 1% with occasionally up to 2%) of total energy due to viscous effect [10]. Viscous boundary layer can occupy substantial portion of expander section causing a severe performance degradation in micronozzle. Traditional nozzle design approach based on inviscid flow theory is quite satisfactory in designing macro-nozzle but it is not totally applicable when designing micronozzles. However, inclusion of viscous effect would complicate the governing equations. In light of the complication, conventional quasi 1D inviscid theory was used as first approximation in current design. Important findings (e.g. larger micronozzle depth, shorter expander length etc.) from previous numerical studies on micronozzles were adapted into current design to mitigate viscous loss [3, 4]. Fabrication constraints inherited from available equipments were considered and relevant design limits were applied during the design process. Nozzle throat width, wt, was limited to not smaller than 50 µm. Micronozzle depth was limited to 200 µm after compromising between performance degradation due to 3D viscous effect and allowable nozzle throat width [4]. The wt was magnified by a factor of two in order to evaluate the fabrication capability while the nozzle depth was reduced accordingly to maintain a constant area of nozzle throat. IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 55 2.1 Conical Micronozzles Conical nozzles being the simplest nozzle shape are commonly used in micropropulsion system. Conical micronozzle with optimum expander half-angle demonstrates very similar performance in comparing with bell-shaped micronozzle [12]. Quasi 1D analysis [10] based on isentropic flow equations and rocket nozzle theory [13] were combined to evaluate for micronozzle design parameters such as micronozzle geometry and its internal flow characteristics. One of the most important parameters, nozzle expansion ratio, is formulated as 1 1 2 2 1 1 1 2 2 1 2                           k k e M k k e MtA e A (2) The design parameters are summarized in Table 1. Conical micronozzles with expander half-angles, θ, ranging from 10º – 50º were designed. Details of the design process can be found in authors’ previous work [11]. Table 1: Micronozzle design parameters. Design Input Value Design Output Value Chamber pressure 280 kPa Expansion ratio 9.8 Chamber temperature 513 K Nozzle throat width 53.5 µm (107 µm) Nozzle depth 200 µm (100 µm) Nozzle exit width 522.7 µm Nozzle exit Mach no. 3.65 Mass flow rate 4.539 x 10-6 kg/s Specific impulse 106.2878 s Throat Reynolds no. 1207 2.2 Aerospike Micronozzles Aerospike nozzles were first investigated in early 1960s. The nozzle is characterized by no outer walls allowing exhaust gas to expand freely without constraining from walls. It is a promising candidate in achieving high performance propulsion system owing to this unique altitude compensation feature. Current aerospike micronozzles were designed with Angelino method (also known as Approximate method) [14]. The method combines Prandtl-Mayer expansion theory and area-Mach relation (equation 2) in determining spike geometry. Figure 1 shows a typical 2D aerospike nozzle. IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 56 Fig. 1: Two dimensional aerospike nozzle [11]. Characteristic lines cross the nozzle boundary at distance l from the lip, can be computed by   tt AAll /1sin//  or M  (3) where α = µ(M) – ν(M). Based on the same design parameters as conical micronozzles, length of the lines were calculated with equation (3) and converted into X and Y coordinates. Corresponding spike contour was plotted in Fig. 2. Fig. 2: Spike contour calculated with Angelino method. 3. FABRICATION OF MICRONOZZLES 3.1 Materials In this study, a 5 inch silicon wafer (Shin-Etsu (M) Sdn. Bhd.) was used as substrate. SU-8 50 photoresist (PR) (MicroChem Corp.) and propylene-glycol-methyl-ether-acetate IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 57 (PGMEA) were used as coating material and developer respectively. The designs were drawn using AutoCAD software and printed on a transparency film using a laser printer (EPSON T60) to generate a photomask. PDMS (Sylgard 184, Dow Corning) was used in replicating micronozzles from SU-8 master. 3.2 Methods The procedure of soft lithography is illustrated in Fig. 3 and described in detail below. (a) Spin coating (b) Soft bake and exposure (c) Post-exposure bake (d) PGMEA development (e) PDMS molding (f) PDMS peel off (g) Sealing Fig. 3: Schematic of soft lithography in fabricating micronozzles. A. Master Preparation Prior to the process, silicon wafer was cleaned with standard RCA cleaning procedure followed by a pre-baked at 200ºC for 5 minutes to remove any residual moisture on wafer surface [15]. About 5 ml of SU-8 PR (1 ml per inch of wafer) was poured onto the wafer and slow spin at 500 rpm for about 20 seconds to spread PR evenly. This was followed by a fast spin at 1250 rpm for 30 seconds for a coating thickness of 100 µm. The spin coated wafer was then soft-baked at 65ºC for 50 minutes. Soft-baking at lower temperature [16] would reduce stress-related crackings. Baking and cooling process were carried out at a slow rate (5 ºC/min) to avoid thermal stress. For coating thickness of 200 µm, the steps were repeated for additional 100 µm. The polymerization of SU-8 PR was assisted using self-customized UV-A radiation system (Philip TL-K 40W/10R) which radiates highly concentrated radiation between 350 nm and 390 nm. Post-exposure bake was carried out at 55ºC for 20 minutes with a ramping rate of 5 ºC/min. After 1 hour of relaxation time, the wafer was immersed into PGMEA and agitated gently to dissolve the un-crosslink SU-8 PR. Developed wafer was silanized with 1,1,2,2 – tetrahydrooctyl-1-trichlorosilanem (UCT, Inc.) to facilitate mold release. The wafer with remaining patterns is known as SU- 8 master. B. Device Fabrication PDMS prepolymer was mixed thoroughly with its curing agent at 10:1 ratio by volume. The mixture was then placed into vacuum desiccators for degassing to remove IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 58 bubbles introduced during mixing process. Degassed mixture was poured onto SU-8 master and cured at 65ºC for 1 hour. PDMS device was then peeled off from the master. C. Device Sealing Fabricated PDMS device was sealed using hydrophilic bond. A piece of PDMS slide (cutting the leftover PDMS from previous steps into desired shape) and PDMS device was oxidized in freshly prepared strong oxidizing solution, a mixture of 98% H2SO4 and 30% H2O2 at 1:1 ratio. The slide was dipped in the solution for 1 minute while the PDMS device for 15 seconds. Dipping for a longer time may damage the microstructures of the device. Both pieces of PDMS were rinsed with distilled water to remove residual acid and dried with compressed stream of air. The device was placed into contact with PDMS slide and heated at 70ºC for 30 minutes. 4. RESULTS AND DISCUSSION Figure 4 shows selected images of micronozzles printed on transparency film (mask) and their replication on PDMS. All images were investigated under inverted light microscope (Olympus, IX 51) with magnification of 40X and bright field mode. In general, all printed features can be distinguished clearly as micronozzles. The printer is capable in printing the smallest nozzle throat width of 53.5 µm as nozzle throat was seen clearly and not filled by printer ink as shown in Fig. 4. Though, in comparison to wt=107 µm, the resolution is not satisfactory as micronozzle expander walls are not smooth in which some ripples were observed. Images of PDMS micronozzles are shown below of each mask pattern for comparison. For wt=53.5 µm and θ=10º, expander section of replicated micronozzle disintegrates from nozzle throat. The SU-8 structural strength is weak at smaller wt, as shown in Fig. 5(a), crack was observed at nozzle throat. It is believed that the crack further propagates during thermal heating. Thermal stress is induced as a result of mismatch in thermal expansion coefficient of silicon and SU-8, leading to the disintegration. For θ=17.5º and wt=53.5 µm, replicated nozzle loses its resolution, a narrow opening at nozzle throat which may retard or block the flow was observed. The resolution improves for wt=107 µm as no potential flow blockage was observed and wt was found close to that in mask. Similar trend was observed for other expander angles. These results suggest wt=107 µm is the more producible nozzle throat width with current equipments. It is supported by SEM images in Fig. 5(b) and (c), both conical micronozzles with wt=107 µm show acceptable accuracy. Current printing device is obviously lack of precision in printing full length spike geometry especially wt=53.5 µm. For wt=107 µm, geometry precision is acceptable up to 80% of total spike length. Images of PDMS aerospike micronozzle in Fig. 4 further justify the producible nozzle throat width is 107 µm as blockage at nozzle throat was observed for wt=53.5 µm. Current equipments unable to produce spike contour with a sharp spike tip. As shown in Fig. 5(d), SU-8 near the sharp tip was not completely removed during the development. As a consequence, the replicated spike tip for both wt are distorted as shown in Fig. 4. IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 59 Fig. 4: Selected images of micronozzles (from left: conical 10º, conical 17.5º and aerospike) with different nozzle throat width (53.5 µm and 107 µm) observed under inverted light microscope. First row is printed patterns on transparency film; second row is the replicated patterns on PDMS. (a) (b) (c) (d) Fig. 5: SEM (FEI, Quanta400F) images: (a) crack at nozzle throat, θ=10º, wt=53.5 µm; (b) conical micronozzle, θ=17.5º, wt=107 µm; (c) conical micronozzle, θ=40º, wt=107 µm; (d) aerospike micronozzle, wt=107 µm. (a) (b) (c) (d) Fig. 6: Size of (a) PDMS conical micronozzle with θ=17.5º wt=107µm; (b) PDMS aerospike micronozzle, wt=107 µm, compared with Malaysian 10 sen coin (dia. ≈1.9 cm) and their SEM images in (c) and (d) respectively. Conical micronozzle Aerospike micronozzle Wt = 116.6µm Wt = 133µm IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 60 Figure 6 (a) and (b) show photographs of PDMS micronozzles in comparison with a Malaysian 10 sen coin (diameter ≈ 1.9 cm). Figure 6 (c) and (d) are SEM images of PDMS conical (θ=17.5º and wt=107 µm) and aerospike (wt=107 µm) micronozzles. Measurements show conical micronozzle has wt=116.6 µm while the aerospike micronozzle has wt=133.3 µm. Further optimization on fabrication procedures is needed to improve wall smoothness and dimensional accuracy. Higher resolution mask is required in fabrication of micronozzles of smaller wt. Truncated aerospike micronozzle should be considered as alternative if better photomask still cannot eliminate the spike tip problem. 5. CONCLUSION In this research, we demonstrated that micronozzles can be fabricated through soft lithography. The technique is inexpensive and relative simple compared to well- established DRIE technique. Conical and aerospike micronozzles were designed and fabricated in two different nozzle throat widths for comparison purpose. Microscopic inspections suggest 107 µm is the more producible dimension with current equipments. The potential application of PDMS-based micronozzle is limited to cold gas microthruster system due to relatively low decomposition temperature of PDMS. Future research would focus on fabrication optimization and the use of alternative structural materials such as ceramic which can work under high temperature. ACKNOWLEDGEMENT The authors would like to acknowledge MOSTI eScience Fund for funding the research. REFERENCES J. Mueller, C. Marrese, J. Polk, E-U. Yang, A. Green, V. White, D. Bame, I. Chakraborty and S. Vargo, “An Overview of MEMS-based micropropulsion development at JPL.”, Acta Astronautica 52 (9-12): 881-895, 2003 R. Osiander, M.A. G. Darrin, and J.L. Champion, “MEMS and microstructures in aerospace applications.”, CRC Press, US, 2006 W.F. Louisos and D.L. Hitt, “Viscous effects on performance of two-dimensional supersonic linear micronozzles.”, Journal of Spacecraft and Rockets, Vol. 45, No. 4, 706- 715, 2008 W.F. Louisos and D.L. Hitt,“Numerical studies of thrust production in 2D supersonic bell micronozzles.”, 44th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, 2008 IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 61 A. Zilic, D.L. Hitt and A.A. Alexeenko,“Numerical simulations of supersonic flow in a linear aerospike micronozzle.”, 37th AIAA Fluid Dynamics Conference & Exhibit, 2007 R.L. Bayt, A.A. Ayon, and K.S. Breuer, “A Performance Evaluation of MEMS-based Micronozzles.”, 33rd AIAA/ASME/SAE/ASEE Joint Propulsion Conf. & Exhibit, 7-9 July, Seattle, WA, 1997 Y. Xia, and G.M. Whitesides, “Soft lithography”, Annual Review of Material Science 28 (1): 153-184, 1998 C.H. Lin, G.B. Lee, B.W. Chang, and G.L. Chang, “A new fabrication process for ultra-thick microfluidic microstructures utilizing SU-8 PR.”, Jour. of Micromech. and Microeng.12 (5), pp. 590-597, 2002 Louisos, W.F., Alexeenko, A.A. Hitt, D.L. and Zilic, A. (2008) Design considerations for supersonic micronozzles. International Journal of Manufacturing Research 3, no. 1: 80- 113. G.P. Sutton, and O. Biblarz, “Rocket propulsion elements.”, 7th Edition, John Wiley & Sons Inc, 2001 K.H. Cheah, J.K. Chin, and K.S. Koh, “Design and Performance Prediction of a 5mN HAN Monopropellant Microthruster”, Int. Conf. on Mech. and Aero. Eng. 16-18 April, Chengdu, China, 2010 W.F. Louisos, and D.L. Hitt, “Numerical studies of thrust production in 2D supersonic bell micronozzles.”, 44th AIAA/ASME/SAE/ASEE Joint Propulsion Conf. & Exh., 21-23 July, Hartford, CT, 2008 J.D. Anderson, “Modern Compressible flow with historical perspective.” 3rd Edition, New York: McGraw-Hill, 2003 G. Angelino, “Approximate method for plug nozzle design.”, AIAA Journal, Vol. 2 (10), 1834-1835, 1964 M. Bachman, “RCA-1 and RCA-2 silicon wafer cleaning.”, INRF Application Note, 1999 T.A. Anhoj, A.M. Jorgensen, D.A. Zauner, and J. Hubner, “The effect of soft bake temperature on the polymerization of SU-8 photoresist.” Journal of Micromechanics and Microengineering, 16, 1819-1824, 2006 NOMENCLATURE Re Reynolds number Ae nozzle exit area, m 2 At nozzle throat area, m2 ρ density, kg/m3 L characteristic length, m V flow velocity, m/s µ viscosity, kg/m2s m mass flow rate, kg/s A area, m2 IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 62 M Mach number k ratio of specific heats α geometric angle with reference axis l length of characteristic segment, m λ l/lt, non-dimensional length for plane nozzle ε A/At, expansion ratio wt nozzle throat width θ expander half-angle B body force due to sediment particles Ca cavitation number CD drag coefficient Cµ anisotropic Diffusion Coefficient Cs1, C s2, C 1ε, C 2ε RS Coefficients dp particle diameter E erosion rate density FD drag force Fcond cavitation condensation coefficient Fvap cavitation vaporization coefficient g gravitational acceleration of the sediment particles k erosion model constant m particle mass n velocity exponent Ps saturation pressure of water p fluid pressure pv vapour pressure RB bubble radius Rep particle Reynolds number IIUM Engineering Journal, Vol. 12, No. 1, 2011 Cheah et al. 63 Rf cavitation rate under relaxation factor f fluid phase volume fractionݎ rnuc nuclei volume fraction p particulate phase volume fractionݎ Uc fluid velocity Up particle velocity U particle relative velocity Vp particle impact velocity Greek letters β particle mass loadings γ particle impact angle ρf fluid density ρp particle density ρs nucleation site density ρr density ratio σ stress tensor σε turbulence model constant σk turbulence model constant for the k-equation σt surface tension coefficient τp particle response time µ fluid dynamic viscosity ν fluid kinematic viscosity