prior to the invention , existing codes that implement the icrp - 66 respiratory tract model provide deterministic , point estimates of organ or effective dose for a given set of inhalation exposure parameters . the invention provides a computer code , termed luduc , which allows a user to solve the deposition , clearance , and dosimetry components of the icrp - 66 respiratory tract model using stochastic , as opposed to deterministic , sampling of input parameter values for example by using latin hypercube techniques . for each of the 69 input parameters of the model , probability density functions can be assigned , rather than single default values . the code provides not only mean and median estimates of doses under any selected set of exposure conditions , but also information on the total uncertainty ( for example , at the 95 % confidence level ) for those doses . such information on dose uncertainty is extremely useful , both for demonstrating compliance with a regulatory dose limit , and for dose reconstruction analyses , which involve correlating past worker doses with observed biological effects , such as diseases . both organ equivalent doses and whole - body effective doses can be determined for at least 233 potentially inhaled radionuclides , including those emitting alpha - particles , beta - particles , and x - or gamma - ray photons . the luduc code permits the user to select various aspects of the exposure scenario , such as radionuclide , particle size distribution and solubility of the particles in the lungs . the code further provides for predictions that take into account biological variability among individuals of an exposed population , including an individual &# 39 ; s sex , age , exertion level , body height , and body mass index . as shown in the examples described herein , variability in one or more of these parameters can significantly affect the confidence levels of the predicted level of radiation exposure . referring now to fig1 , luduc can be structured into nine modules as indicated in the flowchart . the scenario specification / input module is written in visual basic v6 . 0 and can run in the microsoft windows 2000 / xp environment . this programming language makes full use of the graphical user interface and multitasking capabilities offered by the ms windows environment . the assessment problem requires an exposure scenario to be defined by specifying : ( 1 ) an age and gender of an exposed population group , ( 2 ) a physical exertion level of the group during the exposure ( assumed to be acute for computational purposes ), ( 3 ) an activity - particle size distribution , a particle shape factor , and a particle density , and ( 4 ) an ambient temperature and pressure . the user can optionally specify an ambient activity concentration level , and an exposure duration . otherwise , quantities are assessed per unit exposure . the code can allow for monodisperse particles , a lognormal particle size distribution , a uniform distribution , or a user - supplied histogram ( normalized ) of particle sizes . the source code of modules 1 and 9 is written in visual basic . the source code of modules 2 through 8 is written in fortran and thus is portable to other operating systems . modules 2 through 8 have been compiled / linked using lahey fortran 90 v4 . 50 . the second module embodies latin hypercube sampling ( lhs ) algorithms developed by iman and shortencarier ( 1984 ) at sandia national laboratories , with minor modifications to the source code . the lhs module reads an input file created by the first module , creating a matrix of n input parameter arrays to be utilized by luduc in the n trials to be run in the uncertainty analysis . the third module solves the particle deposition component of the icrp 66 respiratory tract model n times , to generate n predictions of particle deposition fractions in the various regions of the lung . the fourth module solves the clearance component of the respiratory tract model and implements an algorithm described by birchall ( 1986 ) to solve the resulting system of differential equations . the purpose of this module is to solve the clearance model n times to generate n predictions of either the number of nuclear transformations or the transformation rate in source components after a specified time . the fifth module can be a stand - alone program that computes a dose matrix for alpha particles used as input to the dose computation module . the sixth and seventh modules are a set of programs in visual basic and fortran used to run monte carlo n - particle ( mcnp ) quickly and efficiently for monoenergetic beta particles and photons , respectively . this set of programs was also utilized to create absorbed fraction ( beta particles ) and specific absorbed fraction ( photons ) data for 233 radionuclides . these data tables are used by luduc when it runs . the eighth module computes equivalent doses in target tissues and the combined lung dose ( for example , the weighted sum of the regional doses ). this module couples results from clearance model computations with target and source geometries . using data generated by the lhs module for parameters such as target and source dimensions and masses ( based on assigned input distributions ), this module solves the dose model n times , for n values of the various dose quantities . module 8 also estimates the effective dose and considers alpha - particle , beta - particle , and photon emitters . the ninth module gives the numerical results from the deposition , clearance , and dosimetry modules in terms of histograms with their corresponding statistical parameters . these parameters can include the minimum , maximum , median , mean , standard deviation , coefficient of variance , geometric mean , geometric standard deviation , and several percentiles for predicted quantities . these quantities can include ( 1 ) deposition fractions in the et 1 , et 2 , bb , bb , ai , total thoracic , and total respiratory tract , ( 2 ) total radionuclide transformations at time t since exposure for all source components in the respiratory tract , and ( 3 ) equivalent doses to these target tissues . the results of the dose computation module are presented by the ninth module for individual radiation types and their overall contribution to equivalent dose . for short - range particles such as alpha and beta particles , module 9 shows results for the extrathoracic and thoracic target tissues . for photons , this module displays results for 34 target organs and tissues . the results for the lung include the contributions of the alpha particles , beta particles , and photons emitted by the radionuclide selected by the luduc user . this module is written in visual basic and runs in the ms windows environment . finally , this module presents the effective dose in sieverts . if needed , the user can easily obtain the data generated by luduc for storage or further processing in spreadsheet programs . modules 1 through 9 can all be run from a single computational platform in ms windows . modules 2 , 3 , 4 , 5 , 6 , 7 and 8 are run by a shelling command , which temporarily transfers control from windows to ms - dos . for a sample size of n = 1 , this platform was shown to run in less than a second on a 2 . 0 - ghz personal computer system . to start luduc , the user clicks on a luduc icon ( not shown ). when luduc is executed , it displays the version number and authors , and a luduc banner . the user can click on the banner to start luduc , or the program will start automatically in several seconds after the luduc banner is displayed . the main menu appears once the banner disappears ( fig2 ). the main menu shows the main commands in luduc such as open an existing project file , save project file , print input / output files , edit / view input data for uncertainty analysis , run main program , view results of parameter uncertainty analysis and exit luduc . the main menu also displays the following filenames : project file , exposure / scenario data , model input parameters , deposition calculations , clearance calculations , dose calculations ( for alpha - and beta - particles , and gamma - rays ) and total dose calculations . these files store the necessary information to run luduc . the user can modify these filenames by clicking on the specific filename and editing it . the main menu also shows the time and date . selection of the main menu command open an existing project file ( shown in fig2 ) allows the user select an initial project file to run luduc ( fig3 ). the project file ( for example , luduc . prj ) is the filename of the file containing the necessary data for the input parameters and output ( i . e ., results ). to modify and / or verify the input parameters , the user clicks on the main menu command edit / view input data for uncertainty analysis , causing display of a window , edit / view input parameter distributions , shown in fig4 . this window is subdivided into several parts , including an exposure scenario and computational method setup window , a respiratory tract model parameter setup window , and an output quantity selection window . the exposure scenario and computational method setup window includes on the right a summary of the exposure scenario and computational setup . five commands , displayed on the left , allow the user to modify the corresponding data . referring to the upper left of fig4 , by selecting the command select population class , the user is able to view or edit the gender , age , exertion level , and fraction of nose breathers for a specific population group , as shown in fig5 . by selecting the specify aerosol characteristics command ( fig4 ), the user is able to view or edit the aerosol activity size distribution , in a window illustrated in fig6 . distribution choices include lognormal , uniform , and user - supplied histogram . aerosol characteristics that the user can view or edit include an activity median aerodynamic diameter ( in microns ), a geometric standard deviation for diameter , an aerosol shape factor , and an aerosol density ( in g / cm 3 ). referring again to the upper left of fig4 , by activating the command specify exposure conditions , the user can view or edit the atmospheric activity concentration , duration of exposure , ambient atmospheric temperature and pressure , and the radionuclide ( s ) for the specific scenario , as shown in fig7 . to select the radionuclide ( s ), the user clicks on periodic table , to display a window illustrating the periodic table of the elements ( fig8 ). this window allows the user to select the radionuclide ( s ) in two different ways : ( 1 ) by clicking on a desired element in the display of elements or ( 2 ) by selecting a radionuclide from the dropdown radionuclide list . below the periodic table , important information about the selected radionuclide is shown , such as half - life , and number of alpha and beta particles emitted . as examples of the type of information displayed , fig8 shows the selection of the main components of weapons grade plutonium and fig9 shows the selection of iodine - 131 . having selected the radionuclide , the user can input the percent of total activity in the radionuclide checklist ( fig8 and 9 ). by selecting the command specify run parameters ( fig4 ), the user can view or edit the parameter setup , which includes the sampling method ( i . e ., simple random sampling or latin hypercube sampling ), the number of trials or realizations , the distribution truncation level , and the random number seed value . a run setup parameter window is shown in fig1 . the random number seed can be selected by a direct input from the user , by a random number generated by the visual basic , or by computer time . referring now to the middle of the window shown in fig4 , the respiratory tract model parameter setup window allows the user to modify and / or verify the parameters for the deposition , clearance , and dose models . by clicking on the command deposition model parameters , the user can view or edit the quantity , distribution type , and distribution parameters of the deposition model , as shown in fig1 a . by double - clicking on a quantity , such as body height , the user can modify the distribution characteristics , for example as shown in fig1 b . four different distribution types are available , including normal ( gaussian ), lognormal , uniform , and triangular . the user selects the distribution and necessary parameters ( mean and standard deviation for a normal distribution ), then accepts the choice by clicking on the accept command . selecting the command clearance model parameters in the respiratory tract model parameter setup ( middle of fig4 ), the user can view or edit the quantity , distribution type , and distribution parameters of the clearance model . ( see fig1 .) the procedure to modify the distribution parameters of the clearance model is similar to that for modifying the parameters of the deposition model . by selecting the dose model parameters command in the respiratory tract model parameter setup ( fig4 ), the user can view or edit the quantity , distribution type , and distribution parameters of the dose model , as shown in fig1 . the procedure to modify the distribution parameters of the dose model is similar to that for the deposition and clearance models . referring to the lower portion of fig4 , the output quantity selection window allows the user to compute integrated doses or dose rates based on elapsed time ( for example , days ) after exposure to radiation . referring now to fig1 , to execute luduc , the user selects run main program . luduc indicates to the user when it is finished running . by selecting the command view results of parameter uncertainty analysis , the user can see the results of the run under the specified conditions . a typical window showing such results is illustrated in fig1 . a summary of input parameters ( population class , particle size distribution , exposure conditions , and distribution sampling techniques ) appears in the upper portion of the window . to observe the results for the deposition , clearance , and dose models , the user clicks on the corresponding commands , shown in the lower left of fig1 . to return to the previous window , the user can click on a return command . when selecting the command deposition calculations ( fig1 ), the user is able to observe the deposition model results ( deposition fractions in various regions of the respiratory tract ) as shown in fig1 . the user is able to modify the horizontal scale of the plots by clicking on the command abscissa scale . this modification of the horizontal scale applies to all the plots in luduc . to view probabilistic results , the user can click on the individual plots , to expand the graph and to see statistical data , such as that shown in fig1 . the feature of clicking on a plot to view the statistics applies to all plots displayed in the luduc program . displayed statistics include the number of observations ( trials ), minimum , maximum , maximum - to - minimum ratio , median , mean , standard deviation , coefficient of variability , geometric mean , and geometric standard deviation . to return to the previous window ( fig1 ), the user clicks on the plot itself . referring again to fig1 , the user can select the commands clearance calculations : extrathoracic and clearance calculations : thoracic , to view the corresponding plots and statistics for the number of disintegrations or the rate of disintegrations ( fig1 and 19 ). again referring to fig1 , to observe the results from the dosimetry model , the user can select one of the several dose calculations commands . the results for the dose model are divided into three main groups : alpha particles ( sample data shown in fig2 - 22 ), beta particles ( fig2 ), and gamma - rays / x - rays ( fig2 - 26 ). the user is also able to see the median , weighted median , and percentiles by clicking on the corresponding commands in the window shown in fig2 . referring again to fig1 , the dose calculation results for gamma - rays / x - ray are further divided into three main groups ( i . e ., groups a , b , and c ). each of these groups includes a number of target organs or tissues , as can be seen in fig2 - 26 . from the menu shown in fig1 , the user is able to determine which tissues or organs are included in each group by locating the mouse pointer over the commands for dose calculations for organ groups a , b , or c . the user is also able to observe the results for equivalent doses and effective dose by clicking on the corresponding commands shown in the menu in fig1 . examples of the displayed results for tissue equivalent dose ( sv ) and effective dose ( sv ) are shown in fig2 and 28 , respectively . predicting the deposition behavior of aerosols in the respiratory tract is necessary to estimate the fractions of radioactivity that are deposited in each anatomical region of the lungs . these fractions are required in the assessment of health risks associated with the inhalation of radioactive aerosols . methods : a complete respiratory tract deposition methodology based on the icrp 66 respiratory tract model was used in an analysis of plutonium , uranium and americium oxide aerosol particles . lung deposition fractions were estimated as probability distributions to reflect the variability or spread in the deposition values . the methodology was implemented using the luduc computer code . results . the deposition fractions followed a lognormal distribution shape for all exposure scenarios examined . in general , median distribution fractions generated by luduc agreed with the reference deposition fractions of the deterministic computer code ludep . however , the results showed that the particle aerodynamic diameter and the physical exertion level ( sleeping , resting , light exertion , and high exertion ) strongly influenced the deposition uncertainty . estimating respiratory tract clearance rates of radioactive aerosols is essential in the estimation of health risks associated with the inhalation of radioactive aerosols and vapors . accurate methodology of clearance kinetics is required because respiratory tract clearance rates determine not only doses to the respiratory tract tissues , but also doses to other organs following systemic uptake . aerosols deposited in the respiratory tract are cleared to the gastrointestinal tract via the pharynx , the regional lymph nodes via the lymphatic channels , and blood via absorption . in general , the rate at which deposited aerosols are cleared depends on the time elapsed since the deposition of aerosols , the physicochemical form of the aerosols , and the location of the aerosols in the respiratory tract . methods : a detailed respiratory tract clearance methodology based on the ircp 66 respiratory tract model was used to study 241 am , 235 u , 238 u , and 239 pu oxide aerosols . the methodology utilized luduc , a computer code that permits lung clearance rates to be calculated as probability distributions to reveal the spread in clearance rate values . results . the clearance rates had a lognormal distribution shape for all examined exposure conditions . the results showed that clearance uncertainty is highly subject to the physicochemical properties of the aerosols . a complete respiratory tract model for predicting lung dosimetry of inhaled radioactive aerosols involves several component models , including models for particle deposition in the airways , biokinetic clearance , radiological decay of deposited materials , and radiological dose to critical target tissues . each component depends on several parameters , which can vary among members of a population group . methods . a methodology was developed , based on conducting parameter uncertainty analyses , to incorporate parameter uncertainties into the predictions of lung doses . results of previous studies were compiled to recommend distributions representative of parameter uncertainties , and the methodology was implemented using luduc , an interactive computer program . doses resulting from inhalation of uranium and plutonium oxide aerosols with aerodynamic diameters ranging from 0 . 1 to 50 microns were investigated . results . dose distributions followed a lognormal distribution shape for all exposure scenarios examined . median doses for uranium and plutonium oxide generally agreed with reference dose values , providing some level of confidence in the approach using reference man . differences in the predicted dose distributions were small when comparing different age and gender groups from 2 to 35 years of age . a sensitivity analysis of all model parameters within the icrp - 66 respiratory tract model is essential for the estimation of probabilistic dose distribution in lung dosimetry . methods . this analysis was performed to identify those model parameters which most influence model predictions , and to determine the contribution made by parameter variabilities to uncertainties in the model predictions . sensitivity analyses were conducted for adult males , 25 - 34 years old , exposed to 239 puo 2 aerosols at a light exertion level , assuming acute deposition using the rank - transformed dose and deposition data generated by the computer code luduc . the sensitivities of model predictions on input variables were performed by determining the standardized rank regression coefficients ( srrcs ) of selected input variables . based on absolute values of their associated ssrcs , input variables were ranked in increasing values from one , with the most important , i . e ., the most sensitive variable being assigned a rank of one . the data were generated by performing n = 1000 trials using latin hypercube sampling techniques . results . in general , calculated uncertainties generally increase as the particle diameter increases from 0 . 1 to 50 μm . however , the calculated median dose decreases with increasing particle diameter over this same size range . generally , uncertainties in lung and tissue equivalent doses can be modeled by lognormal distributions . sensitivities in dose predictions differed between target tissues and were influenced by particle size , due primarily to dependencies in the deposition model . the ssrcs technique was generally able to explain over 90 % of the variablility in the dose and deposition predictions . for the deposition component of the respiratory tract model , a larger portion of the variability in deposition and dose model predictions was attributable to only a few model parameters . the extrathoracic airways and lymph nodes have not been previously represented explicitly in mathematical models of the human body which are utilized to predict transport of photons internally between source and target organs within the body . the current icpr respiratory tract model assumes that the extrathoracic airways are reasonably approximated by using the thyroid as a surrogate source and target region . consequently , the thyroid replaces the extrathoracic airways and lymph nodes ( et 1 , et 2 , and ln et ) as the emission site or deposition site for photons released from inhaled radioactive particulates . methods . a new mathematical model was created to explicitly represent the extrathoracic airways , as well as other respiratory structures in the thorax of the adult . the model incorporated the revised dosimetric medical internal radiation dose ( mird ) model of the adult head and brain and the oak ridge national laboratories model of the adult male . several modifications were made , to include a number of organs and tissue regions absent from previous models . the resulting mathematical model included an external nose , nasal cavity , nasal sinuses ( frontal , ethmoid , sphenoid , and maxillary ), larynx , pharynx , trachea , main bronchi , and esophagus . the model was implemented into the mcnp radiation transport code to determine specific absorbed fractions . the specific absorbed fractions and the new mathematical phantom were incorporated into the luduc computer program . ( see fig1 .) results . the et 1 , et 2 , and ln et regions represented a more realistic mathematical model of the human respiratory tract tissues , enabling more accurate estimation of uncertainties in dose within the icrp - 66 respiratory tract model for photon emitters . beta - particle uncertainty within the icrp - 66 respiratory tract model impact of uncertainties in electron absorbed fractions on lung dose estimates this analysis was performed to investigate the short - range dosimetry model of the icrp - 66 respiratory tract model whereby probability density functions are assigned for target depths , thicknesses , and masses . methods . the luduc probabilistic computer code was modified to include capability to analyze beta - particle emitters . to create the data files , monte carlo transport simulations were performed for beta particles . luduc was then used to assess regional and total lung doses from inhaled aerosols of 90 sr and 90 y compounds . results . dose uncertainty was found to depend mainly on particle size . for strontium and yttrium compounds of the inhalation class y , the results showed that the spread in lung dose increased by factors of about 10 over the particle size range from 0 . 001 to 10 μm . the ratio of the 95 % to 5 % fractile was relatively constant for particle diameters of 0 . 01 to 0 . 2 μm , i . e ., 10 and 3 for 90 sr and 90 y , respectively . this difference increased to about a factor of 100 as the particle diameter approached 10 μm . this was mainly due to the fact that thoracic doses become low at larger particle sizes because most of the deposition occurs in the extrathoracic region . analysis of uncertainties in the electron absorbed fractions within the icrp - 66 respiratory tract model the uncertainties of beta - particle transport and energy deposition were analyzed . for short - ranged beta particles , critical parameters of dose assessment are based in part on estimates of target tissue depths , thickness , and masses , predominantly within thoracic regions of the respiratory tract . methods . to model uncertainty in doses for beta particles , probability density functions were assigned for target tissue depths , thickness and masses , using luduc , as described in example 6 . unlike the methodology utilized by the icrp - 66 model for alpha particles , in which range - energy relationships are used to account for alpha particle deposition in lung tissues , full monte carlo radiation transport simulations were made for beta particles due to their non - linear pathlengths within these tissues . egs4 code was used in the icrp - 66 model to simulate beta - particle energy deposition and absorbed fractions in lung airways . the complexity of the work was significantly simplified due to the fixed geometry for both target cell depths and thicknesses , and source tissue depths and thicknesses . both of these combinations of distances were varied stochastically . furthermore , in situations in which an intermediate tissue was located between the source and target tissues , the thickness of this intermediate region could be varied . as a result , a new scheme was developed and implemented into the mcnp 4c radiation transport code . results . in the new scheme , the airways of the bb ( bronchial ) and bb ( bronchiolar ) regions are subdivided into thin ( i . e ., 1 μm thick ) cylindrical shells . in general , each shell is considered as a potential source . this methodology enables assessment of regional and total lung dose from inhaled aerosols of beta - particles , such as 90 sr and 90 y compounds . while the above specification contains many specifics , these should not be construed as limitations on the scope of the invention , but rather as examples of preferred embodiments thereof . many other variations are possible . accordingly , the scope of the invention should be determined not by the embodiments illustrated , but by the appended claims and their legal equivalents .
systems , machine readable storage , and methods for solving a respiratory tract model of radiation dosage . steps include selecting input parameters associated with a respiratory tract model , computing a probability density function for the input parameters , and solving the respiratory tract model using the computed probability density functions .