The Representation of Multi-component Adsorption in Reservoir Simulation of CO 2 Sequestration in Coal and Enhanced Coalbed Methane Recovery

Size: px
Start display at page:

Download "The Representation of Multi-component Adsorption in Reservoir Simulation of CO 2 Sequestration in Coal and Enhanced Coalbed Methane Recovery"

Transcription

1 The Representton of Mult-component Adsorpton n Reservor Smulton of CO 2 Sequestrton n Col nd Enhnced Colbed Methne Recovery 59 Zhejun Pn nd LD Connell CSIRO Petroleum, Prvte Bg 1, Clyton South, Vctor, Austrl, 3169 ABSTRACT CO 2 sequestrton n col s potentl mngement opton for greenhouse gs emssons. An ttrctve spect to ths process s tht the CO 2 s dsorbed to the col nd therefore, for sttc reservor pressure condtons, s mmoble, reducng the rsk of CO 2 mgrton to surfce. Another spect to ths s tht the njected CO 2 cn dsplce dsorbed methne ledng to regon wthn the col sem of mxed gs. Therefore, n order to understnd gs mgrton wthn the reservor mult-component dsorpton models re requred. Lngmur bsed models, such s the extended Lngmur nd Idel Adsorbed Soluton, re prt of the reservor smultors currently n use for modelng CO 2 sequestrton n col. However these models re nccurte descrptons of the dsorpton process prtculrly t the hgh pressures ssocted wth sequestrton n deep cols. Therefore new pproches re requred tht cn be mplemented n computtonl effcent mnner for reservor smulton. Ths pper presents smulton work usng severl pproches to representng mult-component dsorpton, mplemented nto the col sem gs reservor smultor SIMED II. The dsorpton models re extended Lngmur, Idel Adsorbed Soluton () nd 2D Equton of Stte (). The s shown to be more ccurte descrpton thn the other models for expermentl observtons of dsorpton on Frutlnd col. The effects of ths greter ccurcy on modelng the sequestrton process s nvestgted usng hypothetcl exmple. Whle smulton results for the extended Lngmur nd re very close, results clculted usng re sgnfcntly dfferent. The dfferences re gretest wth the CO 2 njecton rte. A concluson from ths work s tht nccurte sotherm models cn led to sgnfcnt errors n smulton results nd tht the mplementton of the model provdes greter ccurcy n the representton of the sequestrton process.

2 2 The Representton of Mult-component Adsorpton n Reservor Smulton of CO 2 Sequestrton n Col nd Enhnced Colbed Methne Recovery INTRODUCTION CO 2 sequestrton n unmnble col s potentl mngement opton for greenhouse gs emssons. The fesblty of ths process wll be determned by vrous fctors ncludng the col cpcty, the rte of njecton nd the long term behvour of the sequestered CO 2. To nvestgte these requres representtve models for the reservor behvour of the CO 2. One mportnt process n determnng the reservor behvour s dsorpton nd snce mgrton wll often nvolve the dsplcement of dsorbed methne wth CO 2 nd therefore regon of mxed gs wthn the sem, t s mult-component dsorpton tht needs to be represented. Currently, most of the reservor smultors used for sequestrton n col mplement the Lngmur nd Extended Lngmur models () for pure nd mxed gs dsorpton, respectvely. Snce s not thermodynmclly correct [1] nd subject to lrge errors for mxed gs dsorpton, Idel Adsorbed Soluton () theory ws mplemented n some smultors. However, s for mxture dsorpton only nd requres pure gs dsorpton model, nd normlly the Lngmur model s used. Thus, the nccurcy of the pure gs dsorpton propgtes to the mxed gs dsorpton clculton. Recent work hs proposed new dsorpton models for mult-component dsorpton for the full rnge of pressures, such s the Two-Dmensonl Equton of Stte () pproch. Ths pper presents smulton work usng severl pproches to representng mult-component dsorpton, mplemented nto the col sem gs reservor smultor SIMED II. The dsorpton models re extended Lngmur, Idel Adsorbed Soluton nd. The smulton studes nvestgte the sequestrton process nd compre the results usng these dfferent dsorpton pproches. ADSORPTION THEORY Lngmur nd Extended Lngmur Model The Lngmur model ssumes dsorpton occurs on flt surfce s gven by the knetc theory. At equlbrum, contnul process of bombrdment of molecules onto the surfce nd correspondng evporton of molecules from the surfce mntn zero rte of ccumulton t the surfce. The ssumptons of the Lngmur model re [1]: (1) The surfce s homogeneous; tht s, the dsorpton energy s constnt over ll stes. (2) Adsorpton on the surfce s loclzed, whch mens tht the toms or molecules re dsorbed t defnte, loclzed stes. (3) Ech ste cn ccommodte only one molecule or tom. Bsed on the ssumptons bove, the Lngmur model cn be derved s: ω BP θ = = (1) L 1+ BP Where θ s the frctonl coverge, ω s the mount dsorbed, P s the pressure, B s the Lngmur constnt, nd L s mxmum dsorpton cpcty. When extended to mxture dsorpton ssumng no ntercton mong the dsorbed molecules, the Lngmur model becomes: LBPy ω = (2) 1 B Py + j j j where subscrpt represents component, nd y s the mole frcton n the gs phse. Lngmur model s successful n modelng dsorpton t low or md pressure (typclly up to 8 MP for gs dsorpton on cols). However, the dfference between the model nd the expermentl dt s lrge t hgh pressure rnges. Moreover, the extended Lngmur model () s not thermodynmclly correct unless the mxmum cpcty L s the sme for ech component [1]. Consequently, lrge errors my occur for hgh pressure mult-component dsorpton predctons, prtculrly when the dfference of L s lrge.

