Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local minimum”

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 919 records · Page 51

Projecting climate change in South America using variable‐resolution Community Earth System Model: An application to Chile

Abstract We introduce variable‐resolution enabled Community Earth System Model (VR‐CESM) results simulating historical and future climate conditions at 28 km over South America and 14 km over the Andes. Three 30‐year simulations are performed: a historic (1985–2014), a near future (2030–2059), and an end‐century (2070–2099) simulation under the RCP8.5 scenario. Historic results compare favourably to several temperature and precipitation reanalysis products, though local biases are present, particularly during austral summer. Future simulations highlight broad warming patterns (+3–6°C by end‐century) and heterogeneous precipitation responses across South America that qualitatively agree with prior modelling efforts. Our results reveal that the interaction between temperature and precipitation changes produce shifts in several Köppen–Geiger climates. Notable changes include the near‐elimination of the Andean Tundra or Alpine climates, a 15% decrease in Tropical Rainforests and a Tropical Savannah expansion of 20%. To provide a regionally focused analysis of projected climate change and to illustrate the benefits of variable resolution modelling, we analyse changes in the magnitude and trend in seasonal and daily temperature and precipitation in Chile. We also examined several metrics [e.g., snow water equivalent (SWE), temperatures on wet days, and days below 0°C] to evaluate potential impacts of climate change on the Chilean cryosphere between the end‐of‐century and historic periods, finding wide‐ranging indications of cryospheric decline. These changes are interpreted through reductions in the timing (1–2.5 months earlier peak SWE) and magnitude (200–1,000 mm SWE decreases) of water stored as snow in the Andes, a 10–30% decrease in number of cool season wet days with temperatures below 1°C, and 50–200 fewer days (annually) with minimum temperatures below 0°C. Our aim in producing a high‐resolution dataset of climate projections from VR‐CESM is to support analyses of climate change throughout South America but especially in vulnerable montane regions and to provide additional results for comparison with previous, ongoing, and upcoming modelling efforts.

variable-resolution global climate models↗

Analysis and Design of Crew Sleep Station for ISS

This paper details the analysis and design of the Temporary Sleep Station (TeSS) environmental control system for International Space Station (ISS). The TeSS will provide crewmembers with a private and personal space, to accommodate sleeping, donning and doffing of clothing, personal communication and performance of recreational activities. The need for privacy to accommodate these activities requires adequate ventilation inside the TeSS. This study considers whether temperature, carbon dioxide, and humidity within the TeSS remain within crew comfort and safety levels for various expected operating scenarios. Evaluation of these scenarios required the use and integration of various simulation codes. An approach was adapted for this study, whereby results from a particular code were integrated with other codes when necessary. Computational Fluid Dynamics (CFD) methods were used to evaluate the flow field inside the TeSS, from which local gradients for temperature, velocity, and species concentration such as CO (sub 2) could be determined. A model of the TeSS, containing a human, as well as equipment such as a laptop computer, was developed in FLUENT, a finite-volume code. Other factors, such as detailed analysis of the heat transfer through the structure, radiation, and air circulation from the TeSS to the US Laboratory Aisle, where the TeSS is housed, were considered in the model. A complementary model was developed in G189A, a code which has been used by NASA/JSC for environmental control systems analyses since the Apollo program. Boundary conditions were exchanged between the FLUENT and G189A TeSS models. G189A provides human respiration rates to the FLUENT model, while the FLUENT model provides local convective heat transfer coefficients to G189A model. An additional benefit from using an approach with both a systems simulation and CFD model, is the capability to verify the results of each model by comparison to the results of the other model. The G189A and FLUENT models were used to evaluate various ventilation designs for the TeSS over a range of operating conditions with varying crew metabolic load, equipment operating modes, ventilation flow rates, and with the TeSS doors open and closed. Results from the study were instrumental in the optimization of a design for the TeSS ventilation hardware. A special case was considered where failure of the TeSS ventilation system occurred. In this case, a study was conducted in order to determine the time required for the CO (sub 2) concentration inside the TeSS to increase to ISS limit values under transient conditions. A lumped-capacitance code, SINDA-FLUINT was used in this case to provide accurate predictions of the human reaction to the TeSS cabin conditions including core and skin temperatures and body heat storage. A simple two-dimensional CFD model of a crewmember inside the TeSS was developed in FLUENT in order to determine the volume envelope of the respired air from the human, which maintained a minimum velocity profile. This volume was then used in the SINDA-FLUINT model to facilitate the calculations of CO (sub 2) concentrations, dry bulb temperatures and humidity levels inside the TeSS.

Keener, John F.↗

Evaluating space station applications of automation and robotics technologies from a human productivity point of view

