(*-----zadané hodnoty---------------------------------------------------------------------------------------------------*)
delta1=55;      (*°C*)
delta2=30 ;    (*°C*)

(*pro hi*)
Ms=7;                 (* hmotnostní průtok tekutiny v kg/s*)
μs=0.17;       (*viskozita boční stěny při dané teplotě v kg/(m*s)*)
μw=0.14;         (*viskozita stěn trubek při dané teplotě v(kg/m*s)*)
Cps=3400;         (*měrná tepelná kapacita media kJ/kg*)
ks=16.27;        (*tepelná vodivost materiálu w*(m^(-1))*(K^(-1))*)
(*pro Jc*)
Bc=60;                (* velikost přepážek v % Ds *) 
Ds=0.9;              (*vnější průměr výměníku v m*)
d=0.08;              (*tzv, hydraulický průměr m *)
(*Pro J1*)
Layout=90;         (*dispozice 30°, 45°, nebo 90°*)
Dotl=0.85;       (*vnější průměr nejzašší hranice trubic v m obr 5.26 *)
Ltb=0.8;            (*Vůle mezi nejzazší vzdáleností trubek a otvorem přepážky, jeho hodnota je 0,4 nebo 0,8 *)
Nt=70;                (*Celkový počet trubek*)
Ltp=0.23;          (*rozteč trubic v m podle obrázku 5.2*)    
Lbc=0.41;          (*Lbc je střední hodnota roztečí přepážek podle obr.11)*)         
(*Pro Jb*)
Lpl=0;
Lpp=Ltp;                    (*Lpp je efektivní vzdálenost 2 řad ve směru proudění(podle obr. 9)*)
Nss=4;                       (*počet těsnících pásů(párů) v jedné rozteči*)
(*Pro Js*)
Lta=0.25;                 (*vzdálenost mezi přepážkami v m(Podle obr.12)*)
Lbi=0.40;                 (*rozteč přepážek na vstupní straně v m(podle obr.11)*)
Lb0=0.40;                 (*rozteč přepážek na výstupní straně v m (podle obr.11)*)
(*Pro Jr*)





(*-----Výpočet-----------------------------------------------------------------------------------------------------------*)


Res=((d*Gs)/(μs));    (*Reynoldsovo číslo, d je hydraulický průměr, μs je viskozita boční stěny při dané teplotě, Gs je hmotnostní rychlost*)
If[1<=Res<100,ji=1.73*(Res^(-0.694))];         (*Ideální Colburnův j faktor pro různá reynoldsova čísla *)
If[100<=Res<1000,ji=0.717*(Res^(-0.574))];
If[1000<=Res,ji=0.036*(Res^(-0.346))];
                                    
                                                       
Gs=(Ms/Sm);                                                                 (*Gs je hmotnostní rychlost v kg/((m^2)*s), Sm Průtoková plocha v ose svazku, Ms je hmotnostní průtok tekutiny *)
Sm=Lbc*(Lbb+(Dctl/Ltpeff)*(Ltp-d));           (*Sm Průtoková plocha v ose svazku, Lbc je střední hodnota roztečí přepážek , d je hydraulický průměr, Lbb je vzdálenost stěny od trubek
                                      Dctl je nejzazší vzdálenost středů trubek od středu výměníku*)

Lbb=Ds-Dotl;                                                    (*Lbb je vzdálenost stěny od trubek, Ds je vnější průměr výměníku, Dotl je vnější průměr nejzašší hranice trubic*)

If[ 44<Layout<46,Ltpeff=0.707* Ltp,Ltpeff=Ltp]; (*Layout je dispozice, Ltp je vzdálenost průměrů trubek, Ltpeff je efektivní vzdálenost trubek*)

ϕs=(((((μs)/(μw))))^0.14);                                                    (*Korekční faktor viskozity ((ϕs)^n), μs je viskozita boční stěny při dané teplotě,μw je viskozita stěn trubek při dané teplotě*)

