CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    定义加速度的phi

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

      重力g的phi可以写成phiG=g & mesh.Sf(),其中g是uniformDimensionedVectorField。如果g变成一个空间上不均匀的加速度场,类型变成volVectorField,这时候的加速度的phi该如何写??

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

        :wocao: 那不就跟速度U一样了么?也是这样写啊?

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

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

        W 1 条回复 最后回复 回复 引用
        • W
          wwzhao 教授 @李东岳 最后由 编辑

          @东岳 Sf是surfaceVectorField,不能跟volVectorField直接相乘,需要将volVectorField插值到面上再相乘。

          李东岳 D 2 条回复 最后回复 回复 引用
          • 李东岳
            李东岳 管理员 @wwzhao 最后由 编辑

            @wwzhao 主楼知道的吧?

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

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

            1 条回复 最后回复 回复 引用
            • D
              danao @wwzhao 最后由 编辑

              @wwzhao 那应该怎么写啊啊。。。

              1 条回复 最后回复 回复 引用
              • D
                danao @游客 最后由 编辑

                @东岳 @wwzhao ,抱歉,忘了格式了:threaten:

                volVectorField g
                (
                IOobject
                (
                "g",
                runTime.timeName(),
                mesh,
                IOobject::NO_READ,
                IOobject::NO_WRITE
                ),
                g=关于网格空间坐标xy的函数
                );
                dimensionedVectorField g = fvc::interpolate(g,"g");
                phig=g & mesh.Sf()
                1 条回复 最后回复 回复 引用
                • 李东岳
                  李东岳 管理员 最后由 李东岳 编辑

                  volVectorField g
                  (
                  IOobject
                  (
                  "g",
                  runTime.timeName(),
                  mesh,
                  IOobject::NO_READ,
                  IOobject::NO_WRITE
                  ),
                  g=关于网格空间坐标xy的函数
                  );
                  //dimensionedVectorField g = fvc::interpolate(g,"g");
                  phig = fvc::interpolate(g) & mesh.Sf()
                  

                  你这是做什么?怎么会有非均一重力场

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

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

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

                    @东岳 非惯性坐标系啊,体积力与每个网格的具体位置有关,这样就不用动网格了

                    1 条回复 最后回复 回复 引用
                    • Referenced by  P Prometheus10 
                    • First post
                      Last post