Fichier q1 ayant servis pour le calcul du Rayleigh critique


TALK=T;RUN( 1, 1);VDU=X11-TERM

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

Q1 created by GENERAL menu, Version 2.2, Date 26/03/96

CHAR(CPVNAM);CPVNAM=GENERAL

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

IRUNN = 1 ;LIBREF = 0

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

Group 1. Run Title

TEXT(ESSAI 1 )

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

Group 2. Transience

REAL (hhaut,hbas)

hhaut=0

hbas=550

STEADY = T

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

Groups 3, 4, 5 Grid Information

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

RSET(M,40,20,1,1.000E-05)

* Set overall domain extent:

* xulast yvlast zwlast name

XSI= 2.000E-02;YSI= 1.000E-02;ZSI= 1.000E+00;RSET(D,EAU )

* Set objects: x0 y0 z0

* dx dy dz name

XPO= 0.000E+00;YPO= 0.000E+00;ZPO= 0.000E+00

XSI= 2.000E-02;YSI= 0.000E+00;ZSI= 1.000E+00;RSET(B,INF )

XPO= 0.000E+00;YPO= 1.000E-02;ZPO= 0.000E+00

XSI= 2.000E-02;YSI= 0.000E+00;ZSI= 1.000E+00;RSET(B,SUP )

* Modify default grid

RSET(X,1,40,1.000E+00)

RSET(Y,1,20,1.000E+00)

RSET(Z,1,1,1.000E+00)

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

Group 6. Body-Fitted coordinates

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

Group 7. Variables: STOREd,SOLVEd,NAMEd

ONEPHS = T

NAME( 49) =DEN1

* Solved variables list

SOLVE(P1 ,U1 ,V1 ,H1 )

* Stored variables list STORE(DEN1,VMOD,IMB1)

* Additional solver options

SOLUTN(P1 ,Y,Y,Y,N,N,Y)

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

Group 8. Terms & Devices

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

Group 9. Properties

RHO1 = 9.982E+02

TMP1 = GRND2

TMP1A = 0.000E+00 ;TMP1B = 2.391E-04 ;TMP1C = 0.000E+00

CP1 = 4.182E+03

ENUL = 1.006E-06 ;ENUT = 0.000E+00

PRNDTL(H1 ) = 7.034E+00

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

Group 10.Inter-Phase Transfer Processes

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

Group 11.Initialise Var/Porosity Fields

No PATCHes used for this Group

PATCH (TEMPI,LINVLY,1,NX,1,NY,1,1,1,1) INIT(TEMPI,H1,(hhaut-hbas)/YVLAST,hbas)

RSTGRD = F INIADD = T

RESTRT (ALL)

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

Group 12. Convection and diffusion adjustments

No PATCHes used for this Group

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

Group 13. Boundary & Special Sources

PATCH (BUOYANCY,PHASEM,#1,#NREGX,#1,#NREGY,#1,#NREGZ,#1,#NREGT)

COVAL (BUOYANCY,U1 , FIXFLU , GRND3 )

COVAL (BUOYANCY,V1 , FIXFLU , GRND3 )

PATCH (FIXPRESS,CELL ,10,10,10,10,1,1,1,1)

COVAL (FIXPRESS,P1 , 1.000E+00, 0.000E+00)

COVAL (FIXPRESS,H1 , 0.000E+00, SAME )

PATCH (INF ,SWALL ,#1,#1,#1,#1,#1,#1,1,1)

COVAL (INF ,U1 , GRND2 , 0.000E+00)

COVAL (INF ,H1 , 1.422E-01, hbas)

PATCH (SUP ,NWALL ,#1,#1,#1,#1,#1,#1,1,1)

COVAL (SUP ,U1 , GRND2 , 0.000E+00)

COVAL (SUP ,H1 , 1.422E-01, hhaut)

BUOYA = 0.000E+00 ; BUOYB =-9.810E+00 ; BUOYC = 0.000E+00

BUOYD =-4.304E-08

EGWF = T

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

Group 14. Downstream Pressure For PARAB

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

Group 15. Terminate Sweeps

LSWEEP = 8000

SELREF = T

RESFAC = 1.000E-03

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

Group 16.

Terminate Iterations ************************************************************

Group 17. Relaxation

RELAX(P1 ,LINRLX, 1.000E+00)

RELAX(U1 ,FALSDT, 5.000E-02)

RELAX(V1 ,FALSDT, 5.000E-02)

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

Group 18. Limits

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

Group 19. EARTH Calls To GROUND Station

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

Group 20. Preliminary Printout

ECHO = T

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

Group 21. Print-out of Variables

INIFLD = T

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

Group 22. Monitor Print-Out

IXMON = 2 ;IYMON = 2 ;IZMON = 1

TSTSWP = -1

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

Group 23.Field Print-Out & Plot Control

No PATCHes used for this Group

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

Group 24. Dumps For Restarts

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

real (Ra)

Ra=-buoyd*(hbas-hhaut)*9.81*yvlast**3/(enul**2/prndtl(h1))

mesg(le Rayleigh est egal a :Ra:

MENSAV(S,RELX,DEF,5.0000E-04,1,100)

MENSAV(S,PHSPROP,DEF,200,0,9.9823E+02,1.0060E-06)

MENSAV(S,FLPRP,DEF,LAMINAR,CONSTANT)

m

STOP