Search NASA⌕ Search

SEARCH · Search NASA

Results for “continuous time system estimation”

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 631 records · Page 35

NEXT Long-Duration Test Plume and Wear Characteristics after 16,550 h of Operation and 337 kg of Xenon Processed

The NASA s Evolutionary Xenon Thruster (NEXT) program is developing the next-generation ion propulsion system with significant enhancements beyond the state-of-the-art. The NEXT ion propulsion system provides improved mission capabilities for future NASA science missions to enhance and enable Discovery, New Frontiers, and Flagship-type NASA missions. As part of a comprehensive thruster service life assessment utilizing both testing and analyses, a Long-Duration Test (LDT) was initiated to validate and qualify the NEXT propellant throughput capability to a qualification-level of 450 kg, 1.5 times the mission-derived throughput requirement of 300 kg. This wear test is being conducted with a modified, flight-representative NEXT engineering model ion thruster, designated EM3. As of June 25, 2008, the thruster has accumulated 16,550 h of operation: the first 13,042 h at the thruster full-input-power of 6.9 kW with 3.52 A beam current and 1800 V beam power supply voltage. Operation since 13,042 h, i.e., the most recent 3,508 h, has been at an input power of 4.7 kW with 3.52 A beam current and 1180 V beam power supply voltage. The thruster has processed 337 kg of xenon (Xe) surpassing the NSTAR propellant throughput demonstrated during the extended life testing of the Deep Space 1 flight spare. The NEXT LDT has demonstrated a total impulse of 13.3 106 N s; the highest total impulse ever demonstrated by an ion thruster. Thruster plume diagnostics and erosion measurements are obtained periodically over the entire NEXT throttle table with input power ranging 0.5 to 6.9 kW. Observed thruster component erosion rates are consistent with predictions and the thruster service life assessment. There have not been any observed anomalous erosion and all erosion estimates indicate a thruster throughput capability that exceeds ~750 kg of Xe, an equivalent of 36,500 h of continuous operation at the full-power operating condition. This paper presents the erosion measurements and plume diagnostic results for the NEXT LDT to date with emphasis on the change in thruster operating condition and resulting impact on wear characteristics. Ion optics grid-gap data, both cold and operating, are presented. Performance and wear predictions for the LDT throttle profile are presented.

Herman, Daniel A.↗

The NASA Soil Moisture Active Passive (SMAP) Mission - Algorithm and Cal/Val Activities and Synergies with SMOS and Other L-Band Missions

NASA's Soil Moisture Active Passive (SMAP) mission, planned for launch in late 2014, has as its key measurement objective the frequent, global mapping of near-surface soil moisture and its freeze-thaw state. SMAP soil moisture and freeze/thaw measurements at 10 km and 3 km resolutions respectively, would enable significantly improved estimates of water, energy and carbon transfers between the land and atmosphere. Soil moisture control of these fluxes is a key factor in the performance of atmospheric models used for weather forecasts and climate projections Soil moisture measurements are also of great importance in assessing floods and for monitoring drought. In addition, observations of soil moisture and freeze/thaw timing over the boreal latitudes can help reduce uncertainties in quantifying the global carbon balance. The SMAP measurement concept utilizes an L-band radar and radiometer sharing a rotating 6-meter mesh reflector antenna. The SMAP radiometer and radar flight hardware and ground processing designs are incorporating approaches to identify and mitigate potential terrestrial radio frequency interference (RFI). The radar and radiometer instruments are planned to operate in a 680 km polar orbit, viewing the surface at a constant 40-degree incidence angle with a 1000-km swath width, providing 3-day global coverage. Data from the instruments would yield global maps of soil moisture and freeze/thaw state to be provided at 10 km and 3 km resolutions respectively, every two to three days. Plans are to provide also a radiometer-only soil moisture product at 40-km spatial resolution. This product and the underlying brightness temperatures have characteristics similar to those provided by the Soil Moisture and Ocean Salinity (SMOS) mission. As a result, there are unique opportunities for common data product development and continuity between the two missions. SMAP also has commonalities with other satellite missions having L-band radiometer and/or radar sensors applicable to soil moisture measurement, such as Aquarius, SAO COM, and ALOS-2. The algorithms and data products for SMAP are being developed in the SMAP Science Data System (SDS) Testbed. The algorithms are developed and evaluated in the SDS Testbed using simulated SMAP observations as well as observational data from current airborne and spaceborne L-band sensors including SMOS. The SMAP project is developing a Calibration and Validation (Cal/Val) Plan that is designed to support algorithm development (pre-launch) and data product validation (post-launch). A key component of the Cal/Val Plan is the identification, characterization, and instrumentation of sites that can be used to calibrate and validate the sensor data (Level I) and derived geophysical products (Level 2 and higher). In this presentation we report on the development status of the SMAP data product algorithms, and the planning and implementation of the SMAP Cal/Val program. Several components of the SMAP algorithm development and Cal/Val plans have commonality with those of SMOS, and for this reason there are shared activities and resources that can be utilized between the missions, including in situ networks, ancillary data sets, and long-term monitoring sites.

Njoku, Eni↗

