Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. upc_ngh
    U
    • Profile
    • Following 5
    • Followers 3
    • Topics 19
    • Posts 114
    • Groups 0

    upc_ngh

    @upc_ngh

    123
    Profile views
    114
    Posts
    3
    Followers
    5
    Following
    Joined Last Online

    upc_ngh Unfollow Follow

    Latest posts made by upc_ngh

    • RE: 关于网格高度变化的动网格设置

      @田畔的风 非常感谢两位老师@李东岳 ,感觉久旱逢甘露 ,话不多说,我要开始汲取营养了,后面有相关问题和进展,再向老师们请教:huahua: :huahua: :huahua:

      posted in OpenFOAM
      U
      upc_ngh
    • RE: 关于网格高度变化的动网格设置

      @李东岳 东岳老师,您好,目前我遇到了和上面那个同学类似问题@xiaolin ,具体问题如下:二维内,气固两相流动过程中,固态颗粒和壁面(半圆形壁面)碰撞,粘连形成沉积。目前的想法是当贴近壁面的第一层网格内,固相体积分数大于某值时,壁面的边界条件就向上移动一个单元,就是上面"alpha与point的关系“。目前发现相关动网格编程资料很少,盼望各位老师能给点建议或者参考资料学习一下,提前谢谢各位大神!
      图片1.jpg

      posted in OpenFOAM
      U
      upc_ngh
    • RE: 气液流动过程中温度直接求解问题

      @李东岳 收到,谢谢东岳老师!

      posted in OpenFOAM
      U
      upc_ngh
    • RE: 气液流动过程中温度直接求解问题

      @李东岳 对的李老师,我去年参加过OKS课程,周四晚上有时间的,那就提前谢谢李老师了!:high:

      posted in OpenFOAM
      U
      upc_ngh
    • RE: 气液流动过程中温度直接求解问题

      @李东岳 东岳老师,上次那个程序您帮忙看看,有结果了吗:chouchou: :chouchou: :chouchou:

      posted in OpenFOAM
      U
      upc_ngh
    • RE: 气液流动过程中温度直接求解问题

      @李东岳 哇,谢谢东岳老师,非常感谢感谢:qinqin: :qinqin:

      posted in OpenFOAM
      U
      upc_ngh
    • RE: 气液流动过程中温度直接求解问题

      @李东岳 东岳老师,上次温度方程求解错误主要是由于两个原因,一个是边界条件给错了,还一个是温度方程和动量、连续相方程耦合求解了,经过更改已经可以正常运行。但目前在验证算例时发现和实验值存在较大的误差。
      验证算例如下: 一根长为 12 m, 直径为 1 m的垂直管, 管道上端为入口, 管道下端为出口,出口压力是 1e5Pa, 气液相分别为空气和水,温度均为 25 ℃ .入口处 x = 0, 出口处 x = 12 m.初始状态如图 1(a)所示, 气相相含率 0. 2, 气相速度 Vg=0;液相相含率 0. 8, 液相速度 Vl=10 m/s。使用双欧拉求解。耦合求解动量及连续性方程(参考瞬态可压缩压力泊松方程),收敛后求解相率方程,进行下一时间步。图一是示意图,图二是解析解的气相分数及液相流速,图三为求解出的压力、持气率、液相流速,可以看出波动很大,尤其持气率波动尤其大,已经排除网格密度的原因。能力有限,查找一周的错误了,实在找不到什么别的原因了,还望东岳老师和前辈指点一二。
      RJGAUFG(HTID`(89F@THYEN.png

      RIMYF9`R5I}FSN_L@QJ_0BD.png

      }G682D}15VJO01DI1)BC7AY.png
      气液相动量方程:

      fvVectorMatrix UEqn1
                 (
                      fvm::ddt(alpharho1, U1) + fvm::div(alphaphi1, U1)
                      == -pGrad1-Sg
                 );
                 UEqn1.solve();
      
                 fvVectorMatrix UEqn2
                 (
                      fvm::ddt(alpharho2, U2) + fvm::div(alphaphi2, U2)
                         == -pGrad2-Sl
                  );
                 UEqn2.solve();
      

      连续性方程

      psi1 = 1/Cg/Cg;
                 psi2 = rho2/p;
      
                 volScalarField rAU1(1.0/UEqn1.A());
                 surfaceScalarField rhorAUf1("rhorAUf1", fvc::interpolate(alpharho1*rAU1/rhog));
                 volVectorField HbyA1(constrainHbyA(rAU1*UEqn1.H(), U1, p));
      
                 surfaceScalarField phiHbyA1
                   (
                          "phiHbyA1",
                       (
                         fvc::interpolate(psi1*alpha1/rhog)
                         *fvc::flux(HbyA1)
                       ) 
                    );
      
                 volScalarField alphapsi1("alphapsi1",alpha1*psi1/rhog);
      
         
                 volScalarField rAU2(1.0/UEqn2.A());
      
                 surfaceScalarField rhorAUf2("rhorAUf2", fvc::interpolate(alpharho2*rAU2/rho2));
      
                 volVectorField HbyA2(constrainHbyA(rAU2*UEqn2.H(), U2, p));
      
                 surfaceScalarField phiHbyA2
                   (
                          "phiHbyA2",
                        (
                            fvc::interpolate(psi2*alpha2/rho2)
                            *fvc::flux(HbyA2)
                        ) 
                    );
      
                 volScalarField alphapsi2("alphapsi2",alpha2*psi2/rho2);
      
        
         
      fvScalarMatrix pEqn
              (
                  fvm::ddt(alphapsi1, p)
                + fvm::ddt(alphapsi2, p)
                + fvm::div(phiHbyA1, p)
                + fvm::div(phiHbyA2, p)
                - fvm::laplacian(rhorAUf1, p)
                - fvm::laplacian(rhorAUf2, p)
              );
        pEqn.solve();
        rho1 = p*psi1;
        rho1.correctBoundaryConditions();
      
      surfaceScalarField phi1 = fvc::interpolate(rho1*U1) & mesh.Sf();
         surfaceScalarField phi2 = fvc::interpolate(rho2*U2) & mesh.Sf();  
      
         fvScalarMatrix alpharhoEqn1
          (
              fvm::ddt(rho1,alpha1)
            + fvm::div(phi1,alpha1)
          );
         
      fvScalarMatrix alpharhoEqn2
          (
              fvm::ddt(rho2,alpha2)
            + fvm::div(phi2,alpha2)
          );
      
       
         alpharhoEqn1.solve();
         alpharhoEqn2.solve();
        
        
         alpha1 = alpha1/(alpha1+alpha2);
         alpha2 = 1 - alpha1;
         alpha2.correctBoundaryConditions();
      
      posted in OpenFOAM
      U
      upc_ngh
    • RE: 气液流动过程中温度直接求解问题

      @李东岳 对的,东岳老师,我是先求解动量方程,压力修正方程之后,最后一步再求解温度场,压力这个时候已经是已知的吧,然后就各种发散

      posted in OpenFOAM
      U
      upc_ngh
    • 气液流动过程中温度直接求解问题

      目前在求解两相可压缩流动过程时的温度方程时出现一些问题,请教一下各位老师前辈。求解发现,管道内流体流到处温度急剧下降至0.01K量级,直至发散。试图将不同相除以其密度再相加,问题依旧存在,感觉方程和植入没啥问题,不知道为啥始终求得温度接近0K,烦请各位老师指点一二,谢谢!
      1.png

      alphacp1 = alpha1*cp1*rho1;
      alphacp2 = alpha2*cp2*rho2;
      alphaphicp1 = fvc::interpolate(alphacp1*U1) & mesh.Sf();
      alphaphicp2 = fvc::interpolate(alphacp2*U2) & mesh.Sf();
      
      surfaceScalarField birho1 = fvc::interpolate(alpha1*rho1*U1) & mesh.Sf()/fvc::interpolate(rho1);
      
      surfaceScalarField birho2 = fvc::interpolate(alpha2*rho2*U2) & mesh.Sf()/fvc::interpolate(rho2);
      
       fvScalarMatrix TEqn
          (
              fvm::ddt(alphacp1,T)
            + fvm::ddt(alphacp2,T)
            + fvm::div(alphaphicp1,T)
            + fvm::div(alphaphicp2,T)
            + fvc::div(birho1,p)
            + fvc::div(birho2,p)
            - fvc::ddt(alpha1,p)
            - fvc::ddt(alpha2,p)
          );
           TEqn.solve();
      
      
      posted in OpenFOAM
      U
      upc_ngh
    • RE: fluent求助帖 变径管(PBM模型)的收敛问题

      @李东岳 感谢东岳老师,在去掉PBM模型(其他条件不变时)后,收敛起来很容易,证明应该不是流场的原因。加上PBM后不收敛的原因是什么,还请您指点

      posted in Fluent
      U
      upc_ngh