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中文网

Z

zhaominghao

@zhaominghao
关于
帖子
4
主题
1
群组
0
粉丝
0
关注
0

帖子

最新

  • 关于codeFixedValue边界条件读取临界值问题
    Z zhaominghao

    @李东岳
    当边界不做判断直接叠加之后是这样的
    fc3c9e32-1d39-4a41-b44a-85c2594da99b-image.png

    bottom
    {
    type codedFixedValue;
    value uniform 101325;
    name face1P;
    code #{
    scalarField& field = *this;
    const label id = this->patch().boundaryMesh().findPatchID("bottom");//在网格中找到inlet
    const volScalarField& pp = this->db().objectRegistry::lookupObject<volScalarField>("p");
    const scalarField& ppp=pp.boundaryField()[id].patchInternalField();
    forAll(field, faceI)
    {

             field[faceI]=ppp[faceI]+100;   
                 
    	}
    #};
    }
    

  • 关于codeFixedValue边界条件读取临界值问题
    Z zhaominghao

    @李东岳 在 关于codeFixedValue边界条件读取临界值问题 中说:

    存在几个0数据场

    没看到啊老铁

    如果给定边界条件是field[faceI]=ppp[faceI]+100

    输出的参数是b6d72e5c-ece5-4072-90b2-dda806370f59-image.png


  • 关于codeFixedValue边界条件读取临界值问题
    Z zhaominghao

    @李东岳 在 关于codeFixedValue边界条件读取临界值问题 中说:

    存在几个0数据场

    没看到啊老铁

    0数据直接输出会报错 所以给的是边界值加了一个大气压


  • 关于codeFixedValue边界条件读取临界值问题
    Z zhaominghao

    在使用codeFixedValue边界条件实现出口无梯度条件时间发现,读取边界网格的临近值pp.boundaryField()[id].patchInternalField(),并将该值赋给边界值时,赋值后网格左侧会有读取不到数值的情况,存在几个0数据场。

    将出口临近值加上一个大气压后输出给边界,得到如下结果
    b0a1d6c8-0cf1-4e00-9e8e-bc6c3aca7760-ea306b61227d4717d769eb97fdc061c.png

    网格前两个数值为赋值后的101325,将网格加密后结果如下
    8e3af657-b7c1-4663-a794-7c85f062027a-38e14074dd0a229cc5670d71143c530.png

    读取到的0数值数量增加。

    有没有同学遇到过类似情况,请问出现这种现象是由什么引起的,如何解决?
    感谢,
    附上边界代码如下

    bottom
        {
            type            codedFixedValue;
            value           uniform 101325;
            name            face1P;
            code            #{				
                 
              scalarField& field = *this;
    
    	      const label id = this->patch().boundaryMesh().findPatchID("bottom");//在网格中找到inlet
                  //const label id2 = this->patch().boundaryMesh().findPatchID("face1");//在网格中找到inlet
                  const volScalarField& pp = this->db().objectRegistry::lookupObject<volScalarField>("p");  
                  const scalarField& ppp=pp.boundaryField()[id].patchInternalField();//;.patchInternalField();
    	      forAll(field, faceI)
    	        {                         
                      if (ppp[faceI]>101325)
                      {
                        field[faceI]=ppp[faceI]*0.95+0.05*101325;    
                      }
                      else
                      {
                       field[faceI]=201325   ;
                      }
    		}
    	#};
        }
  • 登录

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