The role that automation, robotics, and artificial intelligence will play in Space Station operations is now beginning to take shape. Although there is only limited data on the precise nature of the payoffs that these technologies are likely to afford there is a general consensus that, at a minimum, the following benefits will be realized: increased responsiveness to innovation, lower operating costs, and reduction of exposure to hazards. Nevertheless, the question arises as to how much automation can be justified with the technical and economic constraints of the program? The purpose of this paper is to present a methodology which can be used to evaluate and rank different approaches to automating the functions and tasks planned for the Space Station. Special attention is given to the impact of advanced automation on human productivity. The methodology employed is based on the Analytic Hierarchy Process. This permits the introduction of individual judgements to resolve the confict that normally arises when incomparable criteria underly the selection process. Because of the large number of factors involved in the model, the overall problem is decomposed into four subproblems individually focusing on human productivity, economics, design, and operations, respectively. The results from each are then combined to yield the final rankings. To demonstrate the methodology, an example is developed based on the selection of an on-orbit assembly system. Five alternatives for performing this task are identified, ranging from an astronaut working in space, to a dexterous manipulator with sensory feedback. Computational results are presented along with their implications. A final parametric analysis shows that the outcome is locally insensitive to all but complete reversals in preference.

Bard, J. F.↗

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗

Hydrogen Biogeochemistry in Anaerobic and Photosynthetic Ecosystems

The simple biochemistry of molecular hydrogen is central to a large number of microbial processes, affecting the interaction of organisms with each other and with the environment. In anoxic sediments, a great majority of microbial redox processes involve hydrogen as a reactant, product or potential by-product. Accordingly, the energetics (thermodynamics) of each of these processes is affected by variations in local H2 concentrations. It has long been established that this effect is important in governing microbe-microbe interactions and there are multiple demonstrations that "interspecies hydrogen transfer" can alter the products of, inhibit/stimulate, or even reverse microbial metabolic reactions. In anoxic sediments, H2 concentrations themselves are thought to be controlled by the thermodynamics of the predominant H2-consuming microbial process. In sediments from Cape Lookout Bight, this relationship quantitatively describes the co-variation of H2 concentrations with temperature (for methanogens and sulfate reducers) and with sulfate concentration (for sulfate reducers). The quantitative aspect is import= for two reasons: 1) it permits the modeling of H2-sensitive biogeochemistry, such as anaerobic methane oxidation or pathways of organic matter remineralization, as a function of environmental controls; 2) for such a relationship to be observed requires that intracellular biochemistry and bioenergetics are being directly expressed in a component of the extracellular medium. H2 could therefore be utilized a non-invasive probe of cellular energetic function in intact microbial ecosystems. Based on the latter principle we have measured down-core profiles of H2 and other relevant physico-chemical parameters in order to calculate the metabolic energy yields (DG) that support microbial metabolism in Cape Lookout Bight sediments. Methanogens in this system apparently function with energy yields significantly smaller than the minimum requirements suggested by pure culture studies. Our recent work has extended the study of hydrogen to cyanobacterial mat communities. The large amounts of reducing power generated during photosynthetic activity carry the potential to contribute a swamping term to the H2 economy of the anaerobic microbial populations within the mat - and thereby to alter the population structure and biogeochemical function of the mat as a whole. In hypersaline microbial mats, we observe a distinct diel cycle in H2 production and a substantial corresponding flux. On an early Earth dominated by microbial mats, this transmission of photosynthetic reducing power may have carried important implications for both biospheric and atmospheric evolution.

Hoehler, Tori M.↗

Using Satellite Imagery to Identify Tornado Damage Tracks and Recovery from the April 27, 2011 Severe Weather Outbreak

Emergency response to natural disasters requires coordination between multiple local, state, and federal agencies. Single, relatively weak tornado events may require comparatively simple response efforts; but larger "outbreak" events with multiple strong, long-track tornadoes can benefit from additional tools to help expedite these efforts. Meteorologists from NOAA's National Weather Service conduct field surveys to map tornado tracks, assess damage, and determine the tornado intensity following each event. Moderate and high resolution satellite imagery can support these surveys by providing a high-level view of the affected areas. Satellite imagery could then be used to target areas for immediate survey or to corroborate the results of the survey after it is completed. In this study, the feasibility of using satellite imagery to identify tornado damage tracks was determined by comparing the characteristics of tracks observed from low-earth orbit to tracks assessed during the official NWS storm survey process. Of the 68 NWS confirmed centerlines, 24 tracks (35.3%) could be distinguished from other surface features using satellite imagery. Within each EF category, 0% of EF-0, 3% of EF-1, 50% of EF-2, 77.7% of EF-3, 87.5% of EF-4 and 100% of EF-5 tornadoes were detected. It was shown that satellite data can be used to identify tornado damage tracks in MODIS and ASTER NDVI imagery, where damage to vegetation creates a sharp drop in values though the minimum EF-category which can be detected is dependent upon the type of sensor used and underlying vegetation. Near-real time data from moderate resolution sensors compare favorably to field surveys after the event and suggest that the data can provide some value in the assessment process.