Prs=((μs*Cps)/ks);                                                    (*Prandtlovo číslo, ks je tepelná vodivost materiálu, Cps je měrná tepelná kapacita media, μs je viskozita boční stěny při dané teplotě*)

hi=((ji*Cps*Gs*ϕs)/(Prs^(2/3)));                                       (*Koeficient ideálního přestupu tepla, ji je Ideální Colburnův faktor,Cps je měrná tepelná kapacita media,Gs je hmotnostní rychlost *)



Fw=((θctl)/(2*Pi))-((Sin[θctl])/(2*Pi));                                   (*Fw je poměr trubic v přepážkovém okně, θctl je vnitřní úhel svíraný přepážkami*)
θctl=2*ArcCos[(Ds/Dctl)*(1-((2*Bc)/100))];       (*θctl je vnitřní úhel svíraný přepážkami, Ds je vnější průměr výměníku, Bc je velikost přepážek v % Ds, 
	                                  Dctl je nejzazší vzdálenost středů trubek od středu výměníku*)

Dctl=Dotl-d;                                                   (*Dctl je nejzazší vzdálenost středů trubek od středu výměníku, Dotl je nejzazší vzdálenost trubek, d je hydraulický průměr*)
Fc=1-2*Fw;                                                      (*Fc je poměr trubic v proudu kapaliny, Fw je poměr trubic v přepážkovém okně*)

Jc=0.55+0.72*Fc;                                         (*Jc je korekční faktor přepážkových oken pro rozsah 15-45%, měl by se pohybovat okolo 1.0, Fc je poměr trubic v proudu kapaliny*)




θds=2*ArcCos[1-((2*Bc)/100)];                                 (*θds je úhel svíraný přepážkami,Bc je velikost přepážek v % Ds*)
Lsb=3.1+0.0004*Ds;                                           (*Ds je vnější průměr výměníku*)
Ssb = Pi*Ds*(Lsb/2)*(((2*Pi-θds )/(2*Pi)));                  (*Ds je vnější průměr výměníku, θds je úhel svíraný přepážkami*)
Stb=(4/Pi)*((((d+Ltb))^2)-(d^2))*Nt-(1-Fw);      (*Fw je poměr trubic v přepážkovém okně*)
(*Dctl= Dotl-d*)                                                   (*Dctl je nejzazší vzdálenost středů trubek od středu výměníku, Dotl je nejzazší vzdálenost trubek, d je hydraulický průměr*)
(*Lbb=Ds-Dotl*)                                                      (*Lbb je vzdálenost stěny od trubek, Ds je vnější průměr výměníku, Dotl je vnější průměr nejzašší hranice trubic*)
(*Sm=Lbc*(Lbb+(Dctl/Ltpeff)*(Ltp-d)) *)                 (*Sm Průtoková plocha v ose svazku *)
rs=(Ssb/(Ssb+Stb));                                                             (*Ssb je oblast úniku mezi pláštěm a přepážkou a Stb je oblast úniku mezi trubicemi a přepážkou *)
r1m=((Ssb+Stb)/Sm);                                                         (*Ssb je oblast úniku mezi pláštěm a přepážkou a Stb je oblast úniku mezi trubicemi a přepážkou,Sm Průtoková plocha v ose svazku *)
J1=0.44*(1-rs)+(1-0.44*(1-rs))*(E^(-2.2*r1m));    (*J1 je korekční faktor pro přepážkový únik, rs a r1m jsou poměrové veličiny rozměrů výměníku *)



If[Res<=100,Cbh=1.25] ;
If[100<Res,Cbh=1.35] ;
 Fsbp=(Sb/Sm) ;                                                                     (*Sb je prostor mezi pláštěm a svazkem trubek v rámci jedné přepážky a Sm Průtoková plocha v ose svazku *)
