[phelps]

Hello everyone! I’ve been learning OpenFoam for the past two weeks and this is my first post so I would like to apologize if I dont explain myself correctly or if I ask silly questions.

I'm trying to study the effect of the wake produced by one vertical axys wind turbine in another turbine next to it (see the attached image). So, I need to achieve a flow induced rotation in both turbines. I'm able to do this with only one turbine using pimpleFoam, AMI and sixDoFRigidBodyMotion (I’m using this Holzmann tutorial as a base case: https://holzmann-cfd.com/community/t...l-wind-turbine), but I dont know how to set up the dynamicMeshDict file for more than one body since sixDoFRigidBodyMotion library is only capable to manage one body (as far as I know).

I've been trying to find a way to do this in OpenFoam 8 for a week, but none of the solutions I have found work for this version, everything seems removed.

Last thing I have tried is rigidBodyMotion which seems to be able to handle multiple uncoupled bodies. Here is my dynamicMeshDict file:

 dynamicFvMesh   dynamicMotionSolverFvMesh;
 
 
 motionSolverLibs   ("librigidBodyMeshMotion.so");
 
 
 motionSolver    rigidBodyMotion;
 
 
 report          on;
 
 
 solver
 {
     type Newmark;
 }
 
 
 rho        rhoInf;
  
 rhoInf     1.2;
   
 accelerationRelaxation 0.4;
  
 bodies
 {
     VAWT
     {
         type            rigidBody;
         parent          root;
 
 
         centreOfMass    (0 0 0.025);
         mass             0.000461158*1500;
         inertia         (0.0053176 0.0053176 0.0103473 6.898e-6 3.54533e-6 3.545e-6);
         transform       (1 0 0 0 1 0 0 0 1) (0 0 0);
  
         joint
         {
                     type Rz;
         }
 
         patches         (VAWT);
         innerDistance   0.6;
         outerDistance   0.7;
     }
      
     VAWT2
     {
         type            rigidBody;
         parent          root;
 
         centreOfMass    (2.5 0 0.025);
         mass             0.000461158*1500;
         inertia         (0.0053176 0.0053176 0.0103473 6.898e-6 3.54533e-6 3.545e-6);
         transform       (1 0 0 0 1 0 0 0 1) (0 0 0);
  
         joint
         {
             type    Rz;
         }
  
         patches         (VAWT2);
         innerDistance   0.6;
         outerDistance   0.7;
     }
 }

The simulation started but after a few time steps, the solver stopped and I got this error:


[1] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
[1] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foamperator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[1] #5 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[1] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
[pedro-HP-Laptop-15s-fq1xxx:11240] *** Process received signal ***
[pedro-HP-Laptop-15s-fq1xxx:11240] Signal: Floating point exception (8)
[pedro-HP-Laptop-15s-fq1xxx:11240] Signal code: (-6)
[pedro-HP-Laptop-15s-fq1xxx:11240] Failing at address: 0x3e800002be8
[pedro-HP-Laptop-15s-fq1xxx:11240] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f23a951d210]
[pedro-HP-Laptop-15s-fq1xxx:11240] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f23a951d18b]
[pedro-HP-Laptop-15s-fq1xxx:11240] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x46210)[0x7f23a951d210]
[pedro-HP-Laptop-15s-fq1xxx:11240] [ 3] /opt/openfoam8/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRK NS_5UListIdEE+0xaf)[0x7f23aa13867f]
[pedro-HP-Laptop-15s-fq1xxx:11240] [ 4] pimpleFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7volMes hEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_11d imensionedIdEERKS8_+0x2f7)[0x556fec2d9c37]
[pedro-HP-Laptop-15s-fq1xxx:11240] [ 5] pimpleFoam(+0x2debd)[0x556fec27eebd]
[pedro-HP-Laptop-15s-fq1xxx:11240] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0x7f23a94fe0b3]
[pedro-HP-Laptop-15s-fq1xxx:11240] [ 7] pimpleFoam(+0x313fe)[0x556fec2823fe]
[pedro-HP-Laptop-15s-fq1xxx:11240] *** End of error message ***
--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 0 on node pedro-HP-Laptop-15s-fq1xxx exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------


I don’t have a clue on how to read this error, so I would be really grateful if someone can tell me what is happening here.

Also, please tell me if I’m using the correct solvers and libraries or if there is a better way to simulate multiple body motion with 6DoF.

Thank you.