1 : * ASPARAM PROCEDUR MAGN 08/12/09 21:15:00 6189 2 : 'DEBPROC' asparam rxt*TABLE ; 3 : 4 : TBT = rxt.'TBT' ; 5 : GEO = rxt.'GEO' ; 6 : * 7 : idim = vale 'DIME'; 8 : DIM3D=FAUX ; 9 : Si(EGA idim 3) ; 10 : DIM3D=VRAI ; 11 : Finsi ; 12 : * MESS 'si VRAI dimension 3D' ; 13 : * LIST DIM3D ; 14 : * 15 : un = rxt.'TIC'.'UN' ; 16 : vn = rxt.'TIC'.'VN' ; 17 : xd = rxt.'TIC'.'XD' ; 18 : xdinj = 'IPOL' 0.0 t (rxt.'scenasp'.'xdinj') ; 19 : xdmin = xdinj * 0.01 ; 20 : ixd=(masq xd 'SUPERIEUR' xdmin ) ; 21 : tf = rxt.'TIC'.'TF' ; 22 : td = rxt.'TIC'.'TD' ; 23 : dd = rxt.'TIC'.'DD' ; 24 : t = rxt.'scenasp'.'t' ; 25 : ddinj = 'IPOL' 0.0 t (rxt.'scenasp'.'ddinj') ; 26 : ddmin = ddinj * 0.01 ; 27 : idd=(masq dd 'SUPERIEUR' ddmin ) ; 28 : dd = rxt.'TIC'.'DD' *idd ; 29 : rhogas = rxt.'TIC'.'RHO' ; 30 : rhovap = rxt.'TIC'.'RVP' ; 31 : mugas = rxt.'TIC'.'Mu' ; 32 : Cpgas = extr (rxt.'TIC'.'Cpm') (dime (rxt.'TIC'.'Cpm')) ; 33 : lamgas = TBT.'lamair' ; 34 : rhodis = rxt.'rod' ; 35 : Cpdis = rxt.'Cpd' ; 36 : $vtf = rxt.'GEO'.'$vtf' ; 37 : Lv = TBT.'Lv' ; 38 : dt = rxt.'DT0' ; 39 : * 40 : *corrections appel procedur rvst2.procedur 41 : * rhosat rhosp = rvsat td ; 42 : rgva = TBT.'Rgvap' ; 43 : Pvap = rgva * rhovap * (tf+273.0) + 0.000005 ; 44 : ind200= td masque 'EGINF' 200. ; 45 : ind20sup = td masque 'EGSUP' 20. ; 46 : indch= ind20sup * ind200 ; 47 : * 48 : ** 49 : td1 = td * indch ; 50 : rhosat1 rhosp1 = rvsat td1 ; 51 : rhosp1 = rhosp1 * indch ; 52 : rhosat1 = rhosat1 * indch ; 53 : ind20inf = td masque 'EGINF' 19.999999999999999 ; 54 : ind0sup = td masque 'EGSUP' 0. ; 55 : indfr = ind20inf * ind0sup ; 56 : td2 = td * indfr ; 57 : rhosat2 rhosp2 = rvst2 td2 ; 58 : rhosat2 = rhosat2 * indfr ; 59 : rhosp2 = rhosp2 * indfr ; 60 : rhosat = rhosat1 + rhosat2 ; 61 : rhosp = rhosp1 + rhosp2 ; 62 : * 63 : *fin correction appel procedur 64 : * 65 : Pt = extr (dime rxt.'TIC'.'PT') rxt.'TIC'.'PT' ; 66 : Dvap = (((tf + 273.0)/273.0) ** 1.75) * 2.2 / Pt ; 67 : * 68 : drho = noel $vtf ((rhodis - rhogas)/rhodis) ; 69 : 'SI' DIM3D ; 70 : gdrx = nomc 'UX' (drho * 0.00) ; 71 : gdry = nomc 'UY' (drho * 0.00) ; 72 : gdrz = nomc 'UZ' (drho * 9.81) ; 73 : rxt.'TIC'.'GRD' = kcht $vtf vect centre (gdrx et gdry et gdrz) ; 74 : 'SINON' ; 75 : gdrx = nomc 'UX' (drho * 0.00) ; 76 : gdrz = nomc 'UY' (drho * 9.81) ; 77 : rxt.'TIC'.'GRD' = kcht $vtf vect centre (gdrx et gdrz) ; 78 : FINSI; 79 : * 80 : cmps = extr un comp ; 81 : ur = (psca (vn - un) (vn - un) cmps cmps) ** 0.5 ; 82 : Re = dd * ur * rhogas * (mugas** -1.) ; 83 : Pr = mugas * Cpgas * (lamgas** -1.) ; 84 : Sc = mugas * (rhogas ** -1.0) * (Dvap ** -1.0) ; 85 : * modif 1 2004 86 : Nu = 2.00 + (0.599 * ((Pr ** 0.33333) * (Re ** 0.5))) ; 87 : * Nu = (2.00 + 0.599 * ((Pr ** 0.33333) * (Re ** 0.5))) ; 88 : * modif 2 2004 89 : Sh = 2.00 + (0.599 * ((Sc ** 0.33333) * (Re ** 0.5))) ; 90 : * Sh = (2.00 + 0.599 * ((Sc ** 0.33333) * (Re ** 0.5))) ; 91 : * 92 : F_D = (1.00 + (0.15 * (Re ** 0.687))) ; 93 : * 94 : fk = 18.0 * mugas * F_D * (dd ** -2.0) ; 95 : fh = 6.0 * lamgas * Nu * (dd ** -2.0) ; 96 : fm = 6.0 * Dvap * Sh * (dd ** -2.0) * (rhovap - rhosat) ; 97 : 98 : 99 : 100 : *fmax = (tf - td) * ((rhodis * Cpdis / dt) - (fh / Cpgas)) / Lv; 101 : *fmin = (-1.0*td) * ((rhodis * Cpdis / dt) - (fh / Cpgas)) / Lv; 102 : * 103 : num = ((rhovap - rhosat)/dt) * ((xd + 1.0e-10) ** -1.0) ; 104 : den = (rhosp * ((xd + 1.0e-10) ** -1.0)) ; 105 : den = (den * Lv / (Cpdis * rhodis)) + 1.0 ; 106 : fmax = abs (num * (den ** -1.0)) ; 107 : fmin = -1.0 * fmax ; 108 : * 109 : f1 = fm masque superieur fmax ; 110 : f2 = fm masque inferieur fmin ; 111 : f3 = (1.0 - f1) * (1.0 - f2) ; 112 : fm = (fmax * f1) + (fmin * f2) + (fm * f3) ; 113 : * fm = fm * 0.0 ; 114 : 115 : kgas = fk * xd * (rhogas ** -1.0) * ixd ; 116 : kdis = fk / rhodis ; 117 : hgas = fh * xd / Cpgas * (rhogas ** -1.0) * ixd ; 118 : hdis = fh / Cpdis / rhodis ; 119 : sdia = fm * dd / 3.0 / rhodis ; 120 : sxd = fm * xd / rhodis ; 121 : srvp = fm * xd * (-1.0) ; 122 : *srho = fm * xd * (-1.0) ; 123 : std = fm * (Lv) / Cpdis / rhodis ; 124 : 125 : 126 : * rxt.'TIC'.'GRD' = kcht $vtf vect centre (gdrx et gdry et gdrz) ; 127 : rxt.'TIC'.'KGAS' = kcht $vtf scal sommet kgas ; 128 : rxt.'TIC'.'KDIS' = kcht $vtf scal sommet kdis ; 129 : rxt.'TIC'.'HGAS' = kcht $vtf scal sommet hgas ; 130 : rxt.'TIC'.'HDIS' = kcht $vtf scal sommet hdis ; 131 : rxt.'TIC'.'SDIA' = noel $vtf (kcht $vtf scal sommet sdia) ; 132 : rxt.'TIC'.'SXD' = noel $vtf (kcht $vtf scal sommet sxd ) ; 133 : *rxt.'TIC'.'SRHO' = noel $vtf (kcht $vtf scal sommet srho) ; 134 : rxt.'TIC'.'SRVP' = noel $vtf (kcht $vtf scal sommet srvp) ; 135 : rxt.'TIC'.'STD' = noel $vtf (kcht $vtf scal sommet std ) ; 136 : rxt.'TIC'.'fm' = kcht $vtf scal sommet fm ; 137 : rxt.'TIC'.'fmax' = kcht $vtf scal sommet fmax ; 138 : 139 : 140 : 141 : * mess 'Mini/Max kdis ' (mini kdis) (maxi kdis) ; 142 : * mess 'Mini/Max hdis ' (mini hdis) (maxi hdis) ; 143 : * mess 'Mini/Max kgas ' (mini kgas) (maxi kgas) ; 144 : * mess 'Mini/Max hgas ' (mini hgas) (maxi hgas) ; 145 : * mess 'Mini/Max fk ' (mini fk ) (maxi fk ) ; 146 : * mess 'Mini/Max fh ' (mini fh ) (maxi fh ) ; 147 : * mess 'Mini/Max fm ' (mini fm ) (maxi fm ) ; 148 : * mess 'Mini/Max srvp ' (mini srvp) (maxi srvp) ; 149 : * mess 'Mini/Max xd ' (mini xd ) (maxi xd ) ; 150 : * mess 'Mini/Max rhogas' (mini rhogas) (maxi rhogas) ; 151 : * mess 'Cpgas' Cpgas 'Lv=' lv ; 152 : * mess 'Mini/Max mugas' (mini mugas) (maxi mugas) ; 153 : * mess 'Mini/Max lamgas' lamgas ; 154 : * mess 'Mini/Max Dvap ' (mini Dvap ) (maxi Dvap ) ; 155 : 156 : 157 : 'FINPROC' ; 158 :
© Cast3M 2003 - All rights reserved.
Disclaimer