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 523 records · Page 29

Ultrasonic Measurements of Temperature Profile and Heat Fluxes in Coal-Fired Power Plants (Final Report)

Many industrial processes are inaccessible or inhospitable to characterization by traditional temperature measurement methods, such as thermocouples, especially over prolonged exposure to harsh environments. Ultrasound is an established characterization technology with diverse applications ranging from medical imaging to therapies to flaw detection to nondestructive evaluation. Ultrasound may characterize solid materials and components noninvasively as a nondestructive evaluation modality and obtain internal measurements of material properties. For example, the speed of ultrasound propagation changes with Young’s modulus and Poisson’s ratio, which can be found from its measurements. Traditional ultrasonic characterization assumes all material properties remain constant with the position. When this assumption holds, a property of interest may be measured by relating it to the speed of ultrasound propagation (or a speed of sound, SOS) and measuring the SOS by timing the ultrasound propagation through a known distance. However, when a property of interest is spatially distributed, the propagation time depends on the SOS changing with the position along the ultrasound propagation path. The multiple temperature distributions may lead to an identical time of flight (TOF). Temperature is one property that impacts the speed of ultrasound and often cannot be assumed to remain constant with the position. Previously, in the context of temperature, we addressed the challenge of ultrasonic characterization of spatially distributed properties by developing a method for measuring segmental temperature distributions (MSTD). This method divides the ultrasonic propagation into segments bound by echogenic features. These features provide ultrasonic interfaces where some energy is reflected toward the receiving transducer, and the rest continues through the medium. The time-of-flight between the echoes reflected from echogenic features characterizes the spatial distribution in the properties of interest in the corresponding segment of the ultrasonic propagation path. This project demonstrated the application of the MSTD method in industrial conditions of the coal-fired power plant. We implemented the MSTD using metals and alloys waveguides, which may be the existing structure for which the temperature distribution is characterized or purposefully designed waveguides added to the structure by welding or other means specifically to quantify thermal properties using the MSTD method. Previous iterations of the MSTD method used ceramic and cementitious waveguides, which significantly attenuate ultrasound. On the other hand, low attenuation in metallic waveguides creates interactions between echogenic features which compilates the signal analysis in the segmental TOF measurements. We have established the WG design principles that minimize the interferences between trailing and primary echoes and, in some cases, eliminate them. The waveguides in which echoes do not interfere improve the timing accuracy and the robustness of ultrasonic measurements of the spatial distributions in material properties. Our emphasis remained on the estimation of the temperature distributions. We have developed general recommendations for designing ultrasonically segmented waveguides with the reduced influence of trailing echoes. Two of our waveguide designs were tested in the industry. The first waveguide was designed for insertion into a combustion zone of the utility-scale coal-fired power plant boiler. The second design allows the characterization of temperature distribution in the direction normal to the boiler’s water wall, a large heat exchanger converting the chemical energy released during combustion to the steam driving the electrical power generation turbines. These waveguides were designed to operate within a restrictive space of thermally insulated water wall and incorporate densely located echogenic features while combatting the influence of trailing echoes. The project has successfully demonstrated the feasibility of using the developed method for accurate, continuous, and robust temperature measurements in extreme environments of power generation and other industrial processes. It, therefore, has achieved its overarching goal of advancing the technology readiness level of the novel Ultrasound Measurements of Segmental Temperature Distribution (US-MSTD) method for real-time measurements of the temperature distribution and heat fluxes closer to commercial availability, developing a prototype multipoint measurement system, and validating its performance on coal-fired utility boilers. The success of this project was achieved in collaboration with the power generator, Rocky Mountain Power, and set the stage for the transfer of this technology from the laboratory to the industry.

01 COAL, LIGNITE, AND PEAT↗

Ion Exchange Processing of AN-107 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System

