CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

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

    OpenFOAM
    5
    13
    1726
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • DY大世界
      DY大世界 最后由 编辑

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

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

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

        2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        DY大世界 1 条回复 最后回复 回复 引用
        • DY大世界
          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;
              }
          }
          
          // ************************************************************************* //
          
          1 条回复 最后回复 回复 引用
          • DY大世界
            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;
                }
            }
            
            // ************************************************************************* //
            
            1 条回复 最后回复 回复 引用
            • 李东岳
              李东岳 管理员 最后由 编辑

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

              2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
              CFD高性能服务器 http://dyfluid.com/servers.html

              DY大世界 2 条回复 最后回复 回复 引用
              • DY大世界
                DY大世界 @李东岳 最后由 编辑

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

                1 条回复 最后回复 回复 引用
                • DY大世界
                  DY大世界 @李东岳 最后由 编辑

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

                  队长别开枪 1 条回复 最后回复 回复 引用
                  • 队长别开枪
                    队长别开枪 教授 @DY大世界 最后由 编辑

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

                    DY大世界 1 条回复 最后回复 回复 引用
                    • DY大世界
                      DY大世界 @队长别开枪 最后由 编辑

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

                      队长别开枪 1 条回复 最后回复 回复 引用
                      • 队长别开枪
                        队长别开枪 教授 @DY大世界 最后由 编辑

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

                        DY大世界 N 2 条回复 最后回复 回复 引用
                        • DY大世界
                          DY大世界 @队长别开枪 最后由 编辑

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

                          1 条回复 最后回复 回复 引用
                          • N
                            nbyjn @队长别开枪 最后由 编辑

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

                            1 条回复 最后回复 回复 引用
                            • Y
                              yfclark 讲师 最后由 编辑

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

                              1 条回复 最后回复 回复 引用
                              • First post
                                Last post