Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mixing index”

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 163 records · Page 9

Applicability of the Effective-Medium Approximation to Heterogeneous Aerosol Particles.

The effective-medium approximation (EMA) is based on the assumption that a heterogeneous particle can have a homogeneous counterpart possessing similar scattering and absorption properties. We analyze the numerical accuracy of the EMA by comparing superposition T-matrix computations for spherical aerosol particles filled with numerous randomly distributed small inclusions and Lorenz-Mie computations based on the Maxwell-Garnett mixing rule. We verify numerically that the EMA can indeed be realized for inclusion size parameters smaller than a threshold value. The threshold size parameter depends on the refractive-index contrast between the host and inclusion materials and quite often does not exceed several tenths, especially in calculations of the scattering matrix and the absorption cross section. As the inclusion size parameter approaches the threshold value, the scattering-matrix errors of the EMA start to grow with increasing the host size parameter and or the number of inclusions. We confirm, in particular, the existence of the effective-medium regime in the important case of dust aerosols with hematite or air-bubble inclusions, but then the large refractive-index contrast necessitates inclusion size parameters of the order of a few tenths. Irrespective of the highly restricted conditions of applicability of the EMA, our results provide further evidence that the effective-medium regime must be a direct corollary of the macroscopic Maxwell equations under specific assumptions.

aerosols↗

A Detailed Characterization of the Saharan Dust Collected During the Fennec Campaign in 2011: In Situ Ground-Based and Laboratory Measurements

Millions of tons of mineral dust are lifted by the wind from arid surfaces and transported around the globe every year. The physical and chemical properties of the mineral dust are needed to better constrain remote sensing observations and are of fundamental importance for the understanding of dust atmospheric processes. Ground-based in situ measurements and in situ filter collection of Saharan dust were obtained during the Fennec campaign in the central Sahara in 2011. This paper presents results of the absorption and scattering coefficients, and hence single scattering albedo (SSA), of the Saharan dust measured in real time during the last period of the campaign and subsequent laboratory analysis of the dust samples collected in two supersites, SS1 and SS2, in Algeria and in Mauritania, respectively. The samples were taken to the laboratory, where their size and aspect ratio distributions, mean chemical composition, spectral mass absorption efficiency, and spectral imaginary refractive index were obtained from the ultraviolet (UV) to the near-infrared (NIR) wavelengths. At SS1 in Algeria, the time series of the scattering coefficients during the period of the campaign show dust events exceeding 3500 Mm (exp -1), and a relatively high mean SSA of 0.995 at 670 nm was observed at this site. The laboratory results show for the fine particle size distributions (particles diameter <5μm and mode diameter at2–3 μm) in both sites a spectral dependence of the imaginary part of the refractive index Im(m) with a bow-like shape, with increased absorption in UV as well as in the shortwave infrared. The same signature was not observed, however, in the mixed particle size distribution (particle diameter<10 μm and mode diameter at 4 μm) in Algeria. Im(m) was found to range from 0.011 to 0.001i for dust collected in Algeria and 0.008 to 0.002i for dust collected in Mauritania over the wavelength range of 350-2500 nm. Differences in the mean elemental composition of the dust collected in the supersites in Algeria and in Mauritania and between fine and mixed particle size distributions were observed from EDXRF measurements, although those differences cannot be used to explain the optical properties variability between the samples. Finally, particles with low-density typically larger than 10 μm in diameter were found in some of the samples collected at the supersite in Mauritania, but these low-density particles were not observed in Algeria.

Saharan dust measured in real time↗

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work uses an automated grid adaptation tool recently implemented in VULCAN-CFD, and explores this tool’s ability to solve highspeed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user-generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flowfield parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Furthermore, motivated by the grid convergence index (GCI) used for structured grid simulations, a grid convergence estimate (GCE) was developed and demonstrated for the grid adaptation. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex high-speed mixing problems with only a fraction of the grid cells required to obtain comparable results on a user-generated structured grid.

hypersonics↗

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work uses an automated grid adaptation tool recently implemented in VULCAN-CFD, and explores this tool’s ability to solve highspeed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user-generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flowfield parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Furthermore, motivated by the grid convergence index (GCI) used for structured grid simulations, a grid convergence estimate (GCE) was developed and demonstrated for the grid adaptation. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex high-speed mixing problems with only a fraction of the grid cells required to obtain comparable results on a user-generated structured grid.

hypersonics↗

Towards a Regolith Maturity Index for Howardites