3 Pn, Connell 3 Idel Adsorbed Soluton The Idel Adsorbed Soluton () model s for mxture dsorpton only, developed by Myers nd Prusntz (1965) [2]. The s n dsorpton nlog to the Roult s Lw for vpor-lqud equlbrum. The ssumptons of the model re [3]: (1) The dsorbed solutons nd the gs phse re del. (2) All ctvty coeffcents n the dsorbed phse re unty. The equlbrum condton for the dsorbed phse nd the gs phse s: ˆ g Py φ = P x φ P 1... ( ) N = where P s the gs pressure of component dsorbed t the sme spredng pressure s the mxture, g ˆφ s the prtl fugcty coeffcent of component n the gs phse, x s the dsorbed phse composton nd φ s the pure gs fugcty coeffcent. The spredng pressure π cn be determned by the Gbbs dsorpton sotherm provdng tht the dsorbed phse s del: π = RT P ω d ln P (4) A where R s the gs constnt, T s the temperture nd A s the surfce re of the dsorbent. The theory ssumes tht the spredng pressure of ech component should be equvlent to the multcomponent dsorpton. When ssumng ech component ccesses the sme surfce re, the followng relton pples: P1 ω 2 ω n ω 1 P P 2 n dp = dp =... = dp (5) P P P Menwhle, the mole frcton constrnt s: n = 1 n x = y = 1 = 1 The totl dsorpton, ω T, s gven by: n 1 x = ω P T = 1 ω ( ) The component dsorpton, ω, s gven by: ω = ω T x To perform mxed-gs dsorpton clcultons, pure-component model s needed. Any purecomponent model, for nstnce the Lngmur model, cn be utlzed n the clculton. However, for mny models, nlytcl ntegrton s unvlble. Thus numercl ntegrton my be used n Equton (5). For hgh pressure dsorpton, fugcty should be used n the gs phse. However, ssumng n del dsorbed phse my be nvld t hgh pressure ledng to lrge errors for these condtons. Two-Dmensonl Equton of Stte The generl 2D nlog cn be wrtten s follows wth n ddtonl prmeter m for dded model flexblty [4]: 2 αω m Aπ + [ 1 ( βω) ] = ωrt 2 (9) 1+ Uβω + W( βω) where A s the specfc surfce re, π s the spredng pressure, ω s the specfc mount dsorbed, nd α nd β re model prmeters. The model coeffcents, U, W, nd m must be specfed to obtn specfc form of the for pplcton. For exmple, n nlog of the vn der Wls (VDW) EOS s obtned by settng m = 1 nd U = W = ; smlrly for the Sove-Redlch-Kwong (SRK) (m = U = 1 nd W (3) (6) (7) (8)

4 4 The Representton of Mult-component Adsorpton n Reservor Smulton of CO 2 Sequestrton n Col nd Enhnced Colbed Methne Recovery = ); the Peng-Robnson (PR) (m = 1, U = 2, nd W = -1); the Eyrng (m = 1/2 nd U = W = ) EOS; nd the Zhou-Gsem-Robnson (ZGR) EOS (m = 1/3 nd U = W = ). [4] One-flud mxng rules were used n the to descrbe mult-component dsorpton: α = x x jαj nd β = x x jβj (1) j j The combnton rules below re chosen becuse they work best for dsorbed phse [4]: α = α + α C / (11) j j ( j )( 1 j ) 2 β β ( 1+ D ) β = (12) j j where C j nd D j re bnry ntercton prmeters. To predct the mxture dsorpton, C j nd D j re set to be zero. To relte the dsorbed phse wth the gs phse, t equlbrum, the chemcl potentl of speces n g the dsorbed phse, µ, s equl to tht n the gs phse, : µ Thus, π g = µ d ln fˆ = π * g nd d µ = dµ (13) P P * d ln fˆ g Integrtng Equton 14, t yelds: ln fˆ * g ( π) ln fˆ ( ) ( ) ( ) g * π = ln fˆ P ln fˆ P ˆ * * g * * ( π ) = π fˆ ( P ) = P At very low pressure, f nd * P fˆ * g ( π) = π fˆ ( P) At very low pressure, 2D del gs lw pples: * * π A = ω RT Henry s constnt s defned s: * ω k = (18) * P So, substtuton nto Equton 16, we obtn: g A fˆ = Ax πφˆ = k RTfˆ For pure-gs dsorpton, Equton 19 becomes: g ωz φ = k f The fugcty for the s: ω 1 ( A ) 1 ln ˆ π φ = dω ln Z RTω ω ω T,M s, n j where ω s the mount dsorbed, s the 2D compressblty fctor, φ s the fugcty coeffcent usng Z g the, f s the fugcty for the gs phse. The 3D PR EOS s normlly used to clculte the gs fugcty. Detled dervton of Equton 21 ws provded by Zhou [3]. µ (14) (15) (16) (17) (19) (2) (21) ADSORPTION COMPARISON Pure CH 4, CO 2 nd ther mxture dsorpton on wet Frutlnd col [5] were used to evlute the models (the CO 2 dsorpton dt up to 8 MP were used). The mesurements were performed t 46.1 C wth pressures up to 12.8 MP. The feed molr compostons for the mxture dsorpton were 8%/2%, 6%/4%, 4%/6% nd 2%/8%.

5 Pn, Connell 5 The pure gs dsorpton dt re used to regress the model prmeters. The objectve functon, S, s used to correlte dt wth the dsorpton models. The functon mnmzes the sum of the squredpercentge devtons n predcted dsorpton: 2 NPTS c e ω ω e = 1 S ω = (22) NPTS c where ω nd e ω re the clculted nd expermentl dsorpton mount, respectvely. NPTS s the number of dt ponts. The percentge verge bsolute devton (%AAD) s used to evlute the results: e 1 ω ω % AAD = (23) NPTS NPTS c e = 1 ω The model evluton results for pure gs dsorpton re summrzed n Tble 1 nd the prmeters for the models re lsted n Tble 2. The model evluton results for the mxture dsorpton re lso lsted n Tble 1. Fgure 1 shows the CH 4 dsorpton results represented wth dfferent models. The Lngmur model underestmtes the pure CH 4 dsorpton t pressures hgher thn 7 MP nd slghtly overestmtes the pure CH 4 dsorpton from 2 MP to 6 MP. As result, both nd overestmte CH 4 component dsorpton below 6 MP nd slghtly overestmte the dsorpton t hgher pressure., on the other hnd, represents the pure CH 4 dsorpton more ccurtely s well s CH 4 component dsorpton. Fgure 2 shows the CO 2 dsorpton results usng dfferent models. In smlr fshon s wth the pure CH 4 dsorpton results, the Lngmur model underestmtes the pure CO 2 dsorpton t pressure hgher thn 6 MP nd slghtly overestmtes the pure CO 2 dsorpton below 6 MP. Snce the pure CO 2 dsorpton dt up to 8 MP were used n the model evluton, the models behvor t pressures hgher thn 8 MP re not known. However, for the CO 2 component dsorpton, predcts the dsorpton better thn or, especlly n the hgh pressure regon. Fgure 3 shows the totl dsorpton results usng dfferent models. Both nd underestmte the dsorpton for pressures hgher thn 7 MP. Fgure 4 shows the compostonl results for mxture dsorpton. predcts the dsorbed phse composton better thn the other two models. SIMULATION RESULTS SIMED II s colbed methne smultor cpble of compostonl dul-porosty reservor smulton. The dsorpton models ntegrted n SIMED II were nd. In ths work, the ws mplemented nto SIMED II. A hypothetcl one-dmensonl enhnced col methne through CO 2 njecton cse ws studed to llustrte how the dfferent dsorpton models ffect the reservor smulton results usng the dsorpton dt presented bove for Frutlnd col. The descrpton of the reservor propertes s lsted n Appendx A. Fgure 5 presents the gs rte clculted wth SIMED usng the dfferent dsorpton models. The CH 4 producton rte decreses fter rechng mxmum nd the dfferent dsorpton models predct smlr rtes tll round dy 4. After tht, the CH 4 producton rte ncreses nd predcts lower rte thn the two other models. Fgure 6 presents the dsorbed concentrtons for ech gs nd sotherm model. In Fgure 6(), t the erly stges of the smulton, the two wells hve no ntercton, thus the CH 4 producton well behves s conventonl prmry colbed methne recovery, where the methne s desorbed from the col by lowerng the sem pressure. Snce the representton of pure CH 4 dsorpton t pressures less thn 5 MP s lmost dentcl for the nd Lngmur models, the ntl

