Skip to content
  • 最新
  • Categories
  • 东岳流体
  • 随机看[请狂点我]
Skins
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • Default (No Skin)
  • No Skin
Collapse
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. foam-extend的fsiFoam并行效率问题

foam-extend的fsiFoam并行效率问题

Scheduled Pinned Locked Moved OpenFOAM
7 Posts 4 Posters 6.1k Views
  • Oldest to Newest
  • Newest to Oldest
  • Most Votes
Reply
  • Reply as topic
Log in to reply
This topic has been deleted. Only users with topic management privileges can see it.
  • Number44N Offline
    Number44N Offline
    Number44
    wrote on last edited by
    #1

    通过将fsiFoam算例beamInCrossFlow网格等比放大至800万流体网格+800万固体网格,其他参数不变,foam4.0版本,我测试了一下fsiFoam的并行效率,结果发现fsiFoam大规模并行效率非常糟糕,规整的六面体网格,200核以上直接跑不动,峰值在50核左右,和串行比加速比不超过25。
    仔细地对fsiFoam的代码计时,发现求解过程正常,根源出现在固体求解器的一行最小二乘法的cell to point interpolation上,源代码如下:

    头文件声明:

    //- Cell-to-point interpolation
    leastSquaresVolPointInterpolation volToPoint_;
    

    .C文件:

    volToPoint_.interpolate(D_, pointD_);
    

    在这个插值函数里面,用的大量的scheduled形式的MPI通信,MPI_Probe函数占用了大量的时间。
    底层的东西我改不了,就不求甚解地直接把这个插值类从最小二乘法的leastSquaresVolPointInterpolation 改成了反距离加权法的volPointInterpolation,编译正常,报错的地方改好了且不参加计算。
    (foam-extend提供了volPointInterpolation、leastSquaresVolPointInterpolation 和edgeCorrectedVolPointInterpolation,第三个的输入参数不一样就没用)
    但是计算结果出了问题,网格变形和原算例的不一致,原算例的网格变形类似于倾斜的平行四边形,改了之后的网格变形像个橡皮泥一样从四边形迎流面中部凹陷,背流面也错位凸起。
    图片传不上去,服务器解析出错,后面再补。

    请问各位大佬,这两个插值类的区别那么大,不能混用么,还是我代码改得有问题?我看foam-extend的icoFsiFoam求解器的算例也是一个类似的二维受限空间平板绕流的流固耦合算例,用的volPointInterpolation就没问题。

    算不准,发个散,报error,没问题!

    Number44N 1 Reply Last reply
  • Number44N Offline
    Number44N Offline
    Number44
    replied to Number44 on last edited by
    #2

    补个图,上面是源代码算的结果,下面是改了volPointInterpolation算的结果。

    compare.jpg

    算不准,发个散,报error,没问题!

    1 Reply Last reply
  • S Offline
    S Offline
    shiyu
    wrote on last edited by
    #3

    Hi, Number44
    我遇到和您类似的问题,不过我的case简单不需要并行,但对时间cost要求很严格,我发现我改动的FSI很花费时间,可能与你说的插值函数有关,不知你有改善效率问题吗?
    其次,想问一下,除了foam-extend里面有FSI tool,基金会的openfoam有类似的吗?比较偏向于使用基金会版本,谢谢

    Number44N 1 Reply Last reply
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    wrote on last edited by
    #4

    @Number44 在 foam-extend的fsiFoam并行效率问题 中说:

    过将fsiFoam算例beamInCrossFlow网格等比放大至800万流体网格+800万固体网格,其他参数不变,foam4.0版本,我测试了一下fsiFoam的并行效率,结果发现fsiFoam大规模并行效率非常糟糕,规整的六面体网格,200核以上直接跑不动,峰值在50核左右,和串行比加速比不超过25

    这看起来像是个bug,可以跟foam-extend官方提bug

    9月CFD算法编程课: http://dyfluid.com/class.html

    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    Number44N 1 Reply Last reply
  • Number44N Offline
    Number44N Offline
    Number44
    replied to shiyu on last edited by
    #5

    @shiyu 在 foam-extend的fsiFoam并行效率问题 中说:

    Hi, Number44
    我遇到和您类似的问题,不过我的case简单不需要并行,但对时间cost要求很严格,我发现我改动的FSI很花费时间,可能与你说的插值函数有关,不知你有改善效率问题吗?
    其次,想问一下,除了foam-extend里面有FSI tool,基金会的openfoam有类似的吗?比较偏向于使用基金会版本,谢谢

    能力有限,改了好几次代码,结果都算不对。
    我只用过OpenFOAM5-7,里面应该都没有fsi相关内容,新的OpenFOAM8不清楚。
    另外,我不清楚你用fsi干什么,我之前用过基于OpenFOAM6的一个IBM方法的流固耦合,还行,MIT的一个博士搞的,叫sdfibm,能算固体大位移,暂时不能算形变。效率和OpenFOAM6差不多。

    算不准,发个散,报error,没问题!

    1 Reply Last reply
  • Number44N Offline
    Number44N Offline
    Number44
    replied to 李东岳 on last edited by
    #6

    @李东岳 在 foam-extend的fsiFoam并行效率问题 中说:

    @Number44 在 foam-extend的fsiFoam并行效率问题 中说:

    过将fsiFoam算例beamInCrossFlow网格等比放大至800万流体网格+800万固体网格,其他参数不变,foam4.0版本,我测试了一下fsiFoam的并行效率,结果发现fsiFoam大规模并行效率非常糟糕,规整的六面体网格,200核以上直接跑不动,峰值在50核左右,和串行比加速比不超过25

    这看起来像是个bug,可以跟foam-extend官方提bug

    好不容易在sourceforge注册了个帐号,在foam-extend下写了个ticket说了一下这个问题,不知道维护大佬的看不看。

    算不准,发个散,报error,没问题!

    1 Reply Last reply
  • 五好青年五 Offline
    五好青年五 Offline
    五好青年
    wrote on last edited by
    #7

    @number44 在 foam-extend的fsiFoam并行效率问题 中说:

    foam-extend的icoFsiFoam

    您好,请问下,下面两者有什么区别呢?
    foam-extend中自带的icoFsiFoam
    跟
    单独的fsiFoam

    I am a CFD machine with no emotions.

    1 Reply Last reply

  • Login

  • Don't have an account? Register

  • Login or register to search.
  • First post
    Last post
0
  • 最新
  • Categories
  • 东岳流体
  • 随机看[请狂点我]