The Hanford Site stores an estimated 56 million gallons of mixed radioactive and chemically hazardous waste in large underground tanks. In support of the Direct Feed Low-Activity Waste (DFLAW) Program for expediting Hanford tank waste supernate treatment, laboratory-scale ion exchange processing using prototypic unit operations was conducted on AN-107 tank waste at the Pacific Northwest National Laboratory Radiochemical Processing Laboratory. This report describes the small-scale ion exchange testing with 13.7 L of diluted and filtered supernate from Tank 241-AN-107 (hereafter referred to as AN-107) at 16 °C (62 °F). One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AN-107 tank waste to meet this criterion, only 0.147% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 678. Testing with AN-107 matched current Tank Side Cesium Removal (TSCR) facility prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitanate (CST) bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 40% Cs breakthrough after processing ~1700 BVs of feed; the 50% Cs breakthrough was extrapolated from the breakthrough data to occur at 1873 BVs. Testing compared to previous AP-101 and AP-107 testing at 16 °C showed ~300 BV increases in volume processed to reach the WAC limit for both lead and lag columns. The increase in capacity was determined to be due to the significantly lower K concentration in the AN-107 compared to the other tank waste matrices. A comparison in breakthrough curves for the three tests indicated slightly slower kinetic behavior in the AN-107, with variations in feed matrices (high organic complexants) likely responsible for the deviation. The Cs effluent from the lag column reached the WAC limit after processing 1097 BVs. Anticipating this breakthrough point, the polish column was preemptively installed around 900 BVs. Cs breakthrough from the lag column began at 500 BVs, reaching 3.06×10 0 µCi/mL, or 2.6 % Cs breakthrough, after processing all 1700 BVs of feed. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ion Exchange Processing of AW-105 Hanford Tank Waste through Crystalline Silicotitanate in a Staged 2- then 3-Column System

The Hanford Site stores an estimated 56 million gallons of mixed radioactive and chemically hazardous waste in large underground tanks. In support of the Direct Feed Low-Activity Waste (DFLAW) Program for expediting Hanford tank waste supernate treatment, laboratory-scale ion exchange processing using prototypic unit operations was conducted on AW-105 tank waste at the Pacific Northwest National Laboratory Radiochemical Processing Laboratory. This report describes the small-scale ion exchange testing with 9.2 L of diluted and filtered supernate from Tank 241-AW-105 (hereafter referred to as AW-105) at 16 °C (62 °F). One of the waste acceptance criteria (WAC) for the Waste Treatment Plant (WTP) Low-Activity Waste Facility is that the waste must contain less than 3.18×10 -5 Ci 137 Cs per mole of Na. For the AW-105 tank waste to meet this criterion, only 0.225% of the influent 137 Cs concentration may be delivered to the WTP; this requires a Cs decontamination factor of 445. Testing with AW-105 matched current Tank Side Cesium Removal (TSCR) facility prototypic operations where a lead-lag configuration was used until the lag column reached the WAC limit, then a polish column was brought online for continued processing in a lead-lag-polish column configuration. Feed was processed at 1.9 bed volumes (BVs) per hour; the flowrate, in terms of contact time with the crystalline silicotitanate (CST) bed, matched the expected flowrate at TSCR. The Cs-decontaminated product was retained for vitrification testing (to be reported separately). The lead column reached 83% Cs breakthrough after processing ~1500 BVs of feed; the 50% Cs breakthrough was interpolated from the breakthrough data and occurred at 1041 BVs. Despite the AW-105 having a significantly higher K concentration (0.55 M compared to 0.10 M), testing compared to previous AP-107 ion exchange column testing at 16 °C showed no difference in BVs processed to reach the WAC on the lead column and only an approximate ~20 BV decrease in volume processed to reach the WAC limit on the lag column. The negligible differences in capacity despite the 5x concentration differences in K was determined to be due to the significantly lower NO3 concentration in the AW-105 supernate compared to the AP-107 tank waste matrix. A comparison in breakthrough curves for the two tests also indicated slightly faster kinetic behavior in the AW-105, with the variations in feed matrices (lower NO3 concentration) likely responsible for the deviation. The Cs effluent from the lag column reached the WAC limit after processing 772 BVs. Anticipating this breakthrough point, the polish column was preemptively installed around 675 BVs. Cs breakthrough from the lag column began at 300 BVs, reaching 1.10×10 1 µCi/mL, or 14.13 % Cs breakthrough, after processing all 1500 BVs of feed. The polish column processed nominally 830 BVs and reached 2.10×10 -1 µCi/mL, or 0.27 % Cs breakthrough at the conclusion of the test. Table S.1 and Figure S.1 summarize the observed column performance and relevant Cs loading characteristics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Environmental life cycle assessment of treatment and management strategies for food waste and sewage sludge

