{UTF-8}
********************************************************
	.no2-n2o4-lechatelier
********************************************************~
		Simulation: Peter Bützer
	|

abführen 1=
	IF THEN ELSE(NO2<=0.01, 0, ka1 )
	~	mol/Microsecond [0,0.01]
	~	konstante Geschwindigkeit, mit der NO2 zu- oder abgeführt wird
	|

abführen 2=
	IF THEN ELSE(N2O4<=0.01, 0 , ka2 )
	~	mol/Microsecond [0.01,0.01]
	~	konstante Geschwindigkeit, mit der N2O4 zu- oder abgeführt wird
	|

NO2= INTEG (
	RückR-abführen 1-HinR,
		0.1)
	~	mol [0,?]
	~		|

Druck=
	1
	~	Dmnl [1,5]
	~	Der Druck wird bei den Mengen als Faktor gerechnet. Begründung: Die \
		Konzentration pro Volumen steigt mit dem Druck. Da das Modell von einem \
		konstanten Volumen ausgeht, müssen die Mengen erhöht werden, um eine \
		höhere Konzentration zu simulieren: c=n/V
	|

ka2=
	0.001
	~	mol/Microsecond [-0.01,0.01]
	~	Geschwindigkeits-Konstante, mit der N2O4 zu- oder abgeführt wird
	|

HinR=
	kf*Q10*Druck*NO2*Druck*NO2
	~	mol/Microsecond [0,?]
	~	Reaktionsgeschwindgkeit
	|

RückR=
	kb*Druck*N2O4
	~	mol/Microsecond [0,?]
	~	Reaktionsgeschwindgkeit
	|

ka1=
	0
	~	mol/Microsecond [-0.01,0.01]
	~	Geschwindigkeits-Konstante, mit der NO2 zu- oder abgeführt wird.
	|

N2O4= INTEG (
	HinR-abführen 2-RückR,
		0.1)
	~	mol [0,?]
	~		|

K=
	N2O4/(NO2*NO2)
	~	1/mol [0,?]
	~		|

kb=
	0.7/21
	~	1/Microsecond [0,?]
	~	Reaktionsgeschwindgkeits-Konstante
	|

kf=
	7
	~	1/Microsecond/mol [0,?]
	~	Reaktionsgeschwindgkeits-Konstante, HWZ der Reaktion ca. 0.1 Mikrosekunde
	|

Basis=
	2.5
	~	Dmnl [1.5,3]
	~	Basis der RGT-Regel, Faktor, um den die Geschwindigkeit bei einer \
		Erhöhung um 10 K erhöht werden muss.
	|

Q10=
	1*Basis^((Referenz Temperatur-Reaktions Temperatur)/T10)
	~	Dmnl [0,?]
	~	RGT-Regel, bei der endothermen Reaktion muss (Temp Kelvin-Referenz \
		Temperatur) gerechnet werden, also der negative Wert von einer exothermen \
		Reaktion.
	|

Referenz Temperatur=
	283
	~	K [0,?]
	~	Auf diese Temperatur ist die Gleichgewichtskonstante bezogen
	|

Reaktions Temperatur=
	298
	~	K [263,373]
	~	Temperatur der Reaktionsmischung
	|

T10=
	10*1
	~	K [10,10]
	~	10 K, Bei einer Erhöhung der Temperatur wird der Faktor (Quotient) \
		berechnet
	|

********************************************************
	.Control
********************************************************~
		Simulation Control Parameters
	|

FINAL TIME  = 120
	~	Microsecond
	~	The final time for the simulation.
	|

INITIAL TIME  = 0
	~	Microsecond
	~	The initial time for the simulation.
	|

SAVEPER  = 
        TIME STEP
	~	Microsecond [0,?]
	~	The frequency with which output is stored.
	|

TIME STEP  = 0.01
	~	Microsecond [0,?]
	~	The time step for the simulation.
	|

