Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Space Weathering in the Fine Size Fractions of Lunar Soils: Soil Maturity Effects

The effects of space weathering on the optical properties of lunar materials have been well documented. These effects include a reddened continuum slope, lowered albedo, and attenuated absorption features in reflectance spectra of lunar soils as compared to finely comminuted rocks from the same Apollo sites. However, the regolith processes that cause these effects are not well known, nor is the petrographic setting of the products of these processes fully understood. A Lunar Soil Characterization Consortium has been formed with the purpose of systematically integrating chemical and mineralogical data with the optical properties of lunar soils. Understanding space-weathering effects is critical in order to fully integrate the lunar sample collection with remotely-sensed data from recent robotic missions (e.g., Lunar Prospector, Clementine, and Galileo) We have shown that depositional processes (condensation of impact-derived vapors, sputter deposits, accreted impact material, e.g., splash glass, spherules, etc.) are a major factor in the modification of the optical surfaces of lunar regolith materials. In mature soils, it is the size and distribution of the nanophase metal in the soil grains that has the major effect on optical properties. In this report, we compare and contrast the space-weathering effects in an immature and a mature soil with similar elemental compositions. For this study, we analyzed <10 micron sieve fractions of two Apollo 17 soils, 79221 (mature, Is/FeO = 81) and 71061 (immature, Is/FeO = 14). Details of the sieving procedures and allocation scheme are given else where. The results of other detailed chemical, mineralogical, and spectroscopic analyses of these soil samples are reported elsewhere. A representative sample of each soil was embedded in low-viscosity epoxy, and thin sections (about 70nm thick) were obtained through ultra microtomy. The thin sections used for these analyses typically contained cross sections of up to 500 individual grains. The thin sections were studied using a JEOL 2010 transmission electron microscope (TEM) equipped with a thin window energy-dispersive X-ray (EDX) spectrometer. An individual thin section was selected from each soil, and for each grain in the section we determined (1) the elemental composition by EDX; (2) whether the grain was crystalline or glassy using electron diffraction and darkfield imaging; (3) the presence or absence of rims and accreted material; and (4) the distribution of nanophase Fe where present. Most of the categories are self-evident; however, we divide the agglutinate derived material into agglutinitic glass (glass with approximately the same composition as the bulk soil that contains nanophase Fe with or without vesicles) and agglutinate fragments, which are composed of crystalline grains and agglutinitic glass. Lithic fragments are defined as polymineralic grains with no glass. Pyroxene grains have been divided into high- and low-Ca groups. As expected, there are a number of differences in the petrography of the <10-microns fractions of 79221 and 71061 given the great difference in their respective maturities, but we focus here on two major distinctions: agglutinate content and the number of grains with micropatina. Slightly over 50% of the particles in 79221 consist of agglutinitic glass and agglutinate fragments, while the remainder are predominantly crystalline mineral grains. The agglutinic glass particles contain abundant nanophase Fe and vesicles. Angular particles are rare, with most showing smooth, rounded exteriors, Of the mineral grains analyzed thus far, over 90% of the grains have amorphous rims that contain nanophase Fe (these rims are believed to have formed by vapor deposition and irradiation effects). The nanophase Fe in these rims probably accounts for a significant fraction of the increase in Is/FeO measured in these size fractions. In addition to the rims, the majority of particles also show abundant accreted material in the form of glass splashes and spherules that also contain nanophase Fe. In stark contrast, the surfaces of the mineral grains in the 71061 sample are relatively prisitine, as only about 14% of the mineral grains in the sample exhibited amorphous rims. Furthermore, the mineral particles are more angular and show greater surface roughness than in the mature sample. Accreted material on particle surfaces is rare. Agglutinitic material is a major component of the 71061 sample; however, nanophase Fe and vesicles are not as well developed as in the 79221 sample. It is now recognized that nanophase Fe is probably the main agent in modifying the optical properties of lunar soil grains. The most important result of this study is the observation that in the fine size fractions of mature soils, nearly every grain has nanophase Fe within 100 run of the particle surface. (Additional Information contained in original)

Keller, L. P.↗
NASA NTRS2002

Ion Storage Tests with the High Performance Antimatter Trap (HiPAT)