6 6 The Representton of Mult-component Adsorpton n Reservor Smulton of CO 2 Sequestrton n Col nd Enhnced Colbed Methne Recovery producton rte results for the two models re lso very smlr. However for Fgure 6(b) t 1 yer the CH 4 dsplced due to CO 2 njecton hs strted to rrve t the producton well nd the producton rte ncreses. The dsorbed concentrtons of CH 4 re hgher for the thn the other two models reflectng the dsorpton sotherm dfferences t these hgh pressures. As result the free gs volumes re less nd thus producton rte lower. Fgure 5(b) shows the gs rte results for 16 yer perod. Gs producton rte ncreses shrply fter dy 41, ndctng CO 2 brekthrough. Gs produced s mnly CO 2 fter dy 5. Fgure 6(e) shows tht the dsorbed CO 2 concentrton becomes hgher ner the producton well, ndctng comng CO 2 brekthrough. Fgure 6(f) shows no CH 4 dsorbed n the reservor nd CO 2 beng produced t dy 548. Fgure 5 lso presents CO 2 njecton rte. All dsorpton models consdered hve mxmum njecton rte t round dy 1. However, the CO 2 njecton rte predcted by s sgnfcntly hgher thn by ether or tll round dy 1. Then the njecton rte predcted by s slghtly hgher thn the other models. The lower CO 2 njecton rte t erly tmes by nd s due tht both models underestmte the CO 2 dsorpton t hgh pressures. Fgure 7 shows the gs phse composton profle usng nd. The results re smlr to the dsorbed gs concentrton becuse they re closely relted to ech other. Fgure 8 shows the cumultve CH 4 produced nd CO 2 njected for 2 yer perod. The results re consstent wth the gs rte results. SUMMARY nd CONCLUSIONS Ths pper hs consdered three pproches to representng dsorpton n reservor smulton of enhnced colbed methne through CO 2 sequestrton. Two of the dsorpton sotherm models, extended Lngmur nd Idel Adsorbed Soluton, re populr choces for reservor smulton. These provde computtonlly effcent clcultons but cn nvolve nccurces prtculrly for multcomponent gs mxtures nd sngle component CO 2. The sotherm model hs been consdered by severl uthors s n pproch to representng dsorpton. The sotherm model hs been shown to provde closer mtch to expermentl observtons on Frutlnd col for wde pressure rnge nd gs component mxture thn the extended Lngmur nd Idel Adsorbed Soluton sotherms. The sotherm model ws mplemented n the compostonl dul porosty smultor SIMED. A hypothetcl cse study ws formulted to nvestgte the reltve effects of the three sotherm models on reservor smulton of CO 2 sequestrton nd enhnced colbed methne recovery. In ths cse study CO 2 njecton well s seprted from CH 4 recovery well by 1m for col sem t n ntl reservor pressure of 8.8 MP, correspondng to roughly 88 m depth. Whle the njecton well opertes t pressures sgnfcntly hgher thn the ntl reservor pressure, the producton well produces methne ntlly through pressure drwdown below the desorpton pressure. In the cse study, njected CO 2 dsplces dsorbed methne whch then mgrtes wy from the njecton well nd towrds the producton well n response to the pressure grdent wthn the sem. In terms of dsorbed gses, three regons re creted wthn the sem; towrds the njecton well there s regon of pure CO 2, round the producton well there s regon of pure CH 4. In between these two regons there s zone of mxed CO 2 -CH 4, trnston regon. The sze of ths trnston regon s determned by the rte of dffuson nto the col mtrx; chrcterzed n the modellng by the desorpton tme constnt. Wthn ths trnston regon mxed gs dsorpton sotherm s requred. Intlly the methne producton rtes re very smlr for the three dsorpton models. Ths s consstent wth the ccurcy wth whch the dsorpton behvour s represented; for ths low pressure, sngle component methne regme ll three sotherm models re good descrptons of the dsorpton behvour. However wth tme ths good greement deterortes, wth the methne rtes clculted usng the beng sgnfcntly lower thn those clculted usng the other two dsorpton models. Ths

7 Pn, Connell 7 cn be ttrbuted to the ntercton of free methne dsplced by the CO 2 njecton rrvng t the producton well. As dscussed bove the CO 2 njecton cretes mxed gs regon where multcomponent dsorpton occurs. The s more ccurte descrpton of ths process thn the nd models. In ddton, the njecton of CO 2 nvolves hgh pressures; dsorpton of CO 2 for ths sngle component hgh pressure regme s lso more ccurtely descrbed usng the. The combnton of these two effects leds to ncresng dfferences n methne producton rte wth tme s the ntercton between the two wells becomes stronger. The most sgnfcnt dfferences ssocted wth the sotherm models re seen n the CO 2 njecton rte. The rte wth bsed clcultons s sgnfcntly hgher thn wth the other two dsorpton pproches, prtculrly t erly tmes. Ths behvour reflects the dfferences n the dsorpton sotherms, where t hgh pressures the CO 2 dsorbed mount, s represented by the, s greter thn tht estmted by the nd models. Snce the pressure behvour s very smlr cross the smultons, the njecton rte s therefore hgher to compenste for the greter mount dsorbed. After rechng mxmum the njecton rte decreses wth tme, n response to the pressure behvour wthn the reservor. Ths pper hs demonstrted tht mproved models for dsorpton cn hve sgnfcnt mpcts on reservor smulton of CO 2 sequestrton nd enhnced colbed methne. The mgntude of ths mpct s functon of pressure nd the role of mult-component gs mgrton. Implementton of the dsorpton sotherm model nto SIMED hs therefore mproved the cpblty of ths smultor to represent CO 2 sequestrton n col nd enhnced colbed methne producton. Appendx A. Descrpton of reservor propertes used n smulton cse study Grd system Crtesn x drecton: see Tble 3 Reservor temperture: 46.1 C Permeblty: 2.5 md Reltve permeblty: see Tble 4 Gs desorpton pressure: 5 MP Intl reservor pressure: 8 MP Desorpton tme constnt: 1 dys for CO 2, 27.5 dys for CH 4 Well loctons: Injecton well: Grd 1 Producton well: Grd 78 Well skn fctor= Bottomhole pressure: kp Injecton pressure: 14 MP

8 8 The Representton of Mult-component Adsorpton n Reservor Smulton of CO 2 Sequestrton n Col nd Enhnced Colbed Methne Recovery References 1. Do, D. D., 1998: Adsorpton Anlyss: Equlbr nd Knetcs: London, Imperl College Press. 2. Myers, A. L., Prusntz, J. M., 1965: Thermodynmcs of Mxed-Gs Adsorpton, AIChE J., V. 11, p Zhou, C., 1994: Modelng nd Predcton of Pure nd Multcomponent Gs Adsorpton, Ph.D Dssertton: Oklhom Stte Unversty, Stllwter, Oklhom. 4. Zhou, C., Hll, F., Gsem, K. A. M., Robnson, Jr., R. L., 1994: Predctng Gs Adsorpton Usng Two-Dmensonl Equtons of Stte, I&EC Reserch, V. 33, p Hll, F., Zhou, C., Gsem, K. A. M., Robnson, Jr., R. L., 1994: Adsorpton of Pure Methne, Ntrogen, nd Crbon Doxde nd Ther Bnry Mxtures on Wet Frutlnd Col, presented t the Estern Regonl Conference & Exhbton, Chrleston.

