24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1297  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 找博士生导师 +3 小代想上岸 2026-05-15 4/200 2026-05-17 09:20 by YuY66
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 7hlccowb3h 2026-05-15 4/200 2026-05-17 07:46 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 8/400 2026-05-17 07:26 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:11 by 11n4dfd8yn
[基金申请] 精华III评审感受-评审感受-评审感受 +16 ferrarichen 2026-05-11 20/1000 2026-05-17 01:10 by 南开小綦
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:55 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 x0mp7owy2b 2026-05-15 5/250 2026-05-17 00:29 by xiangfeng
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 4/200 2026-05-17 00:15 by ue3ir18jc3
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 00:10 by ue3ir18jc3
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +9 一切都是空工 2026-05-12 20/1000 2026-05-16 19:52 by Equinoxhua
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
请填处理意见