Search NASA⌕ Search

SEARCH · Search NASA

Results for “log normal distribution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Fragmentation under the Scaling Symmetry and Turbulent Cascade with Intermittency

Fragmentation plays an important role in a variety of physical, chemical, and geological processes. Examples include atomization in sprays, crushing of rocks, explosion and impact of solids, polymer degradation, etc. Although each individual action of fragmentation is a complex process, the number of these elementary actions is large. It is natural to abstract a simple 'effective' scenario of fragmentation and to represent its essential features. One of the models is the fragmentation under the scaling symmetry: each breakup action reduces the typical length of fragments, r (right arrow) alpha r, by an independent random multiplier alpha (0 < alpha < 1), which is governed by the fragmentation intensity spectrum q(alpha), integral(sup 1)(sub 0) q(alpha)d alpha = 1. This scenario has been proposed by Kolmogorov (1941), when he considered the breakup of solid carbon particle. Describing the breakup as a random discrete process, Kolmogorov stated that at latest times, such a process leads to the log-normal distribution. In Gorokhovski & Saveliev, the fragmentation under the scaling symmetry has been reviewed as a continuous evolution process with new features established. The objective of this paper is twofold. First, the paper synthesizes and completes theoretical part of Gorokhovski & Saveliev. Second, the paper shows a new application of the fragmentation theory under the scale invariance. This application concerns the turbulent cascade with intermittency. We formulate here a model describing the evolution of the velocity increment distribution along the progressively decreasing length scale. The model shows that when the turbulent length scale gets smaller, the velocity increment distribution has central growing peak and develops stretched tails. The intermittency in turbulence is manifested in the same way: large fluctuations of velocity provoke highest strain in narrow (dissipative) regions of flow.

Gorokhovski, M.↗
NASA NTRS2000

Analysis of cell mechanics in single vinculin-deficient cells using a magnetic tweezer

A magnetic tweezer was constructed to apply controlled tensional forces (10 pN to greater than 1 nN) to transmembrane receptors via bound ligand-coated microbeadswhile optically measuring lateral bead displacements within individual cells. Use of this system with wild-type F9 embryonic carcinoma cells and cells from a vinculin knockout mouse F9 Vin (-/-) revealed much larger differences in the stiffness of the transmembrane integrin linkages to the cytoskeleton than previously reported using related techniques that measured average mechanical properties of large cell populations. The mechanical properties measured varied widely among cells, exhibiting an approximately log-normal distribution. The median lateral bead displacement was 2-fold larger in F9 Vin (-/-) cells compared to wild-type cells whereas the arithmetic mean displacement only increased by 37%. We conclude that vinculin serves a greater mechanical role in cells than previously reported and that this magnetic tweezer device may be useful for probing the molecular basis of cell mechanics within single cells. Copyright 2000 Academic Press.

Non-NASA Center↗
NASA NTRS2005

Analysis of Multi-Scale Radiometric Data Collected during the Cold Land Processes Experiment-1 (CLPX-1)

Brightness temperature maps at 18.7 and 37 GHz collected at the Fraser and North Park Meso-Scale Areas during the Cold Land Processes Experiment by the NOAA Polarimetric Scanning Radiometer (PSWA) airborne sensor are analyzed. The Fraser site is mostly covered by forest with a typical snowpack depth of 1 m while North Park has no forest cover and is characterized by patches of shallow snow. We examine histograms of the brightness temperatures at 500 m resolution for both the Fraser and North Park areas. The histograms can be modelled by a log-normal distribution in the case of the Fraser MSA and by a bi-modal distribution in the case of the North Park MSA. The histograms of the brightness temperatures at coarser resolutions are also plotted to study the effects of sensor resolution on the shape of the distribution, on the values of the average brightness temperatures and standard deviations. Finally, the values of brightness temperatures obtained by re-sampling (aggregating) the data at 25 km resolution are compared with the values of the brightness temperatures collected by the Advanced Microwave Scanning Radiometer (AMSR-E) and Special Sensor Microwave/Imager (SSMII) satellite radiometers. The results show that in both areas for sensor footprint larger than 5000 m, the brightness temperatures show a flat distribution and the memory of the initial distribution is lost. The values of the brightness temperatures measured by the satellite radiometers are in good agreement with the values obtained averaging the airborne data, even if some discrepancies occur.

