@李东岳
谢谢李老师，我用以下脚本把inlet改正了过来，目前能够重组网格了。

#!/bin/sh timestep=20 for processor in processor*; do sed -i 's/codedFixedValue/fixedValue/' $processor/$timestep/turbulenceProperties:R sed -i 's/name.*;/ /' $processor/$timestep/turbulenceProperties:R done

之后也能够在paraview中检查。

@李东岳 谢谢李老师，我试试

@李东岳 在 请教各位如何用tecplot360(EX2022R1)处理OpenFOAMv2312中颗粒分布 中说：

还真是。我刚才看了一下。应该是我的虚拟机上面没有装tecplot，所以没找到tec的头文件。

dyfluid@dyfluid-virtual-machine:~/OpenFOAM/OpenFOAM-8/applications/utilities/postProcessing/dataConversion/foamToTecplot360$ wmake wmakeLnIncludeAll: running wmakeLnInclude on dependent libraries: wmakeLnInclude error: base directory /home/dyfluid/OpenFOAM/ThirdParty-8/tecio/tecsrc/ does not exist Making dependency list for source file foamToTecplot360.C could not open file TECIO.h for source file foamToTecplot360.C due to No such file or directory could not open file MASTER.h for source file foamToTecplot360.C due to No such file or directory could not open file GLOBAL.h for source file foamToTecplot360.C due to No such file or directory Making dependency list for source file vtkMesh.C Making dependency list for source file tecplotWriter.C could not open file TECIO.h for source file tecplotWriter.C due to No such file or directory could not open file MASTER.h for source file tecplotWriter.C due to No such file or directory could not open file GLOBAL.h for source file tecplotWriter.C due to No such file or directory g++ -std=c++11 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -O3 -DNoRepository -ftemplate-depth-100 -I/home/dyfluid/OpenFOAM/ThirdParty-8/tecio/tecsrc/lnInclude -I/home/dyfluid/OpenFOAM/OpenFOAM-8/src/lagrangian/basic/lnInclude -I/home/dyfluid/OpenFOAM/OpenFOAM-8/src/finiteVolume/lnInclude -I/home/dyfluid/OpenFOAM/OpenFOAM-8/src/dynamicMesh/lnInclude -I/home/dyfluid/OpenFOAM/OpenFOAM-8/src/meshTools/lnInclude -IlnInclude -I. -I/home/dyfluid/OpenFOAM/OpenFOAM-8/src/OpenFOAM/lnInclude -I/home/dyfluid/OpenFOAM/OpenFOAM-8/src/OSspecific/POSIX/lnInclude -fPIC -c tecplotWriter.C -o /home/dyfluid/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/applications/utilities/postProcessing/dataConversion/foamToTecplot360/tecplotWriter.o In file included from tecplotWriter.C:26: tecplotWriter.H:39:10: fatal error: TECIO.h: 没有那个文件或目录 39 | #include "TECIO.h" | ^~~~~~~~~ compilation terminated. make: *** [/home/dyfluid/OpenFOAM/OpenFOAM-8/wmake/rules/General/transform:26：/home/dyfluid/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/applications/utilities/postProcessing/dataConversion/foamToTecplot360/tecplotWriter.o] 错误 1 dyfluid@dyfluid-virtual-machine:~/OpenFOAM/OpenFOAM-8/applications/utilities/posttProcessing/tProcessing/dataCotProttProcessing/dataConversion/foamToTecplot360$

应该确实是，我看了一下时间，我openfoam-8装的时候电脑上已经有装tecplot2023了。但openfoam-7装的时候没有，也许就是编译记录的Log文件里就没有foamToTec的原因。

@dyb 我的理解哈，这个视频我看过的。这里其实就是常规的合成湍流入口，只不过作者给它套了个风资源里面常见的名词（阵风），很有可能它根本就不是我们经常理解的阵风的概念。

并没有把阵风加到湍流里面这种说法，阵风本来就是风流大幅变化的一种叫法而已。既然做LES，时间尺度肯定已经远小于阵风所在的秒级尺度了。

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

torch.set_flush_denormal(True)

https://discuss.pytorch.org/t/training-time-gets-slower-and-slower-on-cpu/145483/4

我在计算固体冲压发动机内两相流颗粒二次燃烧时，C颗粒燃烧采用openfoam自带的模型，B颗粒的燃烧自己编写的表面燃烧模型，但计算时总是在空气和颗粒接触不久之后发散。MaxCo为1时，自适应时间步长大概在2e-6s左右，多次尝试后存在以下几种情况：
（1）两种颗粒的Cloud都启用表面燃烧时，不管如何调整species的Cp等物性参数，均会在0.0068s时出现Energy -> temperature conversion failed to converge，输出的报错信息显示，温度在1000±1K之间来回振荡；
（2）两种颗粒的Cloud都启用表面燃烧，但B颗粒所有反应的速率都设为0时，计算在0.0116s发散，从报错信息来看，应该是ReactingMultiphaseParcel中的calc函数调用后，MomentumParcel的move()函数计算出错；
（3）C颗粒启用表面燃烧模型，B颗粒的表面燃烧模型选择none，此时计算不会发散。

对于情况（1），在将thermo.C文件中的tol_从1e-4改为2e-4后能继续往后计算，到0.0073s时还是发散了。但情况（2）反应速率均为0的情况下，在我理解中应该跟情况（3）是等价的，不知道为什么一个报错，一个却能顺利计算。

@李东岳 明白了，谢谢老师

@李东岳
谢谢老师，一开始我认为是入口的速度梯度计算的不准确，所以用差分方法直接计算。后面对照下看来是接近的。
我认为可能是不可压缩的算法里不考虑严格的能量守恒，所以fb4e0ad7-cdcb-4682-ab71-5ada18146f9d-image.png 直接计算并不为0.
但是上面公式推导中使用高斯方法转换体面积分的时候默认了不可压缩时这项为零。应该写出完整的形式计算0ceb3ad3-73d3-4c9e-81b8-f3f473dfcbbb-image.png

Total viscous power in the domain (volume integral): -1.50451907 Total integrated tau : gradU: 1.48509576 Total integrated viscous power (surface): 0.00384099619

这样计算后数值是接近的，等稳定后误差应该还会减小。

@李东岳 testSimplefoam.zip 谢谢老师，我这边上传了，分别是internalField为边界条件和为0的两个0文件，然后我用的是singleGraph3监测x=2m位置处的速度和湍动能信息。网格数量不大。

@xuebao1989 大概率是网格问题，你就用最常用的数值格式（比如参考商业软件的默认设置），先把网格优化下。

谢谢老师，太牛逼了。

挺好，这玩意之前我也遇到，一直没法解决

收藏一下，感觉有用

@李东岳 好的李老师，感谢感谢

@五好青年 老哥，我也会有这个问题，单核正常，并行发散，这是不是与分配核数有关呀？

@李东岳 是可以不选择的，但怕设置错了，这玩意用的人确实不多啊，老师还有啥推荐的网格软件吗

VIV.zip 这是我配置的openFoam动网格算VIV的算例，从目前的结果来看，程序能正常运行，但程序运行结束以后，圆柱不动，想请各位大佬帮我看看，哪儿出问题了，谢谢大家。

@疏影横斜水清浅 重力以及体积力

@sgggggg 硕士论文我也想拜读一下，可以加个联系方式吗2132306535