导入外界网格造波初始化时入口处有问题

DY大世界

将外界划分好的二维长方形网格导入进OF（1906版）中，用内置的stokesII工具包造波，设置好流体区域以及边界条件后（按照原tutorials里面案例边界条件的设置）。但是不知道为什么，在inlet处出现了如下图所示的现象，水域部分（红色）发生了向上的偏移，导致计算失败。恳请了解的大神给予指教，感谢！

李东岳

看起来是边界条件设置。你手动把边界处的alpha设置为0就行了啊

DY大世界

感谢李老师的回复，我是将alpha.water设置成如下形式，但是就是会出现上述问题，其余边界都没事，唯独入口处。

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            waveAlpha;
        value           uniform 0;
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            zeroGradient;
    }

    frontAndBackPlanes
    {
        type            empty;
    }

    top
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
}

// ************************************************************************* //

DY大世界

@东岳感谢李老师的回复！我是将alpha.water设置成如下形式，但是就是会出现上述问题，其余边界都没事，唯独入口处。

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            waveAlpha;
        value           uniform 0;
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            zeroGradient;
    }

    frontAndBackPlanes
    {
        type            empty;
    }

    top
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
}

// ************************************************************************* //

李东岳

waveAlpha这个在官方版没有。不过类似这种的value值不成问题，是个dummy值，具体值是多少是通过waveAlpha计算的

DY大世界

@东岳要是使用OF1906的blockMesh划分网格就可以正常计算，我用fluent导入的网格但是边界条件都没有改动就会出现这样的问题，很奇怪。李老师会不会就是因为软件的兼容问题不可避免的出错，也不会找的到解决方法呢。

DY大世界

@东岳李老师您好，不好意思又打扰您了。我调试了，问题就是出在了这个waterAlpha上，但是又找不到waterAlpha这个文件，也更无法修改这个文件，请问这个应该怎么解决呢？

队长别开枪

@DY大世界悄悄告诉你，g方向必须是z方向，也就是你的二维网格必须在y-plane。

DY大世界

@队长别开枪十分感谢大神的回复！但是对于从外界导入的二维模型，OF是不是均默认在z方向进行一定的拉伸，因为我发现无论我导入什么样的二维模型进入OF，OF均是在Z方向上进行拉伸厚度的。这点无法改变吗？否则g方向与永远不在z方向上。
QQ截图20191114153838.jpg

队长别开枪

@DY大世界 https://openfoamwiki.net/index.php/TransformPoints

DY大世界

@队长别开枪谢过队长！队长好人！

nbyjn

@队长别开枪
大神
之前看你发过那个开源PLIC求解器
我想问一下
VOF几何方法为什么要多一步界面重构
是为了确定界面位置，
从而在下一步advection的时候用几何信息计算通量吗？
那比如isoAdvetor这种直接构建等值面来advection，不需要上一步的界面位置信息，是不是就不需要这个界面重构的步骤了呢？
谢谢队长！

yfclark

@nbyjn 界面重构是为了消除界面的数值耗散，vof都需要重构，只不过代数重构不需要重构准确的界面位置