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. 关于passive scalar受湍流作用的疑问?

关于passive scalar受湍流作用的疑问?

已定时 已固定 已锁定 已移动 OpenFOAM
4 帖子 3 发布者 4.7k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • chpjz0391C 离线
    chpjz0391C 离线
    chpjz0391
    写于 最后由 编辑
    #1

    最近在把passive scalar项写进simpleFoam来算一个标量场。但是对于湍流对这个标量场的影响有一点疑问。

    首先是标量场的方程,参考scalarTransportFoam方程如下:

                fvm::ddt(T)
              + fvm::div(phi, T)
              - fvm::laplacian(DTt, T)
             ==
                fvOptions(T)
    

    然后在creatFields.H 里面定义湍流对其影响。公式参考如下
    D = D_{molecular} + \frac{\nu_t}{Sc_t}

    代码如下:

    Info<< "Reading diffusivity DT\n" << endl;
    
    dimensionedScalar DT
    (
        transportProperties.lookup("DT")
    );
    
    Info<< "Reading turbulent Schmidt number\n" << endl;
    
    dimensionedScalar Sct
    (
        transportProperties.lookup("Sct")
    );
    
    volScalarField DTt
    (
        IOobject
        (
            "DTt",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        nut/Sct + DT,
        nut.boundaryField().types()
    );
    

    这里面Sct是Schmidt number, DT是 标量场的diffusivity
    这也是目前我正在用的solver
    但是网络还有另一种改法,不需要知道标量场的diffusivity,取而代之的是需要知道Prt 和Pr
    标量方程改成如下:

    volScalarField kappaEff
    (
        "kappaEff",
        turbulence->nu()/Pr + turbulence->nut()/Prt
    );
    
    fvScalarMatrix DTEqn
    (
        //fvm::ddt(DT)
      fvm::div(phi, DT)
      - fvm::laplacian(kappaEff, DT)
      ==
      fvOptions(TS) 
    );
    

    DTEqn.relax();
    DTEqn.solve();

    请问这两种有什么区别?为什么一个需要标量场的diffusivity 另一个不需要?

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

    我懂了。 这里的diffusivity= nu()/Pr 对吧

    D J 2 条回复 最后回复
  • D 离线
    D 离线
    Dingcy
    在 中回复了 chpjz0391 最后由 编辑
    #3

    @chpjz0391 什么方向?溶质传输?还是RTD?

    我希望我得到的少一点,少一点,再少一点......

    1 条回复 最后回复
  • J 离线
    J 离线
    Joe_Bin
    在 中回复了 chpjz0391 最后由 编辑
    #4

    @chpjz0391 嗯 是的 。 其实是一样的 换汤不换药

    1 条回复 最后回复

  • 登录

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