Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Radiative property investigation

Measurements of the bidirectional reflections of NASA paints Z-93 and S-13G were made for 0.5, 1.78, and 2.5 microns, source incidence angles of 0, 30, 60, and 75 deg, and four detector azimuth positions. The bidirectional reflectance distribution was found to be near Lambert for azimuth positions not to close to 180 deg from the source incident azimuth. The reflectance distribution at an azimuth of 180 deg from the source is close to Lambert for near normal source incidents angles, and becomes much larger than Lambert as the source incident angle approach grazing. The reflectance deviates from Lambert as the wavelength increases. Comparison of data for the paints shows Z-93 reflectance is much closer to Lambert than S-13G. When there is a peak in the reflectance distribution it occurred at the specular angle or an angle larger than the specular angle.

Scott, R. L., Jr.↗
DOE OSTI2018

A Framework for Robust Long-Term Voltage Stability of Distribution Systems

Power injection uncertainties in distribution power grids, which are mostly induced by aggressive introduction of intermittent renewable sources, may drive the system away from normal operating regimes and potentially lead to the loss of long-term voltage stability (LTVS). Naturally, there is an ever increasing need for a tool for assessing the LTVS of a distribution system. This paper presents a fast and reliable tool for constructing inner approximations of LTVS regions in multidimensional injection space such that every point in our constructed region is guaranteed to be solvable. Numerical simulations demonstrate that our approach outperforms all existing inner approximation methods in most cases. Furthermore, the constructed regions are shown to cover substantial fractions of the true voltage stability region. In conclusion, the paper will later discuss a number of important applications of the proposed technique, including fast screening for viable injection changes, constructing an effective solvability index and rigorously certified loadability limits.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Rotation curves and dynamical masses of MaNGA barred galaxies

In this paper we analyse a sample of 46 barred galaxies of MaNGA. Our goal is to investigate the stellar kinematics of these galaxies and obtain their rotation curves. Additionally, we aim to derive the total stellar and dynamical masses, as well as the maximum rotation velocity, in order to examine their distributions and scaling relations. Using the Pipe3D dataproducts publicly available we obtained the rotation curves, which were fitted considering two components of an axisymmetric Miyamoto–Nagai gravitational potential. We found a wide range of the maximum rotation velocities (117–340 ${\rm km\, s^{-1}}$), with a mean value of 200 ${\rm km\, s^{-1}}$. In addition we found that the total stellar and dynamical masses are in the range of log(M star /M ⊙ ) = 10.1−11.5, with a mean value of log(M star /M ⊙ ) = 10.8, and log(M dyn /M ⊙ ) = 10.4−12.0, with a mean value of log(M dyn /M ⊙ ) = 11.1, respectively. We found a strong correlation between dynamical mass and maximum velocity, between maximum velocity and magnitude, and between stellar mass and maximum velocity. According to these results, barred galaxies exhibit similar behaviour to that of normal spiral galaxies with respect to these relations, as well as in terms of the distribution of their dynamical mass and maximum rotation velocity. However, we found that the distribution of stellar masses of barred galaxies is statistically different from other samples including non-barred galaxies. Finally, analysing the galaxies that show nuclear activity, we find no difference with the rotation curves of normal galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Test and Validate Distributed Coaxial Cable Sensors for in situ Condition Monitoring of Coal-Fired Boiler Tubes

