ケミカルエンジニアのためのExcelを用いた化学工学計算法

Size: px

Start display at page:

Download "ケミカルエンジニアのためのExcelを用いた化学工学計算法"

れいがみおか
9 years ago
Views:

1 VBA 7.1 f ()= ( f ( )) y = f ()(1) y = f ( )( ) + f ( ) (1) = f ( ) f ( ) (2) 1 n = = y = f() y = () 1 n+1 = n (f( n )f( n ))

2 log()2 145

3 7.2 f ( ) f ( ) (3) (4) ) ( ) ( ) '( = f f f ) ' 1 f = + ) ( f ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) '( ) = = = f f f f f f f f f ( 146

4 X 1 X1 X2: 1 FXf( 1 ) FX1f( ) NX1: DX: 147

5 Sub () ' Dim X As Double, X1 As Double, X2 As Double 'X,X1 Dim N1 As Integer ' N1 EPS =.1 ' X = Cells(3, 3) 'Ecel X X1 = Cells(4, 3) 'Ecel X1 1: FX = X ^ 3-2 * X ^ 2 - X + 2 ' FX FX1 = X1 ^ 3-2 * X1 ^ 2 - X1 + 2 ' FX1 X2 = (X * FX1 - X1 * FX) / (FX1 - FX) ' X2 DX = Abs((X2 - X1) / X1) 'X1 X2 DX N1 = N1 + 1 ' N1 ' Cells(12 + N1, 2) = N1 ' N1 Cells(12 + N1, 3) = X 'X Cells(12 + N1, 4) = X1 'X1 Cells(12 + N1, 5) = X2 ' X2 Cells(12 + N1, 6) = FX ' FX Cells(12 + N1, 7) = FX1 ' FX1 '********************************** ' DX EPS If DX < EPS Then GoTo 2 ' X1 X Else X = X1 ' X1 X X1 = X2 ' X2 X1 GoTo 1 End If '*************************************************** 2: Cells(8, 3) = X1 ' X1 Cells(9, 3) = N1 ' N1 End Sub

'X1 X2 DX N1 = N1 + 1 ' N1 '-------------------- Cells(12 + N1, 2) = N1 ' N1 Cells(12 + N1, 3) = X 'X Cells(12 + N1, 4) = X1 'X1 Cells(12 + N1, 5) = X2 ' X2 Cells(12 + N1, 6) = FX ' FX Cells(12 + N1,

6 f()= f( )f( 1 )f( )f( 1 )< 2 1 y=f()[, 1 ] [, 1 ] 1 1 [, 1 ] c1 c1 ( + 1 )/2 2f( c1 )f( )f( c1 )f( )> = c1 f( c1 ) f( )f( c1 )f( )< 1 = c1 31,2 = 1-1 /.1 f( ) 1 f( c2 ) f( c3 ) f( c1 ) c = c1 c3 c f( c )f( )>c f( c )f( )<c1 149

7 Sub () Dim tin As Range Dim tout As Range Set tin = Range("B1:B1") Set tout = Range("C1:C1") ' = tin.cells(1, 1) 1 = tin.cells(2, 1) f = ^ 3-2 * ^ f1 = 1 ^ 3-2 * 1 ^ If f * f1 > Then MsgBo "" Eit Sub End If Do c = ( + 1) / 2 fc = c ^ 3-2 * c ^ 2 - c + 2 y = f * fc If y = Then Eit Do ElseIf y > Then = c Else 1 = c End If n = n + 1 If n = 5 Then c 1 MsgBo "5" Eit Sub End If Loop Until Abs((1 - ) / 1) <.1 tout.cells(1, 1) = "=" & c tout.cells(2, 1) = "=" & n End Sub 15

cells(2, 1) f = ^ 3-2 * ^ 2 - + 2 f1 = 1 ^ 3-2 * 1 ^ 2-1 + 2 If f * f1 > Then MsgBo "" Eit Sub End If Do c = ( + 1) / 2 fc =

8 7.4 i y i y i f ( i ) f ( i ) e i = y i f ( i ) y=a 1 +a S y e i S y=a 1 +a a a 1 y=a 1 +a S=e i ={y i -(a 1 i +a )} S (S a )= (S a 1 )= a,a S yf () = a + a 1 a a 1 X Y

i +a )} S (S a )= (S a 1 )= a,a S 1 5 2 -yf () = a + a

9 y=a 1 +a y e i Ecel 152

10 -y -y X X Y Y 153

11 X/Y/ (S) (F) 1 R 154

12 L E OK 155

13 3 Sub SAISYO1() Dim NP As Integer Dim XX(1), YY(1) As Double Dim targetin, targetout As Range NP = Cells(2, 2) For i = 1 To NP XX(i) = Cells(3 + i, 2) YY(i) = Cells(3 + i, 3) Net i SX2 = SX = SSX = SY = SXSY = For i = 1 To NP SX2 = SX2 + XX(i) ^ 2 SX = SX + XX(i) 156

14 SY = SY + YY(i) SXSY = SXSY + XX(i) * YY(i) Net i SSX = SX ^ 2 A = (SX2 * SY - SX * SXSY) / (NP * SX2 - SSX) A1 = (-SX * SY + SXSY * NP) / (NP * SX2 - SSX) Cells(2, 5) = "Y=A+A1*X" Cells(3, 5) = "A=" Cells(4, 5) = "A1=" Cells(3, 6) = A Cells(4, 6) = A1 End Sub 157