The NASA/Marshall Space Flight Centers (NASA/MSFC) Propulsion Research Center (PRC) is evaluating an antiproton storage system, referred to as the High Performance Antiproton Trap (HiPAT). This interest stems from the sheer energy represented by matter/antimatter annihilation process with has an energy density approximately 10 order of magnitude above that of chemical propellants. In other terms, one gram of antiprotons contains the equivalent energy of approximately 23 space shuttle external tanks or ET's (each ET contains roughly 740,000 kgs of fuel and oxidizer). This incredible source of stored energy, if harnessed, would be an enabling technology for deep space mission where both spacecraft weight and propulsion performance are key to satisfying aggressive mission requirements. The HiPAT hardware consists of a 4 Tesla superconductor system, an ultra high vacuum test section (vacuum approaching 10(exp -12) torr), and a high voltage confinement electrode system (up to 20 kvolts operation). The current laboratory layout is illustrated. The HiPAT designed objectives included storage of up to 1 trillion antiprotons with corresponding lifetimes approaching 18 days. To date, testing has centered on the storage of positive hydrogen ions produced in situ by a stream of high-energy electrons that passes through the trapping region. However, due to space charge issues and electron beam compression as it passes through the HiPAT central field, current ion production is limited to less then 50,000 ions. Ion lifetime was determined by counting particle populations at the end of various storage time intervals. Particle detection was accomplished by destructively expelling the ions against a micro-channel plate located just outside the traps magnetic field. The effect of radio frequency (RF) stabilization on the lifetime of trapped particles was also examined. This technique, referred to as a rotating wall, made use of a segmented electrode located near the center of the trap on which various phases of a particular frequency were applied. Various experiments were performed illustrating the ability of an RF drive to both prolong and reduced the lifetimes of various ion species depending on the selected frequency. HiPAT is now being reconfigured for testing with an ion source that will provide both positive and negative hydrogen ions from an external source. This ion system shall provide higher fill capacity (order of million of ions per shot), stacking of multiple shots, and injection schemes typical of a realistic antiproton delivery system.

Martin, James J.↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
DOE OSTI2020

Petrophysical and geomechanical properties of gas hydrate-bearing sediments recovered from Alaska North Slope 2018 Hydrate-01 Stratigraphic Test Well

Knowledge of petrophysical and geomechanical properties of gas hydrate-bearing sediments are essential for predicting reservoir response due to gas production from gas hydrate reservoirs. That information will be critical parameters for designing production well completion such as specification of depressurization pump, water storage tank, and mesh size of sand screen. In December 2018, Stratigraphic Test Well Hydrate-01 was drilled in the western part of the Prudhoe Bay Unit, Alaska North Slope as part of the technical planning effort for a future long-term production test by collaborative team of DOE/NETL, USGS, and MH21-S (Boswell et al., 2020, Collett et al., 2020, Okinaka et al., 2020). Data set of logging-while-drilling (LWD) were acquired (Haines et al., 2020, Suzuki et al., 2019) and core sampling depth was determined on-site.Side-wall pressure coring was conducted to recover gas hydrate-bearing sediments from two reservoir sections named Unit-B and Unit-D. A total of 34 cores were successfully recovered by 5 runs of a wire-line deployed pressure corer (CoreVault® System - Halliburton). Core analysis scheme of this project are shown in Figure 1. All cores were quenched in liquid nitrogen while at high pressure in the laboratory of Stratum Reservoir, LLC. at Anchorage (Figure 1, a)). And the cores were removed from the pressure corer autoclave with temperature support of dry ice and stored under liquid nitrogen at atmospheric pressure. 19 damaged cores were processed to index property measurements directly such as grain size, grain density. 4 of another 17 cores were depressurized and trimmed for making a plug to analyze petrophysical properties of host sediments. Unsteady-state permeability test was conducted to obtain relative water permeability to gas and core scale NMR T2 distribution measurement was performed for evaluating pore size distribution at Houston (Figure 1, b)). Remained high quality 13 cores were preserved with gas hydrate for advanced laboratory analysis. National Institute of Advanced Industrial Science and Technology, as a part of the Japanese National Hydrate Research Program (MH21-S, funded by Ministry of Economy, Trade and Industry), received the samples at Sapporo, Japan for advanced core analysis. High-resolution X-ray CT was used to analyze the quality of the samples, which showed undisturbed lithological layers. Cores were lathed into cylindrical shape and distributed for multi property measurements (Figure 1, c)).<p>As a result, sediment from Unit-D is categorized as silty sand at ~37% porosity with ~80% gas hydrate saturation. An average hydration number n = 6.16 was measured by Raman spectroscopy. An average intrinsic permeability of ~400 mD and in situ effective permeability (with hydrate) on the order of ~10 mD. The Unit B recovered cores consisted of well sorted sand at ~40% porosity with ~95% gas hydrate saturation. An average intrinsic permeability of ~1 Darcy and in situ effective permeability on the order of ~30 mD was measured for the Unit B cores. There was a small permeability reduction due to porosity loss with increasing effective stress that simulated consolidation behavior along with depressurization in the highly permeable sandy sediment. The apparent minimum change in porosity and permeability may be caused by the low compressibility of quartz sand grains in the recovered cores. XRD and thermal conductivity analysis also suggested high quartz content. Triaxial compression tests established internal friction angles based on the Mohr-Coulomb's failure criterion, which are 40° for hydrate-bearing sediment and 29.8° for hydrate free sediment.</p>

Yoneda, Jun↗
NASA NTRS2003

Satellite Ocean-Color Validation Using Ships of Opportunity