A consequential life cycle assessment (LCA) was utilized to compare the environmental impacts of food waste and sewage sludge management strategies. The strategies included a novel two-phase anaerobic digestion (AD) system and alternatives including landfill, waste-to-energy, composting, anaerobic membrane bioreactor, and conventional AD (wet continuous stirred-tank reactor [CSTR]). The co-management of food waste with sewage sludge was also considered for the two-phase AD system and for a conventional AD reactor. Further, a multidimensional LCA approach was taken, considering the five-midpoint impact categories of global warming, smog, human health particulate, acidification, and eutrophication estimated using the U.S. EPA Tool for the Reduction and Assessment of Chemical and Other Environmental Impacts. Co-management of food waste and sewage sludge using the novel two-phase AD system was shown to maximize energy recovery and had a net global warming benefit while reducing other environmental impacts compared with the alternative management strategies. It had similar relative environmental advantages across all categories as conventional AD, with the advantage of a smaller physical footprint. However, both approaches featured net environmental burdens when the background electric grid intensity fell below 0.25 kg CO 2 -eq kWh -1 , as could be expected in a decarbonized electric future. Upgrading the biogas produced from AD to renewable natural gas can displace the use of fossil natural gas for other non-electricity energy requirements that are difficult to decarbonize and may extend the time period of significant environmental benefits of utilizing AD for organic waste management. Treatment of the nutrient-rich supernatant generated by the novel two-phase AD system could be an obstacle for utilities with stringent nutrient discharge limits. Future research and full-scale implementation are needed to demonstrate the benefits of the two-phase AD system predicted through this analysis.

54 ENVIRONMENTAL SCIENCES↗

A New Optical Aerosol Spectrometer

An optical particle spectrometer capable of measuring aerosol particle size distributions from 0.02 to 100 micrometers has been developed. This instrument combines several optical methods in one, in-situ configuration; it can provide continuous data collection to encompass the wide dynamic size ranges and concentrations found in studies of modeled planetary atmospheres as well as terrestrial air quality research. Currently, the system is incorporated into an eight liter capacity spherical pressure vessel that is appropriate both for flowthrough and for in-situ particle generation. The optical sizing methods include polarization ratio, The scattering, and forward scattering detectors, with illumination from a fiber-coupled, Argon-ion laser. As particle sizes increase above 0.1 micrometer, a customized electronics and software system automatically shifts from polarization to diffraction-based measurements as the angular scattering detectors attain acceptable signal-to-noise ratios. The number concentration detection limits are estimated to be in the part-per-trillion (ppT by volume) range, or roughly 1000 submicron particles per cubic centimeter. Results from static experiments using HFC134A (approved light scattering gas standard), flow-through experiments using sodium chloride (NaCl) and carbon particles, and dynamic 'Tholin' (photochemical produced particles from ultraviolet (UV)-irradiated acetylene and nitrogen) experiments have been obtained. The optical spectrometer data obtained with particles have compared well with particle sizes determined by electron microscopy. The 'Tholin' tests provided real-time size and concentration data as the particles grew from about 30 nanometers to about 0.8 micrometers, with concentrations ranging from ppT to ppB, by volume. Tests are still underway, to better define sizing accuracy and concentration limits, these results will be reported.

Fonda, Mark↗

A comparison of pulsed and continuous lasers for high-temperature Raman measurements of anhydrite

