东岳流体虚拟机4月更新, 配套ISO镜像同样上传至百度云. 安装方法不再赘述.

@五好青年 大佬，后来者大胆问下这个有解决的好办法或者好思路了么

@李东岳 是的。我又从头到尾排了一边，好像是求解马赫数的方式有点问题，现在还在验证，希望是这个原因

@李东岳 感谢李老师，我准备转向org版本进行研究，采用全新求解器重新实操一下

求助，目前在使用hy2Foam求解器计算类似乘波体外形，已经计算过钝头算例可以得到收敛解。但目前计算外形相对复杂，二维计算，计算域为1*0.55m，网格数量42000。计算过程中经常报错，求问有没有什么调试建议。下面是报错信息。

Mean and max Courant Numbers = 5.29417e-05 0.0491796 deltaT = 3.99074e-11 Time = 2.017990284e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N2, Initial residual = 0.000161695, Final residual = 3.0494e-13, No Iterations 1 DILUPBiCG: Solving for O2, Initial residual = 0.000148597, Final residual = 1.8029e-13, No Iterations 1 DILUPBiCG: Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for Ux, Initial residual = 1.07385e-06, Final residual = 1.49309e-06, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 8.83461e-07, Final residual = 1.19425e-06, No Iterations 1000 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 2.94288e-05, Final residual = 3.15326e+222, No Iterations 400 Attempt to use rho2ReactionThermo out of temperature range 11152 times during this iteration. -1.64159e+236 < Tlow: 200 Phase no 1.0 ExecutionTime = 3564.14 s ClockTime = 3564 s Iteration no 4641 (0.79 s) Mean and max Courant Numbers = 5.41311e-05 0.0489379 deltaT = 4.07806e-11 Time = 2.018031064e-06 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N2, Initial residual = 0.000162296, Final residual = 1.67236e-16, No Iterations 2 DILUPBiCG: Solving for O2, Initial residual = 0.000148383, Final residual = 6.70388e-18, No Iterations 2 DILUPBiCG: Solving for NO, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for N, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for O, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #6 Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:? #8 Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at ??:? #9 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #10 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #11 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam" #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #13 ? in "/home/sy/OpenFOAM/sy-v1706/platforms/linux64GccDPInt32Opt/bin/hy2Foam"

这两个不能叠加使用。

@杜浚平 十分感谢！

大家好！
目前我做一个几何对称的共轭传热模拟，所以我只建模一半后面采用mirrorMesh做镜像，但是遇到问题
bolckMesh→mirrorMesh→splitMeshRegions有问题，划分不了对称后的区域；
bolckMesh→splitMeshRegions→mirrorMesh也有问题，区域对称过来了，但是没有划分不同区域；
具体应该有说明操作，有人有相关经验吗？感谢！

我看你高亮的确实是contour。
paraview我也不会太会高级的操作..

感谢老师们！这里确实是多变量的计算，类似于把单个组分场拆分为很多不同的场来计算，但是输入相同的源项，所以会有两种不同的Y在同一个方程里。程序就是在方程这一步出现了错误，我试试各位老师的方法，多输出一些内容看看哪里出问题了。

@李东岳 还有个问题就是，粒子有时候会分散，有时候就是一坨。李老师有遇到过吗
f622a0de-4015-453e-93fc-2b952bdd661d-image.png 79a28b41-4d27-4440-84dc-6ccfd84fdf7c-image.png 5d210d73-a556-4130-85d7-016a7bc98fb3-image.png

@李东岳 LES和laminar都试过，只要把空化核数降一些就不会发散。

2025年3月，岳子回国教学CFD编程，现场让学生植入湍流模型，方程在黑板一写，学生植入的代码各种不一样，各种bug就出来了，当然岳子也是手到擒来，各种拿捏。

但是有一个问题属实卡住了。学生发现，PDE的求解顺序不一样，导致结果大不相同。

首先这是一个稳态的算法，整体是一个大迭代的过程，因此PDE方程求解顺序会影响结果，但不应该有这么大的影响。然而下面两个图，是k方程与epsilon方程求解顺序互相调换后的的结果。

因为是自己写的代码，所以不方便别人debug。感兴趣可以把OpenFOAM里面的kEpsilon模型的k方程、epsilon方程求解调换一下顺序，会发现原本250步收敛的pitzDaily，现在需要计算2000步（但结果基本差不多）。

因为稳态问题或许会收到求解顺序的影响，但不应该这么大（1、原生kEpsilon模型250步收敛的更改后2000步没收敛；2、或者我们自定义代码的nut的特别大的区别）。

目前我已经找到了解决方法（先求解k且不影响结果），但是原因并不能自圆其说。感兴趣的老铁，欢迎讨论。

捕获.JPG

微信图片_20250419212453.jpg

感谢分享！

谢谢！。

这是算的什么东西？好圆的球。。。

异常一直存在么？会消失么

是个bug，你直接把loss_dataP删了就行，这个应该是我当时要增加数据点用的。目前没有作用了。loss_data可以删了，目前跟loss_pde重复

怎么样，这部分玩的怎么样，有什么感受。PINN跟FVM是不是大不一样

@学流体的小明 在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题 中说：

在 OpenFOAM中partialSlip边界条件中valueFraction的设置问题 中说：
我了解到valueFraction 的type是field 类型，我不知道是否可以定义个类型也为field 类型的
𝛼
，然后赋值给它？

这个思路可以的。
wallShearStress，OpenFOAM有个后处理函数可以生成wallShearStress，它就是存在于patch上的。你可以遍历wall边界的面，先从wall边界上读取wallShearStress，然后有函数找这个面对应的内部网格，也就是你壁面上第一层的网格，这样给$\alpha$幅值就可以了。或者直接自己计算这样可以从wall边界上读取wallShearStress，算一个偏导数再乘以动力粘度就可以了。

！感谢大佬，关注你好久了
OK，那我先试试。

看了各位大佬的评论，我觉得无论是什么样的矩形通道，应该都是用的最小那一边的半槽高度。例如进出口内外均为周期性的平板槽道，它就等价于一个AR=宽/高=无穷的矩形通道，这种情况不可能用水力直径的，应该都是用的最小半槽高度

@李东岳 老师我对这个案例比较感兴趣，也想模拟试一下，谢谢老师