pingu_98/cara
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

update var naming

Browse source
This commit is contained in:
Andrejh 2021-09-10 16:48:17 +02:00
parent 8807942603
commit 6e14f8dbb6
1 changed files with 6 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
cara/results_paper.py
View file

@ -441,17 +441,17 @@ def calculate_cunningham_slip_factor(d:int):
# Deposition factor for a breathing model # Deposition factor for a breathing model
def calculate_deposition_factor(): def calculate_deposition_factor():
exposure_mc = breathing_exposure(activity='Seated', mask='No mask') exposure_mc = breathing_exposure(activity='Heavy exercise', mask='No mask')
exposure_model = exposure_mc.build_model(size=SAMPLE_SIZE) exposure_model = exposure_mc.build_model(size=SAMPLE_SIZE)
Pp = 1000 rho_p = 1000
Uair = 1.8*10**-5 mu_air = 1.8*10**-5
FRC = 0.003 FRC = 0.003
Vt = 0.0004 Vt = 0.0004
g = 9.8 g = 9.8
BRk = exposure_model.exposed.activity.inhalation_rate BRk = exposure_model.exposed.activity.inhalation_rate
diameters = np.linspace(0.3, 30, 200) #particle diameter (multiply later by 10**(-6)) diameters = np.linspace(0.3, 100, 200) #particle diameter (multiply later by 10**(-6))
fractions = [] fractions = []
for d in diameters: for d in diameters:
@ -459,9 +459,9 @@ def calculate_deposition_factor():
cunningham_slip_factor = calculate_cunningham_slip_factor(d) cunningham_slip_factor = calculate_cunningham_slip_factor(d)
f_dep = 0.08 + 0.92 / ( f_dep = 0.08 + 0.92 / (
1 + (4.09*10**-6 * ( 1 + (4.09*10**-6 * (
(((cunningham_slip_factor*Pp*d**2*(BRk/3600))/Uair*FRC)**0.8) + ( (((cunningham_slip_factor*rho_p*d**2*(BRk/3600))/mu_air*FRC)**0.8) + (
0.01*( 0.01*(
((cunningham_slip_factor*g*Pp*d**2*FRC**(2/3))/(Uair*(BRk/3600))**0.4) * ( ((cunningham_slip_factor*g*rho_p*d**2*FRC**(2/3))/(mu_air*(BRk/3600))**0.4) * (
(Vt/FRC)**0.8 (Vt/FRC)**0.8
) )
) )
@ -474,7 +474,6 @@ def calculate_deposition_factor():
fig = plt.figure() fig = plt.figure()
ax = fig.add_subplot(1, 1, 1) ax = fig.add_subplot(1, 1, 1)
ax.plot(diameters, fractions) ax.plot(diameters, fractions)
print(diameters)
plt.ylabel( plt.ylabel(
'Deposition fraction (f$_{dep}$)', fontsize=14) 'Deposition fraction (f$_{dep}$)', fontsize=14)