Error in TFM simulation using Cartesian grid with quadrics

Hello, I am trying to simulate a 2D fluidized bed combustor.

I used quadrics to make the combustor geomerty, and it looks there is no problem in the geometry.

But just after running, the simulation immediately stops with the following error:


Welcome to MFiX - https://mfix.netl.doe.gov
MFiX-GUI version 21.1.4
Ready
Loading cfb_combustor_2d_air_200_test03.mfx from /home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test03
Solver: MFiX-TFM
Reading mesh stats: /home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test03/MESH_STATS.LOG
Loaded cfb_combustor_2d_air_200_test03.mfx
Previous MFiX run is resumable. Reset job to edit model
Ready
Exporting cfb_combustor_2d_air_200_test04
Export success: Successfully exported to:
/home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test04/cfb_combustor_2d_air_200_test04.mfx
Open exported project?
Starting env OMP_NUM_THREADS=2 /home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test03/mfixsolver -s -f /home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test03/cfb_combustor_2d_air_200_test03.mfx
MFiX running: elapsed time 0:00:00
Number of SMP threads: 2
MFiX (21.1.4 ) simulation on host: localhost.localdomai
Run name: CFB_COMBUSTOR_2D_AIR Time: 14:17 Date: 3-31-2022
Project version: 178
Memory required: 9.00 Mb

INFO open_files.f:374
The following mesh file was not found: CFB_COMBUSTOR_2D_AIR_200_TEST03.msh
Mesh generation will be performed and a mesh file will be saved.

The following mesh file was successfully opened: CFB_COMBUSTOR_2D_AIR_200_TEST03.msh
INFO FOR QUADRIC 2 Defining Cone for Cylinder to Cylinder junction Between Quadrics 1 AND 3
QUADRIC: 2 WAS DEFINED AS Y_CONE WITH AN HALF-ANGLE OF 17.650124219930124 DEG.
INFO FOR QUADRIC 4 Defining Cone for Cylinder to Cylinder junction Between Quadrics 3 AND 5
QUADRIC: 4 WAS DEFINED AS Y_CONE WITH AN HALF-ANGLE OF -12.804266065286747 DEG.
INFO FOR QUADRIC 6 Defining Cone for Cylinder to Cylinder junction Between Quadrics 5 AND 7
QUADRIC: 6 WAS DEFINED AS Y_CONE WITH AN HALF-ANGLE OF 16.260204708311946 DEG.
INFO FOR QUADRIC 8 Defining Cone for Cylinder to Cylinder junction Between Quadrics 7 AND 9
QUADRIC: 8 WAS DEFINED AS Y_CONE WITH AN HALF-ANGLE OF -14.036243467926468 DEG.
MFIX WITH CARTESIAN GRID IMPLEMENTATION.
RE-INDEXING IS TURNED OFF.
INITIALIZING VELOCITY NODES…
INTERSECTING GEOMETRY WITH SCALAR CELLS…
INTERSECTING GEOMETRY WITH U-MOMENTUM CELLS…
INTERSECTING GEOMETRY WITH V-MOMENTUM CELLS…
SETTING CUT CELL TREATMENT FLAGS…
COMPUTING INTERPOLATION FACTORS IN U-MOMENTUM CELLS…
COMPUTING INTERPOLATION FACTORS IN V-MOMENTUM CELLS…
COMPUTING 1/DX, 1/DY FOR U-MOMENTUM CELLS…
COMPUTING 1/DX, 1/DY FOR V-MOMENTUM CELLS…
FINDING MASTER CELLS FOR U-MOMENTUM CELLS…
FINDING MASTER CELLS FOR V-MOMENTUM CELLS…
MESH STATISTICS:
NUMBER OF CELLS = 17920
NUMBER OF STANDARD CELLS = 12516 ( 69.84 % of Total)
NUMBER OF CUT CELLS = 1140 ( 6.36 % of Total)
NUMBER OF FLUID CELLS = 13656 ( 76.21 % of Total)
NUMBER OF BLOCKED CELLS = 4264 ( 23.79 % of Total)
CARTESIAN GRID PRE-PROCESSING COMPLETED IN 3.38 SECONDS.
WRITING VTU FILE : CFB_COMBUSTOR_2D_AIR_200_TEST03_MESH.vtu .
… DONE.
Saving Mesh connectivity and flags…
Setting master from flags…
Gather and Write MasterList…
Saving Mesh info for Scalar cells…
Saving Mesh info for U-MOMENTUM cells…
Saving Mesh info for V-MOMENTUM cells…
WROTE BOUNDARY IN VTK FILE : CFB_COMBUSTOR_2D_AIR_200_TEST03_boundary.vtk
PRE_PROCESSING COMPLETE.
next_reca = 513
0 1
1 0
next_reca = 514

INFO set_bc0.f:549
Compressible: IJK_P_g remaining undefined.

t= 0.000000 Wrote
SPx: 1
,2
,3,4,5,6
,7
,8
,9
;
.RES;
Wall Time - Elapsed: 0.12s Est. Remaining: 0.00s
WRITING VTU FILE : BACKGROUND_IC_0000.vtu .

