CEL and fortran for multiphase flow
 

Hi All,

I am trying to use MUSIG model to simulate air-water bubbly flow. CFX is only featured with "Prince and Blanch" model for the calculation of coalescence rate. But it is not applicable in my case. I have to write either a CEL or a fortran routine for the suitable model.

In order to get familiar with CEL and fortran, I use "Prince and Blanch" model for practice. When I use CEL to calculate the coalescence rate in "Prince and Blanch" model, the cfx works well. However, with same equations, when I use fortran routine to calculate the coalescence rate, error occurs:

ERROR #001100279 has occurred in subroutine ErrAction. |
| Message: |
| Floating point exception: Zero divide |
| |
| |
| |
| |
| |
+--------------------------------------------------------------------+

+--------------------------------------------------------------------+
| ERROR #001100279 has occurred in subroutine ErrAction. |
| Message: |
| Stopped in routine FPX: C_FPX_HANDLER |
| |
| |
| |
| |
| |
+--------------------------------------------------------------------+

+--------------------------------------------------------------------+
| An error has occurred in cfx5solve: |
| |
| The ANSYS CFX solver exited with return code 1. No results file |
| has been created. |
+--------------------------------------------------------------------+

Does anyone know what the problem is? It doesn't make sense for me that why CEL and fortran routine have different reults with same equations.

Any help or information is highly appreciated.
Thanks in advance