This project aims to test, validate, and advance the technology readiness level (from TRL5 to TRL7) of a novel low-cost distributed stainless-steel/ceramic coaxial cable sensing (SSC-CCS) technology for in situ monitoring of the boiler tube temperature in existing coal-fired power plants. The novel SSC-CCS sensing technology and associated condition-based monitoring (CBM) software to be demonstrated in this project will lead to an improved understanding of the boiler tube failure mechanisms and a prognostic system to improve the overall performance, reliability, and flexibility of the nation’s coal-fired power plant fleet. A boiler tube monitoring system with distributed coaxial cable temperature sensors and a sensor acquisition system was constructed. The high-temperature coaxial cable sensor with a length of 1.3m was made by using a quartz tube (1mm inner diameter (ID) and 6mm outer diameter (OD)) to concentrically separate a 304 stainless-steel (SS) rod (1mm OD) and SS tube (7.94mm OD and 6.16mm ID). The sensor acquisition system includes a vector network analyzer (VNA), a radio frequency (RF) power amplifier, multiple switches and a USB hub. The distributed stainless-steel quartz coaxial cable sensor (SSQ-CCS) had a linear response to temperature with a resolution uncertainty of σ = 0.77℃. To withstand the harsh conditions of 3,300 steam pressures and 800℃ high temperatures, the sensor was shielded by a protective tube made of the same material as the boiler tube. The protection tube had an OD of 1.5 inches and a thickness of 0.25 inches. In the laboratory tests, the sensor showed good sensitivity and fast response. The drift was bounded between +0.33% and -0.67% during a test at 600℃ for 350 hours, indicating good stability of the sensor. A field test was conducted where four sensors were welded on four superheat tubes (SH-Ts) at a coal-fired power station over 400 days. Conventional thermocouples were welded to the superheater tubes alongside the coaxial cable sensors for the purpose of comparison. Two sensors were capable of distributed sensing, with three multiplexed sensing sections. The other two sensors were single section. During the 400-day test period, the power plant experienced startups and shutdowns. At the steady state operations, the temperature of the boiler tube is about 600℃ (1112°F). The sensors recorded the entire coal-firing processes (start-up, steady state, and shut-down) and the glitch event. A GSM modem and a Watchdog were added to the system to ensure reliable data recording. The GSM modem sent daily messages to plant managers and Clemson team to inform the status of the sensor system. If the system was not normally working, the Watchdog would reboot the system automatically. The new coaxial cable based distributed sensing technology has been proven to be successful in both laboratory and field tests. A comprehensive four-stage multi-physics computational framework has been developed to assist the design, optimization, installation, and operation of SSQ-CCS. With the consideration of various operation conditions, we predict the distributions of flue gas temperatures within coal-fired boilers, the temperature correlation between the boiler tube and SSQ-CCS, and the safety of SSQ-CCS. A conditional-based monitoring system is implemented as well. The computational framework developed in this work can guide the future operation of coal-fired plants and other power plants for the safety prediction of boiler operations.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2016

Exploring the Overabundance of ULXs in Metal- and Dust-Poor Local Lyman Break Analogs

We have studied high-mass X-ray binary (HMXB) populations within two low-metallicity, starburst galaxies, Haro 11 and VV 114. These galaxies serve as analogs to high-redshift (z greater than 2) Lyman break galaxies and, within the larger sample of Lyman break analogs (LBAs), they are sufficiently nearby (less than 87 Mpc) to be spatially resolved by Chandra. Previous studies of the X-ray emission in LBAs have found that the 2-10 keV luminosity per star formation rate (SFR) in these galaxies is elevated, potentially because of their low metallicities (12 + log[O/H] = 8.3-8.4). Theoretically, the progenitors of XRBs forming in lower metallicity environments lose less mass from stellar winds over their lifetimes, producing more massive compact objects (i.e., neutron stars and black holes), and thus resulting in more numerous and luminous HMXBs per SFR. In this paper, we have performed an in-depth study of the only two LBAs that have spatially resolved 2-10 keV emission with Chandra to present the bright end of the X-ray luminosity distribution of HMXBs (L(sub X) approximately greater than 10(exp 39) erg s(exp −1); ultraluminous X-ray sources, ULXs) in these low-metallicity galaxies, based on eight detected ULXs. Compared with the star-forming galaxy X-ray luminosity function (XLF) presented by Mineo et al., Haro 11 and VV 114 host approximately equal to 4 times more L(sub X) greater than 10(exp 40) erg s(exp −1) sources than expected given their SFRs. We simulate the effects of source blending from crowded lower-luminosity HMXBs using the star-forming galaxy XLF and then vary the XLF normalizations and bright-end slopes until we reproduce the observed point source luminosity distributions. We find that these LBAs have a shallower bright-end slope (gamma(sub 2) = 1.90) than the standard XLF (gamma(sub 2) 2.73). If we conservatively assume that the brightest X-ray source from each galaxy is powered by an accreting supermassive black hole rather than an HMXB and eliminate these sources from consideration, the luminosity distribution becomes poorly constrained but does appear to be consistent with a standard XLF.

galaxies: evolution aEuro" galaxies: individual↗
DOE OSTI2018

Intensity Distribution of the Three-Wave Diffraction from Dislocation Epitaxial Layers in the Reciprocal Space

