How to improve the pressure

sun haitao

I have a question. It is a transient problemuration--0.01s. In a rectangle cavity, there is a conductor with Fixed heat flux(1000000w). The cavity only has a outlet and fills with air. The result of pressure and velocity are very small. I have some experimental result, and the pressure should be 2 bar. I don't know how to calculate this problem. I have appended the cheng4.q1 file with the letter.

Best Regards!

cheng4.q1: TALK=T;RUN( 1, 1)

************************************************** **********

Q1 created by VDI menu, Version 3.4, Date 10/08/01 CPVNAM=VDI;SPPNAM=Core ************************************************** ********** IRUNN = 1 ;LIBREF = 0 ************************************************** ********** Group 1. Run Title TEXT(Test ) ************************************************** ********** Group 2. Transience STEADY=F

* Set overall time and no. of steps RSET(U,0.000000E+00,1.000000E-02,10)

* Modify regions ************************************************** ********** Groups 3, 4, 5 Grid Information

* Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,20,20,20)

* Set overall domain extent:

* xulast yvlast zwlast

name XSI= 7.000000E-03; YSI= 2.000000E-03; ZSI= 3.500000E-03 RSET(D,CHAM ) ************************************************** ********** Group 6. Body-Fitted coordinates ************************************************** ********** Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T

* Non-default variable names NAME(147) =WDIS ; NAME(148) =PRPS NAME(149) =TEM1 ; NAME(150) =SPH1

* Solved variables list SOLVE(P1 ,U1 ,V1 ,W1 ,TEM1)

* Stored variables list STORE(SPH1,PRPS,WDIS)

* Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) TURMOD(LVEL)

************************************************** ********** Group 8. Terms & Devices ************************************************** ********** Group 9. Properties SETPRPS(1, 0) RHO1 = 1.189000E+00 PRESS0 = 1.000000E+05 TEMP0 = 2.730000E+02 CP1 = 1.005000E+03 ENUL = 1.544000E-05 DVO1DT = 3.410000E-03 PRNDTL(TEM1) = 1.000000E+00 ************************************************** ********** Group 10.Inter-Phase Transfer Processes ************************************************** ********** Group 11.Initialise Var/Porosity Fields FIINIT(P1 ) = 1.000000E+00 ;FIINIT(WDIS) = 2.000000E-04 FIINIT(TEM1) = 1.000000E+01

No PATCHes used for this Group

INIADD = F ************************************************** ********** Group 12. Convection and diffusion adjustments

No PATCHes used for this Group ************************************************** ********** Group 13. Boundary & Special Sources

No PATCHes used for this Group

EGWF = T ************************************************** ********** Group 14. Downstream Pressure For PARAB ************************************************** ********** Group 15. Terminate Sweeps LSWEEP = 300 RESFAC = 1.000000E-03 ************************************************** ********** Group 16. Terminate Iterations ************************************************** ********** Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) ************************************************** ********** Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06 ************************************************** ********** Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T ASAP = T RADIA = 0.000000E+00 ;RADIB = 0.000000E+00 EMISS = 0.000000E+00 CSG1 ='M' ************************************************** ********** Group 20. Preliminary Printout ECHO = T ************************************************** ********** Group 21. Print-out of Variables OUTPUT(WDIS,Y,N,N,N,N,N) OUTPUT(SPH1,N,N,Y,Y,N,N) ************************************************** ********** Group 22. Monitor Print-Out IXMON = 16 ;IYMON = 15 ;IZMON = 6 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************** ********** Group 23.Field Print-Out & Plot Control NPRINT = 100000 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 10000 ISWPRF = 1 ;ISWPRL = 100000

No PATCHes used for this Group ************************************************** ********** Group 24. Dumps For Restarts NOWIPE = T IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0 CSG1 ='M'

GVIEW(P,0.000000E+00,-1.000000E+00,0.000000E+00) GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00)

> DOM, SIZE, 7.000000E-03, 2.000000E-03, 3.500000E-03
: DOM, MONIT, 5.425000E-03, 1.450000E-03, 9.350000E-04
: DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
: DOM, SNAPSIZE, 1.000000E-02
: DOM, RELAX, 5.000000E-01

