CFD中文网

pressureGradient { alpha alpha; }

类似的，试试

@Dingcy 在 如何输出出口escape的颗粒数量随时间变化的曲线？ 中说：

如何输出出口escape的颗粒数量随时间变化的曲线？之前算单向流的时候我输出出口浓度随时间变化曲线是通过paraview那个plot selection over time, 但是现在这个颗粒好像不能那么整了....:136:
1c05a757-a68a-499f-bfc9-1582a3c49a26-image.png
前辈您好，请问这个操作怎么实现，得到下图:134:

@李东岳 好的好的，谢谢大家

应该是drag太大了，颗粒没弹起来

自问自答，我粗略的整出来了，如图1，要注意画红圈地方的选择:chitang:
746340cd-5bc6-47e8-a92f-b507dd49d9a8-image.png
图1

下面这个代码是对壁面一层网格进行赋值的。你可以参考一下。mag可以。

//- Initialise the near-wall P field to zero forAll(patches, patchi) { const fvPatch& currPatch = patches[patchi]; if (isA<wallFvPatch>(currPatch)) { forAll(currPatch, facei) { label faceCelli = currPatch.faceCells()[facei]; epsilon[faceCelli] = 0.0; G[faceCelli] = 0.0; } }

@李东岳 我现在没有好的算法，知识发现零梯度这种边界条件是可以使得场内的浓度出去，不会在边界处积攒，如果编写算法，应该是和零梯度类似的可以。

LES只能计算瞬态。在湍流计算中，RANS方法全部模化湍流项；LES方法则模化小尺度涡部分，解析大尺度涡部分， 所谓“解析”就是瞬态部分。LES的强制稳态没有意义，因为只计算了很少一部分的模化湍流，而大部分湍流（解析）并没有计算。因此LES的强制稳态结果，可能更接近于强制的层流稳态结果。

patchInternalField表示你说的这些量

@李东岳 是个稳态计算过程，不随时间变化，但是公式中r值更新会影响边界的U的值

@hotairballon1997

在您回帖的链接中好像没有看到关于这个问题的解答，

对呀，要分类讨论，就像你说的，是用 PBM，还是拉格朗日。
我发的那个链接里 MPPIC 和 DPMFoam 中用到的粒子都是粒子包（particle parcel），这种用粒子包去代替粒子的方法我没了解过，但肯定已经不是物理粒子了。所以要看哪个层面，哪种模型的碰撞聚合了。

@dingcy 自问自答了，在tecplot里输出非压缩的视频文件，然后放到ps里，增加分辨率，然后再渲染视频，最后的视频就清楚啦。

@李东岳 好滴好滴，非常感谢，我这就试试:xinxin:

@李东岳 李老师您好！我在用pisoFoam RANS kOmegaSST模型计算圆柱绕流（Re ~10^4）遇到了计算误差大的问题（C_d 误差约10%，St 误差约为30%）。通过阅读文献发现大家对入口k和omega的公式使用出入很大，一时间不确定计算误差大是边界条件问题还是离散和求解的问题。请李老师指教！

@李东岳 我发了帖子之后才发现你之前发过通知了，不是催更，抱歉抱歉~:chitang:

解决了，把1换成1.0就好了，感谢 @麦迪文 师兄的帮助。

@东岳 光看流场的矢量图可能确实看着还凑合对称，但是假如我要在出口得到RTD曲线（停留时间分布曲线），就会得到对称出口明显不同的RTD曲线？

你这个结果看起来很正常，应该没有调用turbulent dispersion force， https://pdfs.semanticscholar.org/ec18/d3e3140efbf5651c07c7e8635402777f7e3b.pdf 我在这里也模拟了这个算例，需要注意的就是调用湍流分散力，

@东岳 ok:huahua:

按照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给注释掉了，我这种做法对吗？