9 Pn, Connell 9 Tbles Tble 1. Adsorpton results - %AAD Pure Gs Adsorpton Lngmur CH CO Mxture Adsorpton CH CO Totl Tble 2. Best ft prmeter vlues for dsorpton sotherms of Frutlnd col Lngmur L (m 3 /ton) B (MP -1 ) CH CO α (br cm 3 g/mmol / mol) β (mmol/g) -1 -ln(k) ln(mmol/g br -1 ) CH CO Tble 3. Grd spcng used n smulton cse study Grd x(m) Grd x(m) Grd x(m) Grd x(m) * x[grd]= x[79-grd] when grd s 4 to 78 Tble 4. Reltve permeblty of wter nd gs used n smulton cse study Wter Sturton Wter phse reltve permeblty Gs phse reltve permeblty

10 1 The Representton of Mult-component Adsorpton n Reservor Smulton of CO 2 Sequestrton n Col nd Enhnced Colbed Methne Recovery Fgures Absolute Adsorpton (m 3 /ton) CH4 Exprmentl Feed Composton: Pure CH 4 8% CH 4 6% CH 4 4% CH 4 2% CH Pressure (MP) Fgure 1. CH 4 dsorpton on wet Frutlnd col t 46.1 C Absolute Adsorpton (m 3 /ton) CO2 Expermentl Feed Composton: Pure CO 2 8% CO 2 6% CO 2 4% CO 2 2% CO Pressure (MP) Fgure 2. CO 2 dsorpton on wet Frutlnd col t 46.1 C

11 Pn, Connell Expermentl Absolute Adsorpton (m 3 /ton) Pure CO 2 8% CO 2 6% CO 2 4% CO 2 2% CO 2 Pure CH Pressure (MP) 14 Fgure 3. Totl dsorpton on wet Frutlnd col t 46.1 C 1 Expermentl.8.6 Feed Composton: 8% CH 4 x CH4 Feed Composton: 6% CH 4.4 Feed Composton: 4% CH 4.2 Feed Composton: 2% CH 4 Low Pressure Hgh Pressure y CH4 Fgure 4. Composton for dsorpton on wet Frutlnd col t 46.1 C

12 12 The Representton of Mult-component Adsorpton n Reservor Smulton of CO 2 Sequestrton n Col nd Enhnced Colbed Methne Recovery Gs (m 3 /dy) 3 CO 2 Injecton Rte 2 1 CH 4 Producton Rte Tme (dy) () 2 yer (73 dys) perod 6 5 Gs (m 3 /dy) CO 2 Injecton Rte 1 Gs Producton Rte Tme (dy) (b) 16 yer (584 dys) perod Fgure 5. Gs rtes t njecton nd producton wells for the hypothetcl smulton exmple, clculted usng SIMED, for the three dsorpton sotherms

13 Pn, Connell Gs Concentrton (m 3 /tonne) CH4 Gs Concentrton (m 3 /tonne) CH () (b) Gs Concentrton (m 3 /tonne) CH4 Gs Concentrton (m 3 /tonne) CH (c) (d) 3 3 Gs Concentrton (m 3 /tonne) CH4 Gs Concentrton (m 3 /tonne) CH (e) (f) Fgure 6. Adsorbed gs concentrtons for SIMED smulton cse study t () Dy 182, (b) Dy 365, (c) Dy 547, (d) Dy 73, (e) Dy 365, (f) Dy 548

14 14 The Representton of Mult-component Adsorpton n Reservor Smulton of CO 2 Sequestrton n Col nd Enhnced Colbed Methne Recovery Gs Phse Composton CH4 Gs Phse Composton CH () (b) Gs Phse Composton CH4 Gs Phse Composton CH (c) Fgure 7. Gs phse composton for SIMED smulton cse study t () Dy 182, (b) Dy 195, (c) Dy 365, (d) Dy 4197 (d) Gs (m 3 ) 3 2 Gs (m 3 ) Tme (dy) () Tme (dy) (b) Fgure 8. Cumultve gs producton for smulton cse study () CH 4 (b) CO 2

Solubilities and Thermodynamic Properties of SO 2 in Ionic

Solubilities and Thermodynamic Properties of SO 2 in Ionic Solubltes nd Therodync Propertes of SO n Ionc Lquds Men Jn, Yucu Hou, b Weze Wu, *, Shuhng Ren nd Shdong Tn, L Xo, nd Zhgng Le Stte Key Lbortory of Checl Resource Engneerng, Beng Unversty of Checl Technology,

More information

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II Mcroeconomc Theory I UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS MSc n Economcs MICROECONOMIC THEORY I Techng: A Lptns (Note: The number of ndctes exercse s dffculty level) ()True or flse? If V( y )

More information

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Chapter Newton-Raphson Method of Solving a Nonlinear Equation Chpter.4 Newton-Rphson Method of Solvng Nonlner Equton After redng ths chpter, you should be ble to:. derve the Newton-Rphson method formul,. develop the lgorthm of the Newton-Rphson method,. use the Newton-Rphson

More information

Chemical Reaction Engineering

Chemical Reaction Engineering Lecture 20 hemcl Recton Engneerng (RE) s the feld tht studes the rtes nd mechnsms of chemcl rectons nd the desgn of the rectors n whch they tke plce. Lst Lecture Energy Blnce Fundmentls F 0 E 0 F E Q W

More information

Chemical Reaction Engineering

Chemical Reaction Engineering Lecture 20 hemcl Recton Engneerng (RE) s the feld tht studes the rtes nd mechnsms of chemcl rectons nd the desgn of the rectors n whch they tke plce. Lst Lecture Energy Blnce Fundmentls F E F E + Q! 0

More information

Investigation phase in case of Bragg coupling

Investigation phase in case of Bragg coupling Journl of Th-Qr Unversty No.3 Vol.4 December/008 Investgton phse n cse of Brgg couplng Hder K. Mouhmd Deprtment of Physcs, College of Scence, Th-Qr, Unv. Mouhmd H. Abdullh Deprtment of Physcs, College

More information

Electrochemical Thermodynamics. Interfaces and Energy Conversion

Electrochemical Thermodynamics. Interfaces and Energy Conversion CHE465/865, 2006-3, Lecture 6, 18 th Sep., 2006 Electrochemcl Thermodynmcs Interfces nd Energy Converson Where does the energy contrbuton F zϕ dn come from? Frst lw of thermodynmcs (conservton of energy):

More information

Name: SID: Discussion Session:

Name: SID: Discussion Session: Nme: SID: Dscusson Sesson: hemcl Engneerng hermodynmcs -- Fll 008 uesdy, Octoer, 008 Merm I - 70 mnutes 00 onts otl losed Book nd Notes (5 ponts). onsder n del gs wth constnt het cpctes. Indcte whether

More information

4. Eccentric axial loading, cross-section core

4. Eccentric axial loading, cross-section core . Eccentrc xl lodng, cross-secton core Introducton We re strtng to consder more generl cse when the xl force nd bxl bendng ct smultneousl n the cross-secton of the br. B vrtue of Snt-Vennt s prncple we

More information

Remember: Project Proposals are due April 11.

Remember: Project Proposals are due April 11. Bonformtcs ecture Notes Announcements Remember: Project Proposls re due Aprl. Clss 22 Aprl 4, 2002 A. Hdden Mrov Models. Defntons Emple - Consder the emple we tled bout n clss lst tme wth the cons. However,

More information

Applied Statistics Qualifier Examination

Applied Statistics Qualifier Examination Appled Sttstcs Qulfer Exmnton Qul_june_8 Fll 8 Instructons: () The exmnton contns 4 Questons. You re to nswer 3 out of 4 of them. () You my use ny books nd clss notes tht you mght fnd helpful n solvng