Oxygen carrier particles (OCPs) serve as fuel oxidant in emerging chemical looping combustion systems. However, chemical looping combustion process optimization is hindered by the lack of online sensors for measurements of OCP oxidation states at temperatures up to 1,000 °C and pressures up to 10 atm. We are investigating Raman spectroscopy as a potential solution, as this technique is known for its ability to provide noninvasive, molecularly specific information in real time in a wide variety of applications. As a case study, Raman spectra from high–temperature (>1,000 °C) anhydrite, a potential OCP, are presented using pulsed and continuous wave excitation at 532 nm. When pulsed excitation is coupled with time–gated detection, background signal from thermal radiation can be significantly reduced from single pulse spectra, outperforming results using an ungated detector. Also, laser heating of the sample is not observed but laser–induced breakdown spectra increasingly appear with increasing temperature and pulse intensity. Compared with pulsed excitation, continuous wave excitation with longer acquisition time, offers higher signal to noise and avoids the risk for laser–induced breakdown spectra but demonstrates minor sample heating. Finally, the data demonstrate practical utility by (a) providing a calibration of anhydrite's ν 1 (1,017 cm –1 ) Raman temperature–dependent band position and by (b) providing estimates of temperature using the ratios of the Stokes/anti–Stokes ν 1 integrated intensities. Furthermore, these results demonstrate the first high–temperature, pulsed Raman spectra from anhydrite while evaluating challenges associated with high–temperature, Raman measurements of OCP materials in general.

20 FOSSIL-FUELED POWER PLANTS↗

Integration of fixed-frequency and FM-CW (frequency-modulated continuous-wave) reflectometers for coincident turbulence measurements on LTX- β (Lithium Tokamak eXperiment- β )

The fixed-frequency and frequency-modulated continuous-wave (FM-CW) reflectometers on LTX-β (Lithium Tokamak eXperiment-β) have been configured to use the same transmission lines and antenna arrays for coincident views of the core and edge plasma. The fixed-frequency channels (13.1–20.5 and 20–40 GHz, tunable between discharges) provide time-resolved measurements of density fluctuations, while the FM-CW channels (13.1–20.2 and 19.5–33.5 GHz) measure the density profile and fluctuations, with high spatial resolution and a sampling rate determined by the frequency sweep interval (5 μs). Data from both reflectometers are synchronously acquired to simultaneously leverage the wide bandwidth and high spatial resolution of the respective systems. Experiments showed that mutual crosstalk interference is momentary and does not diminish the capability of either system. Spectral analysis indicated broad power spectra (several hundreds of kHz) and suggests that the signals from the FM-CW system are consistent with under-sampled fixed-frequency signals. Radial correlations were explored using data from the two reflectometers, as well as from the FM-CW system alone. The core channels showed high levels of agreement between these two comparisons, suggesting that the data from the reflectometers are interchangeable for statistical estimates. For the edge channels, comparisons using data from the FM-CW reflectometer alone showed significant decorrelation due to time lag caused by the finite frequency up-sweep duration. Alternatively, this effect is eliminated when cross-correlating data from the different reflectometers. Finally, these results highlight the advantages of operating the fixed-frequency and FM-CW reflectometers in this manner, where the combined system can overcome the limitations of each separate system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Catalytic Bio-crude Production in a Novel, Short-Contact Time Reactor