OP09O-OP404-9 Wide Field Camera 3 CCD Quantum Efficiency Hysteresis

The HST/Wide Field Camera (WFC) 3 UV/visible channel CCD detectors have exhibited an unanticipated quantum efficiency hysteresis (QEH) behavior. At the nominal operating temperature of -83C, the QEH feature contrast was typically 0.1-0.2% or less. The behavior was replicated using flight spare detectors. A visible light flat-field (540nm) with a several times full-well signal level can pin the detectors at both optical (600nm) and near-UV (230nm) wavelengths, suppressing the QEH behavior. We are characterizing the timescale for the detectors to become unpinned and developing a protocol for flashing the WFC3 CCDs with the instrument's internal calibration system in flight. The HST/Wide Field Camera 3 UV/visible channel CCD detectors have exhibited an unanticipated quantum efficiency hysteresis (QEH) behavior. The first observed manifestation of QEH was the presence in a small percentage of flat-field images of a bowtie-shaped contrast that spanned the width of each chip. At the nominal operating temperature of -83C, the contrast observed for this feature was typically 0.1-0.2% or less, though at warmer temperatures contrasts up to 5% (at -50C) have been observed. The bowtie morphology was replicated using flight spare detectors in tests at the GSFC Detector Characterization Laboratory by power cycling the detector while cold. Continued investigation revealed that a clearly-related global QE suppression at the approximately 5% level can be produced by cooling the detector in the dark; subsequent flat-field exposures at a constant illumination show asymptotically increasing response. This QE "pinning" can be achieved with a single high signal flat-field or a series of lower signal flats; a visible light (500-580nm) flat-field with a signal level of several hundred thousand electrons per pixel is sufficient for QE pinning at both optical (600nm) and near-UV (230nm) wavelengths. We are characterizing the timescale for the detectors to become unpinned and developing a protocol for flashing the WFC3 CCDs with the instrument's internal calibration system in flight. A preliminary estimate of the decay timescale for one detector is that a drop of 0.1-0.2% occurs over a ten day period, indicating that relatively infrequent cal lamp exposures can mitigate the behavior to extremely low levels.

Collins, Nick↗

Modeling and Analysis of Coupled Flow Through Packed Beds as a Diagnostic Tool for Packed Bed Reactor Health - 17172

The US Department of Energy, Office of River Protection, is building the Waste Treatment and Immobilization Plant (WTP) in Richland, Washington, to immobilize the hundreds of millions of liters of radioactive waste stored at the Hanford Site in Washington state. The offgas generated during vitrification of these wastes requires treatment to remove and/or destroy chemical and radioactive hazardous components before release to the atmosphere. In particular, volatile and semi-volatile hydrocarbon components are to be catalytically oxidized to carbon dioxide and water, whereas nitrogen oxides (NOx) are to be catalytically reduced to nitrogen and water. These tasks are accomplished in a series of packed bed reactors. Achieving the requisite Destruction/Removal Efficiency (DRE) for each of the tracked components is highly dependent on the flow characteristics of the offgas through the packed beds. The flow kinetics may be monitored and characterized during the initial warm-up of the system. The gases, initially at 76 deg. C (169 deg. F), (a) flow through one side of a countercurrent heat exchanger to recover heat from the exiting gases (which have already passed through the catalytic beds), (b) are then further heated electrically to approximately 400 deg. C (750 deg. F; the preferred operating temperature of the Thermo Catalytic Oxidizer (TCO)), (c) pass through the TCO where the reduced components of the offgas are oxidized over a supported Pt catalyst bed, (d) are partially cooled to approximately 315 deg. C (600 deg. F) by addition of approximately 25% (by weight) diluted ammonia/air mixture, (e) pass through the Selective Catalytic Reducer (SCR) where the NOx components are converted to N{sub 2} and H{sub 2}O over a supported V2O{sub 5}/TiO{sub 2} catalyst bed, and (f) exit through the hot side of the countercurrent heat exchanger. The electrical heater is not powerful enough to raise the temperature of the offgas to the desired 400 deg. C without the contribution of the heat recovered in the heat exchanger. Moreover, the thermal mass of the beds and housings of the TCO and SCR require a significant input of thermal energy to reach the desired operating temperatures; the better part of a day is required to bring the unit from ambient temperatures to operating temperatures. Therefore, the temperature of the stream during the initial warm-up of the system at each point in the system is dependent upon the thermal mass, flow dispersion, and permeability of the components, quantified as characteristic/representative volumes and residence times. By monitoring the temperature of the offgas as it enters and exits each portion of the system during this warm-up period one may model the system as a series of small Continuously Stirred Tank Reactors (CSTR)s (referred to as 'CSTRlets') to explain the flow kinetics through the catalyst beds. These component CSTRlets are adjustable in size and number, that one may model the flow profile of the gases through the catalyst beds. For example, a single large CSTR will present an effluent temperature profile with a simple exponential decay shape. Concatenating more CSTRlets generate more complex sigmoidal responses. Concatenating one CSTRlet significantly larger than the others in the train also alters the shape of the effluent temperature profile, in a fashion not replicable by a train of equal-sized CSTRlets. The ultimate goal is to accurately emulate the entire system as a representative CSTRlet train. Once accomplished, this characterization may be used in subsequent warm-ups (following shutdown for catalyst maintenance, for example) or by an intentional step change in inlet offgas temperature to detect/diagnose anomalies in flow through the beds (e.g., channeling or partial blockage.) This work presents development of this model through analytical, calculational, and empirical means. - The analytic model is accomplished by Laplace transform analysis of each component to recover component process transfer functions. - The calculational portion is performed via brute force finite element calculation. - Fitting empirical data from the warm-up of the TCO skid allows one to select the hypothetical CSTRlet trains which best fit the TCO and SCR beds for the analytical and finite element calculations. The first portion of this paper presents the development of the analytical (Laplace space) model of the coupled system. Estimated values for the initial temperature, mass, specific heat, electrical heat, and relative quantities of the offgas, catalyst beds, stainless steel housings, and ammonia dilution stream are used to generate expected temperature profiles at several selected junctures of the system. The number and size of the component hypothetical CSTRlets are calculated from the analytical model and are used to predict temperature profiles at the exit of the system and at intermediate junctions between system components (where thermocouples are present in the actual unit.) The finite element calculations employ a series of CSTRlet volumes to produce a similar suite of temperature profiles, for comparison with the analytical model results. Finally, empirical data from the Factory Acceptance Test system warm-up are used to select and compare the results from each modeling strategy to select the CSTRlet train(s) which best model the behavior of the actual system. This work presents a novel method of system component characterization, of value to the WTP plant and similar vitrification facilities as well as chemical process industries unrelated to waste management. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A computational modeling framework for pre-clinical evaluation of cardiac mapping systems

