Bubble and dew point calculations in multicomponent and multireactive mixtures
|
|
- Benedict Jordan
- 6 years ago
- Views:
Transcription
1 Insttuto Tecnologco de Aguascalentes From the SelectedWorks of Adran Bonlla-Petrcolet 2006 Bubble and dew pont calculatons n multcomponent and multreactve mtures Adran Bonlla-Petrcolet Azucena Acosta-Martínez Ulses I. Bravo-Sánchez Juan Gabrel Segova-Hernández, Unversdad de Guanauato Avalable at:
2 A. BONILLA-PETICIOLET et al., Bubble and Dew Pont Calculatons n, Chem. Bochem. Eng. Q. 20 (2) 8 (2006) BubbleandDewPontCalculatonsnMultcomponentandMultreactveMtures A. Bonlla-Petrcolet,* A. Acosta-Martínez, U. I. Bravo-Sánchez, andj. G. Segova-Hernández** Insttuto Tecnológco de Aguascalentes, Departamento de Ingenería Químca, Av. López Mateos 80, CP 20256, Aguascalentes, Ags. Méco Orgnal scentfc paper **Unversdad de Guanauato, Facultad de Cencas Químcas, eceved: August 28, 2005 Nora Alta s/n, C.P. 6050, Guanauato, Gto. Méco Accepted: February, 2006 Bubble and dew pont calculatons are useful n chemcal engneerng and play an mportant role n the study of separaton equpments for non-reactve and reactve mtures. To the best of the authors s knowledge, few methods have been proposed for these calculatons n systems wth several chemcal reactons. The obectve of ths paper s to ntroduce new condtons for performng bubble and dew pont calculatons n reactve mtures. We have developed these condtons based on the applcaton of transformed varables of Ung and Doherty (995). Usng these transformed varables, the soluton space s restrcted to compostons that are already at chemcal equlbrum and by consequence the problem dmenson s also reduced. The relablty and effcency of three equaton-solvng methods are tested and compared usng our equlbrum condtons: a) a smultaneous equaton-solvng approach usng Newton method (SESN), b) an equaton-decouplng approach usng successve substtuton method (EDSS) and c) an optmzaton approach usng the stochastc optmzaton method Smulated Annealng (OSA). Our results ndcated that even for smple reactve systems, bubble and dew pont calculatons are challengng for classcal equaton-solvng methods and requre robust strateges. We conclude that OSA and EDSS methods are relable to locate bubble and dew ponts n reactve systems. EDSS s more effcent than OSA; however, OSA does not need ntal guesses and s more sutable for dffcult problems. Key words: Chemcal equlbrum, phase equlbrum, bubble pont, dew pont, global optmzaton, equaton-solvng method comparson Introducton The descrpton of phase equlbrum n multreacton mtures s an mportant topc nvolved n several ndustral applcatons and t s the bass for desgn, analyss, and synthess of separaton processes. The phase equlbrum problem wth several chemcal reactons s hghly non lnear and multvarable. 2 Durng the last years, there has been a growng nterest for developng new methods to treatment the thermodynamc behavor of mtures under physcal and chemcal equlbrum. 2 Almost all developed methods use the mole fractons and reacton etents as ndependent varables and unknowns of the reactve phase equlbrum problem. Untl now, only a few methods have used theores of varable transformaton wth the am of reducng problem dmensonalty and mprovng the numercal behavor (effcency and relablty) of soluton methods.,9 These knds of algorthms are very attractve for the smulaton of separaton process and favor the study of comple multreactve multcomponent systems. *Correspondng author: (52) et. 27, petrcolet@hotmal.com A specal case of phase equlbrum problems are the bubble and dew pont calculatons. These calculatons are useful n chemcal engneerng and play an mportant role n the study of separaton equpments. 4 Under ths contet, robust, and effcent methods for these calculatons are desrable. To the best of the authors s knowledge, few methods have been developed to calculate bubble and dew ponts n multreactve mtures. 9,2 Ths paper ntroduces a new method for performng dew and bubble pont calculatons n multreactve and multcomponent systems. We have proposed new condtons for these equlbrum problems based on the theory of reactve varables of Ung and Doherty. Also, we have tested and compared three equaton-solvng methods for performng those calculatons, usng our equlbrum condtons. Formulaton of new condtons for bubble and dew pont calculatons n multreactve mtures Ung and Doherty showed that the Gbbs energy n a reactve system behaves as n a non-reac-
3 2 A. BONILLA-PETICIOLET et al., Bubble and Dew Pont Calculatons n, Chem. Bochem. Eng. Q. 20 (2) 8 (2006) tve system f transformed composton varables are used nstead of the conventonal composton varables. Usng transformed compostons, we restrct the soluton space to the compostons that are already n chemcal equlbrum. Ths reduces the problem dmenson by the number of ndependent reactons and makes t possble to represent phase equlbrum n reactve systems n a smlar way as n non-reactve mtures. So, the reactve phase dagrams look smlar to the non-reactve ones and the non-reactve flash algorthms can be easly modfed to account for the equlbrum reactons. 2 Also, Ung and Doherty showed that the chemcal potental follows all the thermodynamc relatonshps of a non-reactve system as long as all the thermodynamc propertes are functons of the transformed composton varables. Consderng ths, equlbrum condtons among two phases of a reactve mture, wth n components and n ndependent chemcal reactons, are met when where T,, nn () T {} 0 (2) T T wth the unversal gas constant, 0 the chemcal potental of the pure component and the transformed chemcal potental of component at phase whch s a functon of the transformed mole fracton. Eq (2) can be epressed n terms of fugacty or actvty coeffcents as follows {} ln ln ( y{}),, nn () T where s the mole fracton of component, s the fugacty coeffcent of component n the mture, s the fugacty coeffcent of the pure component, and s the actvty coeffcent of component, respectvely. Usng Eq () and (), we can deduce the followng relaton for mole fractons at phase equlbrum n n,, (4) Then, we defne the transformed phase equlbrum constant K as {, K },, nn (5) and the transformed mole fractons by 8,9, are gven N v 0,, nn (6) v N where s the column vector of n erence component mole fractons, v s the row vector of stochometrc number of component for each reacton v, s a row vector where each element corresponds to reacton r and t s the sum of the stochometrc number for all components that partcpate n reacton r, and N s a square matr formed from the stochometrc number of the erence components n the n reactons. 8,9, The erence mole fractons are calculated usng Eq (6) and from the equlbrum constants for each reacton K eq by solvng a system of n nonlnear equatons gven by n K a r,, n (7) eq v r where a s the actvty of component and v r s the stochometrc number of component n reacton r, respectvely. When we know the erence mole fractons, the remanng mole fractons are calculated usng Eq (6). We wrte the phase equlbrum condton n terms of transformed varables, usng Eq (4) (6), as follows and K,, nn (8) v v vn ( K ) v N N N (9),, nn (0) The materal balance on component s F. ( ) 0,, nn () and the transformed mole phase fracton s defned as ( v N ) ( v N ) F, (2) where s the conventonal mole fracton of phase whose feasble doman s (0, ) and F s the column vector of n erence component mole
4 A. BONILLA-PETICIOLET et al., Bubble and Dew Pont Calculatons n, Chem. Bochem. Eng. Q. 20 (2) 8 (2006) fractons n the ntal overall composton. Usng Eq (8) (), s gven by nn, F, ( K ),, nn () From the transformed mole fracton restrcton, we have nn K( F, ) 0 (4) K nn F, ( K 0 (5) By applyng the condtons 0 and, we obtan the followng new functons for bubble and dew pont calculatons n multreactve mtures nn f ( K ) 0 (6) bubble f dew F, nn F 0 K (7) In contrast wth other formulatons,,5 our functons are ndependent of reacton etents and have fewer unknowns. Eq (6) and (7) are a functon of temperature or pressure and n n transformed mole fracton of the lqud or vapor phase, respectvely. On the other hand, these equlbrum problems are nonlnear and multvarable so that conventonal numercal methods may present dffcultes, such as poor ntalzaton or dvergence behavor. We have performed a comparson of three approaches to solve our bubble and dew pont functons. In the net secton we descrbe the strateges appled n ths work. Soluton approachs Smultaneous equaton-solvng wth Newton method (SESN) Ths s the classcal approach used for performng flash calculatons n non-reactve and reactve systems snce t s conceptually smple and straghtforward. The popular Newton method s used because t provdes quadratc convergence when the ntal estmates are close to the soluton and t s readly avalable n computer programs. 6 Bubble and dew ponts are calculated by solvng smultaneously the followng n n equatons F K,, nn (8) Fnn f (9) where f s the bubble or dew pont functon, Eq (6) or (7), respectvely. The unknowns are n n transformed mole fractons of vapor or lqud phase and the temperature or pressure. In ths work, the above equatons are solved usng the Numercal ecpes Fortran subroutne NEWT. Equaton-decouplng wth Successve Substtuton method (EDSS) Ths approach s often descrbed n tetbooks on thermodynamcs for flash calculatons n non-reactve mtures and s based on the achford-ce formulaton. Ths method s relable for flash calculatons n non-reactve systems and has a lnear rate of convergence. 6 The procedure used for bubble and dew pont calculatons n reactve systems s as follows: n the nner loop, we have appled a Newton method, usng the Numercal ecpes Fortran subroutne NEWT, to update the temperature or pressure by solvng Eq (6) or (7), whle n the outer loop the transformed mole fractons are updated wth a successve substtuton method usng Eq (8) or (). In ths work, we have not appled an acceleraton technque for ths approach. Optmzaton approach usng Smulated Annealng method (OSA) We have tested an optmzaton technque for solvng our equlbrum condtons. As ndcated by Henderson et al., 7 the formulaton of thermodynamc calculatons for optmzaton problems offers some advantages: a) the use of a robust optmzaton method, b) the possblty of usng a drect optmzaton method whch requres only calculatons of the obectve functon and c) the use of an teratve procedure whose convergence s ndependent on the ntal guess. Thus, the calculaton of bubble and dew pont condtons can be performed by mnmzng the net obectve functon nn 2 2 f ob f F (20) where F s descrbed by Eq (8) for =,, n n and f s the bubble or dew pont functon, respectvely. At the bubble and dew pont condtons, the global mnmum of the obectve functon must be zero, but we assume that a soluton s found when we fnd the transformed mole fractons and temperature or pressure that make the value of the
5 4 A. BONILLA-PETICIOLET et al., Bubble and Dew Pont Calculatons n, Chem. Bochem. Eng. Q. 20 (2) 8 (2006) obectve functon less than 0 6. The stochastc global optmzaton method Smulated Annealng (SA) s used to mnmze our obectve functon. SA s nspred n the thermodynamc process of coolng of molten metals to acheve the lowest free energy state. 8 The SA algorthm presents ease of computatonal mplementaton and, f the values for ts algorthm parameters are properly selected, t can converge to the global optmum ndependently of ntal guesses. Ths optmzaton method has been successfully used n the resoluton of several thermodynamc problems. 7,9 24 It s mportant to note that other phase equlbrum problems (e.g. calculatons of crtcal ponts and homogeneous azeotropes) have been solved usng an optmzaton approach coupled wth the SA method. 7,24 In our calculatons, we have used the algorthm proposed by Corana et al. 25 because of ts hgh relablty n thermodynamc calculatons. In ths algorthm, a tral pont s randomly chosen wthn the step length VM (a vector of length n varables) of a startng pont defned n the feasble doman of optmzaton varables. The functon s evaluated at ths tral pont and ts value s compared to ts value at the ntal pont. The Metropols crteron, 26 wth a parameter called annealng temperature T, s used to accept or reect the tral pont. If the tral pont s accepted, the algorthm moves on from that pont. If t s reected, another pont s chosen nstead for a tral evaluaton. Each element of VM s perodcally adusted so that half of all functon evaluatons n that drecton are accepted. A fall n T s mposed upon the system wth the Tvarable by T + = T T, where s the teraton counter and Ts the temperature reducton factor. Thus, as T declnes, downhll moves are less lkely to be accepted and the percentage of reectons rses. Gven the scheme for the selecton for VM, VM falls. Thus, as T declnes, VM falls and SA focuses upon the most promsng area for optmzaton. A full descrpton of ths algorthm s found n Corana et al., 25 and the FOTAN subroutne mplemented by Goffe et al. 27 s used n ths work. The choce of the coolng schedule s a crucal aspect n the mplementaton of SA because t affects the numercal performance of the optmzaton procedure. Based on prelmnary calculatons, we propose the followng values for the quanttes of SA: T 0 = 0.0, T=0.85andN T =5,whereN T s the teraton number before temperature reducton. esults and dscusson For the present study, four eamples are used to compare the numercal performance of the three equaton-solvng methods. In all eamples, we have consdered a tolerance of 0 6 for the convergence of all methods. All eamples are solved 25 tmes to evaluate the relablty and effcency of solvng methods, usng n each tme dfferent random ntal values for unknowns, T or p. All calculatons were performed on a Processor Intel Centrno 600 MHz wth.00 GB of AM. The performance of the three equaton-solvng methods s tested usng the followng crterons: a) success rate of fndng the bubble or dew pont n the performed calculatons, b) average number of K-value evaluatons and c) computaton tme. Phase stablty of all calculated bubble and dew ponts s evaluated by mnmzng the tangent plane dstance functon for reactve mtures 2 nn TPDF ( {} { F }) (2) The TPDF functon s globally mnmzed usng the Smulated Annealng method. Bonlla-Petrcolet 24 has tested the SA method wth several reactve systems and found that ths method s robust for reactve phase stablty analyss. All calculated bubble and dew ponts are stable. Eample. sobutene + methanol methyl ter-butyl ether Our frst eample s the equlbrum reacton of sobutene () and methanol (2) to produce methyl ter-butyl ether (), A + A 2 A. Ths system has been analyzed by Okasnsk and Doherty, 28 Maer et al. 29 and Hardng and Floudas 0 n the contet of calculaton of reactve homogeneous azeotropes. We have consdered non-dealty n the lqud phase and deal gas n the vapor phase. We use the Wlson equaton for the actvty coeffcent n the lqud phase and the Antone equaton to calculate the propertes of vapor phase. Maer et al. 29 report the parameters for both equatons. The reacton equlbrum constant s assumed to be ndependent of temperature K eq = 0.04, and methyl ter-butyl ether (MTBE) s the erence component so that the transformed mole fractons are and 2 2. We use and temperature as unknowns for all equaton-solvng methods. The ntal values are randomly generated n the feasble domans of (0, ) for and (0.0, 00) C for temperature. These domans are also used n the mplementaton of the OSA method. Bubble and dew pont calculatons are performed at dfferent pressures for a feed of F (0.5, 0.5) and the results are reported n Tab. The numercal performance of tested methods appears n Tab 2. In all calculatons,
6 A. BONILLA-PETICIOLET et al., Bubble and Dew Pont Calculatons n, Chem. Bochem. Eng. Q. 20 (2) 8 (2006) 5 Table Bubble and dew pont calculatons for sobutene + methanol methyl ter-butyl ether. Wlson equaton and deal gas, K eq = 0.04 and F (0.5,0.5). p /atm Dew pont Bubble pont T / C T / C (0.096, ).944 ( , ) (0.0274, ) ( , ) (0.05, ) 5.56 (0.9224, ) (0.0882, 0.968) ( , ) Table 2 Performance of equaton-solvng methods tested n bubble and dew pont calculatons n reactve systems Success rate n percent (No. of K-values evaluatons) Eample No. T or p Bubble pont Dew pont SESN DESS OSA SESN DESS OSA 2atm 6 (7) 00 (56) 00 (2089) 6 (0) 00 (42) 00 (245) 4atm 44 (72) 00 (52) 00 (2409) 6 (54) 00 (59) 00 (229) 6atm 44 (84) 00 (48) 00 (2265) 20 () 00 (02) 00 (245) 8atm 60 (84) 00 (4) 00 (2289) 24 (44) 00 (06) 00 (227) 50 C 68 (55) 00 (7) 96 (42977) 4 (22) 00 (79) 00 (4297) 2 70 C 76 (5) 00 (247) 00 (445) 20 (2) 00 (275) 00 (477) 90 C 68 (59) 00 (208) 00 (4297) 40 (56) 00 (220) 00 (4505) atm 76 (86) 00 (270) 00 (24) 4 (22) 00 (59) 00 (2677) 5atm 28 (50) 00 (29) 00 (86) 2 (29) 00 (420) 00 (22) 0atm 64 (64) 00 (70) 00 (2005) 8 (22) 00 (00) 00 (2449) atm 28 (20) 88 (92) 00 (27) 28 (2) 00 (49) 00 (227) 4 4atm 6 (6) 80 (46) 00 (2257) 28 (24) 00 () 00 (224) 6atm 28 (5) 56 (40) 00 (205) 56 (2) 00 (27) 00 (247) 8atm 28 () 48 (5) 00 (269) 72 (24) 00 (5) 00 (24) Both parameters are calculated based on 25 calculatons wth random ntal values. DESS and OSA methods fnd the bubble and dew ponts wthout problems whle SESN dverged several tmes. We observe that SESN shows more falures n dew pont calculatons. These results ndcated, that even for smple reactve systems, bubble and dew pont calculatons are challengng for classcal equaton-solvng methods and requre robust strateges. The number of K-value evaluatons of solvng-methods s gven by: OSA >>> DESS > SESN. Computaton tme mantans around of 6 s for OSA, 0.0 s for DESS and 0.02 s for SESN, respectvely. Eample 2. sobutene + methanol methyl ter-butyl ether wth n-butane as nert In the second eample, we use the same reacton as before but nclude n-butane (4) as an nert. Ung and Doherty, 9 Hardng and Floudas 0 and Maer et al. 29 have also studed ths reactve equlbrum problem. The transformed mole fractons for ths mture are 2. 2 and
7 6 A. BONILLA-PETICIOLET et al., Bubble and Dew Pont Calculatons n, Chem. Bochem. Eng. Q. 20 (2) 8 (2006) Table Bubble and dew pont calculatons for sobutene + methanol methyl ter-butyl ether wth n-butane as nert. Wlson equaton and deal gas. z (0.,0.,0.4) T / C Dew pont Bubble pont p /atm p / atm (0.400, , 0.4).9 (0.2692, , ) (0.4602, , 0.606) 5.7 (0.4796, , ) (0.4077, , 0.676) 9.02 (0.6826, , ) 4 4 2, where the MTBE s the erence component. Agan, the Wlson and Antone equatons are used for the calculaton of thermodynamc propertes. Maer et al. 29 report the parameters for both equatons. The reacton equlbrum constant s determned from G rs / = T T ln T, where T s n K. For SESN and DESS methods, we use, 2 and T as unknowns whle n the optmzaton approach, the obectve functon s mnmzed wth respect to, 2, 4 and T. We set a feasble doman of (0, ) for all the transformed mole fractons and (0.000, 20) bar for pressure. The ntals values are randomly generated n these domans. For OSA method, at each evaluaton of the obectve functon, the transformed mole fractons are normalzed to unty. We have analyzed an ntal global composton of F (0., 0., 0.4). Tab shows the results of bubble and dew pont calculatons at dfferent temperatures and Tab 2 shows the performance of the three equaton-solvng methods. For ths case, DESS shows 00 % relablty for the locaton of bubble and dew ponts at all tested temperatures, whle OSA method fals only one tme n all calculatons performed. Agan, SESN shows several falures n the locaton of reactve saturaton condtons. For dew pont calculatons, the performance of SESN method s very poor. OSA method shows the hghest computaton tme (22 s) whle SESN and DESS showed a computaton tme of 0.08 and 0. s, respectvely. Eample. acetc acd + sopropanol sopropyl acetate + water Thrd eample s the formaton of sopropyl acetate through esterfcaton reacton of sopropanol and acetc acd: acetc acd () + sopropanol (2) sopropyl acetate () + water (4). Hardng and Floudas 0 and Maer et al. 29 have studed ths system. For the lqud phase, we use actvty coeffcents from the NTL equaton and consder the assocaton of acetc acd n the vapor phase, usng a correcton factor as suggested by Maer et al. 29 The correcton factor for the assocaton s gven by and sat 2 / ( 4kp ) 2 / [ 4kp( 2)] 2 / 2 { [ 42 kp( )] } 2 / ( 2 ){ [ 4kP ( 2 )] } (22) for = 2,, 4(2) sat where p s the saturaton pressure of pure component and the dmerzaton constant s calculated from log 0 k = log 0 k = /T wth T n Kelvn and k and k n Pa. Maer et al. 29 gves the Antone and NTL equaton parameters. The reacton equlbrum constant s ndependent of temperature, K eq = 8.7, and water s the erence component. The transformed mole fractons are gven by 2 4, 2 and We use, 2 and T as 4 unknowns for all methods. The ntal values are randomly generated usng the feasble doman (0, ) for transformed mole fractons and (50, 200) C for temperature. Bubble and dew pont calculatons are performed for a feed of F (0.5, 0., 0.2) at dfferent pressures (see Tab 4and Tab 2) shows the numercal performance of all tested methods. DESS and OSA methods are very relable to fnd the bubble and dew pont condtons and they dd not show any falures n all the calculatons performed. Agan, SESN showed the worst performance and ts numercal behavor s hghly dependent on ntal guesses. The mean number of K-value evaluatons s gven by: SESN < DESS << OSA, whle the computaton tme mantans around of 7 s for OSA, 0.2 s for DESS and 0.0 s for SESN.
8 A. BONILLA-PETICIOLET et al., Bubble and Dew Pont Calculatons n, Chem. Bochem. Eng. Q. 20 (2) 8 (2006) 7 Table 4 Bubble and dew pont calculatons for acetc acd + sopropanol sopropyl acetate + water. NTL equaton and deal gas,k eq = 8.7 and F (0.5,0.,0.2). p / atm Dew pont Bubble pont T / C T / C 94.5 ( , 0.786, 0.842) (0.79, , ) (0.6046, 0.496, ) ( , , ) ( , 0.8, ) (0.442, , ) Eample 4. A +A 2 A and 2A A 4 +A 2 Our last eample s a hypothetcal quaternary system that follows the reactons: A + A 2 A and 2A A 4 + A 2. Ung and Doherty 9 have studed the phase equlbrum behavor of ths mture. We have consdered deal behavor for both lqud and vapor phases. The Antone equaton s used to calculate the vapor pressures of pure components wth the parameters reported by Ung and Doherty. 9 The transformed mole fractons are defned usng A and A 4 as erence components, so they are gven by and The chemcal equlbrum constants are calculated usng K eq ep and 768 T K eq ep , where T T 2 T s gven n Kelvn. We use and temperature as unknowns and ther feasble domans are (0, ) and (27.5, 400) K, respectvely. andom ntal values are generated nsde these domans. Bubble and dew pont calculatons are performed for a feed F (0.5, 0.5) at several pressures. The results of phase equlbrum calculatons appear n Tab 5 and performance of all tested methods are reported n Tab 2. Only, the OSA method showed 00 % relablty to locate the bubble and dew ponts n ths reactve system, however, t showed the hghest computaton tme (. seconds) and the hghest mean number of K-values evaluatons. On the other hand, several falures are reported for the DESS and SESN methods. For non-reactve systems, Mchelsen has ndcated that bubble and dew pont calculatons usng hybrd models rarely present convergence problems. However, our results show that ths doesn t apply for phase equlbrum calculatons n reactve systems. Even for smple reactng systems, SESN and DESS methods present convergence problems and depend on the ntal guesses. The OSA method s very relable to fnd the bubble and dew ponts and ts performance s almost ndependent of ntal values. However, we can not offer a theoretcal guaranty for the locaton of bubble and dew ponts usng ths optmzaton approach. Fnally, we have tested the performance of these equaton-solvng methods wth other reactve systems reported by Okasnsk and Doherty, 28 Hardng and Floudas 0 and Maer et al.; 29 and the results, not reported n ths paper, ndcated that OSA and DESS methods are relable strateges for the calculaton of dew and bubble ponts n multreactve mtures. Conclusons New condtons for bubble and dew pont calculatons n multreactve mtures have been derved n terms of the varable transformaton theory of Ung and Doherty. We have tested and compared the numercal performance of three equaton-solvng meth- Table 5 Bubble and dew pont calculatons for a hypothetcal quaternary mture A +A 2 A and 2A A 4 +A 2. Ideal soluton and deal gas. F (0.5,0.5) p /atm Dew pont Bubble pont T /K T /K ( , ) ( , ) (0.0452, ) 28.7 ( , ) (0.0599, ) 7.9 (0.979, ) (0.06, ) ( , 0.027)
9 8 A. BONILLA-PETICIOLET et al., Bubble and Dew Pont Calculatons n, Chem. Bochem. Eng. Q. 20 (2) 8 (2006) ods n the calculaton of bubble and dew ponts usng the new equlbrum condtons. Our results ndcate that an optmzaton approach usng Smulated Annealng s very relable for ths knd of phase equlbrum calculatons n multreactve systems. However, t shows the hghest computaton tme. In the other hand, the equaton-decouplng method wth successve substtuton also presents a good numercal performance and s effcent. It s clear that the smultaneous equaton solvng approach, usng Newton method, s not sutable for bubble and dew pont calculatons n reactve systems because t can dverge very frequently and depends sgnfcantly on ntal values. Based on relablty, we consder that the optmzaton approach usng Smulated Annealng s the most sutable method and perable choce for these equlbrum calculatons n reactve systems, due to ts acceptable convergence speed, and that ntal guesses are not requred. We wll etend ths procedure for performng these equlbrum calculatons usng equatons of state. ACKNOWLEDGEMENTS Authors acknowledge the fnancal support from CONACYT,Insttuto Tecnológco de Aguascalentes and POMEP. Notaton a actvty, mol dm n number of components K eq chemcal equlbrum constant K transformed phase equlbrum constant N teraton number p pressure gven n atm, atm = 0 25 Pa n reacton number unversal gas constant, J mol K T temperature gven n C or K r v stochometrc number of component transformed mole fracton, mole fracton, F transformed feed composton, F feed composton, transformed phase mole fracton, VM vector length Greek letters, phase at equlbrum chemcal potental, J mol actvty coeffcent, fugacty coeffcent, Inde F feed composton erence component phase eferences. Pérez-Csneros,E. S.,Gan,.,Mchelsen,M. L.,Chem. Eng. Sc. 52 (997) Seder,W. D.,Wdagdo,S.,Flud Phase Equlb 2 (996) 28.. Castllo,J.,Grossmann,I. E.,Comp. Chem. Eng. 5 (98) Caster,M.,asmussen,P.,Fredenslund,A.,Chem. Eng. Sc. 44 (989) Xao,W.,Zhu,K.,Yuan,W.,Chen,H. H.,An algorthm for smultaneous chemcal and phase equlbrum calculaton. AIChE J. 5 (989) Grener,H.,An effcent mplementaton of newton s method for comple nondeal chemcal equlbra. Comp. Chem. Eng. 5 (99) Smth,J. V.,Mssen,. W.,Smth,W..,General optmalty crtera for multphase multreacton chemcal equlbrum. AIChE J. 9 (99) Ung,S.,Doherty,M. F.,AIChE J. 4 (995) Ung,S.,Doherty,M. F.,Chem. Eng. Sc. 50 (995) McDonald,C. M.,Floudas,C. A.,Glopeq: a new computatonal tool for the phase and chemcal equlbrum problem. Comp. Chem. Eng. 2 (997).. Stateva,. P.,Wakeham,W. A.,Ind. Eng. Chem. es. 6 (997) Wasylkewcz,S. K.,Ung,S.,Flud Phase Equlb 75 (2000) 25.. Ung,S.,Doherty M. F.,Chem. Eng. Sc. 50 (995) Eckert,E.,Kubcek,M.,Comp. Chem. Eng. 9 (995) Walas,S. M.,Phase equlbra n chemcal engneerng. Butterworth Publshers, USA, Teh,Y. S.,angaah,G. P.,Trans IChemE 80 (2002) Henderson,N.,Fretas,L.,Platt,G. M.,AIChE J. 50 (2004) Krkpatrck,S.,Gelatt,C.,Vecc,M.,Scence 220 (98) Pan,H.,Froozabad,A.,Comple multphase equlbra calculatons by drect mnmzaton of Gbbs free energy by use of smulated annealng. SPE eservor Eva. Eng. Feb (998) Zhu,Y.,Xu,Z.,Flud Phase Equlb 54 (999) Henderson,N.,Olvera,J..,Ind. Eng. Chem. es. 40 (200) Gomes,J.,Henderson,N.,ocha,M.,Macromol. Theory Smul. 0 (200) angaah,g. P.,Flud Phase Equlb (200) Bonlla-Petrcolet,A.,Calculaton of phase equlbrum n reactve and nonreactve systems. Ph. D. Dssertaton, Insttuto Tecnológco de Celaya, Guanauato, Meco, Corana,A.,Marches,M.,Martn,C.,della,S.,ACM T. Math. Software (987) Metropols,N.,osenbluth,A.,osenbluth,M.,Teller,A., Teller E., J. Chem. Phys. 2 (95) Goffe,B.,Ferrer,G.,ogers,J.,Global optmzaton of statstcal functons wth smulated annealng. J. Econometrcs 60 (994) Okasnsk,M. J.,Doherty,M. F.,AIChE J. 4 (997) Maer. W.,Brennecke J. F.,Stadtherr M. A.,Comp. Chem. Eng. 24 (2000) Hardng,S. T.,Floudas,C. A.,Ind. Eng. Chem. es. 9 (2000) Mchelsen,M. L.,Comp. Chem. Eng. 7 (99) 4.
Prediction of steady state input multiplicities for the reactive flash separation using reactioninvariant composition variables
Insttuto Tecnologco de Aguascalentes From the SelectedWorks of Adran Bonlla-Petrcolet 2 Predcton of steady state nput multplctes for the reactve flash separaton usng reactonnvarant composton varables Jose
More informationA SHORT METHOD TO CALCULATE REACTIVE RESIDUE CURVE MAPS
Dstllaton Absorpton 00 A.B. de Haan, H. Koojman and A. Górak (Edtors) All rghts reserved by authors as per DA00 copyrght notce A SHORT METHOD TO CALCULATE REACTIVE RESIDUE CURVE MAPS M. Carrera-Rodríguez,
More informationCalculation of homogeneous azeotropes in reactive and nonreactive mixtures using a stochastic optimization approach
Insttuto Tecnologco de Aguascalentes From the SelectedWorks of Adran Bonlla-Petrcolet 2009 Calculaton of homogeneous azeotropes n reactve and nonreactve mxtures usng a stochastc optmzaton approach Adran
More informationGibbs Energy Minimization Using Simulated Annealing for Two-phase Equilibrium Calculations in Reactive Systems
A. BONILLA-PETRICIOLET et al., Gbbs Enery Mnmzaton Usn Smulated, Chem. Bochem. En. Q. (3) 85 98 (008) 85 Gbbs Enery Mnmzaton Usn Smulated Annealn for Two-phase Equlbrum Calculatons n Reactve Systems A.
More informationIntroduction 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 information3. Be able to derive the chemical equilibrium constants from statistical mechanics.
Lecture #17 1 Lecture 17 Objectves: 1. Notaton of chemcal reactons 2. General equlbrum 3. Be able to derve the chemcal equlbrum constants from statstcal mechancs. 4. Identfy how nondeal behavor can be
More informationPI-21 Data correlation in homogeneous azeotropic mixtures using NRTL model and stochastic optimization methods. A. Briones-Ramírez
ESAT 009 PI- Data correlaton n homogeneous azeotropc mxtures usng NRTL model and stochastc optmzaton methods C. Gutérrez-Antono (, A. Bonlla-Petrcolet (, J. G. Segova-Hernández (3, (,( 4 A. Brones-Ramírez
More information10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization
10.34 Numercal Methods Appled to Chemcal Engneerng Fall 2015 Homework #3: Systems of Nonlnear Equatons and Optmzaton Problem 1 (30 ponts). A (homogeneous) azeotrope s a composton of a multcomponent mxture
More informationAnalysis and Prediction of Input Multiplicity for the Reactive Flash Separation using Reaction- Invariant Composition Variables
Insttuto Tecnologco de Aguascalentes From the SelectedWors of Adran Bonlla-Petrcolet 2012 Analyss and Predcton of Input Multplcty for the Reactve Flash Separaton usng Reacton- Invarant Composton Varables
More informationINTRODUCTION TO CHEMICAL PROCESS SIMULATORS
INTRODUCTION TO CHEMICAL PROCESS SIMULATORS DWSIM Chemcal Process Smulator A. Carrero, N. Qurante, J. Javaloyes October 2016 Introducton to Chemcal Process Smulators Contents Monday, October 3 rd 2016
More informationDetermination of Structure and Formation Conditions of Gas Hydrate by Using TPD Method and Flash Calculations
nd atonal Iranan Conference on Gas Hydrate (ICGH) Semnan Unersty Determnaton of Structure and Formaton Condtons of Gas Hydrate by Usng TPD Method and Flash Calculatons H. Behat Rad, F. Varamnan* Department
More informationOpen 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 informationEstimation of the composition of the liquid and vapor streams exiting a flash unit with a supercritical component
Department of Energ oltecnco d Mlano Va Lambruschn - 05 MILANO Eercses of Fundamentals of Chemcal rocesses rof. Ganpero Gropp Eercse 8 Estmaton of the composton of the lqud and vapor streams etng a unt
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationCHAPTER 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(1) The saturation vapor pressure as a function of temperature, often given by the Antoine equation:
CE304, Sprng 2004 Lecture 22 Lecture 22: Topcs n Phase Equlbra, part : For the remander of the course, we wll return to the subject of vapor/lqud equlbrum and ntroduce other phase equlbrum calculatons
More informationGeneral Thermodynamics for Process Simulation. Dr. Jungho Cho, Professor Department of Chemical Engineering Dong Yang University
General Thermodynamcs for Process Smulaton Dr. Jungho Cho, Professor Department of Chemcal Engneerng Dong Yang Unversty Four Crtera for Equlbra μ = μ v Stuaton α T = T β α β P = P l μ = μ l1 l 2 Thermal
More informationy i x P vap 10 A T SOLUTION TO HOMEWORK #7 #Problem
SOLUTION TO HOMEWORK #7 #roblem 1 10.1-1 a. In order to solve ths problem, we need to know what happens at the bubble pont; at ths pont, the frst bubble s formed, so we can assume that all of the number
More informationCinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure
nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n
More informationSupplementary 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 informationEnergy, Entropy, and Availability Balances Phase Equilibria. Nonideal Thermodynamic Property Models. Selecting an Appropriate Model
Lecture 4. Thermodynamcs [Ch. 2] Energy, Entropy, and Avalablty Balances Phase Equlbra - Fugactes and actvty coeffcents -K-values Nondeal Thermodynamc Property Models - P-v-T equaton-of-state models -
More informationChapter 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 informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationis 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 informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More information( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.
Chemstry 360 Dr. Jean M. Standard Problem Set 9 Solutons. The followng chemcal reacton converts sulfur doxde to sulfur troxde. SO ( g) + O ( g) SO 3 ( l). (a.) Wrte the expresson for K eq for ths reacton.
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationExercises of Fundamentals of Chemical Processes
Department of Energ Poltecnco d Mlano a Lambruschn 4 2056 MILANO Exercses of undamentals of Chemcal Processes Prof. Ganpero Gropp Exercse 7 ) Estmaton of the composton of the streams at the ext of an sothermal
More information36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to
ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationIf two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9.
9.9 Real Solutons Exhbt Devatons from Raoult s Law If two volatle and mscble lquds are combned to form a soluton, Raoult s law s not obeyed. Use the expermental data n Table 9.3: Physcal Chemstry 00 Pearson
More informationA Self-Consistent Gibbs Excess Mixing Rule for Cubic Equations of State: derivation and fugacity coefficients
A Self-Consstent Gbbs Excess Mxng Rule for Cubc Equatons of State: dervaton and fugacty coeffcents Paula B. Staudt, Rafael de P. Soares Departamento de Engenhara Químca, Escola de Engenhara, Unversdade
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationAdiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram
Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationGlobal 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 informationGrover s Algorithm + Quantum Zeno Effect + Vaidman
Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the
More informationMODELING THE HIGH-PRESSURE BEHAVIOR OF BINARY MIXTURES OF CARBON DIOXIDE+ALKANOLS USING AN EXCESS FREE ENERGY MIXING RULE
Brazlan Journal of Chemcal Engneerng ISSN 0104-6632 Prnted n Brazl Vol. 21, No. 04, pp. 659-666, October - December 04 MODELING THE HIGH-PRESSURE BEHAVIOR OF BINARY MIXTURES OF CARBON DIOXIDE+ALKANOLS
More informationName: SID: Discussion Session:
Name: SID: Dscusson Sesson: Chemcal Engneerng Thermodynamcs 141 -- Fall 007 Thursday, November 15, 007 Mdterm II SOLUTIONS - 70 mnutes 110 Ponts Total Closed Book and Notes (0 ponts) 1. Evaluate whether
More informationMean Field / Variational Approximations
Mean Feld / Varatonal Appromatons resented by Jose Nuñez 0/24/05 Outlne Introducton Mean Feld Appromaton Structured Mean Feld Weghted Mean Feld Varatonal Methods Introducton roblem: We have dstrbuton but
More informationONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00
ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationMarkov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement
Markov Chan Monte Carlo MCMC, Gbbs Samplng, Metropols Algorthms, and Smulated Annealng 2001 Bonformatcs Course Supplement SNU Bontellgence Lab http://bsnuackr/ Outlne! Markov Chan Monte Carlo MCMC! Metropols-Hastngs
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationThe Minimum Universal Cost Flow in an Infeasible Flow Network
Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran
More informationFUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM
Internatonal Conference on Ceramcs, Bkaner, Inda Internatonal Journal of Modern Physcs: Conference Seres Vol. 22 (2013) 757 761 World Scentfc Publshng Company DOI: 10.1142/S2010194513010982 FUZZY GOAL
More informationAn identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites
IOP Conference Seres: Materals Scence and Engneerng PAPER OPE ACCESS An dentfcaton algorthm of model knetc parameters of the nterfacal layer growth n fber compostes o cte ths artcle: V Zubov et al 216
More informationChapter 11: Simple Linear Regression and Correlation
Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests
More informationA MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS
Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationarxiv:quant-ph/ Jul 2002
Lnear optcs mplementaton of general two-photon proectve measurement Andrze Grudka* and Anton Wóck** Faculty of Physcs, Adam Mckewcz Unversty, arxv:quant-ph/ 9 Jul PXOWRZVNDR]QDRODQG Abstract We wll present
More informationNon-Ideality Through Fugacity and Activity
Non-Idealty Through Fugacty and Actvty S. Patel Deartment of Chemstry and Bochemstry, Unversty of Delaware, Newark, Delaware 19716, USA Corresondng author. E-mal: saatel@udel.edu 1 I. FUGACITY In ths dscusson,
More informationReliable Phase Stability Analysis for Asymmetric Models
Relable Phase Stablty Analyss for Asymmetrc Models Gang Xu 1, Wllam D. Haynes 2* and Mark A. Stadtherr 2 1 Invensys/SmSc-Esscor, 26561 Rancho Parkway South, Sute 1, Lake Forest, CA 9263 USA 2 Department
More informationVapor-Liquid Equilibria for Water+Hydrochloric Acid+Magnesium Chloride and Water+Hydrochloric Acid+Calcium Chloride Systems at Atmospheric Pressure
Chnese J. Chem. Eng., 4() 76 80 (006) RESEARCH OES Vapor-Lqud Equlbra for Water+Hydrochlorc Acd+Magnesum Chlorde and Water+Hydrochlorc Acd+Calcum Chlorde Systems at Atmospherc Pressure ZHAG Yng( 张颖 ) and
More informationI 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 informationThe 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 informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationThe ChemSep Book. Harry A. Kooijman Consultant. Ross Taylor Clarkson University, Potsdam, New York University of Twente, Enschede, The Netherlands
The ChemSep Book Harry A. Koojman Consultant Ross Taylor Clarkson Unversty, Potsdam, New York Unversty of Twente, Enschede, The Netherlands Lbr Books on Demand www.bod.de Copyrght c 2000 by H.A. Koojman
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationUniversity of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014
Lecture 16 8/4/14 Unversty o Washngton Department o Chemstry Chemstry 452/456 Summer Quarter 214. Real Vapors and Fugacty Henry s Law accounts or the propertes o extremely dlute soluton. s shown n Fgure
More informationELASTIC 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 informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationLecture 16 Statistical Analysis in Biomaterials Research (Part II)
3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan
More informationLecture 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 informationDesign and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm
Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationMarkov Chain Monte Carlo Lecture 6
where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways
More information: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11:
764: Numercal Analyss Topc : Soluton o Nonlnear Equatons Lectures 5-: UIN Malang Read Chapters 5 and 6 o the tetbook 764_Topc Lecture 5 Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton
More informationNAME and Section No. it is found that 0.6 mol of O
NAME and Secton No. Chemstry 391 Fall 7 Exam III KEY 1. (3 Ponts) ***Do 5 out of 6***(If 6 are done only the frst 5 wll be graded)*** a). In the reacton 3O O3 t s found that.6 mol of O are consumed. Fnd
More informationDesign of Reactive Distillation with Thermal Coupling for the Synthesis of Biodiesel using Genetic Algorithms
19 th European Symposum on Computer Aded Process Engneerng ESCAPE19 J. Jeowsk and J. Thulle (Edtors) 2009 Elsever B.V./Ltd. All rghts reserved. Desgn of Reactve Dstllaton wth Thermal Couplng for the Synthess
More informationConstitutive Modelling of Superplastic AA-5083
TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy
More informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationComputation of Phase Equilibrium and Phase Envelopes
Downloaded from orbt.dtu.dk on: Sep 24, 2018 Computaton of Phase Equlbrum and Phase Envelopes Rtschel, Tobas Kasper Skovborg; Jørgensen, John Bagterp Publcaton date: 2017 Document Verson Publsher's PDF,
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationCurve 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 informationSupporting Information
Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to
More informationOPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming
OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or
More informationON THE MIN-MAX FORMULATION OF MULTIPHASE CHEMICAL EQUILIBRIUM PROBLEM AND STABILITY ANALYSIS
63 PHYSIAL HEMISTRY ON THE MIN-MA FORMULATION OF MULTIPHASE HEMIAL EQUILIBRIUM PROBLEM AND STABILITY ANALYSIS DAN GEANĂ Dept. Appled Physcal hemstry, Unversty "Poltehnca" Bucharest, Spl. Independente 33,
More informationSection 3.6 Complex Zeros
04 Chapter Secton 6 Comple Zeros When fndng the zeros of polynomals, at some pont you're faced wth the problem Whle there are clearly no real numbers that are solutons to ths equaton, leavng thngs there
More informationVQ widely used in coding speech, image, and video
at Scalar quantzers are specal cases of vector quantzers (VQ): they are constraned to look at one sample at a tme (memoryless) VQ does not have such constrant better RD perfomance expected Source codng
More informationModule 3: The Whole-Process Perspective for Thermochemical Hydrogen
"Thermodynamc Analyss of Processes for Hydrogen Generaton by Decomposton of Water" by John P. O'Connell Department of Chemcal Engneerng Unversty of Vrgna Charlottesvlle, VA 2294-4741 A Set of Energy Educaton
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationProblem 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 informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationCOS 521: Advanced Algorithms Game Theory and Linear Programming
COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton
More informationCIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M
CIS56: achne Learnng Lecture 3 (Sept 6, 003) Preparaton help: Xaoyng Huang Lnear Regresson Lnear regresson can be represented by a functonal form: f(; θ) = θ 0 0 +θ + + θ = θ = 0 ote: 0 s a dummy attrbute
More information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationAssignment 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 informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More information