Tedesco, M.↗
NASA NTRS2010

Influence of Averaging Method on the Evaluation of a Coastal Ocean Color Event on the U.S. Northeast Coast

Application of appropriate spatial averaging techniques is crucial to correct evaluation of ocean color radiometric data, due to the common log-normal or mixed log-normal distribution of these data. Averaging method is particularly crucial for data acquired in coastal regions. The effect of averaging method was markedly demonstrated for a precipitation-driven event on the U.S. Northeast coast in October-November 2005, which resulted in export of high concentrations of riverine colored dissolved organic matter (CDOM) to New York and New Jersey coastal waters over a period of several days. Use of the arithmetic mean averaging method created an inaccurate representation of the magnitude of this event in SeaWiFS global mapped chl a data, causing it to be visualized as a very large chl a anomaly. The apparent chl a anomaly was enhanced by the known incomplete discrimination of CDOM and phytoplankton chlorophyll in SeaWiFS data; other data sources enable an improved characterization. Analysis using the geometric mean averaging method did not indicate this event to be statistically anomalous. Our results predicate the necessity of providing the geometric mean averaging method for ocean color radiometric data in the Goddard Earth Sciences DISC Interactive Online Visualization ANd aNalysis Infrastructure (Giovanni).

Acker, James G.↗
NASA NTRS2012

Variability of Black Carbon Deposition to the East Antarctic Plateau, 1800-2000 AD

Refractory black carbon aerosols (rBC) from biomass burning and fossil fuel combustion are deposited to the Antarctic ice sheet and preserve a history of emissions and long-range transport from low- and mid-latitudes. Antarctic ice core rBC records may thus provide information with respect to past combustion aerosol emissions and atmospheric circulation. Here, we present six East Antarctic ice core records of rBC concentrations and fluxes covering the last two centuries with approximately annual resolution (cal. yr. 1800 to 2000). The ice cores were drilled in disparate regions of the high East Antarctic ice sheet, at different elevations and net snow accumulation rates. Annual rBC concentrations were log-normally distributed and geometric means of annual concentrations ranged from 0.10 to 0.18 m cro-g/kg. Average rBC fluxes were determined over the time periods 1800 to 2000 and 1963 to 2000 and ranged from 3.4 to 15.5 m /a and 3.6 to 21.8 micro-g/sq m/a, respectively. Geometric mean concentrations spanning 1800 to 2000 increased linearly with elevation at a rate of 0.025 micro-g/kg/500 m. Spectral analysis of the records revealed significant decadal-scale variability, which at several sites was comparable to decadal ENSO variability.

Bisiaux, M. M.↗
NASA NTRS2012

Accuracy Assessment of Aqua-MODIS Aerosol Optical Depth Over Coastal Regions: Importance of Quality Flag and Sea Surface Wind Speed