There are a variety of difficulties in evaluating clinical cardiac mapping systems, most notably the inability to record the transmembrane potential throughout the entire heart during patient procedures which prevents the comparison to a relevant “gold standard”. Cardiac mapping systems are comprised of hardware and software elements including sophisticated mathematical algorithms, both of which continue to undergo rapid innovation. The purpose of this study is to develop a computational modeling framework to evaluate the performance of cardiac mapping systems. The framework enables rigorous evaluation of a mapping system’s ability to localize and characterize (i.e., focal or reentrant) arrhythmogenic sources in the heart. The main component of our tool is a library of computer simulations of various dynamic patterns throughout the entire heart in which the type and location of the arrhythmogenic sources are known. Our framework allows for performance evaluation for various electrode configurations, heart geometries, arrhythmias, and electrogram noise levels and involves blind comparison of mapping systems against a “silver standard” comprised of computer simulations in which the precise transmembrane potential patterns throughout the heart are known. A feasibility study was performed using simulations of patterns in the human left atria and three hypothetical virtual catheter electrode arrays. Activation times (AcT) and patterns (AcP) were computed for three virtual electrode arrays: two basket arrays with good and poor contact and one high-resolution grid with uniform spacing. The average root mean squared difference of AcTs of electrograms and those of the nearest endocardial action potential was less than 1 ms and therefore appears to be a poor performance metric. In an effort to standardize performance evaluation of mapping systems a novel performance metric is introduced based on the number of AcPs identified correctly and those considered spurious as well as misclassifications of arrhythmia type; spatial and temporal localization accuracy of correctly identified patterns was also quantified. This approach provides a rigorous quantitative analysis of cardiac mapping system performance. Proof of concept of this computational evaluation framework suggests that it could help safeguard that mapping systems perform as expected as well as provide estimates of system accuracy.

59 BASIC BIOLOGICAL SCIENCES↗

Characterization of the Triassic Newark Basin of New York and New Jersey for geologic storage of carbon dioxide

