FORMAC routine for computing the matrix elements of the linearized cross-collision operators of gaseous mixtures
✍ Scribed by Francis J. McCormack
- Publisher
- Elsevier Science
- Year
- 1975
- Tongue
- English
- Weight
- 539 KB
- Volume
- 9
- Category
- Article
- ISSN
- 0010-4655
No coin nor oath required. For personal study only.
✦ Synopsis
A FORMAC program is described which calculates the derivatives of appropriate generating functions to obtain the matrix elements with respect to the Burnett functions of the linearized Boltzmann cross-collision operators of a gaseous mixture. The general procedure discussed, however, is applicable to any linear operator for which the generating function for the matrix elements is known1 LET(ELFELF+FXPANDC(EL*EP*FRD/CFKL*ELMt *FKP*FRKR *FAC CM PP) *F AC C PP.KRP) *EVAL(DERIvCFRI,Y,KRP),y,O) *EVAL(DER!V(FR3,Y,PP-KPP.~,Y,0)))'~i END, END, LETCEL E:ELEI-ELE) ATOMIZE (ElF I) /*WRITE MATOIY ELEMENT IN TERMS OF CUAPMAN-CflWlING OMFr,A-INTFC,PAL 3*/ LET CEL,EgE~PAND(PEPLACECFLF,SP'JT,(1~COST**?)**(I/21)))I LFT(MATELF:fl) LFTCILLOWPfl~CE1E,G),IH~HIGHPflW(ELF,G)), Do T:TNTECEP(TL) TO IMTEGERCIH) BY ?,LET(T"T",k:T/?), I ET ( A =C OF FE C EL. F, C', **I)) LET( JL =L OAP(J W(A , COST) p .IH=HT GMPOW CA, COST)) I JL~IMTFGFPCJL),JH:TNTEGFP(JH) j DO JJN TO JL BY -1,LET(j:"J"), IF .1=JM THEN I ETCDEXPANDCRFPI ACE(A,COST**J, I-OMEGA (J,K ) ) ) ) p EL SF LET(O:FXPAND(RFPLACE(D,CflsT**J,j.OMEGACJ,K))))t DO J:Jt TO JI4,LET(JItJ~),