Coastal regions around the globe are a major source for anthropogenic aerosols in the atmosphere, but the underlying surface characteristics are not favorable for the Moderate Resolution Imaging Spectroradiometer (MODIS) algorithms designed for retrieval of aerosols over dark land or open-ocean surfaces. Using data collected from 62 coastal stations worldwide from the Aerosol Robotic Network (AERONET) from approximately 2002-2010, accuracy assessments are made for coastal aerosol optical depth (AOD) retrieved from MODIS aboard Aqua satellite. It is found that coastal AODs (at 550 nm) characterized respectively by the MODIS Dark Land (hereafter Land) surface algorithm, the Open-Ocean (hereafter Ocean) algorithm, and AERONET all exhibit a log-normal distribution. After filtering by quality flags, the MODIS AODs respectively retrieved from the Land and Ocean algorithms are highly correlated with AERONET (with R(sup 2) is approximately equal to 0.8), but only the Land algorithm AODs fall within the expected error envelope greater than 66% of the time. Furthermore, the MODIS AODs from the Land algorithm, Ocean algorithm, and combined Land and Ocean product show statistically significant discrepancies from their respective counterparts from AERONET in terms of mean, probability density function, and cumulative density function, which suggest a need for future improvement in retrieval algorithms. Without filtering with quality flag, the MODIS Land and Ocean AOD dataset can be degraded by 30-50% in terms of mean bias. Overall, the MODIS Ocean algorithm overestimates the AERONET coastal AOD by 0.021 for AOD less than 0.25 and underestimates it by 0.029 for AOD greater than 0.25. This dichotomy is shown to be related to the ocean surface wind speed and cloud contamination effects on the satellite aerosol retrieval. The Modern Era Retrospective-Analysis for Research and Applications (MERRA) reveals that wind speeds over the global coastal region 25 (with a mean and median value of 2.94 meters per second and 2.66 meters per second, respectively) are often slower than 6 meters per second assumed in the MODIS Ocean algorithm. As a result of high correlation (R(sup 2) greater than 0.98) between the bias in binned MODIS AOD and the corresponding binned wind speed over the coastal sea surface, an empirical scheme for correcting the bias of AOD retrieved from the MODIS Ocean algorithm is formulated and is shown to be effective over the majority of the coastal AERONET stations, and hence can be used in future analysis of AOD trend and MODIS AOD data assimilation.

Anderson, J. C.↗
NASA NTRS2018

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Sagittarius A* (Sgr A* ) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability data set of Sgr A* to date: eight 24 hr epochs of continuous monitoring of Sgr A* at 4.5 μm with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 μm data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 μm data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current fast Fourier transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is b 243 57 82 t = - + minutes (90% credible interval). The PSD has no detectable features on timescales down to 8.5 minutes (95% credible level), which is the ISCO orbital frequency for a dimensionless spin parameter a = 0.92. One light curve measured simultaneously at 2.12 and 4.5 μm during a low flux-density phase gave a spectral index αs = 1.6 ± 0.1 (F µ n s n -a ). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined data sets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a nonthermal electron population from below 20 GHz through the NIR.

accretion↗
NASA NTRS2003

On the Use of the Log-Normal Particle Size Distribution to Characterize Global Rain

Although most parameterizations of the drop size distributions (DSD) use the gamma function, there are several advantages to the log-normal form, particularly if we want to characterize the large scale space-time variability of the DSD and rain rate. The advantages of the distribution are twofold: the logarithm of any moment can be expressed as a linear combination of the individual parameters of the distribution; the parameters of the distribution are approximately normally distributed. Since all radar and rainfall-related parameters can be written approximately as a moment of the DSD, the first property allows us to express the logarithm of any radar/rainfall variable as a linear combination of the individual DSD parameters. Another consequence is that any power law relationship between rain rate, reflectivity factor, specific attenuation or water content can be expressed in terms of the covariance matrix of the DSD parameters. The joint-normal property of the DSD parameters has applications to the description of the space-time variation of rainfall in the sense that any radar-rainfall quantity can be specified by the covariance matrix associated with the DSD parameters at two arbitrary space-time points. As such, the parameterization provides a means by which we can use the spaceborne radar-derived DSD parameters to specify in part the covariance matrices globally. However, since satellite observations have coarse temporal sampling, the specification of the temporal covariance must be derived from ancillary measurements and models. Work is presently underway to determine whether the use of instantaneous rain rate data from the TRMM Precipitation Radar can provide good estimates of the spatial correlation in rain rate from data collected in 5(sup 0)x 5(sup 0) x 1 month space-time boxes. To characterize the temporal characteristics of the DSD parameters, disdrometer data are being used from the Wallops Flight Facility site where as many as 4 disdrometers have been used to acquire data over a 2 km path. These data should help quantify the temporal form of the covariance matrix at this site.

