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. patchInternelField和faceCells对应的值竟然可能不一致

patchInternelField和faceCells对应的值竟然可能不一致

Scheduled Pinned Locked Moved OpenFOAM
3 Posts 2 Posters 2.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.
  • S Offline
    S Offline
    Samuel-Tu
    wrote on last edited by Samuel-Tu
    #1

    今天发现一件奇怪的事,我用patchInternalField找出nut的一个边界紧邻的内部网格上的值,和用faceCells先找到这个紧邻网格的编号,再再nut场中取出这一个值,发现大多数网格两个值是一样的,但是存在一个网格值不一样,代码如下:

    scalarField nutW1=nut.boundaryField()[patchID].patchInternalField();
    
    labelList cellID=mesh.boundary()[patchID].faceCells();
    scalarField nutW2;
    forAll(cellID,i)
    {
        nutW2[i]=nut[cellID[i]];
    }
    

    最后发现nutW1和nutW2几乎所有值一样,但是第一个值,nutW1[0]=0,nutW2[0]是一个正常值。我们知道nut是湍流粘度,一定是大于0的。。不知道为什么会这样。好奇怪。。

    X 1 Reply Last reply
  • X Online
    X Online
    xpqiu 超神
    replied to Samuel-Tu on last edited by
    #2

    @Samuel-Tu 在 patchInternelField和faceCells对应的值竟然可能不一致 中说:

    今天发现一件奇怪的事,我用patchInternalField找出nut的一个边界紧邻的内部网格上的值,和用faceCells先找到这个紧邻网格的编号,再再nut场中取出这一个值,发现大多数网格两个值是一样的,但是存在一个网格值不一样,代码如下:

    scalarField nutW1=nut.boundaryField()[patchID].patchInternalField();
    
    labelList cellID=mesh.boundary()[patchID].faceCells();
    scalarField nutW2;
    forAll(cellID,i)
    {
        nutW2[i]=nut[cellID[i]];
    }
    

    最后发现nutW1和nutW2几乎所有值一样,但是第一个值,nutW1[0]=0,nutW2[0]是一个正常值。我们知道nut是湍流粘度,一定是大于0的。。不知道为什么会这样。好奇怪。。

    试了一下,没法复现这个情况,我这边用这样方式算出来 nutw1和nutw2 是一样的。
    PS: 你代码里面先定义 scalarField nutw2, 然后直接用 nutw2[i] 这样来赋值,会导致段错误的吧,nutw2 没有分配空间的啊,可能你实际用的时候不是这样写的?

    S 1 Reply Last reply
  • S Offline
    S Offline
    Samuel-Tu
    replied to xpqiu on last edited by
    #3

    @xpqiu 对,nutW2是分配了空间了的,大小就是cellID的大小。
    造成这个的原因我怀疑是不是因为的我边界有转折导致的。。我是提取的roof这个边界上对应的第一层网格的值的。
    121.png
    因为当时我做计算的时候除以了nutW1之后计算报错了,我才发现nutW1[0]=0,不能被除。然后才改成使用了nutW2的

    1 Reply Last reply

  • Login

  • Don't have an account? Register

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