https://github.com/OpenFOAM/OpenFOAM-3.0.x/blob/master/src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C

去这里加个源项不就好了？

// Turbulent kinetic energy equation tmp<fvScalarMatrix> kEqn ( fvm::ddt(alpha, rho, k_) + fvm::div(alphaRhoPhi, k_) - fvm::laplacian(alpha*rho*DkEff(), k_) == alpha*rho*G - fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_) - fvm::Sp(alpha*rho*epsilon_/k_, k_) + kSource() // ? + K*k_ );

@东岳 嗯嗯 谢谢东岳老师！

@buaa_hgj 非常感谢

最近在看dynamicRefineFvMesh类，发现他和普通的八叉树网格细化有一些区别，有没有人知道这个算法最早是谁写的？有没有相关的参考文献？

@bestucan 时间太短了，烧不起来，只能强算，往后算

@史浩 ANSYS用的是单元生死技术，帝国理工好像是移动标架法。

@昝磊 @史浩 百年好合！ :laile:

多谢，已更新。

关于你的单核成功，并行发散的问题，我遇到过很多次，但没有一阵见血的答案 :papa:

@东岳 万分感谢！

今天在CFD-online上看到一个帖子，也在问of2.1.1中的cyclic边界条件，看了一下他遇到的问题跟我一样（包括使用的求解器也一样），所以现在认为应该不是我的Case中设置的问题，应该是openfoam2.1.1中cyclic边界条件本身存在bug。打算从一个高版本的of中copy出cyclic边界条件，再自行编译为用户自定义边界条件，并没有成功。目前还没想到更好的解决办法。

@东岳 谢谢您！

早期的图形界面一般在后面，tty6还是7.
前面的都是终端界面。
后来改了，tty1是加载信息的界面，就是开机的时候一堆字码。最终是登录界面的样子。tty2是登录后的图形界面，往后才是终端界面。
而界面切换的工作量还是有的，如果当前系统有什么占用磁盘读写的程序在卡着，切换会相当慢。甚至切到终端界面，打字母都有一两秒延迟。如果看见没反应持续操作，就会给卡死的系统卡的更死。所以，和平解决方案，就是--切一下，然后等。进top杀死高占用，或者直接重启

多谢，已将

将源代码包拷贝到ThirdParty-6文件夹中并解压缩，重命名为ParaView-5.4.0，终端切换到paraview目录下并键入：

改为

将源代码包拷贝到ThirdParty-6文件夹中并解压缩，重命名为ParaView-5.4.0，终端切换到ThirdParty-6目录下并键入：

嗯嗯，矢量分量等值面有了，谢谢！

能看一下你完整的设置文档么，这个按道理小了容易出问题，大了反而不容易吧。

发现了几种方法都可以 一种是在ControlDict文件下加入FunctionObjects,并链接"libfieldFunctionObjects.so"，然后对边界面上的质量通量phi求和，
另一种方法需要编入swak4Foam,9b116a8d-2652-41d4-8fb0-a6674bbdb5b8-image.png

网络问题，重新执行sudo apt-get update 试试

不清楚。我一般通过phaseProperties里面的type pureIsothermalPhaseModel关掉传热。如果不关掉，肯定温度和压力耦合一起了。
我觉得可能跟你的进出口压力边界条件有关系。

@danao 在 源项系数是速度的绝对值 中说：

速度的模也是有单位的

嗯，我也忘了，感谢分享！:xiexie:

速度的模也是有单位的

所以在写代码的各种量都带上量纲，也是一种防止错误发生的方法

勘误：计算所用的我现在用是一个长方形的几何
计算散的我精神已经恍惚了，明明是一个很简单的问题:xinlei: