PISO算法讨论

yjc95

@李东岳
 @李东岳
老师，这是pisoFoam里的速度方程

    fvm::ddt(U) + fvm::div(phi , U)
+  MRF.DDt(U)
+  turbulence->divDevSigma(U)
==
    fvModels.source(U)

solve(UEqn == -fvc::grad(p))

这是icoFoam里的速度方程

    fvm::ddt(U) 
+ fvm::div(phi , U)
-  fvm::laplacian(nu , U)

solve(UEqn == -fvc::grad(p))

比较不难发现，pisoFoam里的速度方程多了MRF.DDt(U)，turbulence->divDevSigma(U)
和== fvModels.source(U)，少了fvm::laplacian(nu , U)。

MRF.DDt(U)是和旋转结构有关，在cavityRAS这个案例中用不到。

1.

turbulence->divDevSigma(U)应该是这个案例中用到的K-epsilon方程，这个方程的具体代码在哪个文件里可以看到呢？fvModels.source(U)就是fvm::laplacian(nu , U)吗？

2.

第一段代码中已经有==了，第二段代码又有一个==，从逻辑上讲，应该将fvModels.source(U)移到等号左边，再等于-fvc::grad(p)，这是C++语法特点吗？

Afr1yne

@yjc95

icoFoam是不能用湍流模型的，fvm::laplacian(nu , U)是动量方程右边的粘性项。pisoFoam用湍流模型把粘性项模化turbulence->divDevSigma(U)。fvModels.source(U)是源项。
2.solve(UEqn == -fvc::grad(p))是momentumPredictor里面的，不开的话是用不上的，我感觉你是不是把它理解成求解了？

yjc95

@Afr1yne 在 PISO算法讨论中说：

@yjc95

icoFoam是不能用湍流模型的，fvm::laplacian(nu , U)是动量方程右边的粘性项。pisoFoam用湍流模型把粘性项模化turbulence->divDevSigma(U)。fvModels.source(U)是源项。
2.solve(UEqn == -fvc::grad(p))是momentumPredictor里面的，不开的话是用不上的，我感觉你是不是把它理解成求解了？

谢谢您的回复。
1.“pisoFoam用湍流模型把粘性项模化turbulence->divDevSigma(U)”具体是咋模化的呢？
2.“fvModels.source(U)是源项”是哪个源项？是体积力、电场力这种？
3.假设我打开了动量预测，这俩==该咋理解呢？

Afr1yne

1.就是湍流模型
2.你可以理解为是个接口，具体是什么根据你的设置，你可以去找下tutorials里面fvModel文件看下
3.去了解下fvVectorMatrix UEqn，不是你想象的这种连等

    fvm::ddt(U) + fvm::div(phi , U)
+  MRF.DDt(U)
+  turbulence->divDevSigma(U)
==
    fvModels.source(U)
== -fvc::grad(p)

请教下这种红字怎么打出来的

yjc95

@Afr1yne
1.我知道是湍流模型，我想在openfoam的代码中找到它
3.我了解过FvVectorMatrix，说它存储了Ax=b的相关信息。我也能猜到不是连等，我只是想看看最底层的定义，但我不知道在哪里找。
在你想标红的字前后加上`

李东岳

turbulence->divDevSigma(U)应该是这个案例中用到的K-epsilon方程，这个方程的具体代码在哪个文件里可以看到呢？

看笔记

“fvModels.source(U)是源项”是哪个源项？是体积力、电场力这种？

所有源项。普适性源项

我了解过FvVectorMatrix，说它存储了Ax=b的相关信息

==是操作符重载，笔记里面写了操作符重载

yjc95

谢谢岳子老师。

yjc95

@李东岳岳子老师，先说明下，我不是在故意挑刺，就是看到有些笔误的地方真心难受，很想把它改过来。
不可压+瞬态里，方程29是方程28+19得到的，等号右边的$\frac{1}{V_{p}}$应该是加了两遍，那么应该是$\frac{2}{V_{p}}$吧，同理，下面的几个方程也是。
不可压+稳态里，方程27以及后面几个方程也是同样的。

李东岳

方程29没问题啊

yjc95

@李东岳老师我看错了，不好意思

CFD中文网

PISO算法讨论

1.

2.