Meneghini, Robert↗
NASA NTRS2008

A New Bond Albedo for Performing Orbital Debris Brightness to Size Transformations

We have developed a technique for estimating the intrinsic size distribution of orbital debris objects via optical measurements alone. The process is predicated on the empirically observed power-law size distribution of debris (as indicated by radar RCS measurements) and the log-normal probability distribution of optical albedos as ascertained from phase (Lambertian) and range-corrected telescopic brightness measurements. Since the observed distribution of optical brightness is the product integral of the size distribution of the parent [debris] population with the albedo probability distribution, it is a straightforward matter to transform a given distribution of optical brightness back to a size distribution by the appropriate choice of a single albedo value. This is true because the integration of a powerlaw with a log-normal distribution (Fredholm Integral of the First Kind) yields a Gaussian-blurred power-law distribution with identical power-law exponent. Application of a single albedo to this distribution recovers a simple power-law [in size] which is linearly offset from the original distribution by a constant whose value depends on the choice of the albedo. Significantly, there exists a unique Bond albedo which, when applied to an observed brightness distribution, yields zero offset and therefore recovers the original size distribution. For physically realistic powerlaws of negative slope, the proper choice of albedo recovers the parent size distribution by compensating for the observational bias caused by the large number of small objects that appear anomalously large (bright) - and thereby skew the small population upward by rising above the detection threshold - and the lower number of large objects that appear anomalously small (dim). Based on this comprehensive analysis, a value of 0.13 should be applied to all orbital debris albedo-based brightness-to-size transformations regardless of data source. Its prima fascia genesis, derived and constructed from the current RCS to size conversion methodology (SiBAM Size-Based Estimation Model) and optical data reduction standards, assures consistency in application with the prior canonical value of 0.1. Herein we present the empirical and mathematical arguments for this approach and by example apply it to a comprehensive set of photometric data acquired via NASA's Liquid Mirror Telescopes during the 2000-2001 observing season.

Mulrooney, Mark K.↗
NASA NTRS1988

Rainfall index over oceans derived from SSM/I data

Ths Special Sensor Microwave/Imager radiometer on board the DMSP satellite measured microwave radiation at 19.35, 22.235, 37.0, and 85.5 GHz with a swath width of 1400 km, providing an opportunity to study global precipitation distributions. A monthly averaged rainfall index was derived using only the 19.35 GHz data. It covers the + or - 50 deg portion of the world with a 5 deg x 5 deg grid. The brightness temperature histogram is analyzed to derive the rainfall distribution. Estimates of the freezing level heights and the water vapor contents are utilized to understand the rainfall distribution better. Log normal and gamma distribution functions are fitted to the data and the derived monthly rainfall index compares favorably with historical precipitation patterns.

Chang, A. T. C.↗
DOE OSTI2016

WE-H-207A-03: The Universality of the Lognormal Behavior of [F-18]FLT PET SUV Measurements

