Evaluating Thermodynamic Properties in LAMMPS

Size: px
Start display at page:

Download "Evaluating Thermodynamic Properties in LAMMPS"

Transcription

1 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle Evaluatng hermodynamc Proertes n LAMMP Davd Keffer Deartment of Materals cence & Engneerng Unversty of ennessee Knoxvlle date begun: March 06 date last udated: March 4 06 able of ontents I. Purose of Document... II. ecfyng the hermodynamc tate... III. Evaluatng hermodynamc Proertes... III.A. Instantaneous Proertes... number of artcles of comonent... mole or atom fracton of comonent x... volume of the system... temerature... ressure... 4 knetc energy KE... 5 otental energy PE... 5 nternal energy U (or total energy E)... 6 enthaly H... 6 III.B. Fluctuaton Proertes (heat caacty sothermal comressblty etc.)... 6 III.. Proertes ncludng Entroy... 8 entroy... 8 free energes A and G... 9 chemcal otental of comonent µ... 9 entroy dfferences... 9

2 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle I. Purose of Document he urose of ths document s to rovde a ractcal ntroducton to the evaluaton of thermodynamc roertes n LAMMP. he notes begn wth some formal theory and conclude wth ractcal mlementaton. From an undergraduate course n thermodynamcs we have a lst of thermodynamc roertes of nterest. Here we revew them n three categores: () nstantaneous roertes mmedately avalable n LAMMP () roertes that can be relably obtaned wth a lttle addtonal work and () roertes nvolvng entroy. II. ecfyng the hermodynamc tate We begn by recallng the Phase Rule that states that the number of thermodynamc degrees of freedom DOF of a system are related to the number of hases (φ) and the number of comonents. DOF φ + () For a sngle comonent sngle hase system DOF meanng that the thermodynamc state s fully secfed by two varables for examle temerature and ressure or and the number densty ρ. In fact any two varables can be used to secfy the state. For examle a state could be secfed by rovdng the enthaly H and the ρ but that s just not very ractcal. he Phase Rule s based on ntensve roertes thus t does not consder the secfcaton of the number of atoms or molecules as one of the state varables. In the smulaton of a sngle comonent sngle hase system we therefore secfy three varables n order to fully secfy the system. Below we lst the four routnely accessble ensembles based on the choce of varables. In the common notaton smulaton volume an extensve varable s secfed but ths does not trouble us snce the number densty s related to the system volume va ρ. E mcrocanoncal ensemble canoncal ensemble sobarc-sothermal ensemble H sobarc-senthalc ensemble In undergraduate thermodynamc courses we ordnarly use the symbol U to reresent the nternal energy whereas n statstcal mechancs the symbol E s used. Here U E. In a multcomonent sngle-hase system the hase rule rovdes that each addtonal seces adds an addtonal degrees of thermodynamc freedom DOF +. ycally ths means that we secfy the number of each tye of atom or molecule whch secfes not but only - degrees of freedom. It s easest to see ths n terms of the ntensve resentaton of the of each seces namely the mole fracton (or equvalently atom fracton) of comonent x. x j j ()

3 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle nce the mole fractons are constraned to sum to unty x j j () only - of the mole fractons are ndeendent. hus there reman two addtonal thermodynamc degrees of freedom such as and or and ρ. As you ncrease the number of hases the thermodynamc degrees of freedom decrease. In a sngle comonent two-hase system DOF. ecfyng one varable say the temerature fully defnes the thermodynamc state. In terms of vaor-lqud equlbrum secfyng the bolng temerature secfes the vaor ressure whch s the ressure at whch the vaor and lqud can co-exst at that temerature. III. Evaluatng hermodynamc Proertes III.A. Instantaneous Proertes he soluton of the classcal equatons of moton results exclusvely n a set of trajectores whch are comrsed of atomc ostons r and veloctes v (or equvalently momenta ) as a functon of tme t. Every thermodynamc roerty s formulated as a functon of these trajectores. nce classcal thermodynamcs s a contnuum theory formulated wthout knowledge of the atomc nature of matter t does not rovde relatonshs between thermodynamc quanttes and atomc trajectores. For that we requre statstcal mechancs sometmes called statstcal thermodynamcs when lmted to equlbrum systems. We make no attemt here to derve any of these relatonshs. We smly resent some relatonshs and lnk them to varables avalable n LAMMP. number of artcles of comonent In the four ensembles descrbed above the number of artcles (atoms or molecules) s held constant. herefore ths extensve varable s always an nut to the smulaton. It s a constant and has no fluctuaton durng the course of the smulaton. mole or atom fracton of comonent x nce the number of artcles s a fxed nut er equaton () the mole fractons are also fxed nuts. volume of the system In the E and ensembles the volume of the smulaton cell s a fxed nut. It does not fluctuate wth tme. In the and H ensembles the volume fluctuates n resonse to the barostat. temerature he temerature s related to the knetc energy KE va the equartton theorem. he knetc energy s related to the veloctes va classcal mechancs.

