From c9cd0ed5facbcfa74cb8f77da11fc6af1abf47cf Mon Sep 17 00:00:00 2001 From: markus Date: Tue, 23 Feb 2021 11:43:40 +0100 Subject: [PATCH] add classroom_model_with_hepa --- cara/montecarlo.py | 33 +++++++++++++++++++++++++++++++++ 1 file changed, 33 insertions(+) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index 6318a2f1..e14f1ec2 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -1085,6 +1085,39 @@ classroom_model = MCExposureModel( ) ) +classroom_model_with_hepa = MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.MultipleVentilation( + ventilations=( + models.SlidingWindow( + active=models.PeriodicInterval(period=120, duration=10), + inside_temp=models.PiecewiseConstant((0, 24), (293,)), + outside_temp=models.PiecewiseConstant((0, 24), (283,)), + window_height=1.6, opening_length=0.6, + ), + models.AirChange(active=models.PeriodicInterval(period=120, duration=120), + air_exch=5) + ) + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 2), (2.5, 4), (5, 7), (7.5, 9))), + masked=False, + virus=MCVirus(halflife=1.1, qID=60), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 2), (2.5, 4), (5, 7), (7.5, 9))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +) + shared_office_model = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=50),