CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    请教一下向量场的系数矩阵长什么样子?

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

      添加隐式源项过程中发现了个问题:为何向量场的系数矩阵(fvVectorMatrix/fvMatrix<vector>)返回的对角元素( Foam::tmp <Foam::scalarField > Foam::fvMatrix<Type>::D() const)依然是一个标量场(scalarField)?

      难道fvVectorMatrix的[x]
      是这样排序的?
      \begin{equation}
      \begin{bmatrix}
      u1 \\
      v2\\
      w3\\
      vdots\\
      ui\\
      vi\\
      wi\\
      \end{bmatrix}
      \end{equation}

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

        公式怎么不显示?测试下

        \begin{equation}
        \begin{bmatrix}
        u_1\\
        u_2\\
        ...\\
        u_n
        \end{bmatrix}
        \end{equation}

        线上CFD课程开始报名:http://www.dyfluid.com/class.html

        CFD高性能服务器 http://dyfluid.com/servers.html

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

          我是不是不能编辑帖子0.0:xinlei:

          1 条回复 最后回复 回复 引用
          • Tong
            Tong @李东岳 最后由 编辑

            你的可以显示,可能是我写的有点问题

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

              才发现,mathJax里面\\需要用\\\\替换

              https://github.com/mathjax/MathJax/issues/1301
              http://docs.mathjax.org/en/latest/tex.html#tex-and-latex-in-html-documents

              线上CFD课程开始报名:http://www.dyfluid.com/class.html

              CFD高性能服务器 http://dyfluid.com/servers.html

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

                就比如UEqn,是三个方向的标量组成的三个方程对应的独立的三个系数矩阵,还是三个方向上的标量组合在一起的?

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

                  是三个方向的标量组成的三个方程对应的独立的三个系数矩阵,
                  你可以理解为$u_1$,$u_2$,$u_3$三个独立的方程

                  线上CFD课程开始报名:http://www.dyfluid.com/class.html

                  CFD高性能服务器 http://dyfluid.com/servers.html

                  Tong 1 条回复 最后回复 回复 引用
                  • Tong
                    Tong @李东岳 最后由 编辑

                    @东岳
                    我一直是这么理解的,可是为什么fvVectorMatrix的D()方法返回的是一个标量场,照理说是不是应该是三个独立的系数矩阵中相同位置对角线元素构成的一个向量场吗?

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

                      因为三个分量离散后的矩阵系数是一样的

                      线上CFD课程开始报名:http://www.dyfluid.com/class.html

                      CFD高性能服务器 http://dyfluid.com/servers.html

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

                        想了下,如果不考虑边界条件的情况下的确是一样的,看了下fvMatrix的父类lduMatrix,里面的的upper diag 和lower的确都是scalarField。

                        但是如果我需要添加一个各向异性的隐式源项三个方向的lduMatrix应该就不一样了吧?

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

                          又研究了一下明白了:如果手动添加各项异性的全隐式源项,
                          需要拷贝一份A()并且构造成为张量,然后添加隐式源项,之后再构造HbyA,这样就能将各向异性的源线用全隐的方式添加到面通量上面了。

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