|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Amir Khan
Join Date: Jan 2023
Posts: 5
Rep Power: 3 ![]() |
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 08:43. |
|
![]() |
![]() |
![]() |
Tags |
dynamicmeshdict, floating point exception, openfoam, pimplefoam, vortex induced vibrations |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
icoFoam floating point exception (8) | leizhao512 | OpenFOAM Running, Solving & CFD | 7 | November 1, 2018 11:43 |
A floating point exception has occurred: floating point exception [Overflow]. | starlight | STAR-CCM+ | 4 | May 4, 2016 09:08 |
A floating point exception - SEM Model | yansheng | STAR-CCM+ | 1 | April 4, 2016 04:57 |
Floating point exception from twoPhaseEulerFoam | openfoammaofnepo | OpenFOAM Running, Solving & CFD | 1 | March 19, 2016 13:56 |
Floating point exception (core dumped) for GAMG solver | yuhou1989 | OpenFOAM Running, Solving & CFD | 2 | March 24, 2015 19:28 |