|
[Sponsors] |
May 18, 2022, 13:13 |
Error with multiSolidBodyMotionSolver
|
#1 |
New Member
Jack
Join Date: May 2022
Posts: 3
Rep Power: 4 |
Hi, I am getting this error when trying to use multiSolidBodyMotionSolver to model a car with spinning wheels. I am using AMIs to rotate the spokes with rotating wall for the tires. I was able to run the simulation with the wheels set as walls without the AMIs. I also got moveDynamicMesh -checkAMI to run without any errors. I went to the scalarField files and checked out the lines that it references, but I can't figure out what it means.
Error: Code:
Starting time loop Courant Number mean: 0.002384573 max: 0.4172231 deltaT = 1.199041e-05 Time = 1.19904e-05 PIMPLE: Iteration 1 AMI: Creating addressing and weights between 12008 source faces and 12083 target faces AMI: Patch source sum(weights) min/max/average = 0.5016563, 2.19114, 0.9984182 AMI: Patch target sum(weights) min/max/average = 0, 2.006371, 0.9324172 AMI: Creating addressing and weights between 12016 source faces and 12093 target faces AMI: Patch source sum(weights) min/max/average = 0.261173, 2.280045, 0.9985811 AMI: Patch target sum(weights) min/max/average = 0, 2.01081, 0.9328772 GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00727951, No Iterations 6 time step continuity errors : sum local = 6.393334e-08, global = -4.66333e-09, cumulative = -4.66333e-09 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.05987776, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.9936112, Final residual = 0.06513882, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.06094046, No Iterations 1 Generating stack trace... Backtrace: ZN10StackTraceC1Ev [0x626c1855+0x25] module: D:\blueCFD-Core-2020\ThirdParty-8\platforms\mingw_w64GccDPInt32\lib\libstack_trace.dll printStack [0x6c30ae5a+0x23a] (printStack.C:94) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll sigFpeHandler [0x6c30b693+0x33] (signals\sigFpe.C:77) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll _gnu_exception_handler [0x408ed2+0x122] (D:\mingwbuild\mingw-w64-crt-git\src\mingw-w64\mingw-w64-crt\crt\crt_handler.c:270) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe _C_specific_handler [0x7ff901117ff8+0x98] module: C:\Windows\System32\msvcrt.dll 0_chkstk [0x7ff902e7209f+0x11f] module: C:\Windows\SYSTEM32\ntdll.dll RtlRaiseException [0x7ff902e21454+0x434] module: C:\Windows\SYSTEM32\ntdll.dll KiUserExceptionDispatcher [0x7ff902e70bce+0x2e] module: C:\Windows\SYSTEM32\ntdll.dll 00divide [0x6c2bbf94+0x24] (fields\Fields\scalarField\scalarField.C:116) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll divide<Foam::fvPatchField, Foam::volMesh> [0x42dbfb+0x9b] (X:\OpenFOAM-8\src\OpenFOAM\fields\FieldFields\scalarFieldField\scalarFieldField.C:91) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe operator/<Foam::fvPatchField, Foam::volMesh> [0x4397a5+0x275] (X:\OpenFOAM-8\src\OpenFOAM\memory\tmp\tmpI.H:146) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe operator/<Foam::fvPatchField, Foam::volMesh> [0x439f50+0x50] (X:\msys64\mingw64\include\c++\10.2.0\bits\basic_string.h:186) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe 0000main [0x45826c+0x371c] (pEqn.H:1) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe __tmainCRTStartup [0x4013c1+0x231] (D:\mingwbuild\mingw-w64-crt-git\src\mingw-w64\mingw-w64-crt\crt\crtexe.c:337) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe mainCRTStartup [0x4014f6+0x16] (D:\mingwbuild\mingw-w64-crt-git\src\mingw-w64\mingw-w64-crt\crt\crtexe.c:221) module: D:\blueCFD-Core-2020\OpenFOAM-8\platforms\mingw_w64GccDPInt32Opt\bin\pimpleFoam.exe BaseThreadInitThunk [0x7ff902c47034+0x14] module: C:\Windows\System32\KERNEL32.DLL RtlUserThreadStart [0x7ff902e22651+0x21] module: C:\Windows\SYSTEM32\ntdll.dll Code:
/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 8 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; location "0"; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 35.7632 0); boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" outlet { type pressureInletOutletVelocity; value uniform (0 0 0); } road { type fixedValue; value uniform (0 35.7632 0); } rtire { type rotatingWallVelocity; origin (0.825969 2.1757 0.288829); axis (1 0 0); omega 114.008; } ftire { type rotatingWallVelocity; origin (0.825969 -0.33559 0.288829); axis (1 0 0); omega 114.008; } inlet { type fixedValue; value uniform (0 35.7632 0); } car { type noSlip; } foil { type noSlip; } walls { type slip; } plane { type symmetry; } fwheel { type movingWallVelocity; value uniform (0 0 0); } rwheel { type movingWallVelocity; value uniform (0 0 0); } fAMI1 { type cyclicAMI; value $internalField; } fAMI2 { type cyclicAMI; value $internalField; } rAMI1 { type cyclicAMI; value $internalField; } rAMI2 { type cyclicAMI; value $internalField; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2112 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // motionSolverLibs ( "libfvMotionSolvers.so" ); dynamicFvMesh dynamicMotionSolverFvMesh; motionSolver multiSolidBodyMotionSolver; fwAMI { solidBodyMotionFunction rotatingMotion; rotatingMotionCoeffs { origin (0.825969 -0.33559 0.288829); axis (1 0 0); omega 114.008; } } rwAMI { motionSolver solidBody; solidBodyMotionFunction rotatingMotion; rotatingMotionCoeffs { origin (0.825969 2.1757 0.288829); axis (1 0 0); omega 114.008; } } Code:
/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 8 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; smoother DICGaussSeidel; cacheAgglomeration no; maxIter 50; tolerance 1e-5; relTol 0.01; } pFinal { $p; tolerance 1e-6; relTol 0; } "pcorr.*" { $p; tolerance 1e-2; relTol 0; } "(U|k|epsilon)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-6; relTol 0.1; } "(U|k|epsilon)Final" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-6; relTol 0; } } PIMPLE { correctPhi yes; nOuterCorrectors 2; nCorrectors 1; nNonOrthogonalCorrectors 0; } relaxationFactors { "(U|k|epsilon).*" .7; } cache { grad(U); } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 8 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) cellLimited Gauss linear 1; } divSchemes { default none; //div(phi,U) Gauss upwind; div(phi,U) Gauss linearUpwind grad(U); div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited corrected 0.33; } interpolationSchemes { default linear; } snGradSchemes { default limited corrected 0.33; } // ************************************************************************* // Last edited by jweron; May 19, 2022 at 13:16. |
|
May 19, 2022, 13:29 |
case file
|
#2 |
New Member
Jack
Join Date: May 2022
Posts: 3
Rep Power: 4 |
I just tried running the simulation with just one rotating zone and it did not work, so I think the issue is with something else. Ill add my case file here:
https://drive.google.com/file/d/1UOS...ew?usp=sharing |
|
May 19, 2022, 13:33 |
image
|
#3 |
New Member
Jack
Join Date: May 2022
Posts: 3
Rep Power: 4 |
Red area is the AMI.
|
|
Tags |
pimplefoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|