Design and Performance Analysis of Micro Proton Exchange
Membrane Fuel Cells
*ZHONG Zhenzhong (ᄲჲᄱ)**, CHEN Junxun (чࢎ༴) and PENG Ronggui (ఔڬ)
Department of Mechanical Engineering, National Chiao Tung University, Hsinchu 300, China
Abstract This study describes a novel micro proton exchange membrane fuel cell (PEMFC) (active area, 2.5 cm2). The flow field plate is manufactured by applying micro-electromechanical systems (MEMS) technology to silicon substrates to etch flow channels without a gold-coating. Therefore, this investigation used MEMS technology for fabrication of a flow field plate and presents a novel fabrication procedure. Various operating parameters, such as fuel temperature and fuel stoichiometric flow rate, are tested to optimize micro PEMFC performance. A single mi-cro PEMFC using MEMS technology reveals the ideal performance of the proposed fuel cell. The optimal power den-sity approaches 232.75 mW·cmˉ1 when the fuel cell is operated at ambient condition with humidified, heated fuel. Keywords micro proton exchange membrane fuel cell, flow field plate, micro-electromechanical systems technology, silicon
1 INTRODUCTION
A proton exchange membrane fuel cell (PEMFC) with hydrogen as the anode fuel has immense promise as a future power system because of its high efficiency and low level of pollution generated. Research into and development of PEMFCs have garnered consid-erable attention. Hu et al. [1] developed a numerical model for a PEMFC and simulated basic transport phenomena as gas-liquid two-phase flows in a work-ing fuel cell. Shao et al. [2] generated a dynamic thermal- transfer model of a PEMFC stack, which was based on energy conservation to improve the temperature control of the PEMFC stack. Zhang et al. [3] studied the degradation mechanism of key components of in a PEMFC membrane electrode assembly (MEA) and identified feasible measures to avoid degradation. Moreover, the novel MEA is particularly suitable for portable electronic power applications. Notably, PEMFCs have characteristics suited to consumer elec-tronic devices. Consequently, the increased demand for a portable and efficient electrical power supply led to the development of micro fuel cell technology. Conventional fuel cells are manufactured using tradi-tional Computer Numeric Control (CNC) technology; however, micro fuel cells are fabricated using Mi-cro-electromechanical Systems (MEMS) technology. The potential applications of micro fuel cells have increased significantly since 2000, driven by a well-defined need of the consumer electronics indus-try for micro PEMFCs that outperform secondary lithium batteries [4].
Early in 2000, Lee et al. [5] applied microfabrica-tion techniques, such as deep silicon etching, pho-tomasked electroplating, physical vapor deposition, anodic bonding, and spin coating, to silicon wafers to create flow channels. They produced milliwatt fuel
cells using novel techniques and materials. Their mi-cro fuel cells have an efficiency of 150 mA·cmˉ2.
Hahn et al. [6] developed a planar PEMFC with a cross sectional area between 1 mm2 and 1 cm2 at the silicon wafer level. They analyzed various patterning technologies to fabricate micro flow fields, and achieved a power density of 80 mW·cmˉ2 during
long-term operation. Yohtaro [7] applied micro fabri-cation technology to preparation processes for micro fuel cells. Yohtaro showed that micro PEMFCs had higher output power than micro direct methanol fuel cells (DMFCs). Hsieh et al. [8] developed a novel de-sign and fabrication process for a micro fuel-cell flow-field plate with a cross sectional area of 5 cm2 and thickness of 800 ȝm. Their novel design had a power density of 25 mW·cmˉ2 at 0.65 V, and a
reli-able and streli-able output power at ambient temperatures. In another study, Hsieh et al. [9] used the SU-8 photo-resist microfabrication process to fabricate micro PEMFC flow structures. The I-V curve of perform-ance shows a modest voltage drop over a very small range of current. Their work contributed to the low-cost mass-production of small flat single fuel cells with a power density of 30 mW·cmˉ2 at 0.35 V.
Chan et al. [10] developed a micro fuel cell using polymeric micromachining. The membrane electrode assembly was embedded in a gold-coated polymethyl methacrylate substrate. The operating conditions in-cluded dead-ending of hydrogen in the anode. Chan et al. showed that the power output by the 3 cm2 fuel cell was 0.947 W when pure hydrogen was fed to the an-ode and pure oxygen was supplied to the cathan-ode at room temperature. When air was utilized on the cath-ode side, power output reached 0.246 W. The design of the base substrate and electrodes performed well. Hydrogen was set to be dead-ended; thus, the utiliza-tion of hydrogen was nearly 100%.
In 2006, Cha et al. [11] studied transport phenomena Received 2008-09-29, accepted 2008-11-27.
* Supported by the National Science Council (NSC 97-2221-E-009-067).
in micro fuel-cell flow channels. The channels were 500 ȝm, 100 ȝm and 20 ȝm wide and constructed us-ing a structural photopolymer. The effects of channel size and the gas diffusion layer (GDL) thickness on performance were examined. A large pressure drop in very small channels improved the convection of air into the GDL and improved fuel cell performance at low current densities. Cha et al. indicated that using a thin GDL improved micro fuel cell performance. Kim et al. [12] investigated the fabrication process and electrochemical characterization of a miniaturized PEMFC with silicon separators. The micro channels on the silicon separator were created by a photolitho-graphic process. The micro channels were 400 ȝm wide and 230 ȝm deep. The single cell with these sili-con separators had a power density of 203 mW·cmˉ2
at 0.6 V, which is similar to that of a conventional non-silicon fuel cell. Kim et al. indicated that a silicon separator can be applied successfully as an alternative to conventional bipolar plates for a micro PEMFC. Hsieh et al. [13] investigated the operational parame-ters of a H2/air micro PEMFC with different flow
con-figurations by impedance spectroscopy. Their study considered a range of operating parameters, such as backpressure and cell temperature, to determine how flow configuration affected micro fuel cell perform-ance. Optimal operating conditions were utilized to test the micro PEMFC stack [14]. Hsieh et al. demon-strated that the effects of operational parameters, in-cluding magnitude, on stack performance are similar to those on a single cell. Xiao et al. [15], who devel-oped a silicon-based PEMFC power system, indicated that the catalyst with a high active surface area was monolithically integrated with the major fuel cell components using the sputtering method. Lee and Chuang [16] employed porous silicon as a gas diffu-sion layer in a micro fuel cell. The platinum (Pt) cata-lyst was deposited on the surface of, and inside, the porous silicon by physical vapor deposition, to im-prove porous silicon conductivity. Porous silicon with Pt catalyst replaced the conventional GDL, and the Pt metal remaining on the rib was used to form a mi-cro-thermal sensor in a single lithographic process. Lee et al. demonstrated that temperature was almost linearly related to resistance. Feng et al. [17] devel-oped a novel method of producing silicon-based mi-crofabricated electrodes with high electrochemical active surface area (EASA) and high catalyst utiliza-tion. Feng et al. showed that the Si-based membrane electrode assembly (MEA) system exhibited a power density of 23.4 mW·cmˉ2. Seyfang et al. [18]
devel-oped a novel, simplified concept for polymer electro-lyte micro-fuel cells. For different flow fields, a maximum power density was ~415 mW·cmˉ2 at a cell
voltage of 425 mV. Zhang et al. [19] developed a novel porous gas diffusion medium with improved thermal and electrical conductivities and controllable porosity fabricated from a metal foil using micro/nano tech-nology in a PEMFC. The performance of the novel GDL was enhanced by applying a microporous layer (MPL) to the GDL, and by enhancing in-plane
trans-port. Renaud et al. [20] studied the hydrodynamic flow characteristics in microfluidic devices. The 2D model combined with Navier-Stokes equations and the no-slip condition were sufficient to describe and cal-culate flows in 20 ȝm deep microchannels etched in silicon; the derived flows were in good agreement with measured flows.
A micro PEMFC is defined as one that generates ˘5 W of electricity. The preferred substrate material in a flow field plate is silicon because of its suitability to MEMS technology. In this study, MEMS technol-ogy is employed to a silicon wafer to embed flow channels in a micro fuel cell.
This study focuses applies MEMS technology to a silicon wafer to manufacture flow channels and tests micro fuel cell performance under different ambient conditions. Various operational conditions are experi-mentally utilized to optimize cell performance. This work develops a novel and simple method for manu-facturing efficient micro PEMFCs.
2 DESIGN AND FABRICATION OF MICRO PEMFC
A typical micro PEMFC has an end plate, gasket, GDL, MEA, flow field plate and current collector. Fig. 1 (a) presents a schematic diagram of a micro PEMFC. Fig. 1 (b) displays the fabricated single micro PEMFC with all components.
(a) Scheme of micro PEMFC
(b) Picture of micro PEMFC Figure 1 Micro PEMFC
The end plate, made of acrylic resin, is 45 mm× 45 mm×13 mm. The gasket, made of silica gel, is 1 mm thick. This gasket isolates and prevents gas leak-age. The 0.4-mm-thick GDL is made of standard car-bon paper (CARBEL CL GDL). The MEA, which is 0.035 mm thick, catalyst loading of the anode and cathode is 0.5 (mg Pt)·cmˉ2, is a commercial product.
A single micro PEMFC is assembled using two flow field plates. Silicon is used for the substrate in the anode and cathode flow field plates. The micro PEMFC is manufactured using MEMS technology. The objective of the silicon wafer etching process is to create micro PEMFC channels. Photolithography is used to transfer flow channel geometric shapes on a mask onto the surface of a silicon wafer. The steps in the photolithography process are wafer cleaning, de-posit silicon nitride, spin coating the resist, soft baking, mask alignment and exposure, development, and hard baking. Fig. 2 shows the fabrication procedure for the silicon wafer. A silicon wafer 1.016 m in diameter is employed as the substrate in this investigation.
Contaminants on the surface of silicon wafers at the start of the MEMS technology process or accumu-lated during processing must be removed during processing to obtain high-performance and highly re-liable semiconductor devices, and prevent equipment contamination, especially under high-temperature oxidation, diffusion and deposition tubes. In this study, RCA (Radio Corporation of America) cleaning is a four-step process for removing organic contaminants, the oxide layer, ions, and heavy metal contaminants from the wafer surface over a 1.016 m, (100)-oriented 750-ȝm-thick silicon wafer.
A spinner is used to dry the wafer after RCA cleaning. A 200-nm silicon nitride layer is deposited by low-pressure chemical vapor deposition (LPCVD). One side of the silicon wafer is polished, and the sili-con substrate is then spin-coated with a resistive layer 7 ȝm thick at various rotating speeds.
The silicon substrate is then soft baked. Soft baking removes most solvents from the photoresist coating. A positive photoresist is then employed. The photoresist is exposed to ultraviolet (UV) light wher-ever the underlying material must be removed. The g-line is used during the exposure process in this study.
After exposure, heat is applied for approximately 23 min for post-exposure baking. Development is done in the vitrics to form micro flow channels.
The photoresist is then hard-baked for about 1 min. After development, the silicon substrate is ex-amined. The flow structure cavities are patterned by wet etching. Reactive ion etching (RIE) is applied to etch a silicon nitride layer 300 ȝm wide. Notably, RIE, an etching technology used during microfabrication, uses chemically reactive plasma to remove materials deposited on wafers. Plasma is generated under low pressure in an electromagnetic field. However, the silicon nitride layer is too thin for the fuel to flow from the inlet to the outlet. Hence, the exposed silicon wafer is etched to pattern micro PEMFC channels in an aqueous solution of 45% KOH at 80°C. Etched channel depth is 500 ȝm. This depth promotes gas uniformity, water management and reduced flow re-sistance. Next, the photoresist is removed from the silicon substrate. The final step employs phosphoric acid to remove the remaining silicon nitride at 180°C.
Three serpentine channels, each 300 ȝm wide and 500 ȝm deep, are etched into the sides of the an-ode and cathan-ode. The fuel holes are drilled into the rear of silicon wafer. The oversize flow field plates are 25 mm×25 mm. A silicon wafer 1.016 m in diameter is used to make four flow field plates. Fig. 3 shows the silicon substrate after etching.
Current collectors provide an electrical pad for transmission of electricity to the external load. Unlike metal, silicon is not considered an ideal current col-lector due to its high electrical resistance. In this in-vestigation, the silicon wafer is only used as a fuel carrier. Therefore, an alternative design in this study uses an MEA and GDL sandwiched between current collectors made of brass foil; this design replaces conventional bipolar plates. Consequently, silicon sub-strates are rarely used for flow field plates as electrons
do not pass through silicon substrates. The advantage of this novel design is that no Au electroplate is quired on the top of the silicon wafer. The only re-quirement is that current collectors be hollowed out to enable fuel to diffuse to the MEA.
In this study of the micro PEMFC fabrication procedure, the MEA is sandwiched between two GDLs, two current collectors, two flow field plates, two gaskets and two end plates. The assembled struc-tures are clamped together with four screws.
3 EXPERIMENTAL
To obtain a micro PEMFC polarization curve, the micro PEMFC test has a gas supply unit, electrical load, power supply, computer, fuel humidifier, fuel temperature controller and mass flow controller. Fig. 4 shows the experimental setup. The gas supply unit provides fuel and the oxidant to the micro fuel cell.
The electrical load is a device that simulates loading on an electrical circuit. Counter to the current source, the electrical load is a current sink. The load current is regulated electronically. The range of measured cur-rents is 015 A and maximum power output is 75 W. The power supply supplies electricity to each compo-nent in this experimental test. A computer regulates all installations orders and evaluates experimental results. The micro PEMFC structure cannot be heated because the end plate is made of acrylic resin. Therefore, a humidifier and temperature controller are used to hu-midify and heat the fuel to improve micro PEMFC performance. The flow rates of hydrogen, oxygen and air are regulated by mass flow controllers. Notably, the flow rate through this mass flow controller cannot be ˘10 cm3·minˉ1 in the anode and 20 cm3·minˉ1 in
the cathode. The micro PEMFC is connected to the electrical load. Polarization curves are plotted after the micro PEMFC reaches a steady state. Various fuel temperatures and fuel flow rates are tested. The fuel flow rates that generate a current of 1 A are calculated based on the theoretical volume. The flow rates of hydrogen and oxygen are 7.6 and 3.8 cm3·minˉ1·Aˉ1,
respectively. Therefore, the flow rates are multiplied by the stoichiometric value of 1.37 to obtain a value of 10.4 cm3·minˉ1·Aˉ1 for the anode, and 1.84 to yield
7.0 cm3·minˉ1·Aˉ1 for the cathode.
4 RESULTS AND DISCUSSION
The objective of this study is to develop a proc-ess for fabricating a micro PEMFC with a simple structure using MEMS technology applied to a silicon substrate. A single micro PEMFC is assembled and used in the polarization curve experiment. Optimal Figure 3 Silicon substrates after etching process
operating conditions are identified.
Figure 5 plots the micro PEMFC polarization curves obtained at fuel temperatures of 30, 40, 50, 60, 70 and 80°C. The experiment is performed at ambient pressure and temperature; pure humidified hydrogen and oxygen are fed into the micro PEMFC. The stoichiometric flow rate is fixed at 3 (H2 is 31.2
cm3·minˉ1 and O
2 is 21.0 cm3·minˉ1). The current
densities are 436.77, 421.82, 402.00, 469.58, 488.27 and 490.83 mA·cmˉ2, respectively, at 0.46 V. The
performance at a fuel temperature of 70°C is similar to that at 80°C. The average power densities are 222.16 and 223.32 mW·cmˉ2, indicating that the best fuel
temperatures for the micro PEMFC are 70 and 80°C in this experiment. Experimental results show that a high fuel temperature increases gas diffusivity and mem-brane conductivity. Therefore, a high fuel temperature induces electrochemical reactions. This experimental result indicates that cell performance is increases as fuel temperature increases. Restated, micro PEMFC performance depends on fuel temperature.
Figure 5 I-V curves at different fuel temperatures ƽ 30°C;Ʒ40°C;Ƶ50°C;ƾ60°C; × 70°C;ƹ80°C
Four operating fuel stoichiometric flow rates are employed to study their effects on micro PEMFC per-formance. The fuel stoichiometric flow rates are 3, 4, 5 and 6 (H2, 31.2, 41.6, 52.0, 62.4 cm3·minˉ1; O2, 21.0,
28.0, 35.0, 42.0 cm3·minˉ1). The experiment is
per-formed at ambient pressure and temperature. Pure humidified hydrogen and oxygen are fed into the mi-cro PEMFC while fuel temperature is maintained at 70°C. Fig. 6 plots experimental results obtained with the four fuel stoichiometric flow rates. When the ini-tial fuel stoichiometric flow rate is 6, cell performance is poorest; that is, current density at 0.46 V is 453.27 mA·cmˉ2. For fuel stoichiometric flow rates of 3, 4
and 5, the current densities at 0.46 V are 468.00, 505.98 and 504.60 mA·cmˉ2, respectively. Therefore,
internal flooding is not obvious. However, a fuel stoichiometric flow rate of 4 yields a performance slightly better than that obtained with other stoichiometric flow rates. Maximum power output with a fuel stoichiometric flow rate of 4 is 232.75 mW·cmˉ2. This experimental result indicates that the
fuel stoichiometric flow rate markedly affects cell performance, i.e., micro fuel cell performance im-proves gradually as the fuel stoichiometric flow rates increase. Increasing gas velocity improves the re-moval of water produced by the cathode. However, the fuel stoichiometric flow rate of 6 dries out the membrane, which increases electrical resistance. Therefore, the fuel stoichiometric flow rate of 6 gen-erates the worst performance.
Figure 7 presents the micro PEMFC polarization behavior when the fuel stoichiometric flow rates are changed from H2/O2 to 3, 4 and 5. The experiment is
performed at ambient pressure and temperature. Pure humidified hydrogen and oxygen are fed into the mi-cro PEMFC. Fuel temperature is controlled at 80°C. As the stoichiometric flow rate increases, cell per-formance worsens. As mentioned, increasing fuel temperature accelerates reactions, and increases the amount of liquid water produced, which adversely affects cell performance. Humidified H2 and O2 are
fed into the micro fuel cell. Accordingly, the amount of liquid water is generated increases when the high Figure 6 Comparison of polarization curves at four oper-ating fuel stoichiometric flow rates at 70°C
ƽ 70°C stoich 3;Ʒ70°C stoich 4;Ƶ70°C stoich 5;ƾ70°C stoich 6
Figure 7 Comparison of polarization curves at three oper-ating fuel stoichiometric flow rates at 80°C
current density. Notably, high humidification can cause excessive water accumulation that can flood the channels. Therefore, when the fuel temperature is 80°C and fuel flow rates are high, cell performance is worse than that when fuel temperature is 70°C. How-ever, when the stoichiometric flow rate of fuel is 4, the reaction is optimal at a fuel temperature of 70°C.
Figure 8 plots polarization curves of pure hu-midified hydrogen and air (70°C), which are used as reactant gases under ambient conditions. The flow of H2/air are controlled with stoichiometric flow rates of
1, 2 and 3 (H2 is 10.4, 20.8, 31.2 cm3·minˉ1 and air is
35.0, 70.0, 105.0 cm3·minˉ1). The performance is
poorer than when pure oxygen is used in the cathode; experimental results demonstrate that current densities are 139.03, 149.90 and 166.65 mA·cmˉ2, respectively,
at 0.46 V. This experiment indicates that a fixed stoichiometric flow rate of 5 optimizes performance of micro PEMFC when the cathode fuel is air.
Figure 8 Comparison of polarization curves at three op-erating fuel (H2/air) stoichiometric flow rates at 70°C
ƽ 70°C H2/air stoich 1;Ʒ70°C H2/air stoich 2; Ƶ70°C H2/air stoich 3
Additionally, experimental results demonstrate that the performance of micro PEMFC exceeds expectations. The optimal current density is 505.98 mA·cmˉ2 at 0.46 V
for a pure hydrogen and pure oxygen stoichiometric flow rate of 4 at a fuel temperature of 70°C.
5 CONCLUSIONS
Portable consumer electronic devices require small, lightweight power supplies with high electricity capacity. Micro PEMFCs meet this requirement. The study employs MEMS technology to etch flow field channels into a silicon substrate. The reaction area of this single micro PEMFC is 2.5 cm2. A single micro PEMFC is successfully fabricated and favorable per-formance achieved. The optimal power density ap-proaches 232.75 mW·cmˉ2 at 0.46 V when the cell is
operated at ambient conditions with humidified and heated fuel. The objective of future work is to assem-ble serial interconnected single micro PEMFCs in a micro PEMFC stack.
ACKNOWLEDGEMENTS
Ted Knoy is appreciated for his editorial assistance. REFERENCES
1 Hu, M.R., Zhu, X.J., Gu, A.Z., “Simulation of the internal transport phenomena for PEM fuel cells with different modes of flow”, Chin. J. Chem. Eng., 12, 1426 (2004).
2 Shao, Q.L., Wei, D., Gao, G.Y., Zhu, X.J., “Dynamic thermal model and temperature control of proton exchange membrane fuel cell stack”, Chin. J. Chem. Eng., 13, 218224 (2005).
3 Zhang, S.S., Yu, H.M., Zhu, H., Hou, J.B., Yi, B.L., Ming, P.W., “Effects of freeze/thaw cycles and gas purging method on polymer electrolyte membrane fuel cells”, Chin. J. Chem. Eng., 14, 802805 (2006). 4 Vielstich, W., Lamm, A., Gasteiger, H. A., Handbook of Fuel Cells,
West Sussex, England (2003).
5 Lee, S.J., Cha, S.W., Liu, Y., O’Hayre, R., Prinz, F.B., “High power-density polymer-electrolyte fuel cells by microfabrication”, Electrochemical Society Proceedings (2000).
6 Hahn, R., Wagner, S., Schmitz, A., Reichl, H., “Development of a planar micro fuel cell with thin film and micro patterning technolo-gies”, Journal of Power Sources, 131, 7378 (2004).
7 Yohtaro, Y., “Application of MEMS technology to micro fuel cells”, Electrochimica Acta, 50, 663666 (2004).
8 Hsieh, S.S., Kuo, J.K., Hwang, C.F., Tsai, H.H., “A novel design and microfabrication for a micro PEMFC”, Microsystem Technolo-gies-micro-and Nanosystems-information Storage and Processing Systems, 10, 121126 (2004).
9 Hsieh, S.S., Hwang, C.F., Kuo, J.K., Tsai, H.H., Yang, S.H., “SU-8 flow field plates for a micro PEMFC”, Journal of Solid State Elec-trochemistry, 9, 121131 (2005).
10 Chan, S.H., Nguyen, N.T., Xia, Z., Wu, Z., “Development of a polymeric micro fuel cell containing laser-micromachined flow channels”, Journal of Micromechanics and Microengineering, 15, 231236 (2005).
11 Cha, S.W., O’Hayre, R., Park, Y.I., Prinz, F.B., “Electrochemical impedance investigation of flooding in micro-flow channels for pro-ton exchange membrane fuel cells”, Journal of Power Sources, 161, 138142 (2006).
12 Kim, J.Y., Kwon, O.J., Hwang, S.M., Kang, M.S., Kim, J.J., “Develop-ment of a miniaturized polymer electrolyte membrane fuel cell with silicon separators”, Journal of Power Sources, 161, 432436 (2006). 13 Hsieh, S.S., Yang, S.H., Feng, C.L., “Characterization of the opera-tional parameters of a H2/air micro PEMFC with different flow fields by impedance spectroscopy”, Journal of Power Sources, 162, 262270 (2006).
14 Hsieh, S.S., Feng, C.L., Huang, C.F., “Development and perform-ance analysis of a H2/air micro PEM fuel cell stack”, Journal of Power Sources, 163, 440449 (2006).
15 Xiao, Z., Yan, G., Feng, C., Chan, P.C.H., Hsing, I.M., “A sili-con-based fuel cell micro power system using a microfabrication technique”, Journal of Micromechanics and Microengineering, 16, 20142020 (2006).
16 Lee, C.Y., Chuang, C.W., “A novel integration approach for com-bining the components to minimize a micro-fuel cell”, Journal of Power Sources, 172, 115120 (2007).
17 Feng, C., Chan, P.C.H., Hsing, I.M., “Catalyzed microelectrode me-diated by polypyrrole/nafion composite film for microfabricated fuel cell applications”, Electrochemistry Communications, 9, 8993 (2007). 18 Seyfang, B.C., Kuhnke, M., Lippert, T., Scherer, G.G., Wokaum, A.,
“A novel, simplified micro-PEFC concept employing glassy carbon micro-structures”, Electrochemistry Communications, 9, 19581962 (2007).
19 Zhang, F.Y., Advani, S.G., Prasad, A.K., “Performance of a metallic gas diffusion layer for PEM fuel cells”, Journal of Power Sources,
176, 293298 (2008).
20 Renaud, L., Malhaire, C., Kleimann, P., Barbier, D., Morin, P., “Theoretical and experimental studies of microflows in silicon mi-crochannels”, Materials Science and Engineering C, 28, 910917 (2008).
