@wsy11 好的，稍等，我提交了

破案了。程序没问题，是我自定义的用于AMR的coded量忘记clear了

@shrine 是只在出口下降了吗，你把出口压强改成零梯度试一下，可以看一下案例里面那个open channel的设置

@Hope 楼主有没有ESDU的document呀，风剖速度和湍流强度以及风谱，没订阅，也搜索不到

各位前辈好，请问我使用toposet创建的新的边界或者面之后，我应该怎么在这个边界或者说面上加上边界层呢，是不是相加边界层只能用snappyHexMesh

@学流体的小明 恩呢，谢谢您的回答。
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: v2212 |
| \ / A nd | Website: www.openfoam.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh dynamicRefineFvMesh;

// How often to refine
refineInterval 1;

// Field to be refinement on
field alpha.water;

// Refine field inbetween lower..upper
lowerRefineLevel 0.001;
upperRefineLevel 0.999;

// If value < unrefineLevel unrefine
unrefineLevel 1;

// Have slower than 2:1 refinement
nBufferLayers 6;

// Refine cells only up to maxRefinement levels
maxRefinement 1;

// Stop refinement if maxCells reached
maxCells 15000000;

// Flux field and corresponding velocity field. Fluxes on changed
// faces get recalculated by interpolating the velocity. Use 'none'
// on surfaceScalarFields that do not need to be reinterpolated.
correctFluxes
(
(phi none)
(nHatf none)
(rhoPhi none)
(alphaPhi_ none)
(ghf none)
(phi0 none)
(dVf_ none)
(alphaPhi0.water none)
(alphaPhiUn none)
);

// Write the refinement level as a volScalarField
dumpLevel true;

// ************************************************************************* //
这是我的dynamicMeshDict文件。
我按您说的方法去试一下。

@Prometheus10
我看了一下那个帖子，可能因为实现方法的差异，两种边界会有细微的差别，但应该不会这么大。

你的其他边界是怎么设置的？symmetry不用特别设置，每一个变量都是sym就行。但是，slip你是具体怎么设置的？速度肯定是slip，其他边界是什么？如果搭配不合适的话，和sym不一致是很正常的事情。

再一个，是不是因为雷诺数特别小，以致实现方法带来的微小区别被放大了？或者说，在低雷诺数下，实现方法带来的差异特别大？建议可以加点流速试试，看一下高流速下是不是还差这么多？

看起来好像你得OpenFOAM没编译完整

wsl2中安装ubuntu18.04后只有64个CPU分配，按照教程新建并修改./wslconfigure文件后，关闭wsl并重启，仍然还是64个核，请问有人遇到到这种情况吗，谢谢！！!图片.png

@coolhhh
大佬您好，我可能忘了介绍我的研究时段了，我研究的是较长时间的模拟，理想模拟时间应该是2-3天；模拟区域也比较大，大概为长宽高为6km6km600m。
1、因为研究的时间较长，所以每个边界的风速肯定是变化的（风速输入用的是每个边界所设置的非均匀分布风塔的观测数据，每个面都有数个观测点，观测数据风速利用timeVaryingMappedFixedValue每隔十分钟输入）因为之前学习了windaroundbuildings这个算例（也是建筑扰流），所以设置求解器为simpleFoam，这个算例风速确实是时不变的，这一点没有考虑到，等换成瞬时求解器试试；
2、目前还没有试过五个入流面，因为我也认为可能不太符合通量守恒，所以先设置了3个速度入口和2个出口，我先做一下大佬在4中提供的建议；
3、LES模拟也考虑过，正在学习，等做完2中试验后尝试一下大佬提供的条件，就不用计算k和epsilon文件的设置了。

在ParaView 中用 Temporal Particles To Pathlines 追踪粒子轨迹，会出现多余的线条是怎么回事，求指导。

abc25b9f-320d-4e32-8d77-4c45ac5ea775-image.png

@coolhhh 我明白了，您不妨试试uniformFixedValue。另外dsmcFoam+的dsmcCloud.C中没有correctBoundaryConditions，我估计是包装在别的函数里了。

在进行模拟时，运行136s出现上述报错：Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
2024-03-15_21-17.png

附上我的算例边界条件，是一个基于真实地形，使用codefixedVal函数实现非均匀边界条件的模拟。1.zip
请哪位大神给点建议，谢谢

TaylorBubble (1) )

