求助：液池热毛细对流算例，已编译好求解器和算例，但就是和文献自编程序不一样

[李东岳]

液池热毛细对流算例，已编译好求解器和算例，但就是和文献自编程序不一样

非常有意思的东西，具体什么对不上？

[麦迪文]

@东岳 老师你看它自由面温度场一块一块的高低分明，我这努力了好几天，写出的边界条件它温度均匀的扩散出去

[麦迪文]

这是按照文献给出的边界条件

tmp<vectorField > marangoniFvPatchVectorField::snGrad() const
{
return this->patch().deltaCoeffs()*(*this - patchInternalField());
}
// Evaluate the field on the patch
void marangoniFvPatchVectorField::evaluate()
{
    //Info << "entering  marangoniFvPatchVectorField::evaluate()" << endl;
    if (!this->updated())
    {
	//Info << "marangoniFvPatchVectorField::evaluate(): calling updatecoeffs" << endl;
        this->updateCoeffs();
    }

     vectorField nHat = this->patch().nf();
     vectorField pif = this->patchInternalField();
     scalarField deltas=this->patch().deltaCoeffs();

 	fvPatchField<vector> gradT =this->patch().lookupPatchField<volVectorField, vector>("gradT");
 	vectorField  gradT_internal = gradT.patchInternalField();
	vectorField gradTplane= transform(I - sqr(nHat),gradT_internal);
 	vectorField pifplane= transform(I - sqr(nHat),pif);
 	fvPatchField<vector> gradC =this->patch().lookupPatchField<volVectorField, vector>("gradC");
 	vectorField  gradC_internal = gradC.patchInternalField();
 	vectorField gradCplane= transform(I - sqr(nHat),gradC_internal);

	vectorField result=pifplane+(tmarangonicoeff_*gradTplane+cmarangonicoeff_*gradCplane)/deltas;

 	vectorField::operator=
 	    (
	     result
 	     );
 

    transformFvPatchVectorField::evaluate();
}

[麦迪文]

哎我去，是不是应该用updatecoff隐式赋值边界条件

[李东岳]

信息量好大..

[麦迪文]

虽然没人关注但是还是做下最近研究的总结，以备后来人作参考少走弯路，经多次验证比较，文献中图
可能为速度场，表面温度场怎么也出不来图中所示效果，代码有点小问题，不过不大总结完毕

[dzw05]

@麦迪文 文献竟然还有这种操作，给作者写个comment吧。

[麦迪文]

近期在研究两相欧拉模型和做实验比较忙，稍后会把openfoam模拟的结果温度云图和速度云图传上来供大家分析下

[LiuZhuang]

@麦迪文 等着看云图

[麦迪文]