我有一个动态模型设置为ODE的(刚性)系统。我目前用CVODE(来自Assimulopython包中的SUNDIALS包)解决了这个问题，一切都很好。我现在想为问题添加一个新的3D散热器(具有与温度相关的热参数)。我的想法是使用现有的FEM或FVM框架为我提供一个接口(interface)，让我可以轻松地提供(t,y)将3Dblock转换为例程，并取回残差y'。原理是使用FEM系统中的方程，而不是求解器。CVODE可以利用稀疏性，但预计组合系统的求解速度比FEM系统自身求解的速度要慢，专为此类问题量身定制。#pseudocodeofaresidualsfunctionforCVODEd