那就初始化把多孔介质中加入这些 specie 啊，初始化成零。 按说也应该初始化，不过有的计算过程默认初始化成零。 估计是没有初始化时，物料运输方程都没有在多孔介质区域离散。就没有计算了。

@李东岳 谢谢老师 回头统计地址的时候我报上

#!/bin/sh # Source tutorial run functions . WM_PROJECT_DIR/bin/tools/RunFunctions # Get application name application="potentialFoam" runApplication blockMesh \cp save/boundary constant/polyMesh/boundary \rm -rf 0 \cp -r 0_org 0 runApplication $application -writep runApplication writeIbMasks 这个是cylinderInChannelPotential中Allrun文件代码，就是比较想知道save文件夹里边有block Mesh，constant文件夹下的ployMesh里也有blockMesh，所以Allrun中的blockMesh命令运行的是哪个文件。然后将save/boundary 拷贝到constant/polyMesh/boundary的意义是什么呢

这是一个非常经典的大涡模拟的算例。简直是打遍天下无敌手。做LES的，一般看到这个经典几何，立马就会啪的一下打开文献，很快啊！左边一个结构块，右边一个结构块。把网格画出来算一算。并且计算结果可以和实验对比的非常好： 在这个算例中，相关条件大体如此： 左、右、前、后均为周期边界条件 上下为壁面，速度无滑移，压力零梯度，湍流采用壁面函数 湍流模型使用SpalartAllmaras 粘度$\nu$为2.643e-6，空气的密度按照1来算 传统功夫点到为止，整个压力-速度耦合迭代1500次就好 网格一共256万。酷炫吊炸纯结构网格。OpenFOAM这面采用16核（自用服务器处理器AMD 7542）计算需要533秒。很快啊！ 有没有人感兴趣用Fluent算算？感兴趣的可以在这里下载网格，说一下CPU配置，以及计算时间，更多信息可再交流 网格链接：https://pan.baidu.com/s/1C_VV6Yai6b63kT6MHKmRhQ 提取码：aywm Xiao and Jenny, A consistent dual-mesh framework for hybrid LES/RANS modeling, Journal of Computational Physics, Volume 231, Issue 4, 20 February 2012, Pages 1848-1865

也是实测。OpenFOAM自带算例。256万网格。虚拟机16核跑了1772s，真系统跑了548s。 虚拟机大约要慢3倍

@zhe 我没有，而且MPPIC调用的是kinematicCloud，stochastic collision model中的很多properties是无法调用的比如下面这一段涉及到thermo的地方就是没有办法调用的。 bool massChanged = collideParcels(dt, p1, p2, m1, m2); if (massChanged) { if (m1 > rootVSmall) { const scalarField X(liquids_.X(p1.Y())); p1.setCellValues(this->owner(), td); p1.rho() = liquids_.rho(td.pc(), p1.T(), X); p1.Cp() = liquids_.Cp(td.pc(), p1.T(), X); p1.sigma() = liquids_.sigma(td.pc(), p1.T(), X); p1.mu() = liquids_.mu(td.pc(), p1.T(), X); p1.d() = cbrt(6.0*m1/(p1.nParticle()*p1.rho()*pi)); } if (m2 > rootVSmall) { const scalarField X(liquids_.X(p2.Y())); p2.setCellValues(this->owner(), td); p2.rho() = liquids_.rho(td.pc(), p2.T(), X); p2.Cp() = liquids_.Cp(td.pc(), p2.T(), X); p2.sigma() = liquids_.sigma(td.pc(), p2.T(), X); p2.mu() = liquids_.mu(td.pc(), p2.T(), X); p2.d() = cbrt(6.0*m2/(p2.nParticle()*p2.rho()*pi)); } } Tut中的aachenbomb 在 stochastic collisionmodel 设置一个banana 就能看到Valid model了。 stochasticcollision model是spraycloud独有的。 而kinematiccloud中的collisionmodel 更多的倾向于是关于软硬球的处理，而且应该用的是deterministic tracking method。 所以如果你要在MPPIC下使用stochasticcollisionmodel的话，就要自己研究怎么添加模型了，基本上相当于自己全程写东西了

