SEARCH · Search NASA

Results for “exponential”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1995

OSSE observations of GX 339-4

The Oriented Scintillation Spectrometer Experiment (OSSE) on the Compton Gamma Ray Observatory (CGRO) observed the Galactic black hole candidate GX 339-4 as a target of oppurtunity in 1991 September, in response to the outburst reported by Burst and Transient Source Experiment (BATSE). We report here on energy spectra in the 50 keV-10 MeV range obtained by OSSE. The source was detected from 50 to 400 keV at a level relative to the Crab Nebula of approximately 30%. The observed spectrum was prescribed reasonably well by a power law with an exponential cutoff; a least-squares fit yielded a photon index of 0.88 +/- 0.05 and a cutoff energy of 68 +/- 2 keV. The addition of a Compton reflection component did not significantly improve the overall fit. An optically thin thermal bremsstrahlung spectrum also provides a good fit, and the thermal Comptonization model of Sunyaev & Titarchuk, while deficient in describing the data above approximately 200 keV, cannot formally be ruled out. A pure power law with reflection does not fit the observed spectrum. During a follow-up observation made in 1991 November the intensity of the source below 100 keV had dropped by more than a factor of 40, and it was no longer detected above approximately 100 keV. The energy spectrum during the November observation could be characterized by a power law with a photon index of 2.3 +/- 0.3; the spectrum was fitted equally well with the same exponentially cutoff power-law model applied to the September observation, reduced in intensity by a factor of approximately 40. During the 1991 September observation, the luminosity in the 50-400 keV band was approximately 2 x 10(exp 37) ergs/s (assuming a distance of 4 kpc), no more than a factor of 5 below the soft X-ray luminosity of GX 339-4 observed in its X-ray high state. The luminosity during the 1991 November observation was approximately 5 x 10(exp 35) ergs/s. Extrapolations of both the exponentially cutoff power-law and Sunyaev-Titarchuk models to the approximately 5-20 keV X-ray band yield flux levels very close to that observed by Ginga during an overlapping interval in 1991 September, when GX 339-4 was reported to be in its low state. This may be one of the strongest indications to date of a direct correspondence between the low X-ray state and gamma-ray outbursts of GX 339-4.

Grabelsky, D. A.↗

NASA NTRS2000

R-Function Relationships for Application in the Fractional Calculus

The F-function, and its generalization the R-function, are of fundamental importance in the fractional calculus. It has been shown that the solution of the fundamental linear fractional differential equation may be expressed in terms of these functions. These functions serve as generalizations of the exponential function in the solution of fractional differential equations. Because of this central role in the fractional calculus, this paper explores various intrarelationships of the R-function, which will be useful in further analysis. Relationships of the R-function to the common exponential function, e(t), and its fractional derivatives are shown. From the relationships developed, some important approximations are observed. Further, the inverse relationships of the exponential function, el, in terms of the R-function are developed. Also, some approximations for the R-function are developed.

Lorenzo, Carl F.↗

NASA NTRS2001

Our Sun V: A Bright Young Sun Consistent with Helioseismology and Warm Temperatures on Ancient Earth and Mars

