fvc::DDt是什么意思？

yfclark

我看到有人算

用了如下的表达式

fvc::DDt(phi,p)

我试着写入代码里，可以运行，但是在网上搜索了半天也没有发现这部分代码的定义。有老哥看到过这个运算的定义吗？
Openfoam 2.3.0

zym604

看上去像是ddt+div
查了一下代码:
https://github.com/OpenFOAM/OpenFOAM-3.0.x/blob/master/src/finiteVolume/finiteVolume/fvc/fvcDDt.C

template<class Type>
tmp<GeometricField<Type, fvPatchField, volMesh> >
DDt
(
    const surfaceScalarField& phi,
    const GeometricField<Type, fvPatchField, volMesh>& psi
)
{
    tmp<GeometricField<Type, fvPatchField, volMesh> > ddtDivPhiPsi
        = fvc::ddt(psi) + fvc::div(phi, psi);

    if (phi.mesh().moving())
    {
        return ddtDivPhiPsi - fvc::div(phi + phi.mesh().phi())*psi;
    }
    else
    {
        return ddtDivPhiPsi - fvc::div(phi)*psi;
    }
}

这是3.0的版本，不确定2.3是不是也是如此。从这个版本里看，就是ddt+div。

至于定义，应该是随体导数（Material derivative）。

bestucan

6.0版的定义。估计之前的难理解，又加了个U。

李东岳

感谢诸位大佬分享，赠送66元漏油奔驰代金券

yfclark

谢谢各位大佬

strive_DUT

@zym604 随体导数最后不是梯度项吗？源代码中两个值作差都是散度呀？这个该怎么理解？

xpqiu

@strive_DUT
根据张量等式
$$
\nabla \cdot ( \mathbf{U} \phi) = \mathbf{U} \cdot \nabla \phi + \phi (\nabla \cdot \mathbf{U})
$$
所以，
$$
\frac{\partial \phi}{\partial t} + \mathbf{U} \cdot \nabla \phi = \frac{\partial \phi}{\partial t} + \nabla \cdot ( \mathbf{U} \phi) - \phi (\nabla \cdot \mathbf{U})
$$
这里等式右边的减号就是二楼贴出来代码里面减号的来源了。

另外，

@bestucan 在 fvc::DDt是什么意思？中说：

6.0版的定义。估计之前的难理解，又加了个U。

这个是来搞笑的么？

strive_DUT

@xpqiu 多谢多谢！这个张量等式我还真不了解！长知识了

bestucan

@xpqiu 是啊，我的签名就声明了我就是来搞笑的滚来滚去，滚去滚来，薛定谔的滚蛋状态布朗滚蛋运动

像您这样的能答上来的人没有出现，只有先搞笑一下，不至于别人问的问题没人搭理。还能刺激一下大家的回答欲。忘了在哪看的研究，一般不愿意回答的问题要是有个错误回答，就会引来很多纠错。因为纠错的成就感更强一些。

李东岳

小弟只能在这里多谢各位捧场了，最近事情太多了，早五晚八的，哎...每天就是在瞎忙..

bestucan

@东岳老师言重了，我从东岳流体上学到很多很多，现在来论坛逛也是回报，不算付出。虽然东岳流体上的内容不收费，那是老师的好，我白占便宜我受不了。

对，我是来论坛打工换东岳流体上的内容的

xpqiu

@bestucan 在 fvc::DDt是什么意思？中说：

@xpqiu 是啊，我的签名就声明了我就是来搞笑的滚来滚去，滚去滚来，薛定谔的滚蛋状态布朗滚蛋运动

像您这样的能答上来的人没有出现，只有先搞笑一下，不至于别人问的问题没人搭理。还能刺激一下大家的回答欲。忘了在哪看的研究，一般不愿意回答的问题要是有个错误回答，就会引来很多纠错。因为纠错的成就感更强一些。