|
@@ -1,7 +1,7 @@
|
|
|
# Hot Banana
|
|
|
#
|
|
|
-# I heard that bananas are radioactive. If they are radioactive, then they
|
|
|
-# radiate energy. How many bananas would you need to power a house?
|
|
|
+# I heard that bananas are radioactive. If they are radioactive, then
|
|
|
+# they radiate energy. How many bananas would you need to power a house?
|
|
|
#
|
|
|
# https://what-if.xkcd.com/158/
|
|
|
|
|
@@ -17,7 +17,9 @@ let occurence_40K = 0.0117 percent
|
|
|
# We can now compute the radioactivity of natural potassium
|
|
|
|
|
|
let decay_rate_40K: Activity = ln(2) / halflife_40K
|
|
|
-let radioactivity_potassium: Activity / Mass = N_A * occurence_40K * decay_rate_40K / molar_mass_40K -> Bq / g
|
|
|
+
|
|
|
+let radioactivity_potassium: Activity / Mass =
|
|
|
+ N_A × occurence_40K × decay_rate_40K / molar_mass_40K -> Bq / g
|
|
|
|
|
|
print(radioactivity_potassium)
|
|
|
|
|
@@ -30,26 +32,34 @@ unit banana
|
|
|
|
|
|
let potassium_per_banana = 451 mg / banana
|
|
|
|
|
|
-let radioactivity_banana: Activity / Banana = potassium_per_banana * radioactivity_potassium -> Bq / banana
|
|
|
+let radioactivity_banana: Activity / Banana =
|
|
|
+ potassium_per_banana × radioactivity_potassium -> Bq / banana
|
|
|
+
|
|
|
print(radioactivity_banana)
|
|
|
|
|
|
# A single 40-K decay releases an energy of
|
|
|
# (https://commons.wikimedia.org/wiki/File:Potassium-40-decay-scheme.svg)
|
|
|
|
|
|
-let energy_per_decay: Energy = 11 percent * 1.5 MeV + 89 percent * 1.3 MeV
|
|
|
+let energy_per_decay: Energy = 11 percent × 1.5 MeV + 89 percent × 1.3 MeV
|
|
|
|
|
|
# Finally: how many bananas do we need to power a single household?
|
|
|
|
|
|
-let power_per_banana: Power / Banana = radioactivity_banana * energy_per_decay -> pW / banana
|
|
|
+let power_per_banana: Power / Banana =
|
|
|
+ radioactivity_banana × energy_per_decay -> pW / banana
|
|
|
+
|
|
|
print(power_per_banana)
|
|
|
|
|
|
unit household
|
|
|
-let power_consumption_household: Power / Household = 3000 kWh per household per year
|
|
|
|
|
|
-let bananas_per_household = power_consumption_household / power_per_banana -> bananas / household
|
|
|
+let power_consumption_household: Power / Household =
|
|
|
+ 3000 kWh per household per year
|
|
|
+
|
|
|
+let bananas_per_household =
|
|
|
+ power_consumption_household / power_per_banana -> bananas / household
|
|
|
+
|
|
|
print(bananas_per_household)
|
|
|
|
|
|
-# TODO: https://what-if.xkcd.com/158/ says this number should be around 300
|
|
|
-# quadrillion, but we only get 0.1 quadrillion. 300 quadrillion times "a couple
|
|
|
-# of picowatt" would be an average power consumption of at least 300 kW /
|
|
|
-# household, which seems … excessive.
|
|
|
+# TODO: https://what-if.xkcd.com/158/ says this number should be around
|
|
|
+# 300 quadrillion, but we only get 0.1 quadrillion. 300 quadrillion
|
|
|
+# times "a couple of picowatt" would be an average power consumption of
|
|
|
+# at least 300 kW / household, which seems … excessive.
|