The three-wave X-ray diffraction in strongly disordered epitaxial layers of GaN and ZnO is experimentally investigated. The charts of the intensity distribution in the reciprocal space are plotted in coordinates q{sub θ} and q{sub ϕ} for the most intensive three-wave combination (1010)/(1011) by means of subsequent θ- and ϕ-scanning. A nontrivial shape of the θ-sections of these contours at a distance from the ϕ center of reflection is revealed; it is different for different samples. For the θ-curves at the center of reflection, we observed a common peak that may be approximated by the Voigt function with a power-low decrease in the intensity at the wings; the decrease law (from–4.5 to–5.0) is found to be considerably greater than that for the similar curves of two-wave diffraction and not depending on the dislocation density and distribution in layers. In some films we observed a coarse-block structure; in addition, it follows from the distribution in the reciprocal space that these blocks are turned with respect to each other around a normal to the surface, which allows us to suggest the existence of low-angle boundaries between them, consisting exclusively of edge dislocations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1990

Pressure and current balance conditions during electron beam injections from spacecraft

Electrostatic charging level of a conducting surface in response to injections of electron beams into space plasma is investigated by means of one-dimensional Vlasov code. Injections of Maxwellian beams into a vacuum shows that the surface can charge up to an electric potential phi sub s greater than W sub b, where W sub b is the average electron beam energy. Since Maxwellian beams have extended trails with electrons having energies greater than W sub b, it is difficult to quantify the charging level in terms of the energies of the injected electrons. In order to quantitatively understand the charging in excess of W sub b, simulations were carried out for water-bag types of beam with velocity distribution functions described by f(V) = A for V sub min approx. less than V approx. less than V sub max and f(V) = O otherwise, where A is a constant making the normalized beam density unity. It is found that V sub max does not directly determine the charging level. The pressure distribution in the electron sheath determines the electric field distribution near the surface. The electric field in turn determines the electrostatic potential of the vehicle. The pressure distribution is determined by the beam parameters such as the average beam velocity and the velocity spread of the beam.

Hwang, K. S.↗
NASA NTRS1990

A solution adaptive grid procedure for an upwind parabolized flow solver

A solution adaptive grid procedure based on an error equi-distribution scheme is developed and applied to a Parabolized Navier-Stokes solver. An improved method for selecting weighting functions is introduced which involves normalizing a combination of flowfield gradients and curvature of a number of dependent variables and then selecting the largest at each point. The scheme re-distributes grid points line-by-line, with grid point motion controlled by forces analogous to tensional and torsional spring forces with the spring constants set equal to the weighting functions. Torsional terms are functions of the grid point positions along neighboring grid lines and provide grid smoothness and stability. A grid-fitting scheme is introduced for external flows in which the number of grid points in the freestream are reduced to a minimum. Results for several problems are presented to demonstrate the improvements obtainable with the solution adaptive grid procedure.

Harvey, Albert D.↗
DOE OSTI2022

Spatial patterns of snow distribution in the sub-Arctic

Abstract. The spatial distribution of snow plays a vital role in sub-Arctic and Arctic climate, hydrology, and ecology due to its fundamental influence on the water balance, thermal regimes, vegetation, and carbon flux. However, the spatial distribution of snow is not well understood, and therefore, it is not well modeled, which can lead to substantial uncertainties in snow cover representations. To capture key hydro-ecological controls on snow spatial distribution, we carried out intensive field studies over multiple years for two small (2017–2019; ∼ 2.5 km2) sub-Arctic study sites located on the Seward Peninsula of Alaska. Using an intensive suite of field observations (> 22 000 data points), we developed simple models of the spatial distribution of snow water equivalent (SWE) using factors such as topographic characteristics, vegetation characteristics based on greenness (normalized different vegetation index, NDVI), and a simple metric for approximating winds. The most successful model was random forest, using both study sites and all years, which was able to accurately capture the complexity and variability of snow characteristics across the sites. Approximately 86 % of the SWE distribution could be accounted for, on average, by the random forest model at the study sites. Factors that impacted year-to-year snow distribution included NDVI, elevation, and a metric to represent coarse microtopography (topographic position index, TPI), while slope, wind, and fine microtopography factors were less important. The characterization of the SWE spatial distribution patterns will be used to validate and improve snow distribution modeling in the Department of Energy's Earth system model and for improved understanding of hydrology, topography, and vegetation dynamics in the sub-Arctic and Arctic regions of the globe.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

THE IMACS CLUSTER BUILDING SURVEY. IV. THE LOG-NORMAL STAR FORMATION HISTORY OF GALAXIES

