Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 请教rhopimpleFoam的一个问题

请教rhopimpleFoam的一个问题

已定时 已固定 已锁定 已移动 OpenFOAM
2 帖子 1 发布者 4.0k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • Y 离线
    Y 离线
    yuan_neu
    写于 最后由 CFD中文网 编辑
    #1
    15 if (pimple.transonic())
       16 {
       17     surfaceScalarField phid
       18     (
       19         "phid",
       20         fvc::interpolate(psi)
       21        *(
       22             fvc::flux(HbyA)
       23           + rhorAUf*fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho)
       24         )
       25     );
       26 
       27     MRF.makeRelative(fvc::interpolate(psi), phid);
       28 
       29     while (pimple.correctNonOrthogonal())
       30     {
       31         fvScalarMatrix pEqn
       32         (
       33             fvm::ddt(psi, p)
       34           + fvm::div(phid, p)
       35           - fvm::laplacian(rhorAUf, p)
       36          ==
       37             fvOptions(psi, p, rho.name())
       38         );
       39 
       40         pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
       41 
       42         if (pimple.finalNonOrthogonalIter())
       43         {
       44             **phi == pEqn.flux();//??**
       45         }
       46     }
       47 }
    

    代码中,为什么对phi的更新是
    phi == pEqn.flux();//?????????
    我查阅了rhosimpleFoam,发现也是一样,实在是弄不明白

    1 条回复 最后回复
  • Y 离线
    Y 离线
    yuan_neu
    写于 最后由 编辑
    #2

    http://blog.sina.com.cn/s/blog_97cadb820102w2h0.html
    在这里找到一些回答

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]