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



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



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



  • 感谢李老师的回复,我是将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;
        }
    }
    
    // ************************************************************************* //
    


  • @东岳 感谢李老师的回复!我是将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计算的



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



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


  • 自由表面模型副教授 OpenFOAM教授

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



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


  • 自由表面模型副教授 OpenFOAM教授



  • @队长别开枪 谢过队长!队长好人!



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


  • OpenFOAM讲师

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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2