CFD中文网

最近涉及到一些化学反应，并不十分理解。在比如reactingTwoPhaseEulerFoam求解器的tutorial，bubbleColumnEvaporatingReacting中，展示了所有物质都在气相状态下的反应，CO^0.93 + H2O^0.24 = CO2^0.69 + H2^1。因为该反应被定义在reactions.gas文件中，反应物和产物的属性也在thermo.gas文件中，所有物质都为气态。

我想问问如果我的反应产物是有液态或者固态参与，比如A(气态) + B(气态) = C(气态) + D(液态)，应该如何设置这个reaction和thermo等文件。有没有相关的教程或者资料可以参考？

直接自己编程写。把每一个网格单元的面网格的值复制到另一套网格里面。

@李东岳 谢谢李老师。

这应该是TJunctionFan这个tutorial吧，这个确实能得到局部的速度增加，但是如果管道足够长的话，管道的入口和出口处仍然是相同的速度和流量。

回复您之前的问题，是inlet，outlet这两个地方。我想实现的结果就是在outlet的速度比inlet速度高，这确实是一个不符合物理的问题。。流量不守恒。没有外力的情况下是不可能实现的。硬要做的话，好像是可以通过添加momentum source实现，但是确实是一个非物理的问题，打扰了。。

@hotairballon1997 再coneinjection中 你可以通过控制flowrate来控制你喷射的parcel的速度

case ftFlowRateAndDischarge: { const scalar A = 0.25*pi*(sqr(dOuter_) - sqr(dInner_)); const scalar massFlowRate = this->massTotal()*flowRateProfile_.value(t)/this->volumeTotal(); const scalar Umag = massFlowRate/(parcel.rho()*Cd_.value(t)*A); parcel.U() = Umag*dirVec; break; }

应该影响不大。这种颗粒应该不多，不会所有的颗粒都恰好那样。所以对alpha的影响也很小。除非颗粒直径非常大。不过直径非常大的话，那就会引起收敛性的问题。

在WallSpringSliderDashpot.C中，有这样几句：

scalar etaN = alpha_*sqrt(p.mass()*kN)*pow025(normalOverlapMag); vector fN_PW = rHat_PW *(kN*pow(normalOverlapMag, b_) - etaN*(U_PW & rHat_PW)); ... scalar etaT = etaN; ... (if条件：颗粒不滑移的条件下) fT_PW = - kT*tangentialOverlap_PW - etaT*USlip_PW;

这里的颗粒碰撞中的切向阻尼系数，etaN，等于法向值，然后分别应用到求碰撞力的切向和法向分量，fN, fT，中。是不是意味着，我在颗粒配置文件中通过alpha设置的回弹系数值e, 就是同时设定了切向和法向的回弹系数，eT, eN, 都等于e？

请问这个理解对不对？

@李东岳 谢谢李老师，这个确实是我想找的功能。

不好意思才看到您的回复，不知道为啥论坛里的消息有时候没有提示，也没有邮件提示。。

@hotairballon1997
kinematicParcel中先看move部分，

初始化各个参数，计算lagrangian 步长。 更新参数，setCellValues，计算dispersion，然后calc方程 在calcVelocity中得到新parcel速度，新parcel速度calcVelocity部分中计算parcel受力得到新的parcel速度，返回到calc中，然后在calc中，根据 1-way 2-way 计算parcel是否对流场产生影响 然后返回到move中。。

https://openfoam.top/sprayVelocity/

所以Q1 在你不说你的模拟是1way 还是2way 是无法解释的
Q2 在上面的链接一定程度上可以解释一些

Explicit contribution, Su specified as a force Implicit coefficient, Sp specified as force/velocity

@samuel-tu 这个可以有。如果可以的话用结构网格是最方便的了，直接在该区域定义一个block，单独进行网格划分

我们当时也是添加各向异性，你可以把这个加到你的dispersion函数里面

const typename TrackCloudType::parcelType& p = static_cast<const typename TrackCloudType::parcelType&>(*this); const point start = p.position();

这个函数上面添加上p.position()然后传入进去，因为td没有position()函数，但是p有

就是 @浪迹天大 说的，一步一步弄，

@hotairballon1997 看着我真诚的眼神，告诉我，你看了我的帖子

or per-case override in the $FOAM_CASE/system/controlDict file

你懂我的意思吧.jpg

@李东岳 在 有没有utility或者force，能使颗粒在RANS模拟中像LES那样发散开来 中说：

我们之前详细的研究过这个东西 目前DPM这面的turbulent dispersion force结果都不太好 你可以尝试一下

想请教一下李老师，在这个stochasticDispersionRAS.C的代码中，dt代表的是什么，并没有看出它的出处在哪。。
从结果上来看，这个dt好像不是模拟中的timestep，它的值对于每个颗都不一样。

// Member Functions //- Update (disperse particles) virtual vector update ( const scalar dt, const label celli, const vector& U, const vector& Uc, vector& UTurb, scalar& tTurb );

codeFixedValue

一个网格一个粒子的话，$\alpha=\frac{v_i}{V_{cell}}$，一个网格多个粒子的话，$\alpha=\frac{\sum v_i}{V_{cell}}$，$v_i$是颗粒体积

在FLUENT里面用欧拉-拉格朗日求解时，欧拉方程里是认为α=1的，所以不适合用于高密度颗粒流的求解，而OpenFOAM中是考虑和计算体积分数α的。

理论上式对的，不过你要确认下在FLUENT里面用欧拉-拉格朗日求解时，欧拉方程里是认为α=1的？

@东岳 好的，谢谢老师～～