The Dawn spacecraft has just arrived at asteroid 4 Vesta, parent of the howardite, eucrite and diogenite (HED) meteorites [1], to begin a yearlong surface study from orbit [2]. As Dawn will view a debris-covered surface, understanding the formation and mixing processes for the debris layer will strongly aid surface data interpretations. Howardites are polymict breccias mainly composed of clasts derived from basaltic (eucritic) and orthopy-roxenitic (diogenitic) parent materials [3]. Some howardites are poorly reworked (fragmental howardites) whilst others have been extensively gardened in an active regolith (regolithic howardites) [4]. The latter may represent an ancient, well-mixed regolith, whilst the former may be from more recent ejecta deposits [4]. Due to environmental differences, regolith development on Vesta differs in detail from that on the Moon [4-6]. We have been developing petrological criteria to apply to howardite thin sections to determine their relative regolithic maturity, which we are fine-tuning with comparison to noble gas data [7, 8]. Whilst we previously emphasized the abundance of reworked clasts (fragmental and impact-melt breccia clasts), this is an imperfect criterion: one howardite with abundant re-worked clasts (EET 99408) shows no evidence of solar wind Ne (SW-Ne), yet, two of our alleged fragmental howardites have clear SW-Ne signatures (LEW 85313, MET 00423) [7, 8]. We are now investigating the diversity in minor and trace element contents of low-Ca pyroxene clasts in howardites as a measure of regolith grade, and will begin analyses of such grains within reworked clasts. Our hypothesis is that regolithic howardites (and the breccia clasts they contain) will show greater diversity because they sampled more diverse diogenitic plutons than fragmental howardites, which formed from ejecta from only a few impacts [e.g. 4]. Our initial LA-ICP-MS work showed ranges in trace element diversity in low-Ca pyroxenes (estimated from the standard error of the mean of analyses), where those howardites considered of medium to high regolithic grade showed greater diversity [9]. Our EMPA results (from a larger howardite suite) show an overall greater diversity in our putative medium to high regolithic grade howardites, though there are exceptions. The greatest diversity is found for paired howardites GRO 95574 and GRO 95581, which were not considered regolithic in our initial study. We will continue investigating avenues to determine regolith maturity in thin section, factoring in bulk rock compositional data, and will coordinate these studies with noble gas results.

Mittlefehldt, David W.↗

Evaluating Model Parameterizations of Submicron Aerosol Scattering and Absorption with in situ Data from ARCTAS 2008

Accurate modeling of the scattering and absorption of ultraviolet and visible radiation by aerosols is essential for accurate simulations of atmospheric chemistry and climate. Closure studies using in situ measurements of aerosol scattering and absorption can be used to evaluate and improve models of aerosol optical properties without interference from model errors in aerosol emissions, transport, chemistry, or deposition rates. Here we evaluate the ability of four externally mixed, fixed size distribution parameterizations used in global models to simulate submicron aerosol scattering and absorption at three wavelengths using in situ data gathered during the 2008 Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) campaign. The four models are the NASA Global Modeling Initiative (GMI) Combo model, GEOS-Chem v9- 02, the baseline configuration of a version of GEOS-Chem with online radiative transfer calculations (called GC-RT), and the Optical Properties of Aerosol and Clouds (OPAC v3.1) package. We also use the ARCTAS data to perform the first evaluation of the ability of the Aerosol Simulation Program (ASP v2.1) to simulate submicron aerosol scattering and absorption when in situ data on the aerosol size distribution are used, and examine the impact of different mixing rules for black carbon (BC) on the results. We find that the GMI model tends to overestimate submicron scattering and absorption at shorter wavelengths by 10-23 percent, and that GMI has smaller absolute mean biases for submicron absorption than OPAC v3.1, GEOS-Chem v9-02, or GC-RT. However, the changes to the density and refractive index of BC in GCRT improve the simulation of submicron aerosol absorption at all wavelengths relative to GEOS-Chem v9-02. Adding a variable size distribution, as in ASP v2.1, improves model performance for scattering but not for absorption, likely due to the assumption in ASP v2.1 that BC is present at a constant mass fraction throughout the aerosol size distribution. Using a core-shell mixing rule in ASP overestimates aerosol absorption, especially for the fresh biomass burning aerosol measured in ARCTAS-B, suggesting the need for modeling the time-varying mixing states of aerosols in future versions of ASP.

Alvarado, Matthew J.↗

Phase separation kinetics in immiscible liquids

