有谁对OpenFOAM+GPU感兴趣么？

[尚善若水]

@李东岳 怪不得呢，不过 exaFoam 好像有基于基金会版本的 GPU 加速，可惜没开源，我找下，记得这个文章我下了

[李东岳]

我目前能找到的，基于新版本openfoam的，应该就是挂petsc那个了。之前有一些老的能gpu算的，都是特别老的openfoam版本，2.3, 2.2这种。

如果挂petsc那个GPU速度比较快的话，可以把ESI那个petsc+gpu移植玩一下，比如让基金会版本的也可以gpu计算，目前我看这种方法最合适。之前rapidcfd他们做的openfoam23x的gpu计算，计算过程中全部通过gpu实现了，代码改动不小，移植乏力

目前也不知道速度怎么样。也得测试一下。

[尚善若水]

• （1）exasim-project 的 GPU 好像是做了一个接口平台 OGL，OpenFOAM Ginkgo Layer，使用Ginkgo的GPU功能；

• （2）米兰理工一个组的工作，好像参与了exaFoam项目，但是代码还没开源
  https://doi.org/10.1016/j.compfluid.2023.106041
  https://onlinelibrary.wiley.com/doi/pdf/10.1002/fld.5297
  https://www.mdpi.com/2226-4310/10/9/792

• （3）Venkat Raman 教授组里也有in-house的GPU
  https://doi.org/10.1016/j.compfluid.2023.105972

• （4）deepflame好像也在做GPU

[Amadeus]

@尚善若水 牛逼了这都能找到....我关注exaFoam一段时间了，还挺期待OpenFOAM能用GPU的

[尚善若水]

@Amadeus 一直在关注这方面，有文章就收集起来了。但是短期内OpenFOAM特别是基金会官方支持GPU不太现实

[Amadeus]

@尚善若水 在 有谁对OpenFOAM+GPU感兴趣么？ 中说：

deepflame好像也在做GPU

24年6月北京OpenFOAM Workshop的时候就听陈炽教授介绍了deepflame，确实是在加入深度学习和GPU

确实org这边用上GPU可能很难了，毕竟ESI那边更有钱?

[尚善若水]

@Amadeus 我觉得主要还是李老师说的 ESI 那面更开放吧，Jasak 就有参与 exaFoam，和 ESI 联系也紧密。基金会感觉主要靠资助，对贡献的代码质量要求也很高，每次更新变化也非常非常大（但是我还是更喜欢用基金会版本😂）

[李东岳]

https://github.com/hpsim/OGL

我看这个是已经出来了？

[thegame]

@李东岳

@李东岳 在 有谁对OpenFOAM+GPU感兴趣么？ 中说：

看错了 @thegame 我以为3060的显卡，原来他们用的A10

这个图里面16c，是16核还是16个CPU？是不是16核？如果16核的话，1个CPU跟4c+4g差不多的速度。莫非是16个CPU，一共16乘以24个核？

后面的4c+4g，是4核+4个GPU？还是4个CPU4个GPU

给我整懵了

我理解的c是CPU。4c 4g也是4个CPU+4个GPU。
这个方法也有一些问题：
这个方法矩阵组装那些是没有GPU加速的，而且涉及到OpenFOAM的LDU矩阵转换为CUDA的CSR会有额外的计算开销。所以甚至用这个方法还会拖慢速度，见下图。我理解这种方法可能只对DNS级别，矩阵很大的计算加速效果较好。对于物理过程复杂的多相流反应流那些加速效果，效果不好。这个方法最大的优势就是不用改代码，只需要外挂库就好。
这套我是在笔记本WSL都可以编，就是按照那个网页的教程。

另：最近狠狠打工，我也想试试李老师的大4090但是真的没时间（我之前试着在国内某超算上弄过，编译通过了，但是调多个GPU计算报错了，所以估摸得专门用一段时间来折腾）

[cccrrryyy]

我可太感兴趣了，有空试试。

[李东岳]

@thegame 多谢反馈。是，我看他们讨论全部GPU加速，以及这种GPU求解矩阵的区别。如果用GPU求解矩阵会导致额外开销的话，那看起来还是全部GPU加速比较好。我看rapidCFD他们应该是全部GPU加速的。就是太老了。还是cuda8的。据我所知fluent那面应该也是全部搞到GPU去了。

@cccrrryyy 大佬感兴趣试试么，试试petsc那种的，我这面有个2000万网格的，我想看看实际效果。

[cccrrryyy]

@李东岳 兴趣大大的，但最近可能要再过个一两周，等开年过了时间就比较空了

[李东岳]

@cccrrryyy 那过段时间联系 老铁

[李东岳]

另外一个就是把rapidcfd重写，目前这个适用于of23，cuda8，但是因为我不懂cuda，所以这个活我也干不了。