The investigation s main objective is to collect from platforms of opportunity (merchant ships, research vessels) concomitant normalized water-leaving radiance and aerosol optical thickness data over the world s oceans. A global, long-term data set of these variables is needed to verify whether satellite retrievals of normalized water-leaving radiance are within acceptable error limits and, eventually, to adjust atmospheric correction schemes. To achieve this objective, volunteer officers, technicians, and scientists onboard the selected ships collect data from portable SIMBAD and Advanced SIMBAD (SIMBADA) radiometers. These instruments are specifically designed for evaluation of satellite-derived ocean color. They measure radiance in spectral bands typical of ocean-color sensors. The SIMBAD version measures in 5 spectral bands centered at 443, 490, 560, 670, and 870 nm, and the Advanced SIMBAD version in 11 spectral bands centered at 350, 380, 412, 443, 490, 510, 565, 620, 670, 750, and 870 nm. Aerosol optical thickness is obtained by viewing the sun disk like a classic sun photometer. Normalized water-leaving radiance, or marine reflectance, is obtained by viewing the ocean surface through a vertical polarizer in a specific geometry (nadir angle of 45o and relative azimuth angle of 135deg) to minimize direct sun glint and reflected sky radiation. The SIMBAD and SIMBADA data, after proper quality control and processing, are delivered to the SIMBIOS project office for inclusion in the SeaBASS archive. They complement data collected in a similar way by the Laboratoire d'Optique Atmospherique of the University of Lille, France. The SIMBAD and SIMBADA data are used to check the radiometric calibration of satellite ocean-color sensors after launch and to evaluate derived ocean-color variables (i.e., normalized water-leaving radiance, aerosol optical thickness, and aerosol type). Analysis of the SIMBAD and SIMBADA data provides information on the accuracy of satellite retrievals of normalized water-leaving radiance, an understanding of the discrepancies between satellite and in situ data, and algorithms that reduce the discrepancies, contributing to more accurate and consistent global ocean color data sets.

Frouin, Robert↗
NASA NTRS2016

Winds Measured by the Rover Environmental Monitoring Station (REMS) During the Mars Science Laboratory (MSL) Rover's Bagnold Dunes Campaign and Comparison with Numerical Modeling Using MarsWRF

A high density of REMS wind measurements were collected in three science investigations during MSL's Bagnold Dunes Campaign, which took place over approx. 80 sols around southern winter solstice (Ls approx. 90deg) and constituted the first in situ analysis of the environmental conditions, morphology, structure, and composition of an active dune field on Mars. The Wind Characterization Investigation was designed to fully characterize the near-surface wind field just outside the dunes and confirmed the primarily upslope/downslope flow expected from theory and modeling of the circulation on the slopes of Aeolis Mons in this season. The basic pattern of winds is 'upslope' (from the northwest, heading up Aeolis Mons) during the daytime (approx. 09:00-17:00 or 18:00) and 'downslope' (from the southeast, heading down Aeolis Mons) at night (approx. 20:00 to some time before 08:00). Between these times the wind rotates largely clockwise, giving generally westerly winds mid-morning and easterly winds in the early evening. The timings of these direction changes are relatively consistent from sol to sol; however, the wind direction and speed at any given time shows considerable intersol variability. This pattern and timing is similar to predictions from the MarsWRF numerical model, run at a resolution of approx. 490 m in this region, although the model predicts the upslope winds to have a stronger component from the E than the W, misses a wind speed peak at approx. 09:00, and under-predicts the strength of daytime wind speeds by approx. 2-4 m/s. The Namib Dune Lee Investigation reveals 'blocking' of northerly winds by the dune, leaving primarily a westerly component to the daytime winds, and also shows a broadening of the 1 Hz wind speed distribution likely associated with lee turbulence. The Namib Dune Side Investigation measured primarily daytime winds at the side of the same dune, in support of aeolian change detection experiments designed to put limits on the saltation threshold, and also appears to show the influence of the dune body on the local flow, though less clearly than in the lee. Using a vertical grid with lower resolution near the surface reduces the relative strength of nighttime winds predicted by MarsWRF and produces a peak in wind speed at approx. 09:00, improving the match to the observed diurnal variation of wind speed, albeit with an offset in magnitude. The annual wind field predicted using this grid also provides a far better match to observations of aeolian dune morphology and motion in the Bagnold Dunes. However, the lower overall wind speeds than observed and disagreement with the observed wind direction at approx. 09:00 suggest that the problem has not been solved and that alternative boundary layer mixing schemes should be explored which may result in more mixing of momentum down to the near-surface from higher layers. These results demonstrate a strong need for in situ wind data to constrain the setup and assumptions used in numerical models, so that they may be used with more confidence to predict the circulation at other times and locations on Mars.

Newman, Claire E.↗
DOE OSTI2025

High-power test of a C-band linear accelerating structure with an RFSoC-based LLRF system

