------------------------------------------------------------------------------------- python,kinetics : -('diff(C,x,1))*u-C*k+('diff(C,x,2))*D=0,numerical,part 3,version 1 peter.vlasschaert@gmail.com,28/01/2023 ------------------------------------------------------------------------------------- mathematical model : reactor ----------------------------- ode :-('diff(C,x,1))*u-C*k+('diff(C,x,2))*D=0 bc1 :at(x=0),Cin = u*C -D*'diff(C,x) bc2 :at(x=L),'diff(C,x) =0 ----------------------------- ╒═════════╤════════════╤═══════════════════╤═══════════════════╕ │ Index │ Distance │ Concentration A │ Concentration B │ ╞═════════╪════════════╪═══════════════════╪═══════════════════╡ │ 0 │ 0 │ 1 │ -0 │ ├─────────┼────────────┼───────────────────┼───────────────────┤ │ 1 │ 1.42857 │ 0.315203 │ 0.684797 │ ├─────────┼────────────┼───────────────────┼───────────────────┤ │ 2 │ 2.85714 │ 0.0993541 │ 0.900646 │ ├─────────┼────────────┼───────────────────┼───────────────────┤ │ 3 │ 4.28571 │ 0.0313134 │ 0.968687 │ ├─────────┼────────────┼───────────────────┼───────────────────┤ │ 4 │ 5.71429 │ 0.00988074 │ 0.990119 │ ├─────────┼────────────┼───────────────────┼───────────────────┤ │ 5 │ 7.14286 │ 0.00308066 │ 0.996919 │ ├─────────┼────────────┼───────────────────┼───────────────────┤ │ 6 │ 8.57143 │ 0.00107819 │ 0.998922 │ ├─────────┼────────────┼───────────────────┼───────────────────┤ │ 7 │ 10 │ 0 │ 1 │ ╘═════════╧════════════╧═══════════════════╧═══════════════════╛ ------------------------------------------------------------------------- distance concentration A concentration B 0 0.000000 1.000000 -0.000000 1 1.428571 0.315203 0.684797 2 2.857143 0.099354 0.900646 3 4.285714 0.031313 0.968687 4 5.714286 0.009881 0.990119 5 7.142857 0.003081 0.996919 6 8.571429 0.001078 0.998922 7 10.000000 0.000000 1.000000