MODSWF.f

Peter

Hi,

I'd like to use the folowing subroutine. I know in theory what it does (heat transfer calculation in buyoncy driven flows). But I do not know, how it acts with the heat transfer calculation in StarCD. Do I get new htran or hcoef? Does the simulation use HWF for calculating the heat flux? And how can I postprocess values like HWF? How can I check, if the subroutine works or if it fails?

Thanks Peter

C************************************************* ************************

SUBROUTINE MODSWF(TBULK,SCBULK,HWF,SMTC) C Modify wall functions C************************************************* ************************ C--------------------------------------------------------------------------* C STAR RELEASE 3.150 * C--------------------------------------------------------------------------* ...

C------------------------------------------------------------------------- C C This subroutine enables the user to specfy heat and mass transfer C wall functions HWF,SMTC(IS) in an arbitrary manner. C C ** Parameters to be returned to STAR: HWF,SMTC(IS) C C------------------------------------------------------------------------- C C Sample coding: Wall heat transfer coefficient for the case of C natural convection C

LOGICAL VERT,HORZ CC---- EXPONENTS FOR MIXED CONVECTION

AN=3.0

THIRD=1./3. C CC---- VOLUMETRIC EXPANSION COEFFICIENT BETA C

GO TO (10,20,30),IDEN C CC---- IDEAL GAS LAW C 10 BETA=1./(0.5*(TW+TBULK(IMAT)))

GO TO 40 C CC---- ISOBARIC DENSITY LAW C 20 GO TO 40 C CC---- SPECIFY SOURCE CODING FOR BETA HERE C 30 CONTINUE 40 CONTINUE C

TDIF=ABS(TW-TBULK(IMAT))

VISK2=VISM**2/DEN**2 CC----GRASHOF NUMBER DIVIDED BY CHARACTERISTIC LENGTH OF THE SOLUTION C DOMAIN AND THERMAL CONDUCTIVITY

G=SQRT(G1**2+G2**2+G3**2)

GRASHF=G*BETA*TDIF*PRANL/VISK2

CONDUC=VISM/PRANL*CPREF C CC---- NATURAL CONVECTION HEAT TRANSFER COEF. C

ALFA=0.113*(GRASHF*PRANL)**THIRD*CONDUC

HWF=(HWF**AN+ALFA**AN)**(1./AN) C-------------------------------------------------------------------------C

RETURN

END C