RTI has developed a catalytic fast pyrolysis (CFP) process to produce biomass-based hydrocarbon replacement fuels. From a technology perspective, this advanced biofuels process produces liquid transportation fuels that can leverage capital expenditures in the existing petroleum distribution infrastructure. From a business perspective, this technology does not face potential market limitations caused by oxygenate blend limitations for light duty vehicle applications. A promising catalyst has been developed with excellent potential for deoxygenation to produce a hydrocarbon-rich biocrude that is more suitable for upgrading with traditional hydroprocessing technology. The role of the catalyst is to control the chemistry during biomass pyrolysis to produce a biocrude in a single step that has lower oxygen content and is more thermally stable than conventional biomass fast pyrolysis oil. Catalyst properties are optimized to minimize gas and coke production and improve catalytic deoxygenation and biocrude yields. Results from this project have demonstrated that biocrudes from catalytic fast pyrolysis contain <20 wt% oxygen, with the organic fraction as low as 12% oxygen. These oxygen numbers are > 50% lower compared to traditional fast pyrolysis oils. The biocrude produced from our CFP process can be upgraded to gasoline range hydrocarbons. A very conservative approach for upgrading included a preliminary stabilization step followed by a hydroprocessing step. The oxygen content in the upgraded product was decreased by 97% (to less than 3wt %) from the oxygen content of the biocrude sample. Additionally, the acid content was reduced by 99% from the acid content of the biocrude sample. In an optimized integrated process, the CFP step and the hydroprocessing (upgrading) step need to be carefully balanced to maximize the input biomass energy that is recovered in the finished biofuels. Hydrogen utilization and carbon recovery during both catalytic biomass pyrolysis and biocrude upgrading need to be optimized to maximize energy recovery in the finished biofuel. A nominal 1 ton per day (TPD) catalytic biomass pyrolysis system was designed based on the laboratory experimental data. The system was designed to achieve the short residence times (0.5–2 sec) and high heat transfer rates required for maximum liquid biocrude yields while optimizing process integration to maintain catalyst activity by continuous regeneration. The catalytic biomass pyrolysis reactor is a continuously circulating single-loop transport reactor design that is flexible enough to allow sensitivity studies around temperature, residence time, biomass feed rate, and catalyst-to-biomass ratio (i.e., catalyst circulation rate) for process optimization. For a nominal 100 lb/hr feed rate of dry biomass, the products include ~25 lb/hr char and coke, ~20 lb/hr permanent gases, and ~55 lb/hr condensable vapors. Pyrolysis temperatures are 350–600°C, and the system is operated close to ambient pressure. Production of the catalyst developed within the project was scaled up to provide enough material (200kg) for bench-scale testing. The 1TPD bench-scale unit was fabricated, installed and commissioned with this project. Biocrude samples from wood were produced in this unit with the CFP catalyst with material balance > 90%. The intermediate biocrude yields and chemical composition compare well with previous laboratory results and a hydroprocessing strategy for upgrading the resulting biocrude was developed. The experimental results obtained within this project were used to develop a process model for a 2000 TPD commercial facility. Detailed flow sheets were developed for all major unit operations in the process and provided a basis for an economic analysis to develop capital and operating costs. Estimates of the biofuel production costs were used to develop internal rates of return for the yields and compositions measured for the current state of technology at the end of this project and for the yields anticipated for an optimized (nth plant) scenario as additional process improvements are made.

09 BIOMASS FUELS↗

Demonstrating Assurance of Model-Based Fault Diagnosis Systems on an Operational Mission

Developers of robotic scientific and commercial spacecraft are trending towards use of onboard autonomous capabilities for responding quickly to dynamic environments and rapidly changing situations. These capabilities need to know the state of the spacecraft’s health. Model-based fault diagnosis (MBFD) is an approach to estimating health by continuously verifying accurate behavior and diagnosing off-nominal behavior. Proper functioning of MBFD depends on 1) the quality of the diagnostic system model that is analyzed and compared to commands and onboard measurements to estimate a system’s health state, and 2) the correct functionality of the diagnosis engine interrogating the model and comparing its analyses to observed system behavior. Our goal is to develop Verification and Validation (V&V) techniques for MBFD to provide future missions sufficient confidence in its functionality and performance to deploy it on the systems they develop. Our work has been focused on infusing the techniques we developed earlier to an operational mission. First, we are constructing diagnostic models of a spacecraft attitude control system and updating our diagnostic engine so they can be demonstrated aboard the Arcsecond Space Telescope Enabling Research in Astrophysics (ASTERIA) mission, an operational spacecraft for which experiments in autonomy are being planned and executed, using the V&V techniques we have previously developed to assure they are both correct and complete. Since it is nearing the end of its life, ASTERIA provides a unique opportunity to demonstrate MBFD since the monitored components are expected to fail. Our demonstration will give system developers additional confidence to make timely, informed MBFD deployment decisions. Second, we will be completing performance assessments of the diagnostic engine/diagnostic model ensemble both on the flight system and ground-based testbeds to gain confidence in MBFD’s ability to run successfully in a spacecraft’s resource-constrained environment without adversely affecting other on-board activities. Finally, we are capturing our experience in preparing this demonstration in a set of checklists and guidance documents. Current practice includes high-level institutional guidance documents and standards, but at a high level of abstraction that does not necessarily address specific MBFD concerns. The purpose of the new checklists is to provide future mission developers clear, unambiguous, procedure-oriented guidance on assuring MBFD. This paper describes our work in these areas. For the first area, we describe the diagnostic models and updated diagnostic engine that will be used for the on-board demonstration. We describe how the V&V techniques we developed earlier are used to assure model and engine correctness and completeness. For the second area, we identify the performance measurement and assessment techniques used to characterize the diagnostic engine and diagnostic models, and discuss the effect of measured performance on overall mission operation. Finally, we present the checklist and guidance documents and describe how they meet the goals of providing system developers with clear, unambiguous, procedure-oriented guidance on MBFD assurance. We show how the techniques we have developed map into those artifacts.

