 |
 |
 |
|
|
|
[Sponsors] | ||||
May 18, 2022, 12: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 12:16. |
|
|
|
|
|
May 19, 2022, 12: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, 12:33
|
image
|
#3 |
|
New Member
Jack
Join Date: May 2022
Posts: 3
Rep Power: 4 |
Red area is the AMI.
|
|
|
|
|
|
 |
| Tags |
| pimplefoam |
|
|
Linear Mode