Cole, Tony A.↗

Using Satellite Imagery to Identify Tornado Damage Tracks and Recovery from the April 27, 2011 Severe Weather Outbreak

Emergency response to natural disasters requires coordination between multiple local, state, and federal agencies. Single, relatively weak tornado events may require comparatively simple response efforts; but larger "outbreak" events with multiple strong, long-track tornadoes can benefit from additional tools to help expedite these efforts. Meteorologists from NOAA's National Weather Service conduct field surveys to map tornado tracks, assess damage, and determine the tornado intensity following each event. Moderate and high resolution satellite imagery can support these surveys by providing a high-level view of the affected areas. Satellite imagery could then be used to target areas for immediate survey or to corroborate the results of the survey after it is completed. In this study, the feasibility of using satellite imagery to identify tornado damage tracks was determined by comparing the characteristics of tracks observed from low-earth orbit to tracks assessed during the official NWS storm survey process. Of the 68 NWS confirmed centerlines, 24 tracks (35.3%) could be distinguished from other surface features using satellite imagery. Within each EF category, 0% of EF-0, 3% of EF-1, 50% of EF-2, 77.7% of EF-3, 87.5% of EF-4 and 100% of EF-5 tornadoes were detected. It was shown that satellite data can be used to identify tornado damage tracks in MODIS and ASTER NDVI imagery, where damage to vegetation creates a sharp drop in values though the minimum EF-category which can be detected is dependent upon the type of sensor used and underlying vegetation. Near-real time data from moderate resolution sensors compare favorably to field surveys after the event and suggest that the data can provide some value in the assessment process.

Molthan, Andrew L.↗

Tropical Cyclone Diurnal Cycle as Observed by TRMM

Using infrared satellite data, previous work has shown a consistent diurnal cycle in the pattern of cold cloud tops around mature tropical cyclones. In particular, an increase in the coverage by cold cloud tops often occurs in the inner core of the storm around the time of sunset and subsequently propagates outward to several hundred kilometers over the course of the following day. This consistent cycle may have important implications for structure and intensity changes of tropical cyclones and the forecasting of such changes. Because infrared satellite measurements are primarily sensitive to cloud top, the goal of this study is to use passive and active microwave measurements from the Tropical Rainfall Measurement Mission (TRMM) Microwave Imager (TMI) and Precipitation Radar (PR), respectively, to examine and better understand the tropical cyclone diurnal cycle throughout a larger depth of the storm's clouds. The National Hurricane Center's best track dataset was used to extract all PR and TMI pixels within 1000 km of each tropical cyclone that occurred in the Atlantic basin between 1998-2011. Then the data was composited according to radius (100-km bins from 0-1000 km) and local standard time (LST; 3-hr bins). Specifically, PR composites involved finding the percentage of pixels with reflectivity greater than or equal to 20 dBZ at various heights (i.e., 2-14 km in increments of 2 km) as a function of radius and time. The 37- and 85- GHz TMI channels are especially sensitive to scattering by precipitation-sized ice in the mid to upper portions of clouds. Hence, the percentage of 37- and 85-GHz polarization corrected temperatures less than various thresholds were calculated using data from all storms as a function of radius and time. For 37 GHz, thresholds of 260 K, 265 K, 270 K, and 275 K were used, and for 85 GHz, thresholds of 200-270 K in increments of 10 K were utilized. Note that convection forced by the interactions of a tropical cyclone with land (e.g., due to frictional convergence) may disrupt the natural convective cycle of a cyclone. Hence, only data pertaining to storms whose centers were greater than 300 km from land were included in the composites. Early results suggest the presence of a diurnal cycle in the PR composites of all Atlantic basin tropical cyclones from a height of 2-12 km from approximately 0-400 km radius, but the cycle is most apparent above 6 km. At a height of 8 km, there is a peak (minimum) in the percentage of PR pixels greater than or equal to 20 dBZ near 0 (21) LST in the inner core with some indication that this signal propagates outward with time. In contrast, the 37- and 85-GHz composites show little indication of a diurnal cycle at any radii, regardless of the threshold used. Ongoing work with this project will involve sub-setting the composites according to storm intensity to see if the diurnal cycle varies with storm strength. Moderate to strong vertical wind shear often leads to asymmetries in tropical cyclone convection and may disrupt the cyclone's natural diurnal cycle. Therefore, wind shear thresholds will be applied to the composites to determine if the diurnal cycle becomes more apparent in a low shear environment. Finally, other work to be completed will involve developing composites for other tropical cyclone basins, including the East Pacific, Northwest Pacific, South Pacific, and Indian Ocean.

tropical cyclones↗

A View of the Lunar Interior Through Lunar Laser Range Analysis

