Mass and Heat Diffusion in Ternary Polymer Solutions: A Classical Irreversible Thermodynamics Approach

Size: px
Start display at page:

Download "Mass and Heat Diffusion in Ternary Polymer Solutions: A Classical Irreversible Thermodynamics Approach"

Transcription

1 Mass and Heat ffuson n Ternary Polymer Solutons: A Classcal Irreversble Thermodynamcs Approach S. Shams Es-hagh and M. Cakmak * epartment of Polymer Engneerng, 5 S. Forge St., The Unversty of Akron, Akron, Oho, , USA ABSTRACT Governng equatons for evoluton of concentraton and temperature n three-component systems were derved n the framework of classcal rreversble thermodynamcs usng Onsager s varatonal prncple and were presented for solvent/solvent/polymer and solvent/polymer/polymer systems. The dervaton was developed from the Gbbs equaton of equlbrum thermodynamcs usng the local equlbrum hypothess, Onsager recprocal relatons and Prgogne s theorem for systems n mechancal equlbrum. It was shown that the detals of mass and heat dffuson phenomena n a ternary system are completely expressed by a 3 3 matrx whose entres are mass dffuson coeffcents (4 entres), thermal dffuson coeffcents ( entres) and three entres that descrbe the evoluton of heat n the system. The entres of the dffuson matrx are related to the elements of Onsager matrx that are bounded by some constrants to satsfy the postve defnteness of entropy producton n the system. All the elements of dffuson matrx were expressed n terms of dervatves of exchange chemcal potentals of the components wth respect to concentraton and temperature. The spnodal curves of ternary polymer solutons were derved from the governng equatons and ther correctness was checked by the Hessan of free energy densty. Moreover, t was proved that settng crossdffuson coeffcents to zero results n a contradcton, and the governng equatons wthout cross-dffuson coeffcents do not express the actual phase behavor of the system. *Correspondng author: cakmak@uakron.edu

2 I. INTROUCTION Mass and heat dffuson n polymer solutons s of great mportance n the technologes related to pantng, coatng, nkjet prntng, plastc flms and producton of electronc devces []. ffuson processes also play a key role n self-assembly [] and mportant areas of study n soft matter physcs [3]. A proper understandng of dffuson phenomena n polymer solutons s requred n order to have a control on the processes nvolvng polymer solutons. Polymer solutons commonly used for ndustral purposes are generally mult-component systems,.e., ternary polymer solutons of one polymer n two solvents or two polymers n one solvent [4]. Most of the expermental and theoretcal studes thus far, have been about dffuson n bnary systems, and despte the broad range of applcatons nvolve ternary polymer solutons [4-4] there s no theoretcal work that formulates thermodynamcally consstent governng equatons for mass and heat dffuson n ternary polymer solutons n a way that all the mass and thermal dffuson coeffcents necessary to descrbe the transport phenomena are derved n terms of thermodynamc varables. By ncreasng the number of components n a mxture and emergence of cross-dffuson effects, the transport phenomena become very complex. As a result, a large number of mass dffuson coeffcents are requred to descrbe the mass transport []. Ths complexty wll ncrease by addng a spatally varyng temperature feld to the system, snce the effect of thermal dffuson and contrbutons of mass fluxes to the heat flow must also be consdered. Mass and heat dffuson are rreversble phenomena and they contrbute to entropy producton n the system. Hence, the governng equatons for mass and heat evoluton n ternary polymer solutons must be derved n the framework of classcal rreversble thermodynamcs n whch the balance of entropy plays a crucal role. An overvew of the governng equatons formulated for mass and

3 heat dffuson n ternary polymer solutons ndcates that these equatons were derved wthout consderng the balance of entropy n the system. Vrentas et al. [5] derved governng equatons for sothermal mass dffuson n ternary polymer solutons based on a theory due to Bearman [6]. They consdered a ternary polymer system contanng one polymer and two solvents. In that system, the concentraton of solvents was much lower than that of the polymer. They set the cross-dffuson coeffcents to zero and estmated the prncpal mass dffuson coeffcents by the self-dffuson coeffcents of the solvents. Usng the same approach based on the statstcal theory of Bearman [6], Shojae et al. [7] developed a model for non-sothermal mass dffuson n a ternary polymer soluton contanng a polymer and two solvents. The heat transfer n ther model was descrbed by the Fourer s law of heat conducton. Alsoy and uda [8] studed the dryng of ternary polymer solutons usng a model n whch the mass dffuson coeffcents were expressed n terms of self-dffuson coeffcents of the solvents. In that model they assumed temperature to be temporally dependent and spatally ndependent. abral et al. [4] modeled the dryng process of a ternary polymer soluton as an sothermal mass dffuson problem. ue to the lack of relable expermental data and the absence of a sutable predctve theory, they set the cross-dffuson coeffcents to zero. In the above-mentoned formulatons, there s no crteron that guarantees the governng equatons are thermodynamcally consstent. An analyss of these models developed for dffuson n mult-component systems ndcates that some of them are not consstent wth the Onsager recprocal relatons [9]. Moreover, n these studes, contrbuton of temperature gradent to mass dffuson n ternary polymer solutons was neglected. Furthermore, the heat transfer was modeled usng the Fourer s law of heat conducton. Ths formulaton does not take nto account the effect of mass dffuson on the heat transfer and the couplng between mass and heat transfer s only due to the temperature dependency of mass dffuson coeffcents. 3

