bluecfd运算出错
-
使用pimpleDyMFoam求解器计算风机所受转矩，使用kωSST模型，运算第一步就出现如下错误：
Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: solidBody Selecting solid-body motion function rotatingMotion Applying solid body motion to cellZone inter PIMPLE: no residual control data found. Calculations will employ 2 corrector loops Reading field p AMI: Creating addressing and weights between 14390 source faces and 12742 target faces AMI: Patch source sum(weights) min/max/average = 0.570448, 1.45816, 1.00126 AMI: Patch target sum(weights) min/max/average = 0.000958814, 1.18605, 0.993037 Reading field U Reading/calculating face flux field phi 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.555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; } No MRF models present Reading/calculating face velocity Uf No finite volume options present Courant Number mean: 0.0006176 max: 0.225383 Starting time loop forces Forces: Not including porosity effects forces forceCoeffs: Not including porosity effects forceCoeffs forceCoeffs: Not including porosity effects Courant Number mean: 0.0006176 max: 0.225383 Time = 0.001 AMI: Creating addressing and weights between 14390 source faces and 12742 target faces AMI: Patch source sum(weights) min/max/average = 0.570366, 1.45766, 1.00147 AMI: Patch target sum(weights) min/max/average = 0, 1.18298, 0.993152 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0220577, No Iterations 4 GAMG: Solving for pcorr, Initial residual = 0.230843, Final residual = 0.0159595, No Iterations 1 GAMG: Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0 GAMG: Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0 GAMG: Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0 GAMG: Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0 time step continuity errors : sum local = 2.90421e-009, global = -5.82674e-012, cumulative = -5.82674e-012 PIMPLE: iteration 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0228862, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0203193, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0226345, No Iterations 1 Generating stack trace... Backtrace: ZN10StackTraceC1Ev [0x705c1465+0x25] module: D:\blueCFD-Core-2017\ThirdParty-5.x\platforms\mingw_w64GccDPInt32\lib\libstack_trace.dll ZN4Foam5error10printStackERNS_7OstreamE [0x18d1c88+0x218] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll ZN4Foam6sigFpe13sigFpeHandlerEi [0x18d2af3+0x33] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll (No symbol) [0x4047cd] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe _C_specific_handler [0x7ffe9fcc7c58+0x98] module: C:\WINDOWS\System32\msvcrt.dll 0_chkstk [0x7ffe9ff8f7dd+0x11d] module: C:\WINDOWS\SYSTEM32\ntdll.dll RtlWalkFrameChain [0x7ffe9fefd856+0x13f6] module: C:\WINDOWS\SYSTEM32\ntdll.dll KiUserExceptionDispatcher [0x7ffe9ff8e70e+0x2e] module: C:\WINDOWS\SYSTEM32\ntdll.dll ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE [0x1877fa4+0x24] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll (No symbol) [0x42c1e7] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe (No symbol) [0x4361db] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe (No symbol) [0x4363dd] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe (No symbol) [0x4546ca] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe (No symbol) [0x4013f7] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe (No symbol) [0x40152b] module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe BaseThreadInitThunk [0x7ffe9d393dc4+0x14] module: C:\WINDOWS\System32\KERNEL32.DLL RtlUserThreadStart [0x7ffe9ff63691+0x21] module: C:\WINDOWS\SYSTEM32\ntdll.dll求教各位大神，这是哪里的设置出现问题？
-
openfoam报错如下：
#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:? #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:? #5 ? at ??:? #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 ? at ??:? -
@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); }
