CFD中文网

@Banbor 现在的网格水平停留在上面图示水平，没有加入建筑物。教程中似乎有个建筑风的案例：
/Tutorials/incompressible/simpleFoam/windAroundBuildings
要不试着和simpleFoam/turbineSiting的地形网格设置拼一下？

@李东岳 code是个好东西啊，得想想怎么拟合曲线咯，谢谢李老师的参考案例啦

function1类型的描述中有这样一句话：

(https://www.openfoam.com/documentation/guides/latest/api/classFoam_1_1Function1.html)
恰好我需要根据一组离散数据点拟合曲线，并将其结果赋值给一个function1对象的变量。
问题：对function1对象如何进行描述中所述函数插值功能？有无参考案例？多谢！

@李东岳 好嘞 ，谢谢李老师

@李东岳 比较显著的是一个角落上没有生成边界层，这个角的地形向下走得比较急，而且我的边界层是绝对高度的，不知道是不是因为这两个原因导致边界层生成时出了边界而出错。
在后续的相对高度边界层添加时，FeatureAngle也无法增加至180，若设为120，则在边界地形凹陷处未能添加上边界层。

工程做法参考就是扩大边界范围，或者后续想办法做边界地形面平滑处理。
在一篇论文中，extrudeMesh方法生成的网格，顶面和地面形状一样，感觉不太好。
后续如果snappyHexMesh没法做得更好，就只好先用这套网格将就将就，再研究别人用blockMesh改的网格方法，或者研究一下将fluent meshing的网格（边界层+多面体的网格质量不错）转过来。

未生成边界层的角：


其他角边界层：

在用snappyHexMesh对复杂地形生成网格时发现下述问题，未能找到根除办法，还请有碰到过同样问题的朋友讨论讨论：

1. 报错信息：

Determining displacement for added points according to pointNormal ...
Detected 0 points with point normal pointing through faces.
Reset displacement at 0 points to average of surrounding points.

displacementMedialAxis : Smoothing using Medial Axis ...
#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::medialAxisMeshMover::calculateDisplacement(Foam::dictionary const&, Foam::Field<double> const&, Foam::List<Foam::snappyLayerDriver::extrudeMode>&, Foam::Field<Foam::Vector<double> >&) at ??:?
#4  Foam::medialAxisMeshMover::move(Foam::dictionary const&, int, Foam::List<int>&) at ??:?
#5  Foam::snappyLayerDriver::addLayers(Foam::layerParameters const&, Foam::dictionary const&, Foam::List<int> const&, int, Foam::decompositionMethod&, Foam::fvMeshDistribute&) at ??:?
#6  Foam::snappyLayerDriver::doLayers(Foam::dictionary const&, Foam::dictionary const&, Foam::layerParameters const&, bool, Foam::decompositionMethod&, Foam::fvMeshDistribute&) at ??:?
#7  ? in "/home/lsfluid/OpenFOAM/OpenFOAM-10/platforms/linux64GccDPInt32Opt/bin/snappyHexMesh"
#8  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? in "/home/lsfluid/OpenFOAM/OpenFOAM-10/platforms/linux64GccDPInt32Opt/bin/snappyHexMesh"
浮点数例外 (核心已转储)

2.处理过程：
更改snappyHexMeshDict中addLayersControls参数featureAngle从180到120，网格如下所示。大部分网格生成边界层，一个角落未生成边界层，若继续增加featureAngle将报上述错误。
3.期望：
* 知道上述错误缘由和其他处理方法；
* 希望能增大featureAngle，使地形曲面能全部添加边界层。

填坑，参考官方t13教程（没记错的话，反正题目有个stl，似乎只有这种方法）

经Gmsh论坛解释，Gmsh不支持该功能，建议用meshlab等软件先进行几何前处理。
弃坑，走人

Gmsh导入的stl文件为网格文件，经gmsh.model.mesh.createGeometry()可以转换为Gmsh的内置几何形式。目前需要对stl面进行切割，但Gmsh仅支持occ内核的布尔运算。所以想要请教一下有什么办法将Gmsh的内置几何转为occ几何，谢谢！

根据t13.geo文件，在创建几何前必须依据角度对离散面进行分类，但我的stl文件近似平面，使用面分类，将角度设置为0，会导致内存暴涨，半天没有反应。所以请教一下各位大神：
1.在不使用面分类的前提下如何将近似平面的stl文件处理成几何曲面？（stl导入gmsh显示为网格，非实体）
2.若必须用面分类，有什么简单可行的操作方法？
谢谢！

填坑：
疑似文件下载过程中会存在网络问题，建议自行下载文件并修改setup.py文件，过程如下:

1. 下载并解压gmsh-4.6.0.tar.gz
   网址：https://pypi.org/project/gmsh/#files

2. 查看自己系统类型
   在python中查看(该代码源自文件中的setup.py的代码)

 import platform
print(platform.system().lower())
print(platform.machine().lower())

3. 手动下载对应的gmsh版本
   网址：https://gmsh.info/bin
   64位windows对应 http://gmsh.info/bin/Windows/gmsh-4.6.0-Windows64-sdk.zip
   文件大小30M
   将下载好的文件放到gmsh-4.6.0.tar.gz 压缩包/gmsh-4.6.0/gmsh-4.6.0/下

4. 修改gmsh-4.6.0.tar.gz 压缩包/gmsh-4.6.0/gmsh-4.6.0/setup.py
   删除第50行的 self._download()
   删除第55行~60行的 _download(self) 函数

5. 在终端切换目录至 X/gmsh-4.6.0/gmsh-4.6.0/ ，输入指令：

python setup.py install

备注：如果出现"the wheel package is not available"的提示，请安装wheel库

pip install wheel

6. 重新开启终端，检查gmsh库是否可用

python
import gmsh

没有异常证明可用

7. 享受gmsh的python API吧，骚年！

直接用python会提示没有gmsh模块
pip install gmsh一直卡在using cached gmsh-4.6.0.tar.gz这一步
不确定gmsh的python接口如何使用

受版权影响，原来习惯用的ICEM不能再使用，对OpenFOAM自带的blockMesh和snappyHexMesh使用水平不够，无法应用于复杂模型的网格划分，所以想请教一下看有没有类似于ICEM、pointwise之类的开源网格软件，谢谢！

请问，simpleFoam中的SIMPLE和SIMPLEC算法有区分吗？
如果有，在什么文件下的什么指令进行区分？
另外，SIMPLER植入了simpleFoam吗？
谢谢！

@东岳 谢谢东岳老师！

@cccrrryyy 哦哦，还有这种操作，谢谢您啦

fvSolution下有一段：

cache
{
    grad(U);
}

不知道什么意思

turbulenceProperties中object为RASProperties或turbulenceProperties有什么区别？

@Samuel-Tu 谢谢您的回答。
plot in line我在使用的时候发现有两个需求无法满足：一是似乎只能一条线；另一个是只能显示而不能输出数据，
我的目标是对流场中的几条线上的物理量进行对比。
最后采用了监控数据的方式，streamline导出的数据不知道为什么有些诡异--（一组出现并行核数次重复数据；三组线没感觉出来，最后的那一大堆不知道什么东东）。probes虽然在输入时需要输很多点，但至少最后能得到想要的数据。

在流场中创建一条线，并把这条线上的物理量（比如速度分量Uy）输出（文本形式）。
非规则形网格这样输出的话，如何给点呢？