24小时热门版块排行榜

返回列表

lscookie

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 2.1
红花: 3
帖子: 52
在线: 21.4小时
虫号: 1273505
注册: 2011-04-22
性别: GG
专业: 固体力学

[求助] 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）和速度矢量图有问题

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

回复此楼