CFD中文网

    CFD中文网
    • 登录
    • 搜索
    • 最新

    OpenFoam中的溃坝算例加入自适应网格计算一定时间步后出现的问题

    OpenFOAM
    3
    6
    1454
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 云
      云翳的黄昏 最后由 编辑

      请问各位大佬有没有曾经遇到过算例计算过程中提示“浮点数例外 (核心已转储)”
      ，谢谢帮忙解答！
      我最近在学习interFOAM求解器相关知识，目前向使用interFOAM求解器加上自适应网格计算两相流问题，在利用OpenFoam自带算例修改无报错后，开始利用interDyMFoam进行求解。但何运算基本在0.3s附近出现下图中的错误中断（of2.3.1和of7中都出现了该报错）
      ![0_1596110438484_2020-07-30 19-57-10屏幕截图.png](正在上传 100%)
      ![0_1596110989816_2020-07-30 20-06-37屏幕截图.png](正在上传 100%)
      而且还有一个问题，但何运行计算到0.3s附近出现该错误，双核并行0.1s，四个核时亦不就中断了。无自适应时算例能正常计算。
      麻烦各位大哥有知道y原因的帮忙回答一下，谢谢!

      云 1 条回复 最后回复 回复 引用
      • 云
        云翳的黄昏 @云翳的黄昏 最后由 李东岳 编辑

        @云翳的黄昏

        OpenFoam-2.3.1
Interface Courant Number mean: 0.0038340104 max: 0.84420783
Courant Number mean: 0.037028394 max: 0.99652771
deltaT = 5.2872478e-05
Time = 0.31246054

Selected 124 cells for refinement out of 273070.
Refined from 273070 to 273938 cells.
Selected 83 split points out of a possible 32348.
Unrefined from 273938 to 273357 cells.
Execution time for mesh.update() = 3.39 s
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.8048258e-06, No Iterations 389
time step continuity errors : sum local = 5.4228827e-10, global = -2.51365e-13, cumulative = 0.00017816493
smoothSolver:  Solving for alpha.water, Initial residual = 0.00055204007, Final residual = 6.4208784e-09, No Iterations 2
Phase-1 volume fraction = 0.13010958  Min(alpha1) = -9.3530875e-08  Max(alpha1) = 1.0000001
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.13010958  Min(alpha1) = -1.3561514e-05  Max(alpha1) = 1.0000079
GAMG:  Solving for p_rgh, Initial residual = 0.033205248, Final residual = 0.00014731475, No Iterations 2
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3   at tensorField.C:?
#4   in "/home/vaseatar/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/bin/interDyMFoam"
#5   in "/home/vaseatar/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/bin/interDyMFoam"
#6   in "/home/vaseatar/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/bin/interDyMFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8   in "/home/vaseatar/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/bin/interDyMFoam"

        上面图片好像没有显示，抱歉！这里添加一下报错误的文本

        OpenFoam-7
Courant Number mean: 0.0742364 max: 0.985473
Interface Courant Number mean: 0.0103739 max: 0.761644
deltaT = 0.000210711
Time = 0.331036

Selected 20 cells for refinement out of 36008.
Refined from 36008 to 36148 cells.
Selected 7 split points out of a possible 4206.
Unrefined from 36148 to 36099 cells.
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.59577e-06, No Iterations 208
time step continuity errors : sum local = 1.48015e-09, global = 1.6711e-12, cumulative = 0.000723778
smoothSolver:  Solving for alpha.water, Initial residual = 0.00142426, Final residual = 1.32652e-09, No Iterations 3
Phase-1 volume fraction = 0.130119  Min(alpha.water) = -4.43817e-11  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.130119  Min(alpha.water) = -3.77094e-05  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0515556, Final residual = 0.0021802, No Iterations 2
time step continuity errors : sum local = 0.000166589, global = -9.96186e-07, cumulative = 0.000722782
DICPCG:  Solving for p_rgh, Initial residual = 0.00314206, Final residual = 0.000140061, No Iterations 10
time step continuity errors : sum local = 1.0632e-05, global = -3.56508e-07, cumulative = 0.000722425
DICPCG:  Solving for p_rgh, Initial residual = 0.000874623, Final residual = 9.92427e-08, No Iterations 170
time step continuity errors : sum local = 7.53188e-09, global = -2.18209e-10, cumulative = 0.000722425
ExecutionTime = 431.97 s  ClockTime = 432 s

Courant Number mean: 0.0750393 max: 0.990053
Interface Courant Number mean: 0.0104568 max: 0.765205
deltaT = 0.000210711
Time = 0.331247

Selected 20 cells for refinement out of 36099.
Refined from 36099 to 36239 cells.
Selected 12 split points out of a possible 4220.
Unrefined from 36239 to 36155 cells.
DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 8.95937e-06, No Iterations 199
time step continuity errors : sum local = 1.55988e-09, global = 1.04349e-11, cumulative = 0.000722425
smoothSolver:  Solving for alpha.water, Initial residual = 0.00142307, Final residual = 1.50653e-09, No Iterations 3
Phase-1 volume fraction = 0.130116  Min(alpha.water) = -4.43604e-11  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.130116  Min(alpha.water) = -4.22756e-05  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 0.0570999, Final residual = 0.00183939, No Iterations 3
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  ? at tensorField.C:?
#4  ? in "/home/vaseatar/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/interFoam"
#5  ? in "/home/vaseatar/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/interFoam"
#6  ? in "/home/vaseatar/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/interFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  ? in "/home/vaseatar/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/interFoam"
        1 条回复 最后回复 回复 引用
        • 云
          云翳的黄昏 最后由 编辑

          这个应该是自适应在二维网格下计算过程中出现的问题，后来尝试把二维算例改成三维算例。就不会出现计算过程中断的问题了。因为我目前还在学习，具体原因尚不清楚。有对这一部分深入研究的可以补充，我也学习一下，谢谢！也帮助我们这些初级学习者快速入门，尽快排除bug

          1 条回复 最后回复 回复 引用
          • 史
            史浩 讲师 最后由 编辑

            这个主要问题是OpenFOAM自带的AMR是基于三维网格的，二维网格就会报错。关于2D自适应网格的资料看这里

            让我们随波逐流

            云 L 2 条回复 最后回复 回复 引用
            • 云
              云翳的黄昏 @史浩 最后由 编辑

              @史浩 恩恩，谢谢大哥！我好好学习一下，有问题再请教！

              1 条回复 最后回复 回复 引用
              • L
                Lisiye @史浩 最后由 编辑

                @史浩 您好，请问对于自适应网格的设置讲解有什么资料可以帮助我学习参数设置的吗

                1 条回复 最后回复 回复 引用
                • First post
                  Last post