Fichier q1 en instationnaire

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

IRUNN = 1 ;LIBREF = 166 ************************************************************

Group 1. Run Title

TEXT(INSTABILITE )

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

RESTRT(ALL)

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

** coeff Rayleigh **

REAL(DDT,DDT2)

REAL(RAYL1,RAYL2)

DDT=.145

DDT2=.15

** coefficients de relaxation **

REAL(coef,FDT,Vref)

coef=100

Vref=0.00001

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

Group 2. Transience

STEADY = F

LSTEP=10

TLAST=50

GRDPWR(T,LSTEP,TLAST,1)

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

Y,-NY,YVLAST,1)

GRDPWR(Z,NZ,ZWLAST,1)

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

Group 6. Body-Fitted coordinates

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

Group 7. Variables: STOREd,SOLVEd,NAMEd

ONEPHS = T

* Non-default variable names

NAME( 45) =DEN1

* Solved variables list

SOLVE(P1 ,U1 ,V1 ,H1 )

* Stored variables list

STORE(IMB1 ,DEN1,VORT,VORT1)

* Additional solver options

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

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

Group 8. Terms & Devices

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

Group 9. Properties

RHO1 = 1.000E+03

CP1 = 4.182E+03

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

PRNDTL(H1 ) = 7.020E+00

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

Group 10.Inter-Phase Transfer Processes

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

Group 11.Initialise Var/Porosity Fields

No PATCHes used for this Group

RSTGRD = F

INIADD = F

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

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

Group 13. Boundary & Special Sources

PATCH (BUOYANCY,PHASEM,1,NX,1,NY,1,1,1,LSTEP)

COVAL (BUOYANCY,U1 , FIXFLU , GRND3 )

COVAL (BUOYANCY,V1 , FIXFLU , GRND3 )

PATCH (FIXPRESS,CELL ,1,1,1,1,1,1,1,LSTEP)

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

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

PATCH (MUR ,NWALL ,1,NX,NY,NY,1,1,1,LSTEP)

COVAL (MUR ,U1 , GRND1 , 0.000E+00)

COVAL (MUR ,H1 , 1.425E-01,0.0)

** Iterations sur les cond initiales **

PATCH (MUR2 ,SWALL ,1,NX,1,1,1,1,1,2)

COVAL (MUR2 ,U1 , GRND1 , 0.000E+00)

COVAL (MUR2 ,H1 , 1.425E-01, 4180*DDT)

** Nouvelles C. Limites **

PATCH (MUR2 ,SWALL ,1,NX,1,1,1,1,1,LSTEP)

COVAL (MUR2 ,U1 , GRND1 , 0.000E+00)

COVAL (MUR2 ,H1 , 1.425E-01, 4180*DDT2)

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

BUOYD =-4.3E-08

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

Group 14. Downstream Pressure For PARAB

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

Group 15. Terminate Sweeps

LSWEEP = 3000

SELREF = T

RESFAC = 1.000E-03

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

Group 16. Terminate Iterations

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

Group 17. Relaxation

RELAX(P1 ,LINRLX, 1.000E+00)

RELAX(U1 ,FALSDT, coef*FDT)

RELAX(V1 ,FALSDT, coef*FDT)

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

Group 18. Limits

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

Group 19. EARTH Calls To GROUND Station

GENK = T

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

Group 20. Preliminary Printout

ECHO = T

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

Group 21. Print-out of Variables

OUTPUT(DEN1,Y,N,Y,N,N,N)

OUTPUT(IMB1,Y,N,Y,N,N,N)

OUTPUT(VORT,Y,Y,Y,Y,Y,Y)

OUTPUT(U1,Y,Y,Y,Y,Y,Y)

OUTPUT(V1,Y,Y,Y,Y,Y,Y)

OUTPUT(H1,Y,Y,Y,Y,Y,Y)

OUTPUT(P1,Y,Y,Y,Y,Y,Y)

IDISPA=1

IDISPB=1

IDISPC=LSTEP

IDISPD=1

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

Group 22. Monitor Print-Out

IXMON = 10 ;IYMON = 10 ;IZMON = 1

TSTSWP = -5

UWATCH = T

USTEER = T

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

Group 23.Field Print-Out & Plot Control No PATCHes used for this Group

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

Group 24. Dumps For Restarts

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

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

MENSAV(S,PHSPROP,DEF,200,0,1.0005E+03,1.0060E-06)

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

RAYL1=9.81*CP1*BUOYD*DDT*YVLAST*YVLAST*YVLAST*PRNDTL(H1)

RAYL1=RAYL1/ENUL/ENUL

RAYL2=9.81*CP1*BUOYD*DDT2*YVLAST*YVLAST*YVLAST*PRNDTL(H1)

MESG(NOUVEAU RAYLEIGH VAUT:RAYL2:

STOP

Retour à la page de garde