4 Recently, we developed a thermodynamcally consstent model of mass and heat dffuson n bnary polymer solutons []. In ths paper, we focus on generalzaton of ths model to ternary polymer solutons. Although ths study ams to develop a theoretcal understandng of dffuson n ternary polymer solutons, the governng equatons derved heren can be used for any ternary system. Snce, the governng equatons have been formulated for the general case of non-sothermal three-component systems, the effect of temperature gradent on the evoluton of concentraton of the components was shown by dervaton of thermal dffuson coeffcents n terms of phenomenologcal coeffcents and dervatves of exchange chemcal potentals of the components wth respect to temperature. Moreover, t was proved that crossdffuson coeffcents play a crucal role n expressng the phase behavor of a ternary system and settng the cross-dffuson coeffcents to zero leads to a contradcton. II. THEORY We consder a non-reactng ternary system n whch non-convectve mass dffuson and heat conducton occur. The system s consdered to be n one phase, far from the crtcal regon of phase separaton. We also assume the system s n mechancal equlbrum. The hypothess of local equlbrum allows the fundamental equaton of classcal thermodynamcs to be vald for every volume element n the system, although the whole system s not n equlbrum []. Usng ths hypothess and Prgogne s theorem for systems n mechancal equlbrum, the rate of change of entropy per unt volume for a three-component system such as solvent/solvent/polymer or solvent/polymer/polymer s gven by (see Appendx A) 4

5 ds Jq,3 J,3J dt T T J J J J J q,3,3,3,3 T T T, () where s, J q, J, J and T are entropy per unt volume, heat flux, mass flux of component, mass flux of component and absolute temperature, respectvely.,3 and,3 are exchange chemcal potentals of components and 3, and and 3, respectvely. Equaton () represents the rate of change of entropy per unt volume of the mxture n terms of dvergence of entropy flux and rate of entropy producton per unt volume of the system, whch can be presented n the blnear form of thermodynamcs forces X,3,,3, T T T T and the conjugated fluxes J J, J, Jq,3 J,3J, T J J J J J. () q,3,3,3,3 T T T As one can see due to a constrant for the fluxes, n a ternary system the mass flux of two components are expressed explctly n the rate of entropy producton. In effect, for a general case of an n -component system, n equaton are requred to descrbe the mass dffuson n the system. The governng equatons can be found usng lnear flux-force relatonshps consderng the Onsager recprocal relatons and also they can be derved from Onsager s varatonal prncple [] by maxmzng J, the dfference between rate of entropy producton n the system and a dsspaton functon presented n terms of thermodynamc forces n the system 3 ( X, X, X 3) Lk X X k, (3) k, 5

6 here L k are phenomenologcal coeffcents. Usng lnear flux-force relatonshps, the fluxes can be wrtten n terms of the forces as shown n Eqs. (4) - (6) whch n matrx presentaton would yeld Eq. (7) J l,3 l,3 l 3 T, T T T J l,3 l,3 l 3 T, T T T Jq,3 J,3J l3,3 l3,3 l33 T, T T T,3 J T l l l3 J l l l. 3,3 T Jq,3 J,3J l3 l3 l 33 T T (4) (5) (6) (7) Entres of matrx L ( lj ) 33 are Onsager s coeffcents and based on Onsager s recprocty relatons, matrx L s a symmetrc matrx [3]. In condtons for whch lnear flux-force relatons are vald, rate of entropy producton takes the quadratc form l,3 l,3 l33 T T T T l,3,3 l3,3 T T T T T l3,3 T. T T (8) Matrx L ( lj ) 33 that satsfes Eq. (8) should be postve defnte and to be so, ts entres should satsfy the condtons 6

7 l l l3 l, ll l, det l l l3. l3 l3 l 33 (9) Snce L s a symmetrc postve defnte matrx, all of ts dagonal entres, l are postve. By consderng l l l l T T T T 3 3 (where = or or 3),,,, () and replacng,3 and,3 wth the rght hand sdes of Eqs. () and (), respectvely, knowng the fact that chemcal potental s a functon of concentraton and temperature, T,,3,3,3,3 T, T T,, () T,,3,3,3,3 T, T T,, () we can recast Eqs. (4) - (6) n the forms J,3,3 T, T,,3,3 T, T,,3,3 T, T T T,, (3) 7

8 J,3,3 T, T,,3,3 T, T,,3,3 T, T T T,, (4),3,3 Jq,3 J,3J T, T,,3,3 T, T,,3,3 3 T. T T T,, (5) be satsfed: In order to preserve the postve defnteness of matrx L, the followng condtons must,,3,. 3 3 (6) After dervng the heat and mass fluxes, usng the defnton of the tme dervatve of enthalpy per unt volume (see Eq. A.4) and mass concentraton for component (see Eq. A.5), we obtan the governng equatons from the dvergence of the fluxes 8

9 d,3,3 dt T, T,,3,3 T, T, T T T,,,3,3 T, (7) d,3,3 dt T, T,,3,3 T, T,,3,3 T, T T T,, (8) dt,3,3 cp,3 J,3J dt T, T,,3,3 T, T,,3,3 3 T, T T T,, (9) where and c P are the mass densty and sobarc specfc heat capacty of the system, respectvely. Usng Eqs. (7) and (), one can express the equaton of heat flux n the form (see Appendx B) J k T, () q 9

