unexpectedly stopping of simulation
 

i am running a simulation of interfoam solver it would have 30 seconds time but dont know the reason why its getting stopped at 2 seconds.


order and list of commands is included in this file


file: https://drive.google.com/file/d/190z...usp=drive_link



Error:
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 1.80784e-09, Final residual = 1.80784e-09, No Iterations 0
Phase-1 volume fraction = 0.00489958 Min(alpha.water) = -0.000396251 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.00489958 Min(alpha.water) = -0.000396251 Max(alpha.water) = 1
[1] [stack trace]
[1] =============
[1] #1 Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
[1] #2 ? in /lib/x86_64-linux-gnu/libc.so.6
[1] #3 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
[1] #4 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
[1] #5 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
[1] #6 Foam::fv::gaussLaplacianScheme<Foam::Vector<double >, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
[1] #7 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
[1] #8 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so
[1] #9 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so
[1] #10 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/bin/interFoam
[1] #11 ? in /lib/x86_64-linux-gnu/libc.so.6
[1] #12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
[1] #13 ? in /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/bin/interFoam
[1] =============
[sahil-Lenovo-ideapad-520-15IKB:05750] *** Process received signal ***
[sahil-Lenovo-ideapad-520-15IKB:05750] Signal: Floating point exception (8)
[sahil-Lenovo-ideapad-520-15IKB:05750] Signal code: (-6)
[sahil-Lenovo-ideapad-520-15IKB:05750] Failing at address: 0x3e800001676
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7f5e6f442520]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 1] /lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c)[0x7f5e6f4969fc]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 2] /lib/x86_64-linux-gnu/libc.so.6(raise+0x16)[0x7f5e6f442476]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 3] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7f5e6f442520]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 4] /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(+0x49cd20)[0x7f5e71a9cd20]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 5] /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(+0xb31a90)[0x7f5e72131a90]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 6] /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(+0xb7adc6)[0x7f5e7217adc6]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 7] /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam2fv20gaussLaplacianSche meINS_6VectorIdEEdE12fvmLaplacianERKNS_14Geometric FieldIdNS_13fvsPatchFieldENS_11surfaceMeshEEERKNS5 _IS3_NS_12fvPatchFieldENS_7volMeshEEE+0x227)[0x7f5e7217b127]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 8] /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(+0xb81f65)[0x7f5e72181f65]
[sahil-Lenovo-ideapad-520-15IKB:05750] [ 9] /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(+0x1101db)[0x7f5e707a81db]
[sahil-Lenovo-ideapad-520-15IKB:05750] [10] /usr/lib/openfoam/openfoam2306/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(+0xb87ec)[0x7f5e707507ec]
[sahil-Lenovo-ideapad-520-15IKB:05750] [11] interFoam(+0x304b0)[0x558fc2e324b0]
[sahil-Lenovo-ideapad-520-15IKB:05750] [12] /lib/x86_64-linux-gnu/libc.so.6(+0x29d90)[0x7f5e6f429d90]
[sahil-Lenovo-ideapad-520-15IKB:05750] [13] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0x80)[0x7f5e6f429e40]
[sahil-Lenovo-ideapad-520-15IKB:05750] [14] interFoam(+0x3f255)[0x558fc2e41255]
[sahil-Lenovo-ideapad-520-15IKB:05750] *** 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 sahil-Lenovo-ideapad-520-15IKB exited on signal 8 (Floating point exception).






please let me know how to correct it i ran similar program (only cad files are changed) Earlier that ran smoothly

It seems that your simulation was crashed. Can you show the solutions just before it stopped?

Hello,
You seem to have a floating point exception error. So you divide by zero somewhere when the solver runs. The first step would be to check your case properly and also boundary conditions.

Thanks

i tried but am unable to solve it can u help i did shared my case link

https://drive.google.com/drive/folde...98?usp=sharing


here you can see the last files generated


please guide further

Hello,
I do not have a straight answer as to why your case crashes.

However, after having taken look at your case, I can suggest the following
1) Decrease the maximum Courant number to 0.5 (even 0.1 to just check)
2) Increase number of outercorr to 3 and check (increase further if needed)
3) Increase number of nAlphaSubCycles
4) If none of the above work, then try to lower your inlet velocity and then see if the crash happens. If it crashes, try steps 1-3 for the low velocity case. If it doesn't crash for low velocity cases then we can try to suggest some modifications in fvSchemes which can be the problem for this issue.

I don't know what exact problem you are simulating, but be careful with the outlet, try to keep it far away from regions of high gradients if possible.
I see that you use inletOutlet or pressureInletOutletVelocity at outlet which is good.

Try to visualize the flow field in Paraview just before it crashes. Then check if there is something weird (anything non-physical) going on.
I would suspect the pressure/velocity field, as your solver crashed after computing phase fractions and before the pressure equation.

Thanks