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. 各位CFDer,请问速度的依时类边界条件与人工湍流如何同时施加呢?

各位CFDer,请问速度的依时类边界条件与人工湍流如何同时施加呢?

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

    如题,我想实现一个依时类边界与人工湍流的耦合,比如以下两种。请问各位大佬该如何实现呢?给个提示也行:chigua:
    依时类边界:

    type            timeVaryingUniformFixedValue;
    fileName        "$FOAM_CASE/150MPa-U";
    outOfBounds     clamp;
    

    人工湍流:

    value             uniform (600 0 0);
    type              decayingTurbulenceInflowGenerator;
    direction         1;
    LField            uniform 0.00011;//integral scale = width
    refField          uniform (600 0 0);
    RField            uniform (0.0729 0 0 0.0729 0 0.0729);
    Umean             (600 0 0);
    
    C 1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    在 中回复了 strive_DUT 最后由 编辑
    #2

    @strive_DUT 这里decayingTurbulenceInflowGenerator是入口随时间变化的湍流边界,还是生成一个三维的湍流初始场?如果是也是随时间变化的入口边界,那以人工湍流为基础,把timeVaryingUniformFixedValue的东西加进来编译,就在最后根据不同坐标点生成速度把两个方法的速度相加

    strive_DUTS 1 条回复 最后回复
  • strive_DUTS 离线
    strive_DUTS 离线
    strive_DUT
    在 中回复了 coolhhh 最后由 编辑
    #3

    @coolhhh 您好!这里的decayingTurbulenceInflowGenerator是生成的随机湍流入口边界,与时间有关。比如它其中一个版本的源码链接文本。

    C 1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    在 中回复了 strive_DUT 最后由 编辑
    #4

    @strive_DUT 不太清楚这里依时类边界跟这个人工湍流具体的区别。我的理解是一般是采用平均速度+湍流施加入口边界,比如以turbulentDFSEMInlet为基准,在代码最后根据每个入口faceI生成速度U[faceI]地方,把想要叠加的平均风速的一起加进去。

    boundaryConditions/turbulentDFSEMInlet/turbulentDFSEMInletFvPatchVectorField.C

            if (singleProc_ || !Pstream::parRun())
            {
                forAll(U, faceI)
                {
                    U[faceI] += c*uDashEddy(eddies_, Cf[faceI]);
                }
            }
    
    strive_DUTS 1 条回复 最后回复
  • strive_DUTS 离线
    strive_DUTS 离线
    strive_DUT
    在 中回复了 coolhhh 最后由 编辑
    #5

    @coolhhh 您提到的这个turbulentDFSEMInlet应该跟InflowGenerator差不多,我看是在V版本中的。但是我的速度分布跟时间相关,我想在这个基础上人为的添加随机扰动。

    C 1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    在 中回复了 strive_DUT 最后由 编辑
    #6

    @strive_DUT turbulentDFSEMInlet这个边界也是随时间变化的,没太理解速度分布跟时间相关与这种随时间变化的湍流有什么区别。最后算U[faceI]时候公式中增加与当前时间变量公式就可以了

    strive_DUTS 1 条回复 最后回复
  • strive_DUTS 离线
    strive_DUTS 离线
    strive_DUT
    在 中回复了 coolhhh 最后由 编辑
    #7

    @coolhhh 我的速度分布是从文件中读取的,射流入口速度是按照下图给出:
    af01f2e3-ccae-4d62-b53d-7b11456c4a5a-image.png

    C 1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    在 中回复了 strive_DUT 最后由 编辑
    #8

    @strive_DUT 那就可以以decayingTurbulenceInflowGenerator的代码为基准,把timeVaryingUniformFixedValue的中怎么读取这条曲线数据的代码加到一块植入,在最后计算U[faceI]时候把这两部分相加

    strive_DUTS 1 条回复 最后回复
  • strive_DUTS 离线
    strive_DUTS 离线
    strive_DUT
    在 中回复了 coolhhh 最后由 编辑
    #9

    @coolhhh 谢谢您!我尝试做一下。

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

    如果每个网格在同个时刻的平均速度是一样的,更简单的方式是按照C程序写法直接读取这条曲线的离散数据,插值出当前时刻的平均速度,并加到decayingTurbulenceInflowGenerator中的U[faceI]

    1 条回复 最后回复

  • 登录

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