Purpose: Log transforming [F-18]FDG PET standardized uptake values (SUVs) has been shown to lead to normal SUV distributions, which allows utilization of powerful parametric statistical models. This study identified the optimal transformation leading to normally distributed [F-18]FLT PET SUVs from solid tumors and offers an example of how normal distributions permits analysis of non-independent/correlated measurements. Methods: Forty patients with various metastatic diseases underwent up to six FLT PET/CT scans during treatment. Tumors were identified by nuclear medicine physician and manually segmented. Average uptake was extracted for each patient giving a global SUVmean (gSUVmean) for each scan. The Shapiro-Wilk test was used to test distribution normality. One parameter Box-Cox transformations were applied to each of the six gSUVmean distributions and the optimal transformation was found by selecting the parameter that maximized the Shapiro-Wilk test statistic. The relationship between gSUVmean and a serum biomarker (VEGF) collected at imaging timepoints was determined using a linear mixed effects model (LMEM), which accounted for correlated/non-independent measurements from the same individual. Results: Untransformed gSUVmean distributions were found to be significantly non-normal (p<0.05). The optimal transformation parameter had a value of 0.3 (95%CI: −0.4 to 1.6). Given the optimal parameter was close to zero (which corresponds to log transformation), the data were subsequently log transformed. All log transformed gSUVmean distributions were normally distributed (p>0.10 for all timepoints). Log transformed data were incorporated into the LMEM. VEGF serum levels significantly correlated with gSUVmean (p<0.001), revealing log-linear relationship between SUVs and underlying biology. Conclusion: Failure to account for correlated/non-independent measurements can lead to invalid conclusions and motivated transformation to normally distributed SUVs. The log transformation was found to be close to optimal and sufficient for obtaining normally distributed FLT PET SUVs. These transformations allow utilization of powerful LMEMs when analyzing quantitative imaging metrics.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1988

Effects of El Chichon volcanic effluents on stratospheric aerosol dynamics

The effects of El Chichon's April 1982 eruption on stratospheric aerosol dynamics are presently discussed in terms of log-normal size distributions over 15-20 km sample altitudes between 30 and 45 deg N over the contiguous U.S. After collection, samples were studied by SEM, and log-normal size distributions were fitted to the data-points obtained. It is found that stratospheric aerosol behavior is explainable by the laws of aerosol mechanics more easily than has been the case for tropospheric aerosol, for which the source-sink relationship is much more complex.

Pueschel, R. F.↗
DOE OSTI2018

Bidisperse and polydisperse suspension rheology at large solid fraction

At the same solid volume fraction, bidisperse and polydisperse suspensions display lower viscosities, and weaker normal stress response, compared to monodisperse suspensions. The reduction of viscosity associated with size distribution can be explained by an increase of the maximum flowable, or jamming, solid fraction. In this work, concentrated or "dense" suspensions are simulated under strong shearing, where thermal motion and repulsive forces are negligible, but we allow for particle contact with a mild frictional interaction with interparticle friction coefficient of 0.2. Aspects of bidisperse suspension rheology are first revisited to establish that the approach reproduces established trends; the study of bidisperse suspensions at size ratios of large to small particle radii (2 to 4) shows that a minimum in the viscosity occurs for zeta slightly above 0.5, where zeta=phi_{large}/phi is the fraction of the total solid volume occupied by the large particles. The simple shear flows of polydisperse suspensions with truncated normal and log normal size distributions, and bidisperse suspensions which are statistically equivalent with these polydisperse cases up to third moment of the size distribution, are simulated and the rheologies are extracted. Prior work shows that such distributions with equivalent low-order moments have similar phi_{m}, and the rheological behaviors of normal, log normal and bidisperse cases are shown to be in close agreement for a wide range of standard deviation in particle size, with standard correlations which are functionally dependent on phi/phi_{m} providing excellent agreement with the rheology found in simulation. The close agreement of both viscosity and normal stress response between bi- and polydisperse suspensions demonstrates the controlling in influence of the maximum packing fraction in noncolloidal suspensions. Microstructural investigations and the stress distribution according to particle size are also presented.

Pednekar, Sidhant↗
NASA NTRS2017

Advanced Characterisation of Aerosol Size Properties from Measurements of Spectral Optical Depth Using the GRASP Algorithm