The Newark Basin is a Triassic-aged rift basin underlying densely populated, industrialized sections of New York, New Jersey and Pennsylvania. The Basin is an elongate half-graben encompassing an area of more than 7,510 square-kilometers (2,900 square-miles), and could represent a key storage component for commercial scale management of carbon dioxide emissions via geologic sequestration. The project team first acquired published reports, surface and subsurface maps, and seismic data, which formed the basis for a three-dimensional model framework for the northern end of the Basin incorporating stratigraphic, hydrologic, and water quality data. Field investigations included drilling, coring, and logging of two stratigraphic test borings in Clarkstown, NY (Exit 14 Tandem Lot Well No. 1), drilled to a depth of 2,099 meters (6,885 feet); and Palisades, NY (Lamont-Doherty Earth Observatory Test Well No. 4) drilled to a depth of 549 meters (1,802 feet). Two two-dimensional seismic reflection data lines arrayed perpendicularly were acquired by Schlumberger/WesternGeco to help characterize the structure and stratigraphy and as part of pre-drilling field screening activities for the deep stratigraphic borehole. A total of 47 meters (155 feet) of continuous whole core was recovered from the Tandem Lot boring from depths of 1,393 meters (4,570 feet) to 1,486 meters (4,877 feet). Twenty-five horizontal rotary cores were collected in mudstones and sandstones in the surface casing hole and fifty-two cores were taken in various lithologies in the deep borehole. Rotary core plugs were analyzed by Weatherford Laboratories for routine and advanced testing. Rotary core plug trim end thin sections were evaluated by the New York State Museum for mineralogical analysis and porosity estimation. Using core samples, Lawrence Berkley National Laboratory designed and completed laboratory experiments and numerical modeling analyses to characterize the dissolution and reaction of carbon dioxide with formation brine and minerals, and resulting effects on injection rate, pressure, effective storage volume, and carbon dioxide migration within a prospective sandstone reservoir. $$Three potential porous and permeable sandstone units were identified in the Passaic Formation at the New York State Thruway Exit 14 location. Potential Flow Unit 1, at a depth of 643 meters (2,110 feet) to 751 meters (2,465 feet); Potential Flow Unit 2 at a depth of 853 meters (2,798 feet) to 1,000 meters (3,280 feet); and Potential Flow Unit 3, at a depth of 1,114 meters (3,655 feet) to 1,294 meters (4,250 feet). Reactive transport simulations of interactions between carbon dioxide, brine and formation minerals were carried out to evaluate changes in formation water chemistry, mineral precipitation and dissolution reactions, and any potential resulting effects on formation permeability. The experimental and modeling analyses suggest that mineral precipitation and dissolution reactions (within the target formation) are not expected to lead to significant changes to the underground hydrologic system over time frames (~30 years) typically relevant for carbon dioxide injection operations. Key findings of this basin characterization study include an estimate of carbon dioxide storage capacity in the Newark Basin. Assuming an average porosity of twelve percent and an aquifer volume of 6.1E+12 meters3, calculated ranges of likely storage capacity range from 1.9 – 20.2 gigatonnes under high temperature (low carbon dioxide density) conditions; and 2.9 – 30.2 gigatonnes under low temperature (low carbon dioxide density) conditions. Intra-basin faulting, geometry of the Palisades Sill, and the presence of altered meta-sediments above and below the Sill, increase potential compartmentalization within the basin. A structural/stratigraphic trap type may occur where porous/permeable sediments are cross-cut by the Palisades Sill. Potential injection intervals are present within the Stockton Formation of the Newark Basin. Additional porous/permeable intervals may be present within sandstones of the Passaic Formation, increasing projected storage capacity. Deeper wedges of strata are likely present in the deeper portions of the basin in southern New York and into northern New Jersey. Abundant mudstones are present in the Passaic, Lockatong, and Stockton Formations. These intervals have the requisite petrophysical properties to form effective primary and secondary containment intervals to industrial-scale sequestration of carbon dioxide in the Newark Basin. Hydro-thermally altered meta-sediments in the region immediately surrounding the top and base of the Palisades Sill is devoid of porosity/permeability and forms an additional effective lateral/vertical sealing cap rock.

58 GEOSCIENCES↗

Summary Describing Integration of ERM Methodology into Supervisory Control Framework with Software Package Documentation; Advanced Reactor Technology Milestone: M4AT-16PN2301052

This project involved the development of enhanced risk monitors (ERMs) for active components in Advanced Reactor (AdvRx) designs by integrating real-time information about equipment condition with risk monitors. Health monitoring techniques in combination with predictive estimates of component failure based on condition and risk monitors can serve to indicate the risk posed by continued operation in the presence of detected degradation. This combination of predictive health monitoring based on equipment condition assessment and risk monitors can also enable optimization of maintenance scheduling with respect to the economics of plant operation. This report summarizes PNNL’s multi-year project on the development and evaluation of an ERM concept for active components while highlighting FY2016 accomplishments. Specifically, this report provides a status summary of the integration and demonstration of the prototypic ERM framework with the plant supervisory control algorithms being developed at Oak Ridge National Laboratory (ORNL), and describes additional case studies conducted to assess sensitivity of the technology to different quantities. Supporting documentation on the software package to be provided to ONRL is incorporated in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation of Small Satellite Dynamics Simulation Modules using ASTERIA Flight Data

ASTERIA (Arcsecond Space Telescope Enabling Research in Astrophysics) was a CubeSat space telescope that operated in low-Earth orbit, having been deployed from the International Space Station in 2017. The spacecraft has achieved sub-arcsecond pointing stability and millikelvin thermal stability over 20-minute observations. A key enabling feature of the ASTERIA mission is the low level of pointing error that has been demonstrated under ASTERIA’s fine pointing control mode. Prior to launch, the ASTERIA mission performed analysis and simulations to estimate the in-flight pointing performance. This analysis used reaction wheel models provided by JPL’s Small Satellite Dynamics Testbed (SSDT), which had performed Kistler table testing to characterize the jitter caused by the reaction wheels of the Blue Canyon Technologies XACT attitude control unit within ASTERIA. These models were subsequently incorporated into the SSDT simulation. The main source of jitter, the high-frequency attitude disturbance over the camera exposure time, on ASTERIA is the set of rotating reaction wheels that are used for maintaining fine pointing during observations. The reaction wheels impart disturbance forces and torques continuously that cause unwanted motion during imaging. Increased jitter consequently leads to blurring of the image. Therefore, to benefit all SmallSat missions devoted to photometric or spectroscopic astrophysics applications, an on-orbit flight data acquisition and jitter testing campaign was performed in an attempt to help validate the SSDT’s simulation models. This paper describes the process used to validate several simulation models and attempt to characterize the ASTERIA jitter by analyzing the size of the resulting spot size for the target stars as a function of wheel speed, having commanded four wheel speeds for each of three stars of known brightness. After summarizing the on-orbit observation and jitter level measurement process, this paper compares the results obtained from the flight operational environment to the corresponding set of simulation jitter levels. Though the flight data obtained during this experiment was not sufficient to validate the simulation’s jitter model, other models have been validated with flight data, including the magnetic field, orbit propagation, sun position, and long-duration orbit decay models, which may be used by all other projects that use the SSDT’s simulation to increase system capability knowledge.

