codedfixedvalue使用问题
-
@soulx7 感觉代码应该是这样的
const fvPatch& boundaryPatch = patch(); const vectorField& Cf = boundaryPatch.Cf(); vectorField& field = *this;
参考的这两个帖子
https://cfd-china.com/topic/4765/关于cyclic-boundary-condition的一些问题?_=1671895451562 -
@wangfei9088 你好,更改后问题依然存在,
看报错提示是激活某个东西失败,是不是用这个边界条件需要一些准备条件,不能直接使用?
-
@wangfei9088 分号加上去了,还是一样的报错提示。它报错的代码部分其实是name那里,这一行我去掉它提示缺少字符,加上就提示激活失败,是不是版本问题导致的?我用的是com官网的v2012版本
-
@wangfei9088 你好,现在代码如下图所示:
inlet1 { type codedFixedValue; value uniform (0 0 0); name field; code #{ const fvPatch& boundaryPatch = patch(); const vectorField& Cf = boundaryPatch.Cf(); vectorField& field = *this; const scalar u1=3; const scalar u2=2; const scalar L=0.00048; forAll(Cf,faceI) { const scalar y = Cf[faceI][1]; field[faceI] = vecotr(0,u1+u2*tanh((y-0.05)/L),0); } operator==(field); #}; } inlet2 { type codedFixedValue; value uniform (0 0 0); name field; code #{ const fvPatch& boundaryPatch = patch(); const vectorField& Cf = boundaryPatch.Cf(); vectorField& field = *this; const scalar u1=3; const scalar u2=2; const scalar L=0.00048; forAll(Cf,faceI) { const scalar y = Cf[faceI][1]; field[faceI] = vecotr(0,u1+u2*tanh((y-0.05)/L),0); } operator==(field); #}; } inlet3 { type codedFixedValue; value uniform (0 0 0); name field; code #{ const fvPatch& boundaryPatch = patch(); const vectorField& Cf = boundaryPatch.Cf(); vectorField& field = *this; const scalar u1=3; const scalar u2=2; const scalar L=0.00048; forAll(Cf,faceI) { const scalar y = Cf[faceI][1]; field[faceI] = vecotr(0,u1+u2*tanh((y-0.05)/L),0); } operator==(field); #}; } inlet4 { type codedFixedValue; value uniform (0 0 0); name field; code #{ const fvPatch& boundaryPatch = patch(); const vectorField& Cf = boundaryPatch.Cf(); vectorField& field = *this; const scalar u1=3; const scalar u2=2; const scalar L=0.00048; forAll(Cf,faceI) { const scalar y = Cf[faceI][1]; field[faceI] = vecotr(0,u1+u2*tanh((y-0.05)/L),0); } operator==(field); #}; } outlet { type zeroGradient; } fixedWalls { type zeroGradient; } frontAndBack { type empty; }
现在报错提示也发生了变化,变成了
,看红字高亮部分是说"vector"没有声明,是代码里漏掉了什么东西吗?
-
@wangfei9088 你好,哈哈昨天晚上搞得太晚神智有点不清了,刚改了编译通过了,多谢