The relatively warm temperatures required on early Earth and Mars have been difficult to account for with warming from greenhouse gases. A slightly more massive young Sun would be brighter than predicted by the standard solar model, simultaneously resolving this problem for both Earth and Mars. We computed high-precision solar models with seven initial masses, from Mi = 1.01 to 1.07 solar mass - the latter being the maximum permitted if the early Earth is not to lose its water via a moist greenhouse effect. The relatively modest early mass loss that is required remains consistent with observational limits on mass loss from young stars and with estimates of the past solar wind obtained from lunar rocks. We considered three types of mass loss rates: (1) a reasonable choice of a simple exponential decline, (2) an extreme step-function case that gives the maximum effect consistent with observations, and (3) the radical case of a linear decline which is inconsistent with the solar wind mass loss estimates from lunar rocks. Our computations demonstrated that mass loss leaves a fingerprint oil the Sun's internal structure large enough to be detectable with helioseismic observations. All of our mass-losing solar models were consistent with the helioseismic observations; in fact, our preferred mass-losing cases were in marginally better agreement with the helioseismology than the standard solar model was, although this difference was smaller than the effects of other uncertainties in the input physics and in the solar composition. Mass loss has only a relatively minor effect on the predicted lithium depletion; the major portion of the solar lithium depletion must still be due to rotational mixing. Thus the modest mass loss cases considered here cannot be ruled out by observed lithium depletions. For the three mass loss types considered, the preferred initial masses were 1.07 solar mass for the exponential case and 1.04 solar mass for the step-function and linear cases; all of these provided high enough solar fluxes at Mars 3.8 Gyr ago to be consistent with the existence of liquid water. For a more massive early Sun, the planets would have had to be closer to the young Sun in order to end up in their present orbits; the orbital radii of the planets would vary inversely with the solar mass. Both of these effects contribute to the fact that the early solar flux at the planets would have been considerably higher than that of the standard solar model at that time. In fact, the 1.07 solar mass exponential case has a flux at birth 5% higher than the present solar flux, while the radical 1.04 solar mass linear case has a nearly constant flux over the first 3 Gyr only about 10% lower than at present. The early solar evolution would be in the opposite direction in the H-R diagram to that of the standard Sun.

Sackmann, I.-Juliana↗

NASA NTRS2002

A Parametric Study of Fine-scale Turbulence Mixing Noise

The present paper is a study of aerodynamic noise spectra from model functions that describe the source. The study is motivated by the need to improve the spectral shape of the MGBK jet noise prediction methodology at high frequency. The predicted spectral shape usually appears less broadband than measurements and faster decaying at high frequency. Theoretical representation of the source is based on Lilley's equation. Numerical simulations of high-speed subsonic jets as well as some recent turbulence measurements reveal a number of interesting statistical properties of turbulence correlation functions that may have a bearing on radiated noise. These studies indicate that an exponential spatial function may be a more appropriate representation of a two-point correlation compared to its Gaussian counterpart. The effect of source non-compactness on spectral shape is discussed. It is shown that source non-compactness could well be the differentiating factor between the Gaussian and exponential model functions. In particular, the fall-off of the noise spectra at high frequency is studied and it is shown that a non-compact source with an exponential model function results in a broader spectrum and better agreement with data. An alternate source model that represents the source as a covariance of the convective derivative of fine-scale turbulence kinetic energy is also examined.

Khavaran, Abbas↗

NASA NTRS2002

A Self-Replicating Ligase Ribozyme

A self-replicating molecule directs the covalent assembly of component molecules to form a product that is of identical composition to the parent. When the newly formed product also is able to direct the assembly of product molecules, the self-replicating system can be termed autocatalytic. A self-replicating system was developed based on a ribozyme that catalyzes the assembly of additional copies of Itself through an RNA-catalyzed RNA ligation reaction. The R3C ligase ribozyme was redesigned so that it would ligate two substrates to generate an exact copy of itself, which then would behave in a similar manner. This self-replicating system depends on the catalytic nature of the RNA for the generation of copies. A linear dependence was observed between the initial rate of formation of new copies and the starting concentration of ribozyme, consistent with exponential growth. The autocatalytic rate constant was 0.011 per min, whereas the initial rate of reaction in the absence of pre-existing ribozyme was only 3.3 x 10(exp -11) M per min. Exponential growth was limited, however, because newly formed ribozyme molecules had greater difficulty forming a productive complex with the two substrates. Further optimization of the system may lead to the sustained exponential growth of ribozymes that undergo self-replication.

Paul, Natasha↗

NASA NTRS2004

Changes in the X-Ray Emission from the Magnetar Candidate 1E 2259+586 During its 2002 Outburst

