我在计算的课题是欠膨胀气液两相流内的激波,想看看滑移速度的影响,有没有可能计算的场景的瞬态特征导致了计算的不稳定啊~
Lubing
帖子
-
drift-flux model中相对速度计算问题 -
drift-flux model中相对速度计算问题这个模型是不是存在一个问题,就是滑移速度在计算时除了和加速度相关,与雷诺数相关,似乎与别的参数没有什么关系。这就会导致随着加速度的增加,滑移速度急剧增加,这里是否还缺少能量上的约束以避免出现不合理的结果呢
-
drift-flux model中相对速度计算问题把源代码放上来,我用的是v2312版本的OpenFOAM,tutorials里面有参考算例
compressibleDriftFoam.zip -
drift-flux model中相对速度计算问题在将这个方程植入到OpenFOAM中之后,出现了明显的相对速度的震荡,惯性加速度最大能到10^5m/s2。植入时直接使用- fvc::ddt(U) - (U & fvc::grad(U))进行计算,使用的时1D的气液激波管问题做的验证,此时无滑移时的速度为20m/s左右,
给定滑移速度时计算的速度场没有震荡,使用的是OpenFOAM中的general方法。
想了解下出现这种问题的可能原因是什么呢~ -
drift-flux model中相对速度计算问题谢谢李老师,我先试试~
-
drift-flux model中相对速度计算问题fluent中提供了不同的滑移速度的计算方法,我想尝试将其移植到OpenFOAM中去,在处理滑移速度的时候遇到了问题:相对速度计算时需要用到相对雷诺数,相对雷诺数计算时需要用相对速度,有点套娃的感觉。这里在计算时需要使用迭代计算的方法么,有没有资料可以参考呢~ 先行谢过各位大佬~
-
OpenFOAM膨胀喷管边界条件设置关联的问题,线没有关联好,放到了wedge类中了
-
compressibleInterFoam中植入drift flux的过程是否正确我先用一维气液激波管试试看看会不会有显著的震荡,谢谢李老师详细的解答~
-
compressibleInterFoam中植入drift flux的过程是否正确了解了,密度梯度那里是个点乘,v1m·grad(rho1),在这种情况下是不是也有可能出现有激波时的震荡啊,因为直接使用了密度的梯度~
-
compressibleInterFoam中植入drift flux的过程是否正确了解了,谢谢李老师~为什么右边的密度梯度是错的呢~
-
compressibleInterFoam中植入drift flux的过程是否正确我大概理解李老师的意思了,我对OpenFOAM中的MULES的原理理解的不透彻,李老师的意思应该是直接基于第一个方程:
直接求解此方程,漂移对流项放到前面去。
我这里处理的时候是想同时使用这两个方程:
最终构造成
,把漂移的影响分开到了源项和对流项中,需要同时修改alphaEqn和alphaSuSp。您觉得这样处理有必要么~
73348540-9b27-449d-9432-26ec0c6fe88b-alphaEqn.H d6ed1b9c-4078-46e3-b212-f7ed1662c6c0-alphaSuSp.H -
compressibleInterFoam中植入drift flux的过程是否正确是的,所以我参照着compressibleInterFoam中的方法对方程做了调整以使最终形式编程\partial alpha/ \partial t+ alpha\div vm的形式
-
compressibleInterFoam中植入drift flux的过程是否正确采用这样处理的原因而不是直接使用$\partial(\alpha*\rho)/\partial t$是想直接套用compressibleInterFoam中的MULES方法,保证alpha的有界性
-
compressibleInterFoam中植入drift flux的过程是否正确可压的场景下,vm和v1m前面的系数不同,我计划参照compressibleInterFoam中的方程的形式来处理组分输运方程,组分输运方程最终的输运形式如下:
完成这个方程的处理之后,按照compressibleInterFoam中的方法,能够得到alpha1Phi和alpha2Phi,刚好能够用于pEqn中的参数构造。想请问下这样是否存在问题~ -
OpenFOAM膨胀喷管边界条件设置是的,二维和三维的能量耗散是不同的,所以结果上会有很大的差异。另外,你可能要重新画个数模以满足OpenFOAM中的wedge类型的使用要求
-
compressibleInterFoam中植入drift flux的过程是否正确这里α1表示分散相,α2表示连续项
-
OpenFOAM膨胀喷管边界条件设置这里应该是有一个明显的设置错误,喷管是轴对称结构,但是你这里算的时候使用的是二维平面结构,应该改为wedge结构
-
compressibleInterFoam中植入drift flux的过程是否正确最近在基于compressibleInterFoam添加漂移通量,想请大佬们帮忙看看过程上是否存在明显的问题:组分输运方程如下:
参照compressibleInterFoam中的处理方法,可以整理为:
两式相加得到div v_m:
在第三个方程中添加±α1 div v_m,相对速度的影响拆分成两部分处理,一部分作为源项参照dgdt的处理方法,一部分作为通量处理:
按照这样的处理方法是否存在明显的问题呢?有没有更好的处理方案呢~ 请各位大佬提提意见~ -
openFoam并行计算的计算效率如何提升@学流体的小明 谢谢啦
-
openFoam并行计算的计算效率如何提升最近在用自己编的求解器运行计算一个300万网格的算例,想请假下怎么从编程的角度提升计算效率,又相关的教程么~
我在网上看到有大佬针对输出显示的优化的介绍:https://www.topcfd.cn/18658/
除了这种还有别的可以参考的么,非常感谢~