CFD中文网

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

@cfd-china
恩恩，:laughing:

@him 谢谢解答，另外一般商业软件的定常求解器是用的时间推进来求解的吗？如果是时间推进的定常求解器，这种不稳定的振荡可否理解为流场随时间的变化？

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

各位专家好：
在多孔介质模拟的过程中，我发现了一个问题，虽然在发表文章中省略了。但是确实是一个直观的问题，我的模型，RANS方法中，不管采用哪一种，在多孔介质的近壁面处的拟合的速度与实测值总是具备一定的差距。这种差距我分析认为是在较高阻力降的多孔介质环境下，气流在流经多孔介质表面时，一部分流体会因为高阻降而产生流向变化，形成比较复杂的近壁面流场。而在这个近壁面流场中，因为Fluent默认是采用的多孔介质域或者是多孔跳跃面，应该是对这层近壁面没有处理的，依然按照充分发展层处理，造成这层近壁面的流速（甚至是流向）与实测值都有较大区别。
虽然能够较好的拟合实时压降（Endo based equations），但是如果在近壁面流场不符合实际的话，那么很显然对于最后的表面覆尘的预算难以符合实际情况。不知道大家对于这种多孔介质的近壁面流是否有过研究，或者能够有比较好的解决思路？欢迎大家提出宝贵的想法和建议，感谢，再次感谢！
下图为提出新模型对于滤管外壁的模拟和实测结果对比：
0_1544421734755_compare.png

做静电吸附还是范德华吸附？如果是范德华吸附，可以写一个概率模型。

但是RC插值没有离散，因为它是显式计算源项 -- “速度插值+高阶压力导数项”。

torch.set_flush_denormal(True)

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

@wangfei9088 哈哈 多谢支持多谢支持 :chigua:

@bestucan 谢谢您！学习了

网格加密到Y+<1，计算结果变化不大，对流通量用二阶MUSCL, AUSM格式

如果使用standard KE，计算二维流场，在某个高雷诺数下把Cd, Cl, St都算对，那么说明碰对了，意思就是好运。这并不能说明standard KE可以算对漩涡脱离。

0_1543372305704_捕获.JPG

https://www.ansys.com/-/media/ansys/corporate/resourcelibrary/article/aa-v2-i2-coupling-momentum-and-continuity.pdf

刚看了下，Fluent那面也只是速度和压力耦合求解。能量抛在了外面。

但是另一个文档速度压力能量放在了一起 ： http://www.bakker.org/dartmouth06/engs150/05-solv.pdf

速度压力耦合求解在高超那面用的多，如果分析N-S方程的数学特征（双曲抛物），你需要把u, v, w, p当成矢向量，然后求特征值特征向量。接下来求解这个矩阵就是耦合算法。我目前只是了解到思路，没亲自植入过。可

@李东岳 在 请教各位如何用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的原因。

@cccrrryyy 好的，谢谢，我试试

从特征来看，是双曲形方程。是步进形式的问题。

补充一下，不知道这些图中的这些数据点是如何提取的，谢谢大家

@东岳 这个边界条件已经整合到exteranlHeatFluxTemperature里面了，算的也是温度梯度，也就是$q=-\kappa \frac{\partial T}{\partial n}$，这种边界按照FDS技术手册里面应该是DNS的计算方法。
LES计算时给定的方法是$q=h \Delta T$,$h$是对流换热系数。openfoam里面的边界条件要么是固定值要么是梯度，这种差值不知道在of里面能不能实现

@西湖冷月 这就是天秀

@bestucan 确实可以，谢谢您！