The kinetics of phase separation in the succinonitrile-water system are being investigated. Experiments involve initial physical mixing of the two immiscible liquids at a temperature above the consolute, decreasing the temperature into the miscibility gap, followed by iamging of the resultant microstructure as it evolves with time. Refractive index differences allow documentation of the changing microstructures by noninvasive optical techniques without the need to quench the liquid structures for analysis.

Ng, Lee H.↗

Phase separation kinetics in immiscible liquids

The kinetics of phase separation in the succinonitrile-water system are being investigated. Experiments involve initial physical mixing of the two immiscible liquids at a temperature above the consolute, decreasing the temperature into the miscibility gap, followed by imaging of the resultant microstructure as it evolves with time. Refractive index differences allow documentation of the changing microstructures by noninvasive optical techniques without the need to quench the liquid structures for analysis.

Sadoway, D. R.↗

Linear mixing model applied to AVHRR LAC data

A linear mixing model was applied to coarse spatial resolution data from the NOAA Advanced Very High Resolution Radiometer. The reflective component of the 3.55 - 3.93 microns channel was extracted and used with the two reflective channels 0.58 - 0.68 microns and 0.725 - 1.1 microns to run a Constraine Least Squares model to generate vegetation, soil, and shade fraction images for an area in the Western region of Brazil. The Landsat Thematic Mapper data covering the Emas National park region was used for estimating the spectral response of the mixture components and for evaluating the mixing model results. The fraction images were compared with an unsupervised classification derived from Landsat TM data acquired on the same day. The relationship between the fraction images and normalized difference vegetation index images show the potential of the unmixing techniques when using coarse resolution data for global studies.

Holben, Brent N.↗

Validation of an All-Pressure Fluid Drop Model: Heptane Fluid Drops in Nitrogen

Despite the fact that supercritical fluids occur both in nature and in industrial situations, the fundamentals of their behavior is poorly understood because supercritical fluids combine the characteristics of both liquids and gases, and therefore their behavior is not intuitive. There are several specific reasons for the lack of understanding: First, data from (mostly optical) measurements can be very misleading because regions of high density thus observed are frequently identified with liquids. A common misconception is that if in an experiment one can optically identify "drops" and "ligaments", the observed fluid must be in a liquid state. This inference is incorrect because in fact optical measurements detect any large change (i.e. gradients) in density. Thus, the density ratio may be well below Omicron(10(exp 3)) that characterizes its liquid/gas value, but the measurement will still identify a change in the index of refraction providing that the change is sudden (steep gradients). As shown by simulations of supercritical fluids, under certain conditions the density gradients may remain large during the supercritical binary fluids mixing, thus making them optically identifiable. Therefore, there is no inconsistency between the optical observation of high density regions and the fluids being in a supercritical state. A second misconception is that because a fluid has a liquid-like density, it is appropriate to model it as a liquid. However, such fluids may have liquid-like densities while their transport properties differ from those of a liquid. Considering that the critical pressure of most fuel hydrocarbons used in Diesel and gas turbine engines is in the range of 1.5 - 3 MPa, and the fact that the maximum pressure attained in these engines is about 6 Mps, it is clear that the fuel in the combustion chamber will experience both subcritical and supercritical conditions. Studies of drop behavior over a wide range of pressures were performed in the past, however none of these studies identified the crucial differences between the subcritical and supercritical behavior. In fact, in two of these studies, it was found that the subcritical and supercritical behavior is similar as the drop diameter decreased according to the classical d(exp 2)-law over a wide range of pressures and drop diameters. The present study is devoted to the exploration of differences in fluid-behavior characteristics under subcritical and supercritical conditions in the particular case of heptane fluid drops in nitrogen; these substances were selected because of the availability of experimental observations for model validation.

Harstad, K.↗

Observations of Ocean Primary Productivity Using MODIS

