Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “moderator”

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 577 records · Page 32

Effect of Nose Length, Fuselage Length, and Nose Fineness Ratio on the Longitudinal Aerodynamic Characteristics of Two Complete Models at High Subsonic Speeds

An investigation has been made of the effects of nose length, fuselage length, and nose fineness ratio on the static longitudinal aerodynamic characteristics of an airplane model with a swept wing and low tail and of a second model with a highly tapered wing of moderate sweep and a T-tail. The tests were conducted in the Langley high-speed 7- by 10-foot tunnel at Mach numbers from 0.60 to 0.92. The nose and body cross sections were circular. For either the model with the swept wing and low tail or the model with the highly tapered wing of moderate sweep and the T-tail, the effects of forebody changes amounted primarily to rotations of the pitching-moment curves (changes in static margin) over the test ranges of angle of attack and Mach number. For the range of body shapes investigated the longitudinal stability at low lift is decreased by an increase in nose length or in fuselage length or by a reduction in nose fineness ratio when the fuselage length is held constant. In general, the stability for all model configurations showed substantially the same variation with changes in forebody area moment. The forebody changes did not alter the angle of attack at which an unstable break occurred in the moment contribution of the T-tail but did alter somewhat the magnitude of the instability.

Goodson, Kenneth W.↗

Optical Detection of Lightning from Space

Optical sensors have been developed to detect lightning from space during both day and night. These sensors have been fielded in two existing satellite missions and may be included on a third mission in 2002. Satellite-hosted, optically-based lightning detection offers three unique capabilities: (1) the ability to reliably detect lightning over large, often remote, spatial regions, (2) the ability to sample all (IC and CG) lightning, and (3) the ability to detect lightning with uniform (i.e., not range-dependent) sensitivity or detection efficiency. These represent significant departures from conventional RF-based detection techniques, which typically have strong range dependencies (biases) or range limitations in their detection capabilities. The atmospheric electricity team of the NASA Marshall Space Flight Center's Global Hydrology and Climate Center has implemented a three-step satellite lightning research program which includes three phases: proof-of-concept/climatology, science algorithm development, and operational application. The first instrument in the program, the Optical Transient Detector (OTD), is deployed on a low-earth orbit (LEO) satellite with near-polar inclination, yielding global coverage. The sensor has a 1300 x 1300 sq km field of view (FOV), moderate detection efficiency, moderate localization accuracy, and little data bias. The OTD is a proof-of-concept instrument and its mission is primarily a global lightning climatology. The limited spatial accuracy of this instrument makes it suboptimal for use in case studies, although significant science knowledge has been gained from the instrument as deployed.

Boccippio, Dennis J.↗

Influence of Subpixel Scale Cloud Top Structure on Reflectances from Overcast Stratiform Cloud Layers

Recent observational studies have shown that satellite retrievals of cloud optical depth based on plane-parallel model theory suffer from systematic biases that depend on viewing geometry, even when observations are restricted to overcast marine stratus layers, arguably the closest to plane parallel in nature. At moderate to low sun elevations, the plane-parallel model significantly overestimates the reflectance dependence on view angle in the forward-scattering direction but shows a similar dependence in the backscattering direction. Theoretical simulations are performed that show that the likely cause for this discrepancy is because the plane-parallel model assumption does not account for subpixel, scale variations in cloud-top height (i.e., "cloud bumps"). Monte Carlo simulation, comparing ID model radiances to radiances from overcast cloud field with 1) cloud-top height variation, but constant cloud volume extinction; 2) flat tops but horizontal variations in cloud volume extinction; and 3) variations in both cloud top height and cloud extinction are performed over a approximately equal to 4 km x 4 km domain (roughly the size of an individual GAC AVHRR pixel). The comparisons show that when cloud-top height variations are included, departures from 1D theory are remarkably similar (qualitatively) to those obtained observationally. In contrast, when clouds are assumed flat and only cloud extinction is variable, reflectance differences are much smaller and do not show any view-angle dependence. When both cloud-top height and cloud extinction variations are included, however, large increases in cloud extinction variability can enhance reflectance difference. The reason 3D-1D reflectance differences are more sensitive to cloud-top height variations in the forward-scattering direction (at moderate to low, sun elevations) is because photons leaving the cloud field in that direction experience fewer scattering events (low-order scattering) and are restricted to the topmost portions of the cloud. While reflectance deviations from 1D theory are much larger for bumpy clouds than for flat clouds with variable cloud extinction, differences in cloud albedo are comparable for these two cases.

Loeb, N. G.↗