Laser ranges between observatories on the Earth and retroreflectors on the Moon started in 1969 and continue to the present. Recent range accuracies are 2 cm while earliest ranges are an C, order of magnitude less certain. Four retroreflectors are ranged: three located at the Apollo 11, 14, and 15 sites and one on the Lunakhod 2 rover. Accurate analysis of the range data determines a number of lunar science parameters. The lunar interior variables include a fluid core parameter. The Lunar Laser Ranging effort is reviewed elsewhere. Many parameters are detected through their influence on rotation. Also detected are solid-body tides and accurate selenocentric reflector locations. Determined through the rotation are moment-of-inertia differences, gravitational harmonics, potential Love number, and dissipation effects due to tides and molten core. The rotation of the Moon is not at its minimum energy state; some recently active process has caused free librations. The moment differences contributed to the recent improvement of the Moon's moment of inertia from the Lunar Prospector gravity field. The Love numbers provide bulk elastic properties. Future possibilities for measurement include oblateness of the core-mantle boundary and core moment. A study of dissipation signatures in the rotation determines tidal Q vs. frequency and concludes that the Moon has a molten core. At 1 month the tidal Q is 37 and at 1 yr it is 60. The core radius is < or = 352 km for Fe and < or = 374 km for the Fe-FeS eutectic. The core detection exceeds 3x its uncertainty. The spin of the core is not aligned with the spin of the mantle and torque arises from the velocity difference at the boundary. Yoder's turbulent boundary layer theory is used to compute the radii. The present heat generation from tides and core interaction is minor compared to radiogenic heating. The heating for ancient times is more interesting. Peale and Cassen investigated lunar tidal heating while the lunar orbit expanded due to tides on Earth. Their calculations predate the measurement of Q and should be multiplied by 3.45 to match the lunar-laser-determined Love number and monthly Q. Tidal-heating computations depend on how fast the lunar orbit evolved and whether the tidal dissipation is localized. Neither is known, but plausible assumptions lead to early central region temperature increases of several hundred degrees. Most of the energy is deposited early in the Moon's history. The turbulent boundary layer theory allows a prediction of energy dissipated at the core-mantle boundary during orbit evolution. Under the assumption that the properties of the early core are the same as at present, the energy dissipated by core-mantle interaction is about the same as for tidal dissipation, but it is deposited in a smaller volume. This source of energy is capable of promoting convection in an early fluid core and driving a dynamo. This is a transient phase with a duration depending on the rate of orbit evolution. Plausible assumptions lead to a duration of a few hundred million years. Thus, the remnant magnetization of many lunar rocks is compatible with a brief global magnetic field powered by dynamical energy dissipation. Analysis of the lunar laser ranges is providing information on lunar geophysics. Future data will improve accuracies of present solution parameters, and several more interior effects should be detectable.

Williams, J. G.↗

Chronostratigraphy of talus flatirons and piedmont alluvium along the Book Cliffs, Utah – Testing models of dryland escarpment evolution

Research on the evolution of dryland escarpments and what drives and controls their erosional retreat have been limited by a lack of chronologically constrained records. Prior investigations of escarpments in the southwestern U.S. and globally garner two endmember conceptual models -- one focused on bottom-up baselevel drivers and autogenic variations in piedmont and toeslope erosion processes, the other focused on the role of climate as a top-down forcing mechanism on sediment production and transport from cliffs. In this work, we test these conceptual models along a section of the Book Cliffs in central Utah, where climate has varied strongly and baselevel has fallen significantly over the late Quaternary. Four generations of talus flatirons and piedmont terraces are preserved in the study area we date them by optically stimulated luminescence and 10Be-exposure techniques. OSL depositional ages cluster into four generations of deposits with mean ages and mean standard errors of 117.8 ± 6.7 ka, 83.4 ± 3.3 ka, 52.4 ± 1.8 ka, and 5.6 ± 0.3 ka. 10Be-exposure results from talus boulders protruding from flatirons confirm that the youngest generation of talus is Holocene, with recent rock-fall boulders still being received on the landform top. Older flatiron generations produce incoherent, minimum exposure ages reflecting an increasing dominance of in situ boulder weathering and a steady-state denudation rate of ~45 mm/ky after tens of thousands of years. Based on the chronostratgraphy and correlation of talus and piedmont deposits, we interpret climate as the primary control on escarpment evolution. Depositional ages do not correspond to glacial epochs as historically presumed, but instead to episodes of intermediate and unstable or variable Pleistocene climate during Marine Isotope Stage-3 and 5 a/b. Results support the significance of high climate variability at millennial timescales in driving both sediment production by mass wasting from the cliff as well as sediment storage and reworking along the piedmont. In contrast, wetter climates of glacial epochs apparently caused greater erosion and sediment transport to trunk drainages. Our chronology also reveals a systematic decrease in deposit age with increasing distance up the piedmont from the regional baselevel of the Price River, consistent with transient fluvial incision and terrace abandonment. The Book Cliffs record notably does not correlate with similarly constrained records in Spain or the contrasting record from the Negev, yet they all share the feature of enhanced talus deposition and cliff retreat corresponding to climate instability rather than in response to steady glacial-climate conditions. Thus, we outline a model of escarpment evolution and retreat governed largely by top-down climate-forcing of sediment production and transport during climate variability, with local changes in baselevel and erosional processes on the piedmont being a modifying, secondary factor.

