-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathequations_just_list_internal_format.txt
118 lines (118 loc) · 4.34 KB
/
equations_just_list_internal_format.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
1 :: one = Instantiate ( value , value )
2 :: zero = Instantiate ( value , value )
3 :: onehalf = Instantiate ( value , value )
4 :: to = Instantiate ( t , value )
5 :: te = Instantiate ( t , value )
6 :: p = neg ( ParDiff ( U , V ) )
7 :: T = ParDiff ( U , S )
8 :: chemPot = ParDiff ( U , n )
9 :: H = U - p . V
10 :: A = U - T . S
11 :: G = U + p . V - T . S
14 :: Ue = inv ( charge ) . U
15 :: v_x = ParDiff ( r_x , t )
16 :: v_y = ParDiff ( r_y , t )
17 :: v_z = ParDiff ( r_z , t )
20 :: _Mm = Mm
21 :: Cp = ParDiff ( H , T )
22 :: Cv = ParDiff ( U , T )
23 :: kq_x = inv ( V ) . ParDiff ( U , T ) . v_x
24 :: kq_y = inv ( V ) . ParDiff ( U , T ) . v_y
25 :: kq_z = inv ( V ) . ParDiff ( U , T ) . v_z
26 :: kq = Stack ( kq_x , kq_y , kq_z )
27 :: kc_x = ( Mm | S in N & S | inv ( chemPot ) ) . inv ( V ) . ParDiff ( U , p ) . v_x
28 :: kc_y = ( Mm | S in N & S | inv ( chemPot ) ) . inv ( V ) . ParDiff ( U , p ) . v_y
29 :: kc_z = ( Mm | S in N & S | inv ( chemPot ) ) . inv ( V ) . ParDiff ( U , p ) . v_z
30 :: kc = Stack ( kc_x , kc_y , kc_z )
31 :: kd_x = inv ( chemPot ) . ( v_x : ( inv ( V ) : ParDiff ( U , chemPot ) ) )
32 :: kd_y = inv ( chemPot ) . ( v_y : ( inv ( V ) : ParDiff ( U , chemPot ) ) )
33 :: kd_z = inv ( chemPot ) . ( v_z : ( inv ( V ) : ParDiff ( U , chemPot ) ) )
34 :: kd = Stack ( kd_x , kd_y , kd_z )
35 :: h = H : inv ( n )
36 :: m = _Mm | S in N & S | n
37 :: _m = m
38 :: density = _m . inv ( V )
39 :: T_NK = P_N_NK | N | T
41 :: Ea = Instantiate ( GasConstant . T_NK , value )
42 :: K_NK = Ko : exp ( neg ( Ea ) . inv ( GasConstant . T_NK ) )
44 :: c = inv ( V ) : n
45 :: _c = c
48 :: Ayz = r_y . r_z
49 :: Axz = r_x . r_z
50 :: Axy = r_x . r_y
51 :: _density = density
52 :: _h = h
53 :: _kq_x = kq_x
54 :: _Cv = Cv
55 :: _kq_y = kq_y
56 :: _kq_z = kq_z
57 :: _kq = kq
58 :: _kc_x = kc_x
59 :: _Cp = Cp
60 :: _kc_y = kc_y
61 :: _kc_z = kc_z
62 :: _kc = kc
63 :: _kd_x = kd_x
64 :: _kd_y = kd_y
65 :: _kd_z = kd_z
66 :: _kd = kd
67 :: fV = inv ( _density ) . _kc_x . Ayz . D | N | p
68 :: fnd_AS = Ayz : neg ( _kd_x ) . D_NS_AS | N & S | chemPot
69 :: fnd = Fd | A & S | fnd_AS
70 :: fHd_A = ( F_NS_AS | N & S | _h ) | S in A & S | fnd_AS
71 :: fHd = F_N_A | A | fHd_A
72 :: d = sign ( F_N_A | N | p )
73 :: c_AS = ( onehalf . ( F_NS_AS - d : abs ( F_NS_AS ) ) ) | N & S | c
74 :: fnc_AS = fV : c_AS
75 :: fnc = F_NS_AS | A & S | fnc_AS
76 :: dndt = fnc + fnd + _nProd
77 :: fHc_A = ( F_NS_AS | N & S | _h ) | S in A & S | fnc_AS
78 :: fHc = F_N_A | A | fHc_A
79 :: fw_A = Instantiate ( fHc_A , value )
80 :: fw = F_N_A | A | fw_A
81 :: fq_A_x = Ayz . _kq_x . D | N | T
82 :: fq = F_N_A | A | fq_A_x
83 :: dHdt = fHc + fHd + fq + fw
84 :: Ho = Instantiate ( H , value )
85 :: no = Instantiate ( n , value )
86 :: n = Integral ( dndt :: t in [ to , te ] ) + no
87 :: H = Integral ( dHdt :: t in [ to , te ] )
88 :: additive = Instantiate ( additive , value )
89 :: i = TotalDiff ( charge , t )
90 :: _i = i
91 :: elResist = inv ( _i ) . Ue
92 :: elResist = Instantiate ( elResist , value )
93 :: elConductC = inv ( elResist ) . additive
94 :: _elConductC = elConductC
95 :: Ue = inv ( _elConductC ) . i
96 :: dUedt = one . Ue
97 :: dUedt = Root ( Ue )
98 :: dUedt = Instantiate ( dUedt , zero )
99 :: _i = i
104 :: _additive = additive
105 :: _T = T
106 :: intensities = MixedStack ( p , T , chemPot , c , Ue )
107 :: nTotal = one_S | S in N & S | n
108 :: e = meas - setPoint
109 :: xo = Instantiate ( x , value )
110 :: dxdt = A | N | x + B | A | e
111 :: x = Integral ( I_N_D | D | dxdt :: t in [ to , te ] )
112 :: _additive = additive
113 :: Imeasured = _i
114 :: _Ue = Ue
115 :: UeMeasured = _Ue
116 :: addMeasured = _additive
117 :: RComputed = inv ( Imeasured ) . UeMeasured
118 :: store = MixedStack ( Imeasured , UeMeasured , RComputed , addMeasured )
119 :: setPoint = Instantiate ( setPoint , value )
123 :: c = P_NK . ( P_NS_KS | N & S | _c )
124 :: co = Instantiate ( c , value )
125 :: x = inv ( co ) . c
126 :: y = C | N | x + D . e
127 :: GasConstant = Avogadro . Boltzmann
128 :: N_NK_KS = P_K_NK | K | N_K_KS
129 :: phi = Product ( x ^ ( N_NK_KS ) , K & S )
130 :: nProd = V | N | ( P_N_NK | N & K | ( ( K_NK . phi ) . ( P_NS_KS | K & S | N_NK_KS ) ) )
131 :: _nProd = nProd
132 :: fnd_AS = neg ( DiffCoeff ) . ( Dd_NS_AS | N & S | c )
133 :: V = m . inv ( _density )