Python

ERI
Cinética Reactores Químicos

Cinética

Reactores Ideais

Como Instalar

Exemplos


  • Python Computations in Science and Engineering

    • Exemplos Precedência

    Gráficos

    from pylab import *
    x=[1,2,3]
    y=[2.3,3.5,10.0]
    plot(x,y)
    plot(x,y,'o--')
    xlabel('x')
    ylabel('y')
    xlim(0,4)
    ylim(0,15)
    title('$F(x)=\\frac{x^2}{y_i}$') show()
    close('all')
    x=arange(0,10,0.2)
    plot(x,sin(x),'gd:',label='$sin(x)$')
    plot(x,cos(x),'rq:',label='cos(x)')
    legend(loc=0)
    show()
    grid('on')
  • Voltar

    • Matemática Simbólica

    from sympy import *
    x=Symbol('x')
    integrate(1/x,x)    		 log(x)
    integrate(1/(1-x),(x,0,x) -log(-x + 1)
    integrate(1/(1-x),(x,0,0.9) 2.30258509299405
    solve(x**2-2*x-1,x) [1 + sqrt(2), -sqrt(2) + 1]
    diff(x**2+2*x+4) 2*x + 2 
    # -*- coding: utf-8 -*-
from __future__ import division
from sympy import *
from pylab import *
from sympy.solvers import solve
close('all')
##################################################################
# Integração analitica de de Equações Diferenciais
from fractions import Fraction
C=Symbol('C')
k=Symbol('k')
Co=Symbol('Co')
t=Symbol('t')
n=[-1,0,1/2,1,1.5,2]
figure(facecolor='white')
ax1=axes(frameon=False)
ax1.get_xaxis().tick_bottom()
ax1.axes.get_yaxis().set_visible(False)
ax1.axes.get_xaxis().set_visible(False)
annotate('$n \qquad  \qquad t=%s \qquad \qquad \qquad C(t)=f(t)$'% '\int_{Co}^C \\
frac{dC}{kC^n}',xy=(0.1,1))
ax1.add_artist(Line2D((0, 0.9), (0.95, 0.95), color='black', linewidth=2))
ax1.add_artist(Line2D((0, 0.9), (0.05, 0.05), color='black', linewidth=2))
for i in n:
    s=integrate(1/(-k*C**i),(C,Co,C))
    annotate('$%.1f\qquad  \qquad   t=%s$'%(i,latex(s)),xy=(0.1,(i+1)/4+0.1))
    s=s-t
    s=solve(s,C)
    if i==-1:
        annotate('$C(t)=%s$'%latex(s[1]),xy=(0.6,(i+1)/4+0.1))
    else:
        annotate('$C(t)=%s$'%latex(s[0]),xy=(0.6,(i+1)/4+0.1))
show()
  • Voltar

    • Matemática Numérica

    # -*- coding: utf-8 -*-
from __future__ import division
from scipy.integrate import odeint
from pylab import *

#Simulação dinâmica: Enchimento de um tanque cilindrico
#Parâmetros
A=1			#m2
Qe=1 	        	#m3s-1
#Variáveis
ho=0 	        	#m
t=arange(0,100,10)
#Equações
dh=lambda h,t:Qe/A
#Resolução
h=odeint(dh,ho,t)
#Representação
figure()
xlabel('tempo (s)')
ylabel('Altura de líquido (m)')
z=lambda x: ho+Qe/A*x
plot(t,h,'wo',label=u"$Solução\ Numérica$")
plot(t,z(t),'b--',label=u'$Solução\ Analítica$')
legend(loc=2)
show()

    # -*- coding: utf-8 -*-
from __future__ import division
from scipy.integrate import *
from pylab import *
result1=quad(lambda x: 1/(1-x),0,0.9)
result2=quad(lambda x: 1/(1-x)**2,0,0.9)
figure(facecolor='white')
ax1=axes(frameon=False)
ax1.get_xaxis().tick_bottom()
ax1.axes.get_yaxis().set_visible(False)
ax1.axes.get_xaxis().set_visible(False)
annotate('$\int_0^{0,9}\\frac{1}{1-x}\ dx=%.2f \\pm %.3e$'%(result1[0],result1[1]),xy=(0.1,0.9))
annotate('$\int_0^{0,9}\\frac{1}{(1-x)^2}\ dx=%.2f \\pm %.3e$'%(result2[0],result2[1]),xy=(0.1,0.75))
show()
  • Voltar

  • Voltar
    • -->