CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    如何导出sprayFoam能量方程中的各项参数?

    OpenFOAM
    3
    5
    447
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 香柏树
      香柏树 最后由 编辑

      目前想分析下能量方程中各项参数(如非定常项、对流项、亚格子项,化学反应源项和扩散项)在总能量中所占比例,针对下面sprayFoam中求解能量方程的EEqn.H文件,请教大家采用什么方法能导出各项参数到新时刻文件夹中呢?
      EEqn.H

      {
          volScalarField& he = thermo.he();
      
          fvScalarMatrix EEqn
          (
              fvm::ddt(rho, he) + mvConvection->fvmDiv(phi, he)
            + fvc::ddt(rho, K) + fvc::div(phi, K)
            + (
                  he.name() == "e"
                ? fvc::div
                  (
                      fvc::absolute(phi/fvc::interpolate(rho), U),
                      p,
                      "div(phiv,p)"
                  )
                : -dpdt
              )
            - fvm::laplacian(turbulence->alphaEff(), he)
           ==
              rho*(U&g)
            + parcels.Sh(he)
            + radiation->Sh(thermo, he)
            + Qdot
            + fvOptions(rho, he)
          );
      
          EEqn.relax();
      
          fvOptions.constrain(EEqn);
      
          EEqn.solve();
      
          fvOptions.correct(he);
      
          thermo.correct();
          radiation->correct();
      
          Info<< "T gas min/max   " << min(T).value() << ", "
              << max(T).value() << endl;
      }
      
          drhoedt=fvm::ddt(rho, he);
      
      
      李东岳 1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 @香柏树 最后由 李东岳 编辑

        @香柏树 你可以把这些项声明一些场出来,然后显性离散计算一下。比如时间项:

        volScalarField ddt
            (
                IOobject
                (
                    "rho",
                    runTime.timeName(),
                    mesh,
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                thermo.rho()
            );
        
        计算之后
        
        ddt = fvc::ddt(rho, he);
        

        CFD高性能服务器 http://dyfluid.com/servers.html

        香柏树 1 条回复 最后回复 回复 引用
        • 香柏树
          香柏树 @李东岳 最后由 编辑

          @李东岳 在 如何导出sprayFoam能量方程中的各项参数? 中说:

          ddt = fvc::ddt(rho, he);

          感谢李老师回复。我按照您说的先声明了时间项后进行计算后,显性离散计算“ddt = fvc::ddt(rho, he);”这块是怎么实现的呢?是在sprayFoam.C文件中加入这句么?还是在后处理中进行啊?初学openfoam有点不太明白

          Z 1 条回复 最后回复 回复 引用
          • Z
            zhouxu @香柏树 最后由 编辑

            @香柏树 在能量方程计算之后,或者当前时间步计算结束之前

            香柏树 1 条回复 最后回复 回复 引用
            • 香柏树
              香柏树 @zhouxu 最后由 编辑

              @zhouxu 谢谢,我尝试下。

              1 条回复 最后回复 回复 引用
              • First post
                Last post