@李东岳 收到 多谢李老师

@AppleKiller 我知道一个开放数据集，但是里边单气泡这部分没开放，我现在密码过期了看不到是u否开放了，可以问问他们http://141.44.132.124/spp1740/index.html要是开放了也@我一下

各位大佬好，我最近在尝试模拟球体在井筒中随流体流动并封堵射孔点，模拟的球直径为15mm，井筒直径为110mm，孔眼采用螺旋布孔，将六个射孔点设置为出口，每个孔眼的直径为12mm，使用清水，流速设置为4m/s。模拟采用全解析，网格划分大小在3mm左右，在twoSpheresGlowinskiMPI案例基础上进行修改。现在发现球体堵住出口后，球体速度显示为0，被堵出口处流速降低，但仍有流速，想请教一下各位大佬有可能是哪方面的原因呢。以下是我目前的模拟结果
75f0fe5d-d5ea-499c-aeff-49b56b7540fe-image.png

算出来的结果激波和边界层都很稳定，激波不会来回移动的很明显。这个是代码naca0012_ma-0.7_re-3_alpha-5_sa.zip

@wangfei9088 多谢指教，已解决

@lihaoyang 我得虚拟机不需要安装这些。里面都有了啊。

看起来好像没成功。不知道你运行了几次。

@李东岳 这样做确实是可行的，我也尝试下修改源代码试试，把扩散项加进去

@李东岳
谢谢李老师，我用以下脚本把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 大概率是网格问题，你就用最常用的数值格式（比如参考商业软件的默认设置），先把网格优化下。

谢谢老师，太牛逼了。

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