Normal conducting linear particle accelerators consist of multiple rf stations with accelerating structure cavities. Low-level rf (LLRF) systems are employed to set the phase and amplitude of the field in the accelerating structure and to compensate for the pulse-to-pulse fluctuation of the rf field in the accelerating structures with a feedback loop. The LLRF systems are typically implemented with analog rf mixers, heterodyne-based architectures, and discrete data converters. There are multiple rf signals from each of the rf stations, so the number of rf channels required increases rapidly with multiple rf stations. With a large number of rf channels, the footprint, component cost, and system complexity of the LLRF hardware will increase significantly. To meet the design goals of being compact and affordable for future accelerators, we have designed the next-generation LLRF (NG-LLRF) with a higher integration level based on RFSoC technology. The NG-LLRF system samples rf signals directly and performs rf mixing digitally. Further, the NG-LLRF has been characterized in loopback mode to evaluate the performance of the system and has also been tested with a standing-wave accelerating structure, a prototype for the Cool Copper Collider (C 3 ) with a peak rf power level up to 16.45 MW. The loopback test demonstrated amplitude fluctuation below 0.15% and phase fluctuation below 0.15°, which are considerably better than the requirements of C 3 . The rf signals from the different stages of the accelerating structure at different power levels are measured by the NG-LLRF, which will be critical references for the control algorithm designs. The NG-LLRF also offers flexibility in waveform modulation, so we have used rf pulses with various modulation schemes, which could be useful for controlling some of the rf stations in accelerators. In this paper, the high-power test results at different stages of the test setup will be summarized, analyzed, and discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

HFTS-2 Completions Design and State-of-the-Art Diagnostics Results

Various completion schemes were tested as part of the Hydraulic Fracture Test Site-2 (HFTS-2) project. The wells were instrumented with permanent fiber optic (FO) cable, and unique near-wellbore and far-field region surveys were acquired to evaluate spatial and temporal hydraulic fracture characteristics. The state-of-the-art diagnostics included substantial monitoring of both the well stimulation and longer-term production performance. The objective was to evaluate and improve stimulation distribution effectiveness (SDE) for different completion designs based on the diagnostics results. A hydraulic fracturing (HF) profile (i.e., proppant/fluid allocation among clusters) was obtained using near-wellbore Distributed Acoustic Sensing (DAS) data acquired during stimulation. Statistical analysis provided a quantitative assessment of SDE for the different completion designs. A novel workflow was developed to assess time-lapse fracture property changes at the cluster and stage levels based on Distributed Strain Sensing (DSS) data acquired during the production phase. These time-lapse fracture property changes were compared with the HF profile. Based on the statistical analysis, Rate Per Cluster (RPC) during proppant slurry placement appears to be a primary completion design variable, rather than the typical Rate Per Perf (RPP). If RPC was not maintained at a certain level, the likelihood of uneven distribution increased dramatically. The stage configurations that created low RPC while maintaining high RPP did not show improved SDE. The strain change results based on DSS data during the production phase indicated near-wellbore fracture property changes (e.g., aperture changes during well shut-in and flow periods) and aligned well with near-wellbore DAS signals acquired during stimulation. More importantly, fracture property changes based on DSS data acquired during two tests over a 7-month period show high consistency and good alignment with the near-wellbore HF profile. This workflow quantitatively evaluates SDE for different completion designs using diagnostic results and could help optimize completion strategy for future development projects. Here, the analysis and interpretation of results demonstrate the technical feasibility of using permanent fiber optic cable for stimulation and production monitoring and indic

58 GEOSCIENCES↗
DOE OSTI2020

Considerations for using Privacy Preserving Machine Learning Techniques for Safeguards

In international nuclear safeguards, the International Atomic Energy Agency (IAEA) is tasked with inspecting and verifying nuclear facilities and their activities. Data analytics and machine learning to support inspections require large amounts of data that nuclear facility operators may consider proprietary or sensitive, so the IAEA may not have full access. Allowing computation over private data without compromising its security therefore has value for safeguards inspections and analysis. Privacy-preserving machine learning (PPML) consists of security-focused techniques that allow data analytics and machine learning algorithms to run on sensitive data without revealing it. This includes ideas like homomorphic encryption (HE), secure multiparty computation (SMPC), and secure enclaves. HE allows algorithms and mathematical operations to be conducted directly on the encrypted data instead of first decrypting it. With SMPC, multiple entities collaboratively compute over distributed data such that no party is able to directly view any others’ original data. Secure enclaves allow computation to take place in a separate and heavily blocked-off section of a CPU. Techniques like these allow for several potential use cases in which the security of data is essential. With SMPC, machine learning models can be trained over the input data from multiple entities, resulting in a model that all users can benefit from without leaking the input data from any particular entity. With SMPC or a zero-knowledge proof (ZKP), an algorithm returning some single answer or truth value can be run on someone else’s data without ever needing to see that data, potentially allowing for verification or proof of some underlying question. HE can allow for outsourcing computation on data to a hostile or untrusted environment. Although most of the research in this field resides within the health and financial domains, tools from PPML may have similar applications in nuclear safeguards. Allowing the IAEA to compute over proprietary information, such as process models and raw sensor data using PPML techniques, provides the baseline for running complex analytics without needing direct unencrypted access to the underlying data, maintaining its privacy. Important limitations to consider for these techniques include the efficiency and level of security required. The security of HE and SMPC come at the cost of speed—the significant amount of overhead means that algorithms implemented in these protocols and encryption schemes are slower than when run on plaintext. Additionally, several important parameters determine what techniques or protocols are used based on the security requirements. SMPC protocols may need to be selected for resistance against a party that attempts to deviate from the protocol to distort the result or gain access to additional information, and a protocol secure against these attacks may further increase the overhead of the algorithm.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

An Adaptable Power System with Software Control Algorithm

