I want to simulate the liquid steel flow in a Ladle stirred by Ar injeted from the bottem. But the result is wrong! The below is my fc file:

>> CFX4

#CALC

RHOL = 1000.; /* Properties */

RHOG = 1.28;

VISL = 1.0E-3;

VISG = 1.0E-5;

SRFTEN = 0.025;

BUBDIA = 2.0E-3;

U1 = 0.0; /* Inlet */

U2 = 0.0;

V1 = 0.0;

V2 = 0.0;

W1 = 0.0;

W2 = 0.0;

VF2 = 0.01;

VF1 = 1.0-VF2;

DIAM = 0.1;

#ENDCALC

>> OPTIONS

NUMBER OF PHASES 2

BUOYANT FLOW

TURBULENT FLOW

END

>> MODEL DATA

>> PHYSICAL PROPERTIES

>> BUOYANCY PARAMETERS

GRAVITY VECTOR 0.0 -9.8 0.

BUOYANCY REFERENCE DENSITY #RHOG

END

>> FLUID PARAMETERS

PHASE NAME 'PHASE1'

VISCOSITY #VISL

DENSITY #RHOL

END

>> FLUID PARAMETERS

PHASE NAME 'PHASE2'

VISCOSITY #VISG

DENSITY #RHOG

END

>> TURBULENCE PARAMETERS

>> TURBULENCE MODEL

PHASE NAME 'PHASE1'

TURBULENCE MODEL 'K-EPSILON'

PARTICLE INDUCED TURBULENCE 'NONE'

END

>> TURBULENCE MODEL

PHASE NAME 'PHASE2'

TURBULENCE MODEL 'K-EPSILON'

PARTICLE INDUCED TURBULENCE 'NONE'

END

>> TURBULENCE CONSTANTS

CMUPI 0.6

>> MULTIPHASE PARAMETERS

>> PHASE DESCRIPTION

PHASE NAME 'PHASE1'

LIQUID

CONTINUOUS

END

>> PHASE DESCRIPTION

PHASE NAME 'PHASE2'

GAS

DISPERSE

MEAN DIAMETER #BUBDIA

END

>> MULTIPHASE MODELS

>> MOMENTUM

INTER PHASE TRANSFER

IPSAC

END

>> INTER PHASE TRANSFER MODELS

>>MOMENTUM

FIRST PHASE NAME 'PHASE1'

SECOND PHASE NAME 'PHASE2'

SURFACE TENSION COEFFICIENT #SRFTEN

>>PARTICLE DRAG MODEL

FLOW REGIME 'DISTORTED'

END

>>NON DRAG FORCES

LIFT FORCE

LIFT COEFFICIENT 0.03

WALL LUBRICATION FORCE

LUBRICATION COEFFICIENTS -0.01 0.05

TURBULENT DISPERSION FORCE

DISPERSION COEFFICIENT 0.1

END

>> AMBIENT VARIABLES

PHASE NAME 'PHASE1'

U VELOCITY #U1

V VELOCITY #V1

W VELOCITY #W1

VOLUME FRACTION #VF1

END

>> AMBIENT VARIABLES

PHASE NAME 'PHASE2'

U VELOCITY 0.0

V VELOCITY 0.0

W VELOCITY 0.0

VOLUME FRACTION #VF2

END

>> SOLVER DATA

>> PROGRAM CONTROL

MAXIMUM NUMBER OF ITERATIONS 500

MASS SOURCE TOLERANCE 1.0E-10

END

>> EQUATION SOLVERS

PRESSURE 'AMG'

END

>> UNDER RELAXATION FACTORS

U VELOCITY 6.5000E-01

V VELOCITY 6.5000E-01

W VELOCITY 6.5000E-01

PRESSURE 1.0000E+00

VOLUME FRACTION 0.6000E+00

VISCOSITY 1.0000E+00

K 3.0000E-01

EPSILON 3.0000E-01

>> FALSE TIMESTEPS

U VELOCITY 1.0E-1

V VELOCITY 1.0E-1

W VELOCITY 1.0E-1

END

>> MODEL BOUNDARY CONDITIONS

>> INLET BOUNDARIES

PATCH NAME 'INLET'

PHASE NAME 'PHASE1'

U VELOCITY #U1

V VELOCITY #V1

VOLUME FRACTION #VF1

>> INLET BOUNDARIES

PATCH NAME 'INLET'

PHASE NAME 'PHASE2'

NORMAL VELOCITY 23.1

VOLUME FRACTION #VF2

>> WALL BOUNDARIES

PATCH NAME 'WALLB'

PHASE NAME 'PHASE1'

TAUX 0.0

TAUY 0.0

END

>> WALL BOUNDARIES

PATCH NAME 'WALLOUT'

PHASE NAME 'PHASE2'

TAUX 0.0

TAUY 0.0

END
:> STOP

What's the wrong?

Could you be more specific? Do you have a preliminary error or is the result unphysical?

At least replace :>STOP with >>STOP

Regs, Gert-Jan

Oh, sorry! the ":" was a mistake made by me when post this message.

My project is a Ladle(a cylinder) in which Liquid Steel is. AIR(or Ar)is injetted into the ladle from a little inlet on the bottom(the Velocity is about 20~30 m/s). And the top of the ladle(above liquid steel) is also air(quiescent).

Regs,