Nonisothermal model added.

__pycache__/

 No newline at end of file

__pycache__/

 No newline at end of file

#!/usr/bin/env python3

# -*- coding: utf-8 -*-

"""

@author: Stefan Hiemer

"""

import numpy as np

import matplotlib.pyplot as plt

from forces import F_newt_isoth,F_newt_vft,F_newt_eyr

if __name__=="__main__":

    km = 0.163 # m s /K

    Th = 473 # K

    Tm = 424 # K

    rho_s = 1250 # kg m**(-3)

    rho_m = 1000 # kg m**(-3)

    h = 93600 # J/mol

    Cs = 1800 # J/ mol K 

    T0 = 293 # K

    alpha=56.31 # degree

    Ro= 0.00197/2 # m

    Ri= 0.0002 # m

    Lc= 0.001 # m 

    Q = 80.574 # J/mol

    R = 8.314 # J /(Kmol)

    eps = Q/R # K

    Tv = 0 # K

    eta0 = 3.546 * 10**3 # Pas 

    etac = eta0 * np.exp(Q/(R*Th))

    #

    fig,ax = plt.subplots(1,1)

    # create grid of velocities

    Us = np.linspace(0,0.006,101)[1:]

    # calculate isothermal

    F = F_newt_isoth(km=km, Th=Th, Tm=Tm, rho_s=rho_s, rho_m=rho_m, h=h,

                     Cs=Cs, T0=T0, Us=Us, eta0=eta0*np.exp(Q/(R*(Th+Tm)/2)), 

                     alpha=alpha, Ro=Ro, Ri=Ri, etac=etac, Lc=Lc,

                     debug=False)

    # plot isoth. plot

    ax.plot(Us*10**3,F,label="isoth.")

    # arrhenius

    F = F_newt_eyr(km=km, Th=Th, Tm=Tm, rho_s=rho_s, rho_m=rho_m, h=h, Cs=Cs, 

                   T0=T0, Us=Us, eta0=eta0, Q=Q, R=R, alpha=alpha, Ro=Ro, Ri=Ri, 

                   etac=etac, Lc=Lc, debug=False)

    # plot arrh.

    ax.plot(Us*10**3,F,label="arrh.")

    # vogel fulcher tamann 

    F = F_newt_vft(km, Th, Tm, rho_s, rho_m, h, Cs, T0,

                   Us, eta0, eps, Tv, alpha, Ro, Ri, etac, Lc,

                   debug=False)

    # plot vft.

    ax.plot(Us*10**3,F,label="vft.")

    ax.legend()

    ax.set_xlabel(r"$U_{s} (m/s)$")

    ax.set_ylabel(r"$U_{F} (N)$")

    ax.set_yscale("log")

    plt.show()

 No newline at end of file