@huangyuhui723 我明白你的意思了。
@huangyuhui723 在 关于在求解器中使用热物理库中的函数 中说:
通过其它计算更新了T之后想对he赋值
如果是求解TEqn,可以是得到温度T后更新he。OpenFOAM里的求解器一般都是求解EEqn,得到he后通过thermo.correct()更新T。
如果已知边界温度T后要对he赋值,应该在求解EEqn之前对he边界赋值,然后求解EEqn。大致代码是这样:
volScalarField& he = thermo.he(); //****************************************//加这一段 forAll(p.boundaryField(), patchi) { forAll(p.boundaryField()[patchi], facei) { he.boundaryFieldRef()[patchi][facei] = function(p.boundaryField()[patchi][facei], T.boundaryField()[patchi][facei]); } } //****************************************// fvScalarMatrix EEqn ( ...... );这个function函数就是T与he的关系,类似上面,混合物的he通过各组分的显焓的质量分数加权得到,显焓用janaf里的关系确定。
供参考。