Numerical Simulation of Three-Dimensional Boattail Afterbody Flowfields

The thin shear-layer approximations of the three-dimensional, compressible Navier-Stokes equations are solved for subsonic, transonic, and supersonic now over axisymmetric boattail bodies at moderate angles of attack. The plume is simulated by a solid body configuration identical to those used In experimental tests. An implicit algorithm of second-order accuracy is used to solve the equations on the ILLIAC 4 computer. The turbulence is expressed by an algebraic model applicable to three-dimensional flowfields with moderate separation. The formulation used is attractive in its independence of boundary-layer parameters. Such a simple model, however, is incapable of supporting detailed quantitative descriptions of complex shear flows. Never-the-less, good qualitative comparisons are found with three different sets of experimental date. Quantitative improvement will depend on improved turbulence transport descriptions.

Deiwert, G. S.↗

Heat Acclimation and Water-Immersion Deconditioning: Responses to Exercise

Simulated subgravity conditions, such as bed rest and water immersion, cause a decrease in a acceleration tolerance (3, 4), tilt tolerance (3, 9, 10), work capacity (5, 7), and plasma volume (1, 8-10). Moderate exercise training performed during bed rest (4) and prior to water immersion (5) provides some protection against the adverse effects of deconditioning, but the relationship between exercise and changes due to deconditioning remains unclear. Heat acclimation increases plasma and interstitial volumes, total body water, stroke volume (11), and tilt tolerance (6) and may, therefore, be a more efficient method of ameliorating deconditioning than physical training alone. The present study was undertaken to determine the effects of heat acclimation and moderate physical training, performed in cool conditions, on water-immersion deconditioning.

Shvartz, E.↗

TMI Rain Rate Estimation Over Land and Ocean Utilizing Convective and Stratiform Discrimination

Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct local minima in a regional map containing a Mesoscale Convective System (MCS). This is because of relatively small footprint size (approximately 5.5 km) and strong extinction properties in this channel of the TMI. A map of rain rate for that region, deduced from simultaneous measurements made by the Precipitation Radar (PR) on board the TRMM satellite, reveals that these T85 minima, produced by scattering, correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we infer from TMI data the presence of three different kinds of thunderstorms or Cbs. They are young, mature, and decaying Cbs that have a scale of about 20 km on the average. Two parameters enable us to infer these three kinds of Cbs objectively: a) the magnitude of scattering depression deduced from local T85 minima and b) the mean horizontal gradient of T85 around such minima. Knowing the category of a given Cb, we can estimate the rain rate associated with it. Such estimation is done with the help of relationships linking T85 minimum to rain rate in each Cb type. Similarly, a weak background rain rate in all the areas where T85 is less than 260 K is deduced with another relationship linking T85 to rain rate. In our rain retrieval model, this background rain constitutes the stratiform rain where the Cbs are absent. Initially, these relationships are optimized or tuned utilizing the PR and TMI data of a few MCS events. After such tuning, the model is applied to independent MCS cases. The areal distribution of light (1-10 mm/hr), moderate (10-20 mm/hr), and intense (> 20 mm/hr) rain rates are retrieved satisfactorally. Accuracy in the estimates of the light, moderate and intense rain areas and the mean rain rates associated with such areas in these independent MCS cases is on the average about 15%. Taking advantage of this ability of our retrieval method, one could derive the latent heat input into the atmosphere over the 760 km wide swath of the TMI radiometer in the tropics.

Prabhakara, C.↗

A Priori Analyses of Three Subgrid-Scale Models for One-Parameter Families of Filters

The decay of isotropic turbulence a compressible flow is examined by direct numerical simulation (DNS). A priori analyses of the DNS data are then performed to evaluate three subgrid-scale (SGS) models for large-eddy simulation (LES): a generalized Smagorinsky model (M1), a stress-similarity model (M2), and a gradient model (M3). The models exploit one-parameter second- or fourth-order filters of Pade type, which permit the cutoff wavenumber k(sub c) to be tuned independently of the grid increment (delta)x. The modeled (M) and exact (E) SGS-stresses are compared component-wise by correlation coefficients of the form C(E,M) computed over the entire three-dimensional fields. In general, M1 correlates poorly against exact stresses (C < 0.2), M3 correlates moderately well (C approx. 0.6), and M2 correlates remarkably well (0.8 < C < 1.0). Specifically, correlations C(E, M2) are high provided the grid and test filters are of the same order. Moreover, the highest correlations (C approx.= 1.0) result whenever the grid and test filters are identical (in both order and cutoff). Finally, present results reveal the exact SGS stresses obtained by grid filters of differing orders to be only moderately well correlated. Thus, in LES the model should not be specified independently of the filter.