Measuring the magnitude and variability of oceanic net primary productivity (NPP) represents a key advancement toward our understanding of the dynamics of marine ecosystems and the role of the ocean in the global carbon cycle. MODIS observations make two new contributions in addition to continuing the bio-optical time series begun with Orbview-2's SeaWiFS sensor. First, MODIS provides weekly estimates of global ocean net primary productivity on weekly and annual time periods, and annual empirical estimates of carbon export production. Second, MODIS provides additional insight into the spatial and temporal variations in photosynthetic efficiency through the direct measurements of solar-stimulated chlorophyll fluorescence. The two different weekly productivity indexes (first developed by Behrenfeld & Falkowski and by Yoder, Ryan and Howard) are used to derive daily productivity as a function of chlorophyll biomass, incident daily surface irradiance, temperature, euphotic depth, and mixed layer depth. Comparisons between these two estimates using both SeaWiFS and MODIS data show significant model differences in spatial distribution after allowance for the different integration depths. Both estimates are strongly dependence on the accuracy of the chlorophyll determination. In addition, an empirical approach is taken on annual scales to estimate global NPP and export production. Estimates of solar stimulated fluorescence efficiency from chlorophyll have been shown to be inversely related to photosynthetic efficiency by Abbott and co-workers. MODIS provides the first global estimates of oceanic chlorophyll fluorescence, providing an important proof of concept. MODIS observations are revealing spatial patterns of fluorescence efficiency which show expected variations with phytoplankton photo-physiological parameters as measured during in-situ surveys. This has opened the way for research into utilizing this information to improve our understanding of oceanic NPP variability. Deriving the ocean bio-optical properties places severe demands on instrument performance (especially band to band precision) and atmospheric correction. Improvements in MODIS instrument characterization and calibration over the first 16 mission months have greatly improved the accuracy of the chlorophyll input fields and FLH, and therefore the estimates of NPP and fluorescence efficiency. Annual estimates now show the oceanic NPP accounts for 40-50% of the global total NPP, with significant interannual variations related to large scale ocean processes. Spatial variations in ocean NPP, and exported production, have significant effects on exchange of CO2 between the ocean and atmosphere. Further work is underway to improve both the primary productivity model functions, and to refine our understanding of the relationships between fluorescence efficiency and NPP estimates. We expect that the MODIS instruments will prove extremely useful in assessing the time dependencies of oceanic carbon uptake and effects of iron enrichment, within the global carbon cycle.

Esaias, Wayne E.↗

Multi-spatio-temporal scales PIV in a turbulent buoyant jet discharging in a linearly stratified environment

Time-resolved particle image velocimetry is implemented with a camera array and several laser sheets; this results in a multi-spatio-temporal scale system that is modular and flexible. The setup is optimized to capture global flow features, while locally resolving in space and time near the Kolmogorov scale. The apparatus relies extensively on machine vision CMOS cameras; they are inexpensive and stream to computer hard drive with virtually continuous time-resolved records (up to one hour for the current system). This leads to statistically converged data and also helps in minimizing down time between experimental runs. Details of the implementation and design of experiment are reported. The system instruments a vertical buoyant jet discharging in a linearly stratified environment. Refractive index matched aqueous solutions enable precise optical deployment. The density difference is 3% and the fluids have similar dynamic viscosity. The jet Reynolds number is 2.00 x 10 4 and is above the mixing transition. Such flows are typically challenging to instrument and few velocity field data are available to date. The velocity statistics are reported as well as first insights gained from the campaign.

42 ENGINEERING↗

Wind shear and wet and dry thermodynamic indices as predictors of thunderstorm motion and severity and application to the AVE 4 experimental data

Two types of parameters are computed and mapped for use in assessing their individual merits as predictors of occurrence and severity of thunderstorms. The first group is comprised of equivalent potential temperature, potential temperature, water vapor mixing ratio, and wind speed. Equivalent potential temperature maxima and strong gradients of equivalent potential temperature at the surface correlate well with regions of thunderstorm activity. The second type, comprised of the energy index, shear index, and energy shear index, incorporates some model dynamics of thunderstorms, including nonthermodynamic forcing. The energy shear index is found to improve prediction of tornadic and high-wind situations slightly better than other indices. It is concluded that further development and refinement of nonthermodynamic aspects of predictive indices are definitely warranted.

Connell, J. R.↗

Biocybernetic system evaluates indices of operator engagement in automated task

A biocybernetic system has been developed as a method to evaluate automated flight deck concepts for compatibility with human capabilities. A biocybernetic loop is formed by adjusting the mode of operation of a task set (e.g., manual/automated mix) based on electroencephalographic (EEG) signals reflecting an operator's engagement in the task set. A critical issue for the loop operation is the selection of features of the EEG to provide an index of engagement upon which to base decisions to adjust task mode. Subjects were run in the closed-loop feedback configuration under four candidate and three experimental control definitions of an engagement index. The temporal patterning of system mode switching was observed for both positive and negative feedback of the index. The indices were judged on the basis of their relative strength in exhibiting expected feedback control system phenomena (stable operation under negative feedback and unstable operation under positive feedback). Of the candidate indices evaluated in this study, an index constructed according to the formula, beta power/(alpha power + theta power), reflected task engagement best.