A low cost, flexible and modular spacecraft power system design was developed in response to a call for an architecture that could accommodate multiple missions in the small to medium load range. Three upcoming satellites will use this design, with one launch date in 1999 and two in the year 2000. The design consists of modular hardware that can be scaled up or down, without additional cost, to suit missions in the 200 to 600 Watt orbital average load range. The design will be applied to satellite orbits that are circular, polar elliptical and a libration point orbit. Mission unique adaptations are accomplished in software and firmware. In designing this advanced, adaptable power system, the major goals were reduction in weight volume and cost. This power system design represents reductions in weight of 78 percent, volume of 86 percent and cost of 65 percent from previous comparable systems. The efforts to miniaturize the electronics without sacrificing performance has created streamlined power electronics with control functions residing in the system microprocessor. The power system design can handle any battery size up to 50 Amp-hour and any battery technology. The three current implementations will use both nickel cadmium and nickel hydrogen batteries ranging in size from 21 to 50 Amp-hours. Multiple batteries can be used by adding another battery module. Any solar cell technology can be used and various array layouts can be incorporated with no change in Power System Electronics (PSE) hardware. Other features of the design are the standardized interfaces between cards and subsystems and immunity to radiation effects up to 30 krad Total Ionizing Dose (TID) and 35 Mev/cm(exp 2)-kg for Single Event Effects (SEE). The control algorithm for the power system resides in a radiation-hardened microprocessor. A table driven software design allows for flexibility in mission specific requirements. By storing critical power system constants in memory, modifying the system code for other programs is simple. These constants can be altered also by ground command, or in response to an anomolous event. All critical power system functions have backup hardware functions to prevent a software or computer glitch from propagating. A number of battery charge control schemes can be implemented by selecting the proper control terms in the code. The architecture allows the design engineer to tune the system response to various system components and anticipated load profiles without costly alterations. A design trade was made with the size, weight and power dissipation of the electronics versus the performance of the power bus to load variations. Linear, fine control is maintained with a streamlined electronics design. This paper describes the hardware design as well as the software control algorithm. The challenges of closing the system control loop digitally is discussed. Control loop margin and power system performance is presented. Lab measurements are shown and compared to the system response of a hardware model running actual flight software.

Castell, Karen↗
NASA NTRS2020

A Global End-Member Approach to Derive aCDOM(440) from Near-Surface Optical Measurements

This study establishes an optical inversion scheme for deriving the absorption coefficient of colored (or chromophoric, depending on the literature) dissolved organic material (CDOM) at the 440 nm wavelength, which can be applied to global water masses with near-equal efficacy. The approach uses a ratio of diffuse attenuation coefficient spectral end members, i.e., a short and long wavelength pair. The global perspective is established by sampling "extremely" clear water plus a generalized extent in turbidity and optical properties that each span three decades of dynamic range. A unique data set was collected in oceanic, coastal, and inland waters (as shallow as 0.6 m) from the North Pacific Ocean, the Arctic Ocean, Hawaii, Japan, Puerto Rico, and the east and west coasts of the United States. The data were partitioned using subjective categorizations to define a validation quality subset of conservative water masses, i.e., the inflow and outflow of properties constrain the range in the gradient of a constituent, plus 15 subcategories of water masses that were not evolving conservatively. The dependence on subcategories was confirmed with an objective methodology based on cluster analysis techniques. The latter defined five distinct classes with validation quality data present in all classes, but which also decreased in percent composition as a function of increasing class number and optical complexity. Four different algorithms based on different validation quality end members were validated with accuracies of 1.–6.2 %, wherein pairs with the largest spectral span were most accurate. Although algorithm accuracy decreased with the inclusion of more subcategories containing non-conservative water masses, changes to the algorithm fit were small when a preponderance of subcategories were included. The high accuracy for all end-member algorithms was the result of data acquisition and data processing improvements, e.g., increased vertical sampling resolution to less than 1mm and a boundary constraint to mitigate wave focusing effects, respectively. An independent evaluation with a historical database confirmed the consistency of the algorithmic approach and its application to quality assurance, e.g., to flag data outside expected ranges, identify suspect spectra, and objectively determine the in-water extrapolation interval by converging agreement for all applicable end-member algorithms. The legacy data exhibit degraded performance (as 44 % uncertainty) due to a lack of high-quality near-surface observations, especially for clear waters wherein wave-focusing effects are problematic. The novel optical approach allows the in situ estimation of an in-water constituent in keeping with the accuracy obtained in the laboratory.

Stanford B Hooker↗
NASA NTRS2018

Acceleration Modes and Transitions in Pulsed Plasma Accelerators

