Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. Tens
    T
    • Profile
    • Following 1
    • Followers 4
    • Topics 20
    • Posts 158
    • Groups 1

    Tens

    @Tens

    讲师

    34
    Profile views
    158
    Posts
    4
    Followers
    1
    Following
    Joined Last Online

    Tens Unfollow Follow
    讲师

    Latest posts made by Tens

    • RE: MPPICInterFoam 中如何通过alphac得到alpha1和alpha2?

      :chitang:

      posted in OpenFOAM
      T
      Tens
    • RE: MPPICInterFoam 中如何通过alphac得到alpha1和alpha2?

      @zhe 对于欧拉场来讲,是单纯的两相流,因为考虑了颗粒体积,所以连续相的体积只有alphac*Vcell,即Vc,所以连续相的计算是在这个基础上进行的,这也是为什么控制方程中都要乘alphac,可能结合单纯的MPPICFoam和InterFoam更好理解,只是将InterFoam的处理运用到MPPICFoam的连续相中

      posted in OpenFOAM
      T
      Tens
    • RE: MPPICInterFoam 中如何通过alphac得到alpha1和alpha2?

      @zhe 体积和占比并不是一个概念

      posted in OpenFOAM
      T
      Tens
    • RE: MPPICInterFoam 中如何通过alphac得到alpha1和alpha2?

      @zhe 对alpha1和alpha2来说,分母不是Vcell,而是Vc
      d39cbc61-6364-43e0-a9d5-a6080d1b453d-image.png

      posted in OpenFOAM
      T
      Tens
    • RE: MPPICInterFoam 中如何通过alphac得到alpha1和alpha2?

      @zhe alpha1一直都不是占cell的体积分数,只是占连续相的体积分数,不是alpha1+alpha2=alphac,而是再次将alphac看做单位1,alpha1+alpha2=1
      对于你的例子,只有粒子和水,那就是alphac=0.8,alpha1=1
      2160b90a-c5db-40fe-ad9b-9253977bcce1-image.png

      posted in OpenFOAM
      T
      Tens
    • RE: MPPICInterFoam中静置情况下水流沸腾

      颗粒和水的密度比大吗,是不是有可能是颗粒在浮力作用下运动引起的波动呢

      posted in OpenFOAM
      T
      Tens
    • RE: dsmcFoam输出哪些流场文件是由什么控制的

      @雨罢_清宵 虽然不知道这个文件是干什么的,但是从名字上可以看出来是和氩气Ar有关的,那其他案例中是否组分一样呢,组分不同当然输出文件会不同了

      posted in OpenFOAM
      T
      Tens
    • RE: dsmcFoam输出哪些流场文件是由什么控制的

      @雨罢_清宵 输出的流场信息不是由fvSolution来设置的,而是在创建这个变量的时候设置的,比如

      volScalarField rho
          (
              IOobject
              (
                  "rho",
                  runTime.timeName(),
                  mesh,
                  IOobject::NO_READ,
                  IOobject::AUTO_WRITE   //输出,NO_WRITE则不输出
              ),
              thermo.rho()
          );
      
      posted in OpenFOAM
      T
      Tens
    • RE: MPPICInterFoam 中如何通过alphac得到alpha1和alpha2?

      @jasper-0
      互相交流,以前从没这么认真的看过这么底层的东西,这次也学到了很多:chitang:

      posted in OpenFOAM
      T
      Tens
    • RE: MPPICInterFoam 中如何通过alphac得到alpha1和alpha2?

      @jasper-0
      0e85dbdf-ee3c-449b-ac56-b0e366bf4679-image.png
      懂了,现在就是说是红色部分的问题对吧,我看了代码中这部分是通过下面代码求的

      fvc::surfaceIntegrate(psiIf, phiPsi)
      

      而里面的phiPsi则对应的这一部分
      87b34b23-8ec8-40f0-bb9e-001511b06be6-image.png
      通过追踪

      alphaPhi10 = talphaPhi1Un;
      
      //又查到
      
      tmp<surfaceScalarField> talphaPhi1Un
              (
                  fvc::flux
                  (
                      phiCN(),
                      cnCoeff*alpha1 + (1.0 - cnCoeff)*alpha1.oldTime(),
                      alphaScheme
                  )
                + fvc::flux
                  (
                     -fvc::flux(-phir, alpha2, alpharScheme),
                      alpha1,
                      alpharScheme
                  )
              );
      

      而其中的phiCN则

          tmp<surfaceScalarField> phiCN(alphaPhic);
      
          // Calculate the Crank-Nicolson off-centred volumetric flux
          if (ocCoeff > 0)
          {
              phiCN = cnCoeff*alphaPhic + (1.0 - cnCoeff)*alphaPhic.oldTime();
          }
      

      如果不采用CN格式phiCN就等于alphaPhic,又查到

       alphacf = fvc::interpolate(alphac);
      alphaPhic = alphacf*phi;
      

      所以本质上图一中红色部分还是通过alphac*alpha1得到的

      posted in OpenFOAM
      T
      Tens