Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution secondary”

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 433 records · Page 24

Evidence of galaxy assembly bias in SDSS DR7 galaxy samples from count statistics

We present observational constraints on the galaxy–halo connection, focusing particularly on galaxy assembly bias from a novel combination of counts-in-cylinders statistics, P(N CIC ), with the standard measurements of the projected two-point correlation function w p (r p ), and number density n gal of galaxies. We measure n gal , w p (r p ), and P(N CIC ) for volume-limited, luminosity-threshold samples of galaxies selected from SDSS DR7, and use them to constrain halo occupation distribution (HOD) models, including a model in which galaxy occupation depends upon a secondary halo property, namely halo concentration. We detect significant positive central assembly bias for the M r < -20.0 and M r <-19.5 samples. Central galaxies preferentially reside within haloes of high concentration at fixed mass. Positive central assembly bias is also favoured in the M r < -20.5 and M r < -19.0 samples. We find no evidence of central assembly bias in the M r < -21.0 sample. We observe only a marginal preference for negative satellite assembly bias in the M r < -20.0 and M r < -19.0 samples, and non-zero satellite assembly bias is not indicated in other samples. Our findings underscore the necessity of accounting for galaxy assembly bias when interpreting galaxy survey data, and demonstrate the potential of count statistics in extracting information from the spatial distribution of galaxies, which could be applied to both galaxy–halo connection studies and cosmological analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

Short-range order and its impact on the CrCoNi medium-entropy alloy

Traditional metallic alloys are mixtures of elements in which the atoms of minority species tend to be distributed randomly if they are below their solubility limit, or to form secondary phases if they are above it. The concept of multiple-principal-element alloys has recently expanded this view, as these materials are single-phase solid solutions of generally equiatomic mixtures of metallic elements. This group of materials has received much interest owing to their enhanced mechanical properties 1,2,3,4,5 . They are usually called medium-entropy alloys in ternary systems and high-entropy alloys in quaternary or quinary systems, alluding to their high degree of configurational entropy. However, the question has remained as to how random these solid solutions actually are, with the influence of short-range order being suggested in computational simulations but not seen experimentally 6,7 . In this paper we report the observation, using energy-filtered transmission electron microscopy, of structural features attributable to short-range order in the CrCoNi medium-entropy alloy. Increasing amounts of such order give rise to both higher stacking-fault energy and hardness. These findings suggest that the degree of local ordering at the nanometre scale can be tailored through thermomechanical processing, providing a new avenue for tuning the mechanical properties of medium- and high-entropy alloys.

36 MATERIALS SCIENCE↗

Coannular supersonic ejector nozzles

The nozzles described exhibit a flow field which is supersonic except for the initial flow region, and the secondary mass flow is typically about five percent of the primary core flow. The features to improve the accuracy of the performance calculations are discussed. A special calculation is made to get as realistic a sonic line as possible for this geometry, using an analysis developed by Brown. The mixing between the secondary and core flows is treated to account for entrainment of the secondary flow into core. Both of these phenomena directly affect the pressure distribution on the shroud and therefore, the thrust that the nozzle produces. The importance of using a realistic sonic line and a mixing analysis is stressed.

Bishop, A. R.↗

The elemental abundances of hydrogen through nickel in the low energy cosmic rays

The relative abundances of the elements H through Ni in the galactic cosmic rays have been measured in the energy range 70-280 MeV/nucleon with the University of Chicago cosmic ray telescope on board the satellite IMP-8 from January 1973 to September 1978. Cosmic ray source abundances have been derived by extrapolating the measured composition back to the source. A key factor in the propagation calculation is the use of a pathlength distribution and a solar modulation level shown to be consistent with the secondary to primary ratios and their energy dependence below about 1 GeV/n.

Garcia-Munoz, M.↗

Experiments on supersonic turbulent flow development in a square duct

The nature of supersonic, turbulent, adiabatic-wall flow in a square duct is investigated experimentally over a development length of x/D between 0 and 20 for a uniform flow, Mach 3.9 condition at the duct inlet. Initial discussion centers on the duct configuration itself, which was designed specifically to minimize wave effects and nozzle-induced distortion in the flow. Total pressure contours and local skin friction coefficient distributions are presented which show that the flow develops in a manner similar to that observed for the incompressible case. In particular, undulations exist in total pressure contours within the cross plane and in transverse skin friction coefficient distributions, which are indicative of the presence of a well-defined secondary flow superimposed upon the primary flow. The results are analyzed to show that local law-of-the-wall behavior extends well into the corner region, which implies that wall functions conventionally applied in two-equation type turbulence models, when suitably defined for compressible flow, may also be applied to supersonic streamwise corner flows.