We present here a simple model for the star formation history (SFH) of galaxies that is successful in describing both the star formation rate density (SFRD) over cosmic time, as well as the distribution of specific star formation rates (sSFRs) of galaxies at the current epoch, and the evolution of this quantity in galaxy populations to a redshift of z = 1. We show first that the cosmic SFRD is remarkably well described by a simple log-normal in time. We next postulate that this functional form for the ensemble is also a reasonable description for the SFHs of individual galaxies. Using the measured sSFRs for galaxies at z {approx} 0 from Paper III in this series, we then construct a realization of a universe populated by such galaxies in which the parameters of the log-normal SFH of each galaxy are adjusted to match the sSFRs at z {approx} 0 as well as fitting, in ensemble, the cosmic SFRD from z = 0 to z = 8. This model predicts, with striking fidelity, the distribution of sSFRs in mass-limited galaxy samples to z = 1; this match is not achieved by other models with a different functional form for the SFHs of individual galaxies, but with the same number of degrees of freedom, suggesting that the log-normal form is well matched to the likely actual histories of individual galaxies. We also impose the sSFR versus mass distributions at higher redshifts from Paper III as constraints on the model, and show that, as previously suggested, some galaxies in the field, particularly low mass galaxies, are quite young at intermediate redshifts. As emphasized in Paper III, starbursts are insufficient to explain the enhanced sSFRs in intermediate redshift galaxies; we show here that a model using only smoothly varying log-normal SFHs for galaxies, which allows for some fraction of the population to have peak star formation at late times, does however fully explain the observations. Finally, we show that this model, constrained in detail only at redshifts z < 1, also produces the main sequence of star-formation observed at 1.5 < z < 2.5, again suggesting that the log-normal SFHs are a close approximation to the actual histories of typical galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Measurement of ν ¯ μ charged-current single π - production on hydrocarbon in the few-GeV region using MINERvA

The antineutrino scattering channel $\bar{\nu}_{\mu} \,\text{CH} \rightarrow \mu^{+} \,\pi^{-} \,X$(nucleon(s)) is analyzed in the incident energy range 1.5 to 10 GeV using the MINERvA detector at Fermilab. Differential cross sections are reported as functions of $\mu^{+}$ momentum and production angle, $\pi^{-}$ kinetic energy and production angle, and antineutrino energy and squared four-momentum transfer. Distribution shapes are generally reproduced by simulations based on the GENIE, NuWro, and GiBUU event generators, however GENIE (GiBUU) overestimates (underestimates) the cross-section normalizations by 8% (10%). Comparisons of data with the GENIE-based reference simulation probe conventional treatments of cross sections and pion intranuclear rescattering. The distribution of non-track vertex energy is used to decompose the signal sample into reaction categories, and cross sections are determined for the exclusive reactions $\mu^{+} \pi^{-} n$ and $ \mu^+ \pi^{-} p$. A similar treatment applied to the published MINERvA sample $\bar{\nu}_{\mu} \,\text{CH} \rightarrow \mu^{+} \,\pi^{0} \,X$(nucleon(s)) has determined the $\mu^{+} \pi^{0} n$ cross section, and the latter is used with $\sigma(\pi^{-} n)$ and $\sigma(\pi^{-} p)$ to carry out an isospin decomposition of $\bar{\nu}_{\mu}$-induced CC($\pi$). The ratio of magnitudes and relative phase for isospin amplitudes $A_{3}$ and $A_{1}$ thereby obtained are: $R^{\bar{\nu}} = 0.99 \pm 0.19$ and $\phi^{\bar{\nu}} = 93^{\circ} \pm 7^{\circ}$. Our results are in agreement with bubble chamber measurements made four decades ago.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Measuring ion velocity distribution functions through high-aspect ratio holes in inductively coupled plasmas

Several issues associated with plasma etching of high aspect ratio structures originate from the ions' bombardment of the sidewalls of the feature. The off normal angle incident ions are primarily due to their temperature at the sheath edge and possibly to charging effects. We have measured the ion velocity distribution function (IVDF) at the wafer surface in an industrial inductively coupled plasma reactor by using multigrid retarding field analyzers (RFA) in front of which we place 400 μm thick capillary plates with holes of 25, 50, and 100 μm diameters. The RFA then probes IVDF at the exit of the holes with Aspect Ratios (AR) of 16, 8, and 4, respectively. The results show that the ion flux dramatically drops with the increase in AR. By comparing the measured IVDF with an analytical model, we concluded that the ion temperature is 0.27 eV in our plasma conditions. The charging effects are also observed and are shown to significantly reduce the ion energy at the bottom of the feature but only with a “minor” effect on the ion flux and the shape of the IVDF.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Streamwise flow and heat transfer distributions for jet array impingement with crossflow

