各位前辈们好：
我在YEqn.H 中定义了所有组分的一个反应速率

forAll(Y, i) { ROP_[i].field() = -(reaction->R(Y[i])()).source(); }

在求解时能正常输出：
QQ图片20220704180349.png

但是当我直接定义 H2的反应速率：

ROP_[H2].field() = -(reaction->R(Y[H2])()).source();

编译的时候就报错：

/home/zhu/OpenFOAM/OpenFOAM-2.4.0/src/OpenFOAM/lnInclude/GeometricField.C:188:1: note: candidate expects 4 arguments, 1 provided In file included from AMRreactingLMFoam.C:51:0: createFields.H:287:13: error: expected unqualified-id before ‘(’ token thermo. (reaction->R(Y[1])()).source()() ^ In file included from AMRreactingLMFoam.C:153:0: YEqn.H:30:14: error: ‘H2’ was not declared in this scope ROP_[H2].field() = -(reaction->R(Y[H2])()).source(); ^ In file included from AMRreactingLMFoam.C:54:0: /home/zhu/OpenFOAM/OpenFOAM-2.4.0/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable] scalar maxDeltaT =

显示H2没有定义，请问各位前辈在求解器中如何表示H2？

这个我还没有涉猎到，感觉很屌的样子:chouchou:

@fireztw 嗯，那看起来还需要两相流求解器。
不清楚，如果哪位大佬有类似的sci可以发来瞅瞅

@vbcwl 谢谢！学习一下，实现过程中有问题再向你请教

@李东岳 明白了，谢谢东岳老师:146:

请问各位老师openfoam中使用les方法计算槽道充分发展湍流使用那个lesmodel比较准确呢。

b3552765-7cd8-4842-bd50-b117acee71d9-image.png

我更新了一下，昨天那个花圈的有笔误

厉害厉害，SOWFA只支持openfoam2.4？

你那个图调节Q的大小，阈值变化一下。另外可能还需要计算一段时间。

没没没，我是说SOWFA这个软件挺好，SOWFA这个软件植入了一些新功能，风场的新功能，不是我自己植入的

@hy1112006 我有点记不清了，现在是可以打开的。可能是路径中的dyfluid-7 中的-不易识别？你可以试一试换一个路径编译，或者改成dyfluid_7试一试

@jack-qq 千万网格跑瞬态最好上好几百个核心 千万网格稳态好收敛的话 弄100 200个核也行

@沐风 您好，我最近也在研究可以计算凝固问题的模型，初步的模拟计划就是水滴在低温平板上逐渐凝固成冰，最好考虑水和冰的密度差异，请问楼主现在可以计算了吗？有点问题想请教一下，非常感谢！

各位前辈好，我现在需要计算H2的反应速率，我在 createFields.H 中进行了如下定义:

PtrList<volScalarField> ROP_(Y.size()); forAll(Y, i) { const word RRname = "ROP." + Y[i].name(); ROP_.set ( i, new volScalarField ( IOobject ( RRname, runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh, dimensionedScalar("zero", dimMass/dimVolume/dimTime, 0.0) ) ); }

随后我在 YEqn.H 中定义了ROP如何进行计算:

forAll(Y, i) { ROP_[i].field() = -(reaction->R(Y[i])()).source()/mesh.V(); }

他这样就输出了所有组分的一个反应速率，但是目前我仅仅想得到氢气的反应速率，并输出H2在每一个时间步上所有网格下数据的一个累加值，对此，我在 AMRreactingLMFoam.H 中这样写：

// ------------------------------------------------------------------------- // fileName myFile = "THRR.dat"; //OFstream myStream("THRR.dat"); std:: ofstream myStream(myFile.c_str(), ios_base::app); // ------------------------------------------------------------------------- // fileName File = "ROP.H2.dat"; //OFstream myStream("THRR.dat"); std:: ofstream Stream(File.c_str(), ios_base::app);

THRR是之前输出的总热释放率，在EEqn.H 中它是这样定义的：

// for THRR recording: THRR = 0.0; forAll(dQ, cellI) { THRR += dQ[cellI]; }

目前我遇到的主要是两个问题：

1.如何单独输出H2的反应速率，当我直接输出ROP.H2时显示H2没有定义

2.对每个时间步下所有网格的H2反应速率进行累加并输出应该如何写代码?

希望能得到大家的解答，谢谢！

@星星星星晴 最近可能工作会需要，又从尘封的硬盘里拷了出来，既熟悉又陌生，感觉连ParaView的用法都生疏了:136:

@bestucan 想偷懒来着

@chengan-wang 我用sprayfoam，你说的算的准不准完全要看你应用什么样的模型，你这问题让人摸不到头脑该怎么回答。。。楼上大佬说了，什么软件都能算得出来，你要是有闲心，matlab或者python都能自己写点东西出来算

@李东岳 好的，谢谢老师，我目前主要是想对多组分气体（氧气和氮气）的稳态流动扩散进行求解，想着基于twomixingLiquidfoam简单修改修改，因为我对MULES研究的不多，如果是稳态计算的话，我可否直接在相方程里面添加对流项fvm::div(phi,alpha1)，将求解alpha1的方程写成：
fvm::div(phi,alpha1) - fvm::laplacian
(
volScalarField("Dab", Dab + alphatab*turbulence->nut()),
alpha1
)
从而删去MULES部分呢？

各位大佬，请教一下，我如果要用双流体模型计算气固两相流，那是早期版本(of6等)里的twophaseEulerFoam还是最新版of9里的multiphaseEulerFoam合适？

alpha方程的对流项是通过MULES来算的。因此区分开了。

@bestucan 感谢前辈回复，已经解决了:146:

有个问题想要请教下各位老师，如图，twoliquidmixingFoam求解器里面这个扩散方程里面为啥没有对流项fvm::div(phi,alpha)？