CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    interFoam的曲率计算

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

      void Foam::interfaceProperties::calculateK()
      {
          const fvMesh& mesh = alpha1_.mesh();
          const surfaceVectorField& Sf = mesh.Sf();
      
          // Cell gradient of alpha
          const volVectorField gradAlpha(fvc::grad(alpha1_, "nHat"));
      
          // Interpolated face-gradient of alpha
          surfaceVectorField gradAlphaf(fvc::interpolate(gradAlpha));
      
          // gradAlphaf -=
          //    (mesh.Sf()/mesh.magSf())
          //   *(fvc::snGrad(alpha1_) - (mesh.Sf() & gradAlphaf)/mesh.magSf());
      
          // Face unit interface normal
          surfaceVectorField nHatfv(gradAlphaf/(mag(gradAlphaf) + deltaN_));
          // surfaceVectorField nHatfv
          // (
          //     (gradAlphaf + deltaN_*vector(0, 0, 1)
          //    *sign(gradAlphaf.component(vector::Z)))/(mag(gradAlphaf) + deltaN_)
          // );
          correctContactAngle(nHatfv.boundaryFieldRef(), gradAlphaf.boundaryField());
      
          // Face unit interface normal flux
          nHatf_ = nHatfv & Sf;
      
          // Simple expression for curvature
          K_ = -fvc::div(nHatf_);
      
          // Complex expression for curvature.
          // Correction is formally zero but numerically non-zero.
          /*
          volVectorField nHat(gradAlpha/(mag(gradAlpha) + deltaN_));
          forAll(nHat.boundaryField(), patchi)
          {
              nHat.boundaryField()[patchi] = nHatfv.boundaryField()[patchi];
          }
      
          K_ = -fvc::div(nHatf_) + (nHat & fvc::grad(nHatfv) & nHat);
          */
      }
      

      interFoam在计算曲率的时候应该是调用了这个函数, 我想问一下这里被注释掉的这段complex expression什么时候可以用? 有没有原始文献可以参考。

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

        或许可以把方程写出来对比一下嘛?

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

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

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