CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新
    1. 主页
    2. WYing
    W
    • 资料
    • 关注 0
    • 粉丝 1
    • 主题 18
    • 帖子 42
    • 群组 0

    WYing

    @WYing

    20
    资料浏览
    42
    帖子
    1
    粉丝
    0
    关注
    注册时间 最后登录

    WYing 取消关注 关注

    WYing 发布的最新帖子

    • 使用sixDoFMotionSolver计算多个patches

      各位老师好,

      我用openfoam4.x - pimpleDyMFoam - sixDoFMotionSolver计算刚体运动,我的几何由body1和body2两部分组成,两个body都是刚体,相连并同时运动,因此也可以看做一个刚体。我目前将其看做一个刚体,因此在dynmaicMeshDict文件中写作:

      dynamicFvMesh dynamicMotionSolverFvMesh;
      
      motionSolverLibs ( "libsixDoFRigidBodyMotion.so" );
      
      solver sixDoFRigidBodyMotion;
      
      sixDoFRigidBodyMotionCoeffs
      
      {
      
      ...
      
      patches (BODY);
      
      ...
      }
      

      计算后能得到刚体上的升阻力,但是我要如何根据现有结果分别得到两个body上的升阻力呢?(如果我理解正确的话,BODY1上的气动力 + BODY2上的气动力 = BODY上的气动力)

      还是说只能将其看做两个刚体并重新计算,即dynamicMeshDict中改为patches (BODY1 BODY2);

      谢谢!

      发布在 OpenFOAM
      W
      WYing
    • RE: fvSolution中PIMPLE算法的参数设置问题

      @李东岳 谢谢李老师的回复!抱歉我的上个残差图绘制有问题,正确的应该是这样

      Residuals_Ur5.png

      从上到下依次是速度残差,压力残差和湍流统计量残差。

      上面那个code是我随意找的,我目前使用的code是:

      PIMPLE
      {
          correctPhi          yes;
          nOuterCorrectors    2;
          nCorrectors         2;
          nNonOrthogonalCorrectors 1;
          turbOnFinalIterOnly  false;
      }
      
      relaxationFactors
      {
          fields
          {
              p                  0.3;
          }
          equations
          {
              "(U|k|omega)"      0.7;
              "(U|k|omega)Final" 1.0;
          }
      }
      

      我主要做了以下修改:
      (1)调整 nOuterCorrectors=2,结果几乎不变
      (2)turbOnFinalIterOnly=false,即每次PIMPLE loop都计算湍流,结果几乎不变
      (3)去掉全部的松弛因子,结果几乎不变。

      在模拟中,CFL保持在0.5-1.5之间,但是和实验结果对的不好。李老师觉得还有可能是哪里出了问题呢?

      发布在 OpenFOAM
      W
      WYing
    • fvSolution中PIMPLE算法的参数设置问题

      各位老师好,我关于fvSolution中PIMPLE算法的参数设置有几点疑问。下面是我找到的一个例子,有以下问题:

      1. nOuterCorrector表示pimple loop的次数,并用residualControl来控制U和p的收敛,请问这里需不需要控制湍流统计量(k,epsilon,omega)的收敛?
      2. moveMeshOuterCorrectors 是weakly(默认)或者strongly coupled FSI,但是找不到更多详尽资料介绍这个参数的含义,请问这个参数对结果一般有什么影响?
      3. consistent 是采取SIMPLEC算法,一般什么情况下使用这个参数?
      PIMPLE
      {
          momentumPredictor        yes;
          correctPhi               yes;
          nOuterCorrectors         1000;
          nCorrectors              1;
          nNonOrthogonalCorrectors 1;
          pRefCell                 0;
          pRefValue                0;
          moveMeshOuterCorrectors  yes;
          consistent               yes;
          
          residualControl
          {
              U 1e-6;
              p  1e-6;
          }  
      }
      
      relaxationFactors
      {
          fields
          {
              p                  0.3;
          }
          equations
          {
              "(U|k|omega)"      0.7;
              "(U|k|omega)Final" 1.0;
          }
      }
      

      Residuals.png 另外,我有一个2D结果湍流统计量k的残差逐渐增大,不太清楚是什么原因造成的,能请各位老师帮忙看下吗?我猜测大概率是PIMPLE算法设置的问题。。。谢谢!

      发布在 OpenFOAM
      W
      WYing
    • RE: 计算钝体涡激振动时发散(pimpleDyMFoam)

      @cresendo 在 计算钝体涡激振动时发散(pimpleDyMFoam) 中说:

      songtao.chen@sjtu.edu.cn

      谢谢!非常感谢您的建议!

      发布在 OpenFOAM
      W
      WYing
    • RE: 计算钝体涡激振动时发散(pimpleDyMFoam)

      @cresendo 我重新查看了source code,aRel的松弛方式确实如您所说,感谢指出错误!我对比了不同aRel的结果,发现Rel = 0.4-0.8都可以正常运行且结果几乎相同,但是aRel >0.8之后会发散。我不确定为什么wikidict上说太小的aRel会造成不准确的响应,至少我的结果并不受aRel取值的影响。相反,恰当的aRel可以使本来发散的case收敛。我之前没有指定过aRel,也就是取默认值1.0,因此比较好奇为什么您会设定aRel=0.4而不是其他值,是根据模拟的经验吗?谢谢!

      发布在 OpenFOAM
      W
      WYing
    • RE: 计算钝体涡激振动时发散(pimpleDyMFoam)

      @cresendo 我按照您的建议,调整了accelerationRelaxation 0.4,计算可以正常进行,CFL数小于1, 流场看上去也正常。

      但有一个小疑问,通过查询dynamicMeshDict文件的官方描述,得知accelerationRelaxation直接减小刚体的加速度,比如在某时间步solver求出的加速度为10m/s^2,如果设置accelerationRelaxation=0.4,那么实际应用的加速为4m/s^2,这是一个很大的松弛。而官方描述中建议的accelerationRelaxation范围是0.9-1,并且提到 “Be careful with this accelerationRelaxation. Too low of a value will mean that the Body does not respond to the fluid forces correctly”。因此我担心设置较小的加速度松弛会不会对结果产生太大影响?因为没有实验结果作为参照,所以无法判断结果的准确性。不知道您一般是如何处理这种误差的?谢谢!

      发布在 OpenFOAM
      W
      WYing
    • RE: 计算钝体涡激振动时发散(pimpleDyMFoam)

      @cresendo 请问您说的dynamicMeshDict里的松弛因子是 accelerationRelaxation 0.9 和 accelerationDamping 0.95 吗?

      发布在 OpenFOAM
      W
      WYing
    • RE: 计算钝体涡激振动时发散(pimpleDyMFoam)

      @cresendo 谢谢您的建议!去掉p松弛之后,算例确实可以计算下去,但是CFL保持在20左右,且速度场和压力场也明显有问题。所以还是不太对,不过还是感谢您的帮助!

      发布在 OpenFOAM
      W
      WYing
    • RE: 计算钝体涡激振动时发散(pimpleDyMFoam)

      @李东岳 谢谢李老师的回复!我尝试使用高版本的OF计算,发现相同设置下算例可以正常计算,不确定这是不是低版本的OF涉及动网格时存在bug。不过有个case在所有版本OF下都发散。。。能麻烦李老师有空的话帮忙看一下吗?算例文件case.zip (大约在time=1.43发散)十分感谢!!

      发布在 OpenFOAM
      W
      WYing
    • RE: 计算钝体涡激振动时发散(pimpleDyMFoam)

      @李东岳 感谢您的回复! crash.zip 这是我的case文件,由于大小限制,文件夹内没有包含最后一个时间步,但是大概在time=1.8所有停止计算。欢迎您提出的任何意见或建议!

      发布在 OpenFOAM
      W
      WYing