Skip to content
  • 最新
  • Categories
  • 东岳流体
  • 随机看[请狂点我]
Skins
  • 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

  • Default (No Skin)
  • No Skin
Collapse
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 方程离散!fvm:div到底output出什么?

方程离散!fvm:div到底output出什么?

Scheduled Pinned Locked Moved OpenFOAM
4 Posts 2 Posters 5.9k Views
  • Oldest to Newest
  • Newest to Oldest
  • Most Votes
Reply
  • Reply as topic
Log in to reply
This topic has been deleted. Only users with topic management privileges can see it.
  • siboS Offline
    siboS Offline
    sibo
    wrote on last edited by 李东岳
    #1

    大家好,
    想请教一个基础问题。
    根据这篇post,http://dyfluid.com/divsnGrad.html,
    fvm:div(Phi, U)应该生成的是对一个cell所有面的流率与速度乘积的加和。按照Guass theorem推到下来的确是这样。
    但是那么对应的ddt(rho,U)就应该是rho, U, V的乘积是么?因为如果做了对Phi, U 项做了加和且没有除以体积的话,那其他项都应该有对体积的积分啊。如果加入source term,也应该对source term进行体积积分。请问大家是这样么?

    因为之前一个师兄写的solver里面出现:rho(t+1)=rho(t)+dt*K*fvc::div(rhoFlux) 类似的code,这样来看的话,fvc::div(rhoFlux) 应该生成的是一个对某cell所有面的rhoFlux的加和,但是rho项并没有对cell的体积进行任何操作啊。。。所以这个dimension在等式左右两边就不一样啊。。

    所以我在想div的这个操作符,是不是还加入了除以体积的操作?
    请各位大神指教!
    谢谢啦!

    1 Reply Last reply
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    wrote on last edited by
    #2

    http://www.cfd-china.com/topic/703/fvc-div的前后量纲关系到底是啥/9

    所以我在想div的这个操作符,是不是还加入了除以体积的操作?

    上面的帖子你看看。

    9月CFD算法编程课: http://dyfluid.com/class.html

    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 Reply Last reply
  • siboS Offline
    siboS Offline
    sibo
    wrote on last edited by
    #3

    @李东岳 Exactly what I wanted! Thanks a lot!
    不过还是想请前辈check一下:
    也就是说,
    fvm返回的是fvMatrix, fvm::ddt, fvm:laplacian都会将结果乘以体积。由于fvm::div用gauss theorem做的面的flux加和,所以没有乘以体积;
    fvc返回的GeometricField, 所以在对fvm和fvc进行加减时会重载operator。
    例如,fvc::div 返回的dimension 本来是 fvm::div返回的dimension 再除以了体积V的,但是在一个方程中同时出现fvm和fvc时,会对 fvc::div 项也乘以体积以保证量纲守恒。但是单独进行fvc::div操作时,返回的dimension 仍然再除以了体积V的。
    请问我理解的对么?谢谢!

    1 Reply Last reply
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    wrote on last edited by
    #4

    是的你说的很正确。

    9月CFD算法编程课: http://dyfluid.com/class.html

    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 Reply Last reply

  • Login

  • Don't have an account? Register

  • Login or register to search.
  • First post
    Last post
0
  • 最新
  • Categories
  • 东岳流体
  • 随机看[请狂点我]