fixedFluxPressure边界条件【转载】

转载自：http://blog.sina.com.cn/s/blog_e256415d0102vikh.html

fixedFluxPressure是OpenFOAM较新的一个边界条件，表示边界处压力通量为固定值。在OpenFOAM-2.2.2中，还同时保留了buoyantPressure与fixedFluxPressure两个边界条件，而较新的版本中已经只能使用fixedFluxPressure条件。在CFD论坛中，有详细的描述，（第5，6楼）

http://www.cfd-online.com/Forums/openfoam-solving/82581-i-need-explanations-about-fixedfluxpressure.html

 

我大致总结一下：buoyantPressure设置压力梯度使其适用于浮力流，由于壁面边界通量差值趋于0，因此边界条件近似zeroGradient，buoyantPressure从物理意义上更符合 Archimedes' principle描述的壁面边界；而fixedFluxPressure则是调整压力梯度使得边界通量与速度边界条件指定的保持一致，收敛性更好。此外，fixedFluxPressure边界条件不限于壁面，也可以用于入口与出口，在interFoam中的渠道流案例中，出口设置为fixedFluxPressure。

 

此外，值得注意的是interFoam中的waterChannel范例，计算输出的log文件中会输出各个面的质量通量，如下所示。如果总和不为0，就表示质量通量不守恒，需要检查一下边界条件和初值条件是否设置正确。

 

faceSource inletFlux output:

sum(inlet) for rho*phi = -8.43073

 

faceSource outletFlux output:

sum(outlet) for rho*phi = 8.43846

 

faceSource atmosphereFlux output:

sum(atmosphere) for rho*phi = 0.00525326

 

这些输出是在system/controlDict文件中增加如下代码实现的：

functions

{

inletFlux

{

type faceSource;

functionObjectLibs ("libfieldFunctionObjects.so");

outputControl timeStep;

log true;

// Output field values as well

valueOutput false;

source patch;

sourceName inlet;

operation sum;

 

fields

(

rho*phi

);

}

 

outletFlux

{

$inletFlux;

sourceName outlet;

}

 

atmosphereFlux

{

$inletFlux;

sourceName atmosphere;

}

}