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. CodedFixedValue边界赋值后,noslip壁面变为滑移边界?

CodedFixedValue边界赋值后,noslip壁面变为滑移边界?

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

    今天在做模拟的时候发现,对入口进行CodedFixedValue边界赋值后,noslip壁面的速度不为0,换成FixedValue边界后又正常了,这是为什么呢?
    Uboundary condition

    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volVectorField;
        object      U;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [0 1 -1 0 0 0 0];
    
    internalField   uniform (10 0 0);
    
    boundaryField
    {
        inlet
        {
            type            codedFixedValue;
            value           uniform (10 0 0);
            name            codeinlet;
            code            #{
                const fvPatch& boundaryPatch = this->patch();
    
                vectorField U = boundaryPatch.Cf();
                vectorField u = boundaryPatch.Cf();
                
                scalar zstar = 0.003;
                scalar kappa = 0.41;
                scalar ustar = 0.726;
                scalar ubar  = 10;
    
                forAll(boundaryPatch, faceI)
                {
                //spatial coordinates, type: scalar
                scalar y = boundaryPatch.Cf()[faceI].y();
                            
                U[faceI] = vector((ustar/kappa)*log(1+y/zstar), 0, 0);
                u[faceI] = 0.1*ubar*exp(-0.5*y)*vector(rand() % 2 - 1, rand() % 2 - 1, rand() % 2 - 1);
                U[faceI] = U[faceI]+u[faceI];
                if (U[faceI].x() < 0)
                {
                    U[faceI].x() = -U[faceI].x();
                }
                }
                
                operator==(U+u);
            #};
        }
    
        outlet
        {
            type            zeroGradient;
        }
    
        land
        {
            type            noSlip;
        }
    
        symmetry
        {
            type            slip;
        }
    }
    

    这是计算得到的流场
    51f92c91-e1be-4eb3-95ff-d41e0212ca4d-image.png
    案例文件在这:case.zip

    V 1 条回复 最后回复
  • V 离线
    V 离线
    veen
    在 中回复了 veen 最后由 编辑
    #2

    @veen 破案了,应该是paraview读取case的问题openfoam-post-processing

    1 条回复 最后回复

  • 登录

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