Pulsed plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas, ionizing and accelerating it through the Lorentz body force. Two plasma accelerator types employing this general scheme have typically been studied: the gas-fed pulsed plasma thruster and the quasi-steady magnetoplasmadynamic (MPD) accelerator. The gas-fed pulsed plasma accelerator is generally represented as a completely transient device discharging in approximately 1-10 microseconds. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j (current density) by B (magnetic field) body force, entraining propellant it encounters. This process is sometimes referred to as detonation-mode acceleration because the current sheet representation approximates that of a strong shock propagating through the gas. Acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new current sheet at the breech and depriving the initial sheet of additional acceleration. In the quasi-steady MPD accelerator, the pulse is lengthened to approximately 1 millisecond or longer and maintained at an approximately constant level during discharge. The time over which the transient phenomena experienced during startup typically occur is short relative to the overall discharge time, which is now long enough for the plasma to assume a relatively steady-state configuration. The ionized gas flows through a stationary current channel in a manner that is sometimes referred to as the deflagration-mode of operation. The plasma experiences electromagnetic acceleration as it flows through the current channel towards the exit of the device. A device that had a short pulse length but appeared to operate in a plasma acceleration regime different from the gas-fed pulsed plasma accelerators was developed by Cheng, et al. The Coaxial High ENerGy (CHENG) thruster operated on the 10-microseconds timescales of pulsed plasma thrusters, but claimed high thrust density, high efficiency and low electrode erosion rates, which are more consistent with the deflagration mode of acceleration. Separate work on gas-fed pulsed plasma thrusters (PPTs) by Ziemer, et al. identified two separate regimes of performance. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency (ratio of jet kinetic energy to input electrical energy) as a function of mass bit. In the second regime at very low mass bits (termed Mode II), the efficiency increased with decreasing mass bit. Work by Poehlmann et al. and by Sitaraman and Raja sought to understand the performance of the CHENG thruster and the Mode I / Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes representing two distinct sets of solutions to the relevant conservation laws. These works studied the proposal that, depending upon the values of the various controllable parameters, the accelerator would operate in either the detonation or deflagration mode. In the present work, we propose a variation on the explanation for the differences in performance between the various pulsed plasma accelerators. Instead of treating the accelerator as if it were only operating in one mode or the other during a pulse, we model the initial stage of the discharge in all cases as an accelerating current sheet (detonation mode). If the current sheet reaches the exit of the accelerator before the discharge is completed, the acceleration mode transitions to the deflagration mode type found in the quasi-steady MPD thrusters. This modeling method is used to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗
NASA NTRS1992

Beneficiation of ilmenite from lumar analog

The results reported were obtained on a meteoric eucrite sample called Millbillillie Sample no. 173. Optical microscopy studies of the sample showed it to consist of ilmenite, troilite, and transparent gangue. The transparent gangue consisted of feldspar (anorthite), pyroxenes, olivines, and opaques. Troilite was present in minor quantities. Screen assay analyses of the 30, 100, 200, and 400 US mesh screen fractions showed that minor concentration of titanium occurred in the 200 x 400 and -400 mesh screen fractions. Scanning electron microscopy (SEM) studies of the bulk sample showed the presence of a variety of ilmenite grains, ranging from 50 microns down to less than 1 micron without any evidence of liberation. Electron Diffraction Scans (EDS) confirmed the ratio of Fe to Ti in the ilmenite grains. Dry magnetic separation in a Frantz Isodynamic Separator was found to be effective only at sizes finer than 150 microns (100 US mesh) and more so at 200 mesh (74 microns). In each case, dedusting of the sample to remove -400 mesh (-0.037 microns) fines was required. Liberation size was determined to be 200 mesh and finer. The highest grade concentrate assaying 3.45 percent Ti was produced by magnetic separation of the -200 + 400 mesh screen fraction assaying 0.44 Ti (from a -30 mesh sample) at a current setting of 0.35 AMP. This concentrate contained 21.2 percent of the Ti values in the screen fraction with 2.72 weight percent of feed to test. The results can be projected to a sample stage ground to -200 mesh. Magnetic separation of the 200 + 400 mesh (-0.074 + 0.037 microns) should produce a concentrate accounting for 1.41 weight percent of the feed. This concentrate will analyze 3.45 percent Ti and contain 10.3 percent of the Ti values in the feed. By changing the Frantz Magnetic Separator settings, a lower grade concentrate analyzing 0.98 percent Ti can be produced at an increased recovery of 25.4 percent. The concentrate weight will be 11.7 percent of the feed. It must be emphasized that improved grades and recoveries can be obtained with the -400 mesh fines. However, beneficiation of these extremely fine materials is not possible in a practical process scheme.

Ramadorai, G.↗
NASA NTRS1994

Deflections of anisotropic sandwich beams with variable face sheets and core thicknesses