Pruett, C. David↗

Influence of Transport on Two-Dimensional Model Simulation. Tracer Sensitivity to 2-D Model Transport: Long Lived Tracers - 1

In this study, we examine the sensitivity of long lived tracers to changes in the base transport components in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric mean age, and increased the rate of removal of material from the stratosphere. Increasing the stratospheric K(sub yy) increased the mean age due to the greater recycling of air parcels through the middle atmosphere, via the residual circulation, before returning to the troposphere. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and produces significantly younger ages throughout the stratosphere. Simulations with very small tropical stratospheric K(sub yy) decreased the globally averaged age of air by as much as 25% in the middle and upper stratosphere, and resulted in substantially weaker vertical age gradients above 20 km in the extratropics. We found only very small stratospheric tracer sensitivity to the magnitude of the horizontal mixing across the tropopause, and to the strength of the mesospheric gravity wave drag and diffusion used in the model. We also investigated the transport influence on chemically active tracers and found a strong age-tracer correlation, both in concentration and calculated lifetimes. The base model transport gives the most favorable overall comparison with a variety of inert tracer observations, and provides a significant improvement over our previous 1995 model transport. Moderate changes to the base transport were found to provide modest agreement with some of the measurements. Transport scenarios with residence times ranging from moderately shorter to slightly longer relative to the base case simulated N2O lifetimes that were within the observational estimates of Volk et al. [1997]. However, only scenarios with rather fast transport rates were comparable with the Volk et al. estimates of CFCl3 lifetimes. This is inconsistent with model-measurement comparisons of mean age in which the base model or slightly slower transport rates compared the most favorably with balloon SF6 data. For all comparisons shown, large transport changes away from the base case resulted in simulations that were outside the range of measurements, and in many cases, far outside this range.

Fleming, Eric L.↗

Punctured Parallel and Serial Concatenated Convolutional Codes for BPSK/QPSK Channels

As available bandwidth for communication applications becomes scarce, bandwidth-efficient modulation and coding schemes become ever important. Since their discovery in 1993, turbo codes (parallel concatenated convolutional codes) have been the center of the attention in the coding community because of their bit error rate performance near the Shannon limit. Serial concatenated convolutional codes have also been shown to be as powerful as turbo codes. In this dissertation, we introduce algorithms for designing bandwidth-efficient rate r = k/(k + 1),k = 2, 3,..., 16, parallel and rate 3/4, 7/8, and 15/16 serial concatenated convolutional codes via puncturing for BPSK/QPSK (Binary Phase Shift Keying/Quadrature Phase Shift Keying) channels. Both parallel and serial concatenated convolutional codes have initially, steep bit error rate versus signal-to-noise ratio slope (called the -"cliff region"). However, this steep slope changes to a moderate slope with increasing signal-to-noise ratio, where the slope is characterized by the weight spectrum of the code. The region after the cliff region is called the "error rate floor" which dominates the behavior of these codes in moderate to high signal-to-noise ratios. Our goal is to design high rate parallel and serial concatenated convolutional codes while minimizing the error rate floor effect. The design algorithm includes an interleaver enhancement procedure and finds the polynomial sets (only for parallel concatenated convolutional codes) and the puncturing schemes that achieve the lowest bit error rate performance around the floor for the code rates of interest.

Acikel, Omer Fatih↗

Supersonic Aftbody Closure Wind-Tunnel Testing, Data Analysis, and Computational Results

This paper reports on the model, test, and results from the Langley Supersonic Aftbody Closure wind tunnel test. This project is an experimental evaluation of the 1.5% Technology Concept Aircraft (TCA) aftbody closure model (Model 23) in the Langley Unitary Plan Wind Tunnel. The baseline TCA design is the result of a multidisciplinary, multipoint optimization process and was developed using linear design and analysis methods, supplemented with Euler and Navier-Stokes numerical methods. After a thorough design review, it was decided to use an upswept blade attached to the forebody as the mounting system. Structural concerns dictated that a wingtip support system would not be feasible. Only the aftbody part of the model is metric. The metric break was chosen to be at the fuselage station where prior aft-sting supported models had been truncated. Model 23 is thus a modified version of Model 20. The wing strongback, flap parts, and nacelles from Model 20 were used, whereas new aftbodies, a common forebody, and some new tails were fabricated. In summary, significant differences in longitudinal and direction stability and control characteristics between the ABF and ABB aftbody geometries were measured. Correcting the experimental data obtained for the TCA configuration with the flared aftbody to the representative of the baseline TCA closed aftbody will result in a significant reduction in longitudinal stability, a moderate reduction in stabilizer effectiveness and directional stability, and a moderate to significant reduction in rudder effectiveness. These reductions in the stability and control effectiveness levels of the baseline TCA closed aftbody are attributed to the reduction in carry-over area.