This study evaluates the potential of using aerosol optical depth (τ (sub a)) measurements to characterise the microphysical and optical properties of atmospheric aerosols. With this aim, we used the recently developed GRASP (Generalized Retrieval of Aerosol and Surface Properties) code for numerical testing of six different aerosol models with different aerosol loads. The direct numerical simulations (self-consistency tests) indicate that the GRASP-AOD retrieval provides modal aerosol optical depths (fine and coarse) to within 0.01 of the input values. The retrieval of the fine-mode radius, width and volume concentration are stable and precise if the real part of the refractive index is known. The coarse-mode properties are less accurate, but they are significantly improved when additional a priori information is available. The tests with random simulated errors show that the uncertainty in the bimodal log-normal size distribution parameters increases as the aerosol load decreases. Similarly, the reduction in the spectral range diminishes the stability of the retrieved parameters. In addition to these numerical studies, we used optical depth observations at eight AERONET locations to validate our results with the standard AERONET inversion products. We found that bimodal log-normal size distributions serve as useful input assumptions, especially when the measurements have inadequate spectral coverage and/or limited accuracy, such as moon photometry. Comparisons of the mode median radii between GRASP-AOD and AERONET indicate average differences of 0.013µm for the fine mode and typical values of 0.2–0.3µm for the coarse mode. The dominant mode (i.e. fine or coarse) indicates a 10 % difference in mode radii between the GRASP-AOD and AERONET inversions, and the average of the difference in volume concentration is around 17 % for both modes. The retrieved values of the fine-mode τ(sub a)(500) using GRASP-AOD are generally between those values obtained by the standard AERONET inversion and the values obtained by the AERONET spectral deconvolution algorithm (SDA), with differences typically lower than 0.02 between GRASP-AOD and both algorithms. Finally, we present some examples of application of GRASP-AOD inversion using moon photometry and the airborne PLASMA sun photometer during the ChArMEx summer 2013 campaign in the western Mediterranean.

Torres, Benjamin↗
NASA NTRS1986

The frequency spectrum of C sub n sup 2 from MST radar data

In a recent study (Nastrom et al., 1986), the variability of the refractivity turbulence structure constant, C sub n(2), was examined using observations from the stratosphere troposphere/mesosphere stratosphere troposphere (ST/MST) radar at Poker Flat, Alaska, and Platteville, Colorado. Variations of C2 with height, season, and weather conditions were examined. Also, the autocorrelation function and the frequency distribution of C sub n(2) were studied, and it was shown that C sub n(2), follows a log-normal frequency distribution. One of the more tentative results given in that paper is a first look at the spectrum of log C sub n(2), as a function of frequency at Poker Flat. This spectrum appears to obey a power law relation with frequency, P(F) approx. F(k), with k near -5/3 at periods between about 4 hours and 6 days, and with k near -1 at shorter periods. Power law behavior of a spectrum often helps us to infer the underlying dynamics which give rise to this spectrum, and it is thus of some concern to establish further confidence in the spectral shape. The purpose here is to address these questions.

Nastrom, G. D.↗
NASA NTRS1978

Aerosol size distributions obtained by inversion of spectral optical depth measurements

Columnar aerosol size distributions have been inferred by numerically inverting particulate optical depth measurements as a function of wavelength. An inversion formula which explicitly includes the magnitude of the measurement variances is derived and applied to optical depth measurements obtained in Tucson with a solar radiometer. It is found that the individual size distributions of the aerosol particles (assumed spherical), at least for radii greater than or approximately equal to 0.1 micron, fall into one of three distinctly different categories. Approximately 50% of all distributions examined thus far can best be represented as a composite of a Junge distribution plus a distribution of relatively monodispersed larger particles centered at a radius of about 0.5 micron. Scarcely 20% of the distributions yielded Junge size distributions, while 30% yielded relatively monodispersed distributions of the log-normal or gamma distribution types. A representative selection of each of these types will be presented and discussed. The sensitivity of spectral attenuation measurements to the radii limits and refractive index assumed in the numerical inversion will also be addressed.

King, M. D.↗
NASA NTRS2000

The BATSE Gamma-Ray Burst E-Peak Distribution

