CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    关于非均一扩散系数的扩散项离散问题

    Algorithm
    2
    4
    1069
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 肖艺
      肖艺 最后由 编辑

      1ed57deb-404c-4705-98ba-70691f34a6fb-image.png
      fbbd9b40-b208-43a8-a080-ee0c298ab491-image.png
      采用二阶精度有限体积法,扩散系数是关于空间位置的函数,面心扩散系数直接给定扩散系数的解析值,扩散项用中心差分格式,时间项用一阶欧拉隐式,计算结果是变量在全域的总积分值会逐渐减少,令扩散系数为常数时计算结果是守恒的。
      这种扩散系数应该怎么离散啊。。。
      大佬们指点下:zoule:

      沧海一声笑

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

        已经有$\alpha$解析解了还需要离散么?

        守恒是什么守恒?全局的$\phi$守恒?

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

        1 条回复 最后回复 回复 引用
        • 肖艺
          肖艺 最后由 编辑

          单元体积乘单元变量值(在全计算域所有单元求和),这个积分值和初始场的积分值相等,既是这里指的守恒

          沧海一声笑

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

            如果不考虑时间项
            \begin{equation}
            \int^{f+1}_ {f} \frac{\partial}{\partial x}\left(\alpha \frac{\partial \phi}{\partial x}\right)\mathrm{d}x=\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f}
            \end{equation}
            把后面两项$\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_{f}$把所有的面连起来,看起来是守恒的,剩下的只有边界条件

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

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