CFD中文网

大佬问题解决了吗，看你贴的报错信息貌似是代码里面这段直接”=“出错了？
/home/dong/foam/foam-extend-5.0/src/foam/lnInclude/PtrList.C:54:18: error: invalid conversion from ‘Foam::Field<double>’ to ‘Foam::CoeffField<double>’ [-fpermissive]
54 | ptrs_[i] = (a[i]).clone().ptr();

@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