 Timestamp:
 Jan 31, 2007, 5:02:46 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/newlanguage/eml/heat_exchangers/HeatExchangerSimplified.mso
r146 r147 30 30 PARAMETERS 31 31 outer PP as Plugin (Brief="External Physical Properties"); 32 HE as Plugin (Brief="STHE Calculations",File="heatex");33 32 outer NComp as Integer (Brief="Number of Components"); 34 33 M(NComp) as molweight (Brief="Component Mol Weight"); … … 553 552 VARIABLES 554 553 555 R as positive; 556 P as positive; 557 Pc as positive; 558 Rho as positive; 559 #Phi as positive; 554 R as positive; 555 P as positive; 556 Pc as positive; 557 Rho as positive; 558 Phi as positive; 559 lambdaN as positive; 560 lambda1 as positive; 560 561 561 562 EQUATIONS … … 564 565 Rho*(1P*R) = (1P); 565 566 566 " Capacity Ratio for LMTD Correction Fator"567 "R: Capacity Ratio for LMTD Correction Fator" 567 568 R*(Outlet.Cold.T  Inlet.Cold.T ) = (Inlet.Hot.TOutlet.Hot.T); 568 569 569 " Non  Dimensional Variable for LMTD Correction Fator"570 "P: Non  Dimensional Variable for LMTD Correction Fator" 570 571 P*(Inlet.Hot.T Inlet.Cold.T)= (Outlet.Cold.TInlet.Cold.T); 571 572 … … 575 576 "Temperature Difference at Outlet" 576 577 DTL = Outlet.Hot.T  Inlet.Cold.T; 577 578 579 switch LMTDcorrection 580 581 case "Bowmann": 582 583 lambdaN =1; 584 lambda1 =1; 585 Phi = 1; 586 578 587 if R equal 1 579 588 580 589 then 581 590 Pc*(2P)= P; 582 583 591 Fc= (sqrt(2)*Pc)/((1Pc)*ln( abs( ( 2Pc*0.585786)/( 2Pc*3.414214)))); 584 592 585 593 else 586 594 Pc = (sqrt(abs(( 1P*R)/(1P)))1)/(sqrt(abs(( 1P*R)/(1P)))R); 587 Fc =sqrt(R*R+1)*ln(abs((1Pc*R)/(1Pc)))/((1R)*ln( abs( ( 2Pc*(R+1sqrt(R*R+1)))/ ( 2Pc*(R + 1 + sqrt(R*R+1)))))); 588 589 end 590 591 592 end 593 595 Fc = sqrt(R*R+1)*ln(abs((1Pc*R)/(1Pc)))/((1R)*ln( abs( ( 2Pc*(R+1sqrt(R*R+1)))/ ( 2Pc*(R + 1 + sqrt(R*R+1)))))); 596 597 end 598 599 case "Fakeri": 600 Pc = P; 601 Phi = (sqrt(((Inlet.Hot.T  Outlet.Hot.T)*(Inlet.Hot.T Outlet.Hot.T))+((Outlet.Cold.T  Inlet.Cold.T)*(Outlet.Cold.T  Inlet.Cold.T))))/(2*((Inlet.Hot.T + Outlet.Hot.T)( Inlet.Cold.T + Outlet.Cold.T))); 602 603 if Rho equal 1 604 605 then 606 lambdaN = 1; 607 lambda1 = 1; 608 Fc = (2*Phi )/(ln(abs((1+Phi )/(1Phi )))); 609 610 else 611 612 lambdaN = (1/ln(sqrt(abs(Rho))))*((2*sqrt(abs(Rho))2)/(sqrt(abs(Rho))+1)); 613 614 lambda1 = (1/ln(abs(Rho)))*((2*Rho2)/(Rho+1)); 615 616 Fc = ((2*Phi *(lambdaN/lambda1))/(ln(abs((1+Phi *(lambdaN/lambda1))/(1Phi *(lambdaN/lambda1))))))*(1/lambdaN); 617 618 end 619 620 621 end 622 623 624 end 594 625 #===================================================================== 595 626 # NTU Method
Note: See TracChangeset
for help on using the changeset viewer.