58 GEOSCIENCES↗

Modeling of a Turbofan Engine with Ice Crystal Ingestion in the NASA Propulsion System Laboratory

The main focus of this study is to apply a computational tool for the flow analysis of the turbine engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) at NASA Glenn Research Center. The PSL has been used to test a highly instrumented Honeywell ALF502R-5A (LF11) turbofan engine at simulated altitude operating conditions. Test data analysis with an engine cycle code and a compressor flow code was conducted to determine the values of key icing parameters, that can indicate the risk of ice accretion, which can lead to engine rollback (un-commanded loss of engine thrust). The full engine aerothermodynamic performance was modeled with the Honeywell Customer Deck specifically created for the ALF502R-5A engine. The mean-line compressor flow analysis code, which includes a code that models the state of the ice crystal, was used to model the air flow through the fan-core and low pressure compressor. The results of the compressor flow analyses included calculations of the ice-water flow rate to air flow rate ratio (IWAR), the local static wet bulb temperature, and the particle melt ratio throughout the flow field. It was found that the assumed particle size had a large effect on the particle melt ratio, and on the local wet bulb temperature. In this study the particle size was varied parametrically to produce a non-zero calculated melt ratio in the exit guide vane (EGV) region of the low pressure compressor (LPC) for the data points that experienced a growth of blockage there, and a subsequent engine called rollback (CRB). At data points where the engine experienced a CRB having the lowest wet bulb temperature of 492 degrees Rankine at the EGV trailing edge, the smallest particle size that produced a non-zero melt ratio (between 3 percent - 4 percent) was on the order of 1 micron. This value of melt ratio was utilized as the target for all other subsequent data points analyzed, while the particle size was varied from 1 micron - 9.5 microns to achieve the target melt ratio. For data points that did not experience a CRB which had static wet bulb temperatures in the EGV region below 492 degrees Rankine, a non-zero melt ratio could not be achieved even with a 1 micron ice particle size. The highest value of static wet bulb temperature for data points that experienced engine CRB was 498 degrees Rankine with a particle size of 9.5 microns. Based on this study of the LF11 engine test data, the range of static wet bulb temperature at the EGV exit for engine CRB was in the narrow range of 492 degrees Rankine - 498 degrees Rankine , while the minimum value of IWAR was 0.002. The rate of blockage growth due to ice accretion and boundary layer growth was estimated by scaling from a known blockage growth rate that was determined in a previous study. These results obtained from the LF11 engine analysis formed the basis of a unique “icing wedge.”

Turbo engine↗

Circuit complexity and functionality: A statistical thermodynamics perspective

Circuit complexity, defined as the minimum circuit size required for implementing a particular Boolean computation, is a foundational concept in computer science. Determining circuit complexity is believed to be a hard computational problem. Recently, in the context of black holes, circuit complexity has been promoted to a physical property, wherein the growth of complexity is reflected in the time evolution of the Einstein-Rosen bridge (“wormhole”) connecting the two sides of an anti-de Sitter “eternal” black hole. Here, we are motivated by an independent set of considerations and explore links between complexity and thermodynamics for functionally equivalent circuits, making the physics-inspired approach relevant to real computational problems, for which functionality is the key element of interest. In particular, our thermodynamic framework provides an alternative perspective on the obfuscation of programs of arbitrary length—an important problem in cryptography—as thermalization through recursive mixing of neighboring sections of a circuit, which can be viewed as the mixing of two containers with “gases of gates.” This recursive process equilibrates the average complexity and leads to the saturation of the circuit entropy, while preserving functionality of the overall circuit. The thermodynamic arguments hinge on ergodicity in the space of circuits which we conjecture is limited to disconnected ergodic sectors due to fragmentation. The notion of fragmentation has important implications for the problem of circuit obfuscation as it implies that there are circuits of same size and functionality that cannot be connected via a polynomial number of local moves. Furthermore, we argue that fragmentation is unavoidable unless the complexity classes NP and coNP coincide, a statement that implies the collapse of the polynomial hierarchy of computational complexity theory to its first level.

Science & Technology - Other Topics↗

The Geoscience Laser Altimeter System (GLAS) for the ICESAT Mission