A sandwich construction consists of a low-density core material with high strength face sheets bounded to the top and bottom surfaces. The construction has been widely used in the aerospace and marine industries due to its outstanding characteristics such as noise absorption, weight minimization, heat insulation, and better bending stiffness. In sandwich structures used in high-performance aircraft, the face sheets are often made of fiber-reinforced composite materials and the core is made of honeycomb. The structures may also have variable thickness so as to satisfy aerodynamic requirements. In the stress analysis, the constant-thickness face sheets are usually considered as membrane and the core is assumed to be inextensible but deformable in the thickness direction. The static behavior of variable-thickness, isotropic and homogeneous sandwich beams was successfully studied by employing a constant-thickness theory but allowing stiffnesses to vary in accordance with local thickness variations. It has been recently found in a refined theory that the analyses based on the constant thickness theory locally can lead to significant errors in structural responses if the sandwich beam is thickness-tapered and the cores are deformable in transverse shear. The errors arise mainly from two factors: (1) the transverse shear components of the membrane forces in the face sheets alter the transverse shears carried by the core; and (2) the face-sheet membrane strains arise from transverse shear deformation of the core. In practice the variable thickness may not only exist in core but also in face sheets. The thickness-variations may even be a type of step function. In this case the transverse shear stress in the face sheets and bending stress in the core should be taken into account in the refined theory mentioned. In the present study, energy principles are employed in deriving governing equations for general bending of anisotropic sandwich beams with variable thickness in both face sheets and cores. Solutions to these equations are based on a finite difference scheme. As an example in application, a simply supported thickness-tapered sandwich beam subject to a concentrated load at its center is considered. Let W' be the maximum deflection of the beam in which face sheets are considered as membrane, while W'' is that based on using the modified refined theory. It is found that W' is always larger than W'', however, the magnitude of (W'- W'') appears to be insensitive to the change of the taper of the beam.

Lu, Chu-Ho↗
NASA NTRS2013

Thermal Performance of ATLAS Laser Thermal Control System Demonstration Unit

The second Ice, Cloud, and Land Elevation Satellite mission currently planned by National Aeronautics and Space Administration will measure global ice topography and canopy height using the Advanced Topographic Laser Altimeter System {ATLAS). The ATLAS comprises two lasers; but only one will be used at a time. Each laser will generate between 125 watts and 250 watts of heat, and each laser has its own optimal operating temperature that must be maintained within plus or minus 1 degree Centigrade accuracy by the Laser Thermal Control System (LTCS) consisting of a constant conductance heat pipe (CCHP), a loop heat pipe (LHP) and a radiator. The heat generated by the laser is acquired by the CCHP and transferred to the LHP, which delivers the heat to the radiator for ultimate rejection. The radiator can be exposed to temperatures between minus 71 degrees Centigrade and minus 93 degrees Centigrade. The two lasers can have different operating temperatures varying between plus 15 degrees Centigrade and plus 30 degrees Centigrade, and their operating temperatures are not known while the LTCS is being designed and built. Major challenges of the LTCS include: 1) A single thermal control system must maintain the ATLAS at 15 degrees Centigrade with 250 watts heat load and minus 71 degrees Centigrade radiator sink temperature, and maintain the ATLAS at plus 30 degrees Centigrade with 125 watts heat load and minus 93 degrees Centigrade radiator sink temperature. Furthermore, the LTCS must be qualification tested to maintain the ATLAS between plus 10 degrees Centigrade and plus 35 degrees Centigrade. 2) The LTCS must be shut down to ensure that the ATLAS can be maintained above its lowest desirable temperature of minus 2 degrees Centigrade during the survival mode. No software control algorithm for LTCS can be activated during survival and only thermostats can be used. 3) The radiator must be kept above minus 65 degrees Centigrade to prevent ammonia from freezing using no more than 135 watts of heater power. 4) The LHP reservoir control heater power is limited to 15 watts with a 70 percent duty cycle. 5) The voltage of the power supply can vary between 26 volts direct current and 34 volts direct current during the spacecraft lifetime. A design analysis shows that a single LTCS can satisfy these requirements. However, shutdown of· the LHP is particularly challenging and the shutdown heater must be wired in series with two reservoir thermostats and two CCHP thermostats at different set points. An LTCS demonstration unit has been tested to verify these performance characteristics experimentally prior to proceeding to the final LTCS design and fabrication. Test results showed that the LHP shutdown scheme would be able to shut down the LHP as designed and the reservoir control heater can maintain the ATLAS mass simulator within the plus or minus 1 degrees Centigrade accuracy under various combinations of the heat load, sink temperature, and power supply voltage.

Ku, Jentung↗
NASA NTRS1993

Analysis of the tropospheric water distribution during FIRE 2

