搜索
 找回密码
 立即注册
查看: 147|回复: 10

[电感耦合等离子体] 等离子体计算收敛的问题

[复制链接]
  • TA的每日心情
    郁闷
    2018-1-15 18:46
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    2

    主题

    17

    帖子

    73

    积分

    注册会员

    Rank: 2

    积分
    73
    发表于 2017-12-25 23:18:05 | 显示全部楼层 |阅读模式
    用等离子体做氢气的电离仿真。
    设置了几个简单的反应。

    但是求解的时候出现了以下问题:

    1、用PARDISO求解器的时候提示找不到一致的初始值,模型结果不收敛

    2、使用MUMPS的时候出现了不收敛的情况,MUMPS 内存分配因子增加到 1.44。
    MUMPS 内存分配因子增加到 1.728……就这样不断增加。如图;

    3、使用SPOOLES求解器,发现也是不收敛,计算了很久,如图:

    搞不懂为什么会出现这样的情况?
    求教!

    调节初始参数效果不好。我想应该是模型的问题,但是无从下手。

    本帖子中包含更多资源

    您需要 登录 才可以下载或查看,没有帐号?立即注册

    x
  • TA的每日心情
    郁闷
    2018-1-15 18:46
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    2

    主题

    17

    帖子

    73

    积分

    注册会员

    Rank: 2

    积分
    73
     楼主 发表于 2017-12-26 11:56:43 | 显示全部楼层
    The time-dependent solver could not modify the initial conditions given to a DAE system to satisfy the stationary equations at the initial time. Make sure the initial values satisfy the equations and boundary conditions. In many cases, this can be achieved by solving for only the algebraic variables using a stationary solver before starting the time-dependent solver
    这是官方的提示,但是最后一句 In many cases, this can be achieved by solving for only the algebraic variables using a stationary solver before starting the time-dependent solver具体怎么进行该操作,求教
  • TA的每日心情

    昨天 08:24
  • 签到天数: 348 天

    [LV.8]以坛为家I

    106

    主题

    1009

    帖子

    1万

    积分

    版主

    Rank: 7Rank: 7Rank: 7

    积分
    13258

    荣誉会员学者

    QQ
    发表于 2017-12-26 14:11:08 | 显示全部楼层
    就是让你好好调初始值,很多人调初始值太保守,就改小一点或大一点,如果不行 继续改
    流场,温度场,低温等离子体
  • TA的每日心情
    郁闷
    2018-1-15 18:46
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    2

    主题

    17

    帖子

    73

    积分

    注册会员

    Rank: 2

    积分
    73
     楼主 发表于 2017-12-26 20:11:46 | 显示全部楼层
    bfadlzt 发表于 2017-12-26 14:11
    就是让你好好调初始值,很多人调初始值太保守,就改小一点或大一点,如果不行 继续改 ...

    您是对的。我后来调整了,收敛了
  • TA的每日心情
    郁闷
    2018-1-15 18:46
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    2

    主题

    17

    帖子

    73

    积分

    注册会员

    Rank: 2

    积分
    73
     楼主 发表于 2017-12-26 20:14:53 | 显示全部楼层
    调整了初始值后,最后收敛了。
    这里的初始值涉及到电子的密度、以及其他各种粒子的密度等,而不是单纯只调整电子的密度和电子能量。
    然而这样又存在的问题是没有一个固定的调节标准,感觉自己在盲目试探。
    调整初值后的收敛图如下。

    本帖子中包含更多资源

    您需要 登录 才可以下载或查看,没有帐号?立即注册

    x
  • TA的每日心情
    郁闷
    2018-1-15 18:46
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    2

    主题

    17

    帖子

    73

    积分

    注册会员

    Rank: 2

    积分
    73
     楼主 发表于 2017-12-26 20:15:43 | 显示全部楼层
    之前我只调节了电子密度和能量,没有调整其他粒子。
  • TA的每日心情

    昨天 08:24
  • 签到天数: 348 天

    [LV.8]以坛为家I

    106

    主题

    1009

    帖子

    1万

    积分

    版主

    Rank: 7Rank: 7Rank: 7

    积分
    13258

    荣誉会员学者

    QQ
    发表于 2017-12-26 21:47:50 | 显示全部楼层
    初值是为了更好的收敛,调试初值是有规律的,比如电压很高,初值就不能太大,否则很快就发散了。
    流场,温度场,低温等离子体
  • TA的每日心情
    郁闷
    2018-1-15 18:46
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    2

    主题

    17

    帖子

    73

    积分

    注册会员

    Rank: 2

    积分
    73
     楼主 发表于 2018-1-3 13:33:19 | 显示全部楼层
    另一个收敛的问题。
    仍然是找不到一致初始值,最后一步不收敛。
    在求解的过程中调初始值发现仍旧很难收敛。
    最后能收敛的时间步发现电子的温度很高。电子的温度变化特别剧烈。
    那么问题来了,怎么调整初始值或者耦合关系来保证收敛,求教!

    本帖子中包含更多资源

    您需要 登录 才可以下载或查看,没有帐号?立即注册

    x
  • TA的每日心情
    开心
    昨天 13:02
  • 签到天数: 64 天

    [LV.6]常住居民II

    6

    主题

    87

    帖子

    373

    积分

    中级会员

    Rank: 3Rank: 3

    积分
    373
    发表于 2018-1-3 20:02:28 | 显示全部楼层
    你好,请问你说的调节初值时的其他粒子是参与反应的物质的原子之类的吗?比如氩气的Ar、Ars之类的?
  • TA的每日心情
    郁闷
    2018-1-15 18:46
  • 签到天数: 7 天

    [LV.3]偶尔看看II

    2

    主题

    17

    帖子

    73

    积分

    注册会员

    Rank: 2

    积分
    73
     楼主 发表于 2018-1-15 18:49:56 | 显示全部楼层
    我名字 发表于 2018-1-3 20:02
    你好,请问你说的调节初值时的其他粒子是参与反应的物质的原子之类的吗?比如氩气的Ar、Ars之类的? ...

    有电子、离子以及中性粒子,都是可以调初始值的

    发表回复

    您需要登录后才可以回帖 登录 | 立即注册

    本版积分规则

    10