哈哈哈，没事，谢谢东岳老师！

差不多，U[cellID].x(), U[cellID].y(), 之类

是这样的，我想在求解过程中使用python做一些处理，要把流场数据传到python，处理后在传回，关于在of求解过程中调用python，我找到了一篇帖子：OpenFOAM 与 python 代码之间的通信
但是volScalarField这样一个类，python肯定不能识别，请问我如何能单独把流场的值作为double数组提取出来？
如果不能把整个流场值变成一个double数组，能每个cell挨个提取出来也是极好的。
我尝试去代码里找volScalarField的结构，但模板套模板实在没读明白，希望各位老师大牛们能不吝赐教。

@wwzhao 嗯嗯，明白了，非常感谢您的解答！

哦哦，感谢
如果是吸收到压力梯度项里去了，那需要对压力梯度项或压力方程做修正吗？

在构建湍流模型时，雷诺应力根据Boussinesq假设可写为

S-A模型中，建立了涡粘系数nut的输运方程，可由此得到nut，那么上式中的k是怎么处理的？
看的几个资料中都没有提及这个部分，请各位老师解个惑

是我应该感谢cfd中文网的各位老师，对我的帮助很大.

另外，我测试的这个构造函数和析构函数只承担了new和delete一个double[]的任务，如果是功能更复杂的构造函数，就不好说了。

是的，调用1000000次构造函数和1次，带来的差别不大

@Zephyr 嗯嗯，感谢，我去尝试一些库

谢谢您的回答！我试了下不使用动态分配而用固定大小的数组p[10000]来代替，确实提高了一些效率，从10S+变成了7S+。
但我这里确实需要动态分配，因为我希望做成一个有一定通用性的矩阵类来储存各种流场信息，我无法在编程时就知道需要多大内存，也许这次计算只要几千网格，也许下次就需要几千万网格。