A Helmholtz energy equation of state for calculating the thermodynamic properties of fluid mixtures

Similar documents
PREDICTION OF VAPOR-LIQUID EQUILIBRIA OF BINARY MIXTURES USING QUANTUM CALCULATIONS AND ACTIVITY COEFFICIENT MODELS

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

Chapter Business Statistics: A First Course Fifth Edition. Learning Objectives. Correlation vs. Regression. In this chapter, you learn:

Block-Based Compact Thermal Modeling of Semiconductor Integrated Circuits

Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions

COV. Violation of constant variance of ε i s but they are still independent. The error term (ε) is said to be heteroscedastic.

ENGI 3423 Simple Linear Regression Page 12-01

Functions of Random Variables

Objectives of Multiple Regression

CHAPTER VI Statistical Analysis of Experimental Data

ENGI 4421 Propagation of Error Page 8-01

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Bayesian Classification. CS690L Data Mining: Classification(2) Bayesian Theorem: Basics. Bayesian Theorem. Training dataset. Naïve Bayes Classifier

Chapter Statistics Background of Regression Analysis

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

Lecture 8: Linear Regression

Third handout: On the Gini Index

Beam Warming Second-Order Upwind Method

Journal of Chemical and Pharmaceutical Research, 2014, 6(7): Research Article

ESS Line Fitting

