TITLE Temperature gradient experiment using DE v.6 w/ new MIXT segs !Total helium fraction is fixed: !integral over duct is used as a target at end of model. !->C:\Documents and Settings\092710\My Documents\ChemicalThermoac\GradT\tempgradTyp.out !Created@08:28:08 15-Jul-2009 with DeltaEC version 6.3b2!under win32, using Win 5.1.2600 (Service Pack 3) under Python DeltaEC. !Total helium fraction is fixed: !integral over duct is used as a target at end of model. !--------------------------------- 0 --------------------------------- BEGIN Acoustics targeted at the endpoints. !Segments 25 cm length, as if 1 m total length. !225 C in middle of duct. !Measured temperatures are used to set the temperature boundary conditions. 8.0000E+04 a Mean P Pa 200.00 b Freq Hz 294.00 c TBeg K 2438.6 d |p| Pa 0.0000 e Ph(p) deg 7.1114E-05 f |U| m^3/s G 190.10 g Ph(U) deg G -22.887 h Htot W G 0.0000 i Ndot mol/s 0.60806 j nL G HeAr Gas type !--------------------------------- 1 --------------------------------- RPN Geometry of duct. a) dia. in in. A) x-sec B) perim 0.1310 a G or T 8.6956E-06 A ChngeMe 1.0453E-02 B ChngeMe 2 pi * inp * 2 / 2.54e-2 * pi inp * inp * 2.54e-2 * 2.54e-2 * 4 / !--------------------------------- 2 --------------------------------- RPN a = GasA/A, A = Atot, B = Radius (m) 4.8000E-02 a G or T 1.8116E-04 A ChngeMe 1.6637E-03 B ChngeMe 1a 2 / 2.54e-2 * 1A inp / !--------------------------------- 3 --------------------------------- RPN Gas displacement, |x1| 0.0000 a G or T 6.5079E-03 A ChngeMe U1 mag 1A / w / !--------------------------------- 4 --------------------------------- MIXTCIRC 1st qtr duct sameas 2A a Area m^2 2507.3 A |p| Pa sameas 2a b GasA/A 43.700 B Ph(p) deg 0.2500 c Length m 9.0169E-05 C |U| m^3/s sameas 2B d radius m -131.79 D Ph(U) deg -22.887 E Htot W -0.11269 F Edot W 0.60806 G nL Beg 0.56046 H nL End 294.00 I TBeg K copper Solid type 377.40 J TEnd K !--------------------------------- 5 --------------------------------- RPN Heat input/loss to meet temperature target -9.5301 a G or T G -32.417 A ChngeMe H2k inp + =H2k !--------------------------------- 6 --------------------------------- RPN Temperature boundary condition 377.40 a G or T =6A 377.40 A ChngeMe Tm !--------------------------------- 7 --------------------------------- MIXTCIRC 2nd qtr duct sameas 2A a Area m^2 2908.1 A |p| Pa sameas 2a b GasA/A 84.848 B Ph(p) deg sameas 4c c Length m 1.0584E-04 C |U| m^3/s sameas 2B d radius m -96.63 D Ph(U) deg -32.417 E Htot W -0.15384 F Edot W 0.56046 G nL Beg 0.51402 H nL End 377.40 I TBeg K sameas 4 Solid type 497.30 J TEnd K !--------------------------------- 8 --------------------------------- RPN Set Temperature 497.30 a G or T =8A 497.30 A ChngeMe Tm !--------------------------------- 9 --------------------------------- RPN Heater flux !Heat into the duct at midpoint 65.893 a G or T G 33.476 A ChngeMe H2k inp + =H2k !--------------------------------- 10 --------------------------------- RPN |p1|, re(P1), im(P1) 3000.0 a G or T 2908.1 A ChngeMe 261.14 B ChngeMe 2896.3 C ChngeMe p1 imag p1 real p1 mag !--------------------------------- 11 --------------------------------- RPN Impedance at midplane. 1.0000 a G or T 1.0152 A ChngeMe p1 mag U1 mag / rho / a / 1A * !--------------------------------- 12 --------------------------------- RPN Phase of impedance 180.00 a G or T 181.48 A ChngeMe p1 arg U1 arg - !--------------------------------- 13 --------------------------------- RPN Conc. at middle of duct. (Kludge, should not be const here.) 0.5000 a G or T 0.51402 A ChngeMe nL !--------------------------------- 14 --------------------------------- MIXTCIRC 3rd qtr duct sameas 2A a Area m^2 3386.3 A |p| Pa sameas 2a b GasA/A 127.26 B Ph(p) deg sameas 4c c Length m 1.1132E-04 C |U| m^3/s sameas 2B d radius m -52.459 D Ph(U) deg 33.476 E Htot W -0.18848 F Edot W 0.51402 G nL Beg 0.44957 H nL End 497.30 I TBeg K sameas 4 Solid type 372.60 J TEnd K !--------------------------------- 15 --------------------------------- RPN Heat input/loss to meet temperature target -13.12 a G or T G 20.356 A ChngeMe H2k inp + =H2k !--------------------------------- 16 --------------------------------- RPN Temperature boundary condition 372.60 a G or T =16A 372.60 A ChngeMe Tm !--------------------------------- 17 --------------------------------- MIXTCIRC 4th (and final) qtr duct sameas 2A a Area m^2 4073.1 A |p| Pa sameas 2a b GasA/A 179.06 B Ph(p) deg sameas 4c c Length m 1.1368E-04 C |U| m^3/s sameas 2B d radius m -6.4843 D Ph(U) deg 20.356 E Htot W -0.23042 F Edot W 0.44957 G nL Beg 0.37974 H nL End 372.60 I TBeg K sameas 4 Solid type 297.50 J TEnd K !--------------------------------- 18 --------------------------------- RPN Temperature Boundary Condition 297.50 a G or T =18A 297.50 A ChngeMe Tm !--------------------------------- 19 --------------------------------- RPN Lower speaker real pressure -4072.6 a G or T =19A -4072.6 A ChngeMe p1 real !--------------------------------- 20 --------------------------------- RPN Lower speaker imag pressure 66.500 a G or T =20A 66.500 A ChngeMe p1 imag !--------------------------------- 21 --------------------------------- RPN B/C: make avg of endpoints 0.5 0.5000 a G or T 0.4939 A ChngeMe 17H 0j + 2 / !--------------------------------- 22 --------------------------------- RPN dTm/dx 1 !Temperature gradient over 1st half of duct 0.0000 a G or T -406.6 A ChngeMe 0c 8a - 2 * !--------------------------------- 23 --------------------------------- RPN dn/dx 1 !Concentration gradient for light component over 1st half of duct 0.0000 a G or T 0.18807 A ChngeMe 0j 7H - 2 * !--------------------------------- 24 --------------------------------- RPN dTm/dx 2 0.0000 a G or T 399.60 A ChngeMe 8a 18a - 2 * !--------------------------------- 25 --------------------------------- RPN dn/dx 2 0.0000 a G or T 0.26856 A ChngeMe 7H 17H - 2 * !--------------------------------- 26 --------------------------------- RPN Integral of n_L dV 0.0000 a G or T 7.5480 A ChngeMe 0j 17H + 3.2 * 0j 17H + 14H 7H + 4H + 2 * + 2 / 2.174 * + !--------------------------------- 27 --------------------------------- RPN Calculated n_L for entire device 0.5000 a G or T =27A 0.5000 A ChngeMe 26A 3.2 2 * 2.174 4 * + / ! The restart information below was generated by a previous run ! and will be used by DeltaEC the next time it opens this file. guessz 0f 0g 0h 0j 5a 9a 15a xprecn -3.7375E-10 -4.8586E-03 -2.0955E-04 -1.3948E-06 6.4691E-06 1.7866E-04 -1.1722E-05 targs 6a 8a 16a 18a 19a 20a 27a ! Plot start, end, and step values. May be edited if you wish. ! Outer Loop: | Inner Loop .