[amirkhan5997]

Hello CFD community,


I hope this message finds you all well. I am currently facing a challenging issue in my VIV simulation past a cylinder and would greatly appreciate any insights or suggestions.


In brief, I have been working on running VIV simulations using pimpleFoam, but I consistently encounter a floating point exception error after a few iterations. Despite trying various troubleshooting steps, including changing solvers, adjusting mesh properties, and modifying relaxation factors, the problem persists.



I suspect that the rapid divergence of acceleration and velocity of the cylinder might be contributing to the crash. For a more detailed analysis, I have attached the sixDoFRigidBodyMotionState file, containing relevant information about the acceleration and velocity during the simulation.

I have already sought advice from my professor and changed computers to rule out system-specific issues, but the problem persists. If anyone has faced a similar challenge or has any suggestions on how to resolve this issue, your insights would be immensely valuable.


Here is the case file for your reference. I appreciate your time and assistance in advance.



https://drive.google.com/drive/folde...jD?usp=sharing

 #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::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  ? at ??:?
#8  Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#9  ? in /usr/lib/openfoam/openfoam2006/platforms/linux64GccDPInt32Opt/bin/pimpleFoam
#10  ? in /lib/x86_64-linux-gnu/libc.so.6
#11  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#12  ? in /usr/lib/openfoam/openfoam2006/platforms/linux64GccDPInt32Opt/bin/pimpleFoam
Floating point exception

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1912                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


motionSolverLibs (sixDoFRigidBodyMotion);

motionSolver    sixDoFRigidBodyMotion;


displacementLaplacianCoeffs
{
    diffusivity     uniform 1;
}

dynamicFvMesh       dynamicMotionSolverFvMesh;

dynamicOversetFvMeshCoeffs
{
//    layerRelax 0.3;
}

sixDoFRigidBodyMotionCoeffs
{

    accelerationRelaxation 1.0;
    accelerationDamping 1.0;

    patches         (cylinder);
    innerDistance   1.0;
    outerDistance   1.1;

    mass            0.03575;
   centreOfMass    (-0.016 0.0 0.0);
   momentOfInertia (1.144E-08 1.144E-08 1.144E-08);

   g               (0 -9.81 0);
   rho             rhoInf;
   rhoInf          14817.1;
   report          on;

   solver
   {
       type symplectic;
   }

   constraints
   {

    yPlan
    {
       sixDoFRigidBodyMotionConstraint   plane;
       normal ( 0 1 0); 
    }    

    yAxis
    {
        sixDoFRigidBodyMotionConstraint orientation;
        centreOfRotation                (0 1 0);
    }
}


restraints
{
    spring
    {
        sixDoFRigidBodyMotionRestraint linearSpring;

        anchor          (-0.016 0.0 0.0);
        refAttachmentPt (-0.016 0.0 0.0);
        stiffness       102.97;
        damping         0.00;
        restLength      0;
    }



}

}
// ************************************************************************* //

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2006                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      binary;
    class       dictionary;
    arch        "LSB;label=32;scalar=64";
    location    "0.002/uniform";
    object      sixDoFRigidBodyMotionState;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

centreOfRotation ( -0.0157166 0 -2.80571e-25 );

orientation     ( 1 0 0 0 1 0 0 0 1 );

velocity        ( 88504.4 0 -2.05583e-15 );

acceleration    ( 1.77008e+08 -88.517 -4.11165e-12 );

angularMomentum ( 0 0 0 );

torque          ( 6.92773e-14 -1.65298e-13 -0.0450204 );


// ************************************************************************* //