CFD中文网

@李东岳 东岳老师您好，我后面设置了一个算例，Re=100，边界条件如下所示：5f6212ed-6a29-41f1-a3ec-a4cb870860a9-image.png ，如果把重力加速度设置为0，计算正常，但是加入重力后，计算到后面就会发散，从云图上的感受是速度堵在出口没办法出去，下面是发散前我截出来的某个时刻的云图：
de9f481a-edf9-45fd-b60a-402d4997321e-9715855a0b8e5273d54b9d49de72baf.png 30d0e1d0-8728-4008-a76b-633adc2be9f1-6c476442c75aacaf9c5e83671ab3e27.png 7b9f80ae-537b-4b5b-a8e3-bab63dd0e9b9-e7e1cf13e2b5c28073eec641983f49e.png a896c296-fce5-41ff-b7f0-283d94ef4036-ac5d8ed87ba20cdbe62a8b20db36a7c.png
能否请教一下老师，帮我分析一下问题可能出在哪呢？:xiezuoye:

@李东岳 好的，我明白啦，感谢李老师~::xinxin:

@东岳 非常感谢东岳老师！按照您的指导，我设置了新的边界条件immersedboundaryEnthalpyFvPatchField，其他部分的设置目前都正常，但是想要使用浸没温度边界上的函数来更新焓值时会报错，我的程序部分如下所示：

const label patchi = patch().index(); wordList type = thermo.T().boundaryField().types(); if(isA<immersedBoundaryFvPatchScalarField>(type[patchi])) { immersedBoundaryFvPatchScalarField& TIb = thermo.T().boundaryField()[patchi]; scalarField & Tref = TIb.refValue(); }

之后程序就会报错，

error: invalid initialization of reference of type ‘Foam::immersedBoundaryFvPatchScalarField& {aka Foam::immersedBoundaryFvPatchField<double>&}’ from expression of type ‘const Foam::fvPatchField<double>’ thermo.T().boundaryField()[patchi];

应该是无法将基类转换为派生类，想和老师请教一下遇到这个问题要怎么解决呢，麻烦老师了！！！:xiexie:

//k rho_*epsilon_/k_ //epsilon rho_*Cmu_*sqr(k_)/epsilon_ 为k_和epsilon_接近0，导致除法崩溃，所以使用bound（）

可以试试增大epsilon 减小k；连续性残差也偏大，可以试试fvSchemes散度项二阶换一阶。

@Jato 我的是在OpenFOAM-2.4.0版本下编译的，没有出现你这个问题