@队长别开枪 专业 学习了学习了 我看朋友用blender做过一些美术动画，看起来还是挺舒适的 目测不太适合工程

@同学博 是的，基于 freeCAD，专门给 OF 开发的前处理工具

有偿求助！看了些论文，只看懂在多孔介质区域设置为非平衡项。，具体的参数又是如何推导出来呢？

@李东岳 谢谢李老师！我发现问题了，这个报错其实就是最后的correctADDStress的问题，所以看这个报错提示就是看最后一个信息，错误的原因是我用wclean去移除编译的库的时候其实没有移除，所以后面不管我怎么改，都是调用的原来的库，解决办法就是手动删除openfoam-7的platform中生成的自己编译的库，重新编译就好了

这种情况，我一会再测试一下Crc和Cti是否参与了计算，谢谢老师百忙之中抽空解答我的疑惑，谢谢

can@x260-TC ~/.local/share/ansys_inc/v181/fluent/fluent18.1.0/src %grep -R 'NV_DS' ./ ./turbulence/turb.h: d10 = NV_DST2(xnode0,xnode);\ ./pbns/flow.h:#define NV_DS(a,EQ,x0,x1,x2,S,s)ND_VEC((a)[0]EQ(x0)S(s),(a)[1]EQ(x1)S(s),(a)[2]EQ(x2)S(s)) ./pbns/flow.h:#define NV_DST2(x1,x2)ND_OP(+,SQR((x2)[0]-(x1)[0]),SQR((x2)[1]-(x1)[1]),SQR((x2)[2]-(x1)[2])) ./pbns/flow.h:#define NV_DST(x1,x2)sqrt(NV_DST2(x1,x2)) ./pbns/flow.h:# define N3V_DS NV_DS NV_DS 由ND_VEC定义 ND_VEC的定义在这个文件里 ./main/global.h:# define ND_VEC(x,y,z)x,y,z ./main/global.h:# define ND_VEC(x,y,z)x,y "."开头就是ansys_inc/v181/fluent/fluent18.1.0/src这里开头，这里就是你的ansys安装的位置。

可能是超算上有不只一个节点，所以节点间通过IP地址交流信息。 而linux系统里有一套“权限”和“环境变量”的继承规则。你加载 openfoam 环境变量的的那句应该是放在家目录的.bashrc了（就是source ~/OpenFOAM/etc/bashrc，大概是这句） 而 miprun 再通过IP调用其他节点的时候，这个网络通信的用户权限和环境变量，应该没有继承你的用户名的环境变量；而是系统的环境变量。所以不包含你的用户目录下.bashrc里的东西。 把你家目录那句 source 放到 /etc/profile 试试 。这是系统范围的环境变量。得是绝对路径：/home/<用户名>/OpenFOAM/…… 当然也可能得在各个节点上都配上环境变量、装上 openfoam，这细致活得慢慢调…… 部分求解器出现上述状况（编译成功，但是不能算），还有一部分求解器完全正常。 没用过多节点的，但问题是出在环境变量上没跑了。 到这个路径下 ~/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin 这都是可运行的二进制文件。 直接./<文件名>，比如 ./icoFoam 就知道是没编译好还是环境变量没配置好了。

俺觉得你说滴对 希望我有一天也能像你这么强

一年多了还是没有解决这个算例，不知道有没有人做过相类似的，重新描述下这个算例，如图所示一个后台阶流动，但是是在图中左下（0,0,0）位置为原点绕Z轴旋转，主要流动参数， 雷诺数$Re_b=\frac{U_b(H-h)}{\nu}=5600$，其中$U_b$是体积平均速度； 旋转数$U_b=\frac{\Omega(H-h)}{U_b}$，$\Omega$是角速度； 台阶前段槽道流部分壁面剪切雷诺数$Re_\tau=\frac{u_\tau(H-h)}{\nu}=360$，这个与槽道流DNS结果是完全一致的，也就是说台阶前是完全发展的槽道流。目前存在以下问题： 1、求解器选择，SRF类还是需要自己修改求解器？ 2、合理的进口条件，不旋转的相对好给，但是旋转的case，速度进口边界条件和初始值不知道如何给合理。 3、无量纲化的流向速度为什么会小于1，如下图所示。以我的经验，槽道流里面最大速度比体积平均速度肯定是大的，也就是$\frac{U_{max}}{U_b}$是大于1的，但是为什么在旋转流动里面就会于1呢？ 参考文献 1、DNS原始文献：Turbulent flow over a backward-facing step. Part 1. Effects of anti-cyclonic system rotation 2、图片来自：Development and Application of Rotation and Curvature Correction，是篇博士论文

