对MRF MixVessel 2D算例的疑问
-
最近玩了一下simpleFoam算例文件夹下面的MixVessel 2D算例,该算例的模型是搅拌器,网格文件中把整个流场都设置成了1个cellzone(rotor),并且在MRFProperties文件中设置了1个相对坐标系(MRF1),我的理解是该坐标系以(0,0,0)为中心,以沿着z轴看是顺时针的方向以104.72rad/s的角速度旋转。
在边界条件设置中,我理解的0/U文件中的相关设置含义是:整个流场和内外的搅拌叶片均随坐标系一起运动,这就很让人困惑了...这是否与实际的搅拌器相悖?还是说我对边界条件的理解有误?
附上MRFProperties文件、0/U文件及paraFoam中显示的搅拌器模型图及计算完毕后的流场矢量图:
MRFProperties:/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object MRFProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // MRF1 { cellZone rotor; active yes; // Fixed patches (by default they 'move' with the MRF zone) nonRotatingPatches (); origin (0 0 0); axis (0 0 1); omega 104.72; } // ************************************************************************* //
0/U:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 4.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { rotor { type noSlip; } stator { type noSlip; } front { type empty; } back { type empty; } } // ************************************************************************* //
搅拌器模型图:
流场矢量图: