Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    怎么编写法向二阶梯度为0的边界条件

    OpenFOAM
    2
    7
    2501
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • Y
      yang山青 last edited by

      在计算中需要用到连续性边界条件,即法向二阶梯度为0。在常用的边界条件中能看到gradientInternalCoeffs()和gradientBoundaryCoeffs(),应该能够通过将这个一阶梯度变为二阶梯度来实现这个边界条件,但是不知道这个函数在哪定义?或者有没有其他的实现这个边界的方法?谢谢

      1 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by

        法向二阶梯度为0

        这不就是法向一阶梯度为固定值么?试试fixedGradient边界条件?已经有了。

        线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        Y 1 Reply Last reply Reply Quote
        • Y
          yang山青 @李东岳 last edited by

          @李东岳 并不是一阶为固定值啊,因为在每个网格上是定梯度,但是在整个边界上并不是所有网格都是相同的梯度。

          1 Reply Last reply Reply Quote
          • 李东岳
            李东岳 管理员 last edited by

            嗯,好像是。好吧,我们缕一缕。

            考虑下面这个方程:
            \begin{equation}
            \frac{\rd}{\rd x}\left(\frac{\rd T}{\rd x}\right)=0
            \end{equation}
            离散后有:
            \begin{equation}
            \left(A\frac{\rd T}{\rd x}\right)_e-\left(A\frac{\rd T}{\rd x}\right)_p=0
            \end{equation}
            考虑边界点$p$,如果是fixedGradient,那么就是
            \begin{equation}
            \left(A\frac{\rd T}{\rd x}\right)_p=A*\mathrm{someValue}
            \end{equation}
            如果是二阶fixedGradient,那么就是
            \begin{equation}
            \left(A\frac{\rd T}{\rd x}\right)_p=A*\mathrm{someValueSec}
            \end{equation}
            其中$\mathrm{someValue}$ (fixedValue)在所有边界面初都一样,$\mathrm{someValueSec}$在每个边界面上都不一样,是这样的吧?

            线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
            CFD高性能服务器 http://dyfluid.com/servers.html

            Y 2 Replies Last reply Reply Quote
            • Y
              yang山青 @李东岳 last edited by

              @李东岳 根据方程来看,是这样的。从物理上可以认为是一个内边界。

              1 Reply Last reply Reply Quote
              • Y
                yang山青 @李东岳 last edited by

                @李东岳 您认为如果在边界周围将网格密度增大,然后使用zeroGradient边界条件能近似看成连续性边界吗?

                1 Reply Last reply Reply Quote
                • 李东岳
                  李东岳 管理员 last edited by

                  在每个网格上是定梯度

                  你可以试试修改fixedGradient边界条件,这里面有一个gradient_成员,默认是均一的需要用户给定,按照你的定义,好像是应该把这个成员设置为和网格边界毗连网格单元的定梯度。

                  线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
                  CFD高性能服务器 http://dyfluid.com/servers.html

                  1 Reply Last reply Reply Quote
                  • First post
                    Last post

                  CFD中文网 | 东岳流体 | 京ICP备15017992号-2
                  论坛登录问题反馈可联系 li.dy@dyfluid.com