TITLE with JOINs and STKSCREENS !->C:\Documents and Settings\092710\My Documents\TippyTBT\travelingWave\DeltaEC-expt-LongThin\ExtraQdot\1215-0138.out !Created@20:11:18 01-Jan-2010 with DeltaEC version 6.3b9!under win32, using Win 5.1.2600 (Service Pack 3) under Python DeltaEC. !Match f, gas, geom of course. !Match pm, complex p1 drive, complex p1 tank. !(guess inertance area and perim; watch these for sanity.) !Match TCold and THot. !Inject Qdot convection in RPN segment 8. !Results: half-stroke; tube temperatures; mid-tube phase. !--------------------------------- 0 --------------------------------- BEGIN Based on pulsetube1 rev1.dwg, !Begin at Endevco location, at the end of the long drill hole. 3.0860E+06 a Mean P Pa 25.000 b Freq Hz 287.63 c TBeg K G 4.9132E+04 d |p| Pa 19.100 e Ph(p) deg 0.0000 f |U| m^3/s 0.0000 g Ph(U) deg helium Gas type !--------------------------------- 1 --------------------------------- DUCT Drill hole to Endevco 1.3000E-05 a Area m^2 Mstr 4.9129E+04 A |p| Pa 1.2781E-02 b Perim m 1a 19.101 B Ph(p) deg 6.7000E-02 c Length m 1.4245E-06 C |U| m^3/s 5.0000E-04 d Srough -75.166 D Ph(U) deg 0.0000 E Htot W ideal Solid type -2.6036E-03 F Edot W !--------------------------------- 2 --------------------------------- BRANCH This is the piston -4.6691E+07 a Re(Zb) Pa-s/m^3 G 4.9129E+04 A |p| Pa 1.5729E+08 b Im(Zb) Pa-s/m^3 G 19.101 B Ph(p) deg 0.0000 c HtotBr W 2.9804E-04 C |U| m^3/s 92.509 D Ph(U) deg 0.0000 E Htot W 2.0905 F Edot W -2.0931 G EdotBr W !--------------------------------- 3 --------------------------------- RPN set correct Hdot 0.0000 a G or T 2.0905 A watts Edot =H2k !--------------------------------- 4 --------------------------------- DUCT entrance region 2.3950E-04 a Area m^2 Mstr 4.9155E+04 A |p| Pa 5.4860E-02 b Perim m 4a 19.091 B Ph(p) deg 2.7000E-02 c Length m 2.8851E-04 C |U| m^3/s 5.0000E-04 d Srough 91.985 D Ph(U) deg 2.0905 E Htot W ideal Solid type 2.0856 F Edot W !--------------------------------- 5 --------------------------------- RPN deg C to Kelvin conversion 11.400 a G or T 284.55 A kelvin inp 273.15 + !--------------------------------- 6 --------------------------------- TX 0.052 inch holes, 91 holes sameas 4a a Area m^2 4.9166E+04 A |p| Pa 0.5206 b GasA/A 19.080 B Ph(p) deg 5.6390E-03 c Length m 2.8716E-04 C |U| m^3/s 6.6040E-04 d radius m 91.950 D Ph(U) deg -3.6639 e HeatIn W G -1.5734 E Htot W sameas 5A f SolidT K =6H 2.0792 F Edot W 287.63 G GasT K copper Solid type 284.55 H SolidT K !--------------------------------- 7 --------------------------------- RPN calcs: want NR >~ 1000 and dp >~ rhov2/2/10 0.0000 a G or T 387.99 A NR1 100.31 B dp 13.697 C rhov2/2 rho 6C 6a / 6b / sto sqrd * 2 / ; p1 12A - mag ; rcl 4 * 9d * rho * mu / !--------------------------------- 8 --------------------------------- RPN insert convect Qdot here, avoiding TX but catching fs 9.8500 a G or T -11.423 A ChngeMe H2k inp - =H2k !--------------------------------- 9 --------------------------------- STKSCREEN 42 mesh, 0.0055 wire, mcmaster sameas 4a a Area m^2 4.9144E+04 A |p| Pa 0.8190 b VolPor 18.966 B Ph(p) deg 6.3500E-03 c Length m 2.9877E-04 C |U| m^3/s 1.5800E-04 d rh m 91.436 D Ph(U) deg 1.0000 e ksFrac -11.423 E Htot W 2.2112 F Edot W 287.63 G TBeg K stainless Solid type 303.36 H TEnd K !--------------------------------- 10 --------------------------------- RPN account for convective heat again sameas 8a a G or T -1.5734 A ChngeMe H2k inp + =H2k !--------------------------------- 11 --------------------------------- JOIN Change Me 4.9144E+04 A |p| Pa 18.966 B Ph(p) deg 2.9828E-04 C |U| m^3/s 91.436 D Ph(U) deg -1.5734 E Htot W 2.2076 F Edot W 303.36 G TBeg K 302.86 H TEnd K !--------------------------------- 12 --------------------------------- RPN p1 here. Also xi 1 0.0000 a G or T (4.6476E+04, 1.5972E+04) A Pa 7.9287E-03 B m U1 mag 13a / w / ; p1 !--------------------------------- 13 --------------------------------- STKDUCT 2.75 inch long, 3/4 od, 0.031 wall sameas 4a a Area m^2 Mstr 4.9171E+04 A |p| Pa 5.4860E-02 b Perim m 13a 18.954 B Ph(p) deg 3.4925E-02 c Length m 2.8751E-04 C |U| m^3/s 4.5500E-05 d WallA m^2 90.405 D Ph(U) deg -1.5734 E Htot W 2.2487 F Edot W 302.86 G TBeg K stainless Solid type 397.31 H TEnd K !--------------------------------- 14 --------------------------------- RPN xi1 and phase, mid tube 0.0000 a G or T 7.6424E-03 A m -71.451 B degr 1.5934E-02 C p1/pm p1 mag pm / p1 arg U1 arg - U1 mag 13a / w / !--------------------------------- 15 --------------------------------- STKDUCT 2.75 inch long, 3/4 od, 0.031 wall sameas 4a a Area m^2 Mstr 4.9192E+04 A |p| Pa 5.4860E-02 b Perim m 15a 18.945 B Ph(p) deg 3.4925E-02 c Length m 2.7665E-04 C |U| m^3/s 4.5500E-05 d WallA m^2 89.319 D Ph(U) deg -1.5734 E Htot W 2.2855 F Edot W 397.31 G TBeg K stainless Solid type 493.97 H TEnd K !--------------------------------- 16 --------------------------------- RPN p1 here. Also xi1. 0.0000 a G or T (4.6528E+04, 1.5971E+04) A p1 7.3537E-03 B xi1 U1 mag 15a / w / ; p1 !--------------------------------- 17 --------------------------------- JOIN Change Me 4.9192E+04 A |p| Pa 18.945 B Ph(p) deg 2.7615E-04 C |U| m^3/s 89.319 D Ph(U) deg -1.5734 E Htot W 2.2813 F Edot W 493.97 G TBeg K 493.32 H TEnd K !--------------------------------- 18 --------------------------------- RPN account for convective heat again sameas 8a a G or T -11.423 A ChngeMe H2k inp - =H2k !--------------------------------- 19 --------------------------------- STKSCREEN Change Me sameas 9a a Area m^2 4.9159E+04 A |p| Pa sameas 9b b VolPor 18.821 B Ph(p) deg sameas 9c c Length m 2.8977E-04 C |U| m^3/s sameas 9d d rh m 88.834 D Ph(U) deg sameas 9e e ksFrac -11.423 E Htot W 2.4345 F Edot W 493.32 G TBeg K stainless Solid type 524.09 H TEnd K !--------------------------------- 20 --------------------------------- RPN account for convective heat again sameas 8a a G or T -1.5734 A ChngeMe H2k inp + =H2k !--------------------------------- 21 --------------------------------- RPN calcs: want NR >~ 1000 and dp >~ rhov2/2/10 0.0000 a G or T 142.88 A NR1 1.6076E+04 B dp 7.6548 C rhov2/2 rho 19C 22a / 22b / sto sqrd * 2 / ; p1 22A - mag ; rcl 4 * 19d * rho * mu / !--------------------------------- 22 --------------------------------- TX same as above sameas 6a a Area m^2 4.9164E+04 A |p| Pa sameas 6b b GasA/A 18.810 B Ph(p) deg sameas 6c c Length m 2.8820E-04 C |U| m^3/s sameas 6d d radius m 88.773 D Ph(U) deg 1.7424 e HeatIn W G 0.16899 E Htot W 523.00 f SolidT K 2.4273 F Edot W 524.09 G GasT K ideal Solid type 525.28 H SolidT K !--------------------------------- 23 --------------------------------- JOIN Change Me 4.9164E+04 A |p| Pa 18.810 B Ph(p) deg 2.8767E-04 C |U| m^3/s 88.773 D Ph(U) deg 0.16899 E Htot W 2.4228 F Edot W 524.09 G TBeg K 523.65 H TEnd K !--------------------------------- 24 --------------------------------- RPN hot temperature target (deg C !) 250.50 a G or T =24A 250.50 A ChngeMe Tm 273.15 - !--------------------------------- 25 --------------------------------- STKDUCT very short, just to get JOIN working like I want it to. sameas 26a a Area m^2 Mstr 4.9164E+04 A |p| Pa 5.4860E-02 b Perim m 25a 18.810 B Ph(p) deg 1.0000E-06 c Length m 2.8767E-04 C |U| m^3/s 4.0000E-05 d WallA m^2 88.773 D Ph(U) deg 0.16899 E Htot W 2.4228 F Edot W 523.65 G TBeg K stainless Solid type 523.65 H TEnd K !--------------------------------- 26 --------------------------------- CONE transition to little tube sameas 6a a AreaI m^2 Mstr 4.9241E+04 A |p| Pa 5.4860E-02 b PerimI m 26a 18.734 B Ph(p) deg 3.1500E-02 c Length m 2.8309E-04 C |U| m^3/s sameas 28a d AreaF m^2 Mstr 88.484 D Ph(U) deg 9.5414E-03 e PerimF m 26d 0.16899 E Htot W 5.0000E-04 f Srough 2.4124 F Edot W ideal Solid type !--------------------------------- 27 --------------------------------- MINOR Change Me sameas 28a a Area m^2 4.8930E+04 A |p| Pa 0.2000 b K+ 17.724 B Ph(p) deg 0.8000 c K- 2.8309E-04 C |U| m^3/s 88.484 D Ph(U) deg 0.16899 E Htot W 2.2823 F Edot W !--------------------------------- 28 --------------------------------- DUCT inertance-resistance 3/16 inch od copper !3/16 refrig tubing, id = 0.128 inch nominal. !Autocad notes show 0.129 inch most probable i.d. !32 inch long !So, 8.432E-06 m2 area, 1.0294E-02 m perim, and 0.8128 m length. 7.2445E-06 a Area m^2 G 5.7836E+04 A |p| Pa 8.9532E-03 b Perim m G 0.90263 B Ph(p) deg 0.8128 c Length m 2.6793E-04 C |U| m^3/s 5.0000E-04 d Srough 88.286 D Ph(U) deg 0.16899 E Htot W ideal Solid type 0.35367 F Edot W !--------------------------------- 29 --------------------------------- RPN area and perim ratios 0.0000 a G or T 0.85917 A ChngeMe 0.86975 B ChngeMe 28b 1.0294e-2 / 28a 8.432e-6 / !--------------------------------- 30 --------------------------------- MINOR Change Me sameas 28a a Area m^2 5.7801E+04 A |p| Pa 0.9000 b K+ 5.6900E-03 B Ph(p) deg 0.2000 c K- 2.6793E-04 C |U| m^3/s 88.286 D Ph(U) deg 0.16899 E Htot W 0.23233 F Edot W !--------------------------------- 31 --------------------------------- COMPLIANCE subtract the coil's volume, and add its surface area !1/4 inch od, 34 inch = 0.8636 m long ---> 27E-6 m3, 0.017 m2 !3/16 inch od, 34 inch = 0.8636 m long --> 15E-6 m3, 0.013 m2 !3/16 inch od, 32 inch = 0.8128 m long --> 14E-6 m3, 0.012 m2 1.2000E-02 a SurfAr m^2 5.7801E+04 A |p| Pa -1.4000E-05 b Volume m^3 5.6900E-03 B Ph(p) deg 2.9256E-04 C |U| m^3/s 89.044 D Ph(U) deg 0.16899 E Htot W ideal Solid type 0.14191 F Edot W !--------------------------------- 32 --------------------------------- COMPLIANCE add the volume of the conical space at base of DUCT 2.1900E-03 a SurfAr m^2 5.7801E+04 A |p| Pa 5.2800E-06 b Volume m^3 5.6900E-03 B Ph(p) deg 2.8324E-04 C |U| m^3/s 89.128 D Ph(U) deg 0.16899 E Htot W ideal Solid type 0.12541 F Edot W !--------------------------------- 33 --------------------------------- DUCT Hot duct 2 inch od 1/8 inch wall 1.5520E-03 a Area m^2 Mstr 5.7806E+04 A |p| Pa 0.13966 b Perim m 33a 5.5478E-03 B Ph(p) deg 0.1016 c Length m 1.3018E-06 C |U| m^3/s 5.0000E-04 d Srough 60.578 D Ph(U) deg 0.16899 E Htot W ideal Solid type 1.8486E-02 F Edot W !--------------------------------- 34 --------------------------------- SURFACE top plate sameas 33a a Area m^2 5.7806E+04 A |p| Pa 5.5478E-03 B Ph(p) deg 1.1579E-06 C |U| m^3/s 78.298 D Ph(U) deg 0.16899 E Htot W ideal Solid type 6.7908E-03 F Edot W !--------------------------------- 35 --------------------------------- STKDUCT little tube to pressure transducer 4.1000E-06 a Area m^2 Mstr 5.7816E+04 A |p| Pa 7.1779E-03 b Perim m 35a 2.4607E-14 B Ph(p) deg 0.1250 c Length m 3.4604E-19 C |U| m^3/s 5.0000E-06 d WallA m^2 -64.129 D Ph(U) deg 0.16899 E Htot W 4.3649E-15 F Edot W 523.65 G TBeg K stainless Solid type 283.00 H TEnd K !--------------------------------- 36 --------------------------------- RPN Endevco at room temperature 283.00 a G or T =36A 283.00 A ChngeMe Tm !--------------------------------- 37 --------------------------------- HARDEND Some heat goes to ambient air here, so ignore Htot 0.0000 a R(1/z) =37G 5.7816E+04 A |p| Pa 0.0000 b I(1/z) =37H 2.4607E-14 B Ph(p) deg 0.0000 c Htot W 3.4604E-19 C |U| m^3/s -64.129 D Ph(U) deg 0.16899 E Htot W 4.3649E-15 F Edot W 1.4912E-15 G R(1/z) -3.0750E-15 H I(1/z) !--------------------------------- 38 --------------------------------- RPN pressure target. 5.7816E+04 a G or T =38A 5.7816E+04 A ChngeMe p1 mag !--------------------------------- 39 --------------------------------- RPN pressure target 0.0000 a G or T =39A 2.4607E-14 A ChngeMe p1 arg !--------------------------------- 40 --------------------------------- RPN Rayleigh number, based on L 0.0000 a G or T 3.7759E+07 A ChngeMe 10 13c 15c + 3 ^ * 15H 15G - * 15H 15G + 2 / / rho sqrd * cp * mu / k / !--------------------------------- 41 --------------------------------- RPN Nusselt number, based on L, a guess based on Incrop. 9.51 0.0000 a G or T 31.504 A ChngeMe 0.2 40A 0.29 ^ * !--------------------------------- 42 --------------------------------- RPN Heat flow estimate, tipped 90 degrees, no acoustics 0.0000 a G or T 1.5265 A watts k 15a * 15H 15G - * 13c 15c + / 41A * !--------------------------------- 43 --------------------------------- RPN ChangeMe 0.0000 a G or T 1.8735E-02 A tankpr 1.5921E-02 B drprat 0d pm / p1 mag pm / !--------------------------------- 44 --------------------------------- RPN ChangeMe 0.0000 a G or T 8871.3 A ChngeMe 13F 43B sqrd / ! The restart information below was generated by a previous run ! and will be used by DeltaEC the next time it opens this file. guessz 0c 2a 2b 6e 22e 28a 28b xprecn 2.6296E-03 401.51 -789.38 -1.0074E-06 -4.4733E-05 -3.0317E-13 3.8024E-08 targs 6f 24a 36a 37a 37b 38a 39a hilite 12B 13G 14A 15H 16B mstr-slave 8 1 -2 4 -2 13 -2 15 -2 25 -2 26 -9 33 -2 35 -2 ! Plot start, end, and step values. May be edited if you wish. ! Outer Loop: | Inner Loop .