CFD中文网

Capture1.PNG
Capture2.PNG

如上2图，白色的点为测量数据点，每个点都有ux,uy,uz
先在inlet 上通过topoSet 专门设置出了一个小的区域，
图1是方形，图2是根据测量区域画出来的一个正方形加两个圆形的区域。

现在最大的问题是，你可以发现topoSet的边界特别明显。。
这种情况应该通过什么办法能解决呢？

为什么用topoSet的原因是如果直接在inlet 上设置为timeVaryingMappedfixedValue 这个边界条件的话，code 会直接插值到最边界。。结果看起来更有问题，topoSet这个解决方法一定程度上规避了interpolation到最边界的情况，但是会造成这个patch 和他原来关联的patch的边界问题。。

求大佬指点迷津

3d408cf3-836a-4aef-bb0b-a370d462f7e5-image.png 如图，左边的上部是大气边界，右边是封闭壁面。如果不在左边画和右边一样的细网格，得到的网格质量很差。如果画了细网格，则计算速度很慢。有更好的处理方式吗

@bestucan 还是那四个字：适者生存

打电话警告怕啥 先用再说

自己写的代码终于可以出结果了，但是用cgns存坐标数据的时候发现只能存一半，很奇怪，代码举例如下
三维的面网格，节点数2x11，因此是1x10个网格, 坐标储存在x,y,z数组中

real, dimesion(2,11) :: x,y,z CALL CG_BASE_WRITE_F(INDEX_FILE, BASENAME, 2, 3, INDEX_BASE, IER) ISIZE(1,1) = 2 ISIZE(2,1) = 11 ISIZE(1,2) = 1 ISIZE(2,2) = 10 ISIZE(:,3) = 0 CALL CG_ZONE_WRITE_F(INDEX_FILE, INDEX_BASE, ZONENAME, ISIZE, Structured, INDEX_ZONE, IER) CALL CG_COORD_WRITE_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, RealSingle, & 'CoordinateX', x, INDEX_COORD, IER) YZ坐标同上

读取时代码为

real :: node(2,11) node = 0.0 call CG_COORD_READ_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, 'CoordinateX', RealSingle, (/1,1/),(/2,11/),node,IER ) write(*,*) node

node的读取结果为

x(1,1) x(2,1) ... x(1,6) 0.0 0.0 ... 0.0

请问有没有大佬能指点一下

@evensun 求解器只负责求解数学式（控制方程），为什么要求解这个控制方程而不是别的控制方程是用户的选择，或者说取决于用户需要的是什么。
如果你需要百分百精确的描述一个非定常的物理问题，那你用稳态求解器原则上来讲是不合理的，但这个不合理不是说解方程这一步不对，而是你没有求解与你的需求相匹配的方程。如果你只需要部分描述这个物理问题，比如你关注时间平均量，那这时候用RANS去解这个问题，这一套下来是合理的。因为你需要的只是时均量，只是这个物理问题的一部分。
从这个角度去思考，你用稳态求解器计算了一个非定常的过程，求解器没有发散，甚至可能还体现出了一些有规律的震荡。这个现象本身其实是有意义的，但你在分析它的时候要非常小心。比如你认为“算得的结果呈现出了涡的发展变化过程”，这个表述是值得商榷的。

再啰嗦一点，关于RANS我一直觉得有个很有意思的角度。我们都知道雷诺数是决定流动的重要因素，我们可以认为控制方程也有一个“数值雷诺数”，它和普通的雷诺数定义类似，只是粘度项应该替换为分子粘度+数值粘度。在RANS中，数值粘度一般高出分子粘度2个量级甚至3个量级，所以尽管物理问题可能雷诺数很大，但数值雷诺数可能只在10的3次方或者4次方量级。10的3次方已经和我们通常认为的临界雷诺数接近了，导致RANS结果看起来是稳定的，就好像层流一样。
说白了，RANS模型对湍流的耗散非常强，把很多非定常的结构都给耗散掉了，留下了定常的部分。但这个耗散能力是有限度的，如果你求解的这个物理问题经过RANS模型后它的数值雷诺数正好处于临界雷诺数附近，那可能就是会有一些小的震荡。
这个湍流模型是自己开发的吗？你可以通过尝试不同的湍流粘度来验证一下，比如量级直接调大10倍，你说的这些震荡还存在吗？

请问有大佬知道如何在Fluent中利用UDF添加Maxwell速度滑移和温度跳跃边界条件吗

@hurricane007 这样做的确可以用。相当于使用linux的命令行，windows的UI。

有个问题想问问大伙，如题

@cccrrryyy 大神，想问问轴向尺寸指的是网格吗还是管道长度呀

周老师，这里我测试发现好像这种结构网格的加密方式对于三维模型可以，对于二维网格，在经过局部加密后再导入到OF里面是会报错的，错误是：--> FOAM FATAL ERROR:
face 11528151 does not belong to cell 210100. Face owner: -1 neighbour: 210099

From Foam::cellShape Foam::extrudedQuadCellShape(Foam::label, const labelList&, const faceList&, const labelList&, const labelList&, Foam::label, Foam::faceList&) in file extrudedQuadCellShape.C at line 105.

FOAM aborting
感觉好像是在z方向拉伸出来的面没有识别到具体的cell中

有大佬用cfx-pre算过二维旋转对称模型吗？我用二维旋转对称算的结果跟实验值相比差了蛮多，但用三维模型算的结果就对的上，但三维模型计算时间太长了，想用二维算。

@jinlinna 感谢！

实际工程上说的压力是什么压，压力表测到的是什么压，还有就是数值模拟计算得到的压力结果我们应该关注什么压呢？或者说从动压静压总压分布图上我们能看出什么东西。谢谢

之前接触过CFD++和STAR CCM+的噪声计算，里面都有专门的模块。openfoam这方面的功能没试过，不过搜到了一个第三方库，里面有源码和一篇文章，你可以参考下：
https://github.com/unicfdlab/libAcoustics
实在不行就把声源面的压力脉动都导出来，按定义写个程序去积分，虽然麻烦了点，但应该也是能做的

我设置了interface，check mesh的时候还是显示Unassigned interface zone detected for这是为什么呢？

@XIN 请问您解决了吗，我也遇到这个问题了

各位大佬，最近想用OF计算固定的多体装置，由于在大域内生成的结构过于粗糙 需在小域里生成 然后用mergeMesh合并到大域里去，但是小域的边界没办法删除
请问各位大佬，有没有什么设置建议？删除边界或指定特殊的边界？

@LuoS 在 snappyHexMesh对复杂地形添加边界层报“浮点数例外 (核心已转储)”错误 中说：

教程中似乎有个建筑风的案例：
/Tutorials/inco

好的谢谢 我试一下

@队长别开枪 谢谢老师指点，我好好研究一下