|
[Sponsors] |
April 3, 2004, 05:50 |
comfiling problem ..
|
#1 |
Guest
Posts: n/a
|
Hi all I set the density based on USER-f(P.T) in morecular property,and i made simple fortran code below. but densit.f doesn't comfiled so that the output of density shows constant value wich i gave as a initial value. what is problem in this case.
regards mspark fortran code SUBROUTINE DENSIT(DEN,DENDP) INCLUDE 'comdb.inc' COMMON/USR001/INTFLG(100) INCLUDE 'usrdat.inc' DIMENSION SCALAR(50) EQUIVALENCE( UDAT12(001), ICTID ) EQUIVALENCE( UDAT11(001), CP ) EQUIVALENCE( UDAT11(005), PR ) EQUIVALENCE( UDAT11(006), P ) EQUIVALENCE( UDAT11(007), T ) EQUIVALENCE( UDAT11(009), SCALAR(01) ) EQUIVALENCE( UDAT11(059), U ) EQUIVALENCE( UDAT11(060), V ) EQUIVALENCE( UDAT11(061), W ) EQUIVALENCE( UDAT11(067), X ) EQUIVALENCE( UDAT11(068), Y ) EQUIVALENCE( UDAT11(069), Z ) DENSITYL1=96693.*T**(-0.7684) DENSITYV1=4.*10.**(-22.)*T**9.2727 DENSITYL=3589.7*(P)**-0.081+0.001/DENSITYL1 DENSITYV=7E-05*(P)**0.9772+0.001/DENSITYV1 SV_F=1/DENSITYL SV_G=1/DENSITYV SV_FG=abs(SV_G-SV_F) REALH=CP*T VT=72.424*(P)**0.1054 HFLUID=7.2051*VT**1.8239 HLATENT=20000.*10000*VT**(-1.2023) QUALITY=(REALH-HLATENT)/HFLUID DEN=1/(SV_F+QUALITY*SV_FG) RETURN END |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
VOF diffuse problem on two fluids problem | Fang Jin | FLUENT | 5 | February 17, 2011 05:38 |
Problem in Modelling Heat Transfer Problem | Deepak R | FLUENT | 1 | December 6, 2007 09:37 |
Problem in cavity flow problem | saad | Main CFD Forum | 4 | November 1, 2007 07:45 |
natural convection problem for a CHT problem | Se-Hee | CFX | 2 | June 10, 2007 06:29 |
problem in solving "wave generation" problem | san | FLUENT | 2 | April 3, 2006 23:37 |