Two-dimensional arrays of circular jets of air impinging on a heat transfer surface parallel to the jet orifice plate are considered. The air, after impingement, is constrained to exit in a single direction along the channel formed by the surface and the jet plate. The downstream jets are subjected to a crossflow originating from the upstream jets. Experimental and theoretical results obtained for streamwise distributions of jet and crossflow velocities are presented and compared. Measured Nusselt numbers resolved to one streamwise hole spacing are correlated with individual spanwise row jet Reynolds numbers and crossflow-to-jet velocity ratios. Correlations are presented for both inline and staggered hole patterns including effects of geometric parameters: streamwise hole spacing, spanwise hole spacing, and channel height, normalized by hole diameter. The physical mechanisms influencing heat transfer coefficients as a function of flow distribution and geometric parameters are also discussed.

Florschuetz, L. W.↗
DOE OSTI2016

Determination of Interfacial Mixing in Tapered Block Polymer Thin Films: Experimental and Theoretical Investigations

Tapered block polymers are an emerging class of macromolecules with unique and diverse self-assembly behavior and properties. Here, we directly examine the manipulation of self-assembled interfaces in poly(isoprene- b -styrene) (I-S)-based block polymers (BPs) by synthesizing non-tapered (I-S), normal tapered (I-IS-S), and inverse tapered (I-SI-S) BPs with controlled monomer segment distributions. We provide the first direct measurements of interfacial mixing for these tapered polymers through X-ray reflectivity (XRR). The density profiles from XRR are compared to results from fluids density functional theory (fDFT) with good agreement. We find that our normal tapered BPs (30 vol % tapering) have similar interfacial mixing to diblock polymers, while our inverse tapered BPs (30 vol % tapering) have much wider interfaces. Additionally, differential scanning calorimetry (DSC) studies elucidate the influence of tapering on the glass transition temperature (T g ) and change of heat capacity (ΔC P ) for each BP phase, and quantitative analysis from ΔC P also indicates enhanced mixing in the inverse tapered I-SI-S BPs. Finally, we investigate the free surface morphologies of these tapered polymers in thin film geometries. The inverse tapered BP form larger island/hole structures likely due to decreased surface elasticity from mixing as a result of the modified interblock interfacial characteristics. These results demonstrate that BPs with similar molecular weights and compositions can exhibit different thermodynamic properties and free surface morphologies in thin film geometries, as influenced by monomer sequence.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Comparisons of angularly and spectrally resolved Bremsstrahlung measurements to two-dimensional multi-stage simulations of short-pulse laser-plasma interactions

A 2-D multi-stage simulation model incorporating realistic laser conditions and a fully resolved electron distribution handoff has been developed and compared to angularly and spectrally resolved Bremsstrahlung measurements from high-Z planar targets. For near-normal incidence and 0.5-1 × 10{sup 20} W/cm{sup 2} intensity, particle-in-cell (PIC) simulations predict the existence of a high energy electron component consistently directed away from the laser axis, in contrast with previous expectations for oblique irradiation. Measurements of the angular distribution are consistent with a high energy component when directed along the PIC predicted direction, as opposed to between the target normal and laser axis as previously measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Cylindrical Kernel Density Estimators For Monte Carlo Neutron Transport Reactor Physics Problems

Kernel Density Estimators (KDEs) have been investigated as an alternative to histogram tallies in Monte Carlo radiation transport simulations for capturing spatially-resolved quantities such as scalar flux and reaction rate densities. KDEs are non-parametric estimators that show potential for producing smooth estimates of quantities of interest with reduced variance compared to traditional histogram tallies. Furthermore, KDEs obtain estimates of the underlying distribution at user-defined points without requiring an underlying mesh, thus making KDEs an attractive potential alternative to histogram tallies. Recent work has introduced the Mean Free Path based KDE (MFP KDE) in order to capture distributions with large curvature at material interfaces, e.g. absorption reaction rates in UO{sub 2} fuel in light water reactors. However, the multivariate MFP KDE produces inaccuracies in the form of localized spikes in problems with material cross sections that contain large resonances. Since the MFP KDE works well in one dimension, it is hypothesized that a different MFP KDE formulation where the normalization coefficient only contains one factor of the cross section would be able to capture these distributions with large curvature without producing spikes in the estimate. This can be accomplished by using a cylindrical KDE whose 2-D support area is wedge-shaped and applying the MFP KDE to the radial component of the cylindrical KDE. These cylindrical KDEs use 2-D (r,θ) cylindrical coordinates with no dependence in axial dimension z. This paper derives the cylindrical KDE and the cylindrical MFP KDE and tests them on a reactor physics pin-cell benchmark problem, producing accurate estimates of flux and reaction rate densities with differences less than 0.5% when compared to a reference histogram. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

