Skip to main content

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);

for

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