Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical simulation”

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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Hybrid Reynolds-Averaged/Large-Eddy Simulations of a Coaxial Supersonic Free-Jet Experiment

Reynolds-averaged and hybrid Reynolds-averaged/large-eddy simulations have been applied to a supersonic coaxial jet flow experiment. The experiment was designed to study compressible mixing flow phenomenon under conditions that are representative of those encountered in scramjet combustors. The experiment utilized either helium or argon as the inner jet nozzle fluid, and the outer jet nozzle fluid consisted of laboratory air. The inner and outer nozzles were designed and operated to produce nearly pressure-matched Mach 1.8 flow conditions at the jet exit. The purpose of the computational effort was to assess the state-of-the-art for each modeling approach, and to use the hybrid Reynolds-averaged/large-eddy simulations to gather insight into the deficiencies of the Reynolds-averaged closure models. The Reynolds-averaged simulations displayed a strong sensitivity to choice of turbulent Schmidt number. The initial value chosen for this parameter resulted in an over-prediction of the mixing layer spreading rate for the helium case, but the opposite trend was observed when argon was used as the injectant. A larger turbulent Schmidt number greatly improved the comparison of the results with measurements for the helium simulations, but variations in the Schmidt number did not improve the argon comparisons. The hybrid Reynolds-averaged/large-eddy simulations also over-predicted the mixing layer spreading rate for the helium case, while under-predicting the rate of mixing when argon was used as the injectant. The primary reason conjectured for the discrepancy between the hybrid simulation results and the measurements centered around issues related to the transition from a Reynolds-averaged state to one with resolved turbulent content. Improvements to the inflow conditions were suggested as a remedy to this dilemma. Second-order turbulence statistics were also compared to their modeled Reynolds-averaged counterparts to evaluate the effectiveness of common turbulence closure assumptions.

Baurle, Robert A.↗
NASA NTRS2016

The XXL Survey I. Scientific Motivations - Xmm-Newton Observing Plan - Follow-up Observations and Simulation Programme

The quest for the cosmological parameters that describe our universe continues to motivate the scientific community to undertake very large survey initiatives across the electromagnetic spectrum. Over the past two decades, the Chandra and XMM-Newton observatories have supported numerous studies of X-ray-selected clusters of galaxies, active galactic nuclei (AGNs), and the X-ray background. The present paper is the first in a series reporting results of the XXL-XMM survey; it comes at a time when the Planck mission results are being finalized. Aims. We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of approx. 5 x 10(exp 15) erg/s/sq cm in the [0.5-2] keV band (completeness limit). The surveys main goals are to provide constraints on the dark energy equation of state from the space-time-distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programs. Methods. We describe the 542 XMM observations along with the associated multi- and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. Results. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. Conclusions. The XXL multi- data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve asa calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps.

X-rays: general – large-scale structure of Unive↗
NASA NTRS2016

Optical and Near-Infrared Observations of SN 2013DX Associated with GRB 130702A

