CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Error with multiSolidBodyMotionSolver

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 18, 2022, 13:13
Default Error with multiSolidBodyMotionSolver
  #1
New Member
 
Jack
Join Date: May 2022
Posts: 3
Rep Power: 4
jweron is on a distinguished road
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
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       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;
    }

}


// ************************************************************************* //
DynamicMeshDict
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;
    }
    }
fvSolution:
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);
}

// ************************************************************************* //
fvSchemes:
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.
jweron is offline   Reply With Quote

Old   May 19, 2022, 13:29
Default case file
  #2
New Member
 
Jack
Join Date: May 2022
Posts: 3
Rep Power: 4
jweron is on a distinguished road
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
jweron is offline   Reply With Quote

Old   May 19, 2022, 13:33
Default image
  #3
New Member
 
Jack
Join Date: May 2022
Posts: 3
Rep Power: 4
jweron is on a distinguished road
Red area is the AMI.
Attached Images
File Type: jpg wheelAMI.jpg (40.0 KB, 14 views)
jweron is offline   Reply With Quote

Reply

Tags
pimplefoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 15:12.