想问一下各位大佬有人做过利用PINN预测建筑群的钝体扰流吗？

我在对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]的报错，我的理解是因为这两个名称不同所导致的。这个问题要怎么去解决？我想的是，或许可以把这三项在方程外部导入到一个临时变量，再输入到方程里，但是不知道怎么用代码去实现。求各位大佬指点。

尝试用openfoam的coded自定义一个cellSet的组分源项，有大佬知道field设置成什么吗？

@李东岳 最近在将低版本OF的自定义边界条件移植到高版本OF的时候也发生了这个问题。串行求解没问题，并行求解报错。在不修改自定义边界条件源码的情况下，解决方案如下：1.在decomposePar前，将controlDict中引用的自定义库注释掉，例如

libs ( "libOpenFOAM.so" //"libUserBC.so" );

2.执行decomposePar，3.取消注释，4.运行求解器，5.重新注释，6.执行reconstructPar。
以上，便可解决该报错。

@wangfei9088 感谢大佬:146: :146:

@学流体的小明 谢谢您哈，能不能加您个联系方式交流一下

问题已解决，y站上没有更新教程，新版本的发行版openfoam会在指定facezone时自动生成主面和从面

@hachikoi 您好，请问您的问题现在解决了吗？我现在也在做和您相似的双浮体模拟，但是在模拟时出现了第一部就会发散的问题，重叠网格挖的洞数量很大。

然后又去pimpleFoam跑了一下，代码改成这样的：

momentumSource { type vectorCodedSource; active yes; name sourceTime; vectorCodedSourceCoeffs { selectionMode all; // cellZone pZone; fields (U); codeInclude #{ #}; codeCorrect #{ // Pout<< "**codeCorrect**" << endl; #}; codeAddSup #{ // Pout<< "**codeAddSup**" << endl; // const vectorField& C = mesh_.C(); const scalarField& V = mesh_.V(); vectorField& Usource = eqn.source(); const vectorField& U = mesh().lookupObject<volVectorField>("U"); // const scalarField& Rho = mesh().lookupObject<volScalarField>("rho"); const scalarField& nu = mesh().lookupObject<volScalarField>("nu"); // const scalarField& magU = mag(U); scalar A = 1e2; scalar B = 1e2; // vector C(0,1e4,0); forAll(V,i) { const scalar x = mesh_.C()[i][0]; const scalar y = mesh_.C()[i][1]; if(x < 0.5 && x > 0 && y < 0.5 && y > 0.45) { Usource[i] += (1e-5 * A + mag(U[i])* B * 0.5 ) * U[i]* V[i]; // Usource = (A * U[i] + B * mag(U[i]) * U[i]) * V[i]; // Usource[i] += - C * V[i]; } } // Info << "***codeAddSup***" << nl; #}; codeSetValue #{ // Pout<< "**codeSetValue**" << endl; #}; // Dummy entry. Make dependent on above to trigger recompilation code #{ $codeInclude $codeCorrect $codeAddSup $codeSetValue #}; } sourceTimeCoeffs { $vectorCodedSourceCoeffs; } }

跑完的速度场是这样：
pimple.gif

就完全处于A和B只能在1e2这个数量级，但凡再大一点就直接浮点溢出了。。。完全不知道是咋回事。。

@biaowai 在 使用mapFields 操作更改网格计算出错 中说：

印象可以修改差值方式或者用mapFieldsPar试试

我用的1706版本，尝试过修改差值方式还是不行。去试试mapFieldsPar。

@wjllll reactingFoam，用的是of的标准燃烧模型

@Amadeus openfoam本身有没有什么api或者怎么设置能够实现从体心速度插值到网格点上，我也在寻找实现方法。

@allanZhong 在某个方向还没有成熟的时候，总会有一些人顶着大佬的帽子来抨击方向的不成熟性，铸造自己的护城河，阻碍技术的发展。就像NLP、CV，发展初期一样有质疑、有诋毁，但是都一点点发展过来，颠覆了传统的认知，再到现在的大模型、deepseek甚至是AGI，都是颠覆传统认知的创新技术。
我们更应该充满信心，随着技术的发展，更多的研究成果的产出，终是会有颠覆性发展的一天的。

@cdc123 现在有看到论文里面会写，用数据辅助pinn训练，可以提升训练速度。

我用的python，用fulidFOAM库和自己写的代码，可以读openfoam网格，可以带编号。然后用pyg的graph存，也能保存网格与网格之间的连通属性。是可以 迁移到C++上面，达到想要的效果的。

@chenboyao 好的感谢大佬，目前我还没有加湍流模型，我先改一下时间格式看看能不能行

最近在用codedsource写个阻力源项，但是写完之后感觉效果都不太好，昨天用paraview看了一眼矢量图，发现用多孔介质模型的区域是没有速度矢量的，但是用codedsource的话从t=0开始就存在速度矢量，然后导致介质满天飞。所以openfoam是怎么实现让多孔介质模型明明作为一个动量源项但是不沾任何速度矢量的。。。

@liujm
问题已解决，添加libs("libfvMotionSolvers.so")，或者参考这里

@lxz-ii 请问一下，只有边界一层网格运动的状况解决了吗？我也遇到动网格文件配置后只有最边界的一层网格在运动，其他网格不动

@李东岳 好的，谢谢老师