We present optical and near-infrared (NIR) light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be E(sub gamma, iso) = 6.4(+1.3/-1.0) x 10(exp 50) erg (1 keV to 10 MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed g'r'i'z' light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves approx. 20% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined SNe Ic, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of approx. 21,000 km/s. We construct a quasi-bolometric (g'r'z'yJ) light curve for SN 2013dx, only the fifth GRB-associated SN with extensive NIR coverage and the third with a bolometric light curve extending beyond (Delta)t > 40 days. Together with the measured photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a Ni-56 mass of M(sub Ni) = 0.37+/- 0.01 Stellar Mass, an ejecta mass of M(sub ej) = 3.1+/- 0.1 Stellar Mass, and a kinetic energy of E(sub K) = (8.2+/- 0.43) x 10(exp 51) erg (statistical uncertainties only), consistent with previous GRB-associated supernovae. When considering the ensemble population of GRB-associated supernovae, we find no correlation between the mass of synthesized Ni-56 and high-energy properties, despite clear predictions from numerical simulations that M(sub Ni) should correlate with the degree of asymmetry. On the other hand, M(sub Ni) clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.

gamma-ray burst: individual (GRB 130702A) – supe↗
NASA NTRS2017

Future ATM Concepts Evaluation Tool (FACET) Interface Control Document

This Interface Control Document (ICD) documents the airspace adaptation and air traffic inputs of NASA's Future ATM Concepts and Evaluation Tool (FACET). Its intended audience is the project manager, project team, development team, and stakeholders interested in interfacing with the system. FACET equips Air Traffic Management (ATM) researchers and service providers with a way to explore, develop and evaluate advanced air transportation concepts before they are field-tested and eventually deployed. FACET is a flexible software tool that is capable of quickly generating and analyzing thousands of aircraft trajectories. It provides researchers with a simulation environment for preliminary testing of advanced ATM concepts. Using aircraft performance profiles, airspace models, weather data, and flight schedules, the tool models trajectories for the climb, cruise, and descent phases of flight for each type of aircraft. An advanced graphical interface displays traffic patterns in two and three dimensions, under various current and projected conditions for specific airspace regions or over the entire continental United States. The system is able to simulate a full day's dynamic national airspace system (NAS) operations, model system uncertainty, measure the impact of different decision-makers in the NAS, and provide analysis of the results in graphical form, including sector, airport, fix, and airway usage statistics. NASA researchers test and analyze the system-wide impact of new traffic flow management algorithms under anticipated air traffic growth projections on the nation's air traffic system. In addition to modeling the airspace system for NASA research, FACET has also successfully transitioned into a valuable tool for operational use. Federal Aviation Administration (FAA) traffic flow managers and commercial airline dispatchers have used FACET technology for real-time operations planning. FACET integrates live air traffic data from FAA radar systems and weather data from the National Weather Service to summarize NAS performance. This information allows system operators to reroute flights around congested airspace and severe weather to maintain safety and minimize delay. FACET also supports the planning and post-operational evaluation of reroute strategies at the national level to maximize system efficiency. For the commercial airline passenger, strategic planning with FACET can result in fewer flight delays and cancellations. The performance capabilities of FACET are largely due to its architecture, which strikes a balance between flexibility and fidelity. FACET is capable of modeling the airspace operations for the continental United States, processing thousands of aircraft on a single computer. FACET was written in Java and C, enabling the portability of its software to a variety of operating systems. In addition, FACET was designed with a modular software architecture to facilitate rapid prototyping of diverse ATM concepts. Several advanced ATM concepts have already been implemented in FACET, including aircraft self-separation, prediction of aircraft demand and sector congestion, system-wide impact assessment of traffic flow management constraints, and wind-optimal routing.

FACET↗
NASA NTRS2020

The Prevalence of Precipitation from Polar Supercooled Clouds

Supercooled clouds substantially impact polar surface energy budgets but large-scale models often underestimate their occurrence, which motivates accurately establishing metrics of basic processes. An analysis of long-term measurements at Utqiaġvik, Alaska, and McMurdo Station, Antarctica, combines lidar-validated use of soundings to identify supercooled cloud layers and colocated ground-based profiling radar measurements to quantify cloud base precipitation. We find that more than 85 % (75 %) of sampled supercooled layers are precipitating over the Arctic (Antarctic) site, with more than 75 % (50 %) precipitating continuously to the surface. Such high frequencies can be reconciled with substantially lesser spaceborne estimates by considering differences in radar hydrometeor detection sensitivity. While ice precipitation into supercooled clouds from aloft is common, we also find that the great majority of supercooled cloud layers without ice falling into them are themselves continuously generating precipitation. Such sustained primary ice formation is consistent with continuous activation of immersion-mode ice nucleating particles (INPs), suggesting that supercooled cloud formation is a principal gateway to ice formation at temperatures greater than ~ −38 °C over polar regions. The prevalence of weak precipitation fluxes is also consistent with supercooled cloud longevity, and with well-observed and widely simulated case studies. An analysis of colocated microwave radiometer retrievals suggests that weak precipitation fluxes can be nonetheless consequential to moisture budgets for supercooled clouds owing to small liquid water paths. Finally, we suggest that these ground-based precipitation rate statistics offer valuable guidance for improving the representation of polar cloud processes in large-scale models.

Polar clouds↗
NASA NTRSDate not supplied

Parallelized Carbon Oxidation and Surface Recession Model in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that often consist of a carbon-based reinforcement and a polymeric matrix. During Earth re-entry, they endure high-temperature oxidation and surface recession. Oxidation is an important mechanism for ablation, sometimes leading to the weakening, spallation, or failure of the oxidized fibers at the surface and in the char layer. However, more details are required including accurate material properties of the fiber microstructure, whether the fibers recede homogenously or localized at pits, and the role of pyrolysis outgassing in order to evaluate the role oxidation plays in the degradation and failure mechanisms of these materials. In this work, we demonstrate a parallelizable oxidation-driven ablation model developed for detailed, large-scale simulations in the DSMC code SPARTA. We also develop robust algorithms which handle the conservation of the surface state after an ablation step. Finally, we verify our model for both simple and more complex chemistry as well as microstructures with reaction statistics, oxidation depth calculations, and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗
NASA NTRS2019

The Diurnal Cycle of Precipitation and Organized Convection in a Set of Global Mesoscale Simulations with the NASA GEOS AGCM

The NASA Global Modeling and Assimilation Office (GMAO) has conducted a series of 40-day nonhydrostatic global simulations with horizontal grid spacing ranging from 200 km to 3 km, as part of the DYAMOND model intercomparison project. The Goddard Earth Observing System (GEOS) model was run with the Grell-Freitas scale-aware convection scheme, which smoothly reduces parameterized deep convection with increasing resolution. Here we evaluate the diurnal cycle and other statistics of precipitation and organized convection as a function of resolution. For validation we use the 0.1 degree IMERG precipitation and 4 km Merged IR brightness temperature datasets, focusing on four regions: the continental United States, Amazonia, the equatorial Indian ocean, and the Maritime Continent. Early results indicate good phase agreement but excessive magnitude of the continental diurnal cycle of precipitation at coarser resolutions, with improved magnitude as resolution increases and the role of parameterization is reduced. Convective cloud clusters are identified with a brightness temperature threshold, and we find realistic numbers of the largest clusters (>10^4 km sq) at all resolutions, while the number of smaller clusters increases with resolution, approaching observations when dx=3 km. The observed diurnal cycle in the cluster size distribution is also reproduced, with realistic magnitude in the highest resolution runs. Precipitation characteristics across cluster sizes are also examined. The results show the potential for global mesoscale simulations as a community science resource, and we invite collaboration to explore these runs in greater detail.

Arnold, Nathan↗
NASA NTRS2022

Combining Machine Learning and Numerical Simulation for High-Resolution PM2.5 Concentration Forecast

Forecasting ambient PM2.5 concentrations with spatiotemporal coverage is key to alerting decision-makers of pollution episodes and preventing detrimental public exposure, especially in regions with limited ground air monitoring stations. The existing methods either rely on chemical transport models (CTMs) to forecast spatial distribution of PM2.5 with nontrivial uncertainty or statistical algorithms to forecast PM2.5 concentration time-series at air monitoring locations without continuous spatial coverage. In this study, we developed a PM2.5 forecast framework by combining the robust Random Forest algorithm with a publicly accessible global CTM forecast product – NASA’s Goddard Earth Observing System “Composition Forecasting” (GEOS-CF), providing spatiotemporally continuous PM2.5 concentration forecasts for the next five days at a 1-km spatial resolution. Our forecast experiment was conducted for a region in Central China including the populous and polluted Fenwei Plain. The forecast for the next two days had overall validation R2 of 0.76 and 0.64, respectively; the R2 was around 0.5 for the following three forecast days. Spatial cross-validation showed similar validation metrics. Our forecast model, with validation normalized mean bias close to zero, substantially reduced the large biases in GEOS-CF. The proposed framework requires minimal computational resources compared to running CTMs at urban scales, enabling near-real-time PM2.5 forecast in resource-restricted environments.

PM2.5↗
NASA NTRSDate not supplied

An Investigation Into Transecting Satellites in Future Space Traffic Management Scenarios

The number of satellites expected to populate the near-Earth space environment is set to dramatically increase in the coming decade as new large constellations are approved and deployed. Current strategies for deploying new batches of these satellites often involve launching into an initial orbit, and then performing apogee raising maneuvers to reach a target altitude. Similarly, end-of-life planning for these constellation satellites can consist of de-orbit burns that lower perigee to permit disposal via re-entry. Both the raising and de-orbiting maneuvers can result in the individual satellites traveling in transecting orbits that have the potential to cross other spacecraft trajectories. While individual large constellations may be able to coexist in separate altitude and inclination bands, having thousands of satellites moving between these bands as new satellites are added and old satellites are removed could pose additional collision risks. Similar concerns have been raised regarding the impact that large numbers of university-class CubeSats might have in terms of their overall collision risk, especially as these satellites typically do not have propulsion systems for active maneuvering. To assess the impact that transecting satellites might have on future space traffic management strategies, this study explored a variety of realistic future scenarios using a high-fidelity simulation tool. The model can simulate the orbit of tens of thousands of resident space objects (RSOs) simultaneously, to include active satellites, debris, rocket bodies, or even future hypothetical satellite constellations, using a realistic force model that incorporates non-spherical gravity, atmospheric drag, and solar radiation pressure, as well as station-keeping. The simulation can be customized to accommodate different methods of calculating the probability of collision, as well as the process for determining probability ellipsoids and screening volumes. This makes it possible to replicate, and compare, different processes used by different spacecraft operators and space situational awareness (SSA) providers. As the model is run forward in time, each conjunction event is recorded, allowing for the analysis of statistics and meta-data related to these events, providing insight into the nature and frequency of potential collisions, such as whether are they active or passive objects, what size are the two objects, and who owns the objects (if known). This information makes it possible to characterize how changes to the status quo affect the number and type of conjunctions that occur, as well as the distributional effects on various types of satellite operators. To assess the general risks that transecting satellites might pose for hypothetical future space object environments, approximately 60,000 new large constellation satellites were considered, in addition to the existing catalog of approximately 7800 known resident space objects (RSOs), over a simulation period of one year. The results indicate that the future space environment will introduce a non-linear increase in conjunction events as the number of RSOs also increase. This will require adjustments to spacecraft fuel budgets in order to conduct the avoidance maneuvers necessary to minimize collision risk, both for existing and new satellites. This increase is due in large part to the higher density of RSOs and the overlap between some constellation orbits. Current catalog objects were shown to require three times more ∆V for collision avoidance (CA) maneuvers in the simulated future environment, and some constellation spacecraft were estimated to devote the majority of their annual ∆V to CA. The impact of small satellites was found to be proportional for the current space environment, and actually decreased in terms of percentage for the future scenario, suggesting that small satellites do not pose an outsized collision risk. Lastly, transecting satellites were found to contribute thousands of additional conjunctions outside of their operational orbit, and may require up to an additional 5% in CA maneuver fuel allocation.

conjunction assessment↗
NASA NTRS2015

Influence of Spanwise Boundary Conditions on Slat Noise Simulations

The slat noise from the 30P/30N high-lift system is being investigated through computational fluid dynamics simulations with the OVERFLOW code in conjunction with a Ffowcs Williams-Hawkings acoustics solver. In the present study, two different spanwise grids are being used to investigate the effect of the spanwise extent and periodicity on the near-field unsteady structures and radiated noise. The baseline grid with periodic boundary conditions has a short span equal to 1/9th of the stowed chord, whereas the other, longer span grid adds stretched grids on both sides of the core, baseline grid to allow inviscid surface boundary conditions at both ends. The results indicate that the near-field mean statistics obtained using the two grids are similar to each other, as are the directivity and spectral shapes of the radiated noise. However, periodicity forces all acoustic waves with less than one wavelength across the span to be two-dimensional, without any variation in the span. The spanwise coherence of the acoustic waves is what is needed to make estimates of the noise that would be radiated from realistic span lengths. Simulations with periodic conditions need spans of at least six slat chords to allow spanwise variation in the low-frequencies associated with the peak of broadband slat noise. Even then, the full influence of the periodicity is unclear, so employing grids with a fine, central region and highly stretched meshes that go to slip walls may be a more efficient means of capturing the spanwise decorrelation of low-frequency acoustic phenomena.

Lockard, David P.↗
NASA NTRS2013

Statistical Issues in Galaxy Cluster Cosmology

The number and growth of massive galaxy clusters are sensitive probes of cosmological structure formation. Surveys at various wavelengths can detect clusters to high redshift, but the fact that cluster mass is not directly observable complicates matters, requiring us to simultaneously constrain scaling relations of observable signals with mass. The problem can be cast as one of regression, in which the data set is truncated, the (cosmology-dependent) underlying population must be modeled, and strong, complex correlations between measurements often exist. Simulations of cosmological structure formation provide a robust prediction for the number of clusters in the Universe as a function of mass and redshift (the mass function), but they cannot reliably predict the observables used to detect clusters in sky surveys (e.g. X-ray luminosity). Consequently, observers must constrain observable-mass scaling relations using additional data, and use the scaling relation model in conjunction with the mass function to predict the number of clusters as a function of redshift and luminosity.

Galaxy↗
NASA NTRS1962

Application of Statistical Filter Theory to the Optimal Estimation of Position and Velocity on Board a Circumlunar Vehicle

Concepts from statistical filter theory are applied to the problem of in-flight determination of the position and velocity of a space vehicle for the purposes of midcourse guidance. The source of information is assumed to be an arbitrary sequence of measurements of any desired set of "observables" (e.g., space angles), the measurements being corrupted by additive errors so that the position and velocity are never known perfectly. A state transition approach is employed which leads naturally to a computational scheme that is readily implemented by a digital computer. The scheme can be regarded as a dynamical time-varying filter which weights the incoming observations in an optimal sense for use in producing an up-to-date optimal estimate of position and velocity. The advantages of the scheme are (1) it provides the best possible estimate (minimum error) based upon ensemble statistics of injection conditions and measurement errors; (2) it is extremely versatile, not requiring adherence to a predetermined observation schedule or reference trajectory; and (3) the computations are sufficiently simple to be practical in an on-board computer. A digital computer simulation of the proposed system is employed to demonstrate the feasibility of an all on-board system and to illustrate the performance attainable in a hypothetical physical situation.

SPACE VEHICLE↗
NASA NTRS1979

Ku-band radar threshold analysis

The statistics of the CFAR threshold for the Ku-band radar was determined. Exact analytical results were developed for both the mean and standard deviations in the designated search mode. The mean value is compared to the results of a previously reported simulation. The analytical results are more optimistic than the simulation results, for which no explanation is offered. The normalized standard deviation is shown to be very sensitive to signal-to-noise ratio and very insensitive to the noise correlation present in the range gates of the designated search mode. The substantial variation in the CFAR threshold is dominant at large values of SNR where the normalized standard deviation is greater than 0.3. Whether or not this significantly affects the resulting probability of detection is a matter which deserves additional attention.

Weber, C. L.↗
NASA NTRS1983

Motion/visual cueing requirements for vortex encounters during simulated transport visual approach and landing

This paper addresses the issues of motion/visual cueing fidelity requirements for vortex encounters during simulated transport visual approaches and landings. Four simulator configurations were utilized to provide objective performance measures during simulated vortex penetrations, and subjective comments from pilots were collected. The configurations used were as follows: fixed base with visual degradation (delay), fixed base with no visual degradation, moving base with visual degradation (delay), and moving base with no visual degradation. The statistical comparisons of the objective measures and the subjective pilot opinions indicated that although both minimum visual delay and motion cueing are recommended for the vortex penetration task, the visual-scene delay characteristics were not as significant a fidelity factor as was the presence of motion cues. However, this indication was applicable to a restricted task, and to transport aircraft. Although they were statistically significant, the effects of visual delay and motion cueing on the touchdown-related measures were considered to be of no practical consequence.

Parrish, R. V.↗
NASA NTRS1982

Compressibility effects in turbulence modeling

Numerical turbulence modeling is discussed with attention given to fluid property variations caused by compressibility in an adiabatic flow. The models are considered in terms of integral quantities expressed by ordinary differential equations and by those formulated as partial differential equations. Compressibility corrections for both integral and partial differential methods are reviewed. Eddy-viscosity models are explored for their capability to characterize the mass-weighted Reynolds stress, which can be accounted for with primitive and/or mass-weighted variables. Compressible flow simulations are currently constrained to low Re and zero mean dilation. The effects of compressibility are defined in wave number space by resolving the Fourier transforms of the velocity vectors into components which are perpendicular and parallel to the wave number vector. Statistical correlations then permit obtaining a value for each contribution.

Rubesin, M. W.↗
NASA NTRS1983

Center of Excellence for Applied Mathematical and Statistical Research in support of development of multicrop production monitoring capability

Efforts in support of the development of multicrop production monitoring capability are reported. In particular, segment level proportion estimation techniques based upon a mixture model were investigated. Efforts have dealt primarily with evaluation of current techniques and development of alternative ones. A comparison of techniques is provided on both simulated and LANDSAT data along with an analysis of the quality of profile variables obtained from LANDSAT data.

Woodward, W. A.↗
NASA NTRS1985

Scaling laws for testing of high lift airfoils under heavy rainfall

The results of studies regarding the effect of rainfall about aircraft are briefly reviewed. It is found that performance penalties on airfoils have been identified in subscale tests. For this reason, it is of great importance that scaling laws be dveloped to aid in the extrapolation of these data to fullscale. The present investigation represents an attempt to develop scaling laws for testing subscale airfoils under heavy rain conditions. Attention is given to rain statistics, airfoil operation in heavy rain, scaling laws, thermodynamics of condensation and/or evaporation, rainfall and airfoil scaling, aspects of splash back, film thickness, rivulets, and flap slot blockage. It is concluded that the extrapolation of airfoil performance data taken at subscale under simulated heavy rain conditions to fullscale must be undertaken with caution.

Bilanin, A. J.↗
NASA NTRS1991

Ordinary chondrites - Multivariate statistical analysis of trace element contents

The contents of mobile trace elements (Co, Au, Sb, Ga, Se, Rb, Cs, Te, Bi, Ag, In, Tl, Zn, and Cd) in Antarctic and non-Antarctic populations of H4-6 and L4-6 chondrites, were compared using standard multivariate discriminant functions borrowed from linear discriminant analysis and logistic regression. A nonstandard randomization-simulation method was developed, making it possible to carry out probability assignments on a distribution-free basis. Compositional differences were found both between the Antarctic and non-Antarctic H4-6 chondrite populations and between two L4-6 chondrite populations. It is shown that, for various types of meteorites (in particular, for the H4-6 chondrites), the Antarctic/non-Antarctic compositional difference is due to preterrestrial differences in the genesis of their parent materials.

Lipschutz, Michael E.↗