Sb=Lbc*(Ds-Dotl+Lpl) ;                                        (*Lbc je střední hodnota roztečí přepážek, Ds je vnější průměr výměníku *)            
rss=(Nss/Ntcc);                                                                     (*Nss je počet těsnících pásů(párů) v jedné rozteči,Ntcc je efektivní počet řad trubic v příčném proudu *)
Ntcc=-(Ds/Lpp)*(1-((2*Bc)/100));                                      (*Ds je vnější průměr výměníku, Lpp je efektivní vzdálenost 2 řad ve směru proudění, Bc je velikost přepážek v % Ds,
                                        Ntcc je efektivní počet řad trubic v příčném proudu *)
Jb=Exp[-Cbh*Fsbp*(1+((Abs[(2*rss)])^(1/3)))];                  (* Fsbp je poměr mezi částí obtékajícího proudu ku proudu hlavnímu, Cbh je rovno 1,25 pro laminární proudění a 
										Reynoldsovo číslo menší než 100,pro turbulentní a přechodné proudění je jeho Cbh rovno 1,35*)


If[Res<=100,n=1,n=0.2];            (*n=1 pro laminární proudění,Res≤100 a n=0.2 pro turbulentní proudění*)
Nb=(Lti/Lbc)-1;                                                          (*Nb je celkový počet přepážek, Lbc je vzdálenost přepážek*)
Lti=Lb0-Lts;                   (*vzdálenost mezi přepážkami v m podle obr 5.58*)
Lts=0.1*Ds;                       (*je tloušťka stěny a její limit je 25 mm, Ds je vnější průměr výměníku*)
Lbc=(Lta/(Nb+1));                         (*Lta je délka trubky mezi stěnami, Nb je celkový počet přepážek*)
Li=(Lbi/Lbc) ;                                (*Lbi je rozteč přepážek na vstupní straně, Lbc je střední hodnota roztečí přepážek *)
L0=(Lb0/Lbc);                              (*Lbc je střední hodnota roztečí přepážek podle, Lb0 je rozteč přepážek na výstupní straně *) 
Js=(((Nb-1)+(((Li))^(1-n))+(L0^(1-n)))/((Nb-1)+(((Li-1))^(1-n))+(L0^(1-n))));            (*Js je korekční faktor nestejnosti vzdáleností přepážek obvykle mezi 0,85 a 1.0*)



Ntcw=(0.8/Lpp)-(((Ds*Bc)/100)-((Ds-Dctl)/2));           (*Ntcw je efektivní počet řad trubek v přepážkovém okně, Ds je vnější průměr výměníku,Lpp je efektivní vzdálenost 2 řad ve směru proudění,
                                    Ds je vnější průměr výměníku, Bc je velikost přepážek v % Ds, Dctl je nejzazší vzdálenost středů trubek od středu výměníku *)
(*Ntcc=-(Ds/Lpp)*(1-((2*Bc)/100)) *)                                     (*Ds je vnější průměr výměníku, Lpp je efektivní vzdálenost 2 řad ve směru proudění, Bc je velikost přepážek v % Ds,
                                        Ntcc je efektivní počet řad trubic v příčném proudu *)
Nc=(Ntcc+Ntcw)*(Nb+1);     (*Nc je celkový počet řad trubek vevýměníku, Nb je celkový počet přepážek, Ntcc je efektivní počet řad trubic v příčném proudu,
							Ntcw je efektivní počet řad trubek v přepážkovém okně*)

Jr=(1.51/(Nc^0.18)) ;                                                           (*Korekční faktor teplotního gradientu laminárního proudění*)
                          
             
hs=hi*Jc*J1*Jb*Js*Jr  (*Koeficient přestupu tepla, 
hi je koeficient ideálního přestupu tepla,
Jc je korekční faktor přepážkových oken
J1 je korekční faktor pro efekt úniku způsobeny přepážkami
Jb je korekční faktor svazků tvořící obtokový proud
Js je korekční faktor nestejnosti vzdáleností přepážkových oken
Jr je korekční faktor pro nepříznivý teplotní spád daný laminárním prouděním
*)

Q=hs*((delta1-delta2)/(Log[(delta1/delta2)]))     (*Bilanční rovnice, hs je koeficient přestupu tepla, delta1 a 2 jsou rozdíly teplot na primární a sekundární straně výměníku*)