@cyw 凝结和气化之差不就已经是体积源了吗,为啥还加vDotcAlpha
飞火流星jyj
帖子
-
-
@cyw 为啥mDot第二项没负号,你少写了吗,还有就是不就是再乘了个体积分数吗?看对不对得对照它的方程看
-
@cyw 代码里只有mDotAlpha()质量源项啊,没有mDot
-
@东岳 在新版本5.0或者V1712的interfoam求解器内都有这个头文件
-
其中zeroField Su;
zeroField Sp;
zeroField divU;
都是零场的作用不清楚,是为了用户以后在这里自行修改吗? -
如何理解最后两项???
是fvm::Sp(vDotvmcAlphal,alpha1)+vDotcAlphal
不是单一的vDotvmcAlphal
,自己看代码好像vDotvmcAlphal
已经代表了m*(1/rhol-alpha1*(1/rhol-1/rhov))
,这不就是公式变形后的最后一项吗?求老师、大神指导 -
fluent中设置是这样的:
Boundary A, B and C, D
are set as periodic boundary conditions. Boundary E is set as no slip
wall. Boundary F is set as opening pressure boundary.
不知道在OF怎样设置好FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 100000; boundaryField { inlet { type cyclic; } gasinlet { type cyclic; } outlet1 { type cyclic; } outlet2 { type cyclic; } left { type cyclic; } right { type cyclic; } TWall { type fixedFluxPressure; } atmosphere { type totalPressure; p0 uniform 100000; U U; phi phi; rho rho; psi none; gamma 1; value uniform 100000; } } dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type cyclic; } gasinlet { type cyclic; } outlet1 { type cyclic; } outlet2 { type cyclic; } left { type cyclic; } right { type cyclic; } TWall { type fixedValue; value uniform (0 0 0); } atmosphere { type pressureInletOutletVelocity; value uniform (0 0 0); } }
体积项用setField,但这样不对,求大神帮忙:big_mouth:
-
@wwzhao 谢谢
-
@小龙 filter 搜一下就可以找到
-
@random_ran 好的
这是没切之前的,网格很正常
slice之后
而且cell变成了正常截面的两倍 -
为啥paraview做了slice后,对应的cell数目变成了两倍,而且选择surface and edge 显示,发现网格中间有斜杠了
-
@mohui 我是基于那个方程进行修改添加一些,只不过它本身存在一些问题,当我设置1的时候温度正常
-
如果不设置nAlphaSubCycles=1的话,会出现超过最大温度,数值很越来越大,不稳定情况,设置了就正常了,增加温度方程的方法网上的interTempFoam,有没有遇到同样的状况,请大神分析一下
-
@李东岳 李老师,是在说paraview里选取网格的point嘛?
-
想知道网格在模型里按照什么顺序来排列的
-
@李东岳
1.volScalarField Source ( IOobject ( "Source", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh );
加了,mesh才能通过编译,
2. Source[celli]是double类型的,是不是类似于对应网格的Source.value(),
3. volScalarField又是定义单位的呢,是不是在速度方程后添加==Source;并且保证量纲统一就可以达到更新速度的目的
4. 我还不是很理解volScalarField和scalarField的区别李老师我问题比较多,看您方便回答哪个,十分感谢
-
@李东岳 李老师我的是3.0.X版本的,
提示:
UEqn.H:68:6: error: ‘Source’ was not declared in this scope
Source 没有指定,这个要怎么创建啊 -
@李东岳 那要怎么写呢,李老师能请简单指点一下:happy:
-
@李东岳 李老师,就是添加源项的速度方程只想对于某种网格求解以此来更新这个网格内的速度
-
对于单个网格如何写速度方程,想更新单个网格的速度,UEqn好像是求解了所有的
-
是不是相中再表示多个成分,目前看VOF模型里还没有
-
@队长别开枪 十分感谢,对于
多个网格的U[celli]的保存,我想用定义一个速度单位的数组存起来,dimensionedvector 好像只是保存一个值吧,这样的数组我不知道该写好,请大神指教 -
@队长别开枪 请问如果在forAll(,celli)循环里,mesh.c()是指循环到的网格的中心嘛,还是所有的,如果不是怎样写循环到的celli的中心坐标:crying: 求大神指教
-
@xpqiu 解决了:happy:
-
@xpqiu版主, 我定义了一个变量然后赋值
U2=U[celli]
计算报错:
Different dimensions for =
dimensions : [0 1 -1 0 0 0 0] = [0 0 0 0 0 0 0]From function dimensionSet::operator=(const dimensionSet&) const in file dimensionSet/dimensionSet.C at line 171.
怎么U[celli]成无单位的量了,不是应该是单位网格的速度
-
@程迪 请问怎么查看是0,1,2还是1,2,3这样的编号:confused:
-
@程迪 好的,一个是polymesh的一个是扩展的FvMesh,那我该用哪个表示出某网格的x,y,z,这该怎么写啊
-
@队长别开枪 好的,十分感谢
-
@队长别开枪 请问您测试过吗,我测试了在
createField.H
加入Info << "\n-Class: cell----------" << endl; cell &cl(cells[100]); Info << "cl.centre(points, faces) = " << cl.centre(points, faces) << endl;
报错如下:
createFields.H:33:14: error: ‘cells’ was not declared in this scope cell &cl(cells[100]); ^ createFields.H:34:52: error: ‘points’ was not declared in this scope Info << "cl.centre(points, faces) = " << cl.centre(points, faces) << endl; ^ createFields.H:34:60: error: ‘faces’ was not declared in this scope Info << "cl.centre(points, faces) = " << cl.centre(points, faces) << endl; ^
怎么修改就好了
-
@李东岳 李老师可以了,3.0.x没问题,生成的数据是怎么排列的呢?
-
@李东岳 李老师我测试了下
Info<< "Reading volume\n" << endl; volScalarField volume ( IOobject ( "volume", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh.V() );
提示错误是:
no matching function for call to ‘Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject, const Foam::DimensionedField<double, Foam::volMesh>&)’ );
没有这种类型
-
@队长别开枪 十分感谢,我回去好好看看:laughing:
-
@赵一铭 赵老师,那网格的坐标呢,x,y,z
-
例如U(i,j,k),请老师指点一下
-
@赵一铭 发现效率很慢啊,赵老师,这个语句出自于of在哪定义的?
-
@赵一铭 谢谢赵老师耐心回答,
forAll(alpha, celli)
这句是怎么解释的 -
-
对于两相流求解器想判断alpha为某值然后两个分支,但if语句不对,不知道怎么写程序语句,求指导:happy:
-
@cfd-china 想要跟饱和压力无关,和饱和温度有关的相变,不知道好改吗
-
@piaston 同求啊
interPhaseChangeFoam相方程
mDot 与 mDotAlpha
interPhaseChangeFoam相方程
新版本的alphaSuSp.H
新版本的alphaSuSp.H
interPhaseChangeFoam相方程
平板降膜边界条件设定(interFoam求解器)
paraview做slice问题
OpenFOAM怎么求对流换热流动中的bulk temperature,对流换热系数和Nusselt数
paraview做slice问题
paraview做slice问题
疑惑对于Adding temperature equation to interFoam
疑惑对于Adding temperature equation to interFoam
如何在paraview中显示mesh编号
如何在paraview中显示mesh编号
如何在求解器编程中编写获取某个网格的速度?
如何在求解器编程中编写获取某个网格的速度?
如何在求解器编程中编写获取某个网格的速度?
如何在求解器编程中编写获取某个网格的速度?
如何在求解器编程中编写获取某个网格的速度?
compressibleInterFoam 耦合多组分传输模型
网格体积
网格体积
如何在求解器编程中编写获取某个网格的速度?
如何在求解器编程中编写获取某个网格的速度?
网格体积
网格体积
网格体积
网格体积
网格体积
网格体积
网格体积
网格体积
如何在求解器编程中编写获取某个网格的速度?
判断alpha==零点几的语句怎么编写
判断alpha==零点几的语句怎么编写
判断alpha==零点几的语句怎么编写
判断alpha==零点几的语句怎么编写
可压缩两相流加入相变后,p方程该如何修改呢???
可压缩两相流加入相变后,p方程该如何修改呢???