10 where J q and k are reduced heat flux and thermal conductvty and are gven by Jq J q J J,3,3, () k T 3. () Equaton () shows thermal conductvty as a functon of phenomenologcal coeffcents and the constrant gven by Eq. (6) guarantees that the thermal conductvty s a postve quantty. There s another way of expressng the equaton of heat dffuson. Ths can be done by replacng the mass fluxes n Eq. (9) usng Eqs. (3) and (4). ong so, the equaton of heat dffuson wll be dt,3,3 cp (,3,3 ) (,3,3 ) dt T, T,,3,3 (,3,3 ) (,3,3 ) T, T,,3,3 (,3,3 ) (,3,3 ) T T T,, (,3,3 3) T. T (3) Equatons (7), (8) and (3) can be wrtten n a compact form as 3 d. dt 3 cpt T (4)

11 We call the 3 3 matrx n Eq. (4) dffuson matrx whose entres j (see Appendx C) descrbe the detals of mass and heat dffuson n the system. Mass dffuson matrx s a matrx wth entres,, and. and are cross-dffuson coeffcents. 3 and 3 are thermal dffuson coeffcents. The three entres n the thrd row of the dffuson matrx descrbe the heat conducton n the system. III. RESULTS AN ISCUSSION In secton, the mass and heat dffuson coeffcents were derved as part of the model n terms of dervatves of exchange chemcal potentals of the components wth respect to concentraton and temperature. Thus, n the model developed heren, there s no need for defnton of dffuson coeffcents. The approach adopted n ths paper for dervng the governng equatons n ternary systems can be compared wth the model developed by Curtss and Brd [4] who derved the generalzed Maxwell-Stefan equatons for the mult-component dffuson. In that model the dffuson coeffcents were not derved. Curtss and Brd [4] consdered two dfferent defntons for mass dffuson coeffcents; zero-dagonal and symmetrc dffusvty defntons and based on these defntons they derved generalzed Maxwell-Stephan equatons. The 3 governng equatons derved n secton can be used for any ternary mxture. The entres of the dffuson matrx wll be dfferent dependng on the chemcal potentals of the components n a partcular system. In order to express the governng equatons for ternary polymer solutons, one needs to have the chemcal potentals of the components. We consder two general ternary polymer solutons; A soluton of one polymer n two solvents,

12 solvent()/solvent()/polymer(3) and a soluton of two polymers n one solvent, solvent()/polymer()/polymer(3). The numbers n parenthess wll be used as subscrpts n equatons to represent components n the system. The chemcal potentals of the components n a ternary polymer soluton can be derved by extenson of Flory-Huggns theory to ternary polymer soluton [5,6]. Then, the exchange chemcal potentals of the components can be derved usng Eqs. (A.) and (A.). A. ffuson matrx for solvent()/solvent()/polymer(3) systems In case of solvent()/solvent()/polymer(3) systems, the exchange chemcal potentals of the components wll be (see Appendx ) M(,3,3 ) ln ln RT N3 N3, M (,3,3 ) ln ln RT N3 N3, (5) (6) where,3 and,3 are,3 3, M N 3 (7),3 3. M N 3 (8) Here, M, M, N 3 and j are chemcal potental of component n ts pure lqud state, molecular weght of component, molecular weght of component, degree of polymerzaton

13 of component 3 (polymer) and the Flory-Huggns nteracton parameter between components and j. Usng Eqs. (5) (8) the dffuson matrx for solvent()/solvent()/polymer(3) systems can be derved (see Appendx E). B. ffuson matrx for solvent()/polymer()/polymer(3) systems In case of solvent()/polymer()/polymer(3) systems, the exchange chemcal potentals of the components wll be (see Appendx ) M(,3,3 ) ln ln RT N3 N3, m (,3,3) ln ln RT N N3 N N3, (9) (3) where m and N are the molecular weght of the repeatng unt of the component (polymer ) and ts degree of polymerzaton, respectvely and,3 and,3 are gven by,3 3, M N 3 (3) N,3 3. M N 3 (3) Usng Eqs. (9) (3) the dffuson matrx for solvent()/polymer()/polymer(3) systems can be derved (see Appendx F). C. Mass and heat dffuson n an n-component system The governng equatons derved for ternary mxtures can easly be generalzed to descrbe mass and heat dffuson n an n -component system. 3

14 gven by The governng equatons for mass and heat dffuson n an n -component system are ( n ) n d. dt ( n ) ( n ) ( n) ( n)( n) ( n) n ( n) ct P n n ( n ) nn T (33) The ( n) ( n ) submatrx shown nsde a box n Eq. (33) s the mass dffuson matrx and the components of the vector n the last column of dffuson matrx are thermal dffuson coeffcents. The Onsager matrx assocated wth the dffuson matrx s an n n symmetrc postve defnte matrx whose entres must satsfy l l l3 l l n l, ll l, det l l l3,, det. l l l l l n nn (34) Ths concept of generalzng governng equatons to an n -component system may become mportant when a soluton of several polymers s prepared n one or several solvents. Moreover, knowng the fact that polymers have a dstrbuton of molecular weght, a soluton of one polymer n one solvent s a mult-component system of polymer chans wth dfferent molecular weghts n that solvent [8].. ervaton of spnodal curve from the governng equatons. Important nformaton about phase behavor of a system can be extracted from the governng equatons for mass dffuson n the system. The spnodal hypersurface for an n - components system s gven by [9] 4

15 ( n) det. ( n) ( n)( n) (35) In effect, the spnodal curve for a system s the hypersurface on whch the mass dffuson matrx becomes sngular. Therefore, for a ternary mxture the spnodal curve s gven by,3,3,3,3 T, T, T, T, det.,3,3,3,3 T, T, T, T, (36) ue to the constrant, the determnant of the mass dffuson matrx wll be,3,3,3,3 T, T, T, T,. (37) Equaton (37) gves the spnodal curve for a general ternary system. Usng the equatons of exchange chemcal potentals of the components, the spnodal curve for solvent()/solvent()/polymer(3) and solvent()/polymer()/polymer(3) systems wll be 3 3 N3( 3) N3( 3) 3 3, N3( 3) (38) 3 3 N3( 3) N N3( 3) 3 3, N3( 3) (39) 5