Mohan, Swati↗

Climate-simulated raceway pond culturing: quantifying the maximum achievable annual biomass productivity of Chlorella sorokiniana in the contiguous USA

Chlorella sorokiniana (DOE 1412) emerged as one of the most promising microalgae strains from the NAABB consortium project, with a remarkable doubling time under optimal conditions of 2.57 hr-1. However, its maximum achievable annual biomass productivity in outdoor ponds in the contiguous United States remained unknown. In order to address this knowledge gap, this alga was cultured in indoor LED-lighted and temperature-controlled raceways in nutrient replete freshwater (BG-11) medium at pH 7 under conditions simulating the daily sunlight intensity and water temperature fluctuations during three seasons in Southern Florida, an optimal outdoor pond culture location for this organism identified by biomass growth modeling. Prior strain characterization indicated that the average maximum specific growth rate (µmax) at 36 ºC declined continuously with pH, with µmax corresponding to 5.92, 5.83, 4.89, and 4.21 day-1 at pH 6, 7, 8, and 9, respectively. In addition, the maximum specific growth rate declined nearly linearly with increasing salinity until no growth was observed above 35 g/L NaCl. In the climate-simulated culturing studies, the volumetric ash-free dry weight-based biomass productivities during the linear growth phase were 57, 69, and 97 mg/L-day for 30-year average light and temperature simulations for January (winter), March (spring), and July (summer), respectively, which corresponds to average areal productivities of 11.6, 14.1, and 19.9 g/m2-day at a constant pond depth of 20.5 cm. The photosynthetic efficiencies (PAR) in the three climate-simulated pond culturing experiments ranged from 4.1 to 5.1%. The annual biomass productivity was estimated as ca. 15 g/m2-day, nearly double the U.S. Department of Energy (DOE) 2015 State of Technology annual cultivation productivity of 8.5 g/m2-day, but this is still significantly below the projected 2022 target of ca. 25 g/m2-day (U.S. DOE, 2016) for economic microalgal biofuel production, indicating the need for additional research in strain biology and system engineering.

annual biomass productivity↗

A Multiple Parameters Biodosimetry Tool with Various Blood Cell Counts - the Hemodose Approach

There continue to be important concerns about the possibility of the occurrence of acute radiation syndromes following nuclear and radiological terrorism or accidents that may result in mass casualties in densely populated areas. To guide medical personnel in their clinical decisions for effective medical management and treatment of the exposed individuals, biological markers are usually applied to examine radiation induced biological changes to assess the severity of radiation injury to sensitive organ systems. Among these the peripheral blood cell counts are widely used to assess the extent of radiation induced bone marrow injury. This is due to the fact that the hematopoietic system is the most vulnerable part of the human body to radiation damage. Particularly, the lymphocyte, granulocyte, and platelet cells are the most radiosensitive of the blood elements, and monitoring their changes after exposure is regarded as a practical and recommended laboratory test to estimate radiation dose and injury. Based upon years of physiological and pathophysiological investigation of mammalian hematopoietic systems, and rigorous coarse-grained bio-mathematical modeling and validation on species from mouse, to dog, monkey, and human, we have developed a set of software tools Hemodose, which can use single or serial granulocyte, lymphocyte, leukocyte, or platelet counts after exposure to estimate absorbed doses of adult victims very rapidly and accurately. Some patient data from historical accidents are utilized as examples to demonstrate the capabilities of these tools as a rapid point-of-care diagnostic or centralized high-throughput assay system in a large-scale radiological disaster scenario. Most significant to the improvement of national and local preparedness of a potential nuclear/radiological disaster, this HemoDose approach establishes robust correlations between the absorbed doses and victim's various types of blood cell counts not only in the early time window (1 or 2 days), but also in the very late phase (up to 4 weeks) after exposure.

Hu, Shaowen↗

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗

Formation of Planetary Systems