16 1 12wt.% 25wt.% 1g F[g]V[g]L[g] w 1 w 2 y 1 y F=V+L (1) w 1 Fy 1 V 1 L (2) w 2 Fy 2 V 2 L (3) y 1 w 1.12w y 1 y F1 (1)(2) V+L1 (1).121V+.25L.25LV12 (2) (1) (2) 159

17 1 1 V 1 =.25 L 12 (4) (5) (5) V 52 [g], L48 [g] wt. SO 3 (SO 3 22wt.%78wt.%) 1g/sSO 3 97wt.% H 2 SO wt.%H 2 SO 4 H 2 SO 4 95wt.% (98.5wt.% H 2 SO 4 ) H 2 SO 4 98SO 3 8S32 (a) 95wt.% H 2 SO 4 (b), y, z, w [g/s] i [g/s] A 16

%) 1g/sSO 3 97wt.% H 2 SO 4 98.5wt.%H 2 SO 4 H 2 SO 4 95wt.

18 A 1 + = i + y (1) (.22)( 1) =.985y (2) (.78)(1)=i (3) B + w = z (4) z (5) (1)(5) A B y z = 1 i w AB AB A.312 B y z = i w = 151, y= 173 z= 353., i= 78 w= 21 a) 95wtH 2 SO 4 151g/sb 98.5wt H SO 4 173g/s 161

8 1 1 1 z = 1 i 78.95.97.985 w AB AB A.312 B.3216 1.312 1.95 1.3216 1.97 1 1 1 y 8.

19 7.6 y = f() [a, b]n, 1, 2, n y, y 1, y 2, y n P, P 1, P 2, P n n S = S 1 + S 2 + S 3 + S n (1) h=1/n S 1 =(h/2)(y +y 1 ), S 2 =(h/2)(y 1 +y 2 ),, S n =(h/2)(y n1 +y n ) (1) S=(h/2)(y +y 1 )+(h/2)(y 1 +y 2 )+(h/2)(y 2 +y 3 )++(h/2)(y n 1+y n ) =(h/2)y + y n + 2(y 1 + y 2 + y n 1) y y = f() P P 1 P 2 P -1 P Pn S 1 S 2 S Sn y y 1 y 2 y 1 y y n a 1 2 n = y 1 b = S 1 1 h y 162

20 Sub fore1() Dim a, b, h, V As Single a = Cells(6, 3) b = Cells(6, 4) h = Cells(6, 6) M = (b- a) / h SS = N = = For N = 1 To M - 1 X = a + N * h SS = SS + 3 * X ^ 2 Net N S =.5 * h * (3 * a ^ * b ^ * SS) Cells(12, 3) = S End Sub y=3 2 =~2 163

21 . 8.1 adsorption 4mg/ 4m 3 6g mg/g/g nw (1) + V[] C [mg/] : W[g] (C -C)V (2) nw(c -C)V (3) C C V n = V = ( C C ) (4) W W V (3) 1 W 1g : n[mol/g] V n g / g) = ( C C W 2 C 4 3 ( 4 ) = ( C 4) 6 164

22 -2/3 Langmuir C n KC n = 1+ KC 1 165

23 g 4mg/ 4m 3 C[mg/] [g/g] 3g g 4mg/ 4m 3 C[mg/] [g/g] 15g

24 8.2 (a) Henry nkc (1) K 1 (a) C n C p (a) (b) (c) (d)bet 1 167

25 (b)langmuir 2 1 [-]r [mols -1 ]a[mols -1 ] 2 Langmuir ra (2) r[mols -1 ](1-)[-]C[mol m -3 ] r b(1-)c (3) b[m 3 s -1 ] ab(1-)c (4) b C θ = (5) a + b C 168

26 n n [molg -1 ] b/a K[m 3 mol -1 ](K)θ = n n KC n = 1+ KC (6) (c)(freundlich) 1(c) nc 1/ (7) n[molg -1 ][-] (d)bet(brunauer-emmett-teller) BET 1(d) 1 BET 3 q = q K c ( 1 c) ( 1 + Kc c) 169

27 3 8.3 Langmuir 12 1 Kn n K 1 n, C, C/n C [molm -3 ] n [molg -1 ] C/n [gm -3 ] (6) 1 K C = 1+ n n KC (9) C C n (1 + KC ) C n KC 1 1 = n K n = C = C (1) 17

28 =1n =1(n K)(1)C/nC =1n =1n K 1 C/nC 1 4 =1/n =.36 =1/(n K)=35.7 n K n =3.27 molg -1 K =.86 m 3 mol -1 n KCn(Cn) 5 = 1 / n 1 n K 4 171

29 n [ molg -1 ] Langmuir n : C [ molm -3 ] 5 12 C n V[m 3 ] C [molm -3 ]W[g] n[molg -1 ]WnW (C -C b )VC b [molm -3 ] nw(c -C b )V (11) n=-v/w(c b -C ) (12) 1 172

30 (12)-V/W C C 1 C 2 C 3 6 n [ molg -1 ] q 1 V W c 1 c C [ molm -3 ] 7 173

31 -V/W-V/W=-1/1=-.1 m 3 g -1 2 mol m -3 (4-2)1=38 mol q 1 V? - W c 1 c

コンピュータ概論

コンピュータ概論 4.1 For Check Point 1. For 2. 4.1.1 For (For) For = To Step (Next) 4.1.1 Next 4.1.1 4.1.2 1 i 10 For Next Cells(i,1) Cells(1, 1) Cells(2, 1) Cells(10, 1) 4.1.2 50 1. 2 1 10 3. 0 360 10 sin() 4.1.2 For