An outburst of more than 80 individual bursts, similar to those seen from Soft Gamma Repeaters (SGRs), was detected from the anomalous X-ray pulsar (AXP) 1E 2259+586 in 2002 June. Coincident with this burst activity were gross changes in the pulsed flux, persistent flux, energy spectrum, pulse profile, and spin-down of the underlying X-ray source. We present Rossi X-Ray Timing Explorer and X-Ray Multi-Mirror Mission observations of 1E 2259+586 that show the evolution of the aforementioned source parameters during and following this episode and identify recovery timescales for each. Specifically, we observe an X-ray flux increase (pulsed and phase-averaged) by more than an order of magnitude having two distinct components. The first component is linked to the burst activity and decays within approx. 2 days, during which the energy spectrum is considerably harder than during the quiescent state of the source. The second component decays over the year following the glitch according to a power law in time with an exponent -0.22 +/- 0.01. The pulsed fraction decreased initially to approx. 15% rms but recovered rapidly to the preoutburst level of approx. 23% within the first 3 days. The pulse profile changed significantly during the outburst and recovered almost fully within 2 months of the outburst. A glitch of size Delta(sib (nu)max) = (4.24 +/- 0.11) x 10(exp -6) was observed in 1E 2259+586, which preceded the observed burst activity. The glitch could not be well fitted with a simple partial exponential recovery. An exponential rise of approx. 20% of the frequency jump with a timescale of approx. 14 days results in a significantly better fit to the data; however, contamination from a systematic drift in the phase of the pulse profile cannot be excluded. A fraction of the glitch (approx. 19%) was recovered in a quasi-exponential manner having a recovery timescale of approx. 16 days. The long-term postglitch spin-down rate decreased in magnitude relative to the preglitch value. The changes in the source properties of 1E 2259+586 during its 2002 outburst are shown to be qualitatively similar to changes seen during or following burst activity in two SGRs, thus further solidifying the common nature of SGRs and AXP's as magnetars. The changes in persistent emission properties of 1E 2259+586 suggest that the star underwent a plastic deformation of the crust that simultaneously impacted the superfluid interior (crustal and possibly core superfluid) and the magnetosphere. Finally, the changes in persistent emission properties coincident with burst activity in 1E 2259+586 enabled us to infer previous burst-active episodes from this and other AXP's. The nondetection of these outbursts by all-sky gamma-ray instruments suggests that the number of active magnetar candidates in our Galaxy is larger than previously thought.

Woods, P. M.↗

NASA NTRS2010

Spacecraft Solar Particle Event (SPE) Shielding: Shielding Effectiveness as a Function of SPE model as Determined with the FLUKA Radiation Transport Code

Analysis of both satellite and surface neutron monitor data demonstrate that the widely utilized Exponential model of solar particle event (SPE) proton kinetic energy spectra can seriously underestimate SPE proton flux, especially at the highest kinetic energies. The more recently developed Band model produces better agreement with neutron monitor data ground level events (GLEs) and is believed to be considerably more accurate at high kinetic energies. Here, we report the results of modeling and simulation studies in which the radiation transport code FLUKA (FLUktuierende KAskade) is used to determine the changes in total ionizing dose (TID) and single-event environments (SEE) behind aluminum, polyethylene, carbon, and titanium shielding masses when the assumed form (i. e., Band or Exponential) of the solar particle event (SPE) kinetic energy spectra is changed. FLUKA simulations have fully three dimensions with an isotropic particle flux incident on a concentric spherical shell shielding mass and detector structure. The effects are reported for both energetic primary protons penetrating the shield mass and secondary particle showers caused by energetic primary protons colliding with shielding mass nuclei. Our results, in agreement with previous studies, show that use of the Exponential form of the event

Koontz, Steve↗

NASA NTRS2010

A Study of Cirrus Ice Particle Size Distribution Using TC4 Observations

An analysis of two days of in situ observations of ice particle size spectra, in convectively generated cirrus, obtained during NASA s Tropical Composition, Cloud, and Climate Coupling (TC4) mission is presented. The observed spectra are examined for their fit to the exponential, gamma, and lognormal function distributions. Characteristic particle size and concentration density scales are determined using two (for the exponential) or three (for the gamma and lognormal functions) moments of the spectra. It is shown that transformed exponential, gamma, and lognormal distributions should collapse onto standard curves. An examination of the transformed spectra, and of deviations of the transformed spectra from the standard curves, shows that the lognormal function provides a better fit to the observed spectra.

Tian, Lin↗

NASA NTRS2005

Improved Compression of Wavelet-Transformed Images