Allen, Jerry↗

Galileo Science Update

Live footage shows Jane Platt, JPL Public Information Office, introducing the moderator of the panel discussion. The moderator introduces the panel members include Bill O'Neil, Project Manager Galileo Primary Mission, Dr. Torrence V. Johnson Galileo Project Scientist, Prof. Ronald Greeley from Arizona State University Galileo Imaging Team, Bob Mitchell Project Manager Galileo Europa Mission, and Dr. Karen Buxbaum Galileo Science Planning Manager. The panelists give the audience information about the Galileo Mission and answers questions from the audience and from Kennedy Space Center. An animation of the Galileo Spacecraft approaching and passing Europa is presented. The panelists mentions High Resolution Images, Detail Gravity studies, Spectral Maps of non-ice materials, Jupiter studies, Callisto studies, Europa studies, and Io studies.

Source record↗

A TRMM Microwave Radiometer Rain Rate Estimation Method with Convective and Stratiform Discrimination

Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct local minima (T85min) in a regional map containing a Mesoscale Convective System (MCS). This is because of relatively small footprint size (approximately 5.5 km) and strong extinction properties in this channel of the TMI. A map of surface rain rate for that region, deduced from simultaneous measurements made by the Precipitation Radar (PR) on board the TRMM satellite, reveals that these T85(sub min), produced by scattering, correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we infer empirically from TMI data the presence of three different kinds of thunderstorms or Cbs. These Cbs are classified as young, mature, and decaying types, and are assumed to have a scale of about 20 km on the average. Two parameters are used to classify these three kinds of Cbs based on the T85 data: a) the magnitude of scattering depression deduced from local T85(sub min) and b) the mean horizontal gradient of T85 around such minima. Knowing the category of a given Cb, we can estimate the rain rate associated with it. Such estimation is done with the help of relationships linking T85min to rain rate in each Cb type. Similarly, a weak background rain rate in all the areas where T85 is less than 260 K is deduced with another relationship linking T85 to rain rate. In our rain retrieval model, this background rain constitutes stratiform rain where the Cbs are absent. Initially, these relationships are optimized or tuned utilizing the PR and TMI data of a few MCS events. After such tuning, the model is applied to independent MCS cases. The areal distribution of light (1-10 mm/hr), moderate (10-20 mm/hr), and intense (>= 20 mm/hr) rain rates are retrieved satisfactorily. Accuracy in the estimates of the light, moderate, and intense rain areas and the mean rain rates associated with such areas in these independent MCS cases is on the average about 15 %. Taking advantage of this ability of our retrieval method, one could derive the latent heat input into the atmosphere over the 760 km wide swath of the TMI radiometer in the tropics.

Prabhakara, C.↗

Energy-Dependent Ionization States of Shock-Accelerated Particles in the Solar Corona

We examine the range of possible energy dependence of the ionization states of ions that are shock-accelerated from the ambient plasma of the solar corona. If acceleration begins in a region of moderate density, sufficiently low in the corona, ions above about 0.1 MeV/amu approach an equilibrium charge state that depends primarily upon their speed and only weakly on the plasma temperature. We suggest that the large variations of the charge states with energy for ions such as Si and Fe observed in the 1997 November 6 event are consistent with stripping in moderately dense coronal. plasma during shock acceleration. In the large solar-particle events studied previously, acceleration occurs sufficiently high in the corona that even Fe ions up to 600 MeV/amu are not stripped of electrons.

Reames, Donald V.↗

Fuel Preheat Effects on Soot-Field Structure in Laminar Gas Jet Diffusion Flames Burning in 0-g and 1-g

An experimental investigation conducted at the 2.2-s drop tower of the NASA Lewis Research Center is presented to quantify the influence of moderate fuel preheat on soot-field structure within 0-g laminar gas jet diffusion flames. Parallel work in 1-g is also presented to delineate the effect of elevated fuel temperatures on soot-field structure in buoyant flames. The experimental methodology implements jet diffusion flames of nitrogen-diluted acetylene fuel burning in quiescent air at atmospheric pressure. Fuel preheat of approximately 100 K in the 0-g laminar jet diffusion flames is found to reduce soot loadings in the annular region, but causes an increase in soot volume fractions at the centerline. In addition, fuel preheat reduces the radial extent of the soot field in 0-g. In 1-g, the same fuel preheat levels have a more moderated influence on soot loadings in the annular region, but are also seen to enhance soot concentrations near the axis low in the flame. The increased soot loadings near the flame centerline, as caused by fuel preheat, are consistent with the hypothesis that preheat levels of approximately 100 K enhance fuel pyrolysis rates. The results show that the growth stage of particles transported along the soot annulus is shortened both in 1-g and 0-g when elevated fuel temperatures are used.

