<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJAppS</journal-id><journal-title-group><journal-title>Open Journal of Applied Sciences</journal-title></journal-title-group><issn pub-type="epub">2165-3917</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojapps.2019.94016</article-id><article-id pub-id-type="publisher-id">OJAppS-91910</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biomedical&amp;Life Sciences</subject><subject> Chemistry&amp;Materials Science</subject><subject> Computer Science&amp;Communications</subject><subject> Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Design and Analysis of MEMS Based Aluminum Nitride (AlN), Lithium Niobate (LiNbO&lt;sub&gt;3&lt;/sub&gt;) and Zinc Oxide (ZnO) Cantilever with Different Substrate Materials for Piezoelectric Vibration Energy Harvesters Using COMSOL Multiphysics Software
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ahmad</surname><given-names>M. Alsaad</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ahmad</surname><given-names>A. Ahmad</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Qais</surname><given-names>M. Al-Bataineh</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Nermeen</surname><given-names>S. Daoud</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mais</surname><given-names>H. Khazaleh</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Physical Sciences, Jordan University of Science and Technology, Irbid, Jordan</addr-line></aff><aff id="aff2"><addr-line>Department of Physics, Faculty of Sciences, University of Hafr Al-Batin, Hafr Al-Batin, Saudi Arabia</addr-line></aff><pub-date pub-type="epub"><day>11</day><month>04</month><year>2019</year></pub-date><volume>09</volume><issue>04</issue><fpage>181</fpage><lpage>197</lpage><history><date date-type="received"><day>12,</day>	<month>March</month>	<year>2019</year></date><date date-type="rev-recd"><day>16,</day>	<month>April</month>	<year>2019</year>	</date><date date-type="accepted"><day>19,</day>	<month>April</month>	<year>2019</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Interest in energy harvesters has grown rapidly over the last decade. The cantilever shaped piezoelectric energy harvesting beam is one of the most employed designs, due to its simplicity and flexibility for further performance enhancement. The research effort in t
  he MEMS Piezoelectric vibration energy harvester designed using three types of cantilever materials, Lithium Niobate (LiNbO<sub>3</sub>), Aluminum Nitride (AlN) and Zinc Oxide (ZnO) with different substrate materials: aluminum, steel and silicon using COMSOL Multiphysics package were designed and analyzed. Voltage, mechanical power and electrical power versus frequency for different cantilever materials and substrates were modeled and simulated using Finite element method (FEM). The resonant frequencies of the LiNbO<sub>3</sub>/Al, AlN/Al and ZnO/Al systems were found to be 187.5 Hz, 279.5 Hz and 173.5 Hz, respectively. We found that ZnO/Al system yields optimum voltage and electrical power values of 8.2 V and 2.8 mW, respectively. For ZnO cantilever on aluminum, steel and silicon substrates, we found the
   
  resonant frequencies to be 173.5 Hz, 170 Hz and 175 Hz, respectively. Interestingly, ZnO/steel yields optimal voltage and electrical power values of 9.83 V and 4.02 mW, respectively. Furthermore, all systems were studied at different differentiate frequencies. We found that voltage and electrical power have increased as the acceleration has increased.
 