A recently developed data-compression method is an adaptive technique for coding quantized wavelet-transformed data, nominally as part of a complete image-data compressor. Unlike some other approaches, this method admits a simple implementation and does not rely on the use of large code tables. A common data compression approach, particularly for images, is to perform a wavelet transform on the input data, and then losslessly compress a quantized version of the wavelet-transformed data. Under this compression approach, it is common for the quantized data to include long sequences, or runs, of zeros. The new coding method uses prefixfree codes for the nonnegative integers as part of an adaptive algorithm for compressing the quantized wavelet-transformed data by run-length coding. In the form of run-length coding used here, the data sequence to be encoded is parsed into strings consisting of some number (possibly 0) of zeros, followed by a nonzero value. The nonzero value and the length of the run of zeros are encoded. For a data stream that contains a sufficiently high frequency of zeros, this method is known to be more effective than using a single variable length code to encode each symbol. The specific prefix-free codes used are from two classes of variable-length codes: a class known as Golomb codes, and a class known as exponential-Golomb codes. The codes within each class are indexed by a single integer parameter. The present method uses exponential-Golomb codes for the lengths of the runs of zeros, and Golomb codes for the nonzero values. The code parameters within each code class are determined adaptively on the fly as compression proceeds, on the basis of statistics from previously encoded values. In particular, a simple adaptive method has been devised to select the parameter identifying the particular exponential-Golomb code to use. The method tracks the average number of bits used to encode recent runlengths, and takes the difference between this average length and the code parameter. When this difference falls outside a fixed range, the code parameter is updated (increased or decreased). The Golomb code parameter is selected based on the average magnitude of recently encoded nonzero samples. The coding method requires no floating- point operations, and more readily adapts to local statistics than other methods. The method can also accommodate arbitrarily large input values and arbitrarily long runs of zeros. In practice, this means that changes in the dynamic range or size of the input data set would not require a change to the compressor. The algorithm has been tested in computational experiments on test images. A comparison with a previously developed algorithm that uses large code tables (generated via Huffman coding on training data) suggests that the data-compression effectiveness of the present algorithm is comparable to the best performance achievable by the previously developed algorithm.

Kiely, Aaron↗

NASA NTRS2003

Quantum Search in Hilbert Space

A proposed quantum-computing algorithm would perform a search for an item of information in a database stored in a Hilbert-space memory structure. The algorithm is intended to make it possible to search relatively quickly through a large database under conditions in which available computing resources would otherwise be considered inadequate to perform such a task. The algorithm would apply, more specifically, to a relational database in which information would be stored in a set of N complex orthonormal vectors, each of N dimensions (where N can be exponentially large). Each vector would constitute one row of a unitary matrix, from which one would derive the Hamiltonian operator (and hence the evolutionary operator) of a quantum system. In other words, all the stored information would be mapped onto a unitary operator acting on a quantum state that would represent the item of information to be retrieved. Then one could exploit quantum parallelism: one could pose all search queries simultaneously by performing a quantum measurement on the system. In so doing, one would effectively solve the search problem in one computational step. One could exploit the direct- and inner-product decomposability of the unitary matrix to make the dimensionality of the memory space exponentially large by use of only linear resources. However, inasmuch as the necessary preprocessing (the mapping of the stored information into a Hilbert space) could be exponentially expensive, the proposed algorithm would likely be most beneficial in applications in which the resources available for preprocessing were much greater than those available for searching.

Zak, Michail↗

NASA NTRS2012

The Ground Flash Fraction Retrieval Algorithm Employing Differential Evolution: Simulations and Applications