16 respectvely. The results derved heren can be checked wth a well-establshed alternatve way n that the spnodal hypersurface of an n -component system can be derved from the determnant of the Hessan matrx of free energy densty of the system [3,3] f f ( n) det. f f ( n) ( n) (4) Therefore, for a ternary mxture the spnodal curve wll be f f f f. (4) Free energy densty functons of solvent()/solvent()/polymer(3) and solvent()/polymer()/polymer(3) systems are gven by [3] f f ln ln ln, (4) N3 ln ln ln, (43) N N3 respectvely. The free energy densty functons are presented n dmensonless form and accordng to Eq. (4) ths does not affect the result of the problem. It s easy to show that the spnodal curves derved from Eq. (4) usng Eqs. (4) and (43) are dentcal wth the spnodal curves derved from the governng equatons. E. The consequence of settng cross-dffuson coeffcents to zero Accordng to the lterature (see for example Ref. 3), due to the lack of data or n order to smplfy the problem the cross-dffuson coeffcents are set to zero. However, t has been observed that droppng cross-dffuson coeffcents resulted n sgnfcant errors n predctons of 6

17 dryng of ternary polymer solutons [3]. In ths secton we dscuss the consequence of gnorng cross-dffuson coeffcents n a ternary system. By droppng the cross-dffuson coeffcent,..e., the mass dffuson matrx for a ternary system wll be presented ether n the form,3 T,,3 T,, (44) or,3 T,,3 T,. (45) ue to the constrant, the spnodal curve of the system wll be,3,3 T, T,, (46) or,3,3 T, T,, (47) respectvely. It s clear that the spnodal curve demonstrated by the governng equatons after reducng the cross-dffuson coeffcents to zero s not the actual spnodal curve of the system. Moreover, Eqs. (46) and (47) lead to a contradcton n that the exchange chemcal potentals of the components,3 and,3 must be a constant or be only a functon of or so that Eqs. 7

18 (46) and (47) are vald. Therefore, one can see that settng the off-dagonal entres of mass dffuson matrx to zero can result n a serous problem n that not only the actual phase behavor of the system wll be altered, but also a contradcton arses. Therefore, one s not elgble to gnore the cross-dffuson coeffcents n the system. It s worth mentonng that the proof provded heren s for a general case and can be consdered for any three-component system. IV. SUMMARY AN CONCLUSIONS 3 governng equatons for mass and heat dffuson n three-component systems were derved n the framework of classcal rreversble thermodynamcs. The formulaton of the governng equatons was developed from the Gbbs equaton usng Onsager s varaton prncple. It was shown that the detals of mass and heat dffuson n a three-component mxture are descrbed by a 3 3 dffuson matrx whose entres are mass dffuson coeffcents, thermal dffuson coeffcents and three entres that descrbe the evoluton of temperature feld n the system. Mass dffuson coeffcents are the entres of mass dffuson matrx, a submatrx whch remans after elmnaton of the thrd row and thrd column of the dffuson matrx. The entres of the dffuson matrx were expressed n terms of phenomenologcal coeffcents and dervatves of exchange chemcal potentals of the components wth respect to concentraton and temperature. The model was formulated for the general case of non-sothermal mass dffuson n three-component systems and can be used for any three-component system, when the chemcal potentals of the components are known as functons of concentraton and temperature. The entres of the dffuson matrx were derved for partcular cases of ternary polymer solutons, solvent/solvent/polymer and solvent/polymer/polymer systems by explotng the Flory-Huggns theory to derve the chemcal potentals of the components. It was shown that the model can be 8

19 easly generalzed to descrbe mass and heat dffuson n an n -component system. In ths general case dffuson matrx has n entres where ( n ) entres are mass dffuson coeffcents, ( n ) entres are thermal dffuson coeffcents and n entres descrbe the evoluton of temperature feld n the system. The spnodal curves of ternary polymer solutons were derved from the entres of the mass dffuson matrx and ther correctness was checked usng the determnant of the Hessan matrx of free energy densty functons of the systems. Moreover, for general case of ternary systems, t was shown that cross-dffuson coeffcents play a crucal role n expressng the phase behavor of the system. In effect, t was proved that settng the cross-dffuson coeffcents to zero leads to a contradcton. ACKNOWLEGMENT We are grateful for fundng enttled Polymde research and Commercalzaton Grant from The State of Oho Thrd fronter program. 9

20 Appendx A. Rate of change of entropy per unt volume of a ternary system Equaton (A.) states the fundamental equaton per unt volume of a mxture; k de Tds dn, (A.) where T, s, e, and n are absolute temperature, entropy per unt volume, enthalpy per unt volume, chemcal potental of component and moles of component per unt volume, respectvely. One can rewrte Eq. (A.), usng the mass concentraton of the components and take a dervatve wth respect to tme to fnd k ds de dc T, dt dt (A.) M dt where c and M are the mass concentraton and molecular weght of component, respectvely. Tme dervatves n Eq. (A.) are substantal tme dervatves gven by d dt v. t (A.3) Here, v s the mean velocty of the components. The tme dervatves of e and c n Eq. (A.) are related to the dvergence of heat and non- convectve mass fluxes, de J, q dt (A.4) dc J. dt (A.5)

