Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • 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

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 用pimpleFOAM做圆柱的涡激振动遇到收敛性问题

用pimpleFOAM做圆柱的涡激振动遇到收敛性问题

已定时 已固定 已锁定 已移动 OpenFOAM
12 帖子 3 发布者 9.2k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 赵 离线
    赵 离线
    赵海盛
    写于 最后由 李东岳 编辑
    #1

    求助啊,大兄弟们!!!错误信息如下:

    PIMPLE: Iteration 1
    Restraint verticalSpring:  attachmentPt - anchor (0 1.14264e-05 0) spring length 1.14264e-05 force (-0 -4.72392e-05 -0)
    6-DoF rigid body motion
        Centre of rotation: (0 1.14788e-05 0.084)
        Centre of mass: (0 1.14788e-05 0.084)
        Orientation: (1 0 0 0 1 0 0 0 1)
        Linear velocity: (0 2.61184e-05 0)
        Angular velocity: (0 0 0)
    GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.0197574, No Iterations 40
    GAMG:  Solving for pcorr, Initial residual = 0.27233, Final residual = 0.0156253, No Iterations 4
    time step continuity errors : sum local = 3.08875e-14, global = 7.0696e-15, cumulative = 9.56635e-09
    smoothSolver:  Solving for Ux, Initial residual = 3.30404e-06, Final residual = 3.30404e-06, No Iterations 0
    smoothSolver:  Solving for Uy, Initial residual = 0.000208158, Final residual = 1.22015e-07, No Iterations 1
    smoothSolver:  Solving for Uz, Initial residual = 2.49119e-05, Final residual = 1.94723e-08, No Iterations 1
    GAMG:  Solving for p, Initial residual = 0.212445, Final residual = 0.00136547, No Iterations 4
    GAMG:  Solving for p, Initial residual = 0.0257973, Final residual = 9.60752e-07, No Iterations 45
    time step continuity errors : sum local = 1.04039e-13, global = -1.47872e-14, cumulative = 9.56634e-09
    PIMPLE: Iteration 2
    smoothSolver:  Solving for Ux, Initial residual = 1.52276e-06, Final residual = 1.52276e-06, No Iterations 0
    smoothSolver:  Solving for Uy, Initial residual = 9.39697e-05, Final residual = 1.19676e-07, No Iterations 1
    smoothSolver:  Solving for Uz, Initial residual = 1.10536e-05, Final residual = 2.00238e-08, No Iterations 1
    GAMG:  Solving for p, Initial residual = 0.25598, Final residual = 0.00164346, No Iterations 4
    [0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    [0] #1  Foam::sigSegv::sigHandler(int) at ??:?
    [0] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
    [0] #3  Foam::lduMatrix::Amul(Foam::Field<double>&, Foam::tmp<Foam::Field<double> > const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char) const at ??:?
    [0] #4  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
    [0] #5  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
    [0] #6  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
    [0] #7  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
    [0] #8  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
    [0] #9  Foam::fvMatrix<double>::solve() in "/home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
    [0] #10  ? in "/home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
    [0] #11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    [0] #12  ? in "/home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
    [111:1838841] *** Process received signal ***
    [111:1838841] Signal: Segmentation fault (11)
    [111:1838841] Signal code:  (-6)
    [111:1838841] Failing at address: 0x3e8001c0ef9
    [111:1838841] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f5a7570c210]
    [111:1838841] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f5a7570c18b]
    [111:1838841] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f5a7570c210]
    [111:1838841] [ 3] /home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam9lduMatrix4AmulERNS_5FieldIdEERKNS_3tmpIS2_EERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEh+0x1dc)[0x7f5a76140c3c]
    [111:1838841] [ 4] /home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0x62)[0x7f5a76166dc2]
    [111:1838841] [ 5] /home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x885)[0x7f5a7616ae35]
    [111:1838841] [ 6] /home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x524)[0x7f5a7616d154]
    [111:1838841] [ 7] /home/abc/OpenFOAM/OpenFOAM-8/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x18b)[0x7f5a783322ab]
    [111:1838841] [ 8] pimpleFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x1e8)[0x55f9e9e962e8]
    [111:1838841] [ 9] pimpleFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x120)[0x55f9e9e96590]
    [111:1838841] [10] pimpleFoam(+0x2e50b)[0x55f9e9e4450b]
    [111:1838841] [11] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f5a756ed0b3]
    [111:1838841] [12] pimpleFoam(+0x313fe)[0x55f9e9e473fe]
    [111:1838841] *** End of error message ***
    --------------------------------------------------------------------------
    Primary job  terminated normally, but 1 process returned
    a non-zero exit code. Per user-direction, the job has been aborted.
    1 条回复 最后回复
  • 同学博同 离线
    同学博同 离线
    同学博
    写于 最后由 编辑
    #2

    这个看不出什么具体问题吧,应该就是发散了,可以检查一下设置,看看边界条件:146: 或者加个松弛试试?

    赵 1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵海盛
    在 中回复了 同学博 最后由 编辑
    #3

    @同学博 已经加过松弛了,速度0.7,压力0.5,还是算不过去

    1 条回复 最后回复
  • F 离线
    F 离线
    Fu
    写于 最后由 编辑
    #4

    二维还是三维呢。。边界条件问题吧一般是

    赵 1 条回复 最后回复
  • 同学博同 离线
    同学博同 离线
    同学博
    写于 最后由 编辑
    #5

    网格量多少,不大的话可以分享一下case,我有时间的时候可以帮忙简单看一看

    赵 2 条回复 最后回复
  • 赵 离线
    赵 离线
    赵海盛
    在 中回复了 同学博 最后由 编辑
    #6

    @同学博 :140: 非常感谢!cylinderWithBase1.zip

    1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵海盛
    在 中回复了 Fu 最后由 编辑
    #7

    @Fu 边界条件挺简单的,应该没啥问题吧。
    U_BC.png
    p_BC.png
    Displacement_BC.png

    1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵海盛
    在 中回复了 同学博 最后由 编辑
    #8

    @同学博 网格有60万,还是挺大的。方便时兄弟帮我看一下我这些文件设置有没有问题吧。计算的问题是“一个带底座的短圆柱在恒定来流下的涡激振动”。

    同学博同 F 2 条回复 最后回复
  • 同学博同 离线
    同学博同 离线
    同学博
    在 中回复了 赵海盛 最后由 编辑
    #9

    @赵海盛 我还没下载 我看你帖子里发的图片你的pointDisplacement设置的有问题吧,直接设置成calculated能振起来吗,我做过简谐振动的,需要在pointDisplacement把边界条件的振动模式设置好。你可以参考一下这篇文章“A tutorial of the sixDofRigidBodyMotion library with multiple bodies”。cy等个后续:140:

    赵 1 条回复 最后回复
  • F 离线
    F 离线
    Fu
    在 中回复了 赵海盛 最后由 编辑
    #10

    @赵海盛 看起来没有问题的。。如果是自激振荡就是caculated,如果是强迫就要oscillatingdisplacement了。。。看边界条件是一个三维的圆形流场吗?。。

    赵 1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵海盛
    在 中回复了 同学博 最后由 编辑
    #11

    @同学博 谢谢兄弟哈,pointDisplacement设置应该没啥问题。现在找到原因了,我其他文件中参数值改了,forceCoeffs中参数没改,晕了。都改完后已经能正常算了。但forceCoeffs不应该只是输出力参数了?

    1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵海盛
    在 中回复了 Fu 最后由 编辑
    #12

    @Fu 是的,大兄弟。我其他文件中改参数了,但forceCoeffs中忘改了,都改后已经可以算了。谢谢哈!

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]