hydrogen peroxide formation in the cathode

Size: px
Start display at page:

Download "hydrogen peroxide formation in the cathode"

Transcription

1 Loughborough Unversty Insttutonal epostory An mpedane model for analy of EIS of polymer eletrolyte fuel ells under hydrogen peroxde formaton n the athode Ths tem was submtted to Loughborough Unverstys Insttutonal epostory by the/an author. Ctaton: CU-MAN S. CEN. and GEEND P.S. 05. An mpedane model for analy of EIS of polymer eletrolyte fuel ells under hydrogen peroxde formaton n the athode. Journal of Eletroanalytal Chemstry 745 pp Metadata eord: Verson: Aepted for publaton Publsher: Elsever B.V. ghts: Ths work s made avalable aordng to the ondtons of the Creatve Commons Attrbuton-NonCommeral-NoDervatves 4.0 Internatonal (CC BY-NC-ND 4.0 lene. Full detals of ths lene are avalable at: Please te the publshed verson.

2 An Impedane Model for Analy of EIS of Polymer Eletrolyte Fuel Cells under ydrogen Peroxde Formaton n the Cathode Samuel Cruz-Manzo* u Chen and Paul Greenwood Department of Aeronautal and Automotve Engneerng Loughborough Unversty Leestershre LE 3TU Unted Kngdom *E-mal: s.ruz-manzo@hotmal.om Abstrat In ths study an mpedane model based on eletrohemal theory onsderng hydrogen peroxde formaton durng a two-step oxygen reduton reaton ( n polymer eletrolyte fuel ells (PEFCs has been developed. To valdate the theoretal treatment eletrohemal mpedane spetrosopy (EIS measurements were arred out n an open-athode 6 m /ar PEFC stak. The results show that ndutve loops at low frequenes of the mpedane spetra are attrbuted to mehansms related to hydrogen peroxde formaton durng. The results also demonstrate that the mehansms durng onsumpton of hydrogen peroxde to form water (seond-step n an be the domnatng proe for loes n the PEFC ompared to the mehansms durng oxygen onsumpton to form hydrogen peroxde (frst-step n. xygen transport lmtatons an be a result of hydrogen peroxde adsorbed onto the surfae of the eletrode whh redues the number of atve stes n the athode atalyst layer for oxygen to reat. Ths study ould support results from other ermental tehnques to dentfy hydrogen peroxde formaton durng the that lmt the performane of PEFCs. Key words: EIS hydrogen peroxde PEFC two-step oxygen reaton mpedane model. Introduton Eletrohemal mpedane spetrosopy (EIS s a powerful ermental tehnque that an be appled n-stu for dagno of polymer eletrolyte fuel ell (PEFC performane. The mpedane response resultng from EIS represents the eletrohemal mehansms of the PEFC n the frequeny doman. It an be shown n a Nyqust plot showng real and magnary omponents the magntude and phase of the mpedane as a funton of frequeny an also be plotted n a Bode plot. The use of eletral equvalent ruts wth the ermental EIS tehnque s a well-establshed methodology to haraterse proees n the PEFC. An eletral equvalent rut an represent an dental mpedane response to that obtaned from the PEFC studed. Eah eletral

3 omponent n the eletral equvalent rut an desrbe a physal proe that takes plae n the PEFC. In the authors prevous study [] an eletral equvalent rut based on eletrohemal theory was developed to smulate and haraterse the frequeny response of the PEFC at any urrent of the polarsaton urve. owever the eletral rut had lmtatons n reprodung EIS measurements wth postve magnary omponents at low frequenes. EIS measurements wth postve magnary omponents at low frequenes are known as an ndutve loop and have been reported n EIS measurements of PEFCs [3]. Makhara et al. [3] suggested that the poble reasons for ths ndutve effet at low frequenes are the sde reaton and ntermedates nvolved n fuel ell reatons. oy et al. [4] developed mpedane models to aount for the reaton mehansms that may be responsble for the ndutve response at low frequenes; the models propose the formaton of hydrogen peroxde and platnum doluton durng the oxygen reduton reaton (. It has been reported [5] that roover of hydrogen to the athode faltates the reaton of oxygen and hydrogen at the athode generatng hydroxyl and hydroperoxyl radals whh reat further to produe hydrogen peroxde at the athode. The hypothe that hydrogen peroxde may be formed at the athode of a fuel ell s supported by the results of Inaba et al. [6]. In ths study an mpedane model onsderng fundamental eletrode theory and onsderng hydrogen peroxde formaton n a two-step has been developed. Modellng onsderatons based on EIS theory were taken nto aount to smplfy the mathematal treatment the hgh frequeny regon related to on dffuson mehansms n porous eletrodes found n the prevous study [] were not onsdered here. The fous of the urrent study s at the low frequeny range where the mehansms relatng to hydrogen peroxde formaton durng the oxygen reduton are represented. The mpedane model onsderng formaton of hydrogen peroxde durng a two-step was valdated aganst EIS measurements arred out n an openathode 6 m PEFC stak.. Eletrohemal mehansms of hydrogen peroxde formaton durng The two-step whh yelds the formaton of hydrogen peroxde ( and the formaton of water ( has been reported n the lterature [47] as: k e ( k e ( The urrent denstes orrespondng to reatons ( and ( have been reported n the lterature [8] as: ( ( b V nfk 3/ (3

4 ( b V nfk 3/ (4 where F s the Faraday onstant k s the rate onstant whh depends on the global order of reaton s the onentraton at the eletrode surfae n s the number of eletrons onsumed durng the reaton fratonal surfae overage of hydrogen peroxde b s the apparent Tafel slope and V s the eletrode potental. The followng aumptons are onsdered to faltate the mathematal treatment: s the. Ths study gnores hydrogen peroxde dffuson n the CCL durng EIS measurements as a frst approxmaton. Inaba et al. [6] reported that when Pt/C atalysts are hghly agglomerated the moleules easly dffuse to nearby Pt/C partles where they are transformed to water. To aount for hydrogen peroxde dffuson durng EIS measurements the effet of partle sze nteratons between Pt partles and the substrate arbon on Pt/C agglomeraton should be studed and onsdered n the modellng work. Ths wll be addreed n future work.. Sethuraman et al. [9] reported that the rate of hydrogen peroxde formaton durng the s nreased wth nreasng proton onentraton. Proton onentraton n the CCL s a funton of the amount of onomer or eletrolyte n the maropore spaes between agglomerates [0]. In ths study t has been onsdered that there s no effet to the opposton of flow of protons n the CCL on hydrogen peroxde formaton. Proton retane n the CCL s negleted. Ths smplfed the mathematal treatment by onsderng that the onentraton of protons aro the CCL s onstant from the CCL- PEM nterfae to the Pt/C atalyst agglomerates. 3. There are three modes of transport of oxygen n the CCL whh an have an effet on ma transport lmtatons: gas-dffuson n the eletrode-pore dolved oxygen dffuson n the lqud water surroundng the agglomerate and dolved oxygen n the onomer phase []. The fnte dffuson dstane for oxygen to reah the reaton stes n the CCL forms a omplated network of mult-phase parallel and seral paths. EIS only measures bulk parameters n the total CCL thkne and reflets a total ma transport retane and a total fnte dffuson dstane for the three modes of oxygen transport n the CCL therefore dffuson of oxygen n the CCL wll be onsdered from the gas dffuson layer (GDL-CCL nterfae to the CCL-PEM nterfae. Ths wll smplfy the mathematal analy and therefore the resultng model wll present parameters ommonly known n the

5 eletrohemal area suh as the arburg mpedane [] whh allows the haratersaton of oxygen transport lmtatons durng EIS measurements. Takng nto aount the frst and seond aumptons prevously defned Eqs. 3 and 4 an be defned as: ( ( b V 3 / nfk (5 ( b V 3/ nfk eq (6 where eq represents the onentraton at equlbrum ondtons as t s aumed that surfae onentraton of protons and hydrogen peroxde s onstant. The rate at whh the transfer of harge ours at equlbrum ondtons s known as the exhange urrent. The exhange urrent denstes n a two-step eletron transfer reaton are desrbed n Appendx A. The exhange urrent denstes for reatons ( and ( at equlbrum ondtons wth ( ( b V 3/ nfk an be defned as: V V eq 0 (7 ( b V 0 nfk (8 3/ where 0 represents the exhange urrent densty s the fratonal surfae overage of hydrogen peroxde under equlbrum ondtons and V eq s the equlbrum potental. Multplyng both terms of Eq. 5 by ( ( b V ( b V eq ( ( and substtutng Eqs. 7 and 8 yelds: [ b ( V V ] 0 and smlarly Eq. 6 by (9 [ b ( V V ] 0 (0 Aumng that the formaton of hydrogen peroxde s governed by the Langmur sothermal t s poble to defne the total urrent and the varaton rate of as [3]: (

6 F d β ( dt where β s the maxmum surfae overage by hydrogen peroxde. Substtutng Eqs. 9 and 0 nto Eqs. and yelds: ( ( [ b ( V V ] [ b ( V V ] 0 (3 F d 0 ( ( [ b ( V V ] [ b ( V V ] b 0 0 (4 eq dt. a mpedane response If a small a perturbaton s appled to the workng eletrode the EIS tehnque allows the use of a lnear equaton to smulate mpedane spetra. A lnear model an be derved usng the Taylor seres anson as reed below: V V (5 To faltate the mathematal treatment the small dsturbane of oxygen from the a perturbaton s not onsdered n Eq. 5. Based on ermental EIS measurements n whh physal parameters are measured aro the entre CCL the hange n oxygen onentraton durng the wll be onsdered from the PEM- CCL nterfae to the CCL-GDL nterfae. The term that defnes the hange of oxygen onentraton from the GDL-CCL nterfae to the PEM-CCL nterfae at the frequeny doman has already been developed n the authors prevous study [] and wll be onsdered n the next seton. It wll be poble to norporate the arburg Impedane whh has been broadly used n the analy of oxygen transport lmtatons n ermental EIS measurements of PEFCs nto the mathematal treatment. The omplex value of an be related to V f Eq. s lnearzed by the followng equaton [4]: Fβ jω V V (6 where ω s the angular frequeny and j s the magnary omponent.

7 Applyng Eqs. 3 and 4 nto Eqs. 5 and 6 and onsderng the oxygen onentraton at the PEM-CCL nterfae yelds: 0 b Fb ( ( 0 b jω ( b η 0 ( b η ( ( b V 0 ( b η ( b η b 0 V ( 0 ( 0 ( b η ( b η (7 0 ( b η ( b η where the subsrpt represents the steady-state η V s the overpotental at steady-state oxygen onentraton at the CCL-PEM nterfae and V (8 s the s the onentraton at the GDL-CCL nterfae. The eletrode surfae overage of hydrogen peroxde at steady-state an be alulated from Eq. 4 onsderng d / dt 0 and 0 ( ( b η eq 0 0 ( ( b ( η b η (9 The eletrode surfae overage of hydrogen peroxde at steady-state represented n Eq. 9 an be reed n terms of steady-state urrents ; ths an be demonstrated n Appendx B. Epelbon and Keddam [5] defned the overage oeffent of adsorbed ntermedate spees at steady-state n a two-step reaton as the rato between the partal urrent of the frst reaton and the sum of the urrent of the frst and seond reatons. The hange n oxygen onentraton aro the CCL resultng from the small a dsturbane of the surfae overage of hydrogen peroxde formaton n the frequeny doman s onsdered neglgble. Therefore the rato of oxygen onentraton between the CCL-PEM and GDL-CCL nterfaes / eq whh are multpled by n Eqs. 7 and 8 (thrd term on the rght hand sde s onsdered to be steady-state. Under suh onsderatons t s poble to defne Eqs. 7 and 8 n terms of harge transfer retanes and steady-state urrents: V [ ] (0

8 Fβ jω V [ ] ( where represents the harge transfer retane: 0 b ( b η ( b 0 ( b η (3 and represents the steady-state urrent: 0 0 ( ( b η ( b η (45. xygen transport lmtatons The oxygen onentraton at the eletrode surfae depends on the dffuson of oxygen through the multphase va parallel and seral paths n the CCL [6]. In the authors prevous study [] the rato between oxygen onentratons at the CCL-PEM and CCL-GDL nterfaes n the frequeny doman was derved to haraterse oxygen transport lmtatons n the mpedane response of PEFCs. Ths rato between oxygen onentratons s onsdered as the EIS tehnque measures the hange n oxygen onentraton n the total CCL. Ths also allows the dervaton of the arburg Impedane whh has been broadly used n the EIS area [] to haraterse oxygen transport lmtatons n the frequeny doman durng PEFC operaton. The rato of oxygen onentraton at the PEM-CCL and GDL-CCL nterfae at the frequeny doman s reed as []: (6 where and s the harge transfer retane for the oxygen reduton whh an be represented through Eq. ( jωt ( jωt tanh (7 s known as the arburg Impedane [] and haraterses oxygen transport lmtatons n EIS measurements of PEFCs wth Tδ (8 z F D *

9 defned as the retane of the dffuson proe and D T δ (9 defned as the tme onstant to dffuse oxygen through the CCL. 3. Eletrohemal Impedane of the PEFC under ydrogen Peroxde Formaton Substtutng Eq. 6 nto Eqs. 0 and gves: [ ] V (30 [ ] V j F ω β (3 Arrangng Eq. 3 for yelds: ( ( V j F ω β (3 Substtutng Eq. 3 nto Eq. 30 yelds: ( ( ( ( V j F ω β (33 The CCL s formed by a double layer struture of the nterfae between dmlar materals.e. the eletrodeeletrolyte nterfae. The eletrode represents the eletron-ondutng network of arbon-supported platnum agglomerates whle the eletrolyte represents the on-ondutng dsperson of Nafon. At ths nterfae an eletr feld s present and has a determnant role n the harge dstrbuton wthn the reatants as well as n the poston and orentaton of the reatant to form the desred produts. Ths double layer an behave lke a apator that s n parallel wth the eletrode reatons; the urrent pang from the eletrode to the eletrolyte an ether take part n the harge transfer reatons or ontrbute to the harge n the apatve effet. The eletrohemal reaton n the eletrode results n an nhomogeneous dstrbuton of harge n the atalyst layer. As a result a non-deal apator (Constant Phase Element [7] defned n the frequeny doman jω has to be onsdered n the theoretal treatment n order to orret for ths nhomogenety:

10 ( ( ( ( ( V j Y V j F P ω ω β (34 where Y represents a parameter related to apatane and supersrpt P represents a parameter to orret the nhomogenety n the dstrbuton of harge. Eq. 34 an be rearranged to re the eletral urrent that enters a rut omposed of retve and apatve elements as shown n Fg. gvng: ( V j Y A P ω φ (35 where ( φ ω τ φ j A F β τ ( ( s I φ and ( ( s I haraterse ndutve loops at low frequenes n the mpedane response of PEFCs. earrangng Eq. 35 yelds the eletrohemal mpedane of the athode aountng for hydrogen peroxde formaton durng the : ( P j Y A ω φ (36 The hydrogen oxdaton reaton ( mehansms ourrng at the anode have been proposed by Malevh et al. [] to be represented as an eletral rut omposed of a retor (harge transfer retane n parallel wth a apator (representng the double layer apatane. takng plae n the anode s a faster and a le omplated reaton sequene than the n the athode and an be represented at the frequeny doman as: aca j a a ω (37 f no ontamnants suh as C are present where a s harge transfer retane and Ca s apatane.

11 The opposton of ons flowng through the PEM and eletrons flowng through the bpolar plate and GDL an be represented by retanes onneted n seres. The overall proe an be represented wth a total retane ohm that aounts for the dfferent retanes n the layers of the PEFC. Therefore the mpedane of the PEFC nludng anode athode and PEM an be represented through the followng equaton: PEFC ohm a (38 where ohm represents the ohm retane of the PEM bpolar plates GDL the seond term on the rght hand sde represents Eq. 37 and the thrd term represents the athode mpedane (Eq. 36. Eq. 38 an be represented through the equvalent eletral rut shown n Fg. For the spef ase when the harge transfer retane related to hydrogen peroxde reaton durng the seond-step s smaller than the sum of harge transfer retane and ma transport retane durng the frst-step ; n Eq. 35 beomes negatve and a apatve loop at the low φ frequeny end of the spetrum s predted nstead of an ndutve loop as shown n Fg Expermental A ommeral pen-cathode 4-ell stak wth a 6 m atve area was used for the ermental tests. The membrane eletrode aembles (MEAs were made of Nafon wth platnum loadngs of 0.4mg/m and arbon blak for the eletrodes. The GDLs were made of arbon felt wth 00 µm wdth. The 4-ell stak onts of open athodes wth two 5V DC fans for oxdant supply (as ar and oolng. The stak onts of bpolar plates made of FU4369 T materal wth a thkne of 5 mm. gh purty hydrogen ( % was used durng the tests to redue fuel mpurtes whh an lmt the performane of the PEFC. pen-cathode fuel ells are ommonly operated wth a dead-ended hydrogen outlet; however under suh a onfguraton a nonsteady response s eted at medum and hgh urrents due to buld-up of water. In ths study the fuel ell stak was run n a through flow mode at the anode. The flow rate of hydrogen n the anode was kept onstant durng the erments wth a stohometry of. The hydrogen suppled was dry. The PEFC stak was operated at an ambent temperature of C and the hydrogen bak preure was held at 0.4 bar(g. A polarsaton urve was reorded pror to mpedane measurements as shown n Fg. 4.

12 EIS measurements were arred out through a #06 onateh frequeny response analyser (FA. EIS measurements were arred out at three dfferent urrent denstes and A/m of the polarsaton urve as shown n Fg. 4. The PEFC stak was run over 60 mnutes at the requred urrent wth no varaton n the voltage to ensure a steady state for EIS measurements. EIS measurements for urrent denstes < 0.875A/m were not poble as the low AC ampltude supermposed onto the DC urrent made t dffult for the FA to dstngush between nose and response. At hgh urrent denstes > A/m the fuel ell stak was not steady for a long perod of tme due to the hgh water onentraton produed by the. The frequeny san was performed from 0 kz down to 0. z wth a 5% AC ampltude of the DC urrent. It has been ommonly onsdered that the oxygen n the fuel ell an be eletrohemally redued dretly to water wthout the formaton of ntermedate spees. The hypothe that a -step s present durng fuel ell operaton s supported by the results from Inaba et al. [6]. The authors revealed that hydrogen peroxde an be formed by a two-eletron reduton pathway and exsts on Pt partles supported on arbon. After hydrogen peroxde s desorbed from the Pt surfae and dffuses away to the bulk of the soluton. If Pt/C atalysts are hghly agglomerated the moleules easly dffuse to nearby Pt/C partles where they are transformed to water by further two-eletron reduton (seres four-eletron reduton or -step. 5. Expermental Valdaton and esults EIS measurements are ommonly analysed by fttng an mpedane model wth ermental data. An nrease n the number of parameters n an equaton that represents the PEFC mpedane response may lead to nreased error n the resultng ftted values wth ermental EIS. A andles rut as reported by Malevh et al. [] was frst ftted to the measured EIS spetra usng Vew.9 software (Srbner Aoates In as shown n Fg. 5. The andles rut was ftted to the EIS measurements through a omplex nonlnear regreon method (Levenberg-Marquardt. Ths allowed the estmaton of Y P ohm and the parameters related to a n Eq. 38. Notng that the andles rut shown n Fg. 5a represents the eletral rut shown n Fg. wthout the parameters. and A φ whh are related to hydrogen peroxde reaton mehansms durng the

13 φ n Eq. 38 s a frequeny dependent parameter therefore t annot be reed as a retor n a ommeral software (e.g. vew. The parameters τ and the ones related to φ from Eq. 38 were estmated usng a Graph User Interfae (GUI developed n Matlab. The use of the GUI for EIS analy has already been demonstrated n the authors prevous study [8]. The parameters τ and from φ were ftted to the EIS measurements usng the GUI to reprodue the ndutve loop at low frequenes as shown n Fg. 6. The GUI allows the fttng of the parameters from Eq. 38 to aheve a good agreement between the ermental and smulated data. The parameter was realulated to aount for the harge transfer retane of oxygen onsumpton durng a two-step. The least squares fttng method was used n order to fnd the best-ft between the model and the measured data. A good qualty ft s obtaned when the sum of the devatons squared (least square error between the smulated and measured mpedane data has a mnmum value for nstane <0.. The parameters of Eq. 38 extrated from the measured EIS data are shown n Table I. Fgure 6 shows that wth the parameters gven n Table I the theoretal model s apable of smulatng the mpedane response of the PEFC stak. The results show good agreement between the measured and smulated data n the Nyqust plot. oy et al. [4] developed an mpedane model to aount for formaton of hydrogen peroxde durng the. The model was ompared wth ermental results and the results show that hydrogen peroxde formaton mehansms ould aount for the ermentally observed low-frequeny ndutve loops n EIS measurements. The omparson between smulated and measured data n the urrent study shows a better agreement than the omparson between smulated and measured data reported by oy et al. [4]. ne of the reasons of fndng a better agreement n ths study s the fat that a onstant phase element whh aounts for nonhomogeneous harge dstrbuton n the atalyst layer was onsdered n the modellng approah. Dfferent to the modellng approah by oy et al. [4] where an deal apator whh represents homogenous harge dstrbuton n the atalyst layer was onsdered. Fg. 6 shows that the dameter of the spetrum dereases from a urrent densty of m to 0.35 A/m and nreases from a urrent densty of 0.35 A/m to A/m. Smlar results were reported by Yuan et al. [9] when EIS measurements were arred out n a 500 sx-ell PEFC stak. Ther results showed that the mpedane spetrum dereased wth nreasng urrent up to 80 amperes and beyond 80 amperes the mpedane spetrum nreased agan. Pagann et al. [0] reported that when the knets of the domnates the ell performane suh as n the low urrent densty range of the polarsaton urve the mpedane spetrum manly

14 represents the harge transfer effet durng the and ts dameter dereases wth nreasng urrent densty. At hgh urrents the nrease n dameter of the spetra wth nreasng urrent densty an be attrbuted to an nrease n oxygen transport lmtatons [9]. ne of the dsadvantages of the EIS tehnque s that multple energy ontrolled proees durng the eletrohemal reaton an be masked n the mpedane plot. The harge transfer retanes related to oxygen reduton and hydrogen peroxde reaton derease when urrent densty s nreased from m to 0.35 A/m as shown n Tab.. The oxygen harge transfer retane s smaller than the hydrogen peroxde harge transfer retane. Tafel slopes related to a two-step wth hydrogen peroxde formaton have been reported n the lterature. The apparent Tafel slope related to oxygen onsumpton to form hydrogen peroxde has been reported to be greater than the apparent Tafel slope related to the hydrogen peroxde onsumpton to form water [4]. A hgh value of apparent Tafel slope (V - dereases the value of the harge transfer retane. Tab. shows that the harge transfer retane of hydrogen peroxde reaton nreases wth nreasng urrent densty from 0.35 A/m to A/m. Ths nrease an be a result due to the fat that Eq. 38 does not take nto aount the varaton of hydrogen peroxde onentraton at the reaton stes. In ths study the surfae onentraton of hydrogen peroxde was onsdered onstant to faltate the mathematal treatment. The nrease of hydrogen peroxde harge transfer retane an be related to a hange of the value of b Tafel slope wth nreasng urrent densty when dffuson effet s present. Further work s requred to develop an analytal reon that onsders dffuson of hydrogen peroxde n EIS measurements as a funton of ell potental. A further valdaton of the mathematal model aganst the modulus and phase angle of the frequeny response of the measured data was arred out. Fgs. 7 8 show the omparson between measured and smulated data n Bode format. Fg. 7 shows that the model s able to reprodue the Bode modulus for all the frequenes. Fg. 8 shows a dsrepany n the phase angle between measured and smulated data at the hghest frequenes ( > 5 kz. razem and Trbollet [] reported n ther study that the Bode modulus and real part omponent of the mpedane plots are relatvely nsenstve to the qualty of the ft of a model to mpedane data. The magnary omponent of the mpedane and Bode phase angle plots are modestly senstve to ft qualty. The on retane n the CCL has been represented n EIS measurements [3] as a 45 o regon at hgh frequenes. At hgh frequenes (Fg. 6 t s learly shown that the 45 straght lne represents the on retane of the atalyst layer. The dsrepany between the measured data and smulated data n phase angle at hgh frequenes

15 (5-0 kz shown n Fg. 8 an be related to the fat that the on retane n the CCL was not onsdered n Eq. 38. A property from EIS measurements s that as the frequeny reahes a value of zero the operatng ondtons of the fuel ell must onverge to steady state ondtons [3]. Ths means that a retve value (no frequeny dependene for the mpedane response s eted as the frequeny reahes zero. Ths retve value s named polarsaton retane and s related to the slope of the E-I (polarsaton urve [4]. oy et al. [4] reported that the presene of sde reatons (e.g. hydrogen peroxde formaton has lttle nfluene on polarsaton urves as these reatons take plae a low rate ompared to and. Also the authors onluded that EIS measurements are more senstve than polarsaton urves wth the sde reatons present. A further study for EIS measurements n a PEFC wth low rate of sde reatons (no ndutve loop and ndung sde reatons (ndutve loop should be arred out. Ths s to study the relaton between polarsaton retane and the mpedane response at frequenes near 0 z for both ondtons. 6. Dsuon The eletral urrent durng a smple eletron-transfer takes part n the harge transfer and smultaneously takes part n the onsumpton of oxygen durng the. Ths proe an be represented through the andles rut where the farada mpedane onts of the harge transfer retane onneted n seres wth the ma transport retane (arburg element []. Durng a two-step eletron transfer the sum of the reatons hydrogen peroxde formaton and water formaton n the two-step yelds the reaton aountng for formaton of water wthout any ntermedates. The resultng eletral rut that aounts for formaton of hydrogen peroxde durng the an be represented through a parallel eletr rut onfguraton between the two-step mehansms as shown n Fg.. In the study of oy et al. [4] eletral ruts were developed to aount for adsorbed spees durng a two-step n whh the mpedane related to the frst step of the s n parallel wth the mpedane related to the seond step of the. agner and Gülzow [5] reported an eletral rut onfguraton to aount for a PEFC durng posonng of the anode wth C n whh the anode harge transfer retane s onneted n parallel wth a relaxaton mpedane desrbng the surfae relaxaton of the eletrode-eletrolyte nterfae. Ther study reported ndutve loops at low frequenes related to C posonng n PEFC eletrodes. gh purty hydrogen ( % was used durng these ermental tests therefore the nfluene of C on EIS measurements an be ruled out n ths work.

16 The harge transfer retane for anode a was one order of magntude le than the harge transfer retanes and related to as shown n Table. Ths demonstrates that the takng plae n the anode s a faster and a le omplated reaton sequene than the n the athode. Commonly t has been aumed durng the analy of EIS measurements that there s no ontrbuton of the anode mpedane to the ell mpedane []. owever t s hghly reommended to onsder the anode mpedane durng the fttng proe of EIS measurements f a referene eletrode s not used wthn the ermental system. Ths would avod unertanly n the resultng data. The value of harge transfer retane ourrng durng the frst step n the (oxygen onsumpton was smaller than the value of the harge transfer retane durng the seond step n the (hydrogen peroxde onsumpton. ydrogen roover from anode to athode has been reported as the mehansm responsble for the formaton of hydrogen peroxde [4]. The model predted that when the harge transfer retane related to hydrogen peroxde reaton the harge transfer retane and ma transport retane durng the seond-step s greater than the sum of durng the frst-step an ndutve loop an be seen at low frequenes of the mpedane spetrum. Ths relaton > was valdated wth ermental measurements as shown n Tab and therefore demonstrates that the knet mehansms durng onsumpton of hydrogen peroxde to form water an be the domnatng proe for loes n the PEFC ompared wth the mehansms durng oxygen reduton. ne of the dsadvantages of the EIS tehnque s that the tme onstants related to the dfferent eletrohemal proees are overlapped n the omplex-mpedane plot [6]. The nrease n oxygen onentraton lmtatons durng the frst step of the an be a result of the ntermedate spees adsorbed onto the eletrode surfae whh redue the number of atve stes n the CCL for oxygen to reat. From Table I t an be demonstrated that the tme onstant proe τ whh relates to the reaton of hydrogen peroxde durng the seond step of the s greater than the tme onstant for oxygen to dffuse through the CCL T. Ths an be demonstrated n the EIS omplex plot where the ndutve loop related to hydrogen peroxde formaton mehansms are present at the lowest frequenes. The proe responsble for the ndutve haratersts at low frequenes has been attrbuted to hydrogen roover through the PEM [5] to form hydrogen peroxde at the athode. ydrogen roover from anode to

17 athode through onventonal graphte for bpolar plates due to ther porous nature ould also be a responsble proe for hydrogen peroxde formaton and ndutve loops n EIS measurements. Fg. 6 shows a slght dfferene n the shape of the ndutve loop from ermental data ompared to the smulated data. Ths dfferene ould reveal that durng fuel ell operaton not only hydrogen peroxde formaton s present but also other eletrohemal mehansms that an ause fuel ell degradaton. It has been proposed [4] that platnum doluton ould also be represented as an ndutve loop at low frequenes. Therefore the results shown n Fg. 6 an lead to the development of a future mpedane model whh ombnes hydrogen peroxde formaton wth platnum doluton. 7. Conlusons An mpedane model based on hydrogen peroxde formaton durng the has been developed to haraterse EIS measurements of PEFCs featurng ndutve loops at low frequenes. The urrent study begns by defnng the eletrohemal equatons representng the rate of harge transfer durng hydrogen peroxde formaton n a two-step. Modellng onsderatons based on EIS theory were taken nto aount to smplfy the mathematal treatment. The equatons defnng the eletrohemal mehansms durng the two-step were transformed from the tme doman to the frequeny doman to smulate the frequeny response of the PEFC and thus the low frequeny ndutve loop. The theoretal treatment s valdated aganst the measured EIS response of an open-athode PEFC stak. The results show that the model an predt the ermental mpedane response presentng an ndutve loop at low frequenes attrbuted to mehansms of hydrogen peroxde reaton. The results also show that eletrohemal mehansms related to the reaton of hydrogen peroxde durng the an lmt the PEFC performane. Ths study ould ast wth other ermental studes to dentfy hydrogen peroxde formaton durng the whh lmt the performane of PEFCs. Aknowledgements The authors thank the Mexan Natonal Counl for Sene and Tehnology (CNACYT for the sponsorshp of the Ph.D researh study of S. Cruz-Manzo (Grant No Appendx A Dard et al. [3] desrbed some examples of the mehansms represented by a two-step eletron transfer reaton wth an adsorbed ntermedate spees suh as athod produton of hydrogen:

18 e (A ads ads e (A eatons A and A an be represented as: K rea tan t e adsorbate (A3 Kr ( rea tant adsorbate e produt K (A4 K r The net urrent for reatons (A3 and (A4 an be reed as: { K o θ K ( θ } FΓ (A5 r { K o ( θ K θ } FΓ (A6 r wth K k ( α FE T and K k ( α FE T oj oj oj / j rj rj rj / where F s the Faraday onstant Γ s the number of eletrosorpton stes per unt area of eletrode surfae and θ s the overage of the eletrode surfae by the adsorbate. The onentraton of the hemal spees s onsdered onstant and s ntrodued nto the rate onstants. At equlbrum ondtons 0 the exhange urrent denstes are defned as: * o eq * r ( q eq 0 FΓK q FΓK (A7 * o * ( qeq F Krq eq 0 FΓK Γ (A8 wth K * k ( α FE T and K k ( α FE T * oj oj oj eq / j rj rj rj eq / Appendx B The steady-state value s alulated from Eq. 4 wth d / dt 0 and gvng: ( ( 0 η 0 ( b η ( b 0 (B Arrangng Eq. B for gves: 0 ( ( b η 0 0 ( ( b η ( b η (B The total urrent for reatons ( and ( n steady-state onsderng Eq. 3 and an be reed as:

19 ( ( (B3 0 ( b η ( b η 0 Substtutng Eq. B nto Eq. B3 yelds: 0 ( 0 ( ( b η 0 ( b η ( b η 0 ( b η (B4 Eqs. (B and (B4 an be reed n terms of steady-state urrents: (B5 (B6 wth 0 0 ( ( b η ( b η Epelbon and Keddam [5] reported a smlar relaton for the total urrent and fratonal surfae overage of hydrogen peroxde at steady-state. Lst of Symbols b apparent tafel slope V Ca double layer apatane n anode F/m oxygen onentraton at the eletrode surfae mol/m 3 eq oxygen onentraton at equlbrum mol/m 3 oxygen onentraton at the CCL-PEM nterfae mol/m 3 * oxygen onentraton at equlbrum ondtons mol/m 3 D F j effetve dffuson oeffent m /se faraday onstant C/mol magnary omponent n mpedane urrent densty A/m 0 exhange urrent densty A/m steady-state urrent densty A/m

20 k P devaton of varable n the steady-state. A/m rate onstant of reaton parameter related to apatane (onstant phase element a harge transfer retane durng Ω.m harge transfer retane Ω.m ohm total ohm retane PEM GDL Plate Ω.m retane for the dffuson proe Ω.m s T t ω laplae doman tme onstant for the dffuson proe se tme se angular frequeny rad/se Y parameter related to apatane ( a mpedane of anode Ω.m mpedane of athode Ω.m mpedane of PEFC Ω.m PEFC arburg mpedane Ω.m se P / Ω. m real part of mpedane magnary part of mpedane Greek β maxmum surfae overage by hydrogen peroxde mol / m fratonal surfae overage of hydrogen peroxde dmensonle fratonal surfae overage of hydrogen peroxde at equlbrum ondtons dmensonle surfae overage at steady-state dmensonle devaton of varable n the steady-state dmensonle δ η fnte dffuson dstane of oxygen transport from CCL-GDL to CCL-PEM nterfae m overpotental n the CCL V

21 η devaton of varable η n the steady-state η V η steady-state overpotental V eferenes [] S. Cruz-Manzo and. Chen J. Eletroanal. Chem (03. [] N. V. Dale M. D. Mann. Salehfar A. M. Dhrde and T. an ASME J. Fuel Cell 7 ( /-0300/0 [3]. Makhara M. F. Mathas and D.. Baker J. Eletrohem. So. 5 A970 (005. [4] S. K. oy M. E. razem and B. Trbollet J. Eletrohem. So. 54 B378 (007. [5] J. Xe D. L. ood III K. L. More P. Atanaov and. L. Borup J. Eletrohem. So. 5 A0 (005. [6] M. Inaba. Yamada J. Tokunaga and A. Tasaka Eletrohem. Sold-State Lett. 7 A474 (004. [7] A. Damjanov and V. Brus Eletrohm. Ata 65 (967. [8] S. K. oy and M. E. razem ECS Trans. 3( 03 (006. [9] Vjay A. Sethuraman John. edner Andrew T. aug Sathya Motupally and Lesa V. Protsalo J. Eletrohem. So. 55 (008 B50-B57 [0] M. Ekerlng A.A. Kornyshev J. Eletroanal. Chem. 453 ( [] D. Malevh E. allop B. A. Peppley J. G. Pharoah K. Karan J. Eletrohem. So. 56 (009 B6-B4. [] N. Fouquet C. Doulet C. Noullant G. Dauphn-Tanguy B. uld-bouamama J. Power Soures 59 ( [3] J.-P. Dard B. Le Gorre and C. Montella J. Eletroanal. Chem (99. [4] I. Epelbon M. Keddam and J. C. Lestrade Faraday Dsu. Chem. So (973. [5] I. Epelbon and M. Keddam J. Eletrohem. So (970. [6] Ekerlng M. and Kornyshev A. A. J Eletroanal. Chem. 475 (999 pp [7] C.. su and F. Mansfeld Corroson 57 ( [8] S. Cruz-Manzo and. Chen J. Eletrohem. So. 60 (03 F09- F5 [9] Yuan X. Sun J. C. Blano S. ang. hang J. and lknson D. P. J. Power Soures 6 (006 pp [0] Pagann V. A. lvera C. L. F. Tanell E. A. Sprnger T. E. and Gonzalez E.. Eletrohm. Ata 43 ( [] razem M. and Trbollet B. Eletrohemal Impedane Spetrosopy ley New Jersey (008. [] S. Cruz-Manzo. Chen and P. ama Journal of Fuel Cell Sene and Tehnology 9 ( [3] G. Maranzana J. Manka. Lottn J. Dllet A. Lambra A. Thomas and S.Dderjean Eletrohm. Ata 83 ( [4] N. agner J. of Appl. Eletrohem. 3 ( [5] N. agner and E. Gülzow J. Power Soures 7 34 (004. [6] S.Cruz-Manzo. Chen and P. ama Int. J. ydrogen Energy 38 (

22 Fgures Fgure. Eletral rut representng eletrohemal mehansms n the reaton wth hydrogen peroxde formaton Fgure. Eletral rut representng the PEFC Fgure 3. Effet of harge transfer retane of hydrogen peroxde reaton on mpedane spetrum a > b <

23 Fgure 4 Polarsaton urve for open-athode stak Fg. 5 Comparson between ermental and smulated data a andles rut b fttng result Fgure 6. Comparson between ermental and smulated data from Eq. 38

24 Fgure 7. Comparson n Bode format (modulus between measured and smulated data Fgure 8. Comparson n Bode format (phase between measured and smulated data

An impedance model for analysis of EIS of polymer electrolyte fuel cells under platinum oxidation and hydrogen peroxide formation in the cathode

An impedance model for analysis of EIS of polymer electrolyte fuel cells under platinum oxidation and hydrogen peroxide formation in the cathode Loughborough Unversty Insttutonal epostory An mpedane model for analyss of EIS of polymer eletrolyte fuel ells under platnum oxdaton and hydrogen peroxde formaton n the athode Ths tem was submtted to Loughborough

More information

Voltammetry. Bulk electrolysis: relatively large electrodes (on the order of cm 2 ) Voltammetry:

Voltammetry. Bulk electrolysis: relatively large electrodes (on the order of cm 2 ) Voltammetry: Voltammetry varety of eletroanalytal methods rely on the applaton of a potental funton to an eletrode wth the measurement of the resultng urrent n the ell. In ontrast wth bul eletrolyss methods, the objetve

More information

Impedance Analysis of Molten Carbonate Fuel Cell

Impedance Analysis of Molten Carbonate Fuel Cell Impedane Analyss of Molten Carbonate Fuel Cell Fnal Projet Report In partal requrement of ECHE 789B ourse Naln Subramanan Course Instrutor: Dr. Branko N. Popov Date submtted: May 5,00 Abstrat A three phase

More information

An Impedance Model for EIS Analysis of Nickel Metal Hydride Batteries

An Impedance Model for EIS Analysis of Nickel Metal Hydride Batteries Ths s a post-refereeng fnal draft. When tng, please refer to the publshed verson: S. ruz-manzo, P. Greenwood and. hen, J. Eletrohem. So. 207 volume 64, ssue 7, A446- A453. An Impedane Model for EIS Analyss

More information

A Theorem of Mass Being Derived From Electrical Standing Waves (As Applied to Jean Louis Naudin's Test)

A Theorem of Mass Being Derived From Electrical Standing Waves (As Applied to Jean Louis Naudin's Test) A Theorem of Mass Beng Derved From Eletral Standng Waves (As Appled to Jean Lous Naudn's Test) - by - Jerry E Bayles Aprl 4, 000 Ths paper formalzes a onept presented n my book, "Eletrogravtaton As A Unfed

More information

A Theorem of Mass Being Derived From Electrical Standing Waves (As Applied to Jean Louis Naudin's Test)

A Theorem of Mass Being Derived From Electrical Standing Waves (As Applied to Jean Louis Naudin's Test) A Theorem of Mass Beng Derved From Eletral Standng Waves (As Appled to Jean Lous Naudn's Test) - by - Jerry E Bayles Aprl 5, 000 Ths Analyss Proposes The Neessary Changes Requred For A Workng Test Ths

More information

Exercise 10: Theory of mass transfer coefficient at boundary

Exercise 10: Theory of mass transfer coefficient at boundary Partle Tehnology Laboratory Prof. Sotrs E. Pratsns Sonneggstrasse, ML F, ETH Zentrum Tel.: +--6 5 http://www.ptl.ethz.h 5-97- U Stoffaustaush HS 7 Exerse : Theory of mass transfer oeffent at boundary Chapter,

More information

Phase Transition in Collective Motion

Phase Transition in Collective Motion Phase Transton n Colletve Moton Hefe Hu May 4, 2008 Abstrat There has been a hgh nterest n studyng the olletve behavor of organsms n reent years. When the densty of lvng systems s nreased, a phase transton

More information

On Adaptive Control of Simulated Moving Bed Plants. Plants Using Comsol s Simulink Interface. Speaker: Marco Fütterer

On Adaptive Control of Simulated Moving Bed Plants. Plants Using Comsol s Simulink Interface. Speaker: Marco Fütterer daptve Smulated Movng ed Plants Usng Comsol s Smulnk Interfae Speaker: Maro Fütterer Insttut für utomatserungstehnk Otto-von-Guerke Unverstät Unverstätsplatz, D-39106 Magdeburg Germany e-mal: maro.fuetterer@ovgu.de

More information

Gravity Drainage Prior to Cake Filtration

Gravity Drainage Prior to Cake Filtration 1 Gravty Dranage Pror to ake Fltraton Sott A. Wells and Gregory K. Savage Department of vl Engneerng Portland State Unversty Portland, Oregon 97207-0751 Voe (503) 725-4276 Fax (503) 725-4298 ttp://www.e.pdx.edu/~wellss

More information

Analytical calculation of adiabatic processes in real gases

Analytical calculation of adiabatic processes in real gases Journal of Physs: Conferene Seres PAPER OPEN ACCESS Analytal alulaton of adabat roesses n real gases o te ths artle: I B Amarskaja et al 016 J. Phys.: Conf. Ser. 754 11003 Related ontent - Shortuts to

More information

Introduction to Molecular Spectroscopy

Introduction to Molecular Spectroscopy Chem 5.6, Fall 004 Leture #36 Page Introduton to Moleular Spetrosopy QM s essental for understandng moleular spetra and spetrosopy. In ths leture we delneate some features of NMR as an ntrodutory example

More information

technische universiteit eindhoven Analysis of one product /one location inventory control models prof.dr. A.G. de Kok 1

technische universiteit eindhoven Analysis of one product /one location inventory control models prof.dr. A.G. de Kok 1 TU/e tehnshe unverstet endhoven Analyss of one produt /one loaton nventory ontrol models prof.dr. A.G. de Kok Aknowledgements: I would lke to thank Leonard Fortun for translatng ths ourse materal nto Englsh

More information

A generic electrical circuit for performance analysis of the fuel cell cathode catalyst layer through electrochemical impedance spectroscopy

A generic electrical circuit for performance analysis of the fuel cell cathode catalyst layer through electrochemical impedance spectroscopy Loughborough University Institutional Repository A generi eletrial iruit for performane analysis of the fuel ell athode atalyst layer through eletrohemial impedane spetrosopy This item was submitted to

More information

ECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 3

ECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 3 C 634 Intermedate M Waves Fall 216 Prof. Davd R. akson Dept. of C Notes 3 1 Types of Current ρ v Note: The free-harge densty ρ v refers to those harge arrers (ether postve or negatve) that are free to

More information

Frequency dependence of the permittivity

Frequency dependence of the permittivity Frequency dependence of the permttvty February 7, 016 In materals, the delectrc constant and permeablty are actually frequency dependent. Ths does not affect our results for sngle frequency modes, but

More information

3D Numerical Analysis for Impedance Calculation and High Performance Consideration of Linear Induction Motor for Rail-guided Transportation

3D Numerical Analysis for Impedance Calculation and High Performance Consideration of Linear Induction Motor for Rail-guided Transportation ADVANCED ELECTROMAGNETICS SYMPOSIUM, AES 13, 19 MARCH 13, SHARJAH UNITED ARAB EMIRATES 3D Numeral Analss for Impedane Calulaton and Hgh Performane Consderaton of Lnear Induton Motor for Ral-guded Transportaton

More information

Controller Design for Networked Control Systems in Multiple-packet Transmission with Random Delays

Controller Design for Networked Control Systems in Multiple-packet Transmission with Random Delays Appled Mehans and Materals Onlne: 03-0- ISSN: 66-748, Vols. 78-80, pp 60-604 do:0.408/www.sentf.net/amm.78-80.60 03 rans eh Publatons, Swtzerland H Controller Desgn for Networed Control Systems n Multple-paet

More information

The calculation of ternary vapor-liquid system equilibrium by using P-R equation of state

The calculation of ternary vapor-liquid system equilibrium by using P-R equation of state The alulaton of ternary vapor-lqud syste equlbru by usng P-R equaton of state Y Lu, Janzhong Yn *, Rune Lu, Wenhua Sh and We We Shool of Cheal Engneerng, Dalan Unversty of Tehnology, Dalan 11601, P.R.Chna

More information

FAULT DETECTION AND IDENTIFICATION BASED ON FULLY-DECOUPLED PARITY EQUATION

FAULT DETECTION AND IDENTIFICATION BASED ON FULLY-DECOUPLED PARITY EQUATION Control 4, Unversty of Bath, UK, September 4 FAUL DEECION AND IDENIFICAION BASED ON FULLY-DECOUPLED PARIY EQUAION C. W. Chan, Hua Song, and Hong-Yue Zhang he Unversty of Hong Kong, Hong Kong, Chna, Emal:

More information

Curvature Dependency of Surface Tension in Multicomponent Systems

Curvature Dependency of Surface Tension in Multicomponent Systems HERMODYNAMICS Curvature Dependeny of Surfae enson n Multomponent Systems Erk Santso Dept. of Chemal Engneerng, North Carolna State Unversty, Ralegh, NC 27695 Abbas Froozabad Reservor Engneerng Researh

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Machine Learning: and 15781, 2003 Assignment 4

Machine Learning: and 15781, 2003 Assignment 4 ahne Learnng: 070 and 578, 003 Assgnment 4. VC Dmenson 30 onts Consder the spae of nstane X orrespondng to all ponts n the D x, plane. Gve the VC dmenson of the followng hpothess spaes. No explanaton requred.

More information

Charged Particle in a Magnetic Field

Charged Particle in a Magnetic Field Charged Partle n a Magnet Feld Mhael Fowler 1/16/08 Introduton Classall, the fore on a harged partle n eletr and magnet felds s gven b the Lorentz fore law: v B F = q E+ Ths velot-dependent fore s qute

More information

GEL 446: Applied Environmental Geology

GEL 446: Applied Environmental Geology GE 446: ppled Envronmental Geology Watershed Delneaton and Geomorphology Watershed Geomorphology Watersheds are fundamental geospatal unts that provde a physal and oneptual framewor wdely used by sentsts,

More information

Theory of Ferromagnetism in Double Perosvkites. Luis Brey CSIC-Madrid F. Guinea CSIC-Madrid S.Das Sarma Univ.Maryland

Theory of Ferromagnetism in Double Perosvkites. Luis Brey CSIC-Madrid F. Guinea CSIC-Madrid S.Das Sarma Univ.Maryland Theory of Ferromagnetsm n Double Perosvktes. Lus Brey CSIC-Madrd F. Gunea CSIC-Madrd S.Das Sarma Unv.Maryland 1 OUTLINE Introduton to Fe based double perovsktes. Chemstry Band struture Ferromagnetsm ndued

More information

Energy-Based State-Space Representation of Modular Multilevel Converters with a Constant Equilibrium Point in Steady-State Operation

Energy-Based State-Space Representation of Modular Multilevel Converters with a Constant Equilibrium Point in Steady-State Operation Energy-ased State-Spae Representaton of Modular Multlevel Converters wth a Constant Equlbrum Pont n Steady-State Operaton Glbert ergna-daz, Jon re Suul, Member IEEE, and Salvatore D'ro bstrat he nternal

More information

Bidimensional Analysis of a Thermoelectric Module using Finite Element Techniques

Bidimensional Analysis of a Thermoelectric Module using Finite Element Techniques Bdmensonal Analyss of a Termoeletr Module usng Fnte Element Tenques *Antono Arenas, Jorge Vázquez, Rafael Palaos Unversdad Pontfa Comllas Esuela Téna Superor de ngenería *Departamento de Fludos y Calor,

More information

Opportunities in Analytical Approaches to Spray Drying of Solid Dosage Forms

Opportunities in Analytical Approaches to Spray Drying of Solid Dosage Forms Opportuntes n Analytal Approahes to Spray Dryng of Sold Dosage Forms Dr. Renhard Vehrng Assoate Professor and George Ford Char n Materals Engneerng Unversty of Alberta, Department of Mehanal Engneerng

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

ALGEBRAIC SCHUR COMPLEMENT APPROACH FOR A NON LINEAR 2D ADVECTION DIFFUSION EQUATION

ALGEBRAIC SCHUR COMPLEMENT APPROACH FOR A NON LINEAR 2D ADVECTION DIFFUSION EQUATION st Annual Internatonal Interdsplnary Conferene AIIC 03 4-6 Aprl Azores Portugal - Proeedngs- ALGEBRAIC SCHUR COMPLEMENT APPROACH FOR A NON LINEAR D ADVECTION DIFFUSION EQUATION Hassan Belhad Professor

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

425. Calculation of stresses in the coating of a vibrating beam

425. Calculation of stresses in the coating of a vibrating beam 45. CALCULAION OF SRESSES IN HE COAING OF A VIBRAING BEAM. 45. Calulaton of stresses n the oatng of a vbratng beam M. Ragulsks,a, V. Kravčenken,b, K. Plkauskas,, R. Maskelunas,a, L. Zubavčus,b, P. Paškevčus,d

More information

JSM Survey Research Methods Section. Is it MAR or NMAR? Michail Sverchkov

JSM Survey Research Methods Section. Is it MAR or NMAR? Michail Sverchkov JSM 2013 - Survey Researh Methods Seton Is t MAR or NMAR? Mhal Sverhkov Bureau of Labor Statsts 2 Massahusetts Avenue, NE, Sute 1950, Washngton, DC. 20212, Sverhkov.Mhael@bls.gov Abstrat Most methods that

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

STK4900/ Lecture 4 Program. Counterfactuals and causal effects. Example (cf. practical exercise 10)

STK4900/ Lecture 4 Program. Counterfactuals and causal effects. Example (cf. practical exercise 10) STK4900/9900 - Leture 4 Program 1. Counterfatuals and ausal effets 2. Confoundng 3. Interaton 4. More on ANOVA Setons 4.1, 4.4, 4.6 Supplementary materal on ANOVA Example (f. pratal exerse 10) How does

More information

ELG4179: Wireless Communication Fundamentals S.Loyka. Frequency-Selective and Time-Varying Channels

ELG4179: Wireless Communication Fundamentals S.Loyka. Frequency-Selective and Time-Varying Channels Frequeny-Seletve and Tme-Varyng Channels Ampltude flutuatons are not the only effet. Wreless hannel an be frequeny seletve (.e. not flat) and tmevaryng. Frequeny flat/frequeny-seletve hannels Frequeny

More information

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit.

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit. Physcs 4B Solutons to Chapter 7 HW Chapter 7: Questons:, 8, 0 Problems:,,, 45, 48,,, 7, 9 Queston 7- (a) no (b) yes (c) all te Queston 7-8 0 μc Queston 7-0, c;, a;, d; 4, b Problem 7- (a) Let be the current

More information

DEVELOPMENT OF TOOLS TO ESTIMATE ACTUAL CORROSION GROWTH RATES OF GAS PIPELINES

DEVELOPMENT OF TOOLS TO ESTIMATE ACTUAL CORROSION GROWTH RATES OF GAS PIPELINES DEVELOPMENT OF TOOLS TO ESTIMATE ACTUAL CORROSION GROWTH RATES OF GAS PIPELINES Fnal Report SwRI Projet.148 Contrat Number: DTPH 56-8-T- Prepared for U.S. Department of Transportaton Ppelne Hazardous Materals

More information

Brander and Lewis (1986) Link the relationship between financial and product sides of a firm.

Brander and Lewis (1986) Link the relationship between financial and product sides of a firm. Brander and Lews (1986) Lnk the relatonshp between fnanal and produt sdes of a frm. The way a frm fnanes ts nvestment: (1) Debt: Borrowng from banks, n bond market, et. Debt holders have prorty over a

More information

PHYSICS 212 MIDTERM II 19 February 2003

PHYSICS 212 MIDTERM II 19 February 2003 PHYSICS 1 MIDERM II 19 Feruary 003 Exam s losed ook, losed notes. Use only your formula sheet. Wrte all work and answers n exam ooklets. he aks of pages wll not e graded unless you so request on the front

More information

Some Results on the Counterfeit Coins Problem. Li An-Ping. Beijing , P.R.China Abstract

Some Results on the Counterfeit Coins Problem. Li An-Ping. Beijing , P.R.China Abstract Some Results on the Counterfet Cons Problem L An-Png Bejng 100085, P.R.Chna apl0001@sna.om Abstrat We wll present some results on the ounterfet ons problem n the ase of mult-sets. Keywords: ombnatoral

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

DOAEstimationforCoherentSourcesinBeamspace UsingSpatialSmoothing

DOAEstimationforCoherentSourcesinBeamspace UsingSpatialSmoothing DOAEstmatonorCoherentSouresneamspae UsngSpatalSmoothng YnYang,ChunruWan,ChaoSun,QngWang ShooloEletralandEletronEngneerng NanangehnologalUnverst,Sngapore,639798 InsttuteoAoustEngneerng NorthwesternPoltehnalUnverst,X

More information

Summary ELECTROMAGNETIC FIELDS AT THE WORKPLACES. System layout: exposure to magnetic field only. Quasi-static dosimetric analysis: system layout

Summary ELECTROMAGNETIC FIELDS AT THE WORKPLACES. System layout: exposure to magnetic field only. Quasi-static dosimetric analysis: system layout Internatonal Workshop on LCTROMGNTIC FILDS T TH WORKPLCS 5-7 September 5 Warszawa POLND 3d approah to numeral dosmetr n quas-stat ondtons: problems and eample of solutons Dr. Nola Zoppett - IFC-CNR, Florene,

More information

Improving the Performance of Fading Channel Simulators Using New Parameterization Method

Improving the Performance of Fading Channel Simulators Using New Parameterization Method Internatonal Journal of Eletrons and Eletral Engneerng Vol. 4, No. 5, Otober 06 Improvng the Performane of Fadng Channel Smulators Usng New Parameterzaton Method Omar Alzoub and Moheldn Wanakh Department

More information

The corresponding link function is the complementary log-log link The logistic model is comparable with the probit model if

The corresponding link function is the complementary log-log link The logistic model is comparable with the probit model if SK300 and SK400 Lnk funtons for bnomal GLMs Autumn 08 We motvate the dsusson by the beetle eample GLMs for bnomal and multnomal data Covers the followng materal from hapters 5 and 6: Seton 5.6., 5.6.3,

More information

Electrostatic Potential from Transmembrane Currents

Electrostatic Potential from Transmembrane Currents Electrostatc Potental from Transmembrane Currents Let s assume that the current densty j(r, t) s ohmc;.e., lnearly proportonal to the electrc feld E(r, t): j = σ c (r)e (1) wth conductvty σ c = σ c (r).

More information

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

Curve Fitting with the Least Square Method

Curve Fitting with the Least Square Method WIKI Document Number 5 Interpolaton wth Least Squares Curve Fttng wth the Least Square Method Mattheu Bultelle Department of Bo-Engneerng Imperal College, London Context We wsh to model the postve feedback

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Implementation of α-qss Stiff Integration Methods for Solving the Detailed Combustion Chemistry

Implementation of α-qss Stiff Integration Methods for Solving the Detailed Combustion Chemistry Proeedngs of the World Congress on Engneerng 2007 Vol II Implementaton of α-qss Stff Integraton Methods for Solvng the Detaled Combuston Chemstry Shafq R. Quresh and Robert Prosser Abstrat Implt methods

More information

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

More information

Role of long range electron - phonon and coulomb interactions in high - T c cuprate superconductors

Role of long range electron - phonon and coulomb interactions in high - T c cuprate superconductors Internatonal Journal of Latest Engneerng Researh and Applatons (IJLERA) ISSN: -77 Volume, Issue, Aprl 7, PP 9-99 Role of long range eletron - phonon and oulomb nteratons n hgh - T uprate superondutors

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics )

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics ) Ismor Fscher, 8//008 Stat 54 / -8.3 Summary Statstcs Measures of Center and Spread Dstrbuton of dscrete contnuous POPULATION Random Varable, numercal True center =??? True spread =???? parameters ( populaton

More information

Assignment 4. Adsorption Isotherms

Assignment 4. Adsorption Isotherms Insttute of Process Engneerng Assgnment 4. Adsorpton Isotherms Part A: Compettve adsorpton of methane and ethane In large scale adsorpton processes, more than one compound from a mxture of gases get adsorbed,

More information

Instance-Based Learning and Clustering

Instance-Based Learning and Clustering Instane-Based Learnng and Clusterng R&N 04, a bt of 03 Dfferent knds of Indutve Learnng Supervsed learnng Bas dea: Learn an approxmaton for a funton y=f(x based on labelled examples { (x,y, (x,y,, (x n,y

More information

Performance Modeling of Hierarchical Memories

Performance Modeling of Hierarchical Memories Performane Modelng of Herarhal Memores Marwan Sleman, Lester Lpsky, Kshor Konwar Department of omputer Sene and Engneerng Unversty of onnetut Storrs, T 0669-55 Emal: {marwan, lester, kshor}@engr.uonn.edu

More information

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017) Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed

More information

CHAPTER 13. Exercises. E13.1 The emitter current is given by the Shockley equation:

CHAPTER 13. Exercises. E13.1 The emitter current is given by the Shockley equation: HPT 3 xercses 3. The emtter current s gen by the Shockley equaton: S exp VT For operaton wth, we hae exp >> S >>, and we can wrte VT S exp VT Solng for, we hae 3. 0 6ln 78.4 mv 0 0.784 5 4.86 V VT ln 4

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

e a = 12.4 i a = 13.5i h a = xi + yj 3 a Let r a = 25cos(20) i + 25sin(20) j b = 15cos(55) i + 15sin(55) j

e a = 12.4 i a = 13.5i h a = xi + yj 3 a Let r a = 25cos(20) i + 25sin(20) j b = 15cos(55) i + 15sin(55) j Vetors MC Qld-3 49 Chapter 3 Vetors Exerse 3A Revew of vetors a d e f e a x + y omponent: x a os(θ 6 os(80 + 39 6 os(9.4 omponent: y a sn(θ 6 sn(9 0. a.4 0. f a x + y omponent: x a os(θ 5 os( 5 3.6 omponent:

More information

4.5. QUANTIZED RADIATION FIELD

4.5. QUANTIZED RADIATION FIELD 4-1 4.5. QUANTIZED RADIATION FIELD Baground Our treatent of the vetor potental has drawn on the onohroat plane-wave soluton to the wave-euaton for A. The uantu treatent of lght as a partle desrbes the

More information

Dr. M. Perumal Professor & Head Department of Hydrology Indian Institute of Technology Roorkee INDIA Co-authors: Dr. B. Sahoo & Dr. C.M.

Dr. M. Perumal Professor & Head Department of Hydrology Indian Institute of Technology Roorkee INDIA Co-authors: Dr. B. Sahoo & Dr. C.M. Dr.. Perumal Professor & Head Department of Hdrolog Indan Insttute of Tehnolog Roorkee INDIA o-authors: Dr. B. Sahoo & Dr... Rao Dr. Dr... Perumal, Professor & & Head, Dept. Dept. of of Hdrolog, I.I.T.

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

AGC Introduction

AGC Introduction . Introducton AGC 3 The prmary controller response to a load/generaton mbalance results n generaton adjustment so as to mantan load/generaton balance. However, due to droop, t also results n a non-zero

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists *

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists * How Strong Are Weak Patents? Joseph Farrell and Carl Shapro Supplementary Materal Lcensng Probablstc Patents to Cournot Olgopolsts * September 007 We study here the specal case n whch downstream competton

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Prediction of Solid Paraffin Precipitation Using Solid Phase Equation of State

Prediction of Solid Paraffin Precipitation Using Solid Phase Equation of State Predton of old Paraffn Preptaton Usng old Phase Equaton of tate Proeedngs of European Congress of Chemal Engneerng (ECCE-6) Copenhagen, 16- eptember 7 Predton of old Paraffn Preptaton Usng old Phase Equaton

More information

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products modelng of equlbrum and dynamc mult-component adsorpton n a two-layered fxed bed for purfcaton of hydrogen from methane reformng products Mohammad A. Ebrahm, Mahmood R. G. Arsalan, Shohreh Fatem * Laboratory

More information

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016 ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have

More information

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results. Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson

More information

Horizontal mergers for buyer power. Abstract

Horizontal mergers for buyer power. Abstract Horzontal mergers for buyer power Ramon Faul-Oller Unverstat d'alaant Llus Bru Unverstat de les Illes Balears Abstrat Salant et al. (1983) showed n a Cournot settng that horzontal mergers are unproftable

More information

A Theorem of Mass Being Derived From Electrical Standing Waves

A Theorem of Mass Being Derived From Electrical Standing Waves A Theorem of Mass Beng Derved From Eletral Standng Waves - by - Jerry E Bayles Ths paper formalzes a onept presented n my book, "Eletrogravtaton As A Unfed Feld Theory", (as well as n numerous related

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

Complement of an Extended Fuzzy Set

Complement of an Extended Fuzzy Set Internatonal Journal of Computer pplatons (0975 8887) Complement of an Extended Fuzzy Set Trdv Jyot Neog Researh Sholar epartment of Mathemats CMJ Unversty, Shllong, Meghalaya usmanta Kumar Sut ssstant

More information

Experimental. SC Solutions, Inc. The COMSOL. cal-mechanicall. wafers. We. for semiconductor. and species transport. model was fabricated,

Experimental. SC Solutions, Inc. The COMSOL. cal-mechanicall. wafers. We. for semiconductor. and species transport. model was fabricated, Expermental Valdaton of Model of Eletro-Chemal-Mehanal Planarzaton (ECMP) of Copper J. L. Ebert*, S. Ghosal, D. de Roover, and A. Emam-Naen SC Solutons, In. *Correspondng author: SC Solutons, In., 161

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

More information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

More information

Statistical Evaluation of WATFLOOD

Statistical Evaluation of WATFLOOD tatstcal Evaluaton of WATFLD By: Angela MacLean, Dept. of Cvl & Envronmental Engneerng, Unversty of Waterloo, n. ctober, 005 The statstcs program assocated wth WATFLD uses spl.csv fle that s produced wth

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

Prediction of the reliability of genomic breeding values for crossbred performance

Prediction of the reliability of genomic breeding values for crossbred performance Vandenplas et al. Genet Sel Evol 217 49:43 DOI 1.1186/s12711-17-318-1 Genets Seleton Evoluton RESERCH RTICLE Open ess Predton of the relablty of genom breedng values for rossbred performane Jéréme Vandenplas

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015 Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want

More information

1) Silicon oxide has a typical surface potential in an aqueous medium of ϕ,0

1) Silicon oxide has a typical surface potential in an aqueous medium of ϕ,0 1) Slcon oxde has a typcal surface potental n an aqueous medum of ϕ, = 7 mv n 5 mm l at ph 9. Whch concentraton of catons do you roughly expect close to the surface? What s the average dstance between

More information

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS 6th ICPT, Sapporo, Japan, July 008 APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS James MAINA Prncpal Researcher, Transport and Infrastructure Engneerng, CSIR Bult Envronment

More information

JAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger

JAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger JAB Chan Long-tal clams development ASTIN - September 2005 B.Verder A. Klnger Outlne Chan Ladder : comments A frst soluton: Munch Chan Ladder JAB Chan Chan Ladder: Comments Black lne: average pad to ncurred

More information

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming EEL 6266 Power System Operaton and Control Chapter 3 Economc Dspatch Usng Dynamc Programmng Pecewse Lnear Cost Functons Common practce many utltes prefer to represent ther generator cost functons as sngle-

More information

Physics 41 Chapter 22 HW Serway 7 th Edition

Physics 41 Chapter 22 HW Serway 7 th Edition yss 41 apter H Serway 7 t Edton oneptual uestons: 1,, 8, 1 roblems: 9, 1, 0,, 7, 9, 48, 54, 55 oneptual uestons: 1,, 8, 1 1 Frst, te effeny of te automoble engne annot exeed te arnot effeny: t s lmted

More information

12 th National Congress on Theoretical and Applied Mechanics September 2013, Saints Constantine and Helena, Varna, Bulgaria

12 th National Congress on Theoretical and Applied Mechanics September 2013, Saints Constantine and Helena, Varna, Bulgaria th Natonal Congress on Theoretal and Appled Mehans 3-6 September 03 Sants Constantne and Helena Varna Bulgara APPLICATION OF GEARING PRIMITIVES TO SKEW AXES GEAR SET SYNTHESIS PART : MATHEMATICAL MODEL

More information

3rd International Conference on Mechanical Engineering and Intelligent Systems (ICMEIS 2015)

3rd International Conference on Mechanical Engineering and Intelligent Systems (ICMEIS 2015) 3rd Internatonal Conferene on Mehanal Engneerng and Intellgent Systems (ICMEIS 15) Researh on Transmsson Lne Loss and Carryng Current Based on Temperature Power Flow Model Le LUO a*, Xngong CHENG b, Xu

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

ECE 422 Power System Operations & Planning 2 Synchronous Machine Modeling

ECE 422 Power System Operations & Planning 2 Synchronous Machine Modeling ECE 422 Power System Operatons & Plannng 2 Synhronous Mahne Moelng Sprng 219 Instrutor: Ka Sun 1 Outlne 2.1 Moelng of synhronous generators for Stablty Stues Synhronous Mahne Moelng Smplfe Moels for Stablty

More information

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors Multple Lnear and Polynomal Regresson wth Statstcal Analyss Gven a set of data of measured (or observed) values of a dependent varable: y versus n ndependent varables x 1, x, x n, multple lnear regresson

More information

Effects of vertical wall and tetrapod weights on wave overtopping in rubble mound breakwaters under irregular wave conditions

Effects of vertical wall and tetrapod weights on wave overtopping in rubble mound breakwaters under irregular wave conditions SNAK, 2014 Int. J. Nav. Arht. Oean Eng. (2014) 6:947~964 http://dx.do.org/10.2478/ijnaoe-2013-0224 pissn: 2092-6782, eissn: 2092-6790 Effets of vertal wall and tetrapod weghts on wave overtoppng n rubble

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information