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,如何设置一个随机力源项(numerical trip)诱导湍流

壁面湍流LES,如何设置一个随机力源项(numerical trip)诱导湍流

已定时 已固定 已锁定 已移动 OpenFOAM
5 帖子 3 发布者 3.3k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • L 离线
    L 离线
    luofq-sysu
    写于 最后由 编辑
    #1
    The source term is to be applied similarly to a trip-wire in a physical experiment, i.e. in a thin strip of cells close to the wall
    The numerical trip adds randomised fluctuations to the velocity in the momentum equation in the defined region, which acts to trip the boundary layer from laminar to turbulent.
    

    请教各位大佬,在OpenFOAM如何实现:在近壁面特定区域的cellZone网格上,增加一个随机脉动的源项,脉动的强度与为自由流速成比例,比如10%

    1 条回复 最后回复
  • 快 离线
    快 离线
    快乐的山炮
    写于 最后由 编辑
    #2

    插个眼 我也想知道怎么做。

    1 条回复 最后回复
  • 学流体的小明学 离线
    学流体的小明学 离线
    学流体的小明 神
    写于 最后由 编辑
    #3

    你看一下fvOptions中的vectorCodedSource,高级编程添加吧

    L 1 条回复 最后回复
  • L 离线
    L 离线
    luofq-sysu
    在 中回复了 学流体的小明 最后由 编辑
    #4

    @学流体的小明 感谢大佬。
    我用的是ESI版本,找到这个codeSource的介绍了。
    https://doc.openfoam.com/2306/tools/processing/numerics/fvoptions/sources/rtm/coded/
    但还没搞懂这个随机量怎么生成和添加

    1 条回复 最后回复
  • 学流体的小明学 离线
    学流体的小明学 离线
    学流体的小明 神
    写于 最后由 编辑
    #5

    代码方面你照猫画虎就行呀,它这个是给scalar场添加源项,你给速度场U添加就行了,该变的地方变一变。我用的可能和这个不一样。

    codedSource
    {
        type            coded;
        selectionMode   all;
    
        fields          (h); //对哪一个场添加源项,你用 U
        name            sourceTime; //名字随便起一个
    
        codeAddSup
        #{
            const Time& time = mesh().time(); //获取时间        
            const scalarField& V = mesh_.V(); //获取网格的体积
            scalarField& heSource = eqn.source(); //获取h的方程的源项的引用,之后在这个scalarField上面加上具体的量就可以了
    
            // Start time
            const scalar startTime = 2.0;
    
            // Retrieve the x component of the cell centres
            const scalarField& cellx = mesh_.C().component(0);//获取网格的x坐标
    
            // Only apply when we have reached the start time
            if (time.value() > startTime)
            {
                // Apply the source
                forAll(cellx, i)
                {
                    // cell volume specific source
                    heSource[i] += 1e5*sin(200*cellx[i])*V[i]; //forAll所有网格的循环,加上这个源项
                };
            }
        #};
    }
    
    1 条回复 最后回复

  • 登录

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