Gamma-ray burst observed by BATSE are found to have approximately the same characteristic energy, denoted as E_p. We examine whether instrumental effects can give rise to this observation. We simulate the derivation of E-p and determine that the values in the BATSE sample are accurate and complete above a minimum fluence. We simulate the triggering of BATSE on gamma-ray bursts, deriving the efficiency of detecting bursts as a function of characteristic energy. From this simulation, we model the observed E_p distribution function expected when the intrinsic distribution function is a power-law. We find that this distribution produces poor fits to the observations. We find that a log-normal intrinsic distribution with a power-law tail gives a good fit to the data. From these fits, we conclude that instrumental effects cannot produce the observed E_p distribution, and that the observed distribution is a consequence of a narrow intrinsic distribution of E_p in gamma-ray bursts.

Brainerd, Jerome J.↗
DOE OSTI2014

Oxygen distribution in tumors: A qualitative analysis and modeling study providing a novel Monte Carlo approach

Purpose: To construct a Monte Carlo (MC)-based simulation model for analyzing the dependence of tumor oxygen distribution on different variables related to tumor vasculature [blood velocity, vessel-to-vessel proximity (vessel proximity), and inflowing oxygen partial pressure (pO{sub 2})]. Methods: A voxel-based tissue model containing parallel capillaries with square cross-sections (sides of 10 μm) was constructed. Green's function was used for diffusion calculations and Michaelis-Menten's kinetics to manage oxygen consumption. The model was tuned to approximately reproduce the oxygenational status of a renal carcinoma; the depth oxygenation curves (DOC) were fitted with an analytical expression to facilitate rapid MC simulations of tumor oxygen distribution. DOCs were simulated with three variables at three settings each (blood velocity, vessel proximity, and inflowing pO{sub 2}), which resulted in 27 combinations of conditions. To create a model that simulated variable oxygen distributions, the oxygen tension at a specific point was randomly sampled with trilinear interpolation in the dataset from the first simulation. Six correlations between blood velocity, vessel proximity, and inflowing pO{sub 2} were hypothesized. Variable models with correlated parameters were compared to each other and to a nonvariable, DOC-based model to evaluate the differences in simulated oxygen distributions and tumor radiosensitivities for different tumor sizes. Results: For tumors with radii ranging from 5 to 30 mm, the nonvariable DOC model tended to generate normal or log-normal oxygen distributions, with a cut-off at zero. The pO{sub 2} distributions simulated with the six-variable DOC models were quite different from the distributions generated with the nonvariable DOC model; in the former case the variable models simulated oxygen distributions that were more similar to in vivo results found in the literature. For larger tumors, the oxygen distributions became truncated in the lower end, due to anoxia, but smaller tumors showed undisturbed oxygen distributions. The six different models with correlated parameters generated three classes of oxygen distributions. The first was a hypothetical, negative covariance between vessel proximity and pO{sub 2} (VPO-C scenario); the second was a hypothetical positive covariance between vessel proximity and pO{sub 2} (VPO+C scenario); and the third was the hypothesis of no correlation between vessel proximity and pO{sub 2} (UP scenario). The VPO-C scenario produced a distinctly different oxygen distribution than the two other scenarios. The shape of the VPO-C scenario was similar to that of the nonvariable DOC model, and the larger the tumor, the greater the similarity between the two models. For all simulations, the mean oxygen tension decreased and the hypoxic fraction increased with tumor size. The absorbed dose required for definitive tumor control was highest for the VPO+C scenario, followed by the UP and VPO-C scenarios. Conclusions: A novel MC algorithm was presented which simulated oxygen distributions and radiation response for various biological parameter values. The analysis showed that the VPO-C scenario generated a clearly different oxygen distribution from the VPO+C scenario; the former exhibited a lower hypoxic fraction and higher radiosensitivity. In future studies, this modeling approach might be valuable for qualitative analyses of factors that affect oxygen distribution as well as analyses of specific experimental and clinical situations.

60 APPLIED LIFE SCIENCES↗