Gessner, F. B.↗

Definition of the Spatial Resolution of X-Ray Microanalysis in Thin Foils

The spatial resolution of X-ray microanalysis in thin foils is defined in terms of the incident electron beam diameter and the average beam broadening. The beam diameter is defined as the full width tenth maximum of a Gaussian intensity distribution. The spatial resolution is calculated by a convolution of the beam diameter and the average beam broadening. This definition of the spatial resolution can be related simply to experimental measurements of composition profiles across interphase interfaces. Monte Carlo calculations using a high-speed parallel supercomputer show good agreement with this definition of the spatial resolution and calculations based on this definition. The agreement is good over a range of specimen thicknesses and atomic number, but is poor when excessive beam tailing distorts the assumed Gaussian electron intensity distributions. Beam tailing occurs in low-Z materials because of fast secondary electrons and in high-Z materials because of plural scattering.

Williams, D. B.↗

Cold-Flow Study of Low Frequency Pressure Instability in Hybrid Rocket Motors

Past experience with hybrid rockets has shown that certain motor operating conditions are conducive to the formation of low frequency pressure oscillations, or flow instabilities, within the motor. Both past and present work in the hybrid propulsion community acknowledges deficiencies in the understanding of such behavior, though it seems probable that the answer lies in an interaction between the flow dynamics and the combustion heat release. Knowledge of the fundamental flow dynamics is essential to the basic understanding of the overall stability problem. A first step in this direction was a study conducted at NASA Marshall Space Flight Center (MSFC), centered around a laboratory-scale two dimensional water flow model of a hybrid rocket motor. Principal objectives included: (1) visualization of flow and measurement of flow velocity distributions: (2) assessment of the importance of shear layer instabilities in driving motor pressure oscillations; (3) determination of the interactions between flow induced shear layers with the mainstream flow, the secondary (wall) throughflow, and solid boundaries; (4) investigation of the interactions between wall flow oscillations and the mainstream flow pressure distribution.

Jenkins, Rhonald M.↗

Effects of Thermal and Mechanical Processing on Microstructures and Desired Properties of Particle-Strengthened Cu-Cr-Nb Alloys

Ternary Cu-Cr-Nb alloys, particularly Cu-8 Cr-4 Nb (in at.%), have demonstrated good thermal stability as well as high strength and conductivity at elevated temperatures. The initial powder material has a bimodal size distribution of Cr2Nb precipitates. Primary Cr2Nb precipitates are approx. 1 micron, and secondary Cr2Nb particles are 30-200 nm. The particle coarsening was analyzed and found to follow LSW-type behavior, This study provides a detailed examination of the stability and strengthening effects of Cr2Nb particles. This investigation also revealed that the primary particles provide direct grain boundary pinning and indirect grain boundary strengthening but virtually no Orowan strengthening. The secondary particles found within grains do provide Orowan strengthening. For extruded material, grain bound-ary strengthening (Hall-Petch effect) accounts for two-thirds of the strength with Orowan effects contributing the remainder. The proven advantages of Cu-Cr-Nb were the motivation to improve these attributes via microstructural refinement. Mechanical milling (MM) of Cu- 4 Cr-2 Nb and Cu-8 Cr-2 Nb produced an increase in hot pressed Vickers hardness of 122% and 96%, respectively. The increase in hardness was more due to Cu grain-size refinement than to Cr,,Nb refinement. This study also demonstrated enhanced stability of MM Cu-4 Cr-2 Nb. Hot pressed 4 h milled Cu-4 Cr-2 Nb experienced only a 22% drop in hardness when annealed at 1273 K for 50 h versus a 30% drop for extruded Cu-8 Cr-4 Nb. The goal of improving the strength and stability of Cu-4 Cr-2 Nb to better than such properties for as- extruded Cu-8 Cr-4 Nb has been met. In addition, a figure-of-merit (FOM) coupling hardness and thermal conductivity was maximized for the case of 4 h milled Cu-4 Cr-2 Nb material. Overall, Cu-Cr-Nb alloys not only possess high strength, conductivity and thermal stability but also can be further developed to improve strength and stability.

Anderson, Kenneth Reed↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optomechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optormechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by, the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Determination of the Three-Dimensional Morphology of ALH84001 and Biogenic MV-1 Magnetite: Comparison of Results from Electron Tomography and Classical Transmission Electron Microscopy