还有一个问题就是我在chemkin反应机理中添加了这个反应：
1.png
对于hv这个光子应该在哪里定义呀？我在chem.inp以及therm.date文件中没有定义光子，也不清楚如何定义，但是采用chemkinToFoam时显示了如下错误：
屏幕截图 2024-03-15 190833.png
hv改成大写以后就成了这样：
1.png
恳请各位老师，大佬解惑，感谢各位了！

:135: :135: :135:

@Prometheus10 对 如果这样理解的话 为什么symmetry计算结果和理论解更相近呢

@tidedrinker @李东岳 找到原因了，是因为我输入命令的问题，我之前输入的是postProcess，应该用 interFoam -postProcess ， 。。。。新手不太会

@李东岳 李老师，最近在做实验，昨天又试验了一下，p场稳定了一些但是还在波动，但是U.air场变得混乱了，李老师有空可帮我看下吗，还有一个问题719e060c-5757-402d-b40a-0d978c2be0a7-tp_of10.zip
b39e6d46-adc8-4f48-b766-281d3efcc0b7-image.png

DPMFoam中的UcEqn.H:

fvVectorMatrix UcEqn
(
fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)

fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc) continuousPhaseTurbulence->divDevRhoReff(Uc)
==
(1.0/rhoc)*cloudSU
);

UcEqn.relax();

volScalarField rAUc(1.0/UcEqn.A());
surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc));

surfaceScalarField phicForces
(
fvc::flux(rAUccloudVolSUSu/rhoc) + rAUcf(g & mesh.Sf())
);

if (pimple.momentumPredictor())
{
solve
(
UcEqn
==
fvc::reconstruct
(
phicForces/rAUcf - fvc::snGrad(p)*mesh.magSf()
)
);
}

///////还有就是这一项代表了什么意思呢////////
surfaceScalarField phicForces
(
fvc::flux(rAUccloudVolSUSu/rhoc) + rAUcf(g & mesh.Sf())
);

@LiuWenchao 好嘞，我去看看，谢谢大佬！

谢谢李老师，那我是不是直接用turbulence->nut()就可以获得混合湍流运动粘度了？我是看到createFields.H定义了这样的指针，对应找到他的 incompressible::turbulenceModel 类，就发现mut和nut是一样的，对于这种多相流的，它实际用的是哪个类，在哪里体现了密度的影响？对openFOAM不太了解，谢谢老师指点
image.png

您好，我在修改了入射条件后依然没有解决这个问题，不知道您有什么解决办法吗

https://github.com/AndreWeiner/of_pytorch_docker
分享一个教程！

@xzdbjx 大佬，可以提供一下编译进Fluent里面的udf吗？ CFD小白，自己写的总是报错:zoule:

@李东岳 多谢李老师，解析写的很清晰，基本弄明白了

没事，找到错误就好。这个方法其实比较笨，因为每次你修改网格，就要重新处理一遍数据。
有的时候会麻烦一些。

祝收敛

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

@上级 大佬，想问下cfMesh上手难吗？我在做相似的研究，想转成cfMesh试一下

@学流体的小明 我也觉得难搞，我重组了某个时间步的网格文件，发现只有一个point文件，可能要一个点一个点的加😫，而且还有各种场文件貌似也要改😂

好的，谢谢李老师！！ 这个问题我纠结了好几天

@李东岳 明白了，谢谢!

@crispy 那你看其他库朗数的结果怎么样，1以下的，是时间步长越大，cd越好么

