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定义入口速度的问题(DSRFG方法)

LES定义入口速度的问题(DSRFG方法)

已定时 已固定 已锁定 已移动 OpenFOAM
82 帖子 16 发布者 163.3k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    在 中回复了 sunss 最后由 编辑
    #37

    sunss邮件发给我了一些非常漂亮的数据图,应该是对算法这面很熟悉,现在卡在openfoam设置这里。你们多分享互相学习。

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

    1 条回复 最后回复
  • X 离线
    X 离线
    xjwang
    在 中回复了 sunss 最后由 编辑
    #38

    @sunss
    这个设置不是很难 你可以看一下这个视频 以及Eddylicious这个软件的document,他们里边也是用的这个边界条件
    https://www.bilibili.com/video/BV1CV411C7Jr

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

    @李东岳
    李老师 突然想起来,您这个图是fluent的图还是openfoam的图?如果是openfoam的图能不能提示一下怎么画出来的?用的什么filter?

    李东岳李 1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    在 中回复了 xjwang 最后由 编辑
    #40

    我那个图是截取 @霜染丹枫 的图,看起来是Fluent

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

    X 1 条回复 最后回复
  • X 离线
    X 离线
    xjwang
    在 中回复了 李东岳 最后由 编辑
    #41

    @李东岳
    哦哦,因为知道您不用fluent 看到您贴出这个图 我还以为paraview的后处理也能这么画!

    霜 1 条回复 最后回复
  • 霜 离线
    霜 离线
    霜染丹枫
    在 中回复了 xjwang 最后由 编辑
    #42

    @xjwang 朋友,是用fluent画的图。我后来更新了结果,在前面的帖子里放上了结果,你看下结果如何?

    X 1 条回复 最后回复
  • S 离线
    S 离线
    sunss
    在 中回复了 xjwang 最后由 编辑
    #43

    @xjwang 在 LES定义入口速度的问题(DSRFG方法) 中说:

    @sunss
    这个设置不是很难 你可以看一下这个视频 以及Eddylicious这个软件的document,他们里边也是用的这个边界条件
    https://www.bilibili.com/video/BV1CV411C7Jr

    你说的这个视频和工具包,我之前都看了,没能找到解决方法。我简化一下问题,假设inlet面上只有四个点,分别是(1 1 0)(1 2 0)(2 1 0)(2 2 0),我要导入两个时间步的速度信息,也就是总共有八个速度的信息,那timeVaryingMappedFixedValue边界条件的point文件和U文件写成下面的形式行不行。point.jpg
    速度信息.jpg

    X 2 条回复 最后回复
  • S 离线
    S 离线
    sunss
    在 中回复了 霜染丹枫 最后由 编辑
    #44

    @霜染丹枫 在 LES定义入口速度的问题(DSRFG方法) 中说:

    @李东岳 如李老师所言:“像是这种速度波动越往后越少的,类似开始有湍流后来无湍流”。我有个疑惑,采用中心差分格式的扩散项可以将扰动向四周传递,迎风格式的对流项具有迁移性,即可以将扰动向下游传递。从上面来看,入口的扰动应该是可以持续向下游传递的。如果几何轴向长度非常长的话,在很远的下游是不是也会保持一定湍流状态?还是由于壁面摩擦力及流体粘性对扰动有削减作用,在下游一直衰减甚至脉动全无,就是这种入口给定方法本身就会出现的现象呢?

    湍流减少的现象应该是湍流生成方法不满足无散度条件,而像满足无散度条件ATSM-L这种方法下游湍流强度不仅不会减少,还会稍微增强。推荐你们去了解一下万嘉伟开发的一款湍流入口生成工具。https://nheri-simcenter.github.io/TinF-Documentation/examples/examples.html
    湍流强度.jpg

    霜 1 条回复 最后回复
  • X 离线
    X 离线
    xjwang
    在 中回复了 sunss 最后由 编辑
    #45

    @sunss
    应该是在constant/boundaryData/inlet文件夹下边有一个points的文件,然后每一个时间步都是一个文件夹,其中包含一个U文件。比如说,我在inlet文件夹下边在terminal中输入ls,会显示:

    0 0.001 0.002 ...... points

    而在每一个时间步文件夹下输入ls,比如在0.001文件夹下,会显示:

    U

    H 1 条回复 最后回复
  • X 离线
    X 离线
    xjwang
    在 中回复了 sunss 最后由 编辑
    #46

    @sunss
    文件格式应该是:
    points:
    .......................................................
    (
    (1 1 0)
    (0 1 0)
    (1 0 0)
    (0 0 0)
    )
    ...........................................................
    U:
    ..........................................................
    (
    (三维速度向量)
    (三维速度向量)
    (三维速度向量)
    (三维速度向量)
    )
    .....................................................
    所以你的设置有两个问题,一个是标出了4或者8这种数量,这个需要删掉,一个是需要把速度分别存储在对应的文件夹下边,不能在一个文件中。还有那些openfoam的星号可以删掉哈哈

    S 1 条回复 最后回复
  • X 离线
    X 离线
    xjwang
    在 中回复了 霜染丹枫 最后由 编辑
    #47

    @霜染丹枫
    不知道你想模拟什么,想和你确认一下是不是模拟ABL,没有得到回复,单凭一张瞬态的速度场图,没法判断你的模拟结果好与坏。

    霜 1 条回复 最后回复
  • 霜 离线
    霜 离线
    霜染丹枫
    在 中回复了 xjwang 最后由 编辑
    #48

    @xjwang 你好,我的研究不是ABL,属于槽道流,但是想用DSRFG这个方法生成湍流入口用于LES计算。

    李东岳李 1 条回复 最后回复
  • 霜 离线
    霜 离线
    霜染丹枫
    在 中回复了 sunss 最后由 编辑
    #49

    @sunss 感谢分享,但是根据huang2010论文中指出该方法严格保证了无散度条件,并且给出了证明过程。

    S 1 条回复 最后回复
  • S 离线
    S 离线
    sunss
    在 中回复了 xjwang 最后由 编辑
    #50

    @xjwang 在 LES定义入口速度的问题(DSRFG方法) 中说:

    @sunss
    文件格式应该是:
    points:
    .......................................................
    (
    (1 1 0)
    (0 1 0)
    (1 0 0)
    (0 0 0)
    )
    ...........................................................
    U:
    ..........................................................
    (
    (三维速度向量)
    (三维速度向量)
    (三维速度向量)
    (三维速度向量)
    )
    .....................................................
    所以你的设置有两个问题,一个是标出了4或者8这种数量,这个需要删掉,一个是需要把速度分别存储在对应的文件夹下边,不能在一个文件中。还有那些openfoam的星号可以删掉哈哈

    多谢!已经成功导入风速数据。

    1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    在 中回复了 霜染丹枫 最后由 编辑
    #51

    @霜染丹枫 所以到现在,你的问题解决了么?

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

    S 霜 2 条回复 最后回复
  • S 离线
    S 离线
    sunss
    在 中回复了 李东岳 最后由 编辑
    #52

    @李东岳 解决了,要感谢@xjwang大佬的帮助。

    1 条回复 最后回复
  • S 离线
    S 离线
    sunss
    在 中回复了 霜染丹枫 最后由 编辑
    #53

    @霜染丹枫 看前面的回复得知你要计算的是槽道流,能量谱应该不需要满足卡门谱,高斯谱可能更合适。TinF这个工具能完美的生成满足无散度条件和高斯谱的湍流(目前这个工具包还在开发中,还不能生成卡门谱的湍流),我有万嘉伟博士给我们讲课的视频,你可以看看。链接:https://pan.baidu.com/s/1iJy7CW6S3iwcDIvzZiAIZg
    提取码:70xw

    X 1 条回复 最后回复
  • 霜 离线
    霜 离线
    霜染丹枫
    在 中回复了 李东岳 最后由 编辑
    #54

    @李东岳 还没有,还在找原因@@

    1 条回复 最后回复
  • 霜 离线
    霜 离线
    霜染丹枫
    写于 最后由 编辑
    #55

    @sunss 感谢分享。我做的是槽道流,能谱是通过槽道流实验得到的(沿轴向的一维谱),通过DSRFG方法生成了入口速度场,得到的入口速度场能谱分布与实验能谱拟合的较好,但是在计算中会出现我上面提到下游脉动减弱的问题。下图是我还原Huang文章中的卡门谱结果,整体看起来方法的实现应该是没有问题。1.jpg
    2.jpg

    X 李东岳李 2 条回复 最后回复
  • X 离线
    X 离线
    xjwang
    在 中回复了 霜染丹枫 最后由 编辑
    #56

    @霜染丹枫
    我做的风场模拟,所以我检察风洞模拟是否正确是通过检察计算域内的风速时程的功率谱是否和目标谱一致,也就是是否和在入口处生成的风速满足同样的能量分布。
    不知道你们槽流最关心的是哪一个数值特征。比如,你可以验证turbulence intensity也可以验证power spectra,但不能通过一个瞬态的流场来判断模拟的好坏。一直不清楚你通过离散能量谱得到的脉动速度为什么要通过观察瞬时速度场来判断模拟的好坏,就算有很强的紊流特性也不能保证这个速度场的能量分布就和你目标谱一样。
    上文中@sunss 同学说的那个tinf,你可以看一下,其中有dfsem, dfm等方法,特别是对于槽流,openfoam的tutorial中有一个RE395的算例(也许把名字记错了),应用的是DFSEM方法。
    如@sunss 同学所说,也许高斯谱更适合你的功率谱,这样的话你可以直接应用fluent自带的那个对应RFG(Smirnov et al,2001)的生成方法,具体名字不记得了,你可以自己去查一下。
    在其他回复中我提到过,DSRFG方法虽然理论上来说是无散度的,但如果应用在数值模拟中,就不是严格的无散度了,需要考虑网格尺寸,所以我一开始就问你网格尺寸和紊流积分尺度的数值。我们完成了一篇这样的论文,希望早点可以和你们分享。

    1 条回复 最后回复

  • 登录

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