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. LES直流槽道边界层模拟,如何得到正则化速度u+以及正则化坐标y+?

LES直流槽道边界层模拟,如何得到正则化速度u+以及正则化坐标y+?

已定时 已固定 已锁定 已移动 OpenFOAM
大涡模拟
125 帖子 14 发布者 116.4k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    写于 最后由 李东岳 编辑
    #92

    我之前一直想做dnsfoam解析,一直没倒出时间。dnsfoam这个跟icofoam啥的唯一的区别就是这个力,这个力的生成算法我还没细看,dnsfoam解析也一直没写出来。我记得是有什么傅里叶变化啥的东西在里面。

    另外新版dnsfoam应该已经没了。变成了一个力源项。

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

    S 1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    在 中回复了 SSSSK 最后由 编辑
    #93

    @SSSSK 确实没看过对均匀各向同性湍流加force的,一般都是多是对类似channel395这种非均匀湍流加force。就是对均匀各向同性湍流加force,要想想是否理论上加force是否就能实现不衰减?或者加的force需要一个很复杂的形式才能实现,而不像channel395只加一个固定的压力梯度。

    S 1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    写于 最后由 编辑
    #94

    对均匀各向同性湍流加force的自保持,如果只是要实现到目标结构位置湍流特性是所想要的,可以对入口湍流乘以放大系数。均匀各向同性湍流衰减曲线的对数率形式接近一条直线,根据衰减曲线估计下在入口处对湍流增强洗漱,应该能调整出目标位置的所需要的湍流特性

    1 条回复 最后回复
  • S 离线
    S 离线
    SSSSK
    在 中回复了 李东岳 最后由 编辑
    #95

    @李东岳 新版的能像dnsfoam一样加白噪声力吗?原版的我稍微看了一下,是通过高斯分布得到白噪声然后再根据维纳过程施加和时间有关的力,最后傅里叶变换变成方程里的force,李老师你说的新版是of11和12吗?

    1 条回复 最后回复
  • S 离线
    S 离线
    SSSSK
    在 中回复了 coolhhh 最后由 编辑
    #96

    @coolhhh 对这种湍流加的力和channel不一样,channel只需要一个固定压力梯度就行,但这个是从波数域加的,需要傅里叶变化,但傅里叶变换里面有个kmesh,就是这个不能并行

    1 条回复 最后回复
  • L 离线
    L 离线
    LESKIN
    在 中回复了 2019201300 最后由 编辑
    #97

    @2019201300文件过期呜呜呜,最近也在研究管道流动问题,想知道大佬怎么做的,可以重新传一下文件嘛

    2 1 条回复 最后回复
  • 2 离线
    2 离线
    2019201300
    在 中回复了 LESKIN 最后由 编辑
    #98

    @LESKIN 这个是个槽道,流向和展向都是周期性边界条件。通过网盘分享的文件:Channel180_Kong.rar
    链接: https://pan.baidu.com/s/1TvKiFbKk-H0n-pm5LBmXyA?pwd=vnfi 提取码: vnfi

    L 1 条回复 最后回复
  • L 离线
    L 离线
    LESKIN
    在 中回复了 2019201300 最后由 编辑
    #99

    @2019201300 在 LES直流槽道边界层模拟,如何得到正则化速度u+以及正则化坐标y+? 中说:

    @LESKIN 这个是个槽道,流向和展向都是周期性边界条件。通过网盘分享的文件:Channel180_Kong.rar
    链接: https://pan.baidu.com/s/1TvKiFbKk-H0n-pm5LBmXyA?pwd=vnfi 提取码: vnfi

    太感谢了

    1 条回复 最后回复
  • 东方白杨东 离线
    东方白杨东 离线
    东方白杨
    写于 最后由 编辑
    #100

    您好,我也是用的上面说的论文里preturb生成初始湍流场方法,可是目前有点问题。一开始初始场的速度型线确实是抛物线,但是算50s就开始变形,算到1000秒速度型线就非常离谱了
    下图是50s的
    f7ab4305d5ac9913f821a7bc952cc13.png
    下图是1000s的

    a69d4a6b2a2d9515062a8bfcf466723.png

    CAS-IET

    1 条回复 最后回复
  • 东方白杨东 离线
    东方白杨东 离线
    东方白杨
    写于 最后由 编辑
    #101

    是我搞错了。。。湍流就应该这样,抛物线型速度型线是层流吧

    CAS-IET

    1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    写于 最后由 编辑
    #102

    是的,你得湍流的没问题,就是时间平均多搞搞,目前不对称

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

    1 条回复 最后回复
  • 东方白杨东 离线
    东方白杨东 离线
    东方白杨
    写于 最后由 编辑
    #103

    @李东岳 嗯嗯,老师您知道Eugene de Villiers论文里的LES数据在哪里能下载吗?我想对比一下

    CAS-IET

    1 条回复 最后回复
  • 东方白杨东 离线
    东方白杨东 离线
    东方白杨
    写于 最后由 编辑
    #104

    一个星期,终于把这个算例搞明白了,纪念一下,感谢上面前辈们的分享!
    image.png

    CAS-IET

    李东岳李 1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    在 中回复了 东方白杨 最后由 李东岳 编辑
    #105

    @东方白杨 我去,厉害啊,老铁,你这个厉害。

    上面帖子当时我没回,后来我就忘了

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

    1 条回复 最后回复
  • 东方白杨东 离线
    东方白杨东 离线
    东方白杨
    写于 最后由 编辑
    #106

    @李东岳 方法都是一样的,只不过算的更多时间更久,疯狂造超算资源(滑稽滑稽)

    CAS-IET

    1 条回复 最后回复
  • 风 离线
    风 离线
    风云5091
    在 中回复了 学流体的小明 最后由 编辑
    #107

    @学流体的小明 大佬,请问你们y+、u+沿某一条轴上的分布曲线是怎么画出来的啊,openfoam自带的后处理函数yPlus貌似只能得到壁面第一层网格上的y+,其余部分都设置的是0:laile:

    学流体的小明学 1 条回复 最后回复
  • 学流体的小明学 离线
    学流体的小明学 离线
    学流体的小明 神
    在 中回复了 风云5091 最后由 编辑
    #108

    @风云5091 我们这些都是槽道流里面的实际坐标,除以特征长度,得到的无量纲化坐标y+。
    你如果只做验证的话,OF输出的yPlus足够了,应该可以在paraview显示出来。
    OF输出的yPlus确实只存在于那个边界patch上。

    风 1 条回复 最后回复
  • 风 离线
    风 离线
    风云5091
    在 中回复了 学流体的小明 最后由 编辑
    #109

    @学流体的小明 谢谢大佬解答。我是做可压缩湍流(平板)边界层内的热化学非平衡效应的,最近要验证openfoam里的求解器,希望能算的准边界层内的流动,应该要用到整个流场里面的u+、y+,我看您上面发了一个网盘链接有直接算摩擦速度的matlab程序,但是链接显示过期了,麻烦能再发下吗,我想学习参考一下。另外大佬之前有用过rhocentralFoam吗,我用的hy2Foam非平衡程序是基于这个求解器开发的,不知道这个求解器算边界层怎么样?

    风 1 条回复 最后回复
  • 风 离线
    风 离线
    风云5091
    在 中回复了 风云5091 最后由 编辑
    #110

    @风云5091 刚刚仔细找了一下,发现大佬后面又发了一下链接,能正常下载了,目前就第二个问题需要大佬解答:papa:

    1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    写于 最后由 编辑
    #111

    根据《无痛苦NS方程笔记》,对数率公式是有个E系数,一般取9.8。《An introduction to computational fluid dynamics: the finite volume method》也是这种写法。
    22bed13e-0ca9-45bf-a614-4900eeba6b40-image.png


    约翰霍普金斯湍流数据库$Re_\tau=1000$结果,包括上文对$Re_\tau=395$结果对数律对比,都是将E=1来考虑。

    https://turbulence.pha.jhu.edu/docs/README-CHANNEL.pdf
    681316e2-d420-46db-92df-3cb952cbd8a5-image.png


    @学流体的小明 @李东岳 ,现在问题是:

    1. 为什么最终画图时候,对数率都取E=1,而不用E=9.8?
    2. 根据约翰霍普金斯湍流数据库,DNS计算$Re_\tau=1000$的对数律,其 k 和 C 值并不是取 $\kappa=0.4, -C=5.5$。换句话说,光滑壁面的对数律系数,是否会随着雷诺数的改变而改变,而不是都是固定的值$\kappa=0.4, -C=5.5$?
    3. 现在有个观点是,of自带算例 channel395,后处理计算$U^+$和$y^+$时,应该用实际计算出来的$u_\tau$,而不能用理论$Re_\tau=395$的反算值。原因是自带算例 channel395 实际模拟得到的$Re_\tau$并不等于395,相当于实际得到的是其他雷诺数的充分发展湍流,用模拟的$u_\tau$计算$U^+$和$y^+$,其剖面在粘性子层依然满足线性关系,但在对数律上公式系数应当类似约翰霍普金斯湍流数据库$Re_\tau=1000$结果重新取值,使得满足另外的对数律关系,而不是沿用$\kappa=0.4, -C=5.5$
    学流体的小明学 C 2 条回复 最后回复

  • 登录

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