Dated at approximately 3.9 billion years of age, carbonate disks, found within fractures of the host rock of Martian meteorite ALH84001, have been interpreted as secondary minerals that formed at low temperature in an aqueous medium. Heterogeneously distributed within these disks are magnetite nanocrystals that are of Martian origin. Approximately one quarter of these magnetites have morphological and chemical similarities to magnetite particles produced by magnetotactic bacteria strain MV-1, which are ubiquitous in aquatic habitats on Earth. Moreover, these types of magnetite particles are not known or expected to be produced by abiotic means either through geological processes or synthetically in the laboratory. The remaining three quarters of the ALH84001 magnetites are likely products of multiple processes including, but not limited to, precipitation from a hydrothermal fluid, thermal decomposition of the carbonate matrix in which they are embedded, and extracellular formation by dissimilatory Fe-reducing bacteria. We have proposed that the origins of magnetites in ALH84001 can be best explained as the products of multiple processes, one of which is biological. Recently the three-dimensional (3-D) external morphology of the purported biogenic fraction of the ALH84001 magnetites has been the subject of considerable debate. We report here the 3-D geometry of biogenic magnetite crystals extracted from MV-1 and of those extracted from ALH84001 carbonate disks using a combination of high resolution classical and tomographic transmission electron microscopy (TEM). We focus on answering the following questions: (1) which technique provides adequate information to deduce the 3-D external crystal morphology?; and, (2) what is the precise 3-D geometry of the ALH84001 and MV-1 magnetites?

Thomas-Keprta, Kathie L.↗

Mechanisms for Diurnal Variability of Global Tropical Rainfall Observed from TRMM

The behavior and various controls of diurnal variability in tropical-subtropical rainfall are investigated using Tropical Rainfall Measuring Mission (TRMM) precipitation measurements retrieved from: (1) TRMM Microwave Imager (TMI), (2) Precipitation Radar (PR), and (3) TMI/PR Combined, standard level 2 algorithms for the 1998 annual cycle. Results show that the diurnal variability characteristics of precipitation are consistent for all three algorithms, providing assurance that TRMM retrievals are providing consistent estimates of rainfall variability. As anticipated, most ocean areas exhibit more rainfall at night, while over most land areas rainfall peaks during daytime ,however, various important exceptions are found. The dominant feature of the oceanic diurnal cycle is a rainfall maximum in late-evening/early-morning (LE-EM) hours, while over land the dominant maximum occurs in the mid- to late-afternoon (MLA). In conjunction with these maxima are pronounced seasonal variations of the diurnal amplitudes. Amplitude analysis shows that the diurnal pattern and its seasonal evolution are closely related to the rainfall accumulation pattern and its seasonal evolution. In addition, the horizontal distribution of diurnal variability indicates that for oceanic rainfall there is a secondary MLA maximum, co-existing with the LE-EM maximum, at latitudes dominated by large scale convergence and deep convection. Analogously, there is a preponderance for an LE-EM maximum over land, co-existing with the stronger MLA maximum, although it is not evident that this secondary continental feature is closely associated with the large scale circulation. The ocean results clearly indicate that rainfall diurnal variability associated with large scale convection is an integral part of the atmospheric general circulation.

Yang, Song↗

Mechanisms for Diurnal Variability of Global Tropical Rainfall Observed from TRMM

The behavior and various controls of diurnal variability in tropical-subtropical rainfall are investigated using Tropical Rainfall Measuring Mission (TRMM) precipitation measurements retrieved from the three level-2 TRMM standard profile algorithms for the 1998 annual cycle. Results show that diurnal variability characteristics of precipitation are consistent for all three algorithms, providing assurance that TRMM retrievals are producing consistent estimates of rainfall variability. As anticipated, most ocean areas exhibit more rainfall at night, while over most land areas, rainfall peaks during daytime; however, important exceptions are noted. The dominant feature of the oceanic diurnal cycle is a rainfall maximum in late-evening-early-morning (LE-EM) hours, while over land the dominant maximum occurs in the mid- to late afternoon (MLA). In conjunction with these maxima are pronounced seasonal variations of the diurnal amplitudes. Amplitude analysis shows that the diurnal pattern and its seasonal evolution are closely related to the rainfall accumulation pattern and its seasonal evolution. In addition, the horizontal distribution of diurnal variability indicates that for oceanic rainfall, there is a secondary MLA maximum coexisting with the LE-EM maximum at latitudes dominated by large-scale convergence and deep convection. Analogously, there is a preponderancy for an LE EM maximum over land coexisting with the stronger MLA maximum, although it is not evident that this secondary continental feature is closely associated with the large-scale circulation. Neither of the secondary maxima exhibit phase behavior that can be considered semidiurnal in nature. Diurnal rainfall variability over the ocean associated with large-scale convection is clearly an integral component of the general circulation. Phase analysis reveals differences in regional and seasonal features of the diurnal cycle, indicating that underlying forcing mechanisms differ from place to place. This is underscored by the appearance of secondary ocean maxima in the presence of large-scale convection, along with other important features. Among these, there are clear-cut differences between the diurnal variability of seasonal rainfall over the mid-Pacific and Indian Ocean Basins. The mid-Pacific exhibits double maxima in spring and winter but only LE-EM maxima in summer and autumn, while the Indian Ocean exhibits double maxima in spring and summer and only an LE-EM maximum in autumn and winter. There are also evident daytime maxima within the major large-scale marine stratocumulus regions off the west coasts of continents. The study concludes with a discussion concerning how the observational evidence either supports or repudiates possible forcing mechanisms that have been suggested to explain diurnal rainfall variability.