ENGI 4421 Joint Probability Distributions Page Joint Probability Distributions [Navidi sections 2.5 and 2.6; Devore sections

Multiple Choice Test. Chapter Adequacy of Models for Regression

ECE 595, Section 10 Numerical Simulations Lecture 19: FEM for Electronic Transport. Prof. Peter Bermel February 22, 2013

hp calculators HP 30S Statistics Averages and Standard Deviations Average and Standard Deviation Practice Finding Averages and Standard Deviations

Module 1 : The equation of continuity. Lecture 5: Conservation of Mass for each species. & Fick s Law

MEASURES OF DISPERSION

General Method for Calculating Chemical Equilibrium Composition

Linear Regression with One Regressor

Simple Linear Regression

Evaluation of uncertainty in measurements

STA 108 Applied Linear Models: Regression Analysis Spring Solution for Homework #1

Summary of the lecture in Biostatistics

STATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS. x, where. = y - ˆ " 1

Analysis of System Performance IN2072 Chapter 5 Analysis of Non Markov Systems

4. Standard Regression Model and Spatial Dependence Tests

Introduction to Matrices and Matrix Approach to Simple Linear Regression

Lecture 3. Sampling, sampling distributions, and parameter estimation

Can we take the Mysticism Out of the Pearson Coefficient of Linear Correlation?

Chemistry 163B Introduction to Multicomponent Systems and Partial Molar Quantities

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Selection of the Optimal Mathematical Model of Multiple Regression in the Ternary Mixture Experiments


12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM

Spreadsheet Problem Solving

Chemistry 163B Introduction to Multicomponent Systems and Partial Molar Quantities

Lecture Notes Types of economic variables

Econometric Methods. Review of Estimation

Analysis of Variance with Weibull Data

Applying the condition for equilibrium to this equilibrium, we get (1) n i i =, r G and 5 i

residual. (Note that usually in descriptions of regression analysis, upper-case

Simple Linear Regression

A New Family of Transformations for Lifetime Data

We have already referred to a certain reaction, which takes place at high temperature after rich combustion.

L5 Polynomial / Spline Curves

THE ROYAL STATISTICAL SOCIETY HIGHER CERTIFICATE

Statistics MINITAB - Lab 5

On modifications of the CSC GM code (R364n)

Chapter 5 Properties of a Random Sample

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

Numerical Simulations of the Complex Modied Korteweg-de Vries Equation. Thiab R. Taha. The University of Georgia. Abstract

Chapter 14 Logistic Regression Models

Journal of Mathematical Analysis and Applications

[ L] υ = (3) [ L] n. Q: What are the units of K in Eq. (3)? (Why is units placed in quotations.) What is the relationship to K in Eq. (1)?

2C09 Design for seismic and climate changes

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Generalization of the Dissimilarity Measure of Fuzzy Sets

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades

Computing Thermal Properties of Natural Gas by Utilizing AGA8 Equation of State

Comparison of Dual to Ratio-Cum-Product Estimators of Population Mean

IFYMB002 Mathematics Business Appendix C Formula Booklet

The equation is sometimes presented in form Y = a + b x. This is reasonable, but it s not the notation we use.

BAL-001-AB-0a Real Power Balancing Control Performance

Lesson 3. Group and individual indexes. Design and Data Analysis in Psychology I English group (A) School of Psychology Dpt. Experimental Psychology

Correlation and Regression Analysis

Introduction to local (nonparametric) density estimation. methods

Lecture 3 Probability review (cont d)

CHAPTER 2. = y ˆ β x (.1022) So we can write

CS 2750 Machine Learning. Lecture 8. Linear regression. CS 2750 Machine Learning. Linear regression. is a linear combination of input components x

Chapter Two. An Introduction to Regression ( )

Ahmed Elgamal. MDOF Systems & Modal Analysis

Midterm Exam 1, section 1 (Solution) Thursday, February hour, 15 minutes

Multi Objective Fuzzy Inventory Model with. Demand Dependent Unit Cost and Lead Time. Constraints A Karush Kuhn Tucker Conditions.

Modeling Thermal Conductivity of Concentrated and Mixed-Solvent Electrolyte Systems

The internal structure of natural numbers, one method for the definition of large prime numbers, and a factorization test

Rock-Paper-Scissors: An Example of a Multi-agent Dynamical System

Rademacher Complexity. Examples

Department of Mechanical Engineering ME 322 Mechanical Engineering Thermodynamics. Ideal Gas Mixtures. Lecture 31

Lecture 5: Interpolation. Polynomial interpolation Rational approximation

MMJ 1113 FINITE ELEMENT METHOD Introduction to PART I

DIFFERENTIAL GEOMETRIC APPROACH TO HAMILTONIAN MECHANICS

An Introduction to. Support Vector Machine

Quantitative analysis requires : sound knowledge of chemistry : possibility of interferences WHY do we need to use STATISTICS in Anal. Chem.?

Chapter 8. Inferences about More Than Two Population Central Values

Investigating Cellular Automata

The number of observed cases The number of parameters. ith case of the dichotomous dependent variable. the ith case of the jth parameter

Module 7: Probability and Statistics

Research on SVM Prediction Model Based on Chaos Theory

MULTIDIMENSIONAL HETEROGENEOUS VARIABLE PREDICTION BASED ON EXPERTS STATEMENTS. Gennadiy Lbov, Maxim Gerasimov

Wu-Hausman Test: But if X and ε are independent, βˆ. ECON 324 Page 1

Transcription:

A Helmholtz eergy equato of state for calculatg the thermodyamc propertes of flud mxtures Erc W. Lemmo, Reer Tller-Roth

Abstract New Approach based o hghly accurate EOS for the pure compoets combed at the reduced T, ρ of the mxture reducg parameters for the T,ρ : composto For smple mxtures wth relatvely smple fuctos very accurate represetato For odeal mxtures modfed reducg fuc.(t, ρ)+departure fuc. 0.% desty, % Cv, % bubble pot pressure

Abstract Two applcatos of mxture model cocepts developed depedetly USA ad Germay over the same tme perod cludg the developmet of dvdual equato for each bary system geeralzato of the model for a wde varety of mxtures Smlartes ad dffereces

Itroducto Dmesoless Helmholtz eergy, α RT Thermodyamc propertes of pure fluds wth hgh accuracy over a wde rage of T, P All thermodyamc propertes from a sgle mathematcal relato A large umber of adjustable coeff. (up to 58 terms reported) Referece equato for the calculato of property tables ad chart accuracy : restrcted oly by the accuracy of avalable measuremets A

Itroducto Due to the complexty Equato wth may coeffcets appled oly to pure fluds Hgh-accuracy EOS developed for oly about 30 pure substaces Oly a few attempts to descrbe mxture propertes usg multparameter EOS

Itroducto Plocer et al. (978) oe-flud theory appled to modfed Beedct-Webb- Rub(mBWR) EOS wth a mxg rule for the pseudo-crtcal temperature VLE ad Ethalpes at hgh pressures Platzer ad Maurer (993) geeralzed the 20-term Beder equato predcto for multcompoets mxture propertes Huber ad Ely (994) exteded correspodg states model based o mbwr EOS

Itroducto Mxture model (pure mxture) based o hgh-accuracy pure flud EOS for the pure compoets pressure-explct EOS smple tegrato of pressure geerally explct Helmholtz eergy adopted for the mxture combato rules for the Helmholtz eergy fucto of mxture compoets (ot mxg rules) ot requred a fxed structure of the pure flud equato sgle mathematcal expresso

Itroducto Geeral form of Helmholtz eergy model depedetly developed by both authors durg overlappg perods of tme Tller-Roth (993) focused o developg a accurate formulato for dvdual bary mxtures for a large amout of accurate expermetal data Lemmo (996) focused o developg a geeralzed model capable of accurate property calculato for a large umber of flud mxtures

Itroducto Purpose addtoal detal ad bacgroud for the models especally deal mxture dfferece betwee the models

Helmholtz eergy EOS for the pure fluds deal part + resdual part reducg parameter T, V : crtcal propertes ( ) ( ) ( ) δ + α δ α δ α,,, r 0 RT A ( ) ( ) [ ] + + + δ+ δ α 2 2 0 exp l l l, m a c c c ( ) ( ) + + δ δ + δ δ α 4 3 3 2 r exp, e m m c c V V T T ρ ρ δ,

Helmholtz eergy Trasformg to the Helmholtz eergy form by tegratg the relato p ρrt + δ α δ r

Mxture Model α of mxture : composto depedece α α 0 + α Ideal gas mxture : aalytcally from fuctos for the deal gas propertes of the pure fluds α 0 α 0 l 0 ( T, V, x) xα ( T, V ) + cosstet wth deal parts of the pure flud equato wth depedet reduced varable, δ r l x l x

Mxture Model Resdual Helmholtz eergy α r l r r (, δ, x ) α (, δ) + Δα (, δ, x) x lear combato of the pure flud resdual eergy departure fucto, Δα r teracto of dfferet speces the mxture Correspodg states prcple shape factors : ratos of the crtcal propertes T r T T h T T, r, r T, T r r T T

Mxture Model

Mxture Model Reducg parameter of the pure flud equato reducg fuctos of the mxture pseudo-crtcal parameter, δ composto depedece as x approaches ( x) T ( x ), δ( x) T T or ( x) T V ( x) ( x), δ( x) δ V V ( x),, V,

Lear mxture model Lear mxture model lear combatos of the pure flud parameters T ( x) xt,, V ( x) departure fucto Δα r omtted o adjustable parameters l l x V,

Lear mxture model : Applcato,,,2-tetrafluoroethae(R-34a) +dfluoromethae(r-32) close behavor of a deal soluto VLE : excellet results (Fg.a) Desty subcrtcal rage (Fg.2a) vapor, lqud desty : good agreemet Desty supercrtcal rage (Fg.2b) larger devato Isochorc heat capacty ±%(about the expermetal ucertaty)

Lear mxture model : Applcato water+ammoa mxture of polar compoets showg large mxg effects, wde two-phase rego VLE (Fg.b) dew curve : well hgh coc. the vapor phasedeal at lower pressures bubble curve : large devato Desty (Fg.2c) -% ~ +3% devato Excess ethalpy (Fg.3c) poor represetato : related to VLE

Lear mxture model : Applcato ethae+carbo doxde azeotropc behavor VLE (Fg.c) correct predcto Isochorc heat capacty (Fg.3b) less accurate >2% devato outsde of the crtcal rego

Lear mxture model : Applcato

Lear mxture model : Applcato

Lear mxture model : Applcato

Ehacemet to the Mxture Model Modfcato of the reducg fuctos T V ( ) x x x jt ( ), j T T + T j ( x) l l l l j x, j T,j, V,j : adjustable parameter x Effect of varyg T,j, V,j VLE of ethae+carbo doxde (Fg.4) ( T,j > V,j ) Cv :.4% devato (Fg.5) j V V, j, j T, j V, j 2 2, ( V + V ),, jj, jj

Ehacemet to the Mxture Model

Ehacemet to the Mxture Model

Ehacemet to the Mxture Model Modfcato by Tler-Roth supplemet the quadratc expresso wth β T, β V T V 2 2 βt ( x) x T, + x2t,2 + 2x2( x2 ) T,2 2 2 βv ( x) x V, + x2v,2 + 2x2( x2 ) V, 2 effect of expoets o the shape (Fg.6) fluece of T,2 or V,2 asymmetrc wth respect to composto Addtoal flexblty for tug of the mxture model

Ehacemet to the Mxture Model

Ehacemet to the Mxture Model Modfcato by Lemmo modfed the lear reducg fucto ζ j,ξ j,β j : modfy the shape of the reducg parameter β j, for ethae+co 2, ζ j -67.4K ( T,j 0.89) + + β ξ + ζ +,, ) ( ) ( l l j j j l l l j j j l x x V x V x x T x T j x x ( )( ) ( )( ),,,,,, + ξ + ζ j V j j j T j j V V T T

Departure fucto Departure fucto Δα r (, δ, x) x Δα (, δ,, ) x j j x x j oly to mxture propertes, o effect at the pure flud lmts magtude of departure fucto geerally oe order smaller tha the resdual α r modelg : regresso aalyss by Wager l l j + a set of of expermetal data lowest stadard devato ba of terms subset of term for Δα j optmzato : modfed Marquardt-Fletcher algorthm

Applcatos by Tller-Roth Itroducto a expoet γ the departure fucto for a bary system Δα r γ ( δ, x ) x ( x ) Δα (, δ, x, x ), 2 2 asymmetrc fluece very effectve for smple system, R-52+R34a reduce the umber of terms the departure fucto appled oly to bary systems to date For water+ammoa, xδ m exp ( e ) ( ) 2 m e δ ad x δ exp δ j j

Applcatos by Tller-Roth Results for R-34a+R-32 fve-term departure fucto Fg.7 ( cf. Fg.2b ) desty at lqud & vapor rego : ±0.% desty at crtcal rego : -6 ~ ±0.4% ( ) ( ) ( ) δ δ + δ δ α Δ γ 5 3 2 exp,, e m m r c c x x x

Applcatos by Tller-Roth

Applcatos by Tller-Roth Results for water+ammoa 4-term departure fucto : 9 adjustable parameters Δα x r (, δ, x) γ ( x ) c δ + c 4 x 2 m δ + 4 6 2 m 4 c δ exp m exp e4 ( δ ) 3 e m e ( δ ) + x c δ exp( δ ) 7 Fg.8 ( cf. Fg.b, Fg.3c ) mproved represetato

Applcatos by Tller-Roth

Applcatos by Lemmo Geeralzed departure fucto Δα Δα r j (, δ, x) x x F Δα (, δ) (, δ) c,, m : fxed the same for all mxture l j + F j, ζ j, ξ j, β j : adjustable parameters from (p,v,t,x) 0.2% desty, % heat capactes 0 c -2% bubble pot pressure at crtcal temp. δ l 5-0% bubble pot pressure at crtcal temp. father apart m j j j

Applcatos by Lemmo Results propae + -butae : Fg.0 less tha 0.04% average absolute devato desty 0.4% average absolute devato bubble pot pressure Ucertaty of the equato reported 0.2% desty % heat capactes Fg. : terary refrgerat mxture

Applcatos by Lemmo

Applcatos by Lemmo

Coclusos EOS for mxture terms of the Helmholtz eergy has bee developed ad appled to several flud systems by two dfferet research groups By Tller-Roth, model of sgle bary systems for whch comprehesve accurate data are avalable By Lemmo, geeralzed model was preseted whch accurately predcts thermodyamc propertes of mxtures Further research for more odeal systems ad complex multcompoet mxtures