CFD中文网

DEM加速比很高
CFD好像一般

@东岳 谢谢！

压差是分析很多工程设备的关键参数，在之前，主要都是针对多孔介质的Ergun模式来处理这一块。
也跟试验对比很好。但是，最近新开发的一个模型，虽然与试验的结果对比比较好，但感觉总是有问题。
于是不知道CFD中静压值到底是怎么得出来的？
对于一块多孔介质好像比较好处理，但是如果这个时候多孔介质表面沉积了灰，这个时候这块的整个压降会升高，然后灰越堆越多整个表面的压降都会升高。然而，这个局部压降和表面灰的分布十分相关，造成整个表面的这个灰分的压降相差一千多Pa。而最后一般的检测压降点，都是压差变松器的两个接口，一般都是设置在容器外壁的。但是如果我模拟检测点和试验一样的时候，这个时候的压差比我的表面计算压差平均值要小得多。（压差是1200多Pa，但是我的多孔介质有800多Pa，加上平均粉尘层的压降有640多Pa，合起来应该有1400多Pa，而不是1200多Pa.）
这就产生了一个问题，这个通过N-S方程求解的压力项到底是怎么处理的?
当他通过一层多孔介质的时候，多孔介质内部有压力源项，而在外部没有，那么当他气体从多孔介质外向内流通的时候，他的这层压力和速度是怎么计算的？如果这个时候我再添加一个源项在多孔介质层，而这个源项在多孔介质外是没有的，那么这层额外的源项是不是也会通过增加流体的静压来表达呢？（通过伯努利方程似乎很好理解，但是N-S是怎么实现的，他的这个\Delta P 是怎么获得的？）

这个是哪国人的名字哟

@mohui 是的，interface处是有间断的，但是one-field形式的控制方程体现不出来间断
压力分布不对，看一下理论解就明白了～

其实这个东西都是很好保存的，如果你的值是网格相关或者是面相关的话，直接利用UDMI就可以很好的实现你的值的随时间保存。
而你所说的EDFINE_RW_FILE是一个把运算数据保存到data文件的宏，但是你必须要注意的是，这个宏的保存是严格顺序的，你怎么样保存顺序，必须再一样的顺序输出。
很多时候，最简单的办法是直接你用C语言输出数据到一个文件中，然后再用ON_DEMAND读取便可。

啊客母~

硕士毕业找到了一些研究所的工作，还是拒了选择读博。因为自己总觉得还有很多东西值得深入地学习。读博是个难得的提升自我的机会，比较enjoy独处的时候可以静下心来思考。生活也不算惨淡，有funding可以不愁吃喝。

@cfd-china 是挺火的，据说那些最火的网红或者游戏主播一年挣几百万甚至上千万。。。

@cfd-china
有钱的话:cheeky: ，可以去构建一些关系嘛。但是这个也有点难。慢慢来嘛:happy:

5.0版本的可以正常安装，但是低版本如2.4.0安装就会报错：

-bash: /public/openfoam/OpenFOAM-2.4.0/etc/bashrc: line 140: syntax error near u nexpected token `;&' -bash: /public/openfoam/OpenFOAM-2.4.0/etc/bashrc: line 140: ` [ "$FOAM_V ERBOSE" -a "$PS1" ] && echo "Sourcing: $1" 1>&2'

@邙山老王 看图

是的你说的很正确。

你不觉得大材小用了，浪费那么多的时间搞科研，不说为人类的进步献一份力把，至少让其学有所用啊，他的那些东西中学根本用不上。 人生价值何在啊。何必当年读博士。

Simulations of gas-liquid dispersions by four-way coupled E-QBMM with aggregation and breakage processes Abstract Results Main References Job: Data processing Test cases debugging and running Writing paper in English if you are willing to Requirements: Benefits: You can get around 30000 yuan bonus from your host (if you write the manuscript). You can use this paper to get promoted faster (if you write the manuscript). Another 10000 yuan from my side If I write the paper, I can support you 40000 yuan as well. However, I dont have much time to write. So I prefer you writing the paper and to be the first author. Project duration

surfaceScalarField cPhi("cPhi", c & mesh.sf()); fvScalarMatrix c ( fvm::ddt(c) + fvc::div(cPhi) == source. );

@李东岳 哦，明白，谢谢李老师

四面体分解确定插值权重啊，用于相间耦合或者场平均（跟用多少耦合有关）。跟重心关联，还可以表示成局部坐标的形式。

带摩擦的运输方程是 u_t = lambda * u_x = eps * u_xx, 通过变化, 我们得到了, 0_1545168405877_3bd01f52-6748-4ece-bac9-e55eeef96e15-image.png ,其中D是导数矩阵, D2 是二阶导数矩阵
然后, 我们对M进行求特征值(EW), 得到了这种图片, 我想知道这是为什么呢?
0_1545168538892_dba5eac1-7474-4ce6-b089-77d3d8903c03-image.png