请教各位大佬，下图公式中的三项怎么在OpenFOAM或者其他后处理软件中导出呢？

[strive_DUT]

如下图红框所示，我想在OpenFOAM中输出这三项所代表的场数据，请问该如何实现呢？跪求

[李东岳]

老铁我记得你上过编程课啊。分别是

fvc::div(U)*w
(fvc::grad(rho) ^ fvc::grad(rho))/sqr(rho)
curl(fvc::div(tau)/rho)

[strive_DUT]

@李东岳 哈哈是的，李老师。我自己提前写了一下，担心不准确，来请教一下各位。
（1）体积膨胀效应中的涡度omega是不是应该定义为fvc::curl(U)呀？
（2）斜压效应中的原式是▽ρ×▽p(密度梯度和压力梯度)，fvc::grad(rho) ^ fvc::grad(p)中间是叉乘吧？
（3）粘性效应中的粘性应力张量tau是不是应该定义为tau=-muEff()*dev(twoSymm(fvc::grad(U_)))呢？
关于涡量和粘性应力张量的定义是要在求解器中重新声明吧？关于上述三个效应的量化，我不确定能否在计算结果中直接导出呢？还请李老师指点。

[尚善若水]

omega 直接 postProcess Ω: vorticity？

[strive_DUT]

@尚善若水 后处理可以直接导出涡度，但是我的计算中涡度是个中间量，所以我写出了计算公式，再把它加入到求解器中了

[尚善若水]

@strive_DUT 明白你的意思了，我最近也在处理相同的问题，也是参考一些后处理中的代码实现。

[李东岳]

你俩这方程打眼一看，简直一模一样。不会是同一个实验室的吧。哇咔咔

[strive_DUT]

@李东岳 是的李老师，都是涡量输运方程，但不是一个实验室哈哈哈