\\\---/// Sketch information - do not modify anything except names
V300  Do not put anything below this section - it will be ignored
*View 1
$192-192-192,0,Times New Roman|12||0-0-0|0-0-0|0-0-255|-1--1--1|-1--1--1|120,120,100,0
10,1,NO2,172,245,37,20,3,3,0,4,0,0,0,0,0-0-0,192-255-192,|12||0-0-0
10,2,N2O4,372,247,38,20,3,3,0,4,0,0,0,0,0-0-0,192-255-255,|12||0-0-0
1,3,5,2,4,0,0,22,0,0,0,-1--1--1,,1|(298,232)|
1,4,5,1,100,0,0,22,0,0,0,-1--1--1,,1|(230,232)|
11,5,252,257,232,5,6,34,3,0,0,1,0,0,0
10,6,HinR,257,246,23,12,40,131,0,0,-1,0,0,0
1,7,9,1,4,0,0,22,0,0,0,-1--1--1,,1|(252,258)|
1,8,9,2,100,0,0,22,0,0,0,-1--1--1,,1|(320,258)|
11,9,300,301,258,5,6,34,3,0,0,1,0,0,0
10,10,RückR,301,276,39,12,40,131,0,0,-1,0,0,0
12,11,0,191,274,30,15,8,0,0,0,3,0,0,0
1,12,1,5,1,0,0,0,0,64,0,-1--1--1,,1|(213,215)|
1,13,2,9,1,0,0,0,0,64,0,-1--1--1,,1|(352,286)|
10,14,kb,311,357,10,11,8,3,0,0,0,0,0,0
10,15,kf,185,184,14,18,8,3,0,0,0,0,0,0
1,16,15,5,1,0,0,0,0,64,0,-1--1--1,,1|(225,183)|
1,17,14,9,1,0,0,0,0,64,0,-1--1--1,,1|(293,312)|
12,18,0,904,245,315,232,3,188,0,0,1,0,0,0
Gleichgewicht
10,19,Reaktions Temperatur,59,161,53,25,8,131,0,0,0,0,0,0
10,20,Referenz Temperatur,90,63,48,22,8,3,0,0,0,0,0,0
10,21,Basis,295,81,24,12,8,3,0,0,0,0,0,0
10,22,Q10,226,149,28,12,8,3,0,0,0,0,0,0
1,23,21,22,1,0,0,0,0,64,0,-1--1--1,,1|(287,119)|
1,24,20,22,1,0,0,0,0,64,0,-1--1--1,,1|(182,91)|
1,25,19,22,1,0,0,0,0,64,0,-1--1--1,,1|(161,136)|
10,26,T10,343,127,19,12,8,3,0,0,0,0,0,0
1,27,26,22,1,0,0,0,0,64,0,-1--1--1,,1|(302,150)|
1,28,22,6,1,0,0,0,0,64,0,-1--1--1,,1|(261,200)|
10,29,K,86,338,9,12,8,3,0,0,0,0,0,0
1,30,1,29,1,0,0,0,0,64,0,-1--1--1,,1|(138,311)|
1,31,2,29,1,0,0,0,0,64,0,-1--1--1,,1|(186,342)|
10,32,Druck,56,246,26,12,8,3,0,0,0,0,0,0
1,33,32,5,1,0,0,0,0,64,0,-1--1--1,,1|(137,195)|
1,34,32,10,1,0,0,0,0,64,0,-1--1--1,,1|(155,295)|
12,35,48,371,382,10,8,0,3,0,0,-1,0,0,0
1,36,38,35,4,0,0,22,0,0,0,-1--1--1,,1|(371,350)|
1,37,38,2,100,0,0,22,0,0,0,-1--1--1,,1|(371,290)|
11,38,48,371,320,8,6,33,3,0,0,4,0,0,0
10,39,abführen 2,423,320,44,12,40,3,0,0,-1,0,0,0
12,40,48,171,384,10,8,0,3,0,0,-1,0,0,0
1,41,43,40,4,0,0,22,0,0,0,-1--1--1,,1|(171,351)|
1,42,43,1,100,0,0,22,0,0,0,-1--1--1,,1|(171,289)|
11,43,48,171,320,8,6,33,3,0,0,4,0,0,0
10,44,abführen 1,221,320,42,12,40,131,0,0,-1,0,0,0
1,45,2,38,1,0,0,0,0,64,0,-1--1--1,,1|(399,288)|
1,46,1,43,1,0,0,0,0,64,0,-1--1--1,,1|(195,288)|
10,47,ka2,446,374,16,12,8,3,0,0,0,0,0,0
10,48,ka1,242,375,16,12,8,3,0,0,0,0,0,0
1,49,48,43,1,0,0,0,0,64,0,-1--1--1,,1|(204,363)|
1,50,47,38,1,0,0,0,0,64,0,-1--1--1,,1|(401,362)|
///---\\\
:GRAPH Gleichgewicht
:TITLE Gleichgewicht
:X-AXIS Time
:SCALE
:VAR NO2
:UNITS mol
:LINE-WIDTH 2
:VAR N2O4
:UNITS mol
:LINE-WIDTH 2
:SCALE
:VAR K
:LINE-WIDTH 2
:L<%^E!@
1:Current.vdf
9:Current
23:0
15:0,0,0,1,0,0
19:100,0
27:2,
34:0,
4:Time
5:abführen 1
24:0
25:119.997
26:119.997