Matlab Codes
Volume of the Aveolus
Va_Alveolus(Radius, Internal Height)
y = 1*pi*height_aveolus_internal*height_aveolus_internal*(3*radius_aveolus-height_aveolus_internal)/3;
Surface Area of the Aveolus
Sa_Alveolus(Radius, Internal Height)
y = 2*pi*radius_aveolus*height_aveolus_internal;
Radius of the Aveolus
r_Alveolus(length of opening radius, Internal Height)
y = ((length_of_opening_radius^2)+(interal_height_alveolus^2))/(2*interal_height_alveolus);
Total volume of Air in the lungs
Vair_Alveolus(Volume of 1 alveolus)
y = Vad + Na*Va_Alveolus;
Maximum Volume of Air in the Alveoulus
Vamax_Alveolus (Vtlc_Alveolus)
y = (Vtlc_Alveolus-Vad)/Na;
Pressure of the cavity wall
Pcw_Alveolus (Vair_Alveolus, Vtlc_Alveolus)
y = 0.071-((log((0.95/((Vair_Alveolus/Vtlc_Alveolus)-0.22))-1))*0.58);
P Hydro
Defined as a constant
Surface tension (Gamma)
gamma_Alveolus (Surface_Area)
So_new(1) = 1e-5;
count = length(Sa);
i = 1:count
Sr(i) = Sa(i) / So_new(i);
if
Sr(i) > 1
So_new(i+1) = Sa(i);
elseif
Sr(i) >= 0.617
So_new(i+1) = So_new(i);
else
So_new(i+1) = (Sa(i) / 0.617);
end
y(i) = (98.9*Sr(i)*Sr(i))-(89*Sr(i))+18,3;
end
end
