474
}
475
476
double Trigidice::iterout(){
477
return(iteration);
478
}
479
480
double Trigidice::tgradout(){
481
QH2 = -RKH(WSAT)*GRADTB+H*(VB*(WSAT-WF)+ Y*VI + ((Y-1)*(WSAT-WF))*VB);
482
return(-1*theta0/eta*QH2*(FSI/RKHI()+(1-FSI)/RKH(WF)));
483
}
484
double Trigidice::heatin(){
485
return(-RKH(WSAT)*GRADTB/QH(1));
486
}
487
double Trigidice::heatout(){
488
QH2 = -RKH(WSAT)*GRADTB+H*(VB*(WSAT-WF)+ Y*VI + ((Y-1)*(WSAT-WF))*VB);
489
return(QH2/QH(1));
490
}
491
492
double Trigidice::waterin(){
493
QW1 = Y*VI + ((Y-1)*(WSAT-WF))*VB;
494
return(QW1*gammaIA*lamda/(eta*nu));
495
}
496
497
double Trigidice::spacing(){
498
return((ZI-ZN)*VI/VB * eta *1000);
499
}
500
501
double Trigidice::thickness(){
502
return(ZI * eta *1000);
503
}
504
505
void Trigidice::look(){
506
cout <<
507
iteration << ","
508
<< counter << ","
509
<< FSI << ","
510
<< (UI*gammaIA/lamda)/kP2P<< ","
511
<< VI * gammaIA*lamda/(eta*nu)* mmpdTomps << ","
512
<< VB * gammaIA*lamda/(eta*nu)* mmpdTomps << ","
513
<< GRADTB*theta0/eta << ",";
514
QW1 = Y*VI + ((Y-1)*(WSAT-WF))*VB;
515
QH1 = -RKH(WSAT)*GRADTB ;
516
QH2 = QH1+H*(VB*(WSAT-WF)+QW1);
517
//
fout << -1.0*theta0/eta*QH2*(VI*RKH(WF)+VB*RKH(0))/((VI+VB)*(RKH(0)*RKH(WF)))
518
cout <<-1.0*theta0/eta*QH2*(FSI/RKHI()+(1-FSI)/RKH(WF))
519
<< ","
520
<< QW1*gammaIA*lamda/(eta*nu) << ","
521
<< (ZI-ZN)*VI/VB * eta *1000 << ","
522
<< ZI * eta *1000
523
<< "\n";
524
}
525
27