CFD中文网

@李东岳 这个案例我用的商软，发现空隙率0.01和0.5差别不大，影响最大的是阻力系数，我试了这两个空隙率，水都能很顺畅的流出来

@东岳 对于$\nabla{\mathbf{U}}$的矩阵形式应该是什么这个问题，从张量的角度来看我站第一种写法，更规范。:huahua:
但是从CFD角度来看，公式这东西本就是规则，制定了规则就是交流用的， 可读性和易理解性还是更为重要的。看个人理解，我觉得并无好坏。:mianmo:

前几天OpenFOAM基金会在群里面发了个消息，说OpenFOAM-9发布啦！然后一片死寂，一个人都没有，跟咱们中国的夸夸群差远了，于是基金会补上了一张图，就是这面这个

tumbleweed.gif

给我笑爆了..

Hello foamers:
如题
比如我想画水体的流场矢量图
那么应该怎样只把液相（体积分数>0.5）的速度场数据提取出来呢？
谢谢！

半新手，看其他论坛上可以实现这样的操作，即把一个复杂的几何体分成多部分分别划分网格之后再进行拟合，但是我自己试了下建块merge后导入到fluent中出错，请问该如何操作呢？不胜感激

主任招博士后不，我给您发个简历

这个问题我解决了，就是在周期性关联的时候，在edit periodic/periodic veteric 时要在中心轴上把所有的veterivic都要关联。即一个点要点两次

按照twoPhaseEulerFoam中的EEqns.H写出能量方程如下（算例文件中能量形式是sensibleInternalEnergy）大神们先看看这么写对吗？

fvScalarMatrix E1Eqn ( fvm::ddt(alpha1, rho1, he1) + fvm::div(alphaRhoPhi1, he1) - fvm::Sp(contErr1, he1) + fvc::ddt(alpha1, rho1, K1) + fvc::div(alphaRhoPhi1, K1) - contErr1*K1 + ( he1.name() == thermo1.phasePropertyName("e") ? fvc::div(fvc::absolute(alphaPhi1, alpha1, U1), p) + p*fvc::ddt(alpha1) : -alpha1*dpdt ) - fvm::laplacian ( fvc::interpolate(alpha1) *fvc::interpolate(thermo1.alphaEff(phase1.turbulence().mut())), he1 ) ); E1Eqn.relax(); E1Eqn -= ( Kh*(thermo2.T() - thermo1.T()) + Kh*he1/Cpv1 - fvm::Sp(Kh/Cpv1, he1) + alpha1*rho1*(U1&g) + fvOptions(alpha1, rho1, he1) );

a8753d25-8ea0-4cf2-8fae-10a71aa5cdfd-image.png
那如果不考虑传热的话，只把导热系数设为0和温度设为一样，这样做对吗？
我为了不考虑传热，是直接把twoPhaseEulerFoam.C中的EEqns.H给注释掉了，我这种做法对吗？

@liujuncfd
本科待遇年薪15w左右，硕士年薪20w左右，博士年薪25w。
https://www.zhipin.com/job_detail/fe4f05de98c982591nN62Ni7GVZT.html?ka=hot-job-1

感谢李老师！:xinxin:

想开个贴讨论下多相流界面追踪类方法目前的前沿水平，例如能算到的最高参数，各类方法都有什么限制等等。目前多相流算法里，计算界面追踪的多相动力学，比如droplet collision, droplet splashing, droplet impact on solid wall，这几个典型案例，目前有用VOF，level-set，phase field，LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态，实验都做出来了，但感觉计算还没跟上。

调研过不下百篇的文献，基本上关键词相关的都翻了一遍，在文献中看到droplet collision大概最好的也就是Re=6500,We=1500，密度比DR=1000左右的，大致是这个级别，像droplet impact有篇08年的用TransAT算到了RE=11000，而最近16年左右的用levelset这些也算到Re=8000，We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟，而且网格也都比较密，计算资源消耗也都比较大，估计都得上超算，当然和实验也都对的上，基本都是到了prompt splattering这样的形态。

想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪？ 因为多相流在算到高RE和We的时候会数值发散，稳定性应该是主要限制，还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况，需要单独使用人工识别薄液层再稳定的方式，对不规则碰撞并不实用。

在第一个时间步内生成了颗粒，怎么把颗粒的位置和尺寸并导出成表，并用表喷射器进行填充

@东岳 谢谢东岳老师，就是下面A beta Ta三个参数求化学反应速率常数，化学反应速率中反应物质浓度次幂直接在方程里面表现出来了哇，没有是默认为0次吗？

网格检查没有错，但是一运行就不行了，不懂啥意思啊。。。

等径管里如果出现激波那就直接堵塞了吧

我在使用solidDisplacementFoam求解器进行热应力计算时，需要对位移向量场重新设置边界条件：其中的x,y方向为zeroGradient边界，z方向为fixedValue边界。请问OpenFOAM里有类似的边界条件或者知道该怎么设置的吗？？

@blgy 谢谢您的回复，
1.第一个问题明白，我继续试一下。
2.对的，这个模型原来的确是采用混合网格，部分采用四面体，部分采用六面体，但是现在想整体都用六面体尝试，所以难度确实比较大。硬件资源应该问题不大，但是只是想用六面体完成，有这个要求
3.周期性网格我还了解的不太多，包括您后面说到的一些，愿闻其详
4.较早前曾使用gambit进行网格划分，但是网格量比较大，之前的计算网格量大概是是过亿了，gambit似乎当初处理不了这么多的网格量，所以才放弃的。
5.确实目前因为icem的网格思路与操作还比较陌生，整体还是比较困难

@东岳 多谢提示，我再研究研究。

我就是这个情况，笔记本会出错，win10系统win7都会出错，至今还没在台式机上发现