@李东岳 李老师，我发现应该是dev2(T(……)）这一小段出的问题，删掉这个转置操作就能编译成功。因为我是在传热求解器的基础上进行修改的，T在createFields.H里声明了是温度场，所以是不是程序将T识别成了温度而进行不了转置操作了呢？

请教各位大佬一个边界条件的选择问题
我的研究区域是一个三维区域，长宽高分别为（6000m*6000m*600m），如下图（俯视图）所示：
08ae4b29-0f40-4ce5-b8b2-db687c3fc3fb-image.png

图中的红色点是在研究区域内的观测站点，每一个观测点位可以得到从地面到600m高度处的风速、气压、温度等数据，
请问应该选择什么边界条件来把这些观测点位的数据输入到openfoam中？

@WYing 在 fvSolution中PIMPLE算法的参数设置问题 中说：

@李东岳 谢谢李老师的回复！抱歉我的上个残差图绘制有问题，正确的应该是这样

Residuals_Ur5.png

从上到下依次是速度残差，压力残差和湍流统计量残差。

上面那个code是我随意找的，我目前使用的code是：

PIMPLE { correctPhi yes; nOuterCorrectors 2; nCorrectors 2; nNonOrthogonalCorrectors 1; turbOnFinalIterOnly false; } relaxationFactors { fields { p 0.3; } equations { "(U|k|omega)" 0.7; "(U|k|omega)Final" 1.0; } }

我主要做了以下修改：
（1）调整 nOuterCorrectors=2，结果几乎不变
（2）turbOnFinalIterOnly=false，即每次PIMPLE loop都计算湍流，结果几乎不变
（3）去掉全部的松弛因子，结果几乎不变。

在模拟中，CFL保持在0.5-1.5之间，但是和实验结果对的不好。李老师觉得还有可能是哪里出了问题呢？

请问同学，你的问题解决了吗

@李东岳 嗯对，感谢李老师！

2、经验教训
不知道说的对不对，前辈们可以建议一下，在新编译软件的时候需要注意，尤其是新加PATH的时候，编译后可能并不会直接出现问题，但是重启电脑将原有.bashrc重新source会引起错误

!感谢大神的回答！

@Shihang-Chen 哇！感谢，确实解决了！

有数据的PINN效果还行，但是完全的无数据的不太行

我不太懂这个。确认一下：你说的这个意思，是说PINN需要CFD算出多套参数的结果，然后这些多套参数和结果就是数据。比如先用CFD算100套进口条件，得出100套结果，把这100套进口参数+结果给到PINN，这样的话效果还行？

绝对到位！嘎嘎一嘎子！

另一种方法：仅添加到UEqn当中，压力方程不做处理，斥力的贡献依靠UEqn.H参与到后续的计算当中。结果就没有压力振荡了。

:146: :146: :146: :146:

Interesting！

好的老师，我算完给您回复。

谢谢岳哥，成了，牛牛牛

只要有浮力就会发散，为什么呢

好的，谢谢李老师，我试一下

0时间步的U，是用湍流生成方法直接生成的，可以计算得到div(U)。然后就可以计算div(phi)： U -> phi = fvc::flux(U) -> div(phi)。测试结果是，这种方式计算的div(U)和div(phi)结果基本一样的，区别是小数点后五六位。此时div(U)=div(phi)=0

那这个算法很屌。不错。

@李东岳 在 CodedFixedValue边界设置问题 中说：

http://dyfluid.com/code.html 来这查一下

好的，谢谢东岳哥

@Calf-Z-DNS 在 用snappyHexMeshDict生成周期性边界问题 中说：

@东岳 谢谢东岳老师
我在cfd-online上面找到一种方法：
下面是原文：
“I encountered the same problem and I could solve it by defining rotational periodicity and then by declaring the corresponding vertices in ICEM as periodic. Then Right Click on Faces in the Blocking tree -> Periodic Faces -> ensure that really every periodic face is colored accordingly. createPatch went fine after this procedure.” ---by Fabian E.
然后我勾选了ICEM里面blocking里的faces的Periodic，再导入openfoam，发现就不会出现上述错误。
但我不太清楚为何这样会消除错误？

请问同学可以告知是怎么设置的吗，icem里面有face pericdic吗？我这里只设置了点对称

velocity.,txt
(
(0 (4 0 0))
(0.5 (4 0.0698 0))
(1 (4 0.0698 0))
(1.5 (4 0.0698 0))
(2 (4 0.0698 0))
(2.5 (4 0.0698 0))
(3 (4 0.0698 0))
(3.5 (4 0.0698 0))
(4 (4 0.1397 0))
(4.5 (4 0.2096 0))
(5 (4 0.2797 0))
(5.5 (4 0.3499 0))
(6 (4 0.4204 0))
(6.5 (4 0.4911 0))
(7 (4 0.5622 0))
(7.5 (4 0.6335 0))
(8 (4 0.7053 0))
(8.5 (4 0.7775 0))
(9 (4 0.8502 0))
(9.5 (4 0.9235 0))
(10 (4 0.9973 0))
(10.5 (4 1.0718 0))
(11 (4 0.9973 0))
(11.5 (4 0.9235 0))
(12 (4 0.8502 0))
(12.5 (4 0.7775 0))
(13 (4 0.7053 0))
(13.5 (4 0.6335 0))
(14 (4 0.5622 0))
(14.5 (4 0.4911 0))
(15 (4 0.4204 0))
(15.5 (4 0.3499 0))
(16 (4 0.2797 0))
(16.5 (4 0.2096 0))
(17 (4 0.1397 0))
(17.5 (4 0.0698 0))

)
入口设置
{
type uniformFixedValue;
uniformValue tableFile;
uniformValueCoeffs
{
file "velocity.txt";
outOfBounds repeat;
}
value $internalField;
}