CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

CFX 4.4 New free surface option

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 3, 2002, 06:27
Default CFX 4.4 New free surface option
  #1
Viatcheslav Anissimov
Guest
 
Posts: n/a
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
  Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
emag beta feature: charge density charlotte CFX 4 March 22, 2011 09:14
transient simulation of a rotating rectangle icesniffer CFX 1 August 8, 2009 07:25
No results for solid domain Gary Holland CFX 10 March 13, 2009 03:30
error using combination of step function xujjun CFX 1 January 15, 2008 16:46
Two-Phase Buoyant Flow Issue Miguel Baritto CFX 4 August 31, 2006 12:02


All times are GMT -4. The time now is 05:32.