Accurate measurements of surface heights and atmospheric backscatter have been demonstrated with the SLA, MOLA and LITE space lidar. Recent MOLA measurements of the Mars surface have 40 cm resolution and have reduced the global uncertainty in Mars topography from a few km to approx. 10 m. GLAS is a next generation lidar being developed as part of NASA's Icesat Mission for Earth orbit . The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, determine the height profiles of the Earth's land topography, and profile the vertical backscatter of clouds and aerosols on a global scale. GLAS will fly on a small dedicated spacecraft in a polar orbit at 598 km altitude with an inclination of 94 degrees. GLAS is scheduled to launch in summer 2001 and to operate continuously for a minimum of 3 years with a goal of 5 years. The primary mission for GLAS is to measure the seasonal and annual changes in the heights of the Greenland and Antarctic ice sheets. GLAS will measure the vertical distance to the ice sheet from orbit with 1064 nm pulses from a Nd:Yag laser at 40 Hz. Each 5 nsec wide laser pulse is used for a single range measurement. When over land GLAS will profile the heights of the topography and vegetation. The GLAS receiver uses a I m diameter telescope and a Si APD detector. The detector signal is sampled by an all digital receiver which records each surface echo waveform with I nsec resolution and a stored echo record lengths of either 200, 400, or 600 samples. Analysis of the echo waveforms within the instrument permits discrimination between cloud and surface echoes. Ground based echo analysis permits precise ranging, determining the roughness or slopes of the surface as well as the vertical distributions of vegetation illuminated by the laser, Errors in knowledge of the laser beam pointing angle can bias height measurements of sloped surfaces. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam centroid to about 8 urad is required to achieve 10 cm height accuracy. GLAS uses a stellar reference system (SRS) to determine the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith whose measurements are combined with a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser pulse is measured with a laser reference system (LRS). Optically measuring each laser far field pattern relative to the star camera and gyroscope permits the angular offsets of each laser pulse to be determined. GLAS will also determine the vertical distributions of clouds and aerosols by measuring atmospheric backscatter profiles at both 1064 and 532 nm. The 1064 nm measurements use an analog detector and profile the height and vertical structure of thicker clouds. Measurements at 532 nm use new highly sensitive photon counting detectors, and measure the height distributions of very thin clouds and aerosol layers. With averaging these can be used to determine the height of the planetary boundary layer. The instrument design and expected performance will be discussed.

Abshire, James B.↗

Design and Field Test of a Mass Efficient Crane for Lunar Payload Handling and Inspection: The Lunar Surface Manipulation System

Devices for lifting, translating and precisely placing payloads are critical for efficient Earthbased construction operations. Both recent and past studies have demonstrated that devices with similar functionality will be needed to support lunar outpost operations. Lunar payloads include: a) prepackaged hardware and supplies which must be unloaded from landers and then accurately located at their operational site, b) sensor packages used for periodic inspection of landers, habitat surfaces, etc., and c) local materials such as regolith which require grading, excavation and placement. Although several designs have been developed for Earth based applications, these devices lack unique design characteristics necessary for transport to and use on the harsh lunar surface. These design characteristics include: a) composite components, b) compact packaging for launch, c) simple in-field reconfiguration and repair, and d) support for tele-operated or automated operations. Also, in contrast to Earth-based construction, where special purpose devices dominate a construction site, a lunar outpost will require versatile devices which provide operational benefit from initial construction through sustained operations. This paper will detail the design of a unique, high performance, versatile lifting device designed for operations on the lunar surface. The device is called the Lunar Surface Manipulation System to highlight the versatile nature of the device which supports conventional cable suspended crane operations as well as operations usually associated with a manipulator such as precise positioning where the payload is rigidly grappled by a tool attached to the tip of the device. A first generation test-bed to verify design methods and operational procedures is under development at the NASA Langley Research Center and recently completed field tests at Moses Lake Washington. The design relied on non-linear finite element analysis which is shown to correlate favorably with laboratory experiments. A key design objective, reviewed in this paper, is the device s simplicity, resulting from a focus on the minimum set of functions necessary to perform payload offload. Further development of the device has the potential for significant mass savings, with a high performance device incorporating composite elements estimated to have a mass less than 3% of the mass of the maximum lunar payload lifted at the tip. The paper will conclude with future plans for expanding the operational versatility of the device.

Doggett, William R.↗

Effect of static local distortions vs. dynamic motions on the stability and band gaps of cubic oxide and halide perovskites

