domain -7.00 9.00 -8.50 7.50 grid 7 2.0e-05 10000 # # 50000 steps = 1 time unit # 2500 steps = 0.05 time units = 1 cycle # # Viscosity = 0.0002, Density = 0.001, units in grams, cm, seconds # # Rescaling: 1 IBIS UNIT LENGTH = airway diameter = 0.3385 cm # viscosity: -m viscosity*diameter # density: -r density*diameter^3 # fluid -m 0.000067 -r 0.000037 # # visual -p binary -b lung -h 512 -w 512 -n 100 graph -n 100 -p line -b lung -d 1 # # -r indicates magnitude of random particle movement # particle -P bolus -r 0.0 # # # spring r -b 1.0e+0 -s 1.0e+1 -t 1.0e+4 spring b -b 1.0e+0 -s 1.0e+1 -t 1.0e+4 # # speaker # line -1.00 4.00 1.00 4.00 -n 40 -A 1 line -1.00 4.10 1.00 4.10 -n 40 -A 2 blink 1 2 # # left branch # polygon -N 11 first -1.50 -3.70 bcurv f 2 -n 60 -0.75 -3.50 -0.50 -2.00 ray -0.50 1.00 -n 60 ray -1.00 3.00 -n 40 ray -1.00 4.00 -n 20 end polygon -N 12 first -1.50 -3.60 bcurv f 2 -n 60 -0.85 -3.50 -0.60 -2.00 ray -0.60 1.00 -n 60 ray -1.10 3.00 -n 40 ray -1.10 4.00 -n 20 end blink 11 12 # # right branch # polygon -N 21 first 1.00 4.00 ray 1.00 3.00 -n 20 ray 0.50 1.00 -n 40 ray 0.50 -2.00 -n 60 bcurv f 3 -n 80 1.50 -3.25 2.50 -3.50 3.20 -3.15 end polygon -N 22 first 1.10 4.00 ray 1.10 3.00 -n 20 ray 0.60 1.00 -n 40 ray 0.60 -2.00 -n 60 bcurv f 3 -n 80 1.50 -3.15 2.50 -3.40 3.20 -3.05 end blink 21 22 # # left divider 1 # polygon -N 31 first -1.50 -4.10 bcurv f 2 -n 20 -1.00 -4.10 -0.80 -3.80 bcurv f 3 -n 60 -1.25 -4.80 -1.50 -5.20 -1.70 -5.40 end polygon -N 32 first -1.50 -4.20 bcurv f 2 -n 20 -1.00 -4.20 -0.90 -3.90 bcurv f 3 -n 60 -1.25 -4.70 -1.50 -5.10 -1.70 -5.30 end blink 31 32 # # left divider 2 # polygon -N 41 first -1.60 -5.65 bcurv f 2 -n 20 -1.40 -5.45 -1.15 -5.35 bcurv f 2 -n 20 -1.20 -5.50 -1.20 -6.00 end polygon -N 42 first -1.60 -5.75 bcurv f 2 -n 20 -1.40 -5.55 -1.25 -5.45 bcurv f 2 -n 20 -1.20 -5.60 -1.30 -6.00 end blink 41 42 # # bottom branch # polygon -N 51 first -0.85 -6.00 bcurv f 3 -n 80 -0.85 -5.50 -0.30 -4.00 0.05 -3.10 bcurv f 4 -n 80 0.75 -3.55 1.50 -3.75 2.10 -3.90 2.70 -4.45 end polygon -N 52 first -0.75 -6.00 bcurv f 3 -n 80 -0.75 -5.50 -0.20 -4.00 0.05 -3.20 bcurv f 4 -n 80 0.75 -3.65 1.40 -3.75 2.00 -3.90 2.60 -4.45 end blink 51 52 # # right divider # polygon -N 61 first 3.00 -4.15 ray 2.75 -3.80 -n 20 ray 3.35 -3.50 -n 20 end polygon -N 62 first 3.00 -4.05 ray 2.85 -3.80 -n 20 ray 3.35 -3.60 -n 20 end blink 61 62