CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新
    1. 主页
    2. 星星星星晴
    星
    • 资料
    • 关注 7
    • 粉丝 13
    • 主题 51
    • 帖子 580
    • 群组 1

    星星星星晴

    @星星星星晴

    979
    资料浏览
    580
    帖子
    13
    粉丝
    7
    关注
    注册时间 最后登录

    星星星星晴 取消关注 关注
    副教授

    星星星星晴 发布的最佳帖子

    • RE: 将2.3版本中自带的kOmega模型改名字后单独编译报错?

      我没用过2.3 不过一般情况下自己改模型的标准流程是:

      1. copy
      2. 修改文件内模型命名,用sed -i可以批量修改,比如原来叫kOmega,你现在可以改为MykOmega,注意看看是不是所有的class都改了,是不是有的地方不应该改,然后sed就给你全局 修改了?
      3. 修改Make中的files和option
      files中注意你的目标文件名是否已经修改为心的 .C文件
      sprayFoamTest.C//原来是sprayFoam
      //是否输出到你的USER_LIBBIN,
      //而且有一个心的命名,如果成功编译在你的USER_APPBIN中会出现一个新的文件,以你设置的下面的名字命名,例如我的就是sprayFoamTest
      EXE = $(FOAM_USER_APPBIN)/sprayFoamTest
      
      
      首行添加引用的地方
      LIB_USER_SRC = $(WM_PROJECT_USER_DIR)/src
      
      EXE_INC = \
          -I. \
          -I../reactingParcelFoam \
          -I$(LIB_SRC)/finiteVolume/lnInclude \
          -I${LIB_SRC}/meshTools/lnInclude \
          -I${LIB_SRC}/sampling/lnInclude \
          -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
          -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \
          -I$(LIB_USER_SRC)/lagrangian/basic/lnInclude \
          -I$(LIB_USER_SRC)/lagrangian/intermediate/lnInclude \
          -I$(LIB_USER_SRC)/lagrangian/spray/lnInclude \
          -I$(LIB_USER_SRC)/lagrangian/distributionModels/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
          -I$(LIB_SRC)/transportModels/compressible/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/properties/liquidProperties/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/properties/liquidMixtureProperties/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/properties/solidProperties/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/properties/solidMixtureProperties/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/thermophysicalFunctions/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/SLGThermo/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/chemistryModel/lnInclude \
          -I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \
          -I$(LIB_SRC)/ODE/lnInclude \
          -I$(LIB_SRC)/regionModels/regionModel/lnInclude \
          -I$(LIB_SRC)/regionModels/surfaceFilmModels/lnInclude \
          -I$(LIB_SRC)/combustionModels/lnInclude
      
      EXE_LIBS = \
          -L$(FOAM_USER_LIBBIN) \ //这个地方也别忘了添加,有可能编译的时候找不到改位置 如果不设置的话
          -lturbulenceModels \
          -lcompressibleTurbulenceModels \
          -llagrangian \
          -llagrangianIntermediate \
          -llagrangianTurbulence \
          -llagrangianSpraytest \  //这边要引用你自己修改的模型的名称
          -lspecie \
          -lcompressibleTransportModels \
          -lfluidThermophysicalModels \
          -lliquidProperties \
          -lliquidMixtureProperties \
          -lsolidProperties \
          -lsolidMixtureProperties \
          -lthermophysicalFunctions \
          -lreactionThermophysicalModels \
          -lSLGThermo \
          -lchemistryModel \
          -lradiationModels \
          -lODE \
          -lregionModels \
          -lsurfaceFilmModels \
          -lcombustionModels \
          -lfiniteVolume \
          -lfvOptions \
          -lmeshTools \
          -lsampling
      

      希望能帮到你

      发布在 OpenFOAM
      星
      星星星星晴
    • 【分享+搬运】自定义非均匀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

      发布在 OpenFOAM
      星
      星星星星晴

    星星星星晴 发布的最新帖子

    • RE: 分享一波Chalmers大学OpenFOAM的课程材料

      @cccrrryyy 暴露年龄之
      楼主好人,一生平安

      发布在 OpenFOAM
      星
      星星星星晴
    • 分享一波Chalmers大学OpenFOAM的课程材料

      课程材料原地址是:
      https://www.tfd.chalmers.se/~hani/kurser/OS_CFD/

      27b721a5-8b98-4495-8a79-2f3af0ef94c9-image.png

      个人认为这里面提供的一些材料很有帮助。

      前两天有一次不知道怎么回事,这个网站登不进去了,于是整个爬了下来,可以解压缩后,打开index.html即可。

      6a790b5a-200b-4de6-be94-60c482d0bbe6-image.png

      链接: https://pan.baidu.com/s/1qvtJRMlIj-V4cAcAX6rlaA 提取码: j9gf

      发布在 OpenFOAM
      星
      星星星星晴
    • RE: 植入拉格朗日粒子的Ueqn中需不需要加phicForces?

      @jasper-0 这就得你自己翻code了,得看你用在算什么,dense的话需要考虑体积比什么的了,看起来这玩意就是和dense的有关

      发布在 OpenFOAM
      星
      星星星星晴
    • RE: 植入拉格朗日粒子的Ueqn中需不需要加phicForces?

      你指的SU应该是动量方程里面的源项,就是2-way coupling中的颗粒对流场的影响,那边是包涵了所有的颗粒对流场的影响的,具体你得去code里面查到底是怎么算的。

      如果你考虑的是颗粒受力的话,你就得去submodel里面看了。

      植入颗粒的话,下面这个教程可能说的更清晰
      https://www.foamacademy.com/wp-content/uploads/2018/03/particles_slides.pdf

      至于你说的phicForces 因为没做DPMFoam 所以不知道具体是什么,不过看起来这个好像用在很dense的颗粒模拟里面吧

      发布在 OpenFOAM
      星
      星星星星晴
    • RE: OpenFoam后处理读取颗粒错误

      @dxl 就是foamtovtk以后,在paraview中读取lagrangefield,然后加一个glyphy

      发布在 OpenFOAM
      星
      星星星星晴
    • RE: 高压下interFoam算液滴运动能量不守恒?

      会不会是质量有问题?sum变成了0.5就不对了吧

      发布在 OpenFOAM
      星
      星星星星晴
    • RE: DPMFoam如何得到粒子的位置

      @vbcwl

      https://www.cfd-online.com/Forums/openfoam-programming-development/241034-number-particles-inside-cell.html

      const polyMesh& pmesh = this->owner().mesh();
      List<DynamicList<typename CloudType::parcelType*>>& cellOccupancy =
          this->owner().cellOccupancy();
      
      // check all cells containing particles
      forAll(pmesh.cells(), celli)
      {
          if ( cellOccupancy[celli].size() == 0 )
          {
              // do something
          }
          else
          {
              // do something else
          }
      }
      
      发布在 OpenFOAM
      星
      星星星星晴
    • RE: 关于欧拉-拉格朗日方法

      @小狗狗 DPM中半径不可以为0

      发布在 Algorithm
      星
      星星星星晴
    • RE: 气液固三项流,耦合interFoam和DPMFoam可以实现吗?

      @小狗狗 你这么说的话 就是质点, 不考虑半径

      发布在 OpenFOAM
      星
      星星星星晴
    • RE: 气液固三项流,耦合interFoam和DPMFoam可以实现吗?

      @小狗狗 我没处理过这样的问题,无法给你准确答案。但是我认为是不可以的,就类似你计算drag的时候同样也要考虑半径,当你考虑空间关系的话,应该也需要将其考虑为真正的小球。比如用拉格朗日模拟气泡壁面的话,肯定也要考虑气泡与避免的距离一样。

      发布在 OpenFOAM
      星
      星星星星晴