Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. Jacobian
    J
    • Profile
    • Following 0
    • Followers 2
    • Topics 6
    • Posts 34
    • Groups 0

    Jacobian

    @Jacobian

    1083
    Profile views
    34
    Posts
    2
    Followers
    0
    Following
    Joined Last Online

    Jacobian Unfollow Follow

    Best posts made by Jacobian

    • RE: paraview并行出错
      1. 终端输入mpirun -np n pvserver (n根据情况调整),保持该终端不关闭
      2. 打开paraFoam,file -> connect -> add server : 其中,Name随意, Server type 用 client/server,Host 用 localhost,port 默认即可。
      3. 加载 .foam文件

      感觉这种方法要比reconstructPar快

      posted in OpenFOAM
      J
      Jacobian

    Latest posts made by Jacobian

    • RE: 粒子与网格归属问题
      把粒子按x坐标排序,i是序号,形成xorder(i,1)存x坐标,xorder(i,2)存粒子编号;
      同理,把粒子按y、z坐标排序,形成yorder(i,1..2),zorder(i,1..2)
      计算网格顶点的xmin,xmax,ymin,ymax,zmin,zmax
      
      loop:循环网格{
            
             找到xmin,xmax在xorder(i,1)中的位置
             同理,处理ymin,ymax,zmin,zmax
             new 一个数组 possible(i)
           
             loop: 处理上面缩小的范围{
                    把以上范围内的粒子编号xorder(i,2),yorder(i,2),
                    zorder(i,2)插入排 序到possible(i)
             }
      
             loop: 确定真正在网格中的粒子{
                   possible(i)中有连续三个相同编号,判断是否真的在网格内,
                   如果真在,就存起来
             }
      }
      

      至于每个网格里的粒子编号怎么存,搞一个数组叫data(i)一个接一个地存上面伪代码得到的编号,另一个数组叫position(j)存单元j在data(i)中的起始位置
      不知道问题理解对不对,抛个砖试试看:chouchou:

      posted in OpenFOAM
      J
      Jacobian
    • RE: rhoSimpleFoam计算翼型发散

      @dzw05 学习了,div(phid,p) 确实不用 bounded 为好
      @东岳 把 z01~02 后面的zip去掉,然后对着 rhoSimple.zip 选择解压

      posted in OpenFOAM
      J
      Jacobian
    • RE: rhoSimpleFoam计算翼型发散

      非气动专业,尝试和大佬讨论,想问一下:
      (1)div(phid,p) 的格式为什么不用 bounded Gauss upwind 呢?经验来看,bounded格式稳定很多;
      (2)grad用最小二乘,以及interpolation用迎风也是为了数值稳定性吧?但貌似当前算例都用中心差分也可以,所以为了精度给改了。
      按照以上两项修改试着跑了几万步,不过不是特别收敛,连续性方程不满足,升阻力系数还在震荡。是不是这个物理问题本身是非定常的,或者是网格太粗糙呢?
      还有就是挺好奇这个网格使用sHM画的吗?像这类混合网格如何处理trailing edge这种特变尖锐的地方?

      posted in OpenFOAM
      J
      Jacobian
    • RE: OpenFOAM中的binaryblock如何读取?

      可能你的操作系统位数和别人不一样。
      试下装一个和你目前系统位数不一样的虚拟机(或者用超算、别人的电脑)
      然后把case搬到这台机器上,在controlDict 里面writeFormat 改成 ascii,终端键入foamFormatConvert看看能不能转换

      posted in OpenFOAM
      J
      Jacobian
    • RE: 请教密度随深度连续变化设置问题

      抱歉说漏了,rho、rho.water、rho.air这三个文件都要放进0文件夹。funkySetFields是能够使用函数表达式的,用来设置场内分布;边界上的值可以用funkySetBoundaryField来设置。在system下添加funkySetFieldsDict和funkySetBoundaryDict来使用。里面的内容,随便举个例子:

      //funkySetFields
      expressions
       (
          pressureWater
          {
           field p_rgh; //需要设置的物理量
           expression "100000 - 99800*pos().y";  //这行写你要的函数
           condition  "(pos().x <= 5) && (pos().x >= 0)"; //起作用的范围、可以注释掉,不影响使用
           keepPatches 1; //填1就好,要不然可能把边界条件改成zeroGradient
          }
       );
      
      //funkySetBoundaryField
      pressureWater
       {
          field p_rgh;
          expressions
          (
              {
                  target value;
                  patchName INLET; //选择边界面
                  //variables "maxY=max(pts().y);";  //可以注释掉,不影响使用
                  expression "100000 - 99800*pos().y"; //这行写你要的函数
              }
               //可以再expression下同时写多个边界的设置
              {
                  target value;
                  patchName OUTLET;
                  //variables "maxY=max(pts().y);";
                  expression "100000 - 99800*pos().y";
              }
          );
      }
      

      更多细节参考这里链接文本,以及swak4foam的tutorials
      另外,你的这个算例里面压力也应当要设置,需要和密度相匹配。

      posted in OpenFOAM
      J
      Jacobian
    • RE: paraview并行出错
      1. 终端输入mpirun -np n pvserver (n根据情况调整),保持该终端不关闭
      2. 打开paraFoam,file -> connect -> add server : 其中,Name随意, Server type 用 client/server,Host 用 localhost,port 默认即可。
      3. 加载 .foam文件

      感觉这种方法要比reconstructPar快

      posted in OpenFOAM
      J
      Jacobian
    • RE: 请教密度随深度连续变化设置问题

      compressibleInterFoam读入rho的方式是READ_IF_PRESENT,也就是说存在rho文件的话,就会被读进去。所以可以在0文件夹下自己添加一个rho,然后用funkySetField 设置你要的初始场。
      至于“如何设置边界密度随空间变化”,compressibleInterFoam是基于压力的求解器,也就是密度由压力决定。所以只要把p_rgh设置对了就可以吧?

      posted in OpenFOAM
      J
      Jacobian
    • RE: OvermultphaseInterDyMFoam

      @maoyanjun_dut 非常感谢你的工作!另外想再请教个问题,如果要把overset用到compressibleInterFoam这类可压缩求解器里面的话,可压缩带来的dgdt,pEqnComp应该怎么处理?

      posted in OpenFOAM
      J
      Jacobian
    • RE: 多相流求解器的Ueqn中 turbulence->divDevRhoReff 在哪段代码求了平均粘性系数?

      谢谢东岳老师考古!讲解非常清楚!

      posted in OpenFOAM
      J
      Jacobian
    • RE: Floating point exception

      如果是除0的话,应该报 #2 Foam::divide at ??:? 这样的错误。
      你的这个错误像是缺少了自己写的库文件。FOAM_USER_APPBIN下能找到吗?试着重新编译一下?

      posted in OpenFOAM
      J
      Jacobian