|
| |
|
Главная Журналы первоначальные оценки. В то же время это поможет выявить некоторые ошибки, которые вы могли допустить при настройке CONDUCT. В любом случае тщательное изучение результатов обогатит ваше понимание лежащих в основе задачи физических явлений. Получайте удовольствие. Моделирование сложных физических явлений с использованием вычислительной программы является увлекательным занятием. С помощью численного моделирования можно удовлетворить свое научное любопытство точно так же, как при изучении физических явлений в лабораторных экспериментах. Такой инструмент, как CONDUCT, дает возможность моделировать широкий диапазон физических явлений. Вы поймете, что это наиболее захватывающий вид деятельности. Приложение 1 ЛИСТИНГ НЕИЗМЕНЯЕМОЙ ЧАСТИ ПРОГРАММЫ CONDUCT Файл «COMMON» сссссссссссссссссссссссссссссссссссссссссссссссссссссссссссс PARAMETER(N1=52,NJ=52,NFMAX=5,NZMAX=1О) CHARACTER*18 TITLE CHARACTER*64 HEADER,PRINTF,PLOTF COMMON F(NI,NJ,NFMAX),ALAM(NI,NJ),GAM(NI,NJ),CON(NI,NJ), 1 AIP(NI,NJ),AIM(NI,NJ),AJP(NI,NJ),AJM(NI,NJ),AP(NI,NJ), 2 IBLOCK(NI,NJ),X(NI),XU{NI),XCV{NI),ARX{NJ),Y(NJ), 3 YV(NJ),YCV(NJ),YCVR(NJ),R(NJ),RV(NJ),SX(NJ),PTX(NI), 4 QTX(NI),PTY{NJ),QTY{NJ),FLUXIl(NJ,NFMAX), 5 FLOXLl(NJ,NFMAX),FLUXJl(NI,NFMAX),FLUXMl(NI,NFMAX), 6 KBCIl(NJ),KBCLl(NJ),KBCJl(NI),KBCMl(NI), 7 FLXCIl(NJ),FLXCLl(NJ),FLXCJl(NI),FLXCMl(NI), 8 FLXPll(NJ),FLXPLl(NJ),FLXPJl(NI),FLXPMl(NI) COMMON/GENL/NF,Ll,L2,L3,Ml,M2,MB,ITER,SMALL,BIG, 1 LAST,TIME,DT,MODE,KORD,KOUT,lUl,I02,KPGR,KSTOP COMMON/NFE/RELAX(NFMAX),KPRINT(NFMAX),KSOLVE(NFMAX), 1 KBLOC(NFMAX),KPLOT(NFMAX),NTIMES(NFMAX),NTC(NFMAX), 2 CRIT(NFMAX) COMMON/TTL/TITLE(NFMAX),HEADER,PRINTF,PLOTF COMMON/EZG/NCVLX,NCVLY,XL,YL,POWERX,POWERY COMMON/ZG/NZX,NZY,NCVX(NZMAX),NCVY(NZMAX),XZONE(NZMAX), 1 YZONE(NZMAX),POWRX(NZMAX),POWRY(NZMAX) DIMENSION SC(NI,NJ),SP(NI,NJ) EQUIVALENCE (CON,SC),(AP,SP) CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC Неизменяемая часть программы CONDUCT CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC PROGRAM CONDUCT COMPUTER PROGRAM CONDUCT BY DR. SUHAS V. PATANKAR COPYRIGHT (C) 1991 INNOVATIVE RESEARCH, INC. SINCLUDE:COMMON CALCULATIONS IN THE GETTING-READY PHASE CALL DEFLT CALL GRID CALL READY CALL BEGIN 10 CONTINUE COME HERE TO START THE ITERATION OR TIME-STEP LOOP CALL OUTPUT IF(KSTOP.NE.0) STOP CALL HEART GO TO 10 END CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC SUBROUTINE DEFRD + * + + * + + + + + + + + * + + + + * + + * + + $INCLUDE:COMMON ENTRY DEFLT COME HERE TO SET THE DEFAULT VALUES С HEADER=USE THE CHARACTER VARIABLE HEADER TO SPECIFY A 1 PROBLEM TITLE PRINTF=PRINTl PLOTF=PLOTl CALL INTA7(KSTOP,0,LAST,5,ITER,0,KORD,2,MODE,1, 1 KPGR,1,K0UT,3) CALL DATA5(SMALL,1.E-2 0,BIG,1.E + 2 0,TIME,0.,DT,1.E+ 20, 1 R(1),0.) CALL DATA2(POWERX,1.,POWERY,1.) DO 10 NZ=1,NZMAX POWRX(NZ)=1. 10 POWRY{NZ)=1. DO 20 N=1,NFMAX CRIT(N)=1.E-5 KSOLVE(N)=0 NTIMES(N)=10 KBLOC(N)=1 RELAX(N)=1. TITLE(N)= KPRINT(N)=0 KPLOT(N)=0 DO 30 1=2,NI FLUXJl(I,N)=0. 30 FLUXMl(I,N)=0. DO 4 0 J=2,NJ FLUXIl(J,N)=0. 40 FLUXLl(J,N)=0. 20 CONTINUE DO 50 J=1,NJ DO 50 1=1,NI CON(I,J)=0. 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 [ 91 ] 92 93 94 95 96 97 98 99 |