21 Accordng to Prgogne s theorem, for systems n mechancal equlbrum an arbtrary frame of reference can be chosen []. The non-convectve mass flux of the component n ths frame of reference whch moves wth the mean velocty v s gven by J v v, (A.6) where, and v are mass densty, volume fracton and the velocty of component, respectvely and the mean velocty v s expressed as k v v. (A.7) It follows from Eqs. (A.6) and (A.7) that the fluxes are not ndependent and ther dependency s gven by k J. (A.8) Rewrtng Eq. (A.) for a three-component system and replacng de dt and dc dt wth the dvergences of the assocated fluxes gven by Eqs. (A.4) and (A.5), we obtan ds T Jq J J J (A.9) dt M M M The subscrpts, and 3 are attrbuted to component, component and component 3 n the system, respectvely. Equaton (A.9) can be reduced after applyng the constrant for fluxes gven by Eq. (A.8) to the form ds T J, q,3 J,3 J (A.) dt where,3 and,3 are the exchange chemcal potentals of the components and are gven by

22 3 3,3, (A.) M M3 3 3,3, (A.) M M3 respectvely. We can rewrte Eq. (A.) by replacng the rght hand sde wth an equvalent form and dvdng both sdes by T, ds Jq,3 J,3J J,3 J,3. (A.3) dt T T T The frst term on the rght hand sde of Eq. (A.3) can be manpulated to an equvalent form so that ds Jq,3 J,3J dt T T J J J J J q,3,3,3,3 T T T. (A.4)

23 Appendx B. ervaton of reduced heat flux Usng Eq. (), and knowng that l l T, l l = T, and 3 3 l3 l3 T the lnear flux-force relatonshps gven by Eq. (7) can be wrtten as J J (B.) T,3,3 T, (B.) T,3,3 T, 3 Jq,3 J,3J,3,3 T, (B.3) T Usng Eqs. (B.) and (B.),,3 and,3 are gven by,3 J J T, T (B.4),3 J J T. T (B.5) By substtutng for,3 and,3 n Eq. (B.3), the heat flux wll be Jq,3 J,3 J T 3 T. (B.6) Equaton (B.6) can be expressed n terms of the reduced heat flux that s gven by Jq J q J J,3,3, (B.7) Therefore, we obtan J q 3 T. T 3 (B.8)

24 Appendx C. Entres of dffuson matrx,3,3 T, T,,,3,3 T, T,,, T T T,3,3 3,,,3,3 T, T,,,3,3 T, T,,, T T T,3,3 3,, ( ) ( ),,3,3 3,3,3,3,3 T, T, ( ) ( ),,3,3 3,3,3,3,3 T, T, ( ) ( ) ( ). T T T,3,3,3,3 3 33,3,3,3,3,, 4

25 Appendx. Exchange chemcal potentals of components n a ternary polymer soluton Exchange chemcal potentals,3 and,3 can be derved usng Eqs. (A.) and (A.) whch can be wrtten n the forms V V,3 3,,3 3, M V3 M V3 (.) where V / V s the rato of molar volumes of the components and wll be consdered equal to the j rato of ther degrees of polymerzaton. Chemcal potentals of the components n a ternary polymer soluton such as solvent()/solvent()/polymer(3) wth respect to ther pure states are gven by RT RT ln , N3 ln , N3 3 3 ln3 N3 N3 RT N N N, (.) (.3) (.4) where, R, N 3 and j are chemcal potentals of component n ts pure lqud state, gas constant, degree of polymerzaton of component 3 (polymer) and the Flory-Huggns nteracton parameter of components and j. Therefore, the exchange chemcal potentals for the solvent()/solvent()/polymer(3) systems wll be M(,3,3) ln ln RT N3 N3, (.5) 5

26 M (,3,3) ln ln RT N3 N3, (.6) where,3 and,3 are,3 3, M N 3 (.7),3 3. M N 3 (.8) Chemcal potentals of the components n a ternary polymer soluton such as solvent()/polymer()/polymer(3) wth respect to ther pure states are gven by ln 3 RT N N3, N ln N 3 RT N3 N N N, N 3 ln3 N3 RT N N N N (.9) (.) (.) Therefore, the exchange chemcal potentals for the solvent()/polymer()/polymer(3) systems wll be M(,3,3 ) ln ln RT N3 N3, (.) 6

27 m (,3,3) ln ln RT N N3 N N3, (.3) where m and N are the molecular weght of the repeatng unt of the component (polymer ) and ts degree of polymerzaton, respectvely and,3 and,3 are gven by,3 3, M N 3 (.4) N,3 3. M N 3 (.5) 7

28 Appendx E. Entres of dffuson matrx for solvent()/solvent()/polymer(3) systems Snce the chemcal potentals of pure components are only functons of temperature [7], the dervatve of wth respect to the concentraton vanshes. Moreover, the nteracton parameters were treated as constants. In the case where a temperature dependency s consdered for nteracton parameters the terms contanng dervatves wth respect to the temperature must be modfed. RT RT, M N3( 3) M N3( 3) RT RT, M N3( 3) M N3( 3) d,3 (,3,3) d,3 (,3,3), dt T dt T T 3 RT RT, M N3( 3) M N3( 3) RT RT, M N3( 3) M N3( 3) d ( ) d ( ), dt T dt T T,3,3,3,3,3,3 3 RT ( ),3,3 3 3 M N3( 3) RT (,3,3 ) 3 3, M N3( 3) 8