More information

Principle Component Analysis

Principle Component Analysis Prncple Component Anlyss Jng Go SUNY Bufflo Why Dmensonlty Reducton? We hve too mny dmensons o reson bout or obtn nsghts from o vsulze oo much nose n the dt Need to reduce them to smller set of fctors

More information

Supporting information

Supporting information upportng nformton for Towrds more useful n vtro toxcty dt wth mesured concentrtons by M.B. Herng R.H.M.M. chreurs F. Busser.T. vn der g B. vn der Burg nd J..M. Hermens Contnng 0 pges 2 tbles nd 2 fgures.

More information

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x) DCDM BUSINESS SCHOOL NUMEICAL METHODS (COS -8) Solutons to Assgnment Queston Consder the followng dt: 5 f() 8 7 5 () Set up dfference tble through fourth dfferences. (b) Wht s the mnmum degree tht n nterpoltng

More information

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Chapter Newton-Raphson Method of Solving a Nonlinear Equation Chpter 0.04 Newton-Rphson Method o Solvng Nonlner Equton Ater redng ths chpter, you should be ble to:. derve the Newton-Rphson method ormul,. develop the lgorthm o the Newton-Rphson method,. use the Newton-Rphson

More information

Quiz: Experimental Physics Lab-I

Quiz: Experimental Physics Lab-I Mxmum Mrks: 18 Totl tme llowed: 35 mn Quz: Expermentl Physcs Lb-I Nme: Roll no: Attempt ll questons. 1. In n experment, bll of mss 100 g s dropped from heght of 65 cm nto the snd contner, the mpct s clled

More information

The Schur-Cohn Algorithm

The Schur-Cohn Algorithm Modelng, Estmton nd Otml Flterng n Sgnl Processng Mohmed Njm Coyrght 8, ISTE Ltd. Aendx F The Schur-Cohn Algorthm In ths endx, our m s to resent the Schur-Cohn lgorthm [] whch s often used s crteron for

More information

4. More general extremum principles and thermodynamic potentials

4. More general extremum principles and thermodynamic potentials 4. More generl etremum prncples nd thermodynmc potentls We hve seen tht mn{u(s, X )} nd m{s(u, X)} mply one nother. Under certn condtons, these prncples re very convenent. For emple, ds = 1 T du T dv +

More information

Definition of Tracking

Definition of Tracking Trckng Defnton of Trckng Trckng: Generte some conclusons bout the moton of the scene, objects, or the cmer, gven sequence of mges. Knowng ths moton, predct where thngs re gong to project n the net mge,

More information

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC

Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC Introducton Rnk One Updte And the Google Mtrx y Al Bernsten Sgnl Scence, LLC www.sgnlscence.net here re two dfferent wys to perform mtrx multplctons. he frst uses dot product formulton nd the second uses

More information

6 Roots of Equations: Open Methods

6 Roots of Equations: Open Methods HK Km Slghtly modfed 3//9, /8/6 Frstly wrtten t Mrch 5 6 Roots of Equtons: Open Methods Smple Fed-Pont Iterton Newton-Rphson Secnt Methods MATLAB Functon: fzero Polynomls Cse Study: Ppe Frcton Brcketng

More information

Lecture 4: Piecewise Cubic Interpolation

Lecture 4: Piecewise Cubic Interpolation Lecture notes on Vrtonl nd Approxmte Methods n Appled Mthemtcs - A Perce UBC Lecture 4: Pecewse Cubc Interpolton Compled 6 August 7 In ths lecture we consder pecewse cubc nterpolton n whch cubc polynoml

More information

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1 Denns Brcker, 2001 Dept of Industrl Engneerng The Unversty of Iow MDP: Tx pge 1 A tx serves three djcent towns: A, B, nd C. Ech tme the tx dschrges pssenger, the drver must choose from three possble ctons:

More information

Many-Body Calculations of the Isotope Shift

Many-Body Calculations of the Isotope Shift Mny-Body Clcultons of the Isotope Shft W. R. Johnson Mrch 11, 1 1 Introducton Atomc energy levels re commonly evluted ssumng tht the nucler mss s nfnte. In ths report, we consder correctons to tomc levels

More information

Bi-level models for OD matrix estimation

Bi-level models for OD matrix estimation TNK084 Trffc Theory seres Vol.4, number. My 2008 B-level models for OD mtrx estmton Hn Zhng, Quyng Meng Abstrct- Ths pper ntroduces two types of O/D mtrx estmton model: ME2 nd Grdent. ME2 s mxmum-entropy

More information

LAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION IN A TWO-LAYERED SLAB

LAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION IN A TWO-LAYERED SLAB Journl of Appled Mthemtcs nd Computtonl Mechncs 5, 4(4), 5-3 www.mcm.pcz.pl p-issn 99-9965 DOI:.75/jmcm.5.4. e-issn 353-588 LAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION

More information

ψ ij has the eigenvalue

ψ ij has the eigenvalue Moller Plesset Perturbton Theory In Moller-Plesset (MP) perturbton theory one tes the unperturbed Hmltonn for n tom or molecule s the sum of the one prtcle Foc opertors H F() where the egenfunctons of

More information

Strategy: Use the Gibbs phase rule (Equation 5.3). How many components are present?

Strategy: Use the Gibbs phase rule (Equation 5.3). How many components are present? University Chemistry Quiz 4 2014/12/11 1. (5%) Wht is the dimensionlity of the three-phse coexistence region in mixture of Al, Ni, nd Cu? Wht type of geometricl region dose this define? Strtegy: Use the

More information

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Supplementary Notes for Chapter 9 Mixture Thermodynamics Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects

More information

523 P a g e. is measured through p. should be slower for lesser values of p and faster for greater values of p. If we set p*

523 P a g e. is measured through p. should be slower for lesser values of p and faster for greater values of p. If we set p* R. Smpth Kumr, R. Kruthk, R. Rdhkrshnn / Interntonl Journl of Engneerng Reserch nd Applctons (IJERA) ISSN: 48-96 www.jer.com Vol., Issue 4, July-August 0, pp.5-58 Constructon Of Mxed Smplng Plns Indexed

More information

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus ESI 34 tmospherc Dnmcs I Lesson 1 Vectors nd Vector lculus Reference: Schum s Outlne Seres: Mthemtcl Hndbook of Formuls nd Tbles Suggested Redng: Mrtn Secton 1 OORDINTE SYSTEMS n orthonorml coordnte sstem

More information

Chemistry 163B Absolute Entropies and Entropy of Mixing

Chemistry 163B Absolute Entropies and Entropy of Mixing Chemstry 163 Wnter 1 Hndouts for hrd Lw nd Entropy of Mxng (del gs, dstngushle molecules) PPENDIX : H f, G f, U S (no Δ, no su f ) Chemstry 163 solute Entropes nd Entropy of Mxng Hº f Gº f Sº 1 hrd Lw

More information

GAUSS ELIMINATION. Consider the following system of algebraic linear equations

GAUSS ELIMINATION. Consider the following system of algebraic linear equations Numercl Anlyss for Engneers Germn Jordnn Unversty GAUSS ELIMINATION Consder the followng system of lgebrc lner equtons To solve the bove system usng clsscl methods, equton () s subtrcted from equton ()

More information

6. Chemical Potential and the Grand Partition Function

