CFD中文网

@yhdthu 大佬，以下是问gpt的回复，不知道对不对：

这个编译错误提示显示在PtrList.C文件的第54行，表明在这一行代码中存在从Foam::Field<double>到Foam::CoeffField<double>的无效转换。 该转换无效的原因可能是因为ptrs_[i]的类型是Foam::CoeffField<double>*，但尝试将Foam::Field<double>类型的对象赋给它是不允许的。 要解决这个问题，您可以尝试以下几种方法： 1. 确保在您的代码中正确使用了正确的类型。检查PtrList.C中的指针类型和涉及到的其他类型是否匹配。 2. 查看代码中是否存在其他地方对ptrs_成员进行了赋值操作，确认这些赋值是否类型匹配。 3. 考虑修改代码以确保在赋值时进行正确的类型转换或者使用适当的类来替换。 您可以检查PtrList.C文件中的相关代码，特别是第54行的内容，以了解更多上下文并尝试解决这个问题。如果需要更多帮助，请提供更多的代码或上下文信息。

@yhdthu 大佬别闹，求赐教:140:

WeChat8db3756fb66113a48305c51893ed1215.png

模型还是很简单的，就是一个柱体+方形，但是我用icem画必须要用两个O grid，网格过渡很不好（至少我看起来很难受），不知有什么其他思路？

如题。

@xpqiu 好的，谢谢邱老师，我再试试😄

可以去看看funkySetFields这个工具

如题，目前的方法是blockMesh先生成方形单层mesh，再利用extrudeMesh生成360的3D网格。
但是出现个问题，里面的面不会自动的成为interior，不太清楚怎么处理。

以下是extrudeMeshDict

constructFrom patch; sourceCase "$FOAM_CASE"; sourcePatches (side1); exposedPatchName side2; extrudeModel sector; nLayers 10;//90; sectorCoeffs { point (0 0 0); axis (0 1 0); angle 20;//360; } flipNormals false; mergeFaces false;

写个类调用下fftw的函数接口就好了

看到个帖子不错，记录一下。关键词：householder transform

@yhdthu 在 如何在domain内加边界条件？ 中说：

域内面积为0的面，然后加边界不？

很抽象，很难想象如何处理的呢

@东岳 是的，这个operator[]不会每个时间步都去读，感觉这样改起来是个无底洞呢。不知是否能设置个贴体坐标，所有的点都相对其运动起来

ESI版本1906，有个bug，在compressibleInterFoam下，将turbulence models换成twoPhaseTransport时，除了laminar外都不能算，直接报错。bug提交在8个月之前，但是并没有人回应，不知其他人有没有研究。

@yhdthu thx

@东岳 如此说来，除了加密网格，没有好的解决方案啦:zoule: :zoule:

为社区做贡献😄

@李东岳
明白了，谢谢李老师。。我另外那问题，还要再麻烦您或其他老师帮指点下

@yhdthu 呵呵，OF的精度问题不止如此，湍流模型如果是后加的，算子分裂可能会导致精度下降到1阶和稳态无法达到的问题。所以Jasak在FE中玩耦合求解是另外的套路。

@yhdthu 前辈能够分享一下解决的经验吗？

@mohui 是的，interface处是有间断的，但是one-field形式的控制方程体现不出来间断
压力分布不对，看一下理论解就明白了～

@李东岳 前辈，我一直用pyFoam来观察残差，我猜的这个bounded variable plot出现是因为计算出现了不稳定，如果这个bounded量计算没过界就不会出现。:confused:

另外观测压力等物理量的残差如下，这个就比较好理解了:cheeky:

0_1508035857600_WechatIMG15956.jpeg