有人用Avizo画网格吗？

@zlwei 好的 感谢回复

@发芽的土豆 也捧个场。只捧场 (我水平低，甚至不知道，你这的“湍流”是强调 turbulence 呢还是 vorticity) “豆”兄好运

@liujunCFD 云岗那边3部 31所什么的

可以在BIOS里面关掉超线程，计算速度感觉快了一点，但是没有详细测试。。。

来，骗！ 来，偷袭！

这是典型的处理方法。先用各种酷炫的格式求hyperbolic的欧拉方程（Roe/AUSM等）。后面的parabolic慢慢搞。

http://c.biancheng.net/view/2226.html this就是当前对象的别称，有哪些成员，就看它所在的函数是哪个对象的成员函数了。 去定义这个对象的头文件查，与之同名的.c文件有各成员函数的定义。

十分感谢，已经解决了

感谢回复，我研究研究。我把体积分转化成了面积分，然后用python calculator进行了求和。感觉结果不太理想

@发芽的土豆 : 感谢赐教，多谢了

@dzw05 你遇到的极值局部点距离给定压力条件的边界接近吗？我今天下午算的时候发现时间步过大时，距离fixedValue p_rgh边界最近的cell中心也会出现温度极值....调小了时间步后就没有了。我算的这段时间一直感觉温度和压强的联系太密切了...

没币发帖

不好判断 但感觉就是设置的问题 如果可以上传算例并且网格数量小的话 我可以给你看看

@shepherd 谢谢！Openfoam Wiki上有of8的编译方式，我照着那个编译好了，of7也编译过，但是不知道为啥超过一定核数就没法计算，把of安装在root下还会有有的东西调用不成功的情况。目前主要是想再编译一下ESI的2006版本，把PETSc用在计算里边，之后还得实现一下跨节点的计算。

电势是相对的，上下机板电势差不变，就像电池正负两极保持电势差一样。不影响是否绝缘。 电荷密度零梯度，仅在极板与空气交界线上。看这种情况：紧挨着交接线两边各有10个电荷，所以零梯度。但是一边诗歌电荷后面站着100个电荷，另一边十个电荷后面站着20个电荷。 不同于流场的变化，仅分析相邻区域即可。电荷之间的相互作用是“超距”的，不相邻也可以用场去作用。 总的来看，电极板之间已有电势，只要有载流子就有电流。就像在电容间加水，就存不住电势差了。 交界处电荷梯度为零可以两边同时增加或者减少，或者增加减少动态平衡。 如果要绝缘，按照物理意义的绝缘最靠谱：没有载流子。 如果强行绝缘，应该是电势梯度为0，没有电势差，就没有电荷受力。但是电势场是场内所有电荷分布造成的结果，应该不能直接设定。 具体怎么实现得看源代码是怎么算电场的了

谢谢李老师！