29 RT ( ),3, M N3( 3) RT (,3,3 ) 3, M N3( 3) d,3 (,3,3 ) ( ) dt T 33,3,3 d,3 (,3,3) (,3,3 3) (,3,3 ). dt T T 9

30 Appendx F. Entres of dffuson matrx for solvent()/polymer()/polymer(3) systems Snce the chemcal potentals of pure components are only functons of temperature [7], the dervatve of wth respect to the concentraton vanshes. Moreover, the nteracton parameters were treated as constants. In the case where a temperature dependency s consdered for nteracton parameters the terms contanng dervatves wth respect to the temperature must be modfed. RT RT, M N3( 3) m N3( 3) RT RT, M N3( 3) m N N3( 3) d,3 (,3,3) d,3 (,3,3), dt T dt T T 3 RT RT, M N3( 3) m N3( 3) RT RT, M N3( 3) m N N3( 3) d ( ) d ( ), dt T dt T T,3,3,3,3,3,3 3 RT ( ),3,3 3 3 M N3( 3) RT (,3,3 ) 3 3, m N3( 3) 3

31 RT ( ),3, M N3( 3) RT (,3,3 ) 3, m N N3( 3) d,3 (,3,3 ) ( ) dt T 33,3,3 d,3 (,3,3) (,3,3 3) (,3,3 ). dt T T 3

32 e T s n t M c J q J NOTATION enthalpy per unt volume absolute temperature entropy per unt volume chemcal potental of component chemcal potental of component n ts pure state moles of component per unt volume tme molecular weght of component mass concentraton of component heat flux non-convectve mass flux of component mass densty of component volume fracton of component v velocty of component v mean velocty exchange chemcal potental of components and j, j rate of entropy producton per unt volume dsspaton functon L Onsager matrx X vector of thermodynamc forces J vector of thermodynamc fluxes,,, phenomenologcal coeffcents mass densty of the system sobarc specfc heat capacty of the system c P J q k j j R f N reduced heat flux thermal conductvty entres of dffuson matrx Flory-Huggns nteracton parameter of components and j gas constant free energy densty functon degree of polymerzaton of component m molecular weght of repeatng unt of component V molar volume of component 3

33 REFERENCES [] K. Ozawa, T. Okuzono, M. o, Japanese Journal of Appled Physcs 45 (6) [] V. K. Vanag, I. R. Epsten, Phys. Chem. Chem. Phys. (9) [3] M. o, Soft Matter Physcs, Oxford Unversty Press, 3. [4] M. abral, L. F. Francs, L. E. Scrven, AIChE Journal 48 () [5] E. B. Nauman, J. Savoca, AIChE Journal 47 () 6-. [6] J. S. Vrentas, C. M. Vrentas, J. Appl. Polym. Sc. 77 () [7] T. P. Lodge, Macromolecules 6 (983) [8] T. Chang, C. C. Han, L. M. Wheeler, T. P. Lodge, Macromolecules (988) [9] W. Brown, P. Zhou, Macromolecules 4 (99) [] P. avs, I. Snook, W. van Megen, B. N. Preston, W.. Comper, Macromolecules 7 (984) [] F. Roby, J.-F. Joanny, Macromolecules 5 (99) [] N. Nemoto, T. Inoue, Y. Makta, Y. Tsunashma, M. Kurata, Macromolecules 8 (985) [3] E. Gessler, A. M. Hecht, J. Chem. Phys. 67 (977) [4] Z. Sun and C. H. Wang, J. Chem. Phys. 6 (997) [5] J. S. Vrentas, J. L. uda, H. C. Lng, J. Appl. Polym. Sc. 3 (985), [6] R. J. Bearman, J. Phys. Chem. 65 (96) [7] S. S. Shojae, W. B. Krantz, A. R. Greenberg, J. Membrane Sc. 94 (994) [8] S. Alsoy, J. L. uda, AIChE Journal 45 (999) [9] J. M. Zelnsk, S. Alsoy, J. Polym. Sc., Part B: Polym. Phys. 39 () [] S. Shams Es-hagh, M. Cakmak, Polymer 79 (5) -8. [] S. R. e Groot, P. Mazur. Nonequlbrum Thermodynamcs. NY: over,

34 [] I. Gyarmat, Nonequlbrum Thermodynamcs Feld Theory and Varatonal Prncples, Sprnger-Verlag, 97. [3] I. Prgogne, Introducton to Thermodynamcs of Irreversble Processes. NY: Interscence, 967. [4] C. F. Curtss, R. Byron Brd, Ind. Eng. Chem. Res. 38 (999) [5] R. L. Scott, J. Chem. Phys. 7 (949) [6] R. L. Scott, J. Chem. Phys. 7 (949) [7] M. o, Introducton to Polymer Physcs, Oxford Unversty Press, New York, 996. [8] M. Kurata, Thermodynamcs of Polymer Solutons, Harwood Academc Publshers, 98. [9] V. Vtaglano, R. Sartoro, S. Scala,. Spaduzz, J. Soluton Chem. 7 (978) [3] V. J. Klenn, Thermodynamcs of Systems Contanng Flexble-Chan Polymers. Amsterdam, Elsever Scence B.V [3] H. Tompa, Polymer Solutons, Butterworths Scentfc Publcatons, 956. [3] P. E. Prce, Jr., I. H. Romdhane, AIChE Journal 49 (3)

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

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

More information

CHEMICAL REACTIONS AND DIFFUSION

CHEMICAL REACTIONS AND DIFFUSION CHEMICAL REACTIONS AND DIFFUSION A.K.A. NETWORK THERMODYNAMICS BACKGROUND Classcal thermodynamcs descrbes equlbrum states. Non-equlbrum thermodynamcs descrbes steady states. Network thermodynamcs descrbes