.
.
DONE.
Time = 0.0000 Dt = 0.10000E-03

Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:

Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:
** Something went wrong while running addr2line. **
** Falling back to a simpler backtrace scheme. **
#26 0x7f9132954545
#27 0x7f9132954cde
#28 0x7f91855483ff
#29 0x7f91330b5133
#30 0x7f91330b5fbc
#31 0x7f91330b74be
#32 0x7f9132f42d3e
#33 0x7f9132f3da68
#34 0x7f9132d1efe7
#35 0x7f9132ecffbb
#36 0x7f9132d21644
#37 0x7f9132d14df6
#38 0x55cbeecb6da7
#39 0x55cbeed156f7
#40 0x55cbeece3767
#41 0x55cbeec91253
#42 0x55cbeece3767
#43 0x55cbeec91253
#44 0x55cbeece3767
#45 0x55cbeecee61b
#46 0x55cbeecb7edd
#47 0x55cbeed456a9
#48 0x55cbeed2da46
#49 0x7f91858e7ea4
#50 0x7f91856108dc
#51 0xffffffffffffffff
Reading mesh stats: /home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test03/MESH_STATS.LOG
Reading mesh boundary: /home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test03/CFB_COMBUSTOR_2D_AIR_200_TEST03_boundary.vtk
/home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test03/mfixsolver: line 3: 23031 Floating point exception(core dumped) env PYTHONPATH="/home/skchoi4091/MFiX/cfb_combustor_2d_air_200_test03":"":${PYTHONPATH:+:$PYTHONPATH} /home/skchoi4091/anaconda3/envs/mfix-21.1.4/bin/python3 -m mfixgui.pymfix “$@”
Previous MFiX run is resumable. Reset job to edit model
MFiX process has stopped
Error: Cmd(‘git’) failed due to: exit code(129)
cmdline: git tag --sort=creatordate
stderr: ‘error: unknown option `sort=creatordate’
usage: git tag [-a|-s|-u ] [-f] [-m |-F ] []
or: git tag -d …
or: git tag -l [-n[]] [–contains ] [–points-at ]
[…]
or: git tag -v …

-l, --list            list tag names
-n[<n>]               print <n> lines of each tag message
-d, --delete          delete tags
-v, --verify          verify tags

Tag creation options
-a, --annotate annotated tag, needs a message
-m, --message
tag message
-F, --file read message from file
-s, --sign annotated and GPG-signed tag
–cleanup how to strip spaces and #comments from message
-u, --local-user
use another key to sign the tag
-f, --force replace the tag if exists
–column[=] show tag list in columns

Tag listing options
–contains print only tags that contain the commit
–points-at print only tags of the object

File “…/gui.py”, line 613, in
btn.pressed.connect(lambda mode=mode: self.change_mode(mode))
File “…/gui.py”, line 1711, in change_mode
vc.update(self.ui.textbrowser_history)
File “…/vc.py”, line 57, in update
report_text = html % gen_report()
File “…/vc.py”, line 26, in gen_report
tags = repo.git.tag(sort=‘creatordate’).split(’\n’)
File “/home/skchoi4091/anaconda3/envs/mfix-21.1.4/lib/python3.8/site-packages/git/cmd.py”, line 585, in
return lambda *args, **kwargs: self._call_process(name, *args, **kwargs)
File “/home/skchoi4091/anaconda3/envs/mfix-21.1.4/lib/python3.8/site-packages/git/cmd.py”, line 1124, in _call_process
return self.execute(call, **exec_kwargs)
File “/home/skchoi4091/anaconda3/envs/mfix-21.1.4/lib/python3.8/site-packages/git/cmd.py”, line 928, in execute
raise GitCommandError(redacted_command, status, stderr_value, stdout_value)

The MFiX version information is as follows:

MFiX version: 21.1.4
Python version: 3.8.10 (default, Jun 4 2021, 15:09:15) [GCC 7.5.0]
Qt Wrapper: PyQt5
Qt Version: 5.14.2
qtpy Version: 1.9.0
Numpy Version: 1.19.3
Nodeworks Version: Unavailable
Flask Version: 1.1.4
Sphinx Version: Unavailable
psutil Version: 5.6.7
VTK Version: 9.0.3
OpenGL Backend: 2
System info: Linux release 3.10.0-1127.8.2.el7.x86_64 running on x86_64
Install location: /home/skchoi4091/anaconda3/envs/mfix-21.1.4/lib/python3.8/site-packages
Default Solver executable:None
Solver source code: /home/skchoi4091/anaconda3/envs/mfix-21.1.4/lib/python3.8/site-packages

Could you please advise me how to resolve this?

Thank you.

cfb_combustor_2d_air_200_test03.mfx (13.9 KB)

I am not able to reproduce. It runs fine for me on Linux with MFIX 21.1.4 . I tried 2 and 32 SMP threads.

Does it run in serial?

It does not run in serial either. Due to the network security of my institute, I was not able to install MFiX with conda command so I have installed with Pip. Could this error be related with the installation method? Other case with a simple geometry without quadrics runs fine.

I have updated the compiler from gcc 4.8.5 to 8.3.1, and now it runs fine.
Thank you!