CFD中文网

我在对coalChemistryFoam改写，想把它和随机场方程进行结合。简单来说是对于连续相用一系列随机场进行计算，颗粒相不用（单独计算）。但这就遇到了一个问题，我结合着输运方程来说明：

for(label i=0; i<Y.size(); i++) { { volScalarField Yi = Ysf[nsf][i]; volScalarField Ymean = Y[i]; WienerTermYi = pow(2*gamma,0.5)*( fvc::grad(Yi) & dW_sf); fvScalarMatrix YiEqn ( fvm::ddt(rho, Yi) + mvConvection->fvmDiv(phi, Yi) - fvm::laplacian(turbulence->muEff()/Sct, Yi) == coalParcels.SYi(i, Ymean) + combustion->R(Ymean) + fvOptions(rho, Ymean) + (rho*WienerTermYi/deltaT) ); YiEqn.relax(); fvOptions.constrain(YiEqn); YiEqn.solve(mesh.solver("Yi")); YiEqn.solve(); fvOptions.correct(Yi); Yi.max(0.0); Ysf[nsf][i] = Yi; Yt += Yi; } }

在这个方程里，Yi指的是随机场的组分项（比如名为N2__01），Ymean指的是平均场的组分项（也就是组分场，名为N2）。颗粒的求解不考虑随机场，来自颗粒的源项被输入到每一个随机场中。但这里就会有一个问题，在计算时，会出现类似incompatible fields for operation [N2] + [N2__01]的报错，我的理解是因为这两个名称不同所导致的。这个问题要怎么去解决？我想的是，或许可以把这三项在方程外部导入到一个临时变量，再输入到方程里，但是不知道怎么用代码去实现。求各位大佬指点。

@李东岳 感谢李老师！原来没有特意写出来的都是同一个:zoule:

是想用cantera算一维火焰做FGM表，用在openfoam计算里ToT

已解决

或许是我的表述不够清晰，我想模拟图一所示的构型，两股不同的喷雾分布从不同的喷雾入口注入，我看到之前DieselFoam(或许是sprayFoam的前身)可以通过指定不同的X来喷注两种液体，这应该在sprayFoam中如何实现呢？

image.png

image.png

倒数第二个图，这种提示都看不出来具体咋回事
倒数第一个图，我挂libtorch的时候用过那个，你那个挂libtorch了么

厉害啊老铁

我的chemistryProperties是这样设置的：

chemistryType { chemistrySolver ode; chemistryThermo psi; } chemistry on; initialChemicalTimeStep 1e-8; EulerImplicitCoeffs { cTauChem 1; equilibriumRateLimiter off; } odeCoeffs { solver seulex; absTol 1e-12; relTol 1e-1; } // ************************************************************************* //

@欧阳 我以前 OpenFOAM-8 用的上面发的代码，好像没问题的话

paraview可以直接对分块的网格后处理

各位老师，之前用的是compressibleInterFoam求解器来计算横射流雾化，后面以期做到燃烧，但由于水平有限，在添加组分方程时便遇到许多问题。
现在发现reactingTwoPhaseEulerFoam，其本身就有YEqns.H，并且据我现在初步了解他可以通过选定reactingPhaseModel相模型来计算反应。

This model represents a phase with multiple species and volumetric reactions.

想问一下如果用reactingTwoPhaseEulerFoam可以做液体横射流雾化燃烧嘛？

目前想用一个类似sprayFoam的求解器，算一个旋流燃烧的算例。在喷入燃料前需要先通入空气，对冷态流场进行计算。但在冷态流场计算时，出现了部分高温、高压的小块出现，同时时间步变的很慢，无法继续正常计算，不知道各位有没有什么建议和指导:140:
18f84d67-316a-46dd-af4d-c6fcaeab642c-image.png
右边绿色区块为inlet，流入空气，仅有流向速度为20m/s（后方圆柱长度约为20mm）。紫色平面包括叶片为wall。
a72a1134-5beb-4e4c-aa34-a9f1909ab6c0-image.p ng
截取了某时间点的温度分布，出现明显高温小块（同时高压），且后续会持续上涨，无法正常计算。

各位老师，晚上好
我参考了$InterAdsFoam: An Open-Source CFD Model for Granular Media–Adsorption Systems with Dynamic Reaction Zones Subject to Uncontrolled Urban Water Fluxes$和$phaseScalarTransport.C$在compressibleInterFoam中添加了组分方程，个人感觉这个组分方程形式也不是非常复杂，但是在求解过程中会出现局部质量分数增加的情况，如下图所示，想请教一下大家有没有什么建议？

volScalarField Deff("Deff",turbulence.muEff()); volScalarField alpharho ( "alpharho", alpha2*rho2 ); surfaceScalarField alphaPhi ( "alphaPhi", linearInterpolate(alpha2*rho2*U) & mesh.Sf() ); fvScalarMatrix alpha2YEqn ( fvm::ddt(alpharho, Yi) + fvm::div(alphaPhi, Yi) - fvm::laplacian ( fvc::interpolate(Deff)*fvc::interpolate(alpha2), Yi ) == fvOptions(alpha2, rho2, Yi) - fvm::ddt(residualAlpha_*rho2, Yi) + fvc::ddt(residualAlpha_*rho2,Yi) ); alpha2YEqn.relax(); fvOptions.constrain(alpha2YEqn); alpha2YEqn.solve(mesh.solver("Yi")); fvOptions.correct(Yi); Yi.max(0.0); //update the results Yi = Yi*alpha2; // should use this guys as Yi Yt += Yi;

O2:
ede448a4-2b25-4d6a-99de-8f78100ea78b-O2.jpeg
alpha:
820eb593-7b64-4668-8589-a37ace97e464-alpha.jpeg
T:
24cd8ab9-c4ea-463c-b7c5-061e76f80f51-T.jpeg
U:ceee3886-0e70-4e37-9112-801c5f0fc52f-u.jpeg

@nbyjn 我扩大了沿射流方向的计算域，现在是已经可以算稳定了（但我不知道是什么原因）
这是600K的结果
f7ee0e69-84ec-4544-be67-a28423a738f8-image.png

@李东岳 是跑求解器的时候出错的

@Rachel0096 在 twoPhaseMixture类修改，编译通过，运行出错 中说：

还是卡在了下面的for循环当中

没瞧出有啥问题。看输出的信息就是遍历网格，计算每个网格内的密度倒数。
如果猜测是卡在for循环里了，这个cellVolMixture函数在heRhoThermo.C文件里也有调用，那里面也有一个for循环，瞅瞅是不是那里边有问题。

好的，李老师，我在回去多调试一下:zoule:

Cp_[cellI]输出一下

@杨英狄

pturb=0.39*（密度*湍动能）

类似这种，可以通过这个实现，不需要复杂的操作

const volScalarField& k = mesh_.lookupObject<volScalarField>("k");

tmp<fv::convectionScheme<scalar> > mvConvection ( fv::convectionScheme<scalar>::New ( mesh, fields, phi, mesh.divScheme("div(phi,Yi_h)") ) ); { reaction->correct(); dQ = reaction->dQ(); label inertIndex = -1; volScalarField Yt(0.0*Y[0]); forAll(Y, i) { if (Y[i].name() != inertSpecie) { volScalarField& Yi = Y[i]; fvScalarMatrix YiEqn ( fvm::ddt(rho, Yi) + mvConvection->fvmDiv(phi, Yi)

这个看起来并不会导致区别。有可能cfd-online的用户没做特别仔细的判定吧。