Stars are observed to be forming within cold regions of the galaxy called molecular clouds. A clump of gas and dust within a molecular cloud can collapse into a rotationally- supported disk orbiting the pre s sure -supported star. Such a disk has the same initial elemental composition as the growing star. At sufficient distances from the central star, it is cool enough for approx. 1 - 2% of this material to be in solid form, either remnant interstellar grains or condensates formed within the disk. During the infall stage, the disk is very active and probably highly turbulent. When the infall slows substantially or stops, the disk becomes more quiescent. The dust grows from micron-sized dust to kilometer-sized planetesimals by physical collisions, possibly aided by collective gravity. The dynamics of larger solid bodies within protoplanetary disks are better characterized. The primary perturbations on the orbits of kilometer-sized and larger planetesimals in protoplanetary disks are mutual gravitational interactions and physical collisions. These interactions lead to accretion (and in some cases erosion and fragmentation) of planetesimals. The most massive planets have the largest gravitationally-enhanced collision cross-sections, and accrete almost everything with which they collide. The size distribution of solid bodies becomes quite skewed, with a few large bodies growing much faster than the rest of the swarm, until they have accumulated most of the small bodies within their gravitational reach. Slower growth continues (at least for solid, earth-like planets) as the eccentricities of planetary embryos are pumped up by long-range mutual gravitational perturbations. As planetary masses increase, they become more efficient at stirring random velocities of neighboring bodies. If sufficiently massive and dense planets exist far enough from the star, they can eject material into interstellar space. In most models, giant planets begin their growth like terrestrial planets, but if they become massive enough before the protoplanetary disk dissipates, then they are able to accumulate substantial amounts of gas. Giant planet growth times predicted by current models are similar to estimates of the lifetime of the gaseous protoplanetary disk, so we don't know whether or not giant planets form in most protoplanetary disks. However, planets that become massive while a substantial amount of gas remains in the disk may migrate into the star as a consequence of their gravitational interactions with the disk.

Lissauer, Jack J.↗

Cleaning at the Edge of Science: NASA's Genesis Mission

As part of NASA's continuing exploration of the origins of our solar system, the California Institute of Technology, Jet Propulsion Laboratory, Lockheed Martin Astronautics, Los Alamos National Laboratory, and the Johnson Space Center are working together to develop the Genesis mission to return solar matter for analysis in terrestrial laboratories. These samples will be used to define a baseline for the chemical and isotopic composition of the solar nebula. Deviations from the baseline resulted as the solar system evolved; thus, providing a tracer for materials incorporated into meteorites, comets and planetary bodies. These differences represent "fossil residues" that provide invaluable insight into how the solar nebula evolved to form the planets. We cannot collect a sample of the Sun as we would for a planet; fortunately, solar material comes to us in the form of the solar wind. Ultrapure materials will be exposed at the Earth-Sun L1, outside the Earth's magnetic influence, where solar wind nuclei will be captured for 2 years before returning to Earth in January 2001. The key challenge to obtaining a good sample of solar wind, uncontaminated by terrestrial atoms, is a clean collection surface in a clean sample canister and clean facilities with which to handle the samples for allocation and future reference. The Johnson Space Center QSQ is responsible for contamination control for the mission, for ensuring the cleanliness of collection surfaces and providing a clean environment for their subsequent handling. The level of cleanliness required is high; at the time of analysis (after sample return), the surface contamination by C, N, O must each be less than 10(exp 15) atoms per centimeter squared and for elements other than C, N, O, the number of atoms per centimeter squared of each surface contaminant shall not exceed the estimated solar wind fluence of the species (varies by element between U at approx. 10 (exp 4) atoms per centimeter squared to Fe, Si, Mg, and Ne at approx. 10(exp 12), atoms per centimeter squared). Typical spacecraft assembly is done in class 10,000 cleanrooms. The final cleaning and reintegration of the Genesis payload canister as well as all sample material handling will be done within a class 10 cleanroom using Dryden suits to protect the collector materials from any human debris. Each component is unique, no standard size, shape, material, or precleaning history. We are developing new final cleaning techniques utilizing ultra-pure water to minimize molecular residues on the hardware components.

Stansbery, Eileen K.↗

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗

Measuring the neutrino-oxygen neutral current quasielastic cross section using the accelerator neutrino neutron interaction experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located on-axis to Fermilab’s Booster Neutrino Beam (BNB). ANNIE is uniquely positioned to perform high-statistics measurements of neutrino-nucleus interactions in water, benefiting from a large neutrino flux due to a short (100-meter) baseline. A central focus of ANNIE’s physics program is the measurement of both charged current (CC) and neutral current (NC) cross sections on water, including neutral current quasielastic (NCQE) and CC-inclusive channels. The NCQE measurement is particularly critical for constraining uncertainties in rare-event searches such as the Diffuse Supernova Neutrino Background (DSNB), where atmospheric $\nu$NCQE interactions constitute a significant and poorly constrained background. This dissertation presents a measurement of the flux-averaged neutrino-oxygen neutral current quasielastic ($\nu$NCQE) cross section using $2.573 \times 10^{20}$~POT of BNB exposure from the 2022 and 2023 beam years. The $\nu$NCQE interaction is identified through the primary $\gamma$-rays produced by nuclear de-excitation of the residual $^{15}$N$^*$ or $^{15}$O$^*$ nucleus following nucleon knockout from $^{16}$O. A dedicated Monte Carlo (MC) re-tuning campaign was conducted using an americium-beryllium (AmBe) calibration source, Michel electrons from stopped muons, and throughgoing dirt muons originating upstream of the detector. This multi-sample approach provided a wide-ranging $\mathcal{O}(\text{MeV})$--$\mathcal{O}(\text{GeV})$ dataset for tuning the simulated detector response, which was subsequently validated against AmBe neutron and Michel electron data for use in the $\nu$NCQE analysis. A dedicated laser calibration campaign was carried out to reduce timing uncertainties across the PMT system, enabling reconstruction of the BNB bunch substructure with sufficient resolution to serve as a background rejection tool. By selecting events in-time with individual neutrino bunches, beam-correlated $\nu$NCQE events are separated from diffuse and accelerator-induced backgrounds, notably skyshine neutrons and externally-originating events, that would otherwise dominate traditional charge-based selections within a small-scale, surface-level, short-baseline detector. A data-driven estimation of the skyshine neutron and external background rates was performed and incorporated into the systematic uncertainty budget. The flux-averaged $\nu$NCQE cross section on oxygen is measured to be $1.57 \pm 0.06\,(\text{stat.})$ $^{+0.91}_{-0.67}\,(\text{syst.})$ $\times 10^{-38}\ \text{cm}^{2}$. A full systematic budget is constructed by propagating uncertainties in the secondary hadronic interaction modeling, background cross section normalizations, detector response, neutrino flux, and the primary $\gamma$-ray emission probabilities from oxygen nuclear de-excitation. An idealized de-excitation model, constructed from existing measurements in the literature is developed to benchmark the predictions of the \textsc{GENIE} event generator. A comparison reveals that \textsc{GENIE} systematically overpredicts the primary $\gamma$-ray emission probability from oxygen de-excitation by a factor of $1.49\times$ for $E_\gamma > 6$~MeV and $3.07\times$ in the $3$--$6$~MeV band. This comparison motivates the dominant systematic uncertainty in this analysis, where a conservative uncertainty of $^{+39.9\%}_{-0\%}$ on the primary $\gamma$-ray signal prediction is assigned. The ANNIE result is consistent with and complementary to existing flux-averaged $\nu$NCQE cross section measurements from T2K and Super-Kamiokande, providing an independent measurement with a different detector, neutrino beam, and analysis methodology. Looking ahead, an upgrade to the ANNIE DAQ infrastructure enabling continuous extended readout will allow a complementary $\nu$NCQE neutron multiplicity measurement, directly relevant to constraining the NCQE background in DSNB searches, competitive with the recent T2K measurement at SK-Gd. The planned Super-SANDI upgrade, deploying a large Water-based Liquid Scintillator (WbLS) volume, will further extend ANNIE's reach to hadronic final states and exclusive NC channels, and enable joint measurements with liquid argon detectors sharing the BNB beamline ahead of DUNE and Hyper-Kamiokande.

Doran, Steven [Iowa State U.]↗

Utility-Scale Solar, 2024 Edition: Empirical Trends in Deployment, Technology, Cost, Performance, PPA Pricing, and Value in the United States [Slides]

Berkeley Lab’s “Utility-Scale Solar, 2024 Edition” presents analysis of empirical plant-level data from the U.S. fleet of ground-mounted photovoltaic (PV), PV+battery, and concentrating solar-thermal power (CSP) plants with capacities exceeding 5 MWAC (PV plants of 5 MWAC or less, including residential rooftop systems, are covered separately in Berkeley Lab’s companion annual report, Tracking the Sun). Key findings from this year’s report include: -18.5 GWAC of new utility-scale PV capacity came online in 2023, bringing cumulative installed capacity to more than 80.2 GWAC across 47 states. Installed costs continued to fall in 2023. Relative to 2022, capacity-weighted averages decreased by 8% to -$\$1.43$/WAC (or $\$1.08$/WDC). Costs, based on a 7.1 GWAC sample of 76 plants completed in 2023, have fallen by 75% (averaging 10% annually) since 2010. Plant-level capacity factors vary widely, from 6% to 36% (on an AC basis), with a sample median of 24%. -Levelized cost of energy (LCOE) of new 2023 projects increased slightly to $\$46$/MWh prior to the application of tax credits but continued to fall to $\$31$/MWh when accounting for federal incentives. PPA prices have largely followed the decline in solar’s LCOE over time, but newly signed longer-term PPA prices have increased since 2021, to an average of $\$35$/MWh (levelized, in 2023 dollars). -Solar’s average energy and capacity value (i.e., ability to offset costs of other power generation sources) across the U.S. was $\$45$/MWh in 2023. Solar’s average market value was lowest in CAISO ($\$27$/MWh), the market with the greatest solar generation share, and highest in ERCOT ($\$67$/MWh). -Newer solar projects had greater market value in 2023 than their generation costs, yielding $\$1.1$ billion in benefits. Projects built in 2022 delivered on average $\$15$/MWh more market value than their costs in 2023. -Solar’s combined value from wholesale electricity markets, public health and climate damage reduction were greater than generation costs and incentives, yielding $\$13.7$ billion in net benefits in 2023. We estimate U.S. health benefits of $\$24$/MWh and reduced global climate damages of $\$101$/MWh. -Adding battery storage is one way to increase the value of solar. Deployment of 52 new PV+battery hybrid plants set a record with 5.3 GW installed in 2023. Our public data file tracks metadata and PPA prices from more than 100 PV+battery hybrid projects that are already online or that have secured offtake arrangements. -Looking ahead, a massive pipeline of at least 1,085 GW of solar capacity dominates the nation’s interconnection queues at the end of 2023. Nearly 571 GW, or 53%, of that total was paired with a battery – in CAISO it was a staggering 98%. Historically only 10% of the requested solar capacity is built. -For more information, and to explore related interactive data visualizations, go to utilityscalesolar.lbl.gov.