文件路径OpenFOAM-7/src/fvMotionSolver/fvMotionSolvers/velocity/laplacian/velocityLaplacianFvMotionSolver.C Foam::tmp<Foam::pointField> Foam::velocityLaplacianFvMotionSolver::curPoints() const { volPointInterpolation::New(fvMesh_).interpolate //就是这里 ( cellMotionU_, pointMotionU_ ); tmp<pointField> tcurPoints ( fvMesh_.points() + fvMesh_.time().deltaTValue()*pointMotionU_.primitiveField() ); twoDCorrectPoints(tcurPoints.ref()); return tcurPoints; } 有没有具体的数学公式，我看了一些书也都没说这些细节的东西。

@veen 多谢大佬！^O^/

湍流模型针对的都是湍流流动，基本都是基于各向同性假设，考虑转捩流动应该用转捩模型吧，我了解不多，不知道您现在是不是有了比较深入的理解

@李东岳 wifi用5G就不会干扰了啊

@史浩 你好，抱歉，之前居然没注意到这个回复，刚刚才发现。现在“导出该球表面积随时间的变化”这个问题应该解决了吧？ 两种方法，一个是save data输出成csv，多个csv合并成一个文件，在excel里处理数据。另一个是plot data over time 也可以 另外，请问，表征形状是不是好的球体这个问题，解决了么？ 如果是少量的液滴，可以求一下质心坐标的偏差，但是你这个图的液滴数量比较多，取坐标什么的应该很麻烦。所以我很好奇这个问题有什么解决办法么？

@史浩 您好，请问对于自适应网格的设置讲解有什么资料可以帮助我学习参数设置的吗

![Y`IW6LQFCDZXZ5RHH%K69V.png**图为剖面图，两部分的网格没有connect，这里有两个内部面，如何将其缝合使网格节点对应 **

取最小网格的最小dt，当做全场的dt

@nanxuan 不好意思你九月份的留言我没注意到。多边形/多面体网格相对传统非结构网格优势就是每个体单元有了更多的面单元/相邻单元，在涉及面单元积分的变量（梯度等）计算上具有优势（https://www.semanticscholar.org/paper/The-advantage-of-polyhedral-meshes-Perić-Ferguson/51ae90047ab44f53849196878bfec4232b291d1c ），其余的还有收敛快等优点（https://www.researchgate.net/publication/292555771_New_element_lops_time_off_CFD_simulations ）。缺点是生成算法复杂，主流多面体网格生成器都是将四面体网格转化为多面体网格，不能很好处理面单元翘曲，不过fluent，cfmesh等通过六面体->四面体->多面体的技术路线缓解了面单元翘曲，但是面对复杂几何外形的时候需要更多控制。我个人开始使用多面体网格的原因纯粹是因为多面体网格好看，可以在女朋友面前装逼 。希望这个解释能帮到你。

@赵一铭 你好，请问对于并行计算的案例该如何用gnuplot实时监测残差呢？

对，看看第一个error，看起来像是编译器的原因，但是貌似你确定编译器没问题

“我们现在做的很多科研成果都是错的”，这句话是有道理的

@红豆沙 fluent meshing画的多面体of是支持的，多面体和六面体核心结合的新型网格形式of不支持，不知道新版的有没有更新

那并行求解器也是人写的…… 用开源不就是为了对计算过程肆意妄为，想干啥干啥 当然能啊，只要想到就能。计算机就是帮助人类做重复无意义的工作的。 你描述的这个其实不用改求解器也行，只是用时间步的数据，算了完再处理也行。编程语言都能处理纯文本数据，或者matlab、excel

TEMOM只展开$(v+v_1)^{1/2}$，任意阶矩展开任意阶都是泰勒的截断误差，是个原式的高阶无穷小，原式也就是$1/2$次方。 但DEMM把$(v+v_1)^k-v^k-v_1^k$也带上展开了，这个截断误差是变的，k越大，误差越大。越高阶矩，误差的量级越大，虽然是都是无穷小。 所以我觉得越高阶这个方法的优势越不明显。当然都是我猜的 在提出DEMM的那篇文章里的数据也有这样的苗头，低阶矩符合很好，高阶的出点问题又修复了。 我把相关资料发给老师。

@李东岳 东岳老师，在kineticTheoryModel.C中有这样一段代码： Foam::tmp<Foam::volScalarField> Foam::RASModels::kineticTheoryModel::k() const { NotImplemented; return nut_; } Foam::tmp<Foam::volScalarField> Foam::RASModels::kineticTheoryModel::epsilon() const { NotImplemented; return nut_; } 但是感觉固相湍流开关打开的话，其湍流黏度完全是按照如下公式计算的： （其中a为固相，b为液相） 也就是固相被看作连续相后不存在真实的k和epsilon哇？