[ehsan_am86]

Hi everyone,

In fact, I have recently developed my solver (Add a class to original twoPhaseEulerFoam (populationBalance class) in OF 2.3.0). My solver works in serial mode works properly but when I want to use that with mpirun or mpiexec it is crashed. Any idea or comments for debugging, even short ones will be pretty much useful and appreciated. Thanks in advance.

P.S: This is short version of log file (complete one is attached):

[3]  in "/home/ehsanam/OpenFOAM/ehsanam-2.3.0/platforms/linux64GccDPOpt/bin/twoPhaseEulerCMFoam"
[cp0233:26977] *** Process received signal ***
[cp0233:26977] Signal: Floating point exception (8)
[cp0233:26977] Signal code:  (-6)
[cp0233:26977] Failing at address: 0x2e5f5200006961
[cp0233:26977] [ 0] /lib64/libc.so.6() [0x3565a326a0]
[cp0233:26977] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x3565a32625]
[cp0233:26977] [ 2] /lib64/libc.so.6() [0x3565a326a0]
[cp0233:26977] [ 3] /home/ehsanam/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam5PBiCG5solveERNS_5FieldIdEERKS2_h+0xaa8) [0x2b309cb346a8]
[cp0233:26977] [ 4] /home/ehsanam/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x11f) [0x2b309b71be2f]
[cp0233:26977] [ 5] twoPhaseEulerCMFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0xf7) [0x45c8e7]
[cp0233:26977] [ 6] twoPhaseEulerCMFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x80) [0x45cb30]
[cp0233:26977] [ 7] /home/ehsanam/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libcompressibleTwoPhaseSystem3.so(_ZN4Foam11classMethod7correctEv+0x123b) [0x2b309a7f787b]
[cp0233:26977] [ 8] twoPhaseEulerCMFoam() [0x4881ce]
[cp0233:26977] [ 9] /lib64/libc.so.6(__libc_start_main+0xfd) [0x3565a1ed5d]
[cp0233:26977] [10] twoPhaseEulerCMFoam() [0x4349a9]
[cp0233:26977] *** End of error message ***
[20]  in "/home/ehsanam/OpenFOAM/ehsanam-2.3.0/platforms/linux64GccDPOpt/bin/twoPhaseEulerCMFoam"
[cp0233:26994] *** Process received signal ***
[cp0233:26994] Signal: Floating point exception (8)
[cp0233:26994] Signal code:  (-6)
[cp0233:26994] Failing at address: 0x2e5f5200006972
[cp0233:26994] [ 0] /lib64/libc.so.6() [0x3565a326a0]
[cp0233:26994] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x3565a32625]
[cp0233:26994] [ 2] /lib64/libc.so.6() [0x3565a326a0]
[cp0233:26994] [ 3] /home/ehsanam/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam5PBiCG5solveERNS_5FieldIdEERKS2_h+0xaa8) [0x2ab939a6d6a8]
[cp0233:26994] [ 4] /home/ehsanam/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x11f) [0x2ab938654e2f]
[cp0233:26994] [ 5] twoPhaseEulerCMFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0xf7) [0x45c8e7]
[cp0233:26994] [ 6] twoPhaseEulerCMFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x80) [0x45cb30]
[cp0233:26994] [ 7] /home/ehsanam/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libcompressibleTwoPhaseSystem3.so(_ZN4Foam11classMethod7correctEv+0x123b) [0x2ab93773087b]
[cp0233:26994] [ 8] twoPhaseEulerCMFoam() [0x4881ce]
[cp0233:26994] [ 9] /lib64/libc.so.6(__libc_start_main+0xfd) [0x3565a1ed5d]
[cp0233:26994] [10] twoPhaseEulerCMFoam() [0x4349a9]
[cp0233:26994] *** End of error message ***
--------------------------------------------------------------------------
mpiexec noticed that process rank 5 with PID 26979 on node cp0233 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
[cp0233:26973] 23 more processes have sent help message help-mpi-runtime.txt / mpi_init:warn-fork
[cp0233:26973] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

[ehsan_am86]

I enhanced twoPhaseEulerFoam solver by including fi-equations to diameter modle (population balance with class method ). In current case, I have 7 classes and consequently, I have 7 equations for fi. Each fi vary between 0 and 1.
In serial mode, I did not have any problem but when I wanted to launch my solver in parallel mode, I faced an error during fi calculations named “floating point exception:”


As following list indicates, I made a lot of efforts to solve it:

1. Increase mesh number
2. Remove layers in all boundaries
3. Increase layer for inlet boundary
4. Use adjustTimeStep
5. Change fi-boundary conditions from 0 to 10e-5
6. Change solver of fi-equations from PBiCG to smoothSolver
7. Remove preconditioner in PBiCG solver for fi-equations
8. Change relaxation factor from 1 to 0.7 and 0.3

9.Change Scheme for div(phia,fi):

1. Gauss vanLeer to Gauss linear
2. Gauss Gauss linear to Gauss upwind (it is solved) ?

Now, I am wondering why upwind scheme solved the bug.