/** Program to calculate vaporization of a block of ice from the thermal radiation of a nuclear blast. See: https://upload.wikimedia.org/wikipedia/commons/4/43/Nuclear_Fireball_Radius_and_Temperature.svg https://de.wikipedia.org/wiki/Kernwaffenexplosion#Feuerball_und_Druckwelle http://www.fourmilab.ch/etexts/www/effects/ */ emissivity = .98 // Good IR emissivity for common materials fireballDiameter = 480 m fireballRadius = fireballDiameter/2 fireballArea = 4 pi fireballRadius^2 detonationHeight = 500 m T = 8000 K targetArea = 1 m^2 targetDist = detonationHeight - fireballRadius targetDensity = water vaporizationEnergy = (80 + 540) calorie/g + (104 degC 1 calorie/g/degC) albedo = .3 t = .2 s totalPower = sigma emissivity T^4 fireballArea powerFraction = targetArea / (4 pi targetDist^2) appliedPower = totalPower powerFraction (1 - albedo) vaporizedMass = appliedPower t / vaporizationEnergy vaporizationDepth = vaporizedMass / targetArea / targetDensity println[vaporizationDepth -> "cm"]