Yang, Song↗

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗

A Chandra X-Ray Survey of Ejecta in the Cassiopeia A Supernova Remnant

We present a survey of the X-ray emitting ejecta in the Cassiopeia A supernova remnant based on an extensive analysis of over 6000 spectral regions extracted on 2.5-10" angular scales using the Chandra 1 Ms observation. We interpret these results in the context of hydrodynamical models for the evolution of the remnant. The distributions of fitted temperature and ionization age are highly peaked and suggest that the ejecta were subjected to multiple secondary shocks. Based on the fitted emission measure and element abundances, and an estimate of the emitting volume, we derive masses for the X-ray emitting ejecta as well as showing the distribution of the mass of various elements over the remnant. The total shocked Fe mass appears to be roughly 0.14 Solar Mass, which accounts for nearly all of the mass expected in Fe ejecta. We find two populations of Fe ejecta, that associated with normal Si-burning and that associated with alpha-rich freeze-out, with a mass ratio of approximately 2:1. Surprisingly, essentially all of this Fe (both components) is well outside the central regions of the SNR, presumably having been ejected by hydrodynamic instabilities during the explosion. We discuss this, and its implications for the neutron star kick.

Hwang, Una↗

Transmission-and-Distribution Dynamic Co-simulation Framework for Distributed Energy Resource Frequency Response

The rapid deployment of distributed energy resources (DERs) in distribution networks has made it challenging to balance the transmission system and stabilize frequency. DERs have the ability to provide frequency regulation services; however, existing frequency dynamic simulation tools - which were developed mainly for the transmission system - lack the capability to simulate distribution network dynamics with high penetrations of DERs. Although electromagnetic transient simulation tools can simulate distribution network dynamics, the computation efficiency limits their use for large-scale transmission-and-distribution (T&D) co-simulation. This paper presents an efficient open-source T&D dynamic co-simulation framework for DER frequency response based on the HELICS platform and off-the-shelf T&D simulators. The challenge of synchronizing the simulation time between the transmission network and the DERs in the distribution network is solved through the detailed modeling of DERs in frequency dynamic models while DER power flow models are also preserved in the distribution networks, thereby respecting local voltage constraints when dispatching DER power for frequency response. DER frequency response (primary and secondary) is simulated in case studies to validate the proposed framework. Last, the accuracy of the proposed co-simulation model is benchmarked, and a large T&D system simulation (2k transmission and 1M distribution nodes) is presented to demonstrate the efficiency and effectiveness of the overall framework.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Confirmation of saturation equilibrium conditions in crater populations

We have continued work on realistic numerical models of cratered surfaces, as first reported at last year's LPSC. We confirm the saturation equilibrium level with a new, independent test. One of us has developed a realistic computer simulation of a cratered surface. The model starts with a smooth surface or fractal topography, and adds primary craters according to the cumulative power law with exponent -1.83, as observed on lunar maria and Martian plains. Each crater has an ejecta blanket with the volume of the crater, feathering out to a distance of 4 crater radii. We use the model to test the levels of saturation equilibrium reached in naturally occurring systems, by increasing crater density and observing its dependence on various parameters. In particular, we have tested to see if these artificial systems reach the level found by Hartmann on heavily cratered planetary surfaces, hypothesized to be the natural saturation equilibrium level. This year's work gives the first results of a crater population that includes secondaries. Our model 'Gaskell-4' (September, 1992) includes primaries as described above, but also includes a secondary population, defined by exponent -4. We allowed the largest secondary from each primary to be 0.10 times the size of the primary. These parameters will be changed to test their effects in future models. The model gives realistic images of a cratered surface although it appears richer in secondaries than real surfaces are. The effect of running the model toward saturation gives interesting results for the diameter distribution. Our most heavily cratered surface had the input number of primary craters reach about 0.65 times the hypothesized saturation equilibrium, but the input number rises to more than 100 times that level for secondaries below 1.4 km in size.

Hartmann, William K.↗