如何看流场中的加速度云图？

DY大世界

在controldict文件中添加了grad(U)字典了，想看加速度场的云图，这里grad(U)有8个magnitude，这是表示什么意思，想看加速度应该选择哪一个？

东岳

@DY大世界在如何看流场中的加速度云图？中说：

想看加速度应该选择哪一个？

gradU是二阶张量，9个分量。
什么是加速度...

DY大世界

@东岳想看一下加速度云图应该怎么设置

DY大世界

@东岳感谢李老师的解答

东岳

额，不好意思，没太看懂你的加速度是什么？加速度是怎么计算的啊？grad(U)不是加速度

史浩

@DY大世界加速度是速度对时间的导数，应该是ddt(U)。grad(U)是空间的导数，相当于速度的梯度

东岳

@史浩谢谢谢谢，我在NS方程笔记里面还写过，自己都给忘了嗯，就是$\bfU$的物质导数

史浩

@东岳你这一提醒，我发现我上面的好像错了。这里说的加速度应该指的是物质导数，物质导数和局部导数有如下关系

我上面写的fvc::ddt(U)只是其中的局部导数。
那问题来了，在OpenFOAM中怎么求物质导数？下面的公式求物质导数看起来有点土

DUDt=fvc::ddt(U)+fvc::div(phi, U)

Yu_Tian

@史浩您好，同求如何查看电场中的电场线云图，是对电场强度Ue求梯度，grad(U),我想问下这个具体操作应该是怎么写进例子中，在哪个文件中补充的具体命令是什么？可不可以给个样板，谢谢

史浩

@Yu_Tian 你说的是电场的streamline是么？是下图这种效果么？

这个直接用电场强度矢量做流线图就可以得到，不用再求导了

Yu_Tian

@史浩对对对，是想做出这个图，但是我采用的openfoam2.1.1版本的interFoamEHD求解器，里面的文件只有这些，其中有电压Ue，没有电场强度E,所以想通过对电压求解，是不是应该求梯度

Yu_Tian

@史浩或者请教一下，这个图里的电场强度是怎么得到的？有没有例子可以分享一下

史浩

@Yu_Tian 电场强度是电势（电压）的梯度

E = fvc::grad(Ue)  // Electric field intensity
j = sigma*E        // Current density, sigma--electric conductivity

Yu_Tian

@史浩感谢回答，受到您的提示后，我采用后处理函数postProcess -func "grad(Ue)",确实得到了电场强度，但是我做了一下电场强度的矢量图，好像方向不太对，如下图所示，背景的云图是电势的分布，上面红色为正电极，下面蓝色为0电位，电场强度的方向不应该是从上指向下面吗？现在方向是不是反了？要不要加一个负号？但是postProcess -func "-grad(Ue)"这种表达又不满足形式要求，这个应该怎么做呢？谢谢

史浩

@Yu_Tian 抱歉，前面的公式少打了一个负号-，不过你已经找到这个错误了。
postProcess这个东西我还没用过，有空研究研究，你看看下面这种书写格式行不行

postProcess -func "grad(-1.0*Ue)"

当时我直接写在程序里面了，把他当成一个结果输出了

Yu_Tian

@史浩感谢回答，但是这样也不行，我也在好奇说这个postProcess函数能不能进行这种加减乘除的运算？或者正确的表达方式是什么样的？不知道您上面提到的那种方式代码是在哪里使用的呢？是不是可以采用您的方式实现呢？

闻久STU

grad(U)是除以距离，应该看ddt吧

DY大世界

@闻久STU 收到~感谢大神指点！

史浩

@Yu_Tian 在createFields.H文件中，添加变量声明

volVectorField js
(
    IOobject
    (
        "js",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    -sigma*fvc::grad(Ue)
);

在求解完Ue后，更新一下js

js=-sigma*fvc::grad(Ue);

我一般这样通过求解器实现的。
也可以用postProcess实现，但是需要具体看一下他的源码，这个我不是很熟悉
放一个我之前写过的一个求边界“axis”上的最大速度，这个是参考之前版本的calc方法写的，后来calc全部整合到postProcess中了，但是在OpenFOAM 6里面，按照calc的思虑一样可以用

#include "timeSelector.H"
#include "calcType.H"

void calc(const argList& args, const Time& runTime, const fvMesh& mesh)
{
    IOobject UHeader
    (
        "U",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ
    );

    if (phiHeader.headerOk())
    {
		volVectorField U(UHeader, mesh);
		forAll(U.boundaryField(), patchi)
		{
			const fvPatchVectorField & pU = U.boundaryField()[patchi];
			const fvPatch & pU_B = pU.patch();
			if(pU_B.name() == "axis")
			{
				Foam::Info<< "mag(U) max : " << max(mag(pU)).value() << Foam::endl;
				break;
			}
		}
	}
	else
	{
		Foam::Info<< "No U exists!" << Foam::endl;
	}

    Info<< "\nEnd\n" << endl;
}

int main(int argc, char *argv[])
{
    #include "setRootCase.H"
    #include "createTime.H"
    Foam::instantList timeDirs = Foam::timeSelector::select0(runTime, args);
    #include "createNamedMesh.H"


    forAll(timeDirs, timeI)
    {
        runTime.setTime(timeDirs[timeI], timeI);
   
        Foam::Info<< "Time = " << runTime.timeName() << Foam::endl;
   
        mesh.readUpdate();
   
        calc(args, runTime, mesh);
   
        Foam::Info<< Foam::endl;
    }
	
    Foam::Info<< "End\n" << Foam::endl;
   
    return 0;
}

Yu_Tian

@史浩好的，感谢大神我仔细研究一下