Ternary ABX 3 perovskites made of corner-sharing BX 6 octahedra have long featured prominently in solid-state chemistry and condensed matter physics. Still, the joint understanding of their two main subgroups—halides and oxides—has not been fully developed. Indeed, unlike the case in simpler compounds having a single, robust repeated motif (“monomorphous”), certain cubic perovskites can manifest a non-thermal (= intrinsic) distribution of local motifs (“polymorphous networks”). Such static deformations can include positional degrees of freedom (e.g., atomic displacements and octahedral tilting) or magnetic moment degrees of freedom in paramagnets. Unlike thermal motion, such static distortions do not time-average to zero, being an expression of the intrinsic symmetry breaking preference of the chemical bonding. Here, the present study compares electronic structure features of oxide and halide perovskites starting from the static polymorphous distribution of motifs described by Density Functional Theory (DFT) minimization of the internal energy, continuing to finite temperature thermal disorder modeled via finite temperature DFT molecular dynamics. We find that (i) different oxide vs. halide ABX 3 compounds adopt different energy-lowering symmetry-breaking modes. The calculated pair distribution function (PDF) of SrTiO 3 from the first-principles agrees with recently measured PDF. (ii) In both oxides and halides, such static distortions lead to band gap blueshifts with respect to undistorted cubic Pm-3m structure. (iii) For oxide perovskites, high-temperature molecular dynamics simulations initiated from the statically distorted polymorphous structures reveal that the thermally-induced distortions can lead to a band gap redshift. (iv) In contrast, for cubic halide perovskite CsPbI 3 , both the intrinsic distortions and the thermal distortions contribute in tandem to band gap blueshift, the former, intrinsic effect being dominant. (v) In the oxide SrTiO 3 and CaTiO 3 (but not in halide) perovskites, octahedral tilting leads to the emergence of a distinct Γ–Γ direct band gap component as a secondary valley minimum to the well-known indirect R–Γ gap. Understanding such intrinsic vs. thermal effects on oxide vs. halide perovskites holds the potential for designing target electronic properties.

36 MATERIALS SCIENCE↗

The Dynamical Role of the Chesapeake Bay on the Local Ozone Pollution Using Mesoscale Modeling—A Case Study

This study investigated the dynamic influence of the Chesapeake Bay (CB) on local ozone (O3) concentration and distribution using a weather forecasting model. The Weather Research and Forecasting model coupled with Chemistry (WRF–Chem) was employed to simulate O3 production and transportation near the CB. Baseline (water) as well as sensitivity (nowater) model experiments of bay circulation were carried out with and without bay water by changing the water surface from water to land (loam). First, the model performance simulating O3 was evaluated using the baseline experiment results and AirNow surface wind and O3 observations. The results showed that the model overestimates surface O3 by up to 20–30%. Further, the comparisons of the baseline and sensitivity experiments revealed higher O3 mixing ratios, primarily due to the resulting bay breeze circulation. These increases, after considering model overestimation, represent a mean bay dynamics circulation-induced contribution of up to 10% at night and 5% during the day. Furthermore, the boundary layer over northern CB, where it is at its narrowest width, was higher (by 1.2 km on average) during daytime due to higher surface temperatures observed. The boundary layer depth difference between the northern, central, and southern regions of the bay leads to a differential in the role of bay circulation dynamics in the observed O3 increase. The relatively wider swath of water surface over southern CB resulted in a lower boundary layer depth and stronger breeze circulation and this circulation contributed to O3 concentrations. Moreover, since the case selected has a minimal bay breeze circulation, the associated surface ozone enhancements represent what is expected at least at a minimum.

Chesapeak Bay↗

