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. SRF/MRF的使用问题?数值模拟中要定义旋转速度。

SRF/MRF的使用问题?数值模拟中要定义旋转速度。

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

    @李东岳 速度是指线速度吗?OF有类似的自带算例吗?

    C 1 条回复 最后回复
  • C 离线
    C 离线
    CFD中文网
    在 中回复了 小龙 最后由 编辑
    #7

    @小龙 OpenFOAM应该不存在相关现成的算例,你尝试过上面的代码了么?能否运行?:big_mouth:

    CFD中国标准用户测试帐号
    目前由徐笑笑登录

    小 2 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 CFD中文网 最后由 编辑
    #8

    @cfd-china 尝试过了,不能运行

    1 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 CFD中文网 最后由 编辑
    #9

    @cfd-china 0_1485330951563_upload-2a6a3824-2ce9-4a18-9d5e-08a9ae5cfa12 麻烦您帮我看一下这样定义的旋转速度对吗,我对OpenFoam学习的还不够深入,感谢您对我的帮助。

    1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵一铭
    写于 最后由 编辑
    #10

    你这个和壁面旋转速度还不一样,你的壁面旋转是那种平的旋转速度吧?像是顶盖驱动流那种平的速度?

    小 1 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 赵一铭 最后由 编辑
    #11

    @赵一铭 嗯对,就是那种类似于顶盖驱动流的平的旋转速度,

    1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵一铭
    写于 最后由 编辑
    #12

    岳哥那个代码提示什么错误?

    小 1 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 赵一铭 最后由 编辑
    #13

    @赵一铭 提示了很多错误,我一会用电脑运行一下,截图您帮我看看,您先看看我在上面截图发的那个对吗,可以用吗?谢谢您了

    赵 1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵一铭
    在 中回复了 小龙 最后由 编辑
    #14

    @小龙
    好的,最好贴文本上来,代码高亮一下。

    小 4 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 赵一铭 最后由 编辑
    #15

    @赵一铭 嗯嗯,谢谢您了,我以前没发过,以后会注意的

    1 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 赵一铭 最后由 编辑
    #16

    @赵一铭 我一会发个文本上来

    1 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 赵一铭 最后由 编辑
    #17

    @赵一铭

    MOVINGWALL
    {
        type            codedFixedValue;
        name    swirl;
    
        code
        #{
            const vector axis(0, 0, 1);
    
            vectorField v(2.0*this->patch().Cf() ^ axis);
            v.replace(vector::Z, 0);
            operator==(v);
        #};
        value           $internalField;
    }
    

    ···

    1 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 赵一铭 最后由 编辑
    #18

    @赵一铭
    这个是李东岳老师给的代码运行的错误

    Create time
    
    Create mesh for time = 0
    
    
    SIMPLE: convergence criteria
        field p	 tolerance 0.01
        field U	 tolerance 0.001
        field "(k|epsilon)"	 tolerance 0.001
    
    Reading field p
    
    --> FOAM Warning : 
        From function Foam::Field<Type>::Field(const Foam::word&, const Foam::dictionary&, Foam::label) [with Type = double; Foam::label = int]
        in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/Field.C at line 317
        Reading "/home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/p.boundaryField.INLET" from line 27 to line 28
        expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
    --> FOAM Warning : 
        From function Foam::Field<Type>::Field(const Foam::word&, const Foam::dictionary&, Foam::label) [with Type = double; Foam::label = int]
        in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/Field.C at line 317
        Reading "/home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/p.boundaryField.OUTLET" from line 32 to line 33
        expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
    AMI: Creating addressing and weights between 1911 source faces and 1911 target faces
    AMI: Patch source sum(weights) min/max/average = 0.999994, 1, 1
    AMI: Patch target sum(weights) min/max/average = 0.924226, 1, 0.999254
    Reading field U
    
    Using dynamicCode for patch MOVINGWALL on field U at line 80 in "/home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL"
    Creating new library in "dynamicCode/dummy/platforms/linux64GccDPInt32Opt/lib/libdummy_998c505417ee368105406475ef07b1aa535fae1b.so"
    Invoking "wmake -s libso /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/dynamicCode/dummy"
    wmakeLnInclude: linking include files to ./lnInclude
    Making dependency list for source file fixedValueFvPatchFieldTemplate.C
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL: In member function ‘virtual void Foam::dummyFixedValueFvPatchVectorField::updateCoeffs()’:
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:93:13: error: ‘v’ was not declared in this scope
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:93:31: error: unable to find numeric literal operator ‘operator""x’
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:93:31: note: use -std=gnu++11 or -fext-numeric-literals to enable more built-in suffixes
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:93:35: error: unable to find numeric literal operator ‘operator""y’
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:93:35: note: use -std=gnu++11 or -fext-numeric-literals to enable more built-in suffixes
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:93:39: error: unable to find numeric literal operator ‘operator""z’
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:93:39: note: use -std=gnu++11 or -fext-numeric-literals to enable more built-in suffixes
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:90:13: warning: unused variable ‘x’ [-Wunused-variable]
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:91:20: warning: unused variable ‘y’ [-Wunused-variable]
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:92:20: warning: unused variable ‘z’ [-Wunused-variable]
    /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL:87:22: warning: unused variable ‘U’ [-Wunused-variable]
    /opt/openfoam4/wmake/rules/General/transform:8: recipe for target 'Make/linux64GccDPInt32Opt/fixedValueFvPatchFieldTemplate.o' failed
    make: *** [Make/linux64GccDPInt32Opt/fixedValueFvPatchFieldTemplate.o] Error 1
    
    
    --> FOAM FATAL IO ERROR: 
    Failed wmake "dynamicCode/dummy/platforms/linux64GccDPInt32Opt/lib/libdummy_998c505417ee368105406475ef07b1aa535fae1b.so"
    
    
    file: /home/xiaolong/OpenFOAM/xiaolong-4.1/run/xuanzhuan/0/U.boundaryField.MOVINGWALL from line 80 to line 84.
    
        From function void Foam::codedBase::createLibrary(Foam::dynamicCode&, const Foam::dynamicCodeContext&) const
        in file db/dynamicLibrary/codedBase/codedBase.C at line 206.
    
    FOAM exiting
    
    1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵一铭
    写于 最后由 编辑
    #19
    code
            #{
                // const vector axis(1, 0, 0);
                const fvPatch& boundaryPatch = this->patch();
                vectorField v = boundaryPatch.Cf();
                forAll(boundaryPatch, faceI)
                {
                    v[faceI] = vector(boundaryPatch.Cf()[faceI].x(),0,0);//修改这个
                }
    
                //vectorField v(1, 1, 1);
                //v.replace(vector::X, 1.0);
                operator==(v);
            #};
    

    试试这个,把速度初始化为网格边界面的x方向的值了即(value, 0, 0),参考这个定义你自己的函数。

    小 1 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 赵一铭 最后由 编辑
    #20

    @赵一铭
    感谢您,我还有两个问题:
    这个是定义了绕x轴的旋转速度吗?
    假如我的旋转速度是1500r/min,我该怎么换算呢,修改哪里呢?

    李东岳李 1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    在 中回复了 小龙 最后由 编辑
    #21

    @小龙 0_1493221439703_webwxgetmsgimg.jpg

    你这个是需要做一些数学运算的,比如图上这些点的速度和坐标点的位置有关。知道方程之后就好办了。coding永远都不是问题,你试试先把方程写出来。

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

    小 1 条回复 最后回复
  • 小 离线
    小 离线
    小龙
    在 中回复了 李东岳 最后由 编辑
    #22

    @李东岳 好的,谢谢您的帮助

    1 条回复 最后回复

  • 登录

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