The Penn State/NCAR mesoscale model, as adapted for use at ARC, was used as a testbed for the development and validation of cloud models for use in General Circulation Models (GCM's). This modeling approach also allows us to intercompare the predictions of the various cloud schemes within the same dynamical framework. The use of the PSU/NCAR mesoscale model also allows us to compare our results with FIRE-II (First International Satellite Cloud Climatology Project Regional Experiment) observations, instead of climate statistics. Though a promising approach, our work to date revealed several difficulties. First, the model by design is limited in spatial coverage and is only run for 12 to 48 hours at a time. Hence the quality of the simulation will depend heavily on the initial conditions. The poor quality of upper-tropospheric measurements of water vapor is well known and the situation is particularly bad for mid-latitude winter since the coupling with the surface is less direct than in summer so that relying on the model to spin-up a reasonable moisture field is not always successful. Though one of the most common atmospheric constituents, water vapor is relatively difficult to measure accurately, especially operationally over large areas. The standard NWS sondes have little sensitivity at the low temperatures where cirrus form and the data from the GOES 6.7 micron channel is difficult to quantify. For this reason, the goals of FIRE Cirrus II included characterizing the three-dimensional distribution of water vapor and clouds. In studying the data from FIRE Cirrus II, it was found that no single special observation technique provides accurate regional distributions of water vapor. The Raman lidar provides accurate measurements, but only at the Hub, for levels up to 10 km, and during nighttime hours. The CLASS sondes are more sensitive to moisture at low temperatures than are the NWS sondes, but the four stations only cover an area of two hundred kilometers on a side. The aircraft give the most accurate measurements of water vapor, but are limited in spatial and temporal coverage. This problem is partly alleviated by the use of the MAPS analyses, a four-dimensional data assimilation system that combines the previous 3-hour forecast with the available observations, but its upper-level moisture analyses are sometimes deficient because of the vapor measurement problem. An attempt was made to create a consistent four-dimensional description of the water vapor distribution during the second IFO by subjectively combining data from a variety of sources, including MAPS analyses, CLASS sondes, SPECTRE sondes, NWS sondes, GOES satellite analyses, radars, lidars, and microwave radiometers.

Westphal, Douglas L.↗
NASA NTRS2000

Power Systems Evaluated for Solar Electric Propulsion Vehicles

Solar electric propulsion (SEP) mission architectures are applicable to a wide range NASA missions including the robotic exploration of the outer planets in the next decade and the human exploration of Mars within the next 2 decades. SEP enables architectures that are very mass efficient with reasonable power levels (1-MW class) aerobrake and cryogenic upper-stage transportation technologies are utilized. In this architecture, the efficient SEP stage transfers the payload from low Earth orbit (LEO) High Energy Elliptical Parking Orbit (HEEPO) within a period of 6 to 12 months. highthrust, cryogenic upper stage and payload then separate from the SEP vehicle for injection to the planetary target, allowing for fast heliocentric trip times. This mission architecture offers a potential reduction in mass to LEO in comparison to alternative all-chemical nuclear propulsion schemes. Mass reductions may allow launch vehicle downsizing enable missions that would have been grounded because of cost constraints. The preceding figure illustrates a conceptual SEP stage design for a human Mars mission. Researchers at the NASA Glenn Research Center at Lewis Field designed conceptual SEP vehicle, conceived the mission architecture to use this vehicle, and analyzed the vehicle s performance. This SEP stage has a dry mass of 35 metric tons (MT), 40 MT of xenon propellant, and a photovoltaic array that spans 110 m, providing power to a cluster of eight 100-kW Hall thrusters. The stage can transfer an 80-MT payload and upper stage to the desired HEEPO. Preliminary packaging studies show this space-station-class SEP vehicle meets the proposed "Magnum" launch vehicle and volume requirements with considerable margin. An SEP vehicle for outer planetary missions, such as the Europa Mapper Mission, would be dramatically smaller than human Mars mission SEP stage. In this mission architecture, the SEP power system with the payload to provide spacecraft power throughout the mission. Several photovoltaic array design concepts were considered for the SEP vehicle power system for the human mission to Mars. These include a space station derivative, a SCARLET (Solar Concentrator Arrays with Refractive Linear Element Technology) derivative, and a hybrid inflatable-deployable thin polymer membrane array with thin-film solar cells (as shown in the concept illustration). This concept is based on a design developed for the Next Generation Space Telescope Sun shield. The array is divided into 16 independent electrical sections with 500-V, negative-grounded solar cell strings. The power system employs a channelized, 500-Vdc power management and distribution (PMAD) architecture with lithium ion batteries for energy storage for vehicle and payload secondary loads (the high-power Hall thrusters do not operate in eclipse periods). The 500-V PMAD voltage permits "direct-drive" thruster operation, greatly reducing the power processing unit size, complexity, and power loss. Similar power system architecture, designs, and technology are assumed for the Europa Mapper Mission SEP vehicle. The primary exceptions are that the photovoltaic array is assumed to consist of two rectangular wings and that the power system rating is 15 kW in Earth orbit and 200 W at Europa. To size the SEP vehicle power system, a dedicated Fortran code was developed to predict detailed power system performance, mass, and thermal control requirements. This code also modeled all the relevant Earth orbit environments; that is, the particulate radiation, plasma, meteoroids and debris, ultraviolet radiation, contamination, and thermal conditions. Analysis results for the Human Mars Mission SEP vehicle show a power system mass of 9-MT and photovoltaic array area of 5800-square meters for the thin-membrane design concept with CuInS2 thin-film cells. Power processing unit input power for a thin-membrane array design with three-junction, amorphous SiGe solar cells is shown in the graph. Power falls off rapidly inhe first weeks of the mission because of light-induced (Staebler-Wronksi) solar cell losses. During the next 200 days, power decreases steadily as the SEP stage spirals through the proton belts and sustains the bulk of the mission radiation damage. Once the vehicle apogee is above approximately four Earth radii, little additional degradation is incurred. From 400 to 800 days, a 1100-km "parking" orbit is maintained to await the next payload transfer opportunity. This orbit is below the main proton belt, and thus, little radiation dose is accumulated during this time period. During the second LEO-to-HEEPO transfer, power degrades somewhat further, but power requirements are still met. In comparison, the Europa Mapper SEP vehicle power system had a mass of 150 kg and a thin membrane array area of 100 square meters.

Kerslake, Thomas W.↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