</p></abstract><kwd-group><kwd>MEMS</kwd><kwd> Piezoelectric</kwd><kwd> Energy Harvester</kwd><kwd> Cantilever</kwd><kwd> Lithium Niobate  (LiNbO&lt;sub&gt;3&lt;/sub&gt;)</kwd><kwd> Aluminum Nitride (AlN)</kwd><kwd> Zinc Oxide (ZnO)</kwd><kwd> Aluminium  Substrate</kwd><kwd> Steel Substrate</kwd><kwd> Silicon Substrate</kwd><kwd> COMSOL</kwd><kwd> Finite Element  Method</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Micro-electro mechanical system (MEMS) devices have a wide area of applications as pressure sensors, accelerometers, gyroscopes etc. [<xref ref-type="bibr" rid="scirp.91910-ref1">1</xref>] . In recent decades, energy harvesting from ambient vibrations of natural environments, such as air flows and human motions, which are obtainable universally and permanently, has attracted much attention of many researches [<xref ref-type="bibr" rid="scirp.91910-ref2">2</xref>] . The investigations in vibration-based energy harvester using piezoelectric transducers have been attracting intensive attention in the research sector of sustainable energy [<xref ref-type="bibr" rid="scirp.91910-ref3">3</xref>] . Piezoelectric energy harvesters (PEH) have the distinctive capacity to convert ambient vibration energy from the surrounding environment into electrical power [<xref ref-type="bibr" rid="scirp.91910-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref6">6</xref>] . A typical PEH system usually consists of two parts: a mechanical structure and an energy harvesting circuit. The coupling between these two parts and the multidisciplinary nature of this field lead to a substantial challenge in modeling PEH system [<xref ref-type="bibr" rid="scirp.91910-ref7">7</xref>] . Piezoelectric materials are widely used in vibration energy harvesters (VEH) as mechanical-to-electrical transducers due to their relatively high power density [<xref ref-type="bibr" rid="scirp.91910-ref8">8</xref>] , scalability and compatibility with conventional integrated circuit technologies [<xref ref-type="bibr" rid="scirp.91910-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref10">10</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref11">11</xref>] .</p><p>Simulation modeling solves real-world problems safely and efficiently. It provides an important method of analysis which is easily verified, communicated, and understood. Across industries and disciplines, simulation modeling provides cherished solutions by giving clear insights into complex systems. Mathematical modeling and simulations have been used to predict the optimized properties and to evaluate the performance of systems similar to our proposed system [<xref ref-type="bibr" rid="scirp.91910-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref13">13</xref>] . COMSOL Multiphysics&#174; (known as FEMLAB before 2005) is a software package designed to investigate a wide range of physical phenomena of different systems using finite element method [<xref ref-type="bibr" rid="scirp.91910-ref14">14</xref>] .</p><p>AlN, LiNbO<sub>3</sub> and ZnO have extensive applications in several fields because of their outstanding piezoelectric and pyroelectric properties, as well as, elastic- and electro-optic effects [<xref ref-type="bibr" rid="scirp.91910-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref17">17</xref>] . Several previous works have studied PEH based on AlN [<xref ref-type="bibr" rid="scirp.91910-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref19">19</xref>] , LiNbO<sub>3</sub> [<xref ref-type="bibr" rid="scirp.91910-ref20">20</xref>] and ZnO [<xref ref-type="bibr" rid="scirp.91910-ref21">21</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref22">22</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref23">23</xref>] Cantilevers due to their significant contributions in the modern technology of MEMS devices.</p><p>In this work, PEH systems were modelled by COMSOL Multiphysics software using AlN, LiNbO<sub>3</sub> and ZnO Cantilevers assembled on Aluminum, Steel and Silicon substrates. The main objective of this work is to reveal the most efficient cantilever/substrate combination that has the optimized electrical and piezoelectric properties. Particularly, we investigate the behavior of the output voltage, power (mechanical and electrical) as well as piezoelectric response a function of frequency and acceleration.</p></sec><sec id="s2"><title>2. Mathematical Model for Piezoelectric Beam</title><p>Cantilever structure is the most common structure used for piezoelectric energy harvester [<xref ref-type="bibr" rid="scirp.91910-ref24">24</xref>] . Cantilever can be defined as a structure with one end fixed and other end free to vibrate, made up of one or more layer of piezoelectric material bonded to an elastic metal in order to increase the sensitivity of the structure and reduce the brittleness of the piezoelectric layer [<xref ref-type="bibr" rid="scirp.91910-ref25">25</xref>] .</p><p>Electric charge is generated when mechanical stress is applied on a piezoelectric material. IEEE standard on piezoelectricity has given different form of piezoelectric constitutive equations. Strain-charge form has been used for cantilever structure and the equations are:</p><p>S = s E T + d i j E &#175; (1)</p><p>D = d i j T + ε T E &#175; (2)</p><p>where: S is mechanical strain, s E is elastic compliance tensor (Pa<sup>−1</sup>), T is mechanical stress vector (N∙m<sup>−2</sup>), E &#175; is electric field vector (V∙m<sup>−1</sup>), D is electrical displacement (C∙m<sup>−2</sup>), ε T is dielectric permittivity tensor (F∙m<sup>−1</sup>) and d i j is electro-mechanical coupling factor (C∙N<sup>−1</sup>), where i is the polarization direction and j is the strain direction.</p><p>There are two approaches used in the design of a piezoelectric harvester: longitudinal and transversal approaches. The first is the longitudinal mode ( d 31 ) where the polarization of the beam is laterally developed in the deposited film. The frequently used is the transversal mode ( d 33 ) where the polarization of the beam is perpendicular to the deposited film.</p><p>The resonant frequency is the most important parameter of a vibration energy harvesting device. It is calculated by using the given equation [<xref ref-type="bibr" rid="scirp.91910-ref26">26</xref>] :</p><p>f n = v n 2 2π 1 L 2 D p m (3)</p><p>where: v n = 1.875 for first mode, m is the mass per unit area and D p is the bending modulus which is a function of Young’s modulus and thickness of the substrate and expressed by:</p><p>D p = E P 2 t p 4 + E s 2 t s 4 + 2 E S E P t p t s ( 2 t p 2 + 2 t s 2 + 3 t p t s ) 12 ( E P t p + E S t s ) (4)</p><p>Hence the variation of resonant frequency is</p><p>f n ∝ 1 L 2 t p t s (5)</p><p>where: L: length of the harvester, t p : thickness of the piezoelectric layer and t s : thickness of the substrate layer. Cantilever oscillates when placed in vibrating environment. The oscillations attain the optimum peak as the vibration frequency of the environment matches the resonance frequency of the cantilever structure, and damps out significantly for all other frequencies. The frequency of the source vibrations present in the environment mostly has frequencies in the range of 50 - 200 Hz. The proof mass lowers the resonance frequency of the cantilever by order of few Hz, which is normally the order of frequency of vibration present in the nature. Proof mass also increases the amount of deflection, hence increasing the stress at the fixed end due to which charge is generated in the cantilever structure. Electrical output voltage is highest when stress is maximum, which is occurs at the resonance frequency [<xref ref-type="bibr" rid="scirp.91910-ref25">25</xref>] .</p></sec><sec id="s3"><title>3. COMSOL Multiphysics Simulation</title><p>To simulate the piezoelectric vibration energy harvester using COMSOL Multiphysics, the piezoelectric devices Multiphysics interface and electrical circuit have to be chosen. Solid mechanics and electrostatics are combined within the piezoelectric devices Multiphysics interface with the constitutive relationships required to model the piezoelectric device.</p><sec id="s3_1"><title>3.1. Geometric Modeling</title><p>2D geometry is considered for the simulations. The cantilever contains aluminum, steel or silicon thin substrate of 40 &#181;m thickness (<xref ref-type="table" rid="table1">Table 1</xref> displays the properties of substrate materials) coated by two layers of AlN, LiNbO<sub>3</sub> or ZnO of 60 &#181;m in thickness‎. The dimensions of the cantilever is (20 &#215; 14) mm, with a mass of aluminum block of dimension (4 &#215; 14 &#215; 1.7) mm on the vibrating end of the cantilever, as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p></sec><sec id="s3_2"><title>3.2. Boundary Settings</title><p>The base of our simulated piezoelectric vibration energy harvester is designed to have one end fixed and the other end is free to vibrate with a mass of aluminum is attached to it. The upper end of the cantilever is taken to be free in response to the applied force. The upper surface of the AlN, LiNbO<sub>3</sub> or ZnO is taken to be ground and the lower surface is connected to external circuit.</p></sec><sec id="s3_3"><title>3.3. Meshing</title><p>Before starting simulation, the piezoelectric vibration energy harvester has be divided into small areas; each is called a “mesh”. In this work, the mesh was</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The properties of substrate materials</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Substrate Materials</th><th align="center" valign="middle" >Density (kg/m<sup>3</sup>)</th><th align="center" valign="middle" >Young’s modulus (10<sup>9</sup> Pa)</th><th align="center" valign="middle" >Poisson’s ratio</th></tr></thead><tr><td align="center" valign="middle" >Aluminum</td><td align="center" valign="middle" >2700</td><td align="center" valign="middle" >70</td><td align="center" valign="middle" >0.33</td></tr><tr><td align="center" valign="middle" >Steel</td><td align="center" valign="middle" >7850</td><td align="center" valign="middle" >200</td><td align="center" valign="middle" >0.30</td></tr><tr><td align="center" valign="middle" >Silicon</td><td align="center" valign="middle" >2329</td><td align="center" valign="middle" >170</td><td align="center" valign="middle" >0.28</td></tr></tbody></table></table-wrap><p>taken to be “Free triangular” with element size parameter ranging between 0.002 and 0.02. <xref ref-type="fig" rid="fig2">Figure 2</xref> illustrates the obtained mesh for a piezoelectric vibration energy harvester.</p></sec><sec id="s3_4"><title>3.4. Simulation Results</title><p>The mass attached to the vibrating end of the piezoelectric vibration energy harvester is simulated to obtain the required mass needed to reach the resonant frequency that yields the maximum voltage, as displayed in <xref ref-type="fig" rid="fig3">Figure 3</xref>(a). <xref ref-type="fig" rid="fig3">Figure 3</xref>(b) demonstrates the obtained electric force at resonant frequency. We found that the resonant frequency of AlN/Al, LiNbO<sub>3</sub>/Al and ZnO/Al are 279.5 Hz, 187.5 Hz and 173.5 Hz, respectively.</p></sec></sec><sec id="s4"><title>4. Results and Discussions</title><p>In this section, we describe and interpret our results on the piezoelectric vibration energy harvester output voltage and power (mechanical and electrical) behavior of the LiNbO<sub>3</sub>, AlN and ZnO cantilevers assembled on different types of substrates (aluminum, steel and silicon) as functions of frequency and acceleration responses.</p><sec id="s4_1"><title>4.1. LiNbO<sub>3</sub>, AlN and ZnO Cantilevers with Aluminum Substrate</title><sec id="s4_1_1"><title>4.1.1. Frequency Response</title><p><xref ref-type="fig" rid="fig4">Figure 4</xref> shows the measured voltage of LiNbO<sub>3</sub>, AlN and ZnO cantilevers assembled on Al substrate as a function of frequency, with a host acceleration of 1 g (g is the acceleration due to gravity) and load resistance of 12 kΩ. As can be seen from <xref ref-type="fig" rid="fig4">Figure 4</xref>, the resonant frequency demonstrated by LiNbO<sub>3</sub>, AlN and ZnO cantilevers is 187.5 Hz, 279.5 Hz and 173.5 Hz, respectively. <xref ref-type="fig" rid="fig4">Figure 4</xref> clearly shows that ZnO cantilever exhibits the maximum voltage value of 8.2 V at a resonant frequency of 173.5 Hz. <xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref> illustrate the calculated</p><p>mechanical power input and electrical output power of the piezoelectric vibration energy harvester designed using LiNbO<sub>3</sub>, AlN and ZnO cantilevers. Obviously, the ZnO cantilever attains the maximum electrical power output of 2.8 mW at a resonant frequency of 173.5 Hz.</p><p>The ZnO/Al cantilever attains the maximum voltage and output electrical power at the resonant frequency in comparison with the other two investigated systems since it exhibits attractive exceptionally high piezoelectric properties. Piezoelectric properties of materials depend on three main parameters: structure, piezoelectric coefficient ( e 13 and e 33 ) and dipole moment. The ZnO and AlN both adopt a wurtzite structure with a space group of (C<sub>6v</sub><sup>4</sup>-P6<sub>3</sub>mc) in their equilibrium phase, thus, exhibiting interesting piezoelectric properties, high voltage and output electrical power at the resonant frequency. On the other hand, LiNbO<sub>3</sub> exhibits a trigonal structure with a space group(R3c), as a result, it attains the minimum voltage and output electrical power at the resonant frequency. We found that AlN piezoelectric coefficients ( e 13 = − 0.580   C ⋅ m − 2 and e 33 = 1.550   C ⋅ m − 2 ), respectively that are close to those of ZnO ( e 13 = − 0.567   C ⋅ m − 2 and e 33 = 1.320   C ⋅ m − 2 ). However, the dipole moment of ZnO is ( μ = 0.345   Debey ) that is much larger than that of AlN ( μ = 0.096   Debey ), as shown in <xref ref-type="table" rid="table2">Table 2</xref>. This is because Zn<sup>2+</sup> cations and O<sup>2−</sup> anions change their position under a stress (force) in a nonhomogenous way. Thus an induced net polarization develops in the material. Moreover, electric dipoles which created by the non-homogenous distribution of Zn<sup>2+</sup> cations and O<sup>2−</sup> anions would be sustained, as long as external stress is applied [<xref ref-type="bibr" rid="scirp.91910-ref27">27</xref>] [<xref ref-type="bibr" rid="scirp.91910-ref28">28</xref>] . The induced polarization caused by the strain at the ZnO/Al interface is large due to the lattice mismatch between ZnO and Al substrate.</p><p>The AlN, and its related alloys are wide (direct) band gap semiconductors with high thermal and mechanical stability. They have attracted great attention for the fabrication of optoelectronic devices operating in the visible and ultraviolet regions at high power and under harsh environmental conditions. The determination of strains at the nanoscale is essential for the development of new electronic devices. The wurtzite hexagonal phase is energetically more stable for AlN and its alloys. This crystalline structure is non-centrosymmetric with a</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Comparative parameters for ZnO, AlN and LiNbO<sub>3</sub> cantilevers</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >ZnO</th><th align="center" valign="middle" >AlN</th><th align="center" valign="middle" >LiNbO<sub>3</sub></th></tr></thead><tr><td align="center" valign="middle" >Structure</td><td align="center" valign="middle" >Wurtzite</td><td align="center" valign="middle" >Wurtzite</td><td align="center" valign="middle" >Trigonal</td></tr><tr><td align="center" valign="middle" >Space Group</td><td align="center" valign="middle" >C<sub>6v</sub><sup>4</sup>-P6<sub>3</sub>mc</td><td align="center" valign="middle" >C<sub>6v</sub><sup>4</sup>-P6<sub>3</sub>mc</td><td align="center" valign="middle" >R<sub>3</sub>c</td></tr><tr><td align="center" valign="middle" >e<sub>13</sub> (C∙m<sup>−2</sup>)</td><td align="center" valign="middle" >−0.567</td><td align="center" valign="middle" >−0.580</td><td align="center" valign="middle" >0.194</td></tr><tr><td align="center" valign="middle" >e<sub>33</sub> (C∙m<sup>−2</sup>)</td><td align="center" valign="middle" >1.320</td><td align="center" valign="middle" >1.550</td><td align="center" valign="middle" >1.309</td></tr><tr><td align="center" valign="middle" >Young’s modulus (GPa)</td><td align="center" valign="middle" >210 [<xref ref-type="bibr" rid="scirp.91910-ref29">29</xref>]</td><td align="center" valign="middle" >344 [<xref ref-type="bibr" rid="scirp.91910-ref30">30</xref>]</td><td align="center" valign="middle" >170 [<xref ref-type="bibr" rid="scirp.91910-ref31">31</xref>]</td></tr><tr><td align="center" valign="middle" >Dipole moment μ (Debey)</td><td align="center" valign="middle" >0.345 [<xref ref-type="bibr" rid="scirp.91910-ref32">32</xref>]</td><td align="center" valign="middle" >0.096 [<xref ref-type="bibr" rid="scirp.91910-ref33">33</xref>]</td><td align="center" valign="middle" >4.40 [<xref ref-type="bibr" rid="scirp.91910-ref34">34</xref>]</td></tr></tbody></table></table-wrap><p>singular polar axis causing the formation of anelectric dipole in the unit cell due to the lack of coincidence of the centre of mass of the negative charge in the N tetrahedrons and the positive charge of the Al atom. The dipolegives rise to a spontaneous polarization. Moreover, in the presence of strain, apiezoelectric polarization is created at the AlN/Al interface due to the lattice mismatch between AlN and the Al substrates that can create fields in the MV/cm range. However, the lattice mismatch in ZnO/Al cantilevers is much larger than lattice mismatch in AlN/Al cantilevers.</p></sec><sec id="s4_1_2"><title>4.1.2. Acceleration Response</title><p><xref ref-type="fig" rid="fig7">Figure 7</xref> displays the calculated voltage of LiNbO<sub>3</sub>, AlN and ZnO cantilevers as a function of host acceleration calculated specifically at the resonant frequency of each of the investigated cantilevers and by keeping the load resistance fixed at 12 kΩ for all cases. Our results indicate that the voltage of the piezoelectric vibration energy harvester increases linearly as the host acceleration increases from 0.25g - 2 g, in excellent agreement with the findings of E. K. Reilly et al. study [<xref ref-type="bibr" rid="scirp.91910-ref35">35</xref>] . In addition, our results clearly demonstrate that ZnO cantilever exhibits larger voltage values than those attained by LiNbO<sub>3</sub> and AlN cantilevers.</p><p><xref ref-type="fig" rid="fig8">Figure 8</xref> and <xref ref-type="fig" rid="fig9">Figure 9</xref> illustrate the calculated mechanical input power and the obtained electrical power output of the piezoelectric vibration energy harvester of the LiNbO<sub>3</sub>, ZnO and AlN cantilevers. Apparently, the ZnO cantilever sustains the best electrical power output at the 173.5 Hz resonant frequency.</p><p>Inspecting the obtained results carefully indicate that among the three different cantilevers investigated, ZnO cantilevers yields better voltage and electrical output power. Now, we turn our attention to focus on investigating the electrical properties of ZnO cantilevers assembled on three different substrates, namely, aluminum, steel and silicon substrates.</p></sec></sec><sec id="s4_2"><title>4.2. ZnO Cantilevers with Aluminum, Steel and Silicon Substrates</title><sec id="s4_2_1"><title>4.2.1. Frequency Response</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref>0 illustrates the calculated voltage of ZnO cantilevers assembled on aluminum, steel and silicon substrates as a function of frequency keeping the host acceleration fixed at 1 g and a predetermined load resistance of 12 kΩ. As can be clearly seen from <xref ref-type="fig" rid="fig1">Figure 1</xref>0, the fundamental resonant frequency of the ZnO cantilevers assembled on aluminum, steel and silicon substrates was approximately 173.5 Hz, 170.0 Hz and 175 Hz, respectively, which is in the vicinity of the resonant frequency used in the experimental work of C.T. Pan et al. for ZnO cantilever on PET substrate [<xref ref-type="bibr" rid="scirp.91910-ref36">36</xref>] . Furthermore, ZnO cantilever assembled</p><p>on steel substrate yields 9.83 V at the resonant frequency. This value is better than the ones obtained using ZnO cantilever assembled on aluminum and silicon substrates.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>1 and <xref ref-type="fig" rid="fig1">Figure 1</xref>2 demonstrate the mechanical power input and the calculated electrical power output of the piezoelectric vibration energy harvester designed using ZnO cantilevers assembled on aluminum, steel and silicon substrates. <xref ref-type="fig" rid="fig1">Figure 1</xref>2 indicates that ZnO cantilevers assembled on steel substrates yields 4.02 of electrical power at resonant frequency 173.5 Hz larger than the electrical output power yielded using ZnO cantilevers assembled on aluminum and silicon substrates.</p></sec><sec id="s4_2_2"><title>4.2.2. Acceleration Response</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref>3 illustrates the calculated voltage of ZnO cantilevers assembled on aluminum, steel and silicon substrates as a function of host acceleration keeping the resonant frequency fixed at a predetermined value and using a constant value of load resistance of 12 kΩ for the three cases. As can be clearly seen from <xref ref-type="fig" rid="fig1">Figure 1</xref>3, The voltage increases linearly as the host acceleration is increased gradually from 0.25 g - 2 g. Apparently, ZnO cantilever assembled on steel substrate yields an optimum value of voltage of approximately 20 V at a resonant frequency of 170 Hz. Whereas ZnO cantilever assembled on aluminum and silicon substrates yields a voltage of 16.4 V, 10.4 V at a resonant frequency of 173.5 Hz, 175 Hz, respectively.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>4 and <xref ref-type="fig" rid="fig1">Figure 1</xref>5 illustrate the mechanical power input and the calculated electrical power output of the ZnO cantilever assembled on aluminum, steel, and silicon substrates. As can be seen from <xref ref-type="fig" rid="fig1">Figure 1</xref>5, the ZnO cantilever assembled on steel substrates yields the maximum electrical power value of approximately 41 mW at a resonant frequency of 170 Hz.</p></sec></sec><sec id="s4_3"><title>4.3. All Cantilever Materials with All Substrate Materials</title><p><xref ref-type="table" rid="table3">Table 3</xref> and <xref ref-type="table" rid="table4">Table 4</xref> summarize the voltage and electrical power output of the three cantilever materials (ZnO, LiNbO<sub>3</sub> and AlN) assembled on aluminum, steel and silicon substrates calculated at a predetermined fixed resonant frequency. The tables summarize our main findings described in details in the text.</p></sec></sec><sec id="s5"><title>5. Conclusions</title><p>In Summary, LiNbO<sub>3</sub>, AlN and zinc oxide (ZnO) cantilever materials assembled on aluminum, steel and silicon substrates have been investigated using COMSOL Multiphysics based on Finite Element Method (FEM).Voltage, mechanical power input and electrical power output versus frequency for each system were calculated and found to change linearly with the frequency.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> The voltage yielded at the predetermined fixed resonant frequency of ZnO, LiNbO<sub>3</sub>, AlN cantilevers assembled on aluminum, steel and silicon substrates</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >ZnO</th><th align="center" valign="middle" >LiNbO<sub>3</sub></th><th align="center" valign="middle" >AlN</th></tr></thead><tr><td align="center" valign="middle" >Aluminum</td><td align="center" valign="middle" >173.5 Hz, 8.20 V</td><td align="center" valign="middle" >187.5 Hz, 1.67 V</td><td align="center" valign="middle" >279.5 Hz, 3.80 V</td></tr><tr><td align="center" valign="middle" >Steel</td><td align="center" valign="middle" >170.0 Hz, 9.82 V</td><td align="center" valign="middle" >183.5 Hz, 1.49 V</td><td align="center" valign="middle" >272.0 Hz, 4.86 V</td></tr><tr><td align="center" valign="middle" >Silicon</td><td align="center" valign="middle" >175.0 Hz, 5.19 V</td><td align="center" valign="middle" >189.0 Hz, 1.00 V</td><td align="center" valign="middle" >281.0 Hz, 4.4 V</td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> The electrical power output yielded at the predetermined fixed resonant frequency of ZnO, LiNbO<sub>3</sub>, AlN cantilevers assembled on aluminum, steel and silicon substrates</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >ZnO</th><th align="center" valign="middle" >LiNbO<sub>3</sub></th><th align="center" valign="middle" >AlN</th></tr></thead><tr><td align="center" valign="middle" >Aluminum</td><td align="center" valign="middle" >173.5 Hz, 2.80 mW</td><td align="center" valign="middle" >187.5 Hz, 0.12 mW</td><td align="center" valign="middle" >279.5 Hz, 0.61 mW</td></tr><tr><td align="center" valign="middle" >Steel</td><td align="center" valign="middle" >170.0 Hz, 4.02 mW</td><td align="center" valign="middle" >183.5 Hz, 0.10 mW</td><td align="center" valign="middle" >272.0 Hz, 0.92 mW</td></tr><tr><td align="center" valign="middle" >Silicon</td><td align="center" valign="middle" >175.0 Hz, 1.12 mW</td><td align="center" valign="middle" >189.0 Hz, 0.04 mW</td><td align="center" valign="middle" >281.0 Hz, 0.81 mW</td></tr></tbody></table></table-wrap><p>The resonant frequency of LiNbO<sub>3</sub>, AlN and ZnO cantilevers are found to be 187.5 Hz, 279.5 Hz and 173.5 Hz, respectively. Interestingly, ZnO cantilever yields the maximum voltage and electrical power values of 8.2 V and 2.8 mW, respectively when examined at resonance frequency of 173.5 Hz for different host accelerations. We attributed this striking result to the dissymmetrical atomic configuration of wurtzite ZnO, in which Zn<sup>2+</sup> cations and O<sup>2−</sup> anions exchange their position under the influence of external load heterogeneously. Simply, the voltage and electrical power output of LiNbO<sub>3</sub>, AlN and ZnO cantilevers were found to increase linearly with the host acceleration.</p><p>Examination of electric properties of different cantilevers assembled on the three different substrates as a function of frequency indicates that the voltage, mechanical power input and electrical power output of all possible combinations of cantilever/substrate changes linearly with frequency. Among all the cantilever/substrate structures studied as a function of frequency, we found that the resonant frequency of ZnO cantilever assembled on aluminum, steel, and silicon substrates to be 173.5 Hz, 170 Hz and 175 Hz, respectively. When examined as a function of host acceleration, our results indicate that ZnO cantilever assembled on steel substrate at resonant frequency of 170 Hz attains maximum voltage and electrical power output of 9.83 V and 4.02 mW, respectively. We hope that our results on different cantilevers investigated in this study could improve the cantilevers applications ranging from aircraft design to architecture, medical diagnostics, nanoscale measurement systems, and forensics.</p></sec><sec id="s6"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s7"><title>Cite this paper</title><p>Alsaad, A.M., Ahmad, A.A., Al-Bataineh, Q.M., Daoud, N.S. and Khazaleh, M.H. (2019) Design and Analysis of MEMS Based Aluminum Nitride (AlN), Lithium Niobate (LiNbO<sub>3</sub>) and Zinc Oxide (ZnO) Cantilever with Different Substrate Materials for Piezoelectric Vibration Energy Harvesters Using COMSOL Multiphysics Software. Open Journal of Applied Sciences, 9, 181-197. https://doi.org/10.4236/ojapps.2019.94016</p></sec></body><back><ref-list><title>References</title><ref id="scirp.91910-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Sharma, S.R. and Pant, B. (2017) Design and Development of Guided Four Beam Cantilever Type MEMS Based Piezoelectric Energy Harvester. Microsystem Technologies, 23, 1751-1759. https://doi.org/10.1007/s00542-016-2940-1</mixed-citation></ref><ref id="scirp.91910-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Firoozy, P., Khadem, S.E. and Pourkiaee, S.M. (2017) Power Enhancement of Broadband Piezoelectric Energy Harvesting Using a Proof Mass and Nonlinearities in Curvature and Inertia. International Journal of Mechanical Sciences, 133, 227-239.  
https://doi.org/10.1016/j.ijmecsci.2017.08.048</mixed-citation></ref><ref id="scirp.91910-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Sun, S. and Peter, W. (2019) Modeling of a Horizontal Asymmetric U-Shaped Vibration-Based Piezoelectric Energy Harvester (U-VPEH). Mechanical Systems and Signal Processing, 114, 467-485. https://doi.org/10.1016/j.ymssp.2018.05.029</mixed-citation></ref><ref id="scirp.91910-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Wei, C. and Jing, X. (2017) A Comprehensive Review on Vibration Energy Harvesting: Modelling and Realization. Renewable and Sustainable Energy Reviews, 74, 1-18. https://doi.org/10.1016/j.rser.2017.01.073</mixed-citation></ref><ref id="scirp.91910-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Ahmed, R., Mir, F. and Banerjee, S. (2017) A Review on Energy Harvesting Approaches for Renewable Energies from Ambient Vibrations and Acoustic Waves Using Piezoelectricity. Smart Materials and Structures, 26, Article ID: 085031.  
https://doi.org/10.1088/1361-665X/aa7bfb</mixed-citation></ref><ref id="scirp.91910-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Toprak, A. and Tigli, O. (2018) Micron Scale Energy Harvesters Using Multiple Piezoelectric Polymer Layers. Sensors and Actuators A: Physical, 269, 412-418.  
https://doi.org/10.1016/j.sna.2017.11.035</mixed-citation></ref><ref id="scirp.91910-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Xiang, H.-J., Zhang, Z.-W., Shi, Z.-F. and Li, H. (2018) Reduced-Order Modeling of Piezoelectric Energy Harvesters with Nonlinear Circuits under Complex Conditions. Smart Materials and Structures, 2, Article ID: 045004.  
https://doi.org/10.1088/1361-665X/aaaf92</mixed-citation></ref><ref id="scirp.91910-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Tang, G., Yang, B., Liu, J.-Q., Xu, B., Zhu, H.-Y. and Yang, C.-S. (2014) Development of High Performance Piezoelectric d33 Mode MEMs Vibration Energy Harvester Based on PMN-PT Single Crystal Thick Film. Sensors and Actuators A: Physical, 205, 150-155. https://doi.org/10.1016/j.sna.2013.11.007</mixed-citation></ref><ref id="scirp.91910-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Elvin, N. and Erturk, A. (2013) Advances in Energy Harvesting Methods. Springer Science &amp; Business Media, Berlin. https://doi.org/10.1007/978-1-4614-5705-3</mixed-citation></ref><ref id="scirp.91910-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Du, S., Jia, Y., Do, C.D. and Seshia, A.A. (2016) An Efficient SSHI Interface with Increased Input Range for Piezoelectric Energy Harvesting under Variable Conditions. IEEE Journal of Solid-State Circuits, 51, 2729-2742.  
https://doi.org/10.1109/JSSC.2016.2594943</mixed-citation></ref><ref id="scirp.91910-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Han, M., Yuan, Q., Sun, X. and Zhang, H. (2014) Design and Fabrication of Integrated Magnetic MEMS Energy Harvester for Low Frequency Applications. Journal of Microelectromechanical Systems, 23, 204-212.  
https://doi.org/10.1109/JMEMS.2013.2267773</mixed-citation></ref><ref id="scirp.91910-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">K&amp;#246;rner, C., Bauerei&amp;#223;, A. and Attar, E. (2013) Fundamental Consolidation Mechanisms during Selective Beam Melting of Powders. Modelling and Simulation in Materials Science and Engineering, 21, Article ID: 085011.  
https://doi.org/10.1088/0965-0393/21/8/085011</mixed-citation></ref><ref id="scirp.91910-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Lee, K. and Fishwick, P.A. (2001) Building a Model for Real-Time Simulation. Future Generation Computer Systems, 17, 585-600.  
https://doi.org/10.1016/j.vacuum.2004.01.052</mixed-citation></ref><ref id="scirp.91910-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">http://www.comsol.com</mixed-citation></ref><ref id="scirp.91910-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Akiyama, M., Nagao, K., Ueno, N., Tateyama, H. and Yamada, T. (2004) Influence of Metal Electrodes on Crystal Orientation of Aluminum Nitride Thin Films. Vacuum, 74, 699-703. https://doi.org/10.1016/j.vacuum.2004.01.052</mixed-citation></ref><ref id="scirp.91910-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Shampa, M. (2014) Preparation of Undoped and Some Doped ZnO Thin Films by Silar and Their Characterization. Bardhaman, New York.</mixed-citation></ref><ref id="scirp.91910-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Mackwitz, P., Rüsing, M., Berth, G., Widhalm, A., Müller, K. and Zrenner, A. (2016) Periodic Domain Inversion in X-Cut Single-Crystal Lithium Niobate Thin Film. Applied Physics Letters, 108, Article ID: 152902.  
https://doi.org/10.1063/1.4946010</mixed-citation></ref><ref id="scirp.91910-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Elfrink, R., Kamel, T., Goedbloed, M., Matova, S., Hohlfeld, D., Van Andel, Y., et al. (2009) Vibration Energy Harvesting with Aluminum Nitride-Based Piezoelectric Devices. Journal of Micromechanics and Microengineering, 19, Article ID: 094005.  
https://doi.org/10.1088/0960-1317/19/9/094005</mixed-citation></ref><ref id="scirp.91910-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Zhao, X., Shang, Z., Luo, G. and Deng, L. (2015) A Vibration Energy Harvester Using AlN Piezoelectric Cantilever Array. Microelectronic Engineering, 142, 47-51.  
https://doi.org/10.1016/j.mee.2015.07.006</mixed-citation></ref><ref id="scirp.91910-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Battista, L., Mecozzi, L., Coppola, S., Vespini, V., Grilli, S. and Ferraro, P. (2014) Graphene and Carbon Black Nano-Composite Polymer Absorbers for a Pyro-Electric Solar Energy Harvesting Device Based on LiNbO3 Crystals. Applied Energy, 136, 357-362. https://doi.org/10.1016/j.apenergy.2014.09.035</mixed-citation></ref><ref id="scirp.91910-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Kumar, B. and Kim, S.-W. (2012) Energy Harvesting Based on Semiconducting Piezoelectric ZnO Nanostructures. Nano Energy, 1, 342-355.  
https://doi.org/10.1016/j.nanoen.2012.02.001</mixed-citation></ref><ref id="scirp.91910-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Song, J., Zhou, J. and Wang, Z.L. (2006) Piezoelectric and Semiconducting Coupled Power Generating Process of a Single ZnO Belt/Wire. A Technology for Harvesting Electricity from the Environment. Nano Letters, 6, 1656-1662.  
https://doi.org/10.1021/nl060820v</mixed-citation></ref><ref id="scirp.91910-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Mahmud, A., Khan, A.A., Voss, P., Das, T., Abdel-Rahman, E. and Ban, D. (2018) A High Performance and Consolidated Piezoelectric Energy Harvester Based on 1D/2D Hybrid Zinc Oxide Nanostructures. Advanced Materials Interfaces, 5, Article ID: 1801167. https://doi.org/10.1002/admi.201801167</mixed-citation></ref><ref id="scirp.91910-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Erturk, A. and Inman, D.J. (2008) On Mechanical Modeling of Cantilevered Piezoelectric Vibration Energy Harvesters. Journal of Intelligent Material Systems and Structures, 19, 1311-1325. https://doi.org/10.1177/1045389X07085639</mixed-citation></ref><ref id="scirp.91910-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Kumari, K. and Khanna, G. (2016) Design and Simulation of Array of Rectangular Micro Cantilevers Piezoelectric Energy Harvester. International Journal of Engineering Research and Applications, 6, 41-49.</mixed-citation></ref><ref id="scirp.91910-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Anton, S. and Sodano, H.A. (2007) A Review of Power Harvesting Using Piezoelectric Materials (2003-2006). Smart Materials and Structures, 16, R1-R21.  
https://doi.org/10.1088/0964-1726/16/3/R01</mixed-citation></ref><ref id="scirp.91910-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Wang, Z.L. and Song, J. (2006) Piezoelectric Nanogenerators Based on Zinc Oxide Nanowire Arrays. Science, 312, 242-246. https://doi.org/10.1126/science.1124005</mixed-citation></ref><ref id="scirp.91910-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Pan, Z.W., Dai, Z.R. and Wang, Z.L. (2001) Nanobelts of Semiconducting Oxides. Science, 291, 1947-1949. https://doi.org/10.1126/science.1058120</mixed-citation></ref><ref id="scirp.91910-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Hu, G., Ma, Y. and Wang, B. (2009) Mechanical Properties and Morphology of Nylon 11/Tetrapod-Shaped Zinc Oxide Whisker Composite. Materials Science and Engineering A, 504, 8-12. https://doi.org/10.1016/j.msea.2008.12.025</mixed-citation></ref><ref id="scirp.91910-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">https://www.memsnet.org/material/aluminumnitridealnbulk/</mixed-citation></ref><ref id="scirp.91910-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">https://www.korth.de/index.php/162/items/19.html</mixed-citation></ref><ref id="scirp.91910-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Nann, T. and Schneider, J. (2004) Origin of Permanent Electric Dipole Moments in Wurtzite Nanocrystals. Chemical Physics Letters, 384, 150-152.  
https://doi.org/10.1016/j.cplett.2003.12.017</mixed-citation></ref><ref id="scirp.91910-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Anota, E.C., Villanueva, M.S. and Cocoletzi, H.H. (2010) Electronic Properties of Group III: A Nitride Sheets by Molecular Simulation. Physica Status Solidi C, 7, 2252-2254. https://doi.org/10.1002/pssc.200983499</mixed-citation></ref><ref id="scirp.91910-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Maggard, P.A., Nault, T.S., Stern, C.L. and Poeppelmeier, K.R. (2003) Alignment of Acentric MoO3F33-Anions in a Polar Material: (Ag3MoO3F3)(Ag3MoO4) Cl. Journal of Solid State Chemistry, 175, 27-33.  
https://doi.org/10.1016/S0022-4596(03)00090-2</mixed-citation></ref><ref id="scirp.91910-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Reilly, E.K., Burghardt, F., Fain, R. and Wright, P. (2011) Powering a Wireless Sensor Node with a Vibration-Driven Piezoelectric Energy Harvester. Smart Materials and Structures, 20, Article ID: 125006.  
https://doi.org/10.1088/0964-1726/20/12/125006</mixed-citation></ref><ref id="scirp.91910-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Pan, C., Liu, Z., Chen, Y. and Liu, C. (2010) Design and Fabrication of Flexible Piezo-Microgenerator by Depositing ZnO Thin Films on PET Substrates. Sensors and Actuators A: Physical, 159, 96-104. https://doi.org/10.1016/j.sna.2010.02.023</mixed-citation></ref></ref-list></back></article>