目前计算一个喷雾燃烧场时，想改变下燃油喷射位置，因此把几何模型的喷嘴结构(5mm圆孔)向上游移动了10mm, map前在system文件夹下创建了mapFieldsDict文件，具体信息如图1所示；然后在对应时刻文件夹下建立了和原始计算场相同时刻的文件（共125个）。在map时使用了1）mapFields /OF6/P25；2）mapFields /OF6/P25 -sourceTime 0.116555 -mapMethod mapNearest -targetRegion region0两种方式均出现phi未能成功map,同时PaSR_kappa和PaSR_Qdot，sprayCloud类文件除了sprayCloudTheta文件，其余均未能map过来，如图2所示。另外在现有map结果基础上试算了下，log和err文件如图3所示，在读取化学反应机理时浮点溢出了，尝试了foamChemistryReader和new format均未能解决。请问这俩问题大家有什么方法解决么？万分感谢！
771f0b91-f3ee-41db-83b0-01decf0d55ac-图片.png
f30a09d4-20a4-45cc-9750-171e698d5778-图片.png

@fangyuanaza 我也遇到这个问题了，貌似是工具里面判定网格中心点距离的问题，太小的网格就提取不到值

你把斜坡做成STL，用snappyHexMesh抠，应该就行了

@myheart 懂的懂的。我这面也经常是，一个bug卡好几天，有个都是几个月，还有按年算的。。。就比如你这个回帖，3行搞定一个问题，其实背后肯定debug好久。辛苦辛苦。

@李东岳 我现在没有好的算法，知识发现零梯度这种边界条件是可以使得场内的浓度出去，不会在边界处积攒，如果编写算法，应该是和零梯度类似的可以。

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

微信图片_20220823200056.jpg 微信图片_20220823200145.jpg 微信图片_20220823200149.jpg
自适应网格达到设置最大的网格数后面的计算就不加密了，
但是我想实现的是网格在压力波传播的附件加密，压力波传播过去后加密过的网格恢复原始大小
请大佬指教是可以怎么解决
（我用的blastFoam求解器）

@李东岳 谢谢李老师！

@星星星星晴 非常感谢！我找找看

@星星星星晴 感谢大佬指点！

看起来像个bug。你可以试试双核的，然后手动把traslationRotation写进processor里面，然后看看报错不。不报错的话可能是bug

@李东岳 谢谢李老师，我再看看

@李东岳 :135: 不好意思，李老师，再打扰您一下，我还有一个问题，我目前使用的是v10版本的multiphaseEulerFoam，我在dragModel里可以通过interface_.dispersed().d()这种来调用直径参数，现在我想在求解器的.c程序里调用它，应该用什么方式呢
我注意到您再小代码里提到过以下代码

meanDiameter { type coded; libs ("libutilityFunctionObjects.so"); name error; codeExecute #{ const volScalarField& d = mesh().lookupObject<volScalarField>("d.alpha.oil"); scalar d32 = d.weightedAverage(mesh().V()).value(); if (Pstream::master()) { std::ofstream file; file.open ("d32", std::ofstream::out | std::ofstream::app); file << mesh().time().timeName() << " " << d32 << "\n"; file.close(); } #}; }

这段我是否能用呢

@李东岳 好的，谢谢老师，有错误信息就好办，看来还是得转战linux:135:

下面这个代码是对壁面一层网格进行赋值的。你可以参考一下。mag可以。

//- Initialise the near-wall P field to zero forAll(patches, patchi) { const fvPatch& currPatch = patches[patchi]; if (isA<wallFvPatch>(currPatch)) { forAll(currPatch, facei) { label faceCelli = currPatch.faceCells()[facei]; epsilon[faceCelli] = 0.0; G[faceCelli] = 0.0; } }

a9956a74-d12c-48a9-8f31-5da7020a5fc6-image.png 各位大佬下午好，我最近在跑一个微流道的圆柱绕流，网格尺寸非常小，然后计算的时间步长选取的也比较小，一个时间步长是5e-5，存储时间为0.5秒存储一次，当时在存储的时候会有小数点后几位出现，各位大佬有知道这是什么原因导致的吗:chouchou:

@李东岳 好的，谢谢:chouchou: 我已经找到对应的黏弹性算例了。

请问您最后是怎样计算磨擦速度的呢，我的算例计算磨擦速度也偏小一点。

@田畔的风 好的，谢谢。我是看到你在后面又把空气密度除了一次，所以存在疑惑。如果你取空气密度为1，那么在数值上应该是没问题的。

@李东岳 老师，求出的数值结果是6e9dde13-6ec5-40ac-98f1-7d6df59514f0-image.png ，还是2c1168c0-20c4-45ef-8e51-45a5f53de754-image.png