compressible two phase flow in CFX4.4
 

following is error message and my command file.

CFX-4 | Solver | forrtl: error (73): floating divide by zero CFX-4 | Solver | Image PC Routine Line Source CFX-4 | Solver | psolvempiNT.exe 004BED98 Unknown Unknown Unknown CFX-4 | Solver | psolvempiNT.exe 0044CC18 Unknown Unknown Unknown CFX-4 | Solver | psolvempiNT.exe 0041DEE6 Unknown Unknown Unknown CFX-4 | Solver | psolvempiNT.exe 004186AC Unknown Unknown Unknown CFX-4 | Solver | psolvempiNT.exe 0040B73B Unknown Unknown Unknown CFX-4 | Solver | psolvempiNT.exe 0040541C Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | Unknown 02290190 Unknown Unknown Unknown CFX-4 | Solver | runsolve4: Job completed CFX-4 | Solver | Output file : /tmp/m07.fo CFX-4 | Solver | CFX-4 | Solver | Finished.

>>CFX4
:>OPTIONS

TURBULENT FLOW

TWO DIMENSIONS

HEAT TRANSFER

COMPRESSIBLE FLOW

LAMINAR FLOW

TRANSIENT FLOW

NUMBER OF PHASES 2
:>MODEL TOPOLOGY
:>CREATE BLOCK

BLOCK NAME 'BLOCK-NUMBER-1'

BLOCK DIMENSIONS 40 20
:>CREATE PATCH

PATCH NAME 'INLET'

BLOCK NAME 'BLOCK-NUMBER-1'

PATCH TYPE 'INLET'

PATCH LOCATION 29 31 20 20 1 1

HIGH J
:>CREATE PATCH

PATCH NAME 'PRESS'

BLOCK NAME 'BLOCK-NUMBER-1'

PATCH TYPE 'PRESSURE BOUNDARY'

PATCH LOCATION 9 11 1 1 1 1

LOW J
:>MODEL DATA
:>TITLE

PROBLEM TITLE 'test'
:>AMBIENT VARIABLES

TEMPERATURE 2.9100E+02
:>PHYSICAL PROPERTIES

>>COMPRESSIBILITY PARAMETERS

PHASE NAME 'PHASE2'

INCOMPRESSIBLE

>>COMPRESSIBILITY PARAMETERS

PHASE NAME 'PHASE1'

WEAKLY COMPRESSIBLE

UNIVERSAL GAS CONSTANT 8.3140E+03

FLUID MOLECULAR WEIGHT 2.8950E+01

REFERENCE PRESSURE 1.0130E+05

>>FLUID PARAMETERS

PHASE NAME 'PHASE2'

VISCOSITY 1.0000E-05

DENSITY 38.0

>>STANDARD FLUID

PHASE NAME 'PHASE1'

FLUID 'AIR'

STANDARD FLUID REFERENCE TEMPERATURE 2.9100E+02

>>HEAT TRANSFER PARAMETERS

PHASE NAME 'PHASE2'

THERMAL CONDUCTIVITY 2.3400E-04

FLUID SPECIFIC HEAT 1* 0.000000E+00

ENTHALPY REFERENCE TEMPERATURE 2.8800E+02

>>MULTIPHASE PARAMETERS

>>PHASE DESCRIPTION

PHASE NAME 'PHASE1'

GAS

CONTINUOUS

MINIMUM VOLUME FRACTION 1.0000E-05

MODIFY EMPTY CELL VELOCITY 0.0000E+00

>>PHASE DESCRIPTION

PHASE NAME 'PHASE2'

SOLID

DISPERSE

MEAN DIAMETER 1.0000E-03

MINIMUM VOLUME FRACTION 1.0000E-05

SOLID PRESSURE

RESTITUTION COEFFICIENT 6.0000E-01

COMPACTION VOLUME FRACTION 6.4000E-01

MODIFY EMPTY CELL VELOCITY 0.0000E+00

>>MULTIPHASE MODELS

>>MOMENTUM

INTER PHASE TRANSFER

SINCE

>>INTER PHASE TRANSFER MODELS

>>MOMENTUM

FIRST PHASE NAME 'PHASE1'

SECOND PHASE NAME 'PHASE2'

>>PARTICLE DRAG MODEL

FLOW REGIME 'ERGUN'

>>TRANSIENT PARAMETERS

>>FIXED TIME STEPPING

TIME STEPS 100* 1.000000E-03

INITIAL TIME 0.0000E+00
:>SOLVER DATA
:>PROGRAM CONTROL

MAXIMUM NUMBER OF ITERATIONS 10

MASS SOURCE TOLERANCE 1.0000E-06
:>CREATE GRID
:>SIMPLE GRID

BLOCK NAME 'BLOCK-NUMBER-1'

DX 40* 3.000000E-03

DY 20* 3.000000E-03

DZ 3.000000E-03

X START 0.0000E+00

Y START 0.0000E+00

Z START 0.0000E+00
:>MODEL BOUNDARY CONDITIONS
:>INLET BOUNDARIES

PHASE NAME 'PHASE1'

PATCH NAME 'INLET'

NORMAL VELOCITY 3.0000E+00

VOLUME FRACTION 9.0000E-01

TEMPERATURE 2.9300E+02
:>INLET BOUNDARIES

PHASE NAME 'PHASE2'

PATCH NAME 'INLET'

NORMAL VELOCITY 2.0000E+00

VOLUME FRACTION 1.0000E-01

TEMPERATURE 2.9300E+02
:>PRESSURE BOUNDARIES

PHASE NAME 'PHASE2'

PATCH NAME 'PRESS'

PRESSURE 0.0000E+00

VOLUME FRACTION 0.0000E+00
:>PRESSURE BOUNDARIES

PHASE NAME 'PHASE1'

PATCH NAME 'PRESS'

PRESSURE 0.0000E+00

VOLUME FRACTION 1.0000E+00
:>STOP