HERO CarbonSAFE Phase 2 Project in the Columbia River Basalt Group

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Basalt CO2 storage has several advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. Associated decreases in pressure and CO2 migration risks, 3. Reduced long-term monitoring requirements with respect to plume tracking, 4. Widespread geographic distribution and, 5. Large storage potential due to thickness, porosity, and CO2 interactions with basalt. And for locations such as the Pacific Northwest, Hawaii, Iceland, India and Japan, basalts may offer the only economically feasible option for local CO2 storage. However, there are limited field-scale assessments of CO2 storage in basalt, and current carbon capture utilization and storage (CCUS) permitting and regulatory frameworks were developed for conventional saline reservoirs. HERO CarbonSAFE is designed to address research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts. The HERO CarbonSAFE project will assess feasibility of developing a commercial-scale (50+ million metric tons of CO2) geological storage complex within the Columbia River Basalt Group (CRBG), a layered continental flood basalt complex that underlies Calpine’s natural gas-fired Hermiston Power Project (HPP) in Hermiston, OR (Figure 1). Under this 2-year CarbonSAFE Phase II project, the HERO team will conduct a data acquisition campaign that includes drilling a stratigraphic well to a total depth of ~1,500 m into the thick layered basalts proximal to HPP. A comprehensive well logging and hydrologic testing program will be augmented with new core collected from flow zones and sealing units, and comprehensive laboratory testing to help refine the kinetic rates of mineralization. The newly acquired information will be integrated with existing data from regional wells to correlate basalt injection zone properties to develop storage hub/commercial-scale models. Using these models, the project team will evaluate injection scenarios to define the technical and economic potential for storing a minimum of 50 million metric tons of CO2 over a 30-year period, along with a robust sensitivity analysis on key parameters governing reservoir viability for sustainable injection over a commercial project lifetime. Specific technical objectives of HERO are: (1) assessing the reservoir response of a series of stacked layered reservoir flowtop sequences occurring in this area of the CRBG to commercial-scale injection volumes; (2) extending prior efforts by the project team to characterize the deep layered basalts encountered in regional studies, to leverage prior investments by U.S. Department of Energy’s (DOE) Carbon Storage program; (3) leveraging DOE’s mineralization characterization efforts to advance model parametrization for commercial scale injection of CO2 in basalts; (4) conducting risk assessments associated with scaling up to commercial storage hub injection goals, while validating DOE’s National Risk Assessment Partnership (NRAP) tools, to identify potential constraints that would prevent the CRBG from serving as a commercial-scale storage complex; (5) developing mitigation plans to address identified risks; (6) developing a commercial-scale injection and monitoring, verification and accounting (MVA) strategy; (7) utilizing computational models to define and minimize, if possible, the Area of Review (AoR) under Class VI regulations; and (8) developing a robust CO2 management strategy for CRBG that also considers a regional source/sink approach that is responsive to stakeholder needs and industrial demand. Specific institutional objectives are: (1) identifying and developing plans to mitigate the nontechnical challenges associated with the build-out of a commercial-scale storage complex within the CRBG with integrated CO2 sources; (2) implementing the community outreach plan; (3) conducting regulatory research, including a survey of issues related to pore space ownership, MVA and long-term assurance of mineralization-based storage, to support an eventual application for a UIC Class VI permit; (4) advancing the project’s plan for CO2 liability management; and (5) continuing to refine and update the project’s economic model. The final objective is the preparation of a comprehensive Site Characterization Plan that draws upon the technical and institutional feasibility assessments to prepare the project for future commercialization efforts.

58 GEOSCIENCES↗

The Geoscience Laser Altimeter System (GLAS) for the ICESAT Mission

The Laser In space Technology Experiment, Shuttle Laser Altimeter and the Mars Observer Laser Altimeter have demonstrated accurate measurements of atmospheric backscatter and Surface heights from space. The recent MOLA measurements of the Mars surface have 40 cm vertical resolution and have reduced the global uncertainty in Mars topography from a few km to about 5 m. The Geoscience Laser Altimeter System (GLAS) is a next generation lidar for Earth orbit being developed as part of NASA's Icesat Mission. The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, establish a grid of accurate height profiles of the Earth's land topography, and profile the vertical backscatter of clouds and aerosols on a global scale. GLAS is being developed to fly on a small dedicated spacecraft in a polar orbit with a 590 630 km altitude at inclination of 94 degrees. GLAS is scheduled to launch in the summer 2001 and to operate continuously for a minimum of 3 years with a goal of 5 years. The primary mission for GLAS is to measure the seasonal and annual changes in the heights of the Greenland and Antarctic ice sheets. GLAS will continuously measure the vertical distance from orbit to the Earth's surface with 1064 nm pulses from a ND:YAG laser at a 40 Hz rate. Each 5 nsec wide laser pulse is used to produce a single range measurement, and the laser spots have 66 m diameter and about 170 m center-center spacings. When over land GLAS will profile the heights of the topography and vegetation. The GLAS receiver uses a 1 m diameter telescope and a Si APD detector. The detector signal is sampled by an all digital receiver which records each surface echo waveform with I nsec resolution and a stored echo record lengths of either 200, 400, or 600 samples. Analysis of the echo waveforms within the instrument permits discrimination between cloud and surface echoes. Ground based echo analysis permits precise ranging, determining the roughness or slopes of the surface as well as the vertical distributions of vegetation illuminated by the laser. Accurate knowledge of the laser beam's pointing angle is needed to prevent height biases when over sloped surfaces. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam's centroid to about 8 urad is needed to achieve 10 cm height accuracy. GLAS uses a stellar reference system (SRS) to determine the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith and a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser is measured pulse relative to the star camera with a laser reference system (LRS). Optically measuring each laser far field pattern relative to the orientation of the star camera and gyroscope permits the precise pointing angle of each laser pulse to be determined. GLAS will also determine the vertical distributions of clouds and aerosols by measuring the vertical profile of laser energy backscattered by the atmosphere at both 1064 and 532 nm. The 1064 nm measurements use the Si APD detector and profile the height and vertical structure of thicker clouds. The measurements at 532 nm use new highly sensitive photon counting, detectors, and measure the height distributions of very thin Clouds and aerosol layers. With averaging these can be used to determine the height of the planetary boundary layer. The instrument design and expected performance will be discussed.

Abshire, James B.↗