Konsur, Bogdan↗

Crew Factors in Flight Operations: A Survey of Fatigue Factors in Regional Airline Operations - 11

This report is the eleventh in a series on the physiological effects of flight operations on flight crews. A 119-question survey was completed by 1,424 flight crewmembers from 26 regional carriers to identify factors contributing to fatigue in regional airline operations. Eighty-nine percent of crewmembers identified fatigue as a moderate or serious concern with 88% reporting that it was a common occurrence and 92% reporting that, when it occurs, fatigue represents a moderate or serious safety issue. However, 86% reported they received no company training addressing fatigue issues. Identified fatigue factors included multiple flight segments, scheduling considerations, varying regulations, and others. The two most commonly cited fatigue factors regarded flying multiple (more than four) segments. Scheduling factors accounted for nine of the ten most common recommendations to reduce fatigue in regional operations. Differing requirements among regulations were cited as contributing to fatigue. Other identified factors were the flight deck environment, automation, and diet. The data suggested specific recommendations, including education of industry personnel about fatigue issues and examination of scheduling practices. Education plays a critical role in any effort to address fatigue. Analyzing scheduling practices and identifying potential improvements may result in reduced fatigue as well as other benefits to operations.

Co, Elizabeth L.↗

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

Comparison of Sunphotometric Measurements During the Fall 1997 ARM Intensive Observation Period

In the Fall of 1997 the Atmospheric Radiation Measurement (ARM) program held an intensive observation period (IOP) to study atmospheric aerosols using in situ and remote sensing techniques at its Southern Great Plains (SGP) site near Lamont, Oklahoma. As part of this experiment five automated, tracking sunphotometers were present to measure total column aerosol optical depth over the three-week period. which included many clear days or parts of days that were clear. The World Meteorological Organization (WMO 1993) has recommended a comparison of tracking sunphotometers to assess the ability of different instruments to arrive at similar aerosol optical depths. It was further recommended that the comparison be staged at a clean mountain site. In fact, this comparison has not occurred, but the comparison that we describe in this paper is representative of what contemporary instruments may accomplish in an environment more typical of sites where aerosols measurements will be required. The measurements were made over the period 15 September to 5 October 1997. The aerosol loading varied from extremely clean to moderately turbid conditions. In the next section the instruments will be described along with a brief explanation of the calibration techniques. The third section contains the results compared graphically on moderately turbid and fairly clean days and in a table representing the whole period. The paper ends with a section of discussion and a summary of the results.

Michalsky, J. J.↗

Geology of the Smythii and Marginis Region of the Moon: Using Integrated Remotely Sensed Data

We characterized the diverse and complex geology of the eastern limb region of the Moon using a trio of remote-sensing data sets: Clementine, Lunar Prospector, and Apollo. On the basis of Clementine-derived iron and titanium maps we classify the highlands into low-iron (3-6 wt % FeO) and high-iron (6-9 wt % FeO) units. The association of the latter with basalt deposits west of Smythii basin suggests that the highland chemical variation is the result of mixing between basalt and highland lithologies. Mare Smythii and Mare Marginis soils are compositionally similar, containing moderate iron (15-18 wt % FeO) and titanium (2.5-3.5 wt % TiO2). Smythii basin, in addition to the basalt deposits, contains an older, moderate-albedo plains unit. Our investigation reveals that the dark basin plains unit has a distinct albedo, chemistry, and surface texture and formed as a result of impact-mixing between highland and mare lithologies in approximately equal proportions. Clementine iron and maturity maps show that swirls along the northern margin of Mare Marginis have the same iron composition as the surrounding nonswirl material and indicate that the swirl material is bright because of its low agglutinate content. Gravity data for the eastern limb show high, positive Bouguer gravity anomalies for areas of thin basalt cover (e.g., Smythii basin and complex craters Joliot, Lomonosov, and Neper). We deduce that the uplift of dense mantle material is the primary (and mare basaltic fill the secondary) source for generating the concentration of mass beneath large craters and basins.

Gillis, Jeffrey J.↗