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

Floating Point Exception

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 4, 2024, 09:20
Default Floating Point Exception
  #1
New Member
 
Amir Khan
Join Date: Jan 2023
Posts: 5
Rep Power: 3
amirkhan5997 is on a distinguished road
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


Code:
 #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
Code:
/*--------------------------------*- 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;
    }



}

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


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

Last edited by amirkhan5997; January 9, 2024 at 09:43.
amirkhan5997 is offline   Reply With Quote

Reply

Tags
dynamicmeshdict, floating point exception, openfoam, pimplefoam, vortex induced vibrations

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
icoFoam floating point exception (8) leizhao512 OpenFOAM Running, Solving & CFD 7 November 1, 2018 12:43
A floating point exception has occurred: floating point exception [Overflow]. starlight STAR-CCM+ 4 May 4, 2016 10:08
A floating point exception - SEM Model yansheng STAR-CCM+ 1 April 4, 2016 05:57
Floating point exception from twoPhaseEulerFoam openfoammaofnepo OpenFOAM Running, Solving & CFD 1 March 19, 2016 14:56
Floating point exception (core dumped) for GAMG solver yuhou1989 OpenFOAM Running, Solving & CFD 2 March 24, 2015 20:28


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