Nikora, Allen↗

Mass Flux in the Ancient Earth-Moon System and Benign Implications for the Origin of Life on Earth

The origin of life on Earth is commonly considered to have been negatively affected by intense impacting in the Hadean, with the potential for the repeated evaporation and sterilization of any ocean. The impact flux is based on scaling from the lunar crater density record, but that record has no tie to any absolute age determination for any identified stratigraphic unit older than approx. 3.9 Ga (Nectaris basin). The flux can be described in terms of mass accretion, and various independent means can be used to estimate the mass flux in different intervals. The critical interval is that between the end of essential crustal formation (approx. 4.4 Ga) and the oldest mare times (approx. 3.8 Ga). The masses of the basin-forming projectiles during Nectarian and early Imbrian times, when the last 15 of the approx.45 identified impact basins formed, can be reasonably estimated as minima. These in sum provide a minimum of 2 x 10(exp 21)g for the mass flux to the Moon during those times. If the interval was 80 million years (Nectaris 3.90 Ga, Orientale 3.82 Ga), then the flux was approx. 2 x 10(exp 13) g/yr over this period. This is higher by more than an order of magnitude than a flux curve that declines continuously and uniformly from lunar accretion to the rate inferred for the older mare plains. This rate cannot be extrapolated back increasingly into pre-Nectarian times, because the Moon would have added masses far in excess of itself in post-crust-formation time. Thus this episode was a distinct and cataclysmic set of events. There are approx. 30 pre-Nectarian basins, and they were probably part of the same cataclysm (starting at approx. 4.0 Ga?) because the crust is fairly intact, the meteoritic contamination of the pre-Nectarian crust is very low, impact melt rocks older than 3.92 Ga are virtually unknown, and ancient volcanic and plutonic rocks have survived this interval. The accretionary flux from approx. 4.4 to approx. 4.0 Ga was comparatively benign. When scaled to Earth, even the late cataclysm does not produce oceane vaporating, globally sterilizing events. The rooted concept that such events took place is based on the extrapolation of a nonexistent lunar record to the Hadean. The Earth from approx. 4.4 to approx. 3.8 Ga was comparatively peaceful, and the impacting itself could have been thermally and hydrothermally beneficial. The origin of life could have taken place at any time between 4.4 and 3.85 Ga, given the current impact constraints, and there is no justification for the claim that life originated (or re-originated) as late as 3.85 Ga in response to the end of hostile impact conditions.

Ryder, Graham↗

Completion of lunar magma ocean solidification at 4.43 Ga