14 SOLAR ENERGY↗

GPS Altimetry

The advent of satellite altimetry has greatly improved our ability to observe global ocean circulation. However, the swath of a single, nadir-viewing satellite altimeter is only a few km and the track spacing is several hundred km to resolve the two-dimensional structure of ocean eddies. Our goal is to increase spatial and temporal coverage by monitoring Global Positioning System (GPS) signals reflected from the ocean. A constellation of spacecraft would each carry a GPS receiver capable of recording 8 reflections simultaneously. The reflections are well distributed in azimuth and elevation and can be tracked continuously while the satellite is in view, and another is then acquired, as illustrated below. The diagram depicts a new approach at altimetry measurements where ocean surface reflected GPS signals are simultaneously tracked and processed in a GPS flight receiver in space. The reflected GPS signals from the ocean must be compared precisely with the direct GPS signals in order to infer the characteristics of the ocean from the combined data set. Understanding the features and accuracy of GPS altimetry measurement is crucial to establishing its suitability for oceanography. Preliminary work has enabled us to theoretically model the signal output of the correlator for a variety of system parameters such as wind speed (sea roughness), receiver height, incidence angle, receiver range and Doppler filter bandwidth and antenna gain. Expected signal-to-noise ratio has been estimated from which we have inferred, to a first approximation, the basic receiver gain requirements for a space-based altimeter and the expected range raw error. In 1998, work on a different task led to the extraction of the first reflected GPS signal observed from a spaceborne receiver during the 1995 Space Transportation System-68 (STS-68) Shuttle Radar Laboratory-2 (SRL-2) high resolution synthetic aperture radar mission. Good comparisons with our signal models have been obtained. Having established that only modest signal-to-noise ratios are obtainable unless very high gain antennas are used, we expect that single measurements might not provide an estimate of sea state parameters as accurate as that obtainable with traditional remote sensing instruments, if costs are limited. Therefore spatial and temporal averaging of many measurements is required. Since the receiving satellite tracks do not repeat, measurements in a given area will be used to refine the local solution as a function of time, and will define the spatial resolution. In order to make this instrument viable in space, we need to detect and process many scattered signals. In 1998 we made a first step at understanding what configurations of antenna gain and orientations capture the largest number of viable signals to be used in the subsequent spatial/temporal averaging process. Additional information is contained in the original.

Lichten, Stephen M.↗

Polarization performance simulation for the GeoXO atmospheric composition instrument: NO2 retrieval impacts

NOAA’s Geostationary Extended Observations (GeoXO) constellation will continue and expand on the capabilities of the current generation of geostationary satellite systems to support US weather, ocean, atmosphere, and climate operations. It is planned to consist of a dedicated atmospheric composition instrument (ACX) to support air quality forecasting and monitoring by providing similar capabilities to missions such as TEMPO (Tropospheric Emission: Monitoring Pollution), currently planned to launch in 2023, and Ozone Monitoring Instrument (OMI), TROPOMI (TROPOspheric Monitoring Instrument), and GEMS (Geostationary Environment Monitoring Spectrometer) currently in operation. As the early phases of ACX development are progressing, design trade-offs are being considered to understand the relationship between instrument design choices and trace gas retrieval impacts. Some of these choices will affect the instrument polarization sensitivity (PS), which can have radiometric impacts on environmental satellite observations. We conducted a study to investigate how such radiometric impacts can affect NO2 retrievals by exploring their sensitivities to time of day, location, and scene type with an ACX instrument model that incorporates PS. The study addresses the basic steps of operational NO2 retrievals: the spectral fitting step and the conversion of slant column to vertical column via the air mass factor (AMF). The spectral fitting step was performed by generating at-sensor radiance from a clear sky scene with a known NO2 amount, the application of an instrument model including both instrument PS and noise, and a physical retrieval. The spectral fitting step was found to mitigate the impacts of instrument PS. The AMF-related step was considered for clear sky and partially cloudy scenes, where instrument PS can lead to errors in interpreting the cloud content, propagating to AMF errors and finally to NO2 retrieval errors. For this step, the NO2 retrieval impacts were small but non-negligible for high NO2 amounts; we estimated that a typical high NO2 amount can cause maximum of 0.25 x 10^15 molecules/cm^2 for a PS of 5%. These simulation capabilities were designed to aid in the development of a GeoXO atmospheric composition instrument that will improve our ability to monitor and understand the Earth’s atmosphere.

trace gases↗