4 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle k KE m v (4) Unlke the volume the temerature fluctuates n all four ensembles. he and ensembles n whch the temerature s secfed are controlled by thermostats that do not kee the temerature exactly at the set ont temerature at all tmes. Instead a roer thermostat allows the temerature to fluctuate about the set ont temerature n a manner that s stll caable of rgorously generatng trajectores that corresond to the ensemble. ressure In three-dmensons the ressure tensor s a symmetrc x tensor wth each comonent defned by a knetc and otental contrbuton KE + (5) β β PE β he knetc contrbuton to the β element of the ressure s KE β mv v β (6) nce the veloctes n dfferent dmensons are not correlated when β the average value of the s zero. However β we have twce the knetc energy n that dmenson. KE β δ β KE β mv (7) If use the equartton theorem relate the knetc energy n terms of the temerature and average over all three dmensons we have the ressure from the deal gas law where the otental energy s zero k KE m v for deal gas only (8) hus the otental energy contrbuton to the ressure can be seen as a devaton from the deal gas. PE β r f β (9) 4

5 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle When β the average value of the PE β s not necessarly zero deendng on the state of stress wthn the system. If the ressure tensor s not of nterest the scalar ressure reorted by LAMMP s the thermodynamc ressure r f (0) mv + Lke the temerature the ressure fluctuates n all four ensembles. he and H ensembles n whch the ressure s secfed are controlled by barostats that do not kee the ressure exactly at the set ont temerature at all tmes. Instead a roer barostat allows the ressure to fluctuate about the set ont ressure n a manner that s stll caable of rgorously generatng trajectores that corresond to the ensemble. knetc energy KE he knetc energy s not a tradtonal thermodynamc varable. However t s an mortant roerty n the atomc-level descrton of matter. he KE s defned by the sum of the knetc energy of all artcles n the system as gven n equaton (4). In ths work the knetc energy s defned exclusvely as a functon of veloctes. Exlctly the knetc energy s not a functon of artcle ostons. If the system does not have net zero momentum then there s a center of mass velocty that exsts n the smulaton. he knetc energy can then be broken nto two arts. Knetc energy of the center of mass of the system and the ecular knetc energy whch s the sum of knetc energy of each artcle based on veloctes relatve to the system center of mass velocty. Imortantly when the equartton theorem s used to calculate temerature t reles on ecular veloctes. ee the notes on thermostats to roerly account for ths effect. otental energy PE he otental energy s not a tradtonal thermodynamc varable. However t s an mortant roerty n the atomc-level descrton of matter. he PE s defned by the sum of all nteractons between all artcles n the system as gven by the nteracton otental. he Lennard-Jones otental Embedded-Atom Method otental and OPL-aa otental are all examles of nteracton otentals that defne the otental energy. In ths work the otental energy s defned exclusvely as a functon of ostons. Exlctly the otental energy s not a functon of artcle veloctes. 5

6 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle nternal energy U (or total energy E) he nternal energy s the sum of the knetc and otental energes. It s also called the total energy. Often n thermodynamcs we here of the excess nternal energy whch subtracts the KE form the total energy. However n general the nternal energy ncludes both contrbutons. U E KE + PE () In the E ensemble the nternal energy should be constant as t s a measure of the conservaton of energy of the classcal equatons of moton. It truth t fluctuates based on the numercal aroxmatons to the soluton of these ordnary dfferental equatons. In the other three ensembles the nternal energy fluctuates. enthaly H he enthaly s defned as the sum of the nternal energy and the roduct of ressure and volume. H U + () All of the quanttes lsted n ths secton can evaluated at each nstant n tme yeldng so called nstantaneous roertes. Of course what s tycally mortant s to reort the average roertes. he average can be comuted n a ost-rocessng mode usng the nstantaneous values rnted out n the LAMMP outut fle. Alternatvely the averages can be comuted n the LAMMP nut fle usng the fx ave/tme command. ee for examle the block averagng examles on the course webste that comute both the average as well as the standard error usng the fx ave/tme command. III.B. Fluctuaton Proertes (heat caacty sothermal comressblty etc.) here are a varety of thermodynamc roertes that are based on the artal dervatve of roerty X wth resect to roerty Y holdng roerty Z constant. ome common thermodynamc roertes of ths form nclude the constant-ressure heat caacty sothermal comressblty β and the thermal exansvty. the H (.a) β (.b) (.c) 6

