free format * ************* * main data * ************* * * numtcr ieos inopt nmat 5 0 1 0 TEST PROBLEM FOR PUMP This is a simple closed loop with a break acting as a pressurizer Flow is driven steady state against a fixed loss coefficient by a pump with constant motor torque * ***************** * namelist data * ***************** * &inopts iadded=50, nhtstr=0, nsend=5000, idiag=0, noair=1, nfrc1=2, ikfac=1, dtstrt=.01 &end * * dstep timet -1 0.000e+00 * stdyst transi ncomp njun ipak 0 1 4 4 1 * epso epss 1.0000e-04 1.0000e-04 * oitmax sitmax isolut ncontr 10 10 0 0 * ntsv ntcb ntcf ntrp ntcp 2 7 10 1 0 * ************************* * component-number data * ************************* * * iorder* 1 20 30 60 e * ************************** * control-parameter data * ************************** * * * signal variables * idsv isvn ilcn icn1 icn2 1 0 0 0 0 *time 110 32 60 2 0 * pump mass flow * * control blocks * * * Pump Torque Table * idcb icbn icb1 icb2 icb3 -1040 101 1 3 0 * gain xmin xmax cbcon1 cbcon2 1.0 0. 1.e12 0. 0. 0. 0.0 s 0.01 23.8 s 10000. 23.8 e * * Pump Mass Flow Set Point (5 kg/sec) * idcb icbn icb1 icb2 icb3 -1045 9 0 0 0 * gain xmin xmax cbcon1 cbcon2 1.00 0. 1.e12 5.0 0. * * Time to Start active control * idcb icbn icb1 icb2 icb3 -1047 9 0 0 0 * gain xmin xmax cbcon1 cbcon2 1.00 0. 1.e12 200.0 0. * * Switch to Check time for Controller * idcb icbn icb1 icb2 icb3 -1048 20 1 -1047 0 * gain xmin xmax cbcon1 cbcon2 1.00 0. 1.e12 0.0 0. * * Set to proper set point at control start time * idcb icbn icb1 icb2 icb3 -1049 22 -1045 110 -1048 * gain xmin xmax cbcon1 cbcon2 1.00 0. 1.e12 0.0 0. * * PI Controller on Pump Torque * idcb icbn icb1 icb2 icb3 -1050 200 110 -1049 0 * Note: change first number in next active card to vary gain * gain xmin xmax cbcon1 cbcon2 100.0 0. 30. 0.0 200. * * First number in the next active card gives the integration time scale * * cbtau cbdt cbwt 0.1 .050 0. * * Set to proper Speed Driver at control start time * idcb icbn icb1 icb2 icb3 -1060 22 -1050 -1040 -1048 * gain xmin xmax cbcon1 cbcon2 1.00 0. 1.e12 0.0 0. * * * trips * * Set a trip that is always off to force pump speed from the control blook * listed in pump variable npmpsd * * ntse ntct ntsf ntdp ntsd 0 0 0 0 0 * idtp isrt iset itst idsg -2 2 0 1 1 * setp(1) setp(2) -1.0000e+00 1.0000e+08 * dtsp(1) dtsp(2) 0.0000e+00 0.0000e+00 * ifsp(1) ifsp(2) 0 0 * ****************** * component data * ****************** * ******* type num id ctitle break 1 1 liquid source * jun1 ibty isat ioff 1 0 0 0 * dxin volin alpin tin pin 1.0000e+00 5.067e-03 0.0000e+00 3.0000e+02 2.0000e+05 * pain concin rbmx poff belv 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 * ******* type num id ctitle pipe 20 20 Flow Resistance Face 4 * ncells nodes jun1 jun2 epsw 10 0 3 4 0.0000e+00 * ichf iconc iacc ipow 1 0 0 0 * iqp3tr iqp3sv nqp3tb nqp3sv nqp3rf * 11 0 0 0 0* * radin th houtl houtv toutl .0127 .002 0.0000e+00 0.0000e+00 3.7300e+02 * toutv 3.7300e+02 * qp3in qp3off rqp3mx qp3scl * 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 * * * dx * f 1.0 e * vol * f 5.0671e-3 e * fa * f 5.0671e-4 e * fric * r02 0.0 17. f 0. e *rv fric * r02 0.0 17. f 0. e * grav * f 0. e * hd * f .0254 e * nff * f -1 e * alp * f 0.0000e+00e * vl * f 0.0000e+00e * vv * f 0.0000e+00e * tl * f 300.00e+00e * tv * f 300.00e+00e * p * f 2.0000e+05e * pa * f 0.0000e+00e *qppp *matid * tw * ******* type num id ctitle pump 60 60 pump 60 * ncells nodes jun1 jun2 epsw 2 0 2 3 0.0000e+00 * ichf iconc ipmpty irp ipm 1 0 3 1 1 * ipmptr ipmpsv npmptb npmpsv npmprf -2 -1060 0 0 0 * npmpmt 2 * iqp3tr iqp3sv nqp3tb nqp3sv nqp3rf * 0 0 0 0 0 * radin th houtl houtv toutl .02540e+00 0.2000e-02 0.0000e+00 0.0000e+00 3.7300e+02 * toutv effmi 3.7300e+02 0.2e+00 * tfr0 tfr1 tfr2 tfr3 tfrb 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 * tfrl0 tfrl1 tfrl2 tfrl3 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 * rhead rtork rflow rrho romega 1.0000e+03 2.5000e+01 5.0000e-03 9.9700e+02 2.0000e+02 * omegan omgoff romgmx omgscl npmpsd 0.0000e+02 -1.0000e+00 2.0000e+03 1.0000e+00 0 * qp3in qp3off rqp3mx qp3scl * 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00* * option 1 * dx * f .1 e * vol * f 5.0671e-5 e * fa * f 5.0671e-4 e * fric * f 0. e * rv fri* f 0. e * grav * f 0. e * hd * f .0254 e * nff * f -1e * alp * f 0.0000e+00e * vl * f 0.0000e+00e * vv * f 0.0000e+00e * tl * f 3.0000e+02e * tv * f 3.0000e+02e * p * f 2.0000e+05e * pa * f 0.0000e+00e *qppp* *matid* * tw * * Pump Torque table (just translates control block output directly to torque) 0.0 0.0 s 1000.0 1000.0 e * ******* type num id ctitle tee 30 30 loop 1 - cold leg section 1 * jcell nodes ichf cost epsw 2 0 1 0.0000e-01 0.0000e+00 * iconc1 ncell1 jun1 jun2 ipow1 0 3 4 2 0 * iqptr1 iqpsv1 nqptb1 nqpsv1 nqprf1 * 0 11 0 0 0* * radin1 th1 houtl1 houtv1 toutl1 1.2700e-02 .2000e-02 0.0000e+00 0.0000e+00 3.7300e+02 * toutv1 3.7300e+02 * qpin1 qpoff1 rqpmx1 qpscl1 * 0.0000e+00 0.0000e+00 1.0000e+06 0.0000e+00* * iconc2 ncell2 jun3 ipow2 0 2 1 0 * iqptr2 iqpsv2 nqptb2 nqpsv2 nqprf2 * 0 11 0 0 0* * radin2 th2 houtl2 houtv2 toutl2 1.2700e-02 .2000e-02 0.0000e+00 0.0000e+00 3.7300e+02 * toutv2 3.7300e+02 * qpin2 qpoff2 rqpmx2 qpscl2 * 0.0000e+00 0.0000e+00 1.0000e+06 0.0000e+00 * * * * dx * f 1.0 e * vol * f 5.0671e-4 e * fa * 5.0671e-4 5.0671e-4 5.0671e-4 5.0671e-4e * fric * f 0. e *rv fric * f 0. e * grav * f 0. e * hd * f .0254 e * nff * -1 1 1 -1 e * alp * f 0.0000e+00e * vl * f 0.0000e+00e * vv * f 0.0000e+00e * tl * f 300.00e+00e * tv * f 300.00e+00e * p * f 2.0000e+05e * pa * f 0.0000e+00e *qppp* *matid* * tw * * * dx * f 1.0 e * vol * f 5.0671e-3 e * fa * f 5.0671e-3 e * fric * f 0. e *rv fric * f 0. e * grav * f 0. e * hd * f .0254 e * nff * -1 f 1 e * alp * f 0.0000e+00e * vl * f 0.0000e+00e * vv * f 0.0000e+00e * tl * f 300.00e+00e * tv * f 300.00e+00e * p * f 2.0000e+05e * pa * f 0.0000e+00e *qppp* *matid* * tw * * end * ****************** * time-step data * ****************** * * dtmin dtmax tend rtwfp 1.0000e-04 0.5000e-01 1.0000e+01 1.0000e+02 * edint gfint dmpint sedint 2.0000e+00 5.0000e-02 1.0000e+02 1.0000e+02 * ****************** * time-step data * ****************** * * endflag -1.0000e+00