[noodleluvmay]

Hello all,

I have been trying to run the rgb_overset tutorial case from foamMooring (https://github.com/hfchen20/foamMooring/tree/master), however, it's giving me floating point exception errors:

[0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #2  ? in /lib64/libpthread.so.0
[0] #3  moordyn::StationaryScheme::Step(double&) in ~/OpenFOAM/xbb16141-v2012/platforms/linux64GccDPInt32Opt/lib/libmoordyn.so.2
[0] #4  moordyn::MoorDyn::icStationary() in ~/OpenFOAM/xbb16141-v2012/platforms/linux64GccDPInt32Opt/lib/libmoordyn.so.2
[0] #5  moordyn::MoorDyn::Init(double const*, double const*, bool) in ~/OpenFOAM/xbb16141-v2012/platforms/linux64GccDPInt32Opt/lib/libmoordyn.so.2
[0] #6  Foam::RBD::restraints::moorDynR2::restrain(Foam::Field<double>&, Foam::Field<Foam::SpatialVector<double> >&, Foam::RBD::rigidBodyModelState const&) const at ??:?
[0] #7  Foam::RBD::rigidBodyModel::applyRestraints(Foam::Field<double>&, Foam::Field<Foam::SpatialVector<double> >&, Foam::RBD::rigidBodyModelState const&) const at ??:?
[0] #8  Foam::RBD::rigidBodySolvers::Newmark::solve(Foam::Field<double> const&, Foam::Field<Foam::SpatialVector<double> > const&) at ??:?
[0] #9  Foam::RBD::rigidBodyMotion::solve(double, double, Foam::Field<double> const&, Foam::Field<Foam::SpatialVector<double> > const&) at ??:?
[0] #10  Foam::rigidBodyMeshMotion::solve() at ??:?
[0] #11  Foam::motionSolver::newPoints() at ??:?
[0] #12  Foam::dynamicMotionSolverListFvMesh::update() at ??:?
[0] #13  Foam::dynamicOversetFvMesh::update() at ??:?
[0] #14  ? at ??:?
[0] #15  __libc_start_main in /lib64/libc.so.6
[0] #16  ? at ??:?
[node033:854632] *** Process received signal ***
[node033:854632] Signal: Floating point exception (8)
[node033:854632] Signal code:  (-6)
[node033:854632] Failing at address: 0x189d000d0a68
[node033:854632] [ 0] /lib64/libpthread.so.0(+0x12d10)[0x1497370acd10]
[node033:854632] [ 1] /lib64/libpthread.so.0(raise+0x10f)[0x1497370acbaf]
[node033:854632] [ 2] /lib64/libpthread.so.0(+0x12d10)[0x1497370acd10]
[node033:854632] [ 3] /users/xbb16141/OpenFOAM/xbb16141-v2012/platforms/linux64GccDPInt32Opt/lib/libmoordyn.so.2(_ZN7moordyn16StationaryScheme4StepERd+0x191)[0x14971d176f51]
[node033:854632] [ 4] /users/xbb16141/OpenFOAM/xbb16141-v2012/platforms/linux64GccDPInt32Opt/lib/libmoordyn.so.2(_ZN7moordyn7MoorDyn12icStationaryEv+0x1414)[0x14971d13de64]
[node033:854632] [ 5] /users/xbb16141/OpenFOAM/xbb16141-v2012/platforms/linux64GccDPInt32Opt/lib/libmoordyn.so.2(_ZN7moordyn7MoorDyn4InitEPKdS2_b+0x18f1)[0x14971d14e511]
[node033:854632] [ 6] /users/xbb16141/OpenFOAM/xbb16141-v2012/platforms/linux64GccDPInt32Opt/lib/librigidBodyMooring.so(_ZNK4Foam3RBD10restraints9moorDynR28restrainERNS_5FieldIdEERNS3_INS_13SpatialVectorIdEEEERKNS0_19rigidBodyModelStateE+0x121c)[0x14971e1eb0ec]
[node033:854632] [ 7] /opt/software/openfoam/gcc-8.5.0/v2012/platforms/linux64GccDPInt32Opt/lib/librigidBodyDynamics.so(_ZNK4Foam3RBD14rigidBodyModel15applyRestraintsERNS_5FieldIdEERNS2_INS_13SpatialVectorIdEEEERKNS0_19rigidBodyModelStateE+0x6c)[0x149729050a1c]
[node033:854632] [ 8] /opt/software/openfoam/gcc-8.5.0/v2012/platforms/linux64GccDPInt32Opt/lib/librigidBodyDynamics.so(_ZN4Foam3RBD16rigidBodySolvers7Newmark5solveERKNS_5FieldIdEERKNS3_INS_13SpatialVectorIdEEEE+0x91)[0x14972906b701]
[node033:854632] [ 9] /opt/software/openfoam/gcc-8.5.0/v2012/platforms/linux64GccDPInt32Opt/lib/librigidBodyDynamics.so(_ZN4Foam3RBD15rigidBodyMotion5solveEddRKNS_5FieldIdEERKNS2_INS_13SpatialVectorIdEEEE+0xe0)[0x149729061140]
[node033:854632] [10] /opt/software/openfoam/gcc-8.5.0/v2012/platforms/linux64GccDPInt32Opt/lib/librigidBodyMeshMotion.so(_ZN4Foam19rigidBodyMeshMotion5solveEv+0xa68)[0x14971bb188c8]
[node033:854632] [11] /opt/software/openfoam/gcc-8.5.0/v2012/platforms/linux64GccDPInt32Opt/lib/libdynamicMesh.so(_ZN4Foam12motionSolver9newPointsEv+0x18)[0x14973924e878]
[node033:854632] [12] /opt/software/openfoam/gcc-8.5.0/v2012/platforms/linux64GccDPInt32Opt/lib/libdynamicFvMesh.so(_ZN4Foam29dynamicMotionSolverListFvMesh6updateEv+0x5f)[0x149738e2ff9f]
[node033:854632] [13] /opt/software/openfoam/gcc-8.5.0/v2012/platforms/linux64GccDPInt32Opt/lib/liboverset.so(_ZN4Foam20dynamicOversetFvMesh6updateEv+0xe)[0x149738b7fbce]
[node033:854632] [14] overInterDyMFoam[0x45620f]
[node033:854632] [15] /lib64/libc.so.6(__libc_start_main+0xe5)[0x149736cfe7e5]
[node033:854632] [16] overInterDyMFoam[0x45ff9e]
[node033:854632] *** End of error message ***

I wondered if anyone have any ideas on what could be causing this, please?

Some extra information:
1. I'm running with OF v2012.
2. I'm running on an HPC, IT helped with the compilation of the code - I believe (?) it should be compiled properly - as the case would run with 'moodyR' restraints, but not 'moorDynR2_pt' or 'moorDynR2_bd' restraint options.
3. I have tried running it in both serial and parallel (4 cores), I get the same errors.
4. I have not modified any tutorial files (i.e. mooring input or dynamicMeshDict). Hence why I am quite confused as to why I'm getting these errors.

For completion, I'm attaching the mooring input file (from the tutorial) here:

lines_v2.txt:

--------------------- MoorDyn v2 Input File ------------------------------------
MoorDyn input file for a moored floating box EsfLOWC (Wu et al., 2018)
----------------------- LINE TYPES ------------------------------------------
TypeName   Diam     Mass/m     EA         BA/-zeta    EI         Cd     Ca     CdAx    CaAx
(name)     (m)      (kg/m)     (N)        (N-s/-)     (N-m^2)    (-)    (-)    (-)     (-)
main       0.003656 0.0607     29         -0.8        0          1.6    1.0    0.05    0.0
---------------------------- BODIES -----------------------------------------------------
ID   Attachment  X0     Y0    Z0    r0      p0     y0     Mass  CG*  I*       Volume  CdA*   Ca
(#)  (-)         (m)    (m)   (m)   (deg)   (deg)  (deg)  (kg)  (m)  (kg-m^2) (m^3)   (m^2)  (-)
1    Coupled     0      0     0     0       0      0      0     0    0        0       0      0
---------------------- CONNECTION PROPERTIES --------------------------------
ID    Type      X       Y       Z       Mass   Volume  CdA    Ca
(#)   (-)       (m)     (m)     (m)     (kg)   (mˆ3)   (m^2)  (-)
1     Fixed     -1.385   0.423  -0.5    0      0       0      0
1     Fixed     -1.385  -0.423  -0.5    0      0       0      0
1     Fixed      1.385   0.423  -0.5    0      0       0      0
1     Fixed      1.385  -0.423  -0.5    0      0       0      0
4     Body1     -0.1     0.1    -0.061  0      0       0      0
4     Body1     -0.1    -0.1    -0.061  0      0       0      0
4     Body1      0.1     0.1    -0.061  0      0       0      0
4     Body1      0.1    -0.1    -0.061  0      0       0      0
---------------------- LINES ----------------------------------------
ID   LineType   AttachA  AttachB  UnstrLen  NumSegs  LineOutputs
(#)   (name)     (#)      (#)       (m)       (-)     (-)
1     main       1        5         1.455     40      p
2     main       2        6         1.455     40      -
3     main       3        7         1.455     40      p
4     main       4        8         1.455     40      -
---------------------- OPTIONS -----------------------------------------
0.00002       dtM           time step to use in mooring integration (s)
3.0e6         kBot          bottom stiffness (Pa/m)
3.0e5         cBot          bottom damping (Pa-s/m)
0.5           WtrDpth       water depth (m)
1.0           dtIC          time interval for analyzing convergence during IC gen (s)
5.0           TmaxIC        max time for ic gen (s)
5.0           CdScaleIC     factor by which to scale drag coefficients during dynamic relaxation (-)
0.001         threshIC      threshold for IC convergence (-)
-------------------------- OUTPUTS --------------------------------
FairTen1 FairTen3
AnchTen1 AnchTen3
L1N1pZ L3N1pZ
------------------------- need this line --------------------------------------

dynamicMeshDict:

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

dynamicFvMesh       dynamicOversetFvMesh;

motionSolverLibs    (rigidBodyMeshMotion);

motionSolver       rigidBodyMotion;

rigidBodyMotionCoeffs
{
    report          on;

    solver
    {
        type Newmark;
        //gamma   0.1;    // Velocity integration coefficient
        //beta    0.1;    // Position integration coefficient
    }

    accelerationRelaxation 0.8;
    
    bodies
    {
        floatingObject
        {
            type            rigidBody;
            parent          root;
            mass            3.148;
            
            centreOfMass    (0 0 0);
            inertia         (0.015 0 0 0.015 0 0.021);
            // Transformation tensor and centre of rotation
            transform       (1 0 0 0 1 0 0 0 1) (0 0 -0.015);
            
            joint
            {
                type            composite;
                joints
                (
                    {
                        type Px;
                    }
                    {
                        type Pz;
                    }
                    {
                        type Ry;
                    }
                );
            }
            
            patches         (floatingObject);
            innerDistance   100.0;
            outerDistance   101.0;
        }
    }
    
    restraints
    {
		/*
        moodyR
        {
            type               moodyR;
            body               floatingObject;
            
            inputFile         "Mooring/boxWu_exPoint.m";
            refAttachmentPt
            (
                (-0.1      0.1    -0.062)
                (-0.1     -0.1    -0.062)
                ( 0.1      0.1    -0.062)
                ( 0.1     -0.1    -0.062)    
            );
        }
		*/
        
        /* Example mooring restraints as defined in libsixDoFMooring
           map3R || moodyR || moorDynR1 || moorDynR2
           
        moorDynR1
        {
            type               moorDynR1;
            body               floatingObject;
            refAttachmentPt
            (
                (-0.1      0.1    -0.062)
                (-0.1     -0.1    -0.062)
                ( 0.1      0.1    -0.062)
                ( 0.1     -0.1    -0.062)
            );
        }
		*/
        
		/*
        moorDynR2_pt
        {
            type               moorDynR2;
            body               floatingObject;
            couplingMode       "POINT"; // optional
            inputFile          "Mooring/lines_v2_point.txt";
            refAttachmentPt
            (
                (-0.1      0.1    -0.062)
                (-0.1     -0.1    -0.062)
                ( 0.1      0.1    -0.062)
                ( 0.1     -0.1    -0.062)  
            );
            writeMooringVTK    true;
            vtkPrefix          "mdv2_pt";
            vtkStartTime       0;
            outerCorrector     3;
        }
		*/

		// /*
        moorDynR2_bd
        {
            type               moorDynR2;
            body               floatingObject;
            couplingMode       "BODY"; // optional
            inputFile          "Mooring/lines_v2.txt";
            bodies
            (
                floatingObject
            );
            writeMooringVTK    true;
            vtkPrefix         "mdv2_bd";
            vtkStartTime       0;
            outerCorrector     3;
        }
		// */

/*
        map3R
        {
            type                map3R;
            body                floatingObject;
            inputFile           "Mooring/esflOWC_4lines.map";
            summaryFile         "Mooring/esflOWC_summary.map";
            writeMooringForces  true;
            outputFile          "Mooring/mapForces.map";
            waterDepth          0.5;

            refAttachmentPt
            (
                (-0.1      0.1    -0.062)
                (-0.1     -0.1    -0.062)
                ( 0.1      0.1    -0.062)
                ( 0.1     -0.1    -0.062)
            );
            
            writeMooringVTK        true;
            nNodes                 10;
            //nodesPerLine           (10 10 10 10);
        }
        */
    }

}

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

Please let me know if any further information would help, or if you have any ideas on the potential cause of this!

Thanks for your help!

[Desire71]

Hi,

I have faced the same error using moordyn. Have you reached any solution for it? I would appreciate it if you could share it if it had worked for you.