The ability to estimate the fraction of ground flashes in a set of flashes observed by a satellite lightning imager, such as the future GOES-R Geostationary Lightning Mapper (GLM), would likely improve operational and scientific applications (e.g., severe weather warnings, lightning nitrogen oxides studies, and global electric circuit analyses). A Bayesian inversion method, called the Ground Flash Fraction Retrieval Algorithm (GoFFRA), was recently developed for estimating the ground flash fraction. The method uses a constrained mixed exponential distribution model to describe a particular lightning optical measurement called the Maximum Group Area (MGA). To obtain the optimum model parameters (one of which is the desired ground flash fraction), a scalar function must be minimized. This minimization is difficult because of two problems: (1) Label Switching (LS), and (2) Parameter Identity Theft (PIT). The LS problem is well known in the literature on mixed exponential distributions, and the PIT problem was discovered in this study. Each problem occurs when one allows the numerical minimizer to freely roam through the parameter search space; this allows certain solution parameters to interchange roles which leads to fundamental ambiguities, and solution error. A major accomplishment of this study is that we have employed a state-of-the-art genetic-based global optimization algorithm called Differential Evolution (DE) that constrains the parameter search in such a way as to remove both the LS and PIT problems. To test the performance of the GoFFRA when DE is employed, we applied it to analyze simulated MGA datasets that we generated from known mixed exponential distributions. Moreover, we evaluated the GoFFRA/DE method by applying it to analyze actual MGAs derived from low-Earth orbiting lightning imaging sensor data; the actual MGA data were classified as either ground or cloud flash MGAs using National Lightning Detection Network[TM] (NLDN) data. Solution error plots are provided for both the simulations and actual data analyses.

Koshak, William↗

NASA NTRS2012

High-latitude E Region Ionosphere-thermosphere Coupling: A Comparative Study Using in Situ and Incoherent Scatter Radar Observations

We present in situ and ground-based measurements of the ratio k of ion cyclotronangular frequency to ion-neutral momentum transfer collision frequency to investigateionosphere-thermosphere (IT) coupling in the auroral E region. In situ observations were obtained by NASA sounding rocket 36.234, which was launched into the nightsideE region ionosphere at 1229 UT on 19 January 2007 from Poker Flat, AK. The payload carried instrumentation to determine ion drift angle and electric field vectors. Neutral winds were measured by triangulating a chemical tracer released from rocket 41.064 launched two minutes later. k is calculated from the rotation of the ion drift angle relative to the E-cross-B drift direction in a frame co-rotating with the payload. Between the altitudes of 118 km and 130 km k increases exponentially with a scale height of 9.3 +/- 0.7 km, deviating from an exponential above 130 km. k = 1 at an altitude z(sub0) of 119.9 +/- 0.5 km. The ratio was also estimated from Poker Flat Incoherent Scatter Radar (PFISR) measurements using the rotation of ion velocity with altitude. Exponential fits to the PFISR measurements made during the flight of 41.064 yield z(sub0) 115.9 +/- 1.2 km and a scale height of 9.1 +/- 1.0 km. Differences between in situ and ground-based measurements show that the E region atmospheric densities were structured vertically and/or horizontally on scales of 1 km to 10 km. There were no signs of ionospheric structure in ion density or ion temperature below scales of 1 km. The observations demonstrate the accuracy with which the in situ and PFISR data may be used as probes of IT coupling.

Burchill, J. K.↗

NASA NTRS2012

Spectrum and Morphology of the Two Brightest Milagro Sources in the Cygnus Region: MGRO J2019+37 and MGRO J2031+41