> OBJ1, NAME, B2
: OBJ1, POSITION, 5.000000E-04, 2.500000E-04, 1.700000E-03
: OBJ1, SIZE, 5.000000E-04, 1.500000E-03, 5.000000E-04
: OBJ1, CLIPART, cube4
: OBJ1, ROTATION24, 1
: OBJ1, TYPE, BLOCKAGE
: OBJ1, MATERIAL, 103
: OBJ1, TIME_LIMITS, ALWAYS_ACTIVE
: OBJ1, HEAT_FLUX, 0.000000E+00, 1.100000E+06

> OBJ2, NAME, PLATE1
: OBJ2, POSITION, 3.500000E-03, 0.000000E+00, 5.000000E-04
: OBJ2, SIZE, 5.000000E-04, 1.500000E-03, 2.000000E-04
: OBJ2, CLIPART, cube14
: OBJ2, ROTATION24, 1
: OBJ2, TYPE, BLOCKAGE
: OBJ2, MATERIAL, 111
: OBJ2, TIME_LIMITS, ALWAYS_ACTIVE
: OBJ2, ADIABATIC, 0.000000E+00, 0.000000E+00

> OBJ3, NAME, B4
: OBJ3, POSITION, 3.500000E-03, 0.000000E+00, 1.000000E-03
: OBJ3, SIZE, 5.000000E-04, 1.500000E-03, 2.000000E-04
: OBJ3, CLIPART, cube14
: OBJ3, ROTATION24, 1
: OBJ3, TYPE, BLOCKAGE
: OBJ3, MATERIAL, 111
: OBJ3, TIME_LIMITS, ALWAYS_ACTIVE
: OBJ3, ADIABATIC, 0.000000E+00, 0.000000E+00

> OBJ4, NAME, B5
: OBJ4, POSITION, 3.500000E-03, 0.000000E+00, 1.500000E-03
: OBJ4, SIZE, 5.000000E-04, 1.500000E-03, 2.000000E-04
: OBJ4, CLIPART, cube14
: OBJ4, ROTATION24, 1
: OBJ4, TYPE, BLOCKAGE
: OBJ4, MATERIAL, 111
: OBJ4, TIME_LIMITS, ALWAYS_ACTIVE
: OBJ4, ADIABATIC, 0.000000E+00, 0.000000E+00

> OBJ5, NAME, B6
: OBJ5, POSITION, 3.500000E-03, 0.000000E+00, 2.000000E-03
: OBJ5, SIZE, 5.000000E-04, 1.500000E-03, 2.000000E-04
: OBJ5, CLIPART, cube14
: OBJ5, ROTATION24, 1
: OBJ5, TYPE, BLOCKAGE
: OBJ5, MATERIAL, 111
: OBJ5, TIME_LIMITS, ALWAYS_ACTIVE
: OBJ5, ADIABATIC, 0.000000E+00, 0.000000E+00

> OBJ6, NAME, B7
: OBJ6, POSITION, 3.500000E-03, 0.000000E+00, 2.500000E-03
: OBJ6, SIZE, 5.000000E-04, 1.500000E-03, 2.000000E-04
: OBJ6, CLIPART, cube14
: OBJ6, ROTATION24, 1
: OBJ6, TYPE, BLOCKAGE
: OBJ6, MATERIAL, 111
: OBJ6, TIME_LIMITS, ALWAYS_ACTIVE
: OBJ6, ADIABATIC, 0.000000E+00, 0.000000E+00

> OBJ7, NAME, OUTLET
: OBJ7, POSITION, 7.000000E-03, 5.000000E-04, 0.000000E+00
: OBJ7, SIZE, 0.000000E+00, 1.000000E-03, 3.500000E-03
: OBJ7, CLIPART, cube12t
: OBJ7, ROTATION24, 1
: OBJ7, TYPE, OUTLET
: OBJ7, PRESSURE, 0.000000E+00
: OBJ7, TEMPERATURE, -1.026000E+04
: OBJ7, COEFFICIENT, 1.000000E+03
: OBJ7, TIME_LIMITS, ALWAYS_ACTIVE

> OBJ8, NAME, B8
: OBJ8, POSITION, 5.000000E-03, 5.000000E-04, 8.000000E-04
: OBJ8, SIZE, 7.000000E-04, 2.000000E-04, 3.500000E-04
: OBJ8, CLIPART, default
: OBJ8, ROTATION24, 1
: OBJ8, TYPE, PROBE
: OBJ8, TEMPERATURE, 0.000000E+00
: OBJ8, PRESSURE, 0.000000E+00
: OBJ8, VELOCITY, 0.000000E+00, 0.000000E+00, 0.000000E+00 STOP