目前,我正在对恒容反应堆-RCM进行建模。我正在努力设定墙的速度。我的教授告诉我,必须使用“体积-时间”历史记录配置壁速。我不知道如何将这两列作为配置文件。
import sys
import cantera as ct
import numpy as np
gas=ct.Solution("gri30.xml")
reac = ct.IdealGasReactor(gas)
env = ct.Reservoir(ct.Solution('air.xml'))
net = ct.Wall(reac, env, velocity=)
netw = ct.ReactorNet([reac])
我有一个这样的文本文件。有一个时间和音量列。
VPRO 0.0000001 1.0000000
VPRO 0.54822008E-02 0.99950299
VPRO 0.65802743E-02 1.0026889
VPRO 0.73418149E-02 0.99627431
VPRO 0.90739698E-02 1.0158893
VPRO 0.96140946E-02 1.0028384
VPRO 0.97804742E-02 1.0070171
VPRO 0.10084646E-01 0.99990454
VPRO 0.10693270E-01 1.0107573
VPRO 0.11517717E-01 1.0099201
VPRO 0.11831626E-01 0.99879979
VPRO 0.29057992E-01 1.0034597
VPRO 0.29139576E-01 1.0001761
VPRO 0.29888237E-01 0.99735919
VPRO 0.30740998E-01 1.0048461
VPRO 0.30806940E-01 1.0092222
VPRO 0.31839109E-01 1.0053230
VPRO 0.32515062E-01 0.99802537
VPRO 0.32889860E-01 0.99232409
VPRO 0.35364048E-01 1.0027345
VPRO 0.35843734E-01 0.99307568
VPRO 0.36135389E-01 0.99638039
VPRO 0.36411413E-01 0.99009293
VPRO 0.37458905E-01 1.0028019
VPRO 0.39592186E-01 0.98974443
VPRO 0.40163110E-01 1.0000900
VPRO 0.43210211E-01 1.0070476
VPRO 0.49542261E-01 1.0026439
VPRO 0.49759985E-01 1.0077359
VPRO 0.51519470E-01 1.0026126
VPRO 0.52322467E-01 0.97408403
VPRO 0.52944216E-01 0.96103576
VPRO 0.54105348E-01 0.93095020
VPRO 0.55044882E-01 0.89119653
VPRO 0.55911692E-01 0.84723777
VPRO 0.56457098E-01 0.81692500
VPRO 0.57248888E-01 0.77123227
VPRO 0.57805937E-01 0.73170784
VPRO 0.58648872E-01 0.66863047
VPRO 0.59157635E-01 0.62958903
VPRO 0.59754123E-01 0.58326969
VPRO 0.59957977E-01 0.56727436
VPRO 0.60453182E-01 0.52899850
VPRO 0.60906861E-01 0.49480953
VPRO 0.61257981E-01 0.46852022
VPRO 0.61704079E-01 0.43503065
VPRO 0.61903909E-01 0.41910567
VPRO 0.62206198E-01 0.39576119
VPRO 0.62507156E-01 0.37251365
VPRO 0.62805821E-01 0.34862793
VPRO 0.63055253E-01 0.32943590
有人可以建议我为此编写代码吗?我在某处读取了这些数据,并将其存储在关键字字典中。然后,我可以使用它来设置壁速度。有人可以帮我吗?