CFD中文网

<a href="#" class="icon icon-home"> 东岳流体|OpenFOAM|CFD <img src="_static/logo.svg" class="logo" alt="Logo"/> </a>

img一行代码缺失，需要加进去

因为分散力模型无法找到，https://zhuanlan.zhihu.com/p/90057423

@李东岳 谢谢！

Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: sixDoFRigidBodyMotion Selecting sixDoFSolver symplectic Translational constraint tensor (1 -0 -0 -0 1 -0 -0 -0 -624) Rotational constraint tensor (0 0 0 0 0 0 0 0 0) PIMPLE: No convergence criteria found PIMPLE: No corrector convergence criteria found Calculations will do 2 corrections PIMPLE: Operating solver in transient mode with 2 outer correctors Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Selecting laminar stress model Stokes No MRF models present No finite volume options present Constructing face velocity Uf Courant Number mean: 0.000625613 max: 0.10552 fieldAverage fieldAverage1: Starting averaging at time 0 forces forces1: Not including porosity effects forces forceCoeffs1: Not including porosity effects forceCoeffs forceCoeffs1: Not including porosity effects Starting time loop Courant Number mean: 0.000625613 max: 0.10552 deltaT = 0.005 Time = 0.005 PIMPLE: Iteration 1 forces forces: Not including porosity effects #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::sixDoFSolvers::symplectic::solve(bool, Foam::Vector<double> const&, Foam::Vector<double> const&, double, double) at ??:? #4 Foam::sixDoFRigidBodyMotion::update(bool, Foam::Vector<double> const&, Foam::Vector<double> const&, double, double) at ??:? #5 Foam::sixDoFRigidBodyMotionSolver::solve() at ??:? #6 Foam::motionSolver::newPoints() at ??:? #7 Foam::dynamicMotionSolverFvMesh::update() at ??:? #8 ? in "/home/z/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam" #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 ? in "/home/z/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/pimpleFoam" 浮点数例外 (核心已转储)

@李东岳 李老师，请问如果求解非守恒型的动量方程，那么对于其中一相为0的情况，是不是就直接求解就好了，不用再做别的处理了，如让相为0的相的速度等于按照模型去求出的速度。

下载个手机版的向日葵，如果手机也上不去大概率就是国外网络问题。或者试试看 Todesk 这个远程软件能用不

@AppleKiller 感谢，我再看看，以前也接触过这些资料，感觉用欧拉法求解器为基础修改的拉格朗日求解器不大好，没怎么深入去学习。我再看看这些思路，感谢

@lyc 你好，我很少遇到。即便遇到都是重启就能解决的那种。

1.JPG
src文件夹下没有fvoptions，但想利用其中的CellSetOption功能，有人知道吗？

@Qing 感谢老师！您说的特别清楚，这对我特别有帮助。再次谢谢您！

没有指导意见，我现在用的sctoch划分，说的是划分后交界面最小，好像也没比simple快

@李东岳 感谢李老师，我试试看能不能画出来

@羽之下 thank you，谢谢你的回复。我记得好像在哪看过说这里引用了一个假设，我决定把这部分系统的看一看，有什么新发现我会继续在这里交流~

@东岳 对，我那个5118是学校的计算集群，我记得我在哪儿看到过说这堆集群都是把超线程关了的，所以不可能使用的线程数量超过物理核心数量；然后睿频的话，应该是说如果你不关他，他默认是开的，但学校集群可能是关了的，你的如果是自己用应该是开着的，就是说如果睿频开着的，CPU温度和功率不超过他的限值的时候他可以运行在更高的频率上。比如你的5120 是14cores 28 threads 2.2G, 但是如果散热比较好功耗不超过某个限制的时候，他是可以整体跑在更高的频率上，这就能更快。
话说我觉得最惊奇的是我的9750H跑6线程居然比服务器CPU慢这么多，理论上说不过去啊，因为我看到跑的时候主频都跑到4G了，服务器才2.3G。得研究一下构架了

感觉DDPM比较适合你的例子，就是计算慢了些，或者用双欧拉

给你一个搜各种书的网站，不需外网。一直用，偶尔令人失望。
https://vdoc.pub/

可以通过TUI实现（在控制台窗口中输入文本命令）

/mesh/modify-zones/copy-move-cell-zone

缺点是先copy再move，可以再用

/mesh/modify-zones/delete-cell-zone

删去移动前的网格域

@悟空空 我也是这样的问题，直接复制并行命令就出错，但输入并行命令就没问题，这是什么原因啊？

其实我只是想实现tecplot里面的ddx函数，但是我找不到tecplot讨论区，所以先发这里了。
我得到的数据变量在储存在节点，而不是cell中心。网格是二维四边形网格。因为是非结构的，所以不能像笛卡尔网格那样直接求差商。我看了tecplot帮助文档里面ddx的实现方法是Moving Least-Squares，说是假定变量关于节点是二次函数，之后最小二乘法拟合出二次函数，然后将二次函数在节点的梯度近似为该变量在节点的梯度。我不明白对于二维的二次函数，tecplot是否考虑了交叉项xy?
如果不考虑的话，那我是否可以偷懒:假设需要求f的梯度，求对x的偏导的时候，完全不管y坐标。对单元的四个节点的x坐标x1，x2，x3，x4，和他们的函数值f1，f2，f3，f4进行一维的二次函数拟合。如果得到的拟合结果是f(x，y)=ax^2+bx+c，那第一个节点的x方向梯度就是把x1带入f'(x)=2ax+b中？

手机打字不太方便，请见谅。

@李东岳 东岳老师您好，我最近发现了一个现象可能造成了我上述提出的问题哈哈哈哈哈，我最开始使用的是DPMFoam默认的的Ergunwenyuforce曳力，我后来换成了sphereforce曳力就能稳定了，而且对于颗粒运动和流场有挺大影响，看来这和我的设置也有一些关系。（最开始图方便就直接使用的默认的模板了，看来还是要根据实际情况来改变）