gusucode.com > 热力学传导matlab源码程序 > Code\HTC.m
function h = HTC(Hsub,H,x) % Three Regime Model for h mfp = 60e-9; % Mean Free Path in Air alpha = 1; % Geometric Correction Factor A = 0.9; % Thermal Accomodation Coeff gamma = 1.4; % Heat Capacity Ratio Pr = 0.707; % Prantl No kb = 8.314; % Boltzmann Const mm = 28.96e-3; % Molecular Mass of Air Tav = 313; % Average Temp for calculation of v rho = 1.12; % Density of air at Tav cp = 1008; % Specific Heat (Const P) of air at Tav kg = 0.03; % Thermal Conductivity of air C = rho*cp; % Volumetric Specific Heat v = sqrt(3*kb*Tav/mm); % RMS velocity of air molecules f = 2*(2 - A)*gamma/(A*(gamma + 1)*Pr); Htot = H - Hsub; h = [ ]; for i = 1:length(x) if (Htot(i)/mfp) <= 1 h(i) = alpha*C*v/( 3*(1 + 2*f) ); elseif (Htot(i)/mfp) < 100 h(i) = ( alpha*kg/Htot(i) )/( 1 + 2*f*mfp/Htot(i) ); else h(i) = alpha*kg/Htot(i); end end