CFD中文网

芋

@18190691340 是的，谢谢老师，我发现是我的0文件夹下的G的参数量纲定义错了，是从别的算例文件那一起考过来的，因为我的计算都是无重力情况下的，所以一直没考虑过是这个文件出了问题，我还以为是求解器中的变量量纲定义错了

芋

@ice_flow 好的谢谢老师，我发现问题了，是我的0文件夹下的G的量纲定义错了，但是之前在用自己改写的求解器计算的时候因为都是考虑无重力的情况，所以一直没想到是这个文件出了问题，非常感谢老师们的指导

芋

@ice_flow 谢谢老师，我今天才看见回复，我按照您的回复对量纲进行了输出

Info << "phiFluid[i] dimensions: " << phiFluid[i].dimensions() << endl;
Info << "alphaEff dimensions: " << alphaEff.dimensions() << endl;
Info << "rad.ST dimensions: " << rad.ST(rhoCp, T).ptr()->dimensions() << endl;

输出结果是这样的，辐射项的量纲不知道为啥输出不出来

phiFluid[i] dimensions: [0 3 -1 0 0 0 0]
alphaEff dimensions: [0 2 -1 0 0 0 0]
rad.ST dimensions: 

--> FOAM FATAL ERROR: 
LHS and RHS of - have different dimensions
     dimensions : [1 0 -3 1 0 0 0] - [1 -1 -3 0 0 0 0]

芋

求助各位大佬，本人刚学习open foam没多久，现在是在多域求解器chtMultiRegionFoam的基础上进行修改，需要把它的原能量方程改写成温度方程，

   while(pimple.correctNonOrthogonal())
   {
      fvScalarMatrix TEqn
    (
         fvm::ddt(T)
         +fvm::div(phiFluid[i],T)
         - fvm::laplacian(alphaEff,T)
       ==
      rad.ST(rhoCp , T)
    );  
          TEqn.relax();
          TEqn.solve();
    }

在不加辐射项的时候算例是可以计算的，但是加上辐射项后会报错（报错信息如下），已经困扰我一个月了，不知道是哪个地方的量纲不对

phiFluid[0] dimensions: [0 3 -1 0 0 0 0]


--> FOAM FATAL ERROR: 
LHS and RHS of - have different dimensions
     dimensions : [1 0 -3 1 0 0 0] - [1 -1 -3 0 0 0 0]


    From function Foam::dimensionSet Foam::operator-(const Foam::dimensionSet&, const Foam::dimensionSet&)
    in file dimensionSet/dimensionSet.C at line 521.

FOAM aborting

CFD中文网

芋圆不圆

帖子