Monte-Carlo simulation of radon equilibrium under varying conditions

Research output: Contribution to ConferencePosterpeer-review


Preliminary results of a theoretical investigation into radon equilibrium under a variety of constraints will be presented, e.g. response to metallic daughter plate-out, response to modulation of radon ingress and phase relationship of alpha activity to cyclically modulated radon ingress. At the core of this investigation is the development of Monte-Carlo simulations of the radon (Rn-222) and thoron (Rn-220) decay chains. Currently, these simulations assume an hypothetical closed cubic metre of atmosphere and the simulations are time-stepped at constant interval. At each time-step, there is a radon or thoron ingress and each nuclei in the decay chains decays probabilistically, with nuclei and activity being aggregated and tabulated at the end of each time-step. Currently, the simulations have been coded in two open-source interpreted mathematical software packages; Scilab and Yorick ( This has required that the decay chains are considered sequentially (in reverse order) within each time-step. However, once the simulations are finalised, these will be compiled to run on multi-core/CPU systems allowing the decay chain calculations to be paralleled at each time-step. This will allow much finer time resolution and generation of a larger statistical base from which to draw conclusions.
Original languageEnglish
Publication statusPublished - 7 Apr 2011
EventEuropean Geosciences Union (EGU) General Assembly - Vienna, Austria
Duration: 4 Apr 2011 → …


ConferenceEuropean Geosciences Union (EGU) General Assembly
Period4/04/11 → …
Internet address


Dive into the research topics of 'Monte-Carlo simulation of radon equilibrium under varying conditions'. Together they form a unique fingerprint.

Cite this