Bird E 1.4-2

bird_e_1_4_2.py

import numpy as np
komponente=np.array(['CO2','O2','N2'])
mm_i=np.array([44.01,32,28.02])
mu=np.array([1462e-7,2031e-7,1754e-7])
x_a=np.array([0.133,0.039,0.828])
phi_ab = np.zeros([komponente.size, komponente.size])
for alpha in range(phi_ab.shape[0]):
    for beta in range(phi_ab.shape[1]):
        phi_ab[alpha, beta] = 1/np.sqrt(8)*(
            1+mm_i[alpha]/mm_i[beta])**(-1/2.)*(
            1+(mu[alpha]/mu[beta])**(1/2.)*(
                mm_i[beta]/mm_i[alpha]
            )**(1/4.)
        )**2
print('phi_alpha_beta')
print(phi_ab)
print('sum(x_alpha phi_alpha_beta)')
print(phi_ab.dot(x_a))
print('Viskosität, g/cm/s')
print('{:g}'.format(
    sum(x_a * mu / phi_ab.dot(x_a))*1e7
)+'(10)^-7 g/cm/s')