TALK=T;RUN(1,1)
text(3D Square Cavity; Various Fluids
#cls
title
DISPLAY
This In-Form case illustrates the setting of properties by
selecting a fluid from those for which property formulae are
assembled in case 089. Their names are listed in this Q1,
First, case 277 is loaded. This is an old library case concerning
steady laminar flow and heat transfer in a cubic cavity with a
moving wall. The material properties were arbitrary round-number
constants.
Then H1, the variable representing enthalpy in 277, is renamed to
TEM1; and PIL settings are made for the wall temperatures.
STORE(ENT1,CP1,RHO1,CP1,ENUL,COND) creates storage for the
properties which In-Form will compute and use.
The relevant In-Form settings are:
1. fluid_name = air
which is one of the fluids in the list.
Note that all that is necessary in order to choose another
fluid is to move the relevant fluid_name setting one space to
the left.
2. (STORED var ENT1 is TEM1*CP1)
which tells PHOENICS how to the new enthalpy cariable.
3, L($089)
which loads case 089
The Q1 contains PHOTON USE commands
ENDDIS
Load($277)
PRESS0
expert=f
GROUP 1. Run title and other preliminaries
TEXT(3D SQUARE CAVITY; various fluids : 761
GROUP 7. Variables stored, solved & named
NAME(H1)=TEM1
STORE(ENT1,CP1,RHO1,ENUL,COND)
PRESS0=1.E5
GROUP 11. Initialization of variable or porosity fields
REAL(THOT,TCOLD); TCOLD=273; THOT=323
FIINIT(TEM1)=0.5*(TCOLD+THOT)
GROUP 13. Boundary conditions and special sources
** South wall, at rest
COVAL(SOUTH,TEM1,1.0,TCOLD)
** North wall, moving
COVAL(MOVING,TEM1,1.0,THOT)
** West wall, at rest
COVAL(WEST,TEM1,1.0,TCOLD)
** East wall, at rest
COVAL(EAST,TEM1,1.0,TCOLD)
Group 15. Terminate Sweeps
ISWC1=1
Group 17. Relaxation
RELAX(U1,FALSDT,0.1)
RELAX(V1,FALSDT,0.1)
RELAX(W1,FALSDT,0.03)
In-Form statements
In-Form statement calculates the enthalpy
INFORM7BEGIN
(STORED var ENT1 is TEM1*CP1)
INFORM7END
char(fluid_name) ! Declared now in o14.htm
Choose a fluid, by moving the chosen fluid_name so that it is the
lowest which starts in the second or third column.
The properties are those which are stored in macro 089.htm
*** GASES AT ATMOSPHERIC PRESSURE ***
fluid_name=Air
fluid_name=Ammonia
fluid_name=Argon
fluid_name=Carbon_Dioxide
fluid_name=Carbon_Monoxide
fluid_name=Helium
fluid_name=Hydrogen
fluid_name=Methane
fluid_name=Nitrogen
fluid_name=Oxygen
fluid_name=superheated_water_vapor
fluid_name=saturated_water_vapor
fluid_name=saturated_vapor_Refrigerant-12
fluid_name=saturated_vapor_Refrigerant-134a
fluid_name=Air_at_20_deg_C_1_atm
fluid_name=Ammonia_at_27_deg_C
fluid_name=Carbon_Dioxide_at_27_deg_C_1_atm
fluid_name=Carbon_Monoxide_at_27_deg_C_1_atm
fluid_name=Helium_at_27_deg_C_1_atm
fluid_name=Hydrogen_at_27_deg_C_1_atm
fluid_name=Nitrogen_at_27_deg_C_1_atm
fluid_name=Oxygen_at_27_deg_C_1_atm
fluid_name=STEAM_at_100_deg_C_1_atm
*** LIQUIDS ***
fluid_name=saturated_water
fluid_name=SAE_5W-30_engine_oil
fluid_name=SAE_10W-30_engine_oil
fluid_name=SAE_20W-20_engine_oil
fluid_name=Ethylene_Glycol
fluid_name=Ethylene_Glycol_50%_by_volume_aqueous_solution
fluid_name=Gasoline
fluid_name=Glycerin
fluid_name=Refrigerant-12
fluid_name=Refrigerant-134a
fluid_name=Therminol_59
fluid_name=Therminol_66
fluid_name=Dowtherm_A
fluid_name=Syltherm_800
fluid_name=FC-72
fluid_name=HFE-7100
fluid_name=Mercury
fluid_name=Water_at_0_deg_C_1_atm
fluid_name=GLASS_Liq_at_1127_deg_C
fluid_name=Coolanol_25_at_25_deg_C_1 atm
fluid_name=Coolanol_45_at_25_deg_C_1_atm
fluid_name=ETHYLENE_GLYCOL_at_27_deg_C
fluid_name=FC75_at_25_deg_C_1_atm
fluid_name=Freon_(refrigerant-12)_at_27_deg_C
fluid_name=Glycerin_at_27_deg_C
fluid_name=Mercury_at_27_deg_C
fluid_name=WATER_at_20_deg_C
------------------------
load fluid-property file
------------------------
text(3D moving-wall; fluid is :fluid_name:..:libref:
l($089)
LIBREF = 762
DISTIL=T
EX(P1)=5.608E-03; EX(U1)=6.330E-02; EX(V1)=4.840E-02
EX(W1)=3.938E-03; EX(TEM1)=2.863E+02; EX(ENT1)=2.894E+05
EX(COND)=2.542E-02; EX(RHO1)=1.215E+00
EX(ENUL)=1.479E-05; EX(CP1)=1.004E+03
PHOTON USE
p;;;;;
vi 1 2 3
gr ou x 1
con tem1 z 5 fi;0.001;gr ou x m;gr ou y 1;gr ou y m
pause
con cp1 z 5 fi;0.001 ;gr ou x m;gr ou y 1;gr ou y m
pause
con enul z 5 fi;0.001;gr ou x m;gr ou y 1;gr ou y m
pause
con rho1 z 5 fi;0.001;gr ou x m;gr ou y 1;gr ou y m
pause
con cond z 5 fi;0.001;gr ou x m;gr ou y 1;gr ou y m
pause
con ent1 z 5 fi;0.001;vec z 5 ;vec y m
gr ou x m;gr ou y 1;gr ou y m
ENDUSE
STOP