6. Chemical Potential and the Grand Partition Function 6. Chemcl Potentl nd the Grnd Prtton Functon ome Mth Fcts (see ppendx E for detls) If F() s n nlytc functon of stte vrles nd such tht df d pd then t follows: F F p lso snce F p F we cn conclude: p In other

More information

INTRODUCTION TO COMPLEX NUMBERS

INTRODUCTION TO COMPLEX NUMBERS INTRODUCTION TO COMPLEX NUMBERS The numers -4, -3, -, -1, 0, 1,, 3, 4 represent the negtve nd postve rel numers termed ntegers. As one frst lerns n mddle school they cn e thought of s unt dstnce spced

More information

International Journal of Pure and Applied Sciences and Technology

International Journal of Pure and Applied Sciences and Technology Int. J. Pure Appl. Sc. Technol., () (), pp. 44-49 Interntonl Journl of Pure nd Appled Scences nd Technolog ISSN 9-67 Avlle onlne t www.jopst.n Reserch Pper Numercl Soluton for Non-Lner Fredholm Integrl

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

Module 2: Rate Law & Stoichiomtery (Chapter 3, Fogler)

Module 2: Rate Law & Stoichiomtery (Chapter 3, Fogler) CHE 309: Chemicl Rection Engineering Lecture-8 Module 2: Rte Lw & Stoichiomtery (Chpter 3, Fogler) Topics to be covered in tody s lecture Thermodynmics nd Kinetics Rection rtes for reversible rections

More information

6.6 The Marquardt Algorithm

6.6 The Marquardt Algorithm 6.6 The Mqudt Algothm lmttons of the gdent nd Tylo expnson methods ecstng the Tylo expnson n tems of ch-sque devtves ecstng the gdent sech nto n tetve mtx fomlsm Mqudt's lgothm utomtclly combnes the gdent

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

Katholieke Universiteit Leuven Department of Computer Science

Katholieke Universiteit Leuven Department of Computer Science Updte Rules for Weghted Non-negtve FH*G Fctorzton Peter Peers Phlp Dutré Report CW 440, Aprl 006 Ktholeke Unverstet Leuven Deprtment of Computer Scence Celestjnenln 00A B-3001 Heverlee (Belgum) Updte Rules

More information

CHM Physical Chemistry I Chapter 1 - Supplementary Material

CHM Physical Chemistry I Chapter 1 - Supplementary Material CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion

More information

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. with respect to λ. 1. χ λ χ λ ( ) λ, and thus:

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. with respect to λ. 1. χ λ χ λ ( ) λ, and thus: More on χ nd errors : uppose tht we re fttng for sngle -prmeter, mnmzng: If we epnd The vlue χ ( ( ( ; ( wth respect to. χ n Tlor seres n the vcnt of ts mnmum vlue χ ( mn χ χ χ χ + + + mn mnmzes χ, nd

More information

Physics 121 Sample Common Exam 2 Rev2 NOTE: ANSWERS ARE ON PAGE 7. Instructions:

Physics 121 Sample Common Exam 2 Rev2 NOTE: ANSWERS ARE ON PAGE 7. Instructions: Physcs 121 Smple Common Exm 2 Rev2 NOTE: ANSWERS ARE ON PAGE 7 Nme (Prnt): 4 Dgt ID: Secton: Instructons: Answer ll 27 multple choce questons. You my need to do some clculton. Answer ech queston on the

More information

Rel Gses 1. Gses (N, CO ) which don t obey gs lws or gs eqution P=RT t ll pressure nd tempertures re clled rel gses.. Rel gses obey gs lws t extremely low pressure nd high temperture. Rel gses devited

More information

Introduction to Numerical Integration Part II

Introduction to Numerical Integration Part II Introducton to umercl Integrton Prt II CS 75/Mth 75 Brn T. Smth, UM, CS Dept. Sprng, 998 4/9/998 qud_ Intro to Gussn Qudrture s eore, the generl tretment chnges the ntegrton prolem to ndng the ntegrl w

More information

LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER

LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER Yn, S.-P.: Locl Frctonl Lplce Seres Expnson Method for Dffuson THERMAL SCIENCE, Yer 25, Vol. 9, Suppl., pp. S3-S35 S3 LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN

More information

Statistics 423 Midterm Examination Winter 2009

Statistics 423 Midterm Examination Winter 2009 Sttstcs 43 Mdterm Exmnton Wnter 009 Nme: e-ml: 1. Plese prnt your nme nd e-ml ddress n the bove spces.. Do not turn ths pge untl nstructed to do so. 3. Ths s closed book exmnton. You my hve your hnd clcultor

More information

PLEASE SCROLL DOWN FOR ARTICLE

PLEASE SCROLL DOWN FOR ARTICLE Ths rtcle ws downloded by:ntonl Cheng Kung Unversty] On: 1 September 7 Access Detls: subscrpton number 7765748] Publsher: Tylor & Frncs Inform Ltd Regstered n Englnd nd Wles Regstered Number: 17954 Regstered

More information

8. INVERSE Z-TRANSFORM

8. INVERSE Z-TRANSFORM 8. INVERSE Z-TRANSFORM The proce by whch Z-trnform of tme ere, nmely X(), returned to the tme domn clled the nvere Z-trnform. The nvere Z-trnform defned by: Computer tudy Z X M-fle trn.m ued to fnd nvere

More information

18.7 Artificial Neural Networks

18.7 Artificial Neural Networks 310 18.7 Artfcl Neurl Networks Neuroscence hs hypotheszed tht mentl ctvty conssts prmrly of electrochemcl ctvty n networks of brn cells clled neurons Ths led McCulloch nd Ptts to devse ther mthemtcl model

More information

Lecture 11. Chapter 7. - Thermodynamic Web - Departure Functions - Review Equations of state (chapter 4, briefly)

Lecture 11. Chapter 7. - Thermodynamic Web - Departure Functions - Review Equations of state (chapter 4, briefly) Lecture 11 Chpter 5 - Thermodnmc We - Deprture unctons - Revew Equtons of stte (chpter 4, refl) Chpter 6 - Equlrum (chemcl potentl) * ure Component * Mtures Chpter 7 - ugct (chemcl potentl fugct equlrum

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

Tilted Plane Symmetric Magnetized Cosmological Models

Tilted Plane Symmetric Magnetized Cosmological Models Tlted Plne Symmetrc Mgnetzed Cosmologcl Models D. D. Pwr # *, V. J. Dgwl @ & Y. S. Solnke & # School of Mthemtcl Scences, Swm Rmnnd Teerth Mrthwd Unversty, Vshnupur, Nnded-0, (Ind) @ Dept. of Mthemtcs,

More information

1 Introduction. 2 ph Calculation

1 Introduction. 2 ph Calculation 1 Introducton The ph, lklnty nd totl norgnc crbon lgorthms tht re ncorported nto WASP come drectly from QUALK & QUALKw (Chpr, Pelleter nd To, 008. Every effort ws mde to nsure tht the mplementton of the

More information

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk Chpter 5 Supplementl Text Mterl 5-. Expected Men Squres n the Two-fctor Fctorl Consder the two-fctor fxed effects model y = µ + τ + β + ( τβ) + ε k R S T =,,, =,,, k =,,, n gven s Equton (5-) n the textook.

More information

Online Appendix to. Mandating Behavioral Conformity in Social Groups with Conformist Members

Online Appendix to. Mandating Behavioral Conformity in Social Groups with Conformist Members Onlne Appendx to Mndtng Behvorl Conformty n Socl Groups wth Conformst Members Peter Grzl Andrze Bnk (Correspondng uthor) Deprtment of Economcs, The Wllms School, Wshngton nd Lee Unversty, Lexngton, 4450

More information

Jens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers

Jens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers Jens Sebel (Unversty of Appled Scences Kserslutern) An Interctve Introducton to Complex Numbers 1. Introducton We know tht some polynoml equtons do not hve ny solutons on R/. Exmple 1.1: Solve x + 1= for

More information

Activator-Inhibitor Model of a Dynamical System: Application to an Oscillating Chemical Reaction System

Activator-Inhibitor Model of a Dynamical System: Application to an Oscillating Chemical Reaction System Actvtor-Inhtor Model of Dynmcl System: Applcton to n Osclltng Chemcl Recton System C.G. Chrrth*P P,Denn BsuP P * Deprtment of Appled Mthemtcs Unversty of Clcutt 9, A. P. C. Rod, Kolt-79 # Deprtment of

More information

13 Design of Revetments, Seawalls and Bulkheads Forces & Earth Pressures

13 Design of Revetments, Seawalls and Bulkheads Forces & Earth Pressures 13 Desgn of Revetments, Sewlls nd Bulkheds Forces & Erth ressures Ref: Shore rotecton Mnul, USACE, 1984 EM 1110--1614, Desgn of Revetments, Sewlls nd Bulkheds, USACE, 1995 Brekwters, Jettes, Bulkheds nd

More information

INTERPOLATION(1) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

INTERPOLATION(1) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek ELM Numercl Anlss Dr Muhrrem Mercmek INTEPOLATION ELM Numercl Anlss Some of the contents re dopted from Lurene V. Fusett, Appled Numercl Anlss usng MATLAB. Prentce Hll Inc., 999 ELM Numercl Anlss Dr Muhrrem

More information

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism CS294-40 Lernng for Rootcs nd Control Lecture 10-9/30/2008 Lecturer: Peter Aeel Prtlly Oservle Systems Scre: Dvd Nchum Lecture outlne POMDP formlsm Pont-sed vlue terton Glol methods: polytree, enumerton,

More information

Solution of Tutorial 5 Drive dynamics & control

Solution of Tutorial 5 Drive dynamics & control ELEC463 Unversty of New South Wles School of Electrcl Engneerng & elecommunctons ELEC463 Electrc Drve Systems Queston Motor Soluton of utorl 5 Drve dynmcs & control 500 rev/mn = 5.3 rd/s 750 rted 4.3 Nm

More information

ORDINARY DIFFERENTIAL EQUATIONS

ORDINARY DIFFERENTIAL EQUATIONS 6 ORDINARY DIFFERENTIAL EQUATIONS Introducton Runge-Kutt Metods Mult-step Metods Sstem o Equtons Boundr Vlue Problems Crcterstc Vlue Problems Cpter 6 Ordnr Derentl Equtons / 6. Introducton In mn engneerng

More information

Part I: Basic Concepts of Thermodynamics

Part I: Basic Concepts of Thermodynamics Prt I: Bsic Concepts o Thermodynmics Lecture 4: Kinetic Theory o Gses Kinetic Theory or rel gses 4-1 Kinetic Theory or rel gses Recll tht or rel gses: (i The volume occupied by the molecules under ordinry

More information

Statistics and Probability Letters

Statistics and Probability Letters Sttstcs nd Probblty Letters 79 (2009) 105 111 Contents lsts vlble t ScenceDrect Sttstcs nd Probblty Letters journl homepge: www.elsever.com/locte/stpro Lmtng behvour of movng verge processes under ϕ-mxng

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

APPROXIMATE INTEGRATION

APPROXIMATE INTEGRATION APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be

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

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero

More information

Characteristics of Fiber Suspension Flow in a Turbulent Boundary Layer

Characteristics of Fiber Suspension Flow in a Turbulent Boundary Layer Chrcterstcs of Fber Suspenson Flow n Turbulent Boundry Lyer Jnzhong Ln, Suhu Shen, Xoke Ku Deprtment of Mechncs Hngzhou, Zheng CHINA Correspondence to: Jnzhong Ln eml: meczln@publc.zu.edu.cn ABSTRACT The

More information

ANALYTICAL AND EXPERIMENTAL STUDY OF CONICAL TELESCOPING SPRINGS WITH NON-CONSTANT PITCH

ANALYTICAL AND EXPERIMENTAL STUDY OF CONICAL TELESCOPING SPRINGS WITH NON-CONSTANT PITCH ANALYTICAL AND EXPEIMENTAL STUDY OF CONICAL TELESCOPING SPINGS WITH NON-CONSTANT PITCH Mnuel Predes Unversté de Toulouse; INSA, UPS, Mnes Alb, ISAE; ICA (Insttut Clément Ader) 5, venue de nguel, F-077

More information

NUMERICAL MODELLING OF A CILIUM USING AN INTEGRAL EQUATION

NUMERICAL MODELLING OF A CILIUM USING AN INTEGRAL EQUATION NUEICAL ODELLING OF A CILIU USING AN INTEGAL EQUATION IHAI EBICAN, DANIEL IOAN Key words: Cl, Numercl nlyss, Electromgnetc feld, gnetton. The pper presents fst nd ccurte method to model the mgnetc behvour

More information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es

More information

Least squares. Václav Hlaváč. Czech Technical University in Prague

Least squares. Václav Hlaváč. Czech Technical University in Prague Lest squres Václv Hlváč Czech echncl Unversty n Prgue hlvc@fel.cvut.cz http://cmp.felk.cvut.cz/~hlvc Courtesy: Fred Pghn nd J.P. Lews, SIGGRAPH 2007 Course; Outlne 2 Lner regresson Geometry of lest-squres

More information

The Study of Lawson Criterion in Fusion Systems for the

The Study of Lawson Criterion in Fusion Systems for the Interntonl Archve of Appled Scences nd Technology Int. Arch. App. Sc. Technol; Vol 6 [] Mrch : -6 Socety of ducton, Ind [ISO9: 8 ertfed Orgnzton] www.soeg.co/st.html OD: IAASA IAAST OLI ISS - 6 PRIT ISS

More information

Title. Phosphate Lattice Loss Simulation. Authors

Title. Phosphate Lattice Loss Simulation. Authors Ttle Phosphte Lttce Loss Smulton Authors Mohmmd Abutyeh, Deprtment of Chemcl Engneerng, Unversty of South Flord, butyeh@ml.usf.edu Scott W. Cmpbell, Deprtment of Chemcl Engneerng, Unversty of South Flord,

More information

NUMERICAL MODELLING OF TWO-DIMENSIONAL HEAVE FOR SLABS- ON-GROUND AND SHALLOW FOUNDATIONS

NUMERICAL MODELLING OF TWO-DIMENSIONAL HEAVE FOR SLABS- ON-GROUND AND SHALLOW FOUNDATIONS 56 TH CANADIAN GEOTECHNICAL CONFERENCE 4 TH JOINT IAH-CNC/CGS CONFERENCE 2003 NAGS CONFERENCE 56 ème CONFÉRENCE CANADIENNE DE GÉOTECHNIQUE 4 ème CONFÉRENCE CONJOINTE AIH-CCN/SCG 2003 NAGS CONFÉRENCE NUMERICAL

More information

Study and modeling on saponification dynamics of the mixture of insect wax and oil-tea camellia seed oil

Study and modeling on saponification dynamics of the mixture of insect wax and oil-tea camellia seed oil Avlble onlne www.jocpr.com Journl of Chemcl nd Phrmceutcl Reserch, 04, 6(4):568-574 Reserch Artcle ISSN : 0975-7384 CODEN(USA) : JCPRC5 Study nd modelng on sponfcton dynmcs of the mxture of nsect wx nd

More information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

Work and Energy (Work Done by a Varying Force)

Work and Energy (Work Done by a Varying Force) Lecture 1 Chpter 7 Physcs I 3.5.14 ork nd Energy (ork Done y Vryng Force) Course weste: http://fculty.uml.edu/andry_dnylov/techng/physcsi Lecture Cpture: http://echo36.uml.edu/dnylov13/physcs1fll.html

More information

Chapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations

Chapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations Cter. Runge-Kutt nd Order Metod or Ordnr Derentl Eutons Ater redng ts cter ou sould be ble to:. understnd te Runge-Kutt nd order metod or ordnr derentl eutons nd ow to use t to solve roblems. Wt s te Runge-Kutt

More information

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems:

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems: TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases

More information

Concept of Activity. Concept of Activity. Thermodynamic Equilibrium Constants [ C] [ D] [ A] [ B]

Concept of Activity. Concept of Activity. Thermodynamic Equilibrium Constants [ C] [ D] [ A] [ B] Conept of Atvty Equlbrum onstnt s thermodynm property of n equlbrum system. For heml reton t equlbrum; Conept of Atvty Thermodynm Equlbrum Constnts A + bb = C + dd d [C] [D] [A] [B] b Conentrton equlbrum

More information

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede Fll Anlss of Epermentl Mesurements B. Esensten/rev. S. Errede Monte Crlo Methods/Technques: These re mong the most powerful tools for dt nlss nd smulton of eperments. The growth of ther mportnce s closel

More information

The heat budget of the atmosphere and the greenhouse effect

The heat budget of the atmosphere and the greenhouse effect The het budget of the tmosphere nd the greenhouse effect 1. Solr rdition 1.1 Solr constnt The rdition coming from the sun is clled solr rdition (shortwve rdition). Most of the solr rdition is visible light

More information

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

SUMMER KNOWHOW STUDY AND LEARNING CENTRE SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18

More information

Demand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert

Demand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert Demnd Demnd nd Comrtve Sttcs ECON 370: Mcroeconomc Theory Summer 004 Rce Unversty Stnley Glbert Usng the tools we hve develoed u to ths ont, we cn now determne demnd for n ndvdul consumer We seek demnd

More information

Lecture 5 Single factor design and analysis

Lecture 5 Single factor design and analysis Lectue 5 Sngle fcto desgn nd nlss Completel ndomzed desgn (CRD Completel ndomzed desgn In the desgn of expements, completel ndomzed desgns e fo studng the effects of one pm fcto wthout the need to tke

More information

I wish to publish my paper on The International Journal of Thermophysics. A Practical Method to Calculate Partial Properties from Equation of State

I wish to publish my paper on The International Journal of Thermophysics. A Practical Method to Calculate Partial Properties from Equation of State I wsh to publsh my paper on The Internatonal Journal of Thermophyscs. Ttle: A Practcal Method to Calculate Partal Propertes from Equaton of State Authors: Ryo Akasaka (correspondng author) 1 and Takehro

More information

Mass Transfer Processes

Mass Transfer Processes Mass Transfer Processes S. Majd Hassanzadeh Department of Earth Scences Faculty of Geoscences Utrecht Unversty Outlne: 1. Measures of Concentraton 2. Volatlzaton and Dssoluton 3. Adsorpton Processes 4.

More information

Equation of State Modeling of Phase Equilibrium in the Low-Density Polyethylene Process

Equation of State Modeling of Phase Equilibrium in the Low-Density Polyethylene Process Equaton of State Modelng of Phase Equlbrum n the Low-Densty Polyethylene Process H. Orbey, C. P. Boks, and C. C. Chen Ind. Eng. Chem. Res. 1998, 37, 4481-4491 Yong Soo Km Thermodynamcs & Propertes Lab.

More information

A Regression-Based Approach for Scaling-Up Personalized Recommender Systems in E-Commerce

A Regression-Based Approach for Scaling-Up Personalized Recommender Systems in E-Commerce A Regresson-Bsed Approch for Sclng-Up Personlzed Recommender Systems n E-Commerce Slobodn Vucetc 1 nd Zorn Obrdovc 1, svucetc@eecs.wsu.edu, zorn@cs.temple.edu 1 Electrcl Engneerng nd Computer Scence, Wshngton

More information

Modeling of CO 2 Cut in CBM Production

Modeling of CO 2 Cut in CBM Production Modelng of CO Cut n CBM Producton A study by: - Kamal Morad - Davd Dunn - ous Mattar Fekete Assocates Inc. Mult-Component Gas n Coalbed Methane Coalbed Methane has dfferent characterstcs from conventonal

More information

Computing a complete histogram of an image in Log(n) steps and minimum expected memory requirements using hypercubes

Computing a complete histogram of an image in Log(n) steps and minimum expected memory requirements using hypercubes Computng complete hstogrm of n mge n Log(n) steps nd mnmum expected memory requrements usng hypercubes TAREK M. SOBH School of Engneerng, Unversty of Brdgeport, Connectcut, USA. Abstrct Ths work frst revews

More information

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:

More information

Fitting a Polynomial to Heat Capacity as a Function of Temperature for Ag. Mathematical Background Document

Fitting a Polynomial to Heat Capacity as a Function of Temperature for Ag. Mathematical Background Document Fttng Polynol to Het Cpcty s Functon of Teperture for Ag. thetcl Bckground Docuent by Theres Jul Zelnsk Deprtent of Chestry, edcl Technology, nd Physcs onouth Unversty West ong Brnch, J 7764-898 tzelns@onouth.edu

More information

Dynamic Power Management in a Mobile Multimedia System with Guaranteed Quality-of-Service

Dynamic Power Management in a Mobile Multimedia System with Guaranteed Quality-of-Service Dynmc Power Mngement n Moble Multmed System wth Gurnteed Qulty-of-Servce Qnru Qu, Qng Wu, nd Mssoud Pedrm Dept. of Electrcl Engneerng-Systems Unversty of Southern Clforn Los Angeles CA 90089 Outlne! Introducton

More information

Predict Global Earth Temperature using Linier Regression

Predict Global Earth Temperature using Linier Regression Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id

More information

Evaluation of Liquefaction Return Period for Bangalore Based on Standard Penetration Test Data: Performance Based Approach

Evaluation of Liquefaction Return Period for Bangalore Based on Standard Penetration Test Data: Performance Based Approach Amercn J. of Engneerng nd Appled Scences 2 (3): 537-543, 2009 ISSN 1941-7020 2009 Scence Publctons Evluton of Lquefcton Return Perod for Bnglore Bsed on Stndrd Penetrton Test Dt: Performnce Bsed Approch

More information