hodnotyParam[Teplota_]:=Module[{out,hustotaVzduch,kinematickaVizkozitaVzduch,tepelnavodivost,mernaTepelnaKapacitaPriStalemTlaku},{hustotaVzduch=Interpolation[{{-100,-50,-20,0,10,20,40,60,80,100,120,140,160,180,200,250,300,350,400,500,600,800,1000},{2.012,1.561,1.376,1.275,1.230,1.188,1.112,1.046,0.986,0.934,0.886,0.843,0.804,0.769,0.736,0.666,0.608,0.559,0.517,0.451,0.399,0.325,0.274}}\[Transpose], InterpolationOrder->1];
tepelnavodivost=Interpolation[{{-100,-50,-20,0,10,20,40,60,80,100,120,140,160,180,200,250,300,350,400,500,600,800,1000},0.01*{1.62,2.04,2.26,2.37,2.45,2.52,2.65,2.80,2.93,3.07,3.2,3.33,3.44,3.57,3.70,3.98,4.29,4.57,4.81,5.4,5.82,6.69,7.62}}\[Transpose]];
mernaTepelnaKapacitaPriStalemTlaku=Interpolation[{{-100,-50,-20,0,10,20,40,60,80,100,120,140,160,180,200,250,300,350,400,500,600,800,1000},{1022,1013,1009,1005,1005,1010,1013,1017,1020,1022,1024,1027,1030,1034,1037,1044,1049,1055,1059,1072,1089,1114,1139}}\[Transpose]];
kinematickaVizkozitaVzduch=Interpolation[{{-100,-50,-20,0,10,20,40,60,80,100,120,140,160,180,200,250,300,350,400,500,600,800,1000},0.000001*{5.86,9.35,11.85,13.49,14.47,15.32,17.27,19.22,21.3,23.34,25.73,27.88,29.98,32.51,35.19,42.04,48.85,56.35,63.83,80.27,98.25,136.9,180.7}}\[Transpose]];
 out={hustotaVzduch[Teplota],kinematickaVizkozitaVzduch[Teplota],tepelnavodivost[Teplota],mernaTepelnaKapacitaPriStalemTlaku[Teplota]}};
out
]

hodnotyParam[50.]

(1/(50+273.))

Pr=0.71;
Gr=(1/(50+273.))*60*9.81*((1^3)/(((18.24*(10^(-6))))^2));

Pr*Gr;
c=0.135;
n=(1/3);

Nu=c*(((Pr*Gr))^n)

α=((Nu*0.0272)/1)

P=α*1*50

α[TP_,To_,Dchar_]:=Module[{Tstr,β,ΔT,ρ,Cp,λ,visk,Pr,Gr,fceNC,αin},
Tstr=((TP+To)/2);
 β=(1/(Tstr+273.15));
ΔT=Abs[TP-To];
ρ=hodnotyParam[Tstr][[1]];
Cp=hodnotyParam[Tstr][[4]];
λ=hodnotyParam[Tstr][[3]];
visk=hodnotyParam[Tstr][[2]];
Pr=((ρ*Cp*visk)/λ);
Gr=((β*ΔT*9.81*(Dchar^3))/(visk^2));
fceNC[souc_]:=1.18*(souc^(1/8))/;souc≤5*(10^2);
fceNC[souc_]:=0.54*(souc^(1/4))/;5*(10^2)≤souc<2*(10^7);
fceNC[souc_]:=0.135*(souc^(1/3))/;souc≥2*(10^7);
αin=((fceNC[Pr*Gr]*λ)/Dchar)
];
λ=0.06;
l=1;
ϵ=1;
To=20;
σ=5.67*(10^(-8));
D1=0.1;
D2=0.3;
T1=80;
Riz=(1/(2*Pi*λ*l))*Log[(D2/D1)];
Qsalani=ϵ*Pi*D2*l*σ*((((TP+273))^4)-(((To+273))^4));
Qkonvekce=α[TP,To,D2]*Pi*D2*l*(TP-To);
rce=Qsalani+Qkonvekce==((T1-TP)/Riz);
FindRoot[rce,{TP,100}]
Qsalani+Qkonvekce/.FindRoot[rce,{TP,100}]
α[TP,To,D2]/.FindRoot[rce,{TP,100}]