python - 试图同时求解微分方程
问题描述
我正在尝试为化学反应器设计构建一个代码,该代码能够解决反应器的压降、转化率和温度问题。所有这些参数都有微分方程,所以我试图在一个函数中定义它们,以便能够使用 ODEINT 对它们进行积分。然而,我构建的函数似乎有一个错误,我无法弄清楚是哪个让我无法集成它。
我遇到的错误:
---------------------------------------------------------------------------
TypeError Traceback (most recent call last)
<ipython-input-32-63c706e84be5> in <module>
1 X0=[0.05,1200,2]
----> 2 y=func(X0,0)
<ipython-input-27-6cfd4fef5ee2> in func(x, W)
8 kp=np.exp(((42311)/(R*T))-11.24)
9 deltah=-42471-1.563*(T-1260)+0.00136*(T**2 -1260**2)- 2,459*10e-7*(T**3-1260**3)
---> 10 ra=k*np.sqrt((1-X)/X)*((0.2-0.11*X)/(1-0.055*X)*(P/P0)-(x/(kp*(1-x)))**2)
11 summ = 57.23+0.014*T-1.94*10e-6*T**2
12 dcp=-1.5625+2.72*10e-3*T-7.38*10e-7*T**2
TypeError: unsupported operand type(s) for -: 'int' and 'list'
这是完整的代码:
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
fi = 0.45
gas_density = 0.054 #density
Pres0 = 2 #pressure
visc = 0.09
U = 10
Ac = 0.0422
T0 = 1400 #and 1200 also
gc = 4.17 * 10**8
bed_density = 33.8
Ta = 1264.6 #wall temp
fa0 = 0.188
def func(x,W):
X = x[0]
T = x[1]
P = x[2]
P0 = 2
R = 0.7302413
k = np.exp((-176008 / T) - (110.1 * np.log(T) + 912.8))
kp = np.exp(((42311) / (R * T)) - 11.24)
deltah = -42471 - 1.563 * (T - 1260) + 0.00136 * (T**2 - 1260**2) - 2,459 * 10e-7 * (T**3 - 1260**3)
ra = k * np.sqrt((1 - X) / X) * ((0.2 - 0.11 * X) / (1 - 0.055 * X) * (P / P0) - (x / (kp * (1 - x)))**2)
summ = 57.23 + 0.014 * T - 1.94 * 10e-6 * T**2
dcp = -1.5625 + 2.72 * 10e-3 * T - 7.38 * 10e-7 * T**2
dxdw = 5.31 * k * np.sqrt((1 - X) / X) * ((0.2 - 0.11 * X) / (1 - 0.055 * X) * (P / P0) - (x / (kp * (1 - x)))**2)
dpdw = (((-1.12 * 10**-8) * (1 - 0.55 * X) * T) / P) * (5500 * visc + 2288)
dtdw = (5.11 * (Ta - T) + (-ra) * deltah) / (fa0 * (summ + x * dcp))
return [dxdw, dpdw, dtdw]
X0 = [0.05, 1200, 2]
y = func(X0, 0)
提前致谢
解决方案
内线
ra=k*np.sqrt((1-X)/X)*((0.2-0.11*X)/(1-0.055*X)*(P/P0)-(x/(kp*(1-x)))**2)
您可能想要使用...X/(kp*(1-X))...
而不是...x/(kp*(1-x))...
(即使用上 X),下x
是列表类型。
如果您想在l
某处使用某个列表变量作为多个值,则将其转换为 numpy 数组la = np.array(l)
并la
在 numpy 矢量化表达式中使用。
推荐阅读
- ios - UITableView.visibleCells.contains(JDHeaderTableViewCell()) 总是返回 false
- javascript - 如何修复“鼠标按住可拖动元素”
- jenkins - 如何在下游作业中读取上游作业发送的参数?
- r - 在 R 中创建分类变量
- android - 如何在java中使用ndk获取证书值
- python-3.x - python定义一个地方到处使用
- javascript - 在 JSDOM 中使用滚动和宽度
- set - 编程语言中集合背后的理论是什么?
- ruby-on-rails - 在 Rails 控制器中拆分参数
- javascript - NodeJS如何将文件转换为字符串