24小时热门版块排行榜    

查看: 1398  |  回复: 1

lscookie

铜虫 (小有名气)

[求助] dpm 模型 求解若干问题 求助

最近一直在研究dpm模型 现在对问题设置完后进行了初步求解,遇到一些问题,希望同在研究dpm的同志们能给些建议, 高手们能指点一下, 非常感谢
下图是一个电解槽阳极及其计算域 。在阳极底部面上定义了file型颗粒,共400个粒子初始速度小于0.004m/s,方向大小随机分布,位置在面内随机分布,直径小于1cm随机排布,温度400k;问题是在静止的连续相电解质流体中有离散相气泡喷出, 研究离散相对连续相流场的影响,以及对比不同模型下离散相逃逸的难易程度。
网格划分后 最小网格尺寸10mm,阳极底部距逃逸面(阳极顶部粉色区域)400mm,考虑重力场,粒子喷射时间足够长end time =40s,初始化(如下图)后,颗粒受力物理模型中仅考虑了Two-Way Turbulence Coupling ,模型中离阳极较近的面及其右侧面是对称边界条件,另外两个侧壁为回弹的wall,reflect条件,顶面是escape边界,底部也是reflect, 采用了两种求解方式:(稳态追踪方式下发现max.num设大了也会出现incomplete或计算很缓慢)
A 非稳态追踪方式的稳态求解:(考虑到瞬态求解时时间步长过小,到粒子escape所用时间步理论上太大),每十步连续相计算后及算离散相轨道,max.num =500000,Step Length Factor   =4,为尽快查看ecsape情况,Particle Time Step Size   =0.2,Number of Time Steps   =1即每十步后喷射一次injection,迭代1500步后,发现以下问题:
1.连续相残差很难收敛,在1e-2到1e-3间震荡,调整松弛因子后仍有这种现象,网格全是结构化网格,检查质量时都比较好;
2.湍流云图显示结果和文献差异较大,速度矢量图显示流场速度最高达3m/s ,文献中最大在0.4m/s附近,不知哪些设置上出了问题



B非稳态追踪方式的瞬态求解:(考虑流动非定常很明显的话稳态求解不太合理)
连续相时间步长设置为0.002(试过0.0005和0.001每时间步内收敛,为了减小计算时间设0.002),每步最大迭代30步,共计算了2000步,发现以下问题:
1.前2000步基本都满足每步收敛,后面开始有出现不收敛情况(主要还是连续相上);
2.前1000步手动停止查看结果同A中情况一样,部分区域湍流云图(Turbulent Kinetic Energy)和速度矢量图有问题

请大家帮忙给些修改意见或指点一下 不胜感激








回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lscookie

铜虫 (小有名气)

哪位 做过dpm的高手帮忙给指点下吧
2楼2011-04-28 15:59:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lscookie 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于代码变化问题,想知道的进来 +14 且听虎啸 2026-08-07 20/1000 2026-08-09 18:03 by 天神眷顾
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +7 Lanmanbaby 2026-08-09 11/550 2026-08-09 17:48 by zhanghaozhu
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见