【分享+搬运】自定义非均匀inlet U



  • 有条件的可以自行油管:https://www.youtube.com/watch?v=K-nAF3qAPTc
    感谢youtuber:Mark Kimber
    侵删

    1. 首先你要有一个自己的进口速度场的公式

      0_1545211435195_151b6040-02bc-4564-b234-278bda56300f-image.png

      或者自行拟合一个,根据实际情况来。

      0_1545211492714_342c9ffa-58c9-4586-a02d-c67091a23998-image.png

    2. 在Terminal中输入

      writeCellCentres
      

      获取网格每个cell的坐标,会生成ccx,ccy,ccz三个文件,内涵每个patch中的每个cell的坐标(x,y,z),找到inlet patch,我的是topWall。3个文件中的topwall(就是我的inlet)的2100代表的是这个patch上又2100个网格。同理对于ccy文件同样找到这个patch的所有cell坐标。分别复制这些数据,生成两个新的文本,叫inletX, inletY。根据实际网格情况,以及坐标轴自行选择用哪2个文件。

      0_1545211861236_d4842d0a-141a-4d61-ba47-720c64fbded6-image.png
      0_1545212222778_b5c53c18-efcc-4912-b38e-f0298d2a48a0-image.png

    3. 进入octave(一个类似matlab的程序,youtuber是这么说的)这一步主要是根据你的速度场的公式,xy坐标,计算出进口的vector。
      0_1545212556614_ba95eede-affa-46b6-8d20-5754edb8bdf5-image.png

      x= load('inletX');   //读取inlet patch的x坐标
      y=load('inletY');    //读取inlet patch的y坐标
      Uz==1.1+0.18988/0.01415/sqrt(3.1415926/2)*exp(-2*(x.^2+y.^2)/0.1415.^2);  //公式
      fid= fopen('inletU','w')// inletU文件要提前新建好,不知道为什么视频中输出结果为5,我的结果一直为3,不过不影响使用
      N=length(x)
      
      for i=1:N
      fprintf(fid,'(%8.4f %8.4f %8.4f) \n',0,0,Uz(i));
      end
      
      fclose(fid)
      
      

      0_1545217781545_7be93b21-67bc-454a-94e0-ebcd256ac948-image.png

      0_1545215405217_03c8ce27-6a6c-4acc-abac-a45bbef4ec72-image.png

      其实第3步完全可以用excel做。。。。。。。

    4. 把得到的vector的值带入到0/U中
      0_1545218022098_fb103f76-1896-4296-a840-d3f675d77abd-image.png

    5. 进入paraview中查看即可

      0_1545218042261_365addff-5093-4e67-9908-c49de0f74690-image.png


  • 管理员

    感谢分享,楼主好人



  • 想咨询下, 这个实现的功能跟setFields有什么区别?比如说这个方法设置的入口面速度值是一直固定的?而etFields是否只是个初始值,但会随时间改变?



  • @aiweimo setfields 应该是定义cell centers 的值吧(这个我不是太确定,如果错误希望各位大神指正),这个应该是定义的入口边界处的,还有那个writeCellCentres现在应该是在postProcessing中了,不能直接用。如果你想定义随时间改变的inlet,可以参考一下timeVaringMappedFixedValue(也许拼错了哈哈)。具体实现的例子,可以参考一个叫eddylicious的工具包,或者查看openfoam的代码,或者去cfdonline上查看对应的例子。


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1