动量方程如何加入加压力梯度源项?
-
目的:求解充分发展槽道流动,进出口设置为cyclic周期性边界。通过动量方程中加入流动的驱动力:压力梯度。
求解器:OpenFOAM v3.0+
其他参考:fvOptions 的demo:simpleFoam\channel395,设置的是平均速度Ubar,还不是驱动压力;
目前遇到的问题:通过下面的这3种方法,加入压力梯度源项后,测试了 kEpsilon 和 kOmega 湍流模型,kEpsilon表现正常,kOmega在槽道中心区域算出诡异的速度分布,分析问题可能出现在两方面:1)加源项的方法,2)kOmega湍流模型本身。具体算例,我会在另一个帖子里说明。这篇帖子向大家请教添加源项的方法,下面是我尝试的几种方法,不知道是不是正确,请大家帮忙批评指正 :)
- 尝试1:
createField 中定义一个压力梯度向量 pGrad,变量类型为volVectorField
volVectorField pGrad ( IOobject ( "pGrad", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh );
simpleFoam 作如下改动:
solve(UEqn() == -fvc::grad(p)-pGrad); // UEqn.H U = HbyA - rAtU()*fvc::grad(p) - rAtU()*pGrad; // pEqn.H
- 尝试2:
同样需要在 createField.H 中定义 pGrad,不同的是换一种方法改动
simpleFoam,在 UEqn 中加入 -pGrad:
tmp<fvVectorMatrix> UEqn ( fvm::div(phi, U) + MRF.DDt(U) + turbulence->divDevReff(U) == fvOptions(U) - pGrad );
- 尝试3:
通过fvOptions加入源项,该方法参考 cfdonline,
这样给动量方程加入了 pGrad=10 的压力梯度
momentumSource { type vectorSemiImplicitSource; active on; selectionMode all; vectorSemiImplicitSourceCoeffs { //volumeMode absolute; // specific volumeMode specific; injectionRateSuSp { U ( (10 0 0) 0); } } }
- 尝试1:
-
@李东岳 谢谢你关注我的问题,我又尝试了槽道层流,上面的源项方法没有问题,个人比较偏向第2种,把压力梯度加入到 UEqn 矩阵中,感觉有点像隐式的求解方法。
-
@东岳 老师您好,我想请教几个问题:
- 楼主加入了一个源项体积矢量pGrad,这个矢量,OpenFOAM是如何识别为压力梯度呢?也就是说OpenFOAM怎么知道pGrad就是▽p呢?楼主上面的场的设置,OpenFOAM应该只认识pGrad是一矢量,而不知道它就是代表着▽p吧?
- 楼主尝试的第三种方法,U ((10 0 0) 0);表示什么意思?(10 0 0)不就表示速度矢量值了吗?
若能够获得老师的回复,我将不胜感激!!
-
-
@d-benjamin 请问按照方案二加入压力梯度后,在算例中应该如何给压力梯度赋值呢?
-
@气固两相流 你好,正如东岳老师所说,只要在0文件夹下有p文件,就会按照给的式子进行赋值了
-
@d-benjamin 嗯嗯,谢谢您的回答,直接赋值fvc::grad ( p)是没有问题的
-
你好,你在尝试1中加压力梯度源项,为什么对UEqn.H和pEqn.H都做了修改?
simple算法求NS方程是分两步,UEqn.H是预估步,pEqn.H包括校正步。既然在预估步就加入了源项,为什么在校正步pEqn.H还要添加源项呢?
另外,对于尝试3 fvoptions,这种方法的源项是加在simple求解的哪一步中了?是和尝试2相同吗? -
@我是河滩 你好,他只修改了U的方程,并没有修改p方程,只是定义了一个新的压力梯度场。他上面三个尝试应该都是一样的。有一点我不是很明白,就是尝试3中的volumeMode,两个关键字absolute和specific有什么具体的区别。如果东岳老师可以帮忙解答一下就好了。@东岳
-
-
4/24