CFD中文网

在使用overPimpleDyMFoam计算旋转机械 时，出现了如下报错：

Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicOversetFvMesh Selecting motion solver: rigidBodyMotion Selecting rigidBodySolver Newmark PIMPLE: Operating solver in PISO mode Reading field p Reading field U Reading/calculating face flux field phi Creating cellMask field to block out hole cells inverseDistance : detected 2 mesh regions zone:0 nCells:29272 voxels:(193 193 1) bb:(-300.001077033 -200.001077033 1.99892296657) (700.001077033 200.001077033 3.00107703343) zone:1 nCells:8145 voxels:(193 193 1) bb:(-50.0001414007 -49.9939730307 1.99985859933) (50.0001414007 49.9720506607 3.00014140067) --> FOAM Warning : From function bool Foam::oversetPolyPatch::master() const in file oversetPolyPatch/oversetPolyPatch.C at line 149 The master overset patch is not the first patch. Generally the first patch should be an overset patch to guarantee consistent operation. Overset analysis : nCells : 37417 calculated : 37186 interpolated : 204 (interpolated from local:204 mixed local/remote:0 remote:0) hole : 27 Creating interpolatedCells field Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kOmegaSST Selecting patchDistMethod meshWave RAS { RASModel kOmegaSST; turbulence on; printCoeffs on; alphaK1 0.85; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.856; gamma1 0.555555555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; decayControl false; kInf 0; omegaInf 0; } Reading/calculating face velocity Uf No MRF models present No finite volume options present Courant Number mean: 0.0195783147649 max: 21.9831544632 Starting time loop Courant Number mean: 0.0195783147649 max: 21.9831544632 Time = 0.005 forces forces: Not including porosity effects --> FOAM FATAL ERROR: Could not find rho:rho From function void Foam::functionObjects::forces::initialise() in file forces/forces.C at line 204. FOAM exiting

使用的动网格文件如下：

motionSolverLibs ("librigidBodyMeshMotion.so"); dynamicFvMesh dynamicOversetFvMesh; motionSolver rigidBodyMotion; // rigidBodyMotionCoeffs - not needed any more // { report on; solver { type Newmark; gamma 0.1; // Velocity integration coefficient beta 0.1; // Position integration coefficient } accelerationRelaxation 1.0; //- prescribedRotation requires some sub-iterations to converge faster // to desired value. nIter 3; bodies { rotor { type rigidBody; parent root; centreOfMass (-0.265 0 0); mass 1; rho 1; inertia (0.001 0 0 0.001 0 0.001); transform (1 0 0 0 1 0 0 0 1) (-0.265 0 0); patches (BLADE1); innerDistance 1000; outerDistance 2000; joint { type Rz; } } } restraints { rotorRotation { type prescribedRotation; body rudder; referenceOrientation (1 0 0 0 1 0 0 0 1); axis (0 0 1); // Axis of rotation omega sine; // Function1 entry omegaCoeffs { frequency 4; amplitude 0.2; scale (0 1 0); level (0 1 0); } } }

但是如果把动网格文件改为：

dynamicFvMesh dynamicOversetFvMesh; dynamicOversetFvMeshCoeffs { // layerRelax 0.3; } solver multiSolidBodyMotionSolver; multiSolidBodyMotionSolverCoeffs { rotor { solidBodyMotionFunction rotatingMotion; rotatingMotionCoeffs { origin (0 0 0); axis (0 0 1); omega 1.08; } } }

就可以正常计算了。请问这是哪里出现错误了，跪求各位大神解答。

openfoam手册中写到mergeMeshes可以将两套网格合并，那么具体是怎么操作的呢？controlDict中需要添加什么信息？

@dzw05
k：

dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.081; boundaryField { inlet { type zeroGradient; } outlet { type zeroGradient; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } blades { type kqRWallFunction; value $internalField; } }

nut

dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedValue; value $internalField; } outlet { type calculated; value uniform 0; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } blades { type nutkWallFunction; value uniform 0; } }

omega

dimensions [0 0 -1 0 0 0 0]; internalField uniform 0.095; boundaryField { outlet { type zeroGradient; } blades { type omegaWallFunction; value $internalField; } inlet { type zeroGradient; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } }

p

dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } blades { type zeroGradient; } }

u

dimensions [0 1 -1 0 0 0 0]; internalField uniform (14 0 0); boundaryField { inlet { type fixedValue; value uniform (14 0 0); } outlet { type zeroGradient; } frontAndBack { type symmetry; } inter2 { type cyclicAMI; } inter1 { type cyclicAMI; } blades { type rotatingWallVelocity; origin (0 0 0); axis (0 0 1); omega 0; } }

boundary

inlet { type patch; nFaces 400; startFace 6377766; } outlet { type patch; nFaces 400; startFace 6378166; } frontAndBack { type symmetry; inGroups 1(symmetry); nFaces 4800; startFace 6378566; } blades { type wall; inGroups 1(wall); nFaces 80334; startFace 6383366; } inter1 { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 13994; startFace 6463700; matchTolerance 0.0001; transform noOrdering; neighbourPatch inter2; } inter2 { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 11572; startFace 6477694; matchTolerance 0.0001; transform noOrdering; neighbourPatch inter1; }

fvSchemes

ddtSchemes { default Euler; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) Gauss linear; } divSchemes { default none; div(phi,U) Gauss linearUpwind grad(U); div(phi,k) Gauss limitedLinear 1; div(phi,omega) Gauss limitedLinear 1; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited corrected 0.5; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; }

fvSolution

solvers { "pcorr.*" { solver GAMG; tolerance 0.02; relTol 0; smoother GaussSeidel; } p { $pcorr; tolerance 1e-7; relTol 0.01; } pFinal { $p; tolerance 1e-7; relTol 0; } "(U|k|omega)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-06; relTol 0.1; } "(U|k|omega)Final" { $U; tolerance 1e-06; relTol 0; } cellDisplacement { solver GAMG; tolerance 1e-5; relTol 0; smoother GaussSeidel; } } PIMPLE { correctPhi yes; nOuterCorrectors 2; nCorrectors 2; nNonOrthogonalCorrectors 5; pRefCell 0; pRefPoint 0; pRefValue 0; } relaxationFactors { fields { p 0.2; } equations { "(U|k|omega)" 0.7; "(U|k|omega)Final" 0.7; } } cache { grad(U); }

@东岳 这个模型是根据实际的物体建的，所以模型重心就直接取了实际物体的重心坐标，没有数学公式。

@fireztw

https://openfoamwiki.net/index.php/Parameter_Definitions_-_dynamicMotionSolverFvMesh

https://forum.cfdac.com/t/topic/660/2

这两个网站可以参考一下~~

据我了解OpenFOAM里面不能同时实现俩个，不过实现应该不难，你需要编程把两个类库结合起来。纯编程问题，因为动网格算法已经有了。

@lv1995 网格这块我不太懂，应该是六面体网格。checkMesh检查出来的非正交性大概50-60，我把非正交修正次数从“0”改为“1”就不发散了，没试过更大的非正交修正系数。看教程上也说一般设置为“0”或“1”。

http://www.cfd-china.com/topic/862/多个网格合并问题/9

@lv1995 不知道是否正确，innerdistance的意思是随着动边界一起运动的网格厚度，不变形，outerdistance是网格变形域的范围。仅供参考。

OF有标准版，什么都不需要改，如果仔细看程序的话。
Wall shear在有fv3的情况下不是特别好而已。

@xpqiu 边界条件刚上来的那个图是什么意思呢？