7 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle he exresson gven above are thermodynamc defntons. tatstcal mechancs rovdes an elegant theoretcal way to evaluate these roertes based on the fluctuatons of roertes. hese fluctuatons are based on the varance or covarance of the nstantaneous thermodynamc roertes dscussed n the revous secton. For examle n the ensemble the constant ressure heat caacty s roortonal to the varance of the enthaly < σ > H (4.a) k he sothermal comressblty s obtaned through the varance of volume β < σ > (4.b) k he thermal exanson coeffcent s based on the covarance of volume and enthaly of the < σ H k ( )( ) > (4.c) he angled brackets ndcate an ensemble average namely an average over both all artcles and tme. ote that ths varance s based on nstantaneous values corresondng to the square of standard devaton of the varables. Exlctly these fluctuatons do not corresond to the square of the standard error generated by the block averagng rocedure. hese fluctuatons seem to rovde an effcent method to obtan these thermodynamc roertes snce we can obtan all these three roertes from one sngle MD smulaton at the desred state ont. he formulae gven n equaton (4) are for the ensemble. Other formulae exst for other ensembles. We could also calculate β and usng a second rocedure namely centered fnte dfference usng average values of roertes over multle smulatons. In ths aroach one runs three smulatons for each thermodynamc roerty. For examle for the heat caacty one runs a smulaton at the state ont defned by () then two other smulatons at (+ε) and (-ε). he choce of ε s made so that one can observe statstcally vald dfferences n the enthaly between the state onts. hus the heat caacty can be determned as ( + ε ) H ( ε ) H (5.a) ε mlarly the sothermal comressblty and thermal exansvty can be exressed as ( ) ( + ε ) ( ε ) β (5.b) ε 7

8 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle ( ) ( + ε ) ( ε ) ε (5.c) A detaled comarson of these aroaches s gven n the followng reference. Wang Q. Keffer D.J. Petrovan. homas J.B. Molecular Dynamcs mulaton of Polyethylene erehthalate Olgomers J. Phys. hem. B. 4() do: 0.0/j90976j. her analyss found that the fnte dfference aroach was markedly more accurate. Any advantage thought to be ganed by effcency of obtanng the result from a sngle smulaton was lost snce n order to get reasonable statstcs from the fluctuatons one would have to run much larger systems for much longer tmes. ertanly a set of three smulatons could be set u n a sngle LAMMP nut fle to evaluate these roertes va the centered-fnte dfference aroach. hermodynamc relatons can also be exloted to evaluate other fluctuaton roertes. For examle the constant volume heat caacty s related to the above three roertes va v v (6.a) β and the sentroc comressblty can be obtaned from the followng relaton v β (6.b) β s III.. Proertes ncludng Entroy entroy here are many reasons to want to know the entroy of a system. here exst some aroxmate ways to evaluate an entroy based on an average confguraton. ee for examle the followng two references. D.. Wallace On the role of densty-fluctuatons n the entroy of a flud J. hem. Phys. 87 (4) (987) A. Baranya D.J. Evans Drect entroy calculaton from comuter smulaton of lquds Phys. Rev. A 40 (7) (989) 87 8 However these efforts are fraught wth suffcent dffculty and errors that they have not been acceted as a standard way to comute entroy for arbtrarly comlex systems as s the case for the nstantaneous thermodynamc roertes. herefore today there s no easly mlemented and straghtforward way to calculate the entroy n a classcal smulaton. 8

9 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle 9 free energes A and G nce the entroy s not easly avalable there s also no drect way to calculate ether the Helmholtz or the Gbbs free energy from a smulaton. U A (7.a) U H G + (7.b) You should be suscous of any artcle that clams to have done so. chemcal otental of comonent µ he chemcal otental of comonent s defned as a thermodynamc artal dervatve j j A G µ (8) here are some sohstcated ways to comute the chemcal otental n a smulaton. We leave these technques to a searate lecture to be resented later n the course. entroy dfferences hose of you wth a fondness for thermodynamcs may recall Maxwell relatons. wo Maxwell relatons most relevant to the calculaton of entroy dfferences are (9.a) (9.b) hese equatons can be manulated and ntegrated to yeld ( ) ( ) d d (0.a) ( ) ( ) d d (0.b)

10 D. Keffer ME 64 Det. of Materals cence & Engneerng Unversty of ennessee Knoxvlle In the frst case one can erform a seres of smulatons all at the same temerature n whch the volume vares from to. One can then comute the value of across ths range usng for examle the centered-fnte dfference aroach descrbed above for each smulaton at ( ). From ths data one can use a numercal ntegraton scheme (such as the traezodal rule) to calculate an aroxmate value of the entroy dfference. ote that ths requres a seres of smulatons. In the second case one can erform a seres of smulatons all at the same temerature n whch the ressure vares from to. One can then comute the value of across ths range usng for examle the centered-fnte dfference aroach descrbed above for at each smulaton ( ). From ths data one can agan use a numercal ntegraton scheme (such as the traezodal rule) to calculate an aroxmate value of the entroy dfference. ote that ths requres a seres of smulatons. Other thermodynamc relatons can also be exloted to calculate entroy dfferences such as v (.a) (.b) As before these equatons can be manulated and ntegrated to yeld ( ) v d d (.a) ( ) ( ) ( ) d d (.b) Agan a seres of smulatons can be run ether at constant volume or constant ressure from to. If the arorate heat caacty s evaluated for each smulaton then the ntegral yelds the change n entroy. Once the change n entroy s known the change n free energy can be comuted drectly from equaton (7) snce the same smulatons generated values for the nternal energy and enthaly. One (of many) examles where such an ntegraton s done (usng equaton.a) s gven n Fgure.a. and Fgure.b of the followng reference. Mcutt.W. Wang Q. Keffer D.J. Ros O. Entroy-drven tructure and Dynamcs n arbon anocrystalltes J. anoartcle Res. 6 (4) 04 artcle # 65 do: 0.007/s

Non-Ideality Through Fugacity and Activity

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

6. Hamilton s Equations

6. Hamilton s Equations 6. Hamlton s Equatons Mchael Fowler A Dynamcal System s Path n Confguraton Sace and n State Sace The story so far: For a mechancal system wth n degrees of freedom, the satal confguraton at some nstant

More information

THERMODYNAMICS. Temperature

THERMODYNAMICS. Temperature HERMODYNMICS hermodynamcs s the henomenologcal scence whch descrbes the behavor of macroscoc objects n terms of a small number of macroscoc arameters. s an examle, to descrbe a gas n terms of volume ressure

More information

Digital PI Controller Equations

Digital PI Controller Equations Ver. 4, 9 th March 7 Dgtal PI Controller Equatons Probably the most common tye of controller n ndustral ower electroncs s the PI (Proortonal - Integral) controller. In feld orented motor control, PI controllers

More information

ME 440 Aerospace Engineering Fundamentals

ME 440 Aerospace Engineering Fundamentals Fall 006 ME 440 Aerosace Engneerng Fundamentals roulson hrust Jet Engne F m( & Rocket Engne F m & F ρ A - n ) ρ A he basc rncle nsde the engne s to convert the ressure and thermal energy of the workng

More information

Chapter 8 Balances on Nonreactive Processes 8.1 Elements of Energy Balances Calculations 8.1a Reference States A Review

Chapter 8 Balances on Nonreactive Processes 8.1 Elements of Energy Balances Calculations 8.1a Reference States A Review Chater 8 Balances on Nonreactve Processes 8.1 Elements of Energy Balances Calculatons 8.1a Reference States A Revew We can never know the absolute values of U and H for a seces at a gven state. t Fortunately,

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

Mechanics Physics 151

Mechanics Physics 151 Mechancs hyscs 151 Lecture Canoncal Transformatons (Chater 9) What We Dd Last Tme Drect Condtons Q j Q j = = j, Q, j, Q, Necessary and suffcent j j for Canoncal Transf. = = j Q, Q, j Q, Q, Infntesmal CT

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

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

