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. 不均匀温度边界 热流密度

不均匀温度边界 热流密度

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

    目前我只会做到在边界某段添加固定温度,整个边界都使用热流密度条件也是可以的,请问如何在边界的某段添加热流密度边界~比如下面这段0.000333<=x<=0.000666使用热流边界 应该如何修改~

    OUTLET
    {
        type            codedFixedValue;
        name            heated;
        value           uniform 300;
        
            codeInclude
            #{
                #include "fvCFD.H"
            #};
    
            codeOptions
            #{
                -I$(LIB_SRC)/finiteVolume/lnInclude \
                -I$(LIB_SRC)/meshTools/lnInclude
            #};
    
    	    //libs needed to visualize BC in paraview
    	    codeLibs
    	    #{
        	-lmeshTools \
        	-lfiniteVolume
    	    #};
    
            code
            #{
                /*const IOdictionary& d = static_cast<const IOdictionary&>
    			(
                    dict.parent().parent()
                );
                const fvMesh& mesh = refCast<const fvMesh>(d.db());
                const label id = mesh.boundary().findPatchID("OUTLET");
                const fvPatch& patch = mesh.boundary()[id];*/
    
                
                scalarField T(patch().size(), 300);	
    
                forAll(T, i)
                {
                    const scalar x = patch().Cf()[i][0];
                    //const scalar y = patch.Cf()[i][1];
                    //const scalar z = patch.Cf()[i][2];
                    if(x>=0.000333||x<=0.000666)
                    {
    	         T[i] = 400;
    	     }
    	     else
    	     {
    	         T[i] = 300;
    	     }    
                }
    
                (*this)==T;//T.writeEntry("", os);
            #};
         
    }
    
    1 条回复 最后回复

  • 登录

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