Crystallization of the lunar magma ocean yielded a chemically unique liquid residuum named KREEP. This component is expressed as a large patch on the near side of the Moon and a possible smaller patch in the northwest portion of the Moon’s South Pole-Aitken basin on the far side. Thermal models estimate that the crystallization of the lunar magma ocean (LMO) could have spanned from 10 and 200 My, while studies of radioactive decay systems have yielded inconsistent ages for the completion of LMO crystallization covering over 160 My. Here, we show that the Moon achieved >99% crystallization at 4,429 ± 76 Ma, indicating a lunar formation age of ~4,450 Ma or possibly older. Using the 176 Lu– 176 Hf decay system (t 1/2 = 37 Gy), we found that the initial 176 Hf/ 177 Hf ratios of lunar zircons with varied U–Pb ages are consistent with their crystallization from a KREEP-rich reservoir with a consistently low 176 Lu/ 177 Hf ratio of 0.0167 that emerged ~140 My after solar system formation. The previously proposed younger model age of ~4.33 Ga for the source of mare basalts (240 My after solar system formation) might reflect the timing of a large impact. Our results demonstrate that lunar magma ocean crystallization took place while the Moon was still battered by planetary embryos and planetesimals leftover from the main stage of planetary accretion. The study of Lu–Hf model ages for samples brought back from the South Pole-Aitken basin will help to assess the lateral continuity of KREEP and further understand its significance in the early history of the Moon.

Science & Technology - Other Topics↗

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗

Orbital elements of 4U 0115+63 and the nature of the hard X-ray transients

Extended SAS 3 timing observations of the hard transient X-ray source 4U 0115+63 are reported, and a definitive measurement of the binary orbit of this transient source is presented. It is shown that this source is in a long orbit (period of approximately 24.3 days) that is moderately eccentric (e about 0.34) and that the mean value of the rate of decrease of the pulse period is consistent with the expected spinup of a rotating neutron star that is accreting from a disk. A distance of about 2.5 kpc is inferred, and the B-star optical counterpart is estimated to have an absolute magnitude of approximately -1.5 and a mass of at least 5 solar masses. It is suggested that the companion is a Be star which does not fill its Roche lobe and that the eccentricity and transient nature of the source result from the large orbital separation. It is proposed that hard X-ray transients as a class are collapsed stars (perhaps all neutron stars) in binary systems that are substantially wider than the more persistent X-ray binaries and that the large orbital separation, the small radius of the companion, or both, result in episodic rather than continuous mass transfer onto the X-ray star.

Rappaport, S.↗

Operational Evaluatioin of Dynamic Weather Routes at American Airlines

Dynamic Weather Routes (DWR) is a search engine that continuously and automatically analyzes inflight aircraft in en route airspace and proposes simple route amendments for more efficient routes around convective weather while considering sector congestion, traffic conflicts, and active Special Use Airspace. NASA and American Airlines (AA) are conducting an operational trial of DWR at the AA System Operations Center in Fort Worth, TX. The trial includes only AA flights in Fort Worth Center airspace. Over the period from July 31, 2012 through August 31, 2012, 45% of routes proposed by DWR and evaluated by AA users - air traffic control coordinators and flight dispatchers - were rated as acceptable as proposed or with some modifications. The wind-corrected potential flying time savings for these acceptable routes totals 470 flying min, and results suggest another 1,500 min of potential savings for flights not evaluated due to staffing limitations. A sector congestion analysis shows that in only two out of 83 DWR routes rated acceptable by AA staff were the flights predicted to fly through a congested sector inside of 30 min downstream of present position. This shows that users considered sector congestion data provided by DWR automation and in nearly all cases did not accept routes through over-capacity sectors. It is estimated that 12 AA flights were given reroute clearances as a direct result of DWR for a total savings of 67 flying min.

air traffic trajectory automation↗

Hybrid Electric Aircraft Thermal Management: Now, New Visions and Future Concepts and Formulation