Cerebral Cortex/physiology↗

High-energy Emission from the Composite Supernova Remnant MSH 15-56

MSH 1556 (G326.3-1.8) is a composite supernova remnant (SNR) that consists of an SNR shell and a displaced pulsar wind nebula (PWN) in the radio. We present XMM-Newton and Chandra X-ray observations of the remnant that reveal a compact source at the tip of the radio PWN and complex structures that provide evidence for mixing of the supernova (SN) ejecta with PWN material following a reverse shock interaction. The X-ray spectra are well fitted by a non-thermal power-law model whose photon index steepens with distance from the presumed pulsar, and a thermal component with an average temperature of 0.55 keV. The enhanced abundances of silicon and sulfur in some regions, and the similar temperature and ionization timescale, suggest that much of the X-ray emission can be attributed to SN ejecta that have either been heated by the reverse shock or swept up by the PWN. We find one region with a lower temperature of 0.3 keV that appears to be in ionization equilibrium.Assuming the Sedov model, we derive a number of SNR properties, including an age of 16,500 yr. Modeling of the gamma-ray emission detected by Fermi shows that the emission may originate from the reverse shock-crushed PWN.

Temim, Tea↗

BOREAS AFM-1 NOAA/ATDD Long-EZ Aircraft Flux data Over the SSA

This data set contains measurements from the Airborne Flux and Meteorology (AFM)-1 National Oceanographic and Atmospheric Administration/Atmospheric Turbulence and Diffusion Division (NOAA/ATDD) Long-EZ Aircraft collected during the 1994 Intensive Field Campaigns (IFCs) at the southern study area (SSA). These measurements were made from various instruments mounted on the aircraft. The data that were collected include aircraft altitude, wind direction, wind speed, air temperature, potential temperature, water mixing ratio, U and V components of wind velocity, static pressure, surface radiative temperature, downwelling and upwelling total radiation, downwelling and upwelling longwave radiation, net radiation, downwelling and upwelling photosynthectically active radiation (PAR), greenness index, CO2 concentration, O3 concentration, and CH4 concentration. There are also various columns that indicate the standard deviation, skewness, kurtosis, and trend of some of these data. The data are stored in tabular ASCII files. The NOAA/ATDD Long-EZ aircraft flux data are available from the Earth Observing System Data and Information System (EOSDIS) Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center (DAAC). The data files are available on a CD-ROM (see document number 20010000884).

Crawford, Timothy L.↗

Mapping Temperate Forest Phenology Using Tower, UAV, and Ground-Based Sensors

Phenology is a distinct marker of the impacts of climate change on ecosystems. Accordingly, monitoring the spatiotemporal patterns of vegetation phenology is important to understand the changing Earth system. A wide range of sensors have been used to monitor vegetation phenology, including digital cameras with different viewing geometries mounted on various types of platforms. Sensor perspective, view-angle, and resolution can potentially impact estimates of phenology. We compared three different methods of remotely sensing vegetation phenology—an unoccupied aerial vehicle (UAV)-based, downward-facing RGB camera, a below-canopy, upward-facing hemispherical camera with blue (B), green (G), and near-infrared (NIR) bands, and a tower-based RGB PhenoCam, positioned at an oblique angle to the canopy—to estimate spring phenological transition towards canopy closure in a mixed-species temperate forest in central Virginia, USA. Our study had two objectives: (1) to compare the above- and below-canopy inference of canopy greenness (using green chromatic coordinate and normalized difference vegetation index) and canopy structural attributes (leaf area and gap fraction) by matching below-canopy hemispherical photos with high spatial resolution (0.03 m) UAV imagery, to find the appropriate spatial coverage and resolution for comparison; (2) to compare how UAV, ground-based, and tower-based imagery performed in estimating the timing of the spring phenological transition. We found that a spatial buffer of 20 m radius for UAV imagery is most closely comparable to below-canopy imagery in this system. Sensors and platforms agree within +/− 5 days of when canopy greenness stabilizes from the spring phenophase into the growing season. We show that pairing UAV imagery with tower-based observation platforms and plot-based observations for phenological studies (e.g., long-term monitoring, existing research networks, and permanent plots) has the potential to scale plot-based forest structural measures via UAV imagery, constrain uncertainty estimates around phenophases, and more robustly assess site heterogeneity.

Jeff W. Atkins↗