CFD中文网

@cyw 凝结和气化之差不就已经是体积源了吗，为啥还加vDotcAlpha

@cyw 为啥mDot第二项没负号，你少写了吗，还有就是不就是再乘了个体积分数吗？看对不对得对照它的方程看

@cyw 代码里只有mDotAlpha（）质量源项啊，没有mDot

@东岳 在新版本5.0或者V1712的interfoam求解器内都有这个头文件

其中zeroField Su；
zeroField Sp；
zeroField divU；
都是零场的作用不清楚，是为了用户以后在这里自行修改吗？

如何理解最后两项？？？
是fvm::Sp(vDotvmcAlphal,alpha1)+vDotcAlphal 不是单一的vDotvmcAlphal，自己看代码好像vDotvmcAlphal已经代表了m*(1/rhol-alpha1*(1/rhol-1/rhov)),这不就是公式变形后的最后一项吗？求老师、大神指导

fluent中设置是这样的：
Boundary A, B and C, D
are set as periodic boundary conditions. Boundary E is set as no slip
wall. Boundary F is set as opening pressure boundary.
不知道在OF怎样设置好

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 100000;

boundaryField
{
    inlet
    {
        type cyclic;
    }
    gasinlet
    {
        type cyclic;
    }
    outlet1
    {
        type cyclic;
    }
    outlet2
    {
        type cyclic;
    }
    left
    {
       type cyclic;
    }
    right
    {
        type cyclic;
    }
    TWall
    {
       type            fixedFluxPressure;
    }
    
    atmosphere
    {
        type            totalPressure;
        p0              uniform 100000;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           1;
        value           uniform 100000;
    }
}
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type cyclic;
    }
    gasinlet
    {
       type cyclic;
    }
   outlet1
    {
       type cyclic;
    }
    outlet2
    {
       type cyclic;
    } 
    left
    {
      type cyclic;
    }
    right
    {
      type cyclic;
    }
    TWall
    {
       type            fixedValue;
       value           uniform (0 0 0);
    }
    
    atmosphere
    {
       type            pressureInletOutletVelocity;
       value           uniform (0 0 0);
    }
}

体积项用setField，但这样不对，求大神帮忙:big_mouth:

@wwzhao 谢谢

@小龙 filter 搜一下就可以找到

@random_ran 好的
这是没切之前的，网格很正常

slice之后

而且cell变成了正常截面的两倍

为啥paraview做了slice后，对应的cell数目变成了两倍，而且选择surface and edge 显示，发现网格中间有斜杠了

@mohui 我是基于那个方程进行修改添加一些，只不过它本身存在一些问题，当我设置1的时候温度正常

如果不设置nAlphaSubCycles=1的话，会出现超过最大温度，数值很越来越大，不稳定情况，设置了就正常了，增加温度方程的方法网上的interTempFoam，有没有遇到同样的状况，请大神分析一下

@李东岳 李老师，是在说paraview里选取网格的point嘛？

想知道网格在模型里按照什么顺序来排列的

@李东岳
1.

volScalarField Source
(
    IOobject
    (
                "Source",
		runTime.timeName(),
		mesh,
		IOobject::NO_READ,
		IOobject::AUTO_WRITE
    ),
    mesh
);

加了，mesh才能通过编译，
2. Source[celli]是double类型的，是不是类似于对应网格的Source.value（），
3. volScalarField又是定义单位的呢，是不是在速度方程后添加==Source；并且保证量纲统一就可以达到更新速度的目的
4. 我还不是很理解volScalarField和scalarField的区别

李老师我问题比较多，看您方便回答哪个，十分感谢

@李东岳 李老师我的是3.0.X版本的，
提示：
UEqn.H:68:6: error: ‘Source’ was not declared in this scope
Source 没有指定，这个要怎么创建啊

@李东岳 那要怎么写呢，李老师能请简单指点一下:happy:

@李东岳 李老师，就是添加源项的速度方程只想对于某种网格求解以此来更新这个网格内的速度

对于单个网格如何写速度方程，想更新单个网格的速度，UEqn好像是求解了所有的