The Cygnus region is a very bright and complex portion of the TeV sky, host to unidentified sources and a diffuse excess with respect to conventional cosmic-ray propagation models. Two of the brightest TeV sources, MGRO J2019+37 and MGRO J2031+41, are analyzed using Milagro data with a new technique, and their emission is tested under two different spectral assumptions: a power law and a power law with an exponential cutoff. The new analysis technique is based on an energy estimator that uses the fraction of photomultiplier tubes in the observatory that detect the extensive air shower. The photon spectrum is measured in the range 1-100 TeV using the last three years of Milagro data (2005-2008), with the detector in its final configuration. An F-test indicates that MGRO J2019+37 is better fit by a power law with an exponential cutoff than by a simple power law. The best-fitting parameters for the power law with exponential cutoff model are a normalization at 10 TeV of 7(sup +5 sub −2) × 10(exp −10)/ s /sq m/ TeV, a spectral index of 2.0(sup +0.5 sub -10), and a cutoff energy of 29(sup +50 sub −16) TeV. MGRO J2031+41 shows no evidence of a cutoff. The best-fitting parameters for a power law are a normalization of 2.1(sup +0.6 sub −0.6) × 10(exp −10)/ s/sq m/ TeV and a spectral index of 3.22(sup +0.23 sub −0.18. The overall flux is subject to a approx.. 30% systematic uncertainty. The systematic uncertainty on the power-law indices is approx. 0.1. Both uncertainties have been verified with cosmic-ray data. A comparison with previous results from TeV J2032+4130, MGRO J2031+41, and MGRO J2019+37 is also presented.

BRIGHTEST↗

NASA NTRS2003

Rapid Acceleration of a Coronal Mass Ejection in the Low Corona and Implications of Propagation

A high-velocity Coronal Mass Ejection (CME) associated with the 2002 April 21 X1.5 flare is studied using a unique set of observations from the Transition Region and Coronal Explorer (TRACE), the Ultraviolet Coronagraph Spectrometer (UVCS), and the Large-Angle Spectrometric Coronagraph (LASCO). The event is first observed as a rapid rise in GOES X-rays, followed by simultaneous conjugate footpoint brightenings connected by an ascending loop or flux-rope feature. While expanding, the appearance of the feature remains remarkably constant as it passes through the TRACE 195 A passband and LASCO fields-of-view, allowing its height-time behavior to be accurately determined. An analytic function, having exponential and linear components, is found to represent the height-time evolution of the CME in the range 1.05-26 R. The CME acceleration rises exponentially to approx. 900 km/sq s within approximately 20-min, peaking at approx.1400 m/sq s when the leading edge is at approx. 1.7 R. The acceleration subsequently falls off as a slowly varying exponential for approx.,90-min. At distances beyond approx. 3.4 R, the height-time profile is approximately linear with a constant velocity of approx. 2400 km/s. These results are briefly discussed in light of recent kinematic models of CMEs.

Gallagher, Peter T.↗

NASA NTRSDate not supplied

Modeling the Worst-Case-Scenario of Soil Drydown for Agricultural Drought Monitoring in Kansas Using NASA Earth Observations and In-situ Parameters

Climatologists in Kansas have observed rapid shifts from wet to dry conditions and rapid intensification of drought conditions for the state in recent years. As a leading state in agricultural production, drought events can lead to decreased crop yields and impact the state’s economy. The exponential decay of soil moisture content is a major consequence of drought and can exacerbate these impacts on agricultural production. This study examined the rate of soil moisture drydown in Kansas to understand and monitor the worst-case-scenario of short-term soil water shortage. Soil Moisture Active Passive (SMAP) L-band Radiometer rootzone soil moisture imagery from March 2015 to June 2019 was used to estimate the minimum and maximum soil moisture capacity throughout the state. In situ soil moisture parameters, such as the exponential decay constant, were approximated from the Kansas Mesonet system. These variables were inputs into an exponential decay model that forecasted gridded rootzone soil moisture in a scenario that assumed no precipitation for a user-identified period of time. The forecasts were output on the order of weeks to capture the rapid shifts in soil moisture content. Another model was used to identify areas in Kansas currently below a given percentage of the relative soil moisture saturation, measured as a function of minimum, maximum, and current soil moisture. This model forecasted the number of days until each cell in the soil moisture grid reached the input percentage of relative soil moisture saturation. The gridded model outputs and workflow were developed in partnership with the Kansas Water Office and Kansas Office of the State Climatologist at Kansas State University, as a part of their drought monitoring efforts.

NASA DEVELOP↗

NASA NTRSDate not supplied

Quantifying Leak Rates and Vacuum Impact on Granular Soil Column Pneumatic Seal

Introduction: All In-Situ Resource Utilization (ISRU) regolith reactor systems that aim to extract gases or vapors must maintain a seal along the flow path of the regolith from the inlet hopper to where extracted gases or vapors are generated. Using a column of regolith at the inlet has been proposed as an alternative to systems using series of mechanical valves to reduce mass and complexity. Preliminary investigations into the use of regolith itself as a barrier to escaping fluids have been conducted, which this research seeks to expand upon. Different vertical regolith simulant plug diameters and heights were tested under vacuum with dry nitrogen to observe leak rates and approximate steady state pressure gradient equilibria. The effect of absolute ambient pressure was also studied by repeating tests at incremental pressures ranging from 10e-5 Torr up to 1 atm to corroborate ambient bench tests from previous work. The impact of regolith dynamics was addressed to emulate continuous-feed regolith reaction processes by measuring leak rates through continuously draining simulant columns. OB-1 lunar highlands regolith simulant was used. Experimental Setup: Regolith simulant columns were constructed using a motorized butterfly valve separating upper acrylic tubes of various diameters from a larger, lower aluminum pipe plenum (See Figure 1). The butterfly valve was cracked open such that gas could freely pass through whereas regolith would bridge and stay suspended as a cohesive plug above in the acrylic portion. These assemblies were loaded into a vacuum chamber, and a compressed dry nitrogen line was routed into the lower plenum. Pressure transducers measured the absolute pressure of the lower plenum, which was compared to absolute chamber pressure to monitor the pressure gradient across the regolith simulant. A precision mass flow controller enabled fine control of plenum pressure. The plenum pressure was quasi-statically increased to set points informed by previous work [1] before shutting off the gas supply to allow the now pressurized gas volume to dissipate through the simulant into the chamber in what was referred to as a ‘leak test.’ In a separate set of tests, the pressure was allowed to quasi-statically build until visual disturbances in the soil were observed via remote camera feeds. Initial leak rate tests were repeated on the disturbed soil to observe changes in leak rates between initial as-poured and post-bubbling soil conditions. This test schedule was applied to simulant columns of three different diameters (approximately 2, 4, and 6 inches) as well as two different heights (around 8 and 16 inches) and was repeated at least three times for each configuration. Plenums were also pressurized while the butterfly valves were opened slightly further to facilitate slow, controlled, constant flowrate draining of the simulant to observe the impact that continuous-feed systems might have on leak rate through the regolith column. Finally, the regolith columns were reset and refilled and static leak rate tests were repeated approximately every decade in ambient pressure from 10e-5 Torr up to ambient conditions to observe the impact of transitioning out of the molecular flow regime on regolith column seal performance. Analysis: Trends with respect to each factor examined were constructed and reported alongside commentary on the most likely responsible physical phenomena. Confidence intervals were applied to the predictive findings of the data by quantifying instrumentation error and accounting for probabilistic impacts of curve fitting the generic exponential function onto pressure decay data. The leak rate is assumed directly proportional to the exponential term coefficient, while the steady-state static pressure gradient is assumed to be the asymptotic limit of each curve fit. Finally, a generic model developed by Ogino et al. [2] and later tuned for ambient testing was re-tuned using vacuum data to estimate leak rates on the lunar surface. Results: Pressure decay rates decreased up to 27-fold under vacuum when compared to ambient pressure, as shown in Figure 2. Simulant column diameter was found to exponentially correlate to both regolith column leak rate and steady-state equilibrium pressure delta. Height was found to be linearly related to leak rate –doubling the regolith column height consistently yielded around a 40% reduction in leak rate. Dynamic regolith draining was found to slightly increase leak rate, although to a lower degree than expected. Conclusion: Because of the transition into molecular gas flow through regolith, vacuum testing is required to properly estimate performance. Regolith columns should be designed as narrow as possible to avoid bridging and minimize leak rate, and as tall as is feasible to both further reduce leak rate and to help prevent bubbling, which has been shown to degrade the quality of the regolith seal. Disrupted soil was found to exhibit higher leak rates than pristine, as-poured soil, indicating the need to consider the impact of recovery from depressurization events. Regolith columns appear equally feasible for both dosed and continuous-feed reactors.

ISRU↗

NASA NTRS2024

The abcd Reliability Growth Model

This paper presents a modification of the well-known Duane-Crow reliability growth model. In the abcd reliability growth model, the initial period of exponential decline of the failure rate in the Duane-Crow model may be followed by a period of constant failure rate. Data often show that an exponential decline in failures is followed by a constant failure rate. If a growth model including only the initial period of exponential decline is applied to increasingly longer failure rate data sets, the data will include longer periods of constant failure rate, and the estimated reliability growth rate will decline from an initially high value down toward zero. Using the Duane-Crow model without extending it to include a possible period of constant failure rate may create the mistaken impression that the initial reliability growth continues forever, but at an ever decreasing rate.

Reliability growth↗

NASA NTRS2024

The abcd Reliability Growth Model

Reliability growth↗