Distribution of Snow and Maximum Snow Water Equivalent Obtained by LANDSAT Data and Degree Day Method

Maximum snow water equivalence and snowcover distribution are estimated using several LANDSAT data taken in snowmelting season over a four year period. The test site is Okutadami-gawa Basin located in the central position of Tohoku-Kanto-Chubu District. The year to year normalization for snowmelt volume computation on the snow line is conducted by year to year correction of degree days using the snowcover percentage within the test basin obtained from LANDSAT data. The maximum snow water equivalent map in the test basin is generated based on the normalized snowmelt volume on the snow line extracted from four LANDSAT data taken in a different year. The snowcover distribution on an arbitrary day in snowmelting of 1982 is estimated from the maximum snow water equivalent map. The estimated snowcover is compared with the snowcover area extracted from NOAA-AVHRR data taken on the same day. The applicability of the snow estimation using LANDSAT data is discussed.

Takeda, K.↗
DOE OSTI2015

Breast dose in mammography is about 30% lower when realistic heterogeneous glandular distributions are considered

Purpose: Current dosimetry methods in mammography assume that the breast is comprised of a homogeneous mixture of glandular and adipose tissues. Three-dimensional (3D) dedicated breast CT (bCT) data sets were used previously to assess the complex anatomical structure within the breast, characterizing the statistical distribution of glandular tissue in the breast. The purpose of this work was to investigate the effect of bCT-derived heterogeneous glandular distributions on dosimetry in mammography. Methods: bCT-derived breast diameters, volumes, and 3D fibroglandular distributions were used to design realistic compressed breast models comprised of heterogeneous distributions of glandular tissue. The bCT-derived glandular distributions were fit to biGaussian functions and used as probability density maps to assign the density distributions within compressed breast models. The MCNPX 2.6.0 Monte Carlo code was used to estimate monoenergetic normalized mean glandular dose “DgN(E)” values in mammography geometry. The DgN(E) values were then weighted by typical mammography x-ray spectra to determine polyenergetic DgN (pDgN) coefficients for heterogeneous (pDgN{sub hetero}) and homogeneous (pDgN{sub homo}) cases. The dependence of estimated pDgN values on phantom size, volumetric glandular fraction (VGF), x-ray technique factors, and location of the heterogeneous glandular distributions was investigated. Results: The pDgN{sub hetero} coefficients were on average 35.3% (SD, 4.1) and 24.2% (SD, 3.0) lower than the pDgN{sub homo} coefficients for the Mo–Mo and W–Rh x-ray spectra, respectively, across all phantom sizes and VGFs when the glandular distributions were centered within the breast phantom in the coronal plane. At constant breast size, increasing VGF from 7.3% to 19.1% lead to a reduction in pDgN{sub hetero} relative to pDgN{sub homo} of 23.6%–27.4% for a W–Rh spectrum. Displacement of the glandular distribution, at a distance equal to 10% of the compressed breast width in the superior and inferior directions, resulted in a 37.3% and a −26.6% change in the pDgN{sub hetero} coefficient, respectively, relative to the centered distribution for the Mo–Mo spectrum. Lateral displacement of the glandular distribution, at a distance equal to 10% of the compressed breast width, resulted in a 1.5% change in the pDgN{sub hetero} coefficient relative to the centered distribution for the W–Rh spectrum. Conclusions: Introducing bCT-derived heterogeneous glandular distributions into mammography phantom design resulted in decreased glandular dose relative to the widely used homogeneous assumption. A homogeneous distribution overestimates the amount of glandular tissue near the entrant surface of the breast, where dose deposition is exponentially higher. While these findings are based on clinically measured distributions of glandular tissue using a large cohort of women, future work is required to improve the classification of glandular distributions based on breast size and overall glandular fraction.

60 APPLIED LIFE SCIENCES↗