More information

4.2 Chemical Driving Force

4.2 Chemical Driving Force 4.2. CHEMICL DRIVING FORCE 103 4.2 Chemcal Drvng Force second effect of a chemcal concentraton gradent on dffuson s to change the nature of the drvng force. Ths s because dffuson changes the bondng n a

More information

Solution Thermodynamics

Solution 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 information

On the symmetric character of the thermal conductivity tensor

On the symmetric character of the thermal conductivity tensor On the symmetrc character of the thermal conductvty tensor Al R. Hadjesfandar Department of Mechancal and Aerospace Engneerng Unversty at Buffalo, State Unversty of New York Buffalo, NY 146 USA ah@buffalo.edu

More information

Lecture. Polymer Thermodynamics 0331 L Chemical Potential

Lecture. Polymer Thermodynamics 0331 L Chemical Potential Prof. Dr. rer. nat. habl. S. Enders Faculty III for Process Scence Insttute of Chemcal Engneerng Department of Thermodynamcs Lecture Polymer Thermodynamcs 033 L 337 3. Chemcal Potental Polymer Thermodynamcs

More information

Entropy generation in a chemical reaction

Entropy generation in a chemical reaction Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy

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

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

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

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Energy, Entropy, and Availability Balances Phase Equilibria. Nonideal Thermodynamic Property Models. Selecting an Appropriate Model

Energy, 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 information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

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

More information

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

( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.

( ) 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 information

Chemical Equilibrium. Chapter 6 Spontaneity of Reactive Mixtures (gases) Taking into account there are many types of work that a sysem can perform

Chemical Equilibrium. Chapter 6 Spontaneity of Reactive Mixtures (gases) Taking into account there are many types of work that a sysem can perform Ths chapter deals wth chemcal reactons (system) wth lttle or no consderaton on the surroundngs. Chemcal Equlbrum Chapter 6 Spontanety of eactve Mxtures (gases) eactants generatng products would proceed

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

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

Lecture 20: Noether s Theorem

Lecture 20: Noether s Theorem Lecture 20: Noether s Theorem In our revew of Newtonan Mechancs, we were remnded that some quanttes (energy, lnear momentum, and angular momentum) are conserved That s, they are constant f no external

More information

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

3. 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 information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

On the non-singularity of the thermal conductivity tensor and its consequences

On the non-singularity of the thermal conductivity tensor and its consequences On the non-sngularty of the thermal conductvty tensor and ts consequences Al R. Hadjesfandar Department of Mechancal and Aerospace Engneerng Unversty at Buffalo, State Unversty of New Yor Buffalo, NY 1426

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

More information

If two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9.

If 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 information

Difference Equations

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

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

Linear Approximation with Regularization and Moving Least Squares

Linear 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

Thermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak

Thermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak Thermodynamcs II Department of Chemcal Engneerng Prof. Km, Jong Hak Soluton Thermodynamcs : theory Obectve : lay the theoretcal foundaton for applcatons of thermodynamcs to gas mxture and lqud soluton

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2 Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

The Symmetries of Kibble s Gauge Theory of Gravitational Field, Conservation Laws of Energy-Momentum Tensor Density and the

The Symmetries of Kibble s Gauge Theory of Gravitational Field, Conservation Laws of Energy-Momentum Tensor Density and the The Symmetres of Kbble s Gauge Theory of Gravtatonal Feld, Conservaton aws of Energy-Momentum Tensor Densty and the Problems about Orgn of Matter Feld Fangpe Chen School of Physcs and Opto-electronc Technology,Dalan

More information

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram

Adiabatic 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 information

Lecture 7: Boltzmann distribution & Thermodynamics of mixing

Lecture 7: Boltzmann distribution & Thermodynamics of mixing Prof. Tbbtt Lecture 7 etworks & Gels Lecture 7: Boltzmann dstrbuton & Thermodynamcs of mxng 1 Suggested readng Prof. Mark W. Tbbtt ETH Zürch 13 März 018 Molecular Drvng Forces Dll and Bromberg: Chapters

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

e - c o m p a n i o n

e - c o m p a n i o n OPERATIONS RESEARCH http://dxdoorg/0287/opre007ec e - c o m p a n o n ONLY AVAILABLE IN ELECTRONIC FORM 202 INFORMS Electronc Companon Generalzed Quantty Competton for Multple Products and Loss of Effcency

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

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

More information

Introduction to Statistical Methods

Introduction to Statistical Methods Introducton to Statstcal Methods Physcs 4362, Lecture #3 hermodynamcs Classcal Statstcal Knetc heory Classcal hermodynamcs Macroscopc approach General propertes of the system Macroscopc varables 1 hermodynamc

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

A Self-Consistent Gibbs Excess Mixing Rule for Cubic Equations of State: derivation and fugacity coefficients

A 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 information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Lecture 5.8 Flux Vector Splitting

Lecture 5.8 Flux Vector Splitting Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form

More information

...Thermodynamics. If Clausius Clapeyron fails. l T (v 2 v 1 ) = 0/0 Second order phase transition ( S, v = 0)

...Thermodynamics. If Clausius Clapeyron fails. l T (v 2 v 1 ) = 0/0 Second order phase transition ( S, v = 0) If Clausus Clapeyron fals ( ) dp dt pb =...Thermodynamcs l T (v 2 v 1 ) = 0/0 Second order phase transton ( S, v = 0) ( ) dp = c P,1 c P,2 dt Tv(β 1 β 2 ) Two phases ntermngled Ferromagnet (Excess spn-up

More information

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure

CinChE 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 information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 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 information

Indeterminate pin-jointed frames (trusses)

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

More information

Supplemental document

Supplemental document Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON 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 information

Name: SID: Discussion Session:

Name: 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 information

NAME and Section No. it is found that 0.6 mol of O

NAME 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 information

Solution Thermodynamics

Solution Thermodynamics CH2351 Chemcal Engneerng Thermodynamcs II Unt I, II www.msubbu.n Soluton Thermodynamcs www.msubbu.n Dr. M. Subramanan Assocate Professor Department of Chemcal Engneerng Sr Svasubramanya Nadar College of

More information

and Statistical Mechanics Material Properties

and Statistical Mechanics Material Properties Statstcal Mechancs and Materal Propertes By Kuno TAKAHASHI Tokyo Insttute of Technology, Tokyo 15-855, JAPA Phone/Fax +81-3-5734-3915 takahak@de.ttech.ac.jp http://www.de.ttech.ac.jp/~kt-lab/ Only for

More information

Numerical Heat and Mass Transfer

Numerical 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 information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

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

More information

Thermodynamics General

Thermodynamics General Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,

More information

Appendix II Summary of Important Equations

Appendix II Summary of Important Equations W. M. Whte Geochemstry Equatons of State: Ideal GasLaw: Coeffcent of Thermal Expanson: Compressblty: Van der Waals Equaton: The Laws of Thermdynamcs: Frst Law: Appendx II Summary of Important Equatons

More information

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva Econ 39 - Statstcal Propertes of the OLS estmator Sanjaya DeSlva September, 008 1 Overvew Recall that the true regresson model s Y = β 0 + β 1 X + u (1) Applyng the OLS method to a sample of data, we estmate

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 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 information

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Electrical double layer: revisit based on boundary conditions

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

More information

Thermodynamics and statistical mechanics in materials modelling II

Thermodynamics and statistical mechanics in materials modelling II Course MP3 Lecture 8/11/006 (JAE) Course MP3 Lecture 8/11/006 Thermodynamcs and statstcal mechancs n materals modellng II A bref résumé of the physcal concepts used n materals modellng Dr James Ellott.1

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-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 information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 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 information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Research Article Green s Theorem for Sign Data

Research Article Green s Theorem for Sign Data Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of

More information

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

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

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Electrostatic Potential from Transmembrane Currents

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

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Temperature. Chapter Heat Engine

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

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

Three-Phase Distillation in Packed Towers: Short-Cut Modelling and Parameter Tuning

Three-Phase Distillation in Packed Towers: Short-Cut Modelling and Parameter Tuning European Symposum on Computer Arded Aded Process Engneerng 15 L. Pugjaner and A. Espuña (Edtors) 2005 Elsever Scence B.V. All rghts reserved. Three-Phase Dstllaton n Packed Towers: Short-Cut Modellng and

More information

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t 8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes

More information

CHEMICAL ENGINEERING

CHEMICAL ENGINEERING Postal Correspondence GATE & PSUs -MT To Buy Postal Correspondence Packages call at 0-9990657855 1 TABLE OF CONTENT S. No. Ttle Page no. 1. Introducton 3 2. Dffuson 10 3. Dryng and Humdfcaton 24 4. Absorpton

More information

Uncertainty and auto-correlation in. Measurement

Uncertainty and auto-correlation in. Measurement Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at

More information

Problem Points Score Total 100

Problem Points Score Total 100 Physcs 450 Solutons of Sample Exam I Problem Ponts Score 1 8 15 3 17 4 0 5 0 Total 100 All wor must be shown n order to receve full credt. Wor must be legble and comprehensble wth answers clearly ndcated.

More information

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA

More information

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity Int. Journal of Math. Analyss, Vol. 6, 212, no. 22, 195-114 Unqueness of Weak Solutons to the 3D Gnzburg- Landau Model for Superconductvty Jshan Fan Department of Appled Mathematcs Nanjng Forestry Unversty

More information

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Hidden Markov Models & The Multivariate Gaussian (10/26/04) CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE III LECTURE - 2 EXPERIMENTAL DESIGN MODELS Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 2 We consder the models

More information

Determination of Structure and Formation Conditions of Gas Hydrate by Using TPD Method and Flash Calculations

Determination 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 information

An Improved Model for the Droplet Size Distribution in Sprays Developed From the Principle of Entropy Generation maximization

An Improved Model for the Droplet Size Distribution in Sprays Developed From the Principle of Entropy Generation maximization ILASS Amercas, 9 th Annual Conference on Lqud Atomzaton and Spray Systems, oronto, Canada, May 6 An Improved Model for the Droplet Sze Dstrbuton n Sprays Developed From the Prncple of Entropy Generaton

More information

Module 3: The Whole-Process Perspective for Thermochemical Hydrogen

Module 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 information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Experment-I MODULE VII LECTURE - 3 ANALYSIS OF COVARIANCE Dr Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Any scentfc experment s performed

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

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse

More information

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980 MT07: Multvarate Statstcal Methods Mke Tso: emal mke.tso@manchester.ac.uk Webpage for notes: http://www.maths.manchester.ac.uk/~mkt/new_teachng.htm. Introducton to multvarate data. Books Chat eld, C. and

More information

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013 ISSN: 2277-375 Constructon of Trend Free Run Orders for Orthogonal rrays Usng Codes bstract: Sometmes when the expermental runs are carred out n a tme order sequence, the response can depend on the run

More information