在两相流代码中添加空化模型的一些问题



  • 原两相流代码是基于有限差分方法,
    速度压力解耦是使用的投影法,
    ad1f6557-6415-4c89-98bf-5a76e78a121f-image.png
    第一步:求中间速度V*2fe2b97a-41e1-4b66-a878-dd55e6de46ee-image.png
    第二步:5be57b52-cb6c-495b-a90c-1df10b105e62-image.png 》》》》24695d15-3822-4515-81ec-d9ee5b7b2c6e-image.png 求压力
    第三步:求解速度
    c034c027-4f48-4d17-8d2f-751eb2a51127-image.png
    VOF 使用的是THINC/WLIC方法,(https://doi.org/10.1016/j.jcp.2007.06.020)通量更新如下
    be9e2423-3067-4f87-9c8c-ef65a73e77c7-image.png
    328abed6-ee01-447f-a346-7091d012e5fa-image.png
    现在想把SchnerrSauer模型加进去,求解的方程是如下四个
    5edb54a5-f49b-4be4-92e9-d4ad9e88e94b-image.png
    我的思路如下:
    首先在某一点,利用时间步n上的值,r,p利用方程(4)计算得到m,把带入方程(1)(2),利于前面提到的投影法的思路计算中间速度,用方程1,3,得到压力泊松方程(方程里面带有m),使用隐式迭代方法计算压力,用压力计算得到速度,然后使用前面提到的THINC 方法更新相分数,这样整体思路有问题吗,
    如果没有,有几个涉及的细节问题,
    1,因为m里面带有体积分数α,压力P,在求解压力泊松方程的时候如果不考虑空化离散的方程如下
    139b1db0-ed6b-495d-9508-ae01a7b0bd3b-image.png ,使用迭代法求解,考虑空化的情况会有m相关的项,m里面涉及的压力p需要如何考虑,就直接用算出来的m放到压力泊松方程里面,然后 迭代求解吗,
    2.同样相分数更新的时候涉及m的项里面也有α,这个如何考虑,也是直接用算出来的m做一个源项计算吗,
    3,在上面提到的THINC 方法中对相分数的跟新使用的是x,y方向交替更新的,加入空化后涉及到源项m,我在交替更新的时候x,y方向,每一个方向考虑0.5*m可以吗,
    问题有点多,如果有对相关问题熟悉或者了解的前辈,希望不吝赐教,非常感谢。
    或者有涉及到求解附带空化模型的N-S方程算法的文章可以分享吗,可能是我找文献方式有问题,没找到关于如何求解附带空化模型的N-S方程算法的文章



  • 看起来很复杂的样子.... :143:


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2