然后又去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 请问一下，只有边界一层网格运动的状况解决了吗？我也遇到动网格文件配置后只有最边界的一层网格在运动，其他网格不动

@李东岳 好的，谢谢老师

@liujm 在 为什么我在用致动线模型算的时候得到的涡量图和文献中不一样 中说：

@卡洛 卡老师，能分享一下边界条件的修改思路吗:140:

怎么个不一样

有没有大佬有openfoam各类文档的总结，我可以用fastgpt或者dify做一个知识库，然后我提供api（用Gemini2.0 flash模型）来让大家使用。

有同学想组队团报课程有优惠可以享受，可以联系vx:wwjc1230_，一块报名

各位老师同学大家好！请问有无六自由度求解器开发的经验，我目前是想开发基于foam extend4.0的浸没边界法下的六自由度求解器，主要用来定义并求解物体从高处自由落体入水运动过程，不知道从何处下手，希望各位老师同学不吝赐教！感激不尽！

使用hy2Foam求解器做RamC模型验证，在驻点线画图时发现异常，组分质量分数结果使用线性坐标为正常曲线，如下所示。6f6a06b3-134f-4d53-be6f-85c18312c808-image.png
但当转换为指数坐标时，曲线出现奇怪波动，请问是为什么呢？7e641883-94cd-4d05-bf9f-0405d67feb7d-image.png

可以的

@coolhhh 厉害，测试代码copy了。不过，这里直接操作就又可以了:wocao:

OpenFOAM-v2206/tutorials/incompressible/pimpleFoam/LES/planeChannel/setups.orig/common/system/controlDict，有一个计算Cf场的函数，可以参考下：

Cf { type coded; libs (utilityFunctionObjects); name Cf; timeStart $/timeStart; writeControl writeTime; codeExecute #{ auto* CfPtr = mesh().getObjectPtr<volScalarField>("Cf"); if (!CfPtr) { Info<< "Create skin-friction coefficient field" << nl; CfPtr = new volScalarField ( IOobject ( "Cf", mesh().time().timeName(), mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh(), dimless ); regIOobject::store(CfPtr); } auto& Cf = *CfPtr; Info<< "Computing skin-friction coefficient field\n" << endl; const auto& tau = mesh().lookupObject<volVectorField>("wallShearStress"); const dimensionedScalar Ubulk(dimVelocity, 17.55); Cf = mag(tau.component(0))/(0.5*sqr(Ubulk)); #}; }