CFX 4.4 New free surface option

Viatcheslav Anissimov

Hi

I have tried to use new option "FREE SURFACE GRID MOVEMENT" of CFX 4.4. The result obtained is different from result with using symmetry plain in place of free surface. Problem is that I could not see the greed change of free surface in postprocessor (Analyse). Is it normal or not?

Any suggestions?

Thanks in advice.

Viatcheslav.

Citrate from command file:

>>CFX4
:>OPTIONS

FREE SURFACE GRID MOVEMENT
:>USER FORTRAN

USRBCS
:>MODEL TOPOLOGY
:>MODIFY PATCH

OLD PATCH NAME 'USER3D_SURFACE_REGION'

NEW PATCH NAME 'SURFACE REGION'
:>MODIFY PATCH

OLD PATCH NAME 'USER2D_FREE_SURFACE'

NEW PATCH NAME 'FREE SURFACE' /*

NEW PATCH TYPE 'SYMMET' */
:>MODIFY PATCH

OLD PATCH NAME 'USER2D_BASE_SURFACE'

NEW PATCH NAME 'BASE SURFACE'
:>MODEL DATA
:>PHYSICAL PROPERTIES

>>BUOYANCY PARAMETERS

GRAVITY VECTOR 0. 0. -9.81

THERMAL EXPANSION COEFFICIENT 1.216E-4

BUOYANCY REFERENCE TEMPERATURE 350.
:>SOLVER DATA
:>PROGRAM CONTROL PRESSURE REFERENCE POSITION 0.000001 2.39 -2.300001 /* AT FREE SURFACE */
:>FALSE TIMESTEPS

U VELOCITY 0.05

V VELOCITY 0.05

W VELOCITY 0.05
:>CREATE GRID
:>FREE SURFACE GRID PARAMETERS

CONVECTION CONTROLLED

FALSE TIME STEP 0.05

Z VERTICAL
:>MODEL BOUNDARY CONDITIONS
:>INLET BOUNDARIES

PATCH NAME 'USER2D_RISER_OUTLET'

NORMAL VELOCITY 0.879221

Citrate from usrbcs.f:

C-- INTERROGATE GETVAR FOR VARIABLE NUMBERS. C

CALL GETVAR('USRBCS','P ',IPRESS) C C SET IPHS = 1 FOR SINGLE PHASE FLOW. C

IPHS = 1 C C USE IPALL TO FIND 1D ADDRESSES OF ALL CELL CENTRES AT FREE SURFACE

CALL IPALL('FREE SURFACE','*','PATCH','CENTRES',

# IPT,NPT,CWORK,IWORK) C LOOP OVER ALL INTERIOR CELLS

DO 110 I=1,NPT C USE ARRAY IPT TO GET ADDRESS

INODE=IPT(I) C SET VARBCS

VARBCS(IPRESS,IPHS,INODE) = 0.0 110 CONTINUE