The global fuel consumption by commercial airlines has increased each year since 2009 and is predicted to reach an all-time high of 97 billion gallons in 2019. There is also an environmental impact from this: CO2 emissions from commercial passenger and freight operations totaled 918 Mt in 2018 (ICCT, 2019), or around 2.5% of global energy-related CO2 emissions. Passenger transport accounted for 81% of the total. Emissions from aviation have grown 32% over the past five years. Coupled with this aspect, there is a continuous and growing need to satisfy ever-growing electrical power needs on commercial and military aircraft. All the armed services (Army, Navy, and Air Force) are continuously trying to enhance UAV (unmanned aerial vehicle) endurance and range across a broad fleet of different aircraft. The commercial Boeing 787 requires about 1.2MWe and that is expected to grow. Current technologies used to supply increased on-board electrical power are generally: 1) “burn more fuel and convert through on-board generators” and 2) use additional heavy (i.e., weight-inefficient) and sometimes unsafe battery systems on-board the aircraft. The aircraft industry is seeking new, innovative ways to satisfy this increasing power demand. One as-yet-untapped power source is the enormous amount of “waste” thermal energy flowing out the jet engine exhaust; some estimates in smaller “by-pass” flow jet engines is several hundreds of kilowatts (e.g., Pratt & Whitney Canada PW545B turbofan). This quantity is much higher in large jet engines associated with commercial aircraft. This large waste thermal energy manifests itself in large temperature differences within the by-pass-flow engine exhaust system relative to outside ambient conditions, because of the actual by-pass engine design configuration. There is strong need to develop thermal technologies and systems that could harness and convert at least a portion of this thermal energy into useful electrical energy to satisfy growing on-board electrical needs. In addition, there is a strong desire within the aircraft and engine manufacturing community to reduce the “carbon footprint” of the industry though reduced fuel usage worldwide. NASA has a robust aircraft electrification program to meet these desires and support industry in its aircraft electrification objectives. This program is integrating thermoacoustic systems, advanced lightweight heat exchanger technology, and advanced heat pipe technology to capture and transport large amounts of engine waste thermal energy for on-board power conversion, advanced heat-pump cooling, and exergy enhancement (i.e., temperature lift). Advanced lightweight heat exchangers are envisioned to capture engine exhaust thermal energy at approximately 673 K and deliver it to efficient thermoacoustic power conversion systems operating at temperature ratios (Thot/Tcold) > 1.6. Advanced heat pipe systems are envisioned to transport thermal energy from low temperature sources, through thermoacoustic heat pumps, to high temperature needs such as wing anti-icing, fuel pre-heating, and combustion air pre-heating. The paper will discuss the current state-of-the-art, objectives, system design architecture, and remaining technical challenges in system formulation within the NASA aircraft electrification program.

Dyson, Rodger↗

Towards real-time monitoring: data assimilated time-lapse full waveform inversion for seismic velocity and uncertainty estimation

SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

58 GEOSCIENCES↗

SPARR: Spiking/Processing Array for Wide Dynamic Range and High Resolution Photonic Sensing

The Spiking/Processing Array (spARR) is a novel photonic focal plane that uses pixels which generate electronic spikes autonomously and without a clock. These spikes feed into a network of digital asynchronous processing elements or DAPES. By building a useful assemblage of DAPES, and connecting them together in the correct way, sophisticated signal processing can be accomplished within the focal plane. Autonomous self-resetting pixels (AsP) enable SPARR to generate electronic response with very small signals--as little as a single photon in the case of Geiger mode avalanche photodiodes to as few as several hundred photons for in-cmos photodetectors. These spiking pixels enable fast detector response, but do not draw as much continuous power as synchronous clocked designs. The spikes emitted by the pixels all have the same magnitude, the information from the scene is effectively encoded into the rate of spikes and the time at which the spike is emitted. The spiking pixels, having converted incident light into electronic spikes, supply the spikes to a network of digital asynchronous processors. These are small state machines which respond to the spikes arriving at their input ports by either remaining unchanged or updating their internal state and possibly emitting a spike on one or more output ports. We show a design that accomplishes the sophisticated signal processing of a Haar spatial wavelet transform with spatial-spectral whitening. We furthermore show how this design results in a data streams which support imaging and transient optical source detection. Two simulators support this analysis: SPICE and sparrow. The CMOS SPICE simulator Cadence provides accurate CMOs design with accounting for effects of circuit parasitics throughout layout, accurate timing, and accurate energy consumption estimates. To more rapidly assess larger networks with more pixels, sparrow is a custom discrete event simulator that supports the non-homogeneous Poisson processes that lie behind photoelectric interaction. Sparrow is a photon-exact simulator that nevertheless performs SPARR system simulator for large-scale systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