CHE-201. I n t r o d u c t i o n t o Chemical E n g i n e e r i n g. C h a p t e r 6. Multiphase Systems

CHE-201. I n t r o d u c t i o n t o Chemical E n g i n e e r i n g. C h a p t e r 6. Multiphase Systems I n t r o d u c t o n t o Chemcal E n g n e e r n g CHE-201 I N S T R U C T O R : D r. N a b e e l S a l m b o - G h a n d e r C h a t e r 6 Multhase Systems Introductory Statement: Phase s a state of

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 151 Lecture 22 Canoncal Transformatons (Chater 9) What We Dd Last Tme Drect Condtons Q j Q j = = j P, Q, P j, P Q, P Necessary and suffcent P j P j for Canoncal Transf. = = j Q, Q, P j

More information

Adsorption: A gas or gases from a mixture of gases or a liquid (or liquids) from a mixture of liquids is bound physically to the surface of a solid.

Adsorption: A gas or gases from a mixture of gases or a liquid (or liquids) from a mixture of liquids is bound physically to the surface of a solid. Searatons n Chemcal Engneerng Searatons (gas from a mxture of gases, lquds from a mxture of lquds, solds from a soluton of solds n lquds, dssolved gases from lquds, solvents from gases artally/comletely

More information

Statistical mechanics handout 4

Statistical mechanics handout 4 Statstcal mechancs handout 4 Explan dfference between phase space and an. Ensembles As dscussed n handout three atoms n any physcal system can adopt any one of a large number of mcorstates. For a quantum

More information

Lecture 3 Examples and Problems

Lecture 3 Examples and Problems Lecture 3 Examles and Problems Mechancs & thermodynamcs Equartton Frst Law of Thermodynamcs Ideal gases Isothermal and adabatc rocesses Readng: Elements Ch. 1-3 Lecture 3, 1 Wllam Thomson (1824 1907) a.k.a.

More information

EP523 Introduction to QFT I

EP523 Introduction to QFT I EP523 Introducton to QFT I Toc 0 INTRODUCTION TO COURSE Deartment of Engneerng Physcs Unversty of Gazante Setember 2011 Sayfa 1 Content Introducton Revew of SR, QM, RQM and EMT Lagrangan Feld Theory An

More information

Hans-Joachim Kretzschmar and Katja Knobloch

Hans-Joachim Kretzschmar and Katja Knobloch Sulementary Backward Equatons for the Industral Formulaton IAPWS-IF of Water and Steam for Fast Calculatons of Heat Cycles, Bolers, and Steam Turbnes Hans-Joachm Kretzschmar and Katja Knobloch Deartment

More information

Lecture 14: Forces and Stresses

Lecture 14: Forces and Stresses The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal

More information

EXCESS FUNCTIONS MATHEMATICAL EXPLANATION

EXCESS FUNCTIONS MATHEMATICAL EXPLANATION EXCESS FUNCTIONS Excess functons... 1 Mathematcal exlanaton... 1 Physcal examle: Molar volume of a ure substance... 2 Chemcal otental for a ure substance... 3 Actvty, actvty coeffcent and fugacty (ure

More information

Statistical Evaluation of WATFLOOD

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

More information

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

A Quadratic Cumulative Production Model for the Material Balance of Abnormally-Pressured Gas Reservoirs F.E. Gonzalez M.S.

A Quadratic Cumulative Production Model for the Material Balance of Abnormally-Pressured Gas Reservoirs F.E. Gonzalez M.S. Formaton Evaluaton and the Analyss of Reservor Performance A Quadratc Cumulatve Producton Model for the Materal Balance of Abnormally-Pressured as Reservors F.E. onale M.S. Thess (2003) T.A. Blasngame,

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

MOLECULAR DYNAMICS ,..., What is it? 2 = i i

MOLECULAR DYNAMICS ,..., What is it? 2 = i i MOLECULAR DYNAMICS What s t? d d x t 2 m 2 = F ( x 1,..., x N ) =1,,N r ( x1 ( t),..., x ( t)) = v = ( x& 1 ( t ),..., x& ( t )) N N What are some uses of molecular smulatons and modelng? Conformatonal

More information

ENV05: METROLOGY FOR OCEAN SALINITY AND ACIDITY

ENV05: METROLOGY FOR OCEAN SALINITY AND ACIDITY Kck-off meetng EMRP ENV5 Ocean Metrology - Setember, PB Berln ENV5: MEROLOGY FOR OCEAN SALINIY AND ACIDIY WP: Extenson of measurement range for thermohyscal arameters (INRM, PB) ISIUO MEROLOGICA Smona

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

A Quadratic Cumulative Production Model for the Material Balance of Abnormally-Pressured Gas Reservoirs F.E. Gonzalez M.S.

A Quadratic Cumulative Production Model for the Material Balance of Abnormally-Pressured Gas Reservoirs F.E. Gonzalez M.S. Natural as Engneerng A Quadratc Cumulatve Producton Model for the Materal Balance of Abnormally-Pressured as Reservors F.E. onale M.S. Thess (2003) T.A. Blasngame, Texas A&M U. Deartment of Petroleum Engneerng

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

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

12. The Hamilton-Jacobi Equation Michael Fowler

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

STATISTICAL MECHANICS

STATISTICAL MECHANICS STATISTICAL MECHANICS Thermal Energy Recall that KE can always be separated nto 2 terms: KE system = 1 2 M 2 total v CM KE nternal Rgd-body rotaton and elastc / sound waves Use smplfyng assumptons KE of

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

Departure Process from a M/M/m/ Queue

Departure Process from a M/M/m/ Queue Dearture rocess fro a M/M// Queue Q - (-) Q Q3 Q4 (-) Knowledge of the nature of the dearture rocess fro a queue would be useful as we can then use t to analyze sle cases of queueng networs as shown. The

More information

Lecture 4. Macrostates and Microstates (Ch. 2 )

Lecture 4. Macrostates and Microstates (Ch. 2 ) Lecture 4. Macrostates and Mcrostates (Ch. ) The past three lectures: we have learned about thermal energy, how t s stored at the mcroscopc level, and how t can be transferred from one system to another.

More information

between standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we

between standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we hermodynamcs, Statstcal hermodynamcs, and Knetcs 4 th Edton,. Engel & P. ed Ch. 6 Part Answers to Selected Problems Q6.. Q6.4. If ξ =0. mole at equlbrum, the reacton s not ery far along. hus, there would

More information

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding. Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

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

More information

Problem 10.1: One-loop structure of QED

Problem 10.1: One-loop structure of QED Problem 10.1: One-loo structure of QED In Secton 10.1 we argued form general rncles that the hoton one-ont and three-ont functons vansh, whle the four-ont functon s fnte. (a Verfy drectly that the one-loo

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

The Dirac Equation for a One-electron atom. In this section we will derive the Dirac equation for a one-electron atom.

The Dirac Equation for a One-electron atom. In this section we will derive the Dirac equation for a one-electron atom. The Drac Equaton for a One-electron atom In ths secton we wll derve the Drac equaton for a one-electron atom. Accordng to Ensten the energy of a artcle wth rest mass m movng wth a velocty V s gven by E

More information

coordinates. Then, the position vectors are described by

coordinates. Then, the position vectors are described by Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

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

More information

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

Review of Thermodynamics

Review of Thermodynamics Revew of hermodynamcs from Statstcal Physcs usng Mathematca James J. Kelly, 1996-2002 We revew the laws of thermodynamcs and some of the technques for dervaton of thermodynamc relatonshs. Introducton Equlbrum

More information

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

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

More information

Isothermal vs. adiabatic compression

Isothermal vs. adiabatic compression Isothermal vs. adabatc comresson 1. One and a half moles of a datomc gas at temerature 5 C are comressed sothermally from a volume of 0.015 m to a volume of 0.0015 m. a. Sketch the rocess on a dagram and

More information

Confidence intervals for weighted polynomial calibrations

Confidence intervals for weighted polynomial calibrations Confdence ntervals for weghted olynomal calbratons Sergey Maltsev, Amersand Ltd., Moscow, Russa; ur Kalambet, Amersand Internatonal, Inc., Beachwood, OH e-mal: kalambet@amersand-ntl.com htt://www.chromandsec.com

More information

Combinational Circuit Design

Combinational Circuit Design Combnatonal Crcut Desgn Part I: Desgn Procedure and Examles Part II : Arthmetc Crcuts Part III : Multlexer, Decoder, Encoder, Hammng Code Combnatonal Crcuts n nuts Combnatonal Crcuts m oututs A combnatonal

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

Global Sensitivity. Tuesday 20 th February, 2018

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

More information

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

Review of Classical Thermodynamics

Review of Classical Thermodynamics Revew of Classcal hermodynamcs Physcs 4362, Lecture #1, 2 Syllabus What s hermodynamcs? 1 [A law] s more mpressve the greater the smplcty of ts premses, the more dfferent are the knds of thngs t relates,

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

Problem Set 9 Solutions

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

More information

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

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

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

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

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

Spring Force and Power

Spring Force and Power Lecture 13 Chapter 9 Sprng Force and Power Yeah, energy s better than orces. What s net? Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi IN THIS CHAPTER, you wll learn how to solve problems

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

Logistic regression with one predictor. STK4900/ Lecture 7. Program

Logistic regression with one predictor. STK4900/ Lecture 7. Program Logstc regresson wth one redctor STK49/99 - Lecture 7 Program. Logstc regresson wth one redctor 2. Maxmum lkelhood estmaton 3. Logstc regresson wth several redctors 4. Devance and lkelhood rato tests 5.

More information

A quote of the week (or camel of the week): There is no expedience to which a man will not go to avoid the labor of thinking. Thomas A.

A quote of the week (or camel of the week): There is no expedience to which a man will not go to avoid the labor of thinking. Thomas A. A quote of the week (or camel of the week): here s no expedence to whch a man wll not go to avod the labor of thnkng. homas A. Edson Hess law. Algorthm S Select a reacton, possbly contanng specfc compounds

More information

THE SUMMATION NOTATION Ʃ

THE SUMMATION NOTATION Ʃ Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the

More information

The Decibel and its Usage

The Decibel and its Usage The Decbel and ts Usage Consder a two-stage amlfer system, as shown n Fg.. Each amlfer rodes an ncrease of the sgnal ower. Ths effect s referred to as the ower gan,, of the amlfer. Ths means that the sgnal

More information

Lecture # 02: Pressure measurements and Measurement Uncertainties

Lecture # 02: Pressure measurements and Measurement Uncertainties AerE 3L & AerE343L Lecture Notes Lecture # 0: Pressure measurements and Measurement Uncertantes Dr. Hu H Hu Deartment of Aerosace Engneerng Iowa State Unversty Ames, Iowa 500, U.S.A Mechancal Pressure

More information

AN ASYMMETRIC GENERALIZED FGM COPULA AND ITS PROPERTIES

AN ASYMMETRIC GENERALIZED FGM COPULA AND ITS PROPERTIES Pa. J. Statst. 015 Vol. 31(1), 95-106 AN ASYMMETRIC GENERALIZED FGM COPULA AND ITS PROPERTIES Berzadeh, H., Parham, G.A. and Zadaram, M.R. Deartment of Statstcs, Shahd Chamran Unversty, Ahvaz, Iran. Corresondng

More information

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson Publcaton 2006/01 Transport Equatons n Incompressble URANS and LES Lars Davdson Dvson of Flud Dynamcs Department of Appled Mechancs Chalmers Unversty of Technology Göteborg, Sweden, May 2006 Transport

More information

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation 1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve

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

Monte Carlo method II

Monte Carlo method II Course MP3 Lecture 5 14/11/2006 Monte Carlo method II How to put some real physcs nto the Monte Carlo method Dr James Ellott 5.1 Monte Carlo method revsted In lecture 4, we ntroduced the Monte Carlo (MC)

More information

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

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

More information

STATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1

STATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1 CHE 54 A. Panagotopoulos STATSTCAL MECHACAL ESEMBLES MCROSCOPC AD MACROSCOPC ARABLES The central queston n Statstcal Mechancs can be phrased as follows: f partcles (atoms, molecules, electrons, nucle,

More information

Physics 240: Worksheet 30 Name:

Physics 240: Worksheet 30 Name: (1) One mole of an deal monatomc gas doubles ts temperature and doubles ts volume. What s the change n entropy of the gas? () 1 kg of ce at 0 0 C melts to become water at 0 0 C. What s the change n entropy

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

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

Taylor series coefficients of the HP-polynomial as an invariant for links in the solid torus

Taylor series coefficients of the HP-polynomial as an invariant for links in the solid torus Al. Math. Inf. c. 7, No. 1, 23-28 (213) 23 Aled Mathematcs & Informaton cences An Internatonal Journal c 213 NP Natural cences Publshng Cor. aylor seres coeffcents of the HP-olynomal as an nvarant for

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

Property calculation I

Property calculation I .0, 3.0, 0.333,.00 Introducton to Modelng and Smulaton Sprng 0 Part I Contnuum and partcle methods Property calculaton I Lecture 3 Markus J. Buehler Laboratory for Atomstc and Molecular Mechancs Department

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 II LECTURE - GENERAL LINEAR HYPOTHESIS AND ANALYSIS OF VARIANCE Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 3.

More information

NUMERICAL SIMULATION OF SUDDEN-EXPANSION PARTICLE-LADEN FLOWS USING THE EULERIAN-LAGRANGIAN APPROACH

NUMERICAL SIMULATION OF SUDDEN-EXPANSION PARTICLE-LADEN FLOWS USING THE EULERIAN-LAGRANGIAN APPROACH THERMAL SCIENCE: Year 2012, Vol. 16, No. 4,. 1005-1012 1005 NUMERICAL SIMULATION OF SUDDEN-EXPANSION PARTICLE-LADEN FLOWS USING THE EULERIAN-LAGRANGIAN APPROACH by Mohamed Al MERGHENI a, b*, Jean-Charles

More information

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

2016 Wiley. Study Session 2: Ethical and Professional Standards Application 6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton

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

Topology optimization of plate structures subject to initial excitations for minimum dynamic performance index

Topology optimization of plate structures subject to initial excitations for minimum dynamic performance index th World Congress on Structural and Multdsclnary Otmsaton 7 th -2 th, June 25, Sydney Australa oology otmzaton of late structures subject to ntal exctatons for mnmum dynamc erformance ndex Kun Yan, Gengdong

More information

PETE 310 Lectures # 24 & 25 Chapter 12 Gas Liquid Equilibrium

PETE 310 Lectures # 24 & 25 Chapter 12 Gas Liquid Equilibrium ETE 30 Lectures # 24 & 25 Chapter 2 Gas Lqud Equlbrum Thermal Equlbrum Object A hgh T, Object B low T Intal contact tme Intermedate tme. Later tme Mechancal Equlbrum ressure essels Vale Closed Vale Open

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

ReaxFF potential functions

ReaxFF potential functions ReaxFF otental functons Suortng nformaton for the manuscrt Nelson, K.D., van Dun, A.C.T., Oxgaard, J.., Deng, W. and Goddard III, W.A. Develoment of the ReaxFF reactve force feld for descrbng transton

More information

Linear system of the Schrödinger equation Notes on Quantum Mechanics

Linear system of the Schrödinger equation Notes on Quantum Mechanics Lnear sstem of the Schrödnger equaton Notes on Quantum Mechancs htt://quantum.bu.edu/notes/quantummechancs/lnearsstems.df Last udated Wednesda, October 9, 003 :0:08 Corght 003 Dan Dll (dan@bu.edu) Deartment

More information

Modeling of diffusion and chemical reactions in ion-exchange resins during swelling and shrinking

Modeling of diffusion and chemical reactions in ion-exchange resins during swelling and shrinking Modelng of dffuson and chemcal reactons n on-exchange resns durng swellng and shrnkng Tuomo ano * and Erkk Paatero, Laeenranta Unversty of Technology, Laboratory of Industral Chemstry, knnarlankatu 34,

More information

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics) CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O

More information

The ChemSep Book. Harry A. Kooijman Consultant. Ross Taylor Clarkson University, Potsdam, New York University of Twente, Enschede, The Netherlands

The 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 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

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

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

Poisson brackets and canonical transformations

Poisson brackets and canonical transformations rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order

More information

Implicit Integration Henyey Method

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

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1 Problem Set 4 Suggested Solutons Problem (A) The market demand functon s the soluton to the followng utlty-maxmzaton roblem (UMP): The Lagrangean: ( x, x, x ) = + max U x, x, x x x x st.. x + x + x y x,

More information

The Robustness of a Nash Equilibrium Simulation Model

The Robustness of a Nash Equilibrium Simulation Model 8th World IMACS / MODSIM Congress, Carns, Australa 3-7 July 2009 htt://mssanz.org.au/modsm09 The Robustness of a Nash Equlbrum Smulaton Model Etaro Ayosh, Atsush Mak 2 and Takash Okamoto 3 Faculty of Scence

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

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

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

More information