dpm 模型 求解若干问题 求助

lscookie

最近一直在研究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）和速度矢量图有问题

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

[ 本帖最后由 lscookie 于 2011-4-26 21:42 编辑 ]