Index of /book-cep/matlab/Likevekt_ved_min_gibbs_fra_Tore_Haug-Warberg

	Name	Last modified	Size	Description
	Parent Directory		-
	MinGibb.m	2010-07-14 06:43	1.5K
	README.html	2010-07-14 06:43	1.5K
	readme.txt	2010-07-14 06:43	1.5K

-----------------------------------
Matlab-program mingibbs: 
Kjemisk likevekt for ideell gass.
-----------------------------------

Likevekten bestemmes med å mimimalisere Gibbs energi for blandingen 
under bibetingelsen om at støkiometrien er oppfylt. 
Skrevet av Tore Haug-Warberg, juli 2004

Programmet trenger følgende input:
	gf (kalt m0  programmet): Gibbs energi for hver komponent ved T og 1 bar [J/mol]
	A: atom-matrise som beskriver støkiometrien
	n0: fødesammensetning [mol]
	T: temperatur [K]
	p: trykk [N/m2]
Output'en er likevekts-sammensetningen n.


Eksempel på bruk:

% Komponenter: CH3COOH, CH3OH, CO, CO2, H2O, H2
% Føde er 50 mol metanol og 50 mol CO2 

gf0 = [-375 -163 -137 -394 -229 0]'*1.e3; % Gf for komponentene ved 298K
f0 = [-433 -201 -111 -394 -242 0]'*1.e3; % Hf for komponentene ved 298K
cp0 = [67 44 29 37 34 29]';
sf0 = -(gf0-hf0)/298.15; % Sf for komponentene
T = 200+273.15; % temperatur [K]
hf = hf0 + cp0*(T-298.15);  % antar her konstant cp
sf = sf0 + cp0*log(T/298.15);
gf  = hf-T*sf % 
A  = [2 1 1 1 0 0; 4 4 0 0 2 2; 2 1 1 2 1 0]; % atommatriser for C, H, O
n0 = [0 50 0 50 0 0]'; % føde
p = 100e5; %  = 100 bar

n=mingibb(n0,T,p,gf,A)

----
n =

   36.6734 eddiksyre
    0.0087 metanol
    0.0203 CO
   26.6242 CO2
   23.3758 H2O
    3.2599 H2

Kommentar:
Ønsket totalreaksjon er: 4 CH3OH + 2 CO2 = 3 CH3COOH  + 2 H2O 
dvs. vi har overskudd CO2 i føden.
Ved likevekt vil nesten all metanol omsettes og det dannes lite H2 og omtrent intet CO.