Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained optimization”

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

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

At least 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evaluating adaptive and predictive power management strategies for optimizing visualization performance on supercomputers

Power is becoming an increasingly scarce resource on the next generation of supercomputers, and should be used wisely to improve overall performance. One strategy for improving power usage is hardware overprovisioning, i.e., systems with more nodes than can be run at full power simultaneously without exceeding the system-wide power limit. With this study, we compare two strategies for allocating power throughout an overprovisioned system – adaptation and prediction – in the context of visualization workloads. While adaptation has been suitable for workloads with more regular execution behaviors, it may not be as suitable on visualization workloads, since they can have variable execution behaviors. This study considers a total of 104 experiments, which vary the rendering workload, power budget, allocation strategy, and node concurrency, including tests processing data sets up to 1 billion cells and using up to 18,432 cores across 512 nodes. Overall, we find that prediction is a superior strategy for this use case, improving performance up to 27% compared to an adaptive strategy.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Advanced timeline systems

The Mission Planning Division of the Mission Operations Laboratory at NASA's Marshall Space Flight Center is responsible for scheduling experiment activities for space missions controlled at MSFC. In order to draw statistically relevant conclusions, all experiments must be scheduled at least once and may have repeated performances during the mission. An experiment consists of a series of steps which, when performed, provide results pertinent to the experiment's functional objective. Since these experiments require a set of resources such as crew and power, the task of creating a timeline of experiment activities for the mission is one of resource constrained scheduling. For each experiment, a computer model with detailed information of the steps involved in running the experiment, including crew requirements, processing times, and resource requirements is created. These models are then loaded into the Experiment Scheduling Program (ESP) which attempts to create a schedule which satisfies all resource constraints. ESP uses a depth-first search technique to place each experiment into a time interval, and a scoring function to evaluate the schedule. The mission planners generate several schedules and choose one with a high value of the scoring function to send through the approval process. The process of approving a mission timeline can take several months. Each timeline must meet the requirements of the scientists, the crew, and various engineering departments as well as enforce all resource restrictions. No single objective is considered in creating a timeline. The experiment scheduling problem is: given a set of experiments, place each experiment along the mission timeline so that all resource requirements and temporal constraints are met and the timeline is acceptable to all who must approve it. Much work has been done on multicriteria decision making (MCDM). When there are two criteria, schedules which perform well with respect to one criterion will often perform poorly with respect to the other. One schedule dominates another if it performs strictly better on one criterion, and no worse on the other. Clearly, dominated schedules are undesireable. A nondominated schedule can be generated by some sort of optimization problem. Generally there are two approaches: the first is a hierarchical approach while the second requires optimizing a weighting or scoring function.

Bulfin, R. L.↗
DOE OSTI2025

PET-FBA: A lightweight enzyme allocation and thermodynamics-constrained flux analysis approach to explore Escherichia coli metabolic adaptation to intracellular acidification

Escherichia coli employs diverse strategies to adapt to acidic environments that disrupt enzyme activity and the thermodynamic feasibility of essential reactions. To understand the impact of pH stress on cell metabolism, we present the PET-FBA (pH-, Enzyme protein allocation-, and Thermodynamics-constrained Flux Balance Analysis) framework. PET-FBA extends genome-scale modeling by integrating enzyme protein costs and reaction Gibbs free energy changes. Additionally, by incorporating pH-dependent enzyme kinetics in response to intracellular acidification, this framework enables the simulation of E. coli's metabolic adjustments across varying external pH levels. The model's accuracy is validated by comparing in silico growth simulations with experimental measurements under both anaerobic and aerobic conditions, as well as in silico gene knockouts of essential genes. By explicitly incorporating pH effects, our model accurately replicates the metabolic shift towards lactate production as the primary fermentation product at low pH in anaerobic conditions. This shift is only predicted when enzyme kinetics are dynamically adjusted as a function of pH. Further analysis revealed that this shift can be attributed to the reduced protein efficiency of the acetyl-CoA branch compared to lactate dehydrogenase under acidic stress, which then becomes crucial for maintaining NAD regeneration and cell growth at low pH. Furthermore, we identified strategies for enhancing cell growth under acidic anaerobic conditions by improving the enzyme activity of lactate dehydrogenase and pyruvate formate lyase, which increases NAD production efficiency and reduces enzyme protein allocation costs. Designed as a lightweight yet versatile framework, PET-FBA enables efficient genome-scale metabolic analysis. Using E. coli as a model system, our framework provides a systematic approach to understanding metabolic responses to environmental stress, pinpointing key metabolic bottlenecks, and identifying potential targets for strain optimization.

42 ENGINEERING↗
DOE OSTI2022

A shared-mobility-based framework for evacuation planning and operations under forecast uncertainty

To meet evacuation needs from carless populations who need personalized assistance to evacuate safely, in this article we propose a ridesharing-based evacuation program that recruits volunteer drivers before a disaster strikes, and then matches volunteer drivers with evacuees once demand is realized. Here we optimize resource planning and evacuation operations under uncertain spatiotemporal demand, and construct a two-stage stochastic mixed-integer program to ensure high demand fulfillment rates. We consider three formulations to improve the number of evacuees served, by minimizing an expected penalty cost, imposing a probabilistic constraint, and enforcing a constraint on the conditional value at risk of the total number of unserved evacuees, respectively. We discuss the benefits and disadvantages of the different risk measures used in the three formulations, given certain carless population sizes and the variety of evacuation modes available. We also develop a heuristic approach to provide quick, dynamic and conservative solutions. We demonstrate the performance of our approaches using five different networks of varying sizes based on regions of Charleston County, South Carolina, an area that experienced a mandatory evacuation order during Hurricane Florence, and utilize real demographic data and hourly traffic count data to estimate the demand distribution.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2020

Revisiting the Effectiveness Of HCHO/NO2 Ratios for Inferring Ozone Sensitivity to Its Precursors Using High Resolution Airborne Remote Sensing Observations in A High Ozone Episode During the KORUS-AQ Campaign

The nonlinear chemical processes involved in ozone production (P(O3)) have necessitated using proxy indicators to convey information about the primary dependence of P(O3) on volatile organic compounds (VOCs) or nitrogen oxides (NOx). In particular, the ratio of remotely sensed columns of formaldehyde (HCHO) to nitrogen dioxide (NO2) has been widely used for studying O3 sensitivity. Previous studies found that the errors in retrievals and the incoherent relationship between the column and the near-surface concentrations are a barrier in applying the ratio in a robust way. In addition to these obstacles, we provide calculational-observational evidence, using an ensemble of 0-D photochemical box models constrained by DC-8 aircraft measurements on an ozone event during the Korea-United States Air Quality (KORUS-AQ) campaign over Seoul, to demonstrate the chemical feedback of NO2 on the formation of HCHO is a controlling factor for the transition line between NOx-sensitive and NOx-saturated regimes. A fixed value (~2.7) of the ratio of the chemical loss of NOx (LNOx) to the chemical loss of HO2+RO2 (LROx) perceptibly differentiates the regimes. Following this value, data points with a ratio of HCHO/NO2 less than 1 can be safely classified as NOx-saturated regime, whereas points with ratios between 1 and 4 fall into one or the other regime. We attribute this mainly to the HCHO-NO2 chemical relationship causing the transition line to occur at larger (smaller) HCHO/NO2 ratios in VOC-rich (VOC-poor) environments. We then redefine the transition line to LNOx/LROx~2.7 that accounts for the HCHO-NO2 chemical relationship leading to HCHO = 3.7 × (NO2 – 1.14 × 1016 molec.cm-2). Although the revised formula is locally calibrated (i.e., requires for readjustment for other regions), its mathematical format removes the need for having a wide range of thresholds used in HCHO/NO2 ratios that is a result of the chemical feedback. Therefore, to be able to properly take the chemical feedback into consideration, the use of HCHO = a × (NO2 – b) formula should be preferred to the ratio in future works. We then use the Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) airborne instrument to study O3 sensitivity in Seoul. The unprecedented spatial (250 × 250 m2) and temporal (~every 2 h) resolutions of HCHO and NO2 observations form the sensor enhance our understanding of P(O3) in Seoul; rather than providing a crude label for the entire city, more in-depth variabilities in chemical regimes are observed that should be able to inform mitigation strategies correspondingly.

KORUS-AQ; formaldehyde; NOx; VOC; OMI; GeoTASO; oz↗
DOE OSTI2020

Conserving Land–Atmosphere Synthesis Suite (CLASS)

Accurate estimates of terrestrial water and energy cycle components are needed to better understand climate processes and improve models’ ability to simulate future change. Various observational estimates are available for the individual budget terms; however, these typically show inconsistencies when combined in a budget. In this work, a Conserving Land–Atmosphere Synthesis Suite (CLASS) of estimates of simultaneously balanced surface water and energy budget components is developed. Individual CLASS variable datasets, where possible, 1) combine a range of existing variable product estimates, and hence overcome the limitations of estimates from a single source; 2) are observationally constrained with in situ measurements; 3) have uncertainty estimates that are consistent with their agreement with in situ observations; and 4) are consistent with each other by being able to solve the water and energy budgets simultaneously. First, available datasets of a budget variable are merged by implementing a weighting method that accounts both for the ability of datasets to match in situ measurements and the error covariance between datasets. Then, the budget terms are adjusted by applying an objective variational data assimilation technique (DAT) that enforces the simultaneous closure of the surface water and energy budgets linked through the equivalence of evapotranspiration and latent heat. Comparing component estimates before and after applying the DAT against in situ measurements of energy fluxes and streamflow showed that modified estimates agree better with in situ observations across various metrics, but also revealed some inconsistencies between water budget terms in June over the higher latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Cognition at the Point of Sensing

Over the last 15 years, compressive sensing techniques have been developed which have the potential to greatly reduce the amount of data collected by systems while preserving the amount of information obtained. A cost of this efficiency is that a computationally-intensive optimization routine must be used to put the sensed data into a form that a person can interpret. At the same time, machine learning techniques have experienced tremendous growth as well. Machines have demonstrated the ability learn how to effectively perform tasks such as detection and classification at speeds much faster than humanly possible. Our goal in this project was to study the feasibility of using compressive sensing systems "at the edge." That is, how can compressive sensing sensors be deployed such that information is created at the remote sensor rather than sending raw data to a central processing location? Studies were performed to analyze whether machine learning could be done on the compressively sensed data in its raw form. If a machine is performing the task, is it possible to do so without putting the data into a human interpretable form? We show that this is possible for some systems, in particular a compressive sensing snapshot imaging spectrometer. Machine learning tasks were demonstrated to be more effective and more robust to noise when the machine learning algorithm worked on data in its raw form. This system is shown to outperform a traditional spectrometer. Techniques for reducing the complexity of the reconstruction routine were also analyzed. Techniques for such as data regularization, deep neural networks, and matrix completion were studied and shown to have benefits over traditional reconstruction techniques. In this project we showed that compressive sensing sensors are indeed feasible at the edge. As always, sensors and algorithms must be carefully tuned to work in the constrained environment. In this project we developed tools and techniques to enable those analyses.

47 OTHER INSTRUMENTATION↗
NASA NTRS1999

Chain Dynamics in Magnetorheological Suspensions

Magnetorheological (MR) suspensions are composed of colloidal particles which acquire dipole moments when subjected to an external magnetic field. At sufficient field strengths and concentrations, the dipolar particles rapidly aggregate to form long chains. Subsequent lateral cross-linking of the dipolar chains is responsible for a rapid liquid-to-solid-like rheological transition. The unique, magnetically-activated rheological properties of MR suspensions make them ideal for interfacing mechanical systems to electronic controls. Additionally, the ability to experimentally probe colloidal suspensions interacting through tunable anisotropic potentials is of fundamental interest. Our current experimental work has focused on understanding the fluctuations of dipolar chains. It has been proposed by Halsey and Toor (HT) that the strong Landau-Peierls thermal fluctuations of dipolar chains could be responsible for long-range attractions between chains. Such interactions will govern the long-time relaxation of MR suspensions. We have synthesized monodisperse neutrally buoyant MR suspensions by density matching stabilized ferrofluid emulsion droplets with D2O. This allows us to probe the dynamics of the dipolar chains using light scattering without gravitational, interfacial, and polydispersity effects to resolve the short-wavelength dynamics of the dipolar chains. We used diffusing wave spectroscopy to measure these dynamics. The particle displacements at short times that show an independence to the field strength, but at long times exhibit a constrained, sub-diffusive motion that slows as the dipole strength is increased. The experiments are in good qualitative agreement with Brownian dynamics simulations of dipolar chains. Although there have been several important and detailed studies of the structure and interactions in MR suspensions, there has not been conclusive evidence that supports or contradicts the HT model prediction that long-range interactions exist between fluctuating chains of dipolar particles. Resolving this issue would contribute greatly to the understanding of these interesting and important materials. We have begun to test the predictions of the HT model by both examining the dynamics of individual chains and by measuring the forces between dipolar chains directly to accurately and quantitatively assess the interactions that they experience. To do so, we employ optical trapping techniques and video-microscopy to manipulate and observe our samples on the microscopic level. With these techniques, it is possible to observe chains that are fluctuating freely in three-dimensions, independent of interfacial effects. More importantly, we are able to controllably observe the interactions of two chains at various separations to measure the force-distance profile. The techniques also allow us to study the mechanical properties of individual chains and chain clusters. Our work to this point has focused on reversibly-formed dipolar chains due to field induced dipoles where the combination of this chaining, the dipolar forces, and the hydrodynamic interactions that dictate the rheology of the suspensions. One can envision, however, many situations where optical, electronic, or rheological behavior may be optimized with magneto-responsive anisotropic particles. Chains of polarizable particles may have the best properties as they can coil and flex in the absence of a field and stiffen and orient when a field is applied. We have recently demonstrated a synthesis of stable, permanent paramagnetic chains by both covalently and physically linking paramagnetic colloidal particles. The method employed allows us to create monodisperse chains of controlled length. We observed the stability, field-alignment, and rigidity of this new class of materials. The chains may exhibit unique rheological properties in an applied magnetic field over isotropic suspensions of paramagnetic particles. They are also useful rheological models as bead-spring systems. These chains form the basis for our current experiments with optical traps.

Gast, A. P.↗
NASA NTRS2011

Ground Vibration Testing Options for Space Launch Vehicles

New NASA launch vehicles will require development of robust systems in a fiscally-constrained environment. NASA, Department of Defense (DoD), and commercial space companies routinely conduct ground vibration tests as an essential part of math model validation and launch vehicle certification. Although ground vibration testing must be a part of the integrated test planning process, more affordable approaches must also be considered. A study evaluated several ground vibration test options for the NASA Constellation Program flight test vehicles, Orion-1 and Orion-2, which concluded that more affordable ground vibration test options are available. The motivation for ground vibration testing is supported by historical examples from NASA and DoD. The approach used in the present study employed surveys of ground vibration test subject-matter experts that provided data to qualitatively rank six test options. Twenty-five experts from NASA, DoD, and industry provided scoring and comments for this study. The current study determined that both element-level modal tests and integrated vehicle modal tests have technical merits. Both have been successful in validating structural dynamic math models of launch vehicles. However, element-level testing has less overall cost and schedule risk as compared to integrated vehicle testing. Future NASA launch vehicle development programs should anticipate that some structural dynamics testing will be necessary. Analysis alone will be inadequate to certify a crew-capable launch vehicle. At a minimum, component and element structural dynamic tests are recommended for new vehicle elements. Three viable structural dynamic test options were identified. Modal testing of the new vehicle elements and an integrated vehicle test on the mobile launcher provided the optimal trade between technical, cost, and schedule.

Patterson, Alan↗
NASA NTRSDate not supplied

High-Temperature MEMS Based Venus Seismometer

Seismology is the method of choice for studying a planet’s interior and assessing current tectonic activity. With a single strategically placed seismometer, or preferably a network of seismometers, a variety of key aspects of a planet can be evaluated. The level and nature of seismicity observed provides a gauge of current geologic, especially tectonic, activity. Where tectonic activity is and is not occurring illuminates what regions and geologic features on the planet are currently active, thus, constraining global geodynamics. Larger distal earthquakes can be used to determine major compositional and structural boundaries within a planet’s interior, providing major constraints on planet formation and interior evolution. Seismometers have been deployed successfully on the Moon and Mars. Only small subsets of the scientific knowledge that can be gained from seismology can be obtained from alternative scientific approaches, such as global gravity solutions or repeat-pass interferometry. By virtue of its similar size to Earth, Venus is likely to be highly seismically active (Fig. 1). Even at Earth-like levels of seismology, returning meaningful data would require observation over a period that is at least 2–3 orders of magnitude longer than the 1–2 h lifetime of previous Soviet landers. Clever insulation could extend the lifetime of standard electronics in a modern lander to perhaps 24 h. While active cooling of electronics for a Venusian lander is plausible, this approach would require considerable technological advancement and will likely come at a high monetary and mass cost. The high density of the Venusian atmosphere at the surface should mean that there is good coupling of seismic energy into the atmosphere, such that seismology using infrasound from a balloon platform is currently being explored, and observations of the upper levels of the atmosphere from orbit has also been proposed [3]. Besides losing all shear wave information with these two methods, calibrating and interpreting seismic information that has been distorted by, and convolved with, atmospheric phenomena will be challenging without existing seismology data collected through surface seismometers (optimally, simultaneous surface seismological observations would be made).

silicon carbide↗
NASA NTRS1994

Control of flow separation in airfoil/wing design applications

Existing aerodynamic design methods have generally concentrated on the optimization of airfoil or wing shapes to produce a minimum drag while satisfying some basic constraints such as lift, pitching moment, or thickness. Since the minimization of drag almost always precludes the existence of separated flow, the evaluation and validation of these design methods for their robustness and accuracy when separated flow is present has not been aggressively pursued. However, two new applications for these design tools may be expected to include separated flow and the issues of aerodynamic design with this feature must be addressed. The first application of the aerodynamic design tools is the design of airfoils or wings to provide an optimal performance over a wide range of flight conditions (multipoint design). While the definition of 'optimal performance' in the multipoint setting is currently being hashed out, it is recognized that given a wide enough range of flight conditions, it will not be possible to ensure a minimum drag constraint at all conditions, and in fact some amount of separated flow (presumably small) may have to be allowed at the more demanding flight conditions. Thus a multipoint design method must be tolerant of the existence of separated flow and may include some controls upon its extent. The second application is in the design of wings with extended high speed buffet boundaries of their flight envelopes. Buffet occurs on a wing when regions of flow separation have grown to the extent that their time varying pressures induce possible destructive effects upon the wing structure or adversely effect either the aircraft controllability or the passenger comfort. A conservative approach to the expansion of the buffet flight boundary is to simply expand the flight envelope of nonseparated flow under the assumption that buffet will also thus be alleviated. However, having the ability to design a wing with separated flow and thus to control the location, extent, and severity of the separated flow regions may allow aircraft manufacturers to gain an advantage in the early design stages of an aircraft, when configuration changes are relatively inexpensive to make. Continuing the work begun last year, an airfoil design package has been modified to provide some control over the existence and extent of flow separation. This package consists of a 2-D Navier-Stokes flow solver which is coupled to the CDISC (constrained direct/iterative surface curvature) design method. The first modification is a prediction method for determining whether separation is likely based solely upon a given pressure distribution. If separation is predicted but is undesirable, the new routines will modify the pressure distribution to alleviate the problem. This new pressure distribution is then used in the design method to generate a new aerodynamic shape. Since separation may be acceptable in some cases, particularly if the separation does not extend to the trailing edge, another added logic estimates the extent of separation based upon a correlation with calculated separated flow cases. If the flow behind a shock induced separation is not predicted to reattach before the trailing edge, the logic weakens the shock strength and otherwise alters the pressure distribution in order to promote reattachment. This later addition is as yet unreliable due to secondary separation effects, but additional work is being pursued to improve the method.

Gally, Thomas A.↗
NASA NTRS2001

EFQPSK Versus CERN: A Comparative Study

This report presents a comparative study on Enhanced Feher's Quadrature Phase Shift Keying (EFQPSK) and Constrained Envelope Root Nyquist (CERN) techniques. These two techniques have been developed in recent times to provide high spectral and power efficiencies under nonlinear amplifier environment. The purpose of this study is to gain insights into these techniques and to help system planners and designers with an appropriate set of guidelines for using these techniques. The comparative study presented in this report relies on effective simulation models and procedures. Therefore, a significant part of this report is devoted to understanding the mathematical and simulation models of the techniques and their set-up procedures. In particular, mathematical models of EFQPSK and CERN, effects of the sampling rate in discrete time signal representation, and modeling of nonlinear amplifiers and predistorters have been considered in detail. The results of this study show that both EFQPSK and CERN signals provide spectrally efficient communications compared to filtered conventional linear modulation techniques when a nonlinear power amplifier is used. However, there are important differences. The spectral efficiency of CERN signals, with a small amount of input backoff, is significantly better than that of EFQPSK signals if the nonlinear amplifier is an ideal clipper. However, to achieve such spectral efficiencies with a practical nonlinear amplifier, CERN processing requires a predistorter which effectively translates the amplifier's characteristics close to those of an ideal clipper. Thus, the spectral performance of CERN signals strongly depends on the predistorter. EFQPSK signals, on the other hand, do not need such predistorters since their spectra are almost unaffected by the nonlinear amplifier, Ibis report discusses several receiver structures for EFQPSK signals. It is observed that optimal receiver structures can be realized for both coded and uncoded EFQPSK signals with not too much increase in computational complexity. When a nonlinear amplifier is used, the bit error rate (BER) performance of the CERN signals with a matched filter receiver is found to be more than one decibel (dB) worse compared to the bit error performance of EFQPSK signals. Although channel coding is found to provide BER performance improvement for both EFQPSK and CERN signals, the performance of EFQPSK signals remains better than that of CERN. Optimal receiver structures for CERN signals with nonlinear equalization is left as a possible future work. Based on the numerical results, it is concluded that, in nonlinear channels, CERN processing leads towards better bandwidth efficiency with a compromise in power efficiency. Hence for bandwidth efficient communications needs, CERN is a good solution provided effective adaptive predistorters can be realized. On the other hand, EFQPSK signals provide a good power efficient solution with a compromise in band width efficiency.

Borah, Deva K.↗
DOE OSTI2020

Comprehensive geoneutrino analysis with Borexino

This paper presents a comprehensive geoneutrino measurement using the Borexino detector, located at Laboratori Nazionali del Gran Sasso (LNGS) in Italy. The analysis is the result of 3262.74 days of data between December 2007 and April 2019. The paper describes improved analysis techniques and optimized data selection, which includes enlarged fiducial volume and sophisticated cosmogenic veto. The reported exposure of ( 1.29 ± 0.05 ) × 10 32 protons × year represents an increase by a factor of two over a previous Borexino analysis reported in 2015. By observing 52.6 - 8.6 + 9.4 ( stat ) - 2.1 + 2.7 ( sys ) geoneutrinos (68% interval) from U 238 and Th 232 , a geoneutrino signal of 47.0 - 7.7 + 8.4 ( stat ) - 1.9 + 2.4 ( sys ) TNU with - 17.2 + 18.3 % total precision was obtained. This result assumes the same Th/U mass ratio as found in chondritic CI meteorites but compatible results were found when contributions from U 238 and Th 232 were both fit as free parameters. Antineutrino background from reactors is fit unconstrained and found compatible with the expectations. The null-hypothesis of observing a geoneutrino signal from the mantle is excluded at a 99.0% C.L. when exploiting detailed knowledge of the local crust near the experimental site. Measured mantle signal of 21.2 - 9.0 + 9.5 ( stat ) - 0.9 + 1.1 ( sys ) TNU corresponds to the production of a radiogenic heat of 24.6 - 10.4 + 11.1 TW (68% interval) from U 238 and Th 232 in the mantle. Assuming 18% contribution of K 40 in the mantle and 8.1 - 1.4 + 1.9 TW of total radiogenic heat of the lithosphere, the Borexino estimate of the total radiogenic heat of the Earth is 38.2 - 12.7 + 13.6 TW , which corresponds to the convective Urey ratio of 0.78 - 0.28 + 0.41 . These values are compatible with different geological predictions, however there is a ~ 2.4 σ tension with those Earth models which predict the lowest concentration of heat-producing elements in the mantle. In addition, by constraining the number of expected reactor antineutrino events, the existence of a hypothetical georeactor at the center of the Earth having power greater than 2.4 TW is excluded at 95% C.L. Particular attention is given to the description of all analysis details which should be of interest for the next generation of geoneutrino measurements using liquid scintillator detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2020

Investigation of Magmatic Activities on Early Mars Using Igneous Mineral Chemistry in Gale Crater, Mars

One objective of rover missions is exploring the geological context of the surroundings. Over the years, igneous petrology and sedimentology have been disconnected, the first investigating magmatic processes and volcanic activities, and the second seeking environmental conditions in the past and assessing the habitability of the planet. Although different, one is related to the other: igneous rocks are altered and broken down, leading to the formation of sedimentary rocks, which can in turn be used to back out the nature of their magmatic source. The Curiosity rover that landed in the 3.7 Gyr old impact crater Gale is traveling through sedimentary rocks. About fifty float rocks have been observed, and several of them with ambiguous texture and composition have been classified as igneous or sedimentary depending on studies such as Jake_M. The composition of several unambiguous igneous rocks has been analyzed [4- 6] but their heterogeneity at a larger instrumental (measurement size < 2 cm) scale prevents the measurement of a bulk composition as performed on Earth. An original approach avoiding these two last issues is to consider igneous mineral chemistry analyzed within igneous and sedimentary rocks to assess magmatic processes that could have formed them. Most Curiosity data are used to explore ancient environmental conditions, and a significant number of compositional analyses are under-explored for constraining magmatic activities. We will present how we can make use of sedimentary data for investigating igneous processes in the vicinity of Gale crater. Geological Context: We focus on the first 750 martian days, corresponding to measurements in a coherent lacustrine sedimentary unit called Bradbury, because all sedimentary rocks were sourced from the same watershed and appear to have a consistent source with minimal alteration [2-3]. Igneous detrital minerals including feldspar and pyroxene, are observed in sedimentary rocks. Monte Carlo models showed that minimal cation loss is observed based on the composition of all Bradbury rocks, implying negligible weathering [3]. Although clay minerals are detected in few rocks [7], chemical compositions of rocks can be explained by a mixture of primary igneous minerals [3]. Variation of composition within Bradbury rocks can be explained by mineral sorting and one distinct source component. While a common magmatic source is suggested, Bradbury sediments likely come from several volcanic eruptions from a single magmatic chamber [9- 10]. The occurrence of alkali minerals like sanidine and K-rich rocks throughout Bradbury supports the presence of a potassic component, likely trachytic, while plagioclase and a mafic composition suggest a basaltic component [8-9]. Instruments: Mineral chemistry can be estimated by three instruments onboard Curiosity. The CheMin instrument enables detection of mineral assemblages using X-ray diffraction (XRD). Using Rietveld refinement, each mineral is identified according to their 1D XRD pattern [11]. Note that distinction between pyroxene minerals is challenging with the CheMin instrument due to overlapping peaks on XRD patterns and low angular resolution of the instrument [12]. Then, using least square regression and optimization algorithms based on unit-cell parameters, mineral chemistry has been estimated by [11]. Plagioclase composition has been estimated using the NaAlSi3O8- CaAl2Si2O8 system and alkali feldspar is based on the NaAlSi3O8-KAlSi3O8 system (stars in Fig.1). Two mudstone samples (John Klein and Cumberland) and one sandstone sample (Windjana) were analyzed by CheMin at Yellowknife Bay and Kimberley, respectively. Figure 1. Ternary diagrams of feldspars (top) and pyroxene (bottom) quadrilateral. Stars correspond to CheMin composition and the gray patches to ChemCam composition. The colored dots are the composition of feldspar and pyroxene that crystallized during fractional crystallization at FMQ+1 of a melt extracted at distinct melting degree during the adiabatic ascent of a primitive mantle composition, without any water (left panels) and with 0.5 wt.% of water (right panels) at distinct pressure. The ChemCam instrument enables the analysis of the chemical compositions of rocks at hundreds of micrometer scale (350-550 μm) using laser induced breakdown spectroscopy (LIBS), which may provide the composition of minerals when they are larger than the beam spot (>550 μm) [13]. Within >5000 LIBS points, we performed a typical stoichiometric filtering allowing us to distinguish 56 feldspar and 10 pyroxene mineral compositions (grey patches in Fig. 1). Finally, the Alpha Particle X-ray Spectrometer (APXS) analyzes the composition of rocks with a 1.6 cm diameter spot size. Monte Carlo mass balance modeling allowed [3] to decipher a feldspar range varying between An30 and An40 (Fig. 1). Discussion: Although there could be a more complex history and other ways to form the whole compositional range of igneous minerals analyzed within the Bradbury formation, we are presenting here simple magmatic pathways commonly occurring on Earth using the thermodynamical softwares pMELTS and rhyoliteMELTS [14]. The objective is to find reasonable igneous processes that produce minerals that parallel the compositions of feldspar and pyroxene analyzed by the Curiosity rover. As commonly observed for mid-ocean ridge basalts, the adiabatic ascent of a primitive mantle composition [15] partially melting at 2 GPa has been modeled, followed by the extraction of a liquid at distinct degrees of partial melting, which undergoes fractional crystallization at an oxygen fugacity +1 log unit above the fayalite-magnetite-quartz (FMQ) buffer within the crust (0.02-0.4 GPa) with H2O = 0- 0.5 wt. %. These latter conditions correspond to those recorded within igneous clasts from the Noachian martian breccia NWA 7034 and paired and within Gale igneous rocks (colored dots in Fig. 1) [16-17]. To check the reliability of these 2-step models, we also tested fractional crystallization at similar conditions (FMQ+1; P=0.02-0.4 GPa; H2O = 0-0.5 wt. %) of starting compositions corresponding to that of magmas with distinct melting degrees obtained from isobaric experiments at 2 GPa [18]. Mineral compositions obtained from both models are similar. As shown on Fig. 1, the whole range of observed feldspar compositions cannot be reproduced by fractionation of one magma only. Indeed, while alkali feldspar and Na-plagioclase likely crystallized from fractional crystallization of a low-degree melt (here <15%), plagioclase and pyroxene can only be formed by fractional crystallization of a higher degree melt (here >19%). The corresponding liquid descent lines are broadly in agreement with compositions estimated by ChemCam corresponding to float igneous rocks (Fig. 2) [4-6]. Figure 2. Silica versus alkali content. Lines show the liquid lines of descent from magmas with distinct degrees of melting. Gray patches represents the composition of Gale igneous rocks [4-6]. Trachytic to rhyolitic magmas crystallize alkali feldspar, and andesite to dacite magmas likely form plagioclase. Therefore, at least two starting magmas at distinct melting degrees, which could easily come from a single mantle source, are necessary to explain the whole compositional range of feldspar and pyroxene analyzed within Bradbury rocks. Conclusion: Because rocks from the Bradbury formation are likely originating from the same magmatic source with minimal weathering as supported by several studies using different approaches, igneous mineral chemistry analyzed by CheMin and ChemCam allows us to back out reasonable magmatic pathways that could have crystallized them. Fractional crystallization of at least two starting magmas originating from distinct melting degrees of a single mantle source can explain the whole range of feldspar and pyroxene composition. Both alkaline and sub-alkaline liquids can be produced, with compositions corresponding to those of the igneous rocks analyzed by ChemCam within the Bradbury formation, highlighting the complexity of Mars magmatism.

Payre, V.↗
DOE OSTI2024

Daytime cooling efficiencies of urban trees derived from land surface temperature are much higher than those for air temperature

Accurately capturing the impact of urban trees on temperature can help optimize urban heat mitigation strategies. Recently, there has been widespread use of remotely sensed land surface temperature (T s ) to quantify the cooling efficiency (CE) of urban trees. However, remotely sensed T s reflects emitted radiation from the surface of an object seen from the point of view of the thermal sensor, which is not a good proxy for the air temperature (T a ) perceived by humans. The extent to which the CEs derived from T s reflect the true experiences of urban residents is debatable. Therefore, this study systematically compared the T s -based CE (CE T s ) with the T a -based CE (CE T a ) in 392 European urban clusters. CE T s and CE T a were defined as the reductions in T s and T a , respectively, for every 1% increase in fractional tree cover (FTC). The results show that the increase in FTC has a substantial impact on reducing T s and T a in most cities during daytime. However, at night, the response of T s and T a to increased FTC appears to be much weaker and ambiguous. On average, for European cities, daytime CE T s reaches 0.075 °C % -1 , which is significantly higher (by an order of magnitude) than the corresponding CE T a of 0.006 °C % -1 . In contrast, the average nighttime CE T s and CE T a for European cities are similar, both approximating zero. Overall, urban trees can lower daytime temperatures, but the magnitude of their cooling effect is notably amplified when using remotely sensed T s estimates compared to in situ T a measurements, which is important to consider for accurately constraining public health benefits. Our findings provide critical insights into the realistic efficiencies of alleviating urban heat through tree planting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Planck 2018 results

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E -mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: f NL local = -0.9 ± 5.1; f NL equil = -26 ± 47; and f NL ortho = -38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarization-only bispectra display a significant improvement in robustness; they can now be used independently to set primordial NG constraints with a sensitivity comparable to WMAP temperature-based results and they give excellent agreement. In addition to the analysis of the standard local, equilateral, and orthogonal bispectrum shapes, we consider a large number of additional cases, such as scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is, however, detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 σ . Beyond estimates of individual shape amplitudes, we also present model-independent reconstructions and analyses of the Planck CMB bispectrum. Our final constraint on the local primordial trispectrum shape is g NL local = (-5.8 ± 6.5) × 10 4 (68% CL, statistical), while constraints for other trispectrum shapes are also determined. Exploiting the tight limits on various bispectrum and trispectrum shapes, we constrain the parameter space of different early-Universe scenarios that generate primordial NG, including general single-field models of inflation, multi-field models (e.g. curvaton models), models of inflation with axion fields producing parity-violation bispectra in the tensor sector, and inflationary models involving vector-like fields with directionally-dependent bispectra. Our results provide a high-precision test for structure-formation scenarios, showing complete agreement with the basic picture of the ΛCDM cosmology regarding the statistics of the initial conditions, with cosmic structures arising from adiabatic, passive, Gaussian, and primordial seed perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Evaluating Global Emission Inventories of Biogenic Bromocarbons

Emissions of halogenated very short-lived substances (VSLS) are poorly constrained. However, their inclusion in global models is required to simulate a realistic inorganic bromine (Bry) loading in both the troposphere, where bromine chemistry perturbs global oxidizing capacity, and in the stratosphere, where it is a major sink for ozone (O3). We have performed simulations using a 3-D chemical transport model (CTM) including three top-down and a single bottom-up derived emission inventory of the major brominated VSLS bromoform (CHBr3) and dibromomethane (CH2Br2). We perform the first concerted evaluation of these inventories, comparing both the magnitude and spatial distribution of emissions. For a quantitative evaluation of each inventory, model output is compared with independent long-term observations at National Oceanic and Atmospheric Administration (NOAA) ground-based stations and with aircraft observations made during the NSF (National Science Foundation) HIAPER Pole-to-Pole Observations (HIPPO) project. For CHBr3, the mean absolute deviation between model and surface observation ranges from 0.22 (38 %) to 0.78 (115 %) parts per trillion (ppt) in the tropics, depending on emission inventory. For CH2Br2, the range is 0.17 (24 %) to 1.25 (167 %) ppt. We also use aircraft observations made during the 2011 Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere (SHIVA) campaign, in the tropical western Pacific. Here, the performance of the various inventories also varies significantly, but overall the CTM is able to reproduce observed CHBr3 well in the free troposphere using an inventory based on observed sea-to-air fluxes. Finally, we identify the range of uncertainty associated with these VSLS emission inventories on stratospheric bromine loading due to VSLS (Br(VSLS/y)). Our simulations show Br(VSLS/y) ranges from approximately 4.0 to 8.0 ppt depending on the inventory. We report an optimized estimate at the lower end of this range (approximately 4 ppt) based on combining the CHBr3 and CH2Br2 inventories which give best agreement with the compilation of observations in the tropics.

bromocarbons↗
DOE OSTI2020

A Custom Machine to Convey Radiologically Impacted Cohesive Soils to the Orion ScanSortSM System - 20382

An established method to reduce the volume of material required to be disposed of as radioactive waste is conveyor-based sorting of potentially radiologically impacted soils. With the Orion ScanSortSM system, potentially impacted soils are conveyed beneath radiation detectors and sorted into above and below criteria bins based on the detectors' response. Confident measurements can only be efficiently achieved when the soil column being conveyed has a reasonably consistent geometry. Highly cohesive soils present a very significant challenge as they tend to clump and adhere to the conveyor hopper, strike-off bar, and belt skirting. This causes voids and valleys in the soil column that affect the measurement geometry and reduce the confidence in the measurement. The reduction in confidence necessitates a longer residence time (i.e. reduction in conveyor belt speed and processing rates) or that the soils in highly unfavorable geometries be dispositioned as impacted or segregated for resurvey (should site logistics support that option). These outcomes may increase the volume of material required to be dispositioned as radiological waste and have a negative impact on project cost and schedule. To minimize impact to sorting operations, Wood designed and built a customized Extruder, to optimize the soil column prior to assay. The design was developed and refined to minimize the volume of highly cohesive soils presented in highly unfavorable geometries in the soil column in order to confidently assess and disposition such soils in a high-production environment. The Extruder is a customized 90 cm wide flat conveyor that has an oversized hopper with a pair of motor-driven rollers that force the soils in the hopper through an opening of adjustable height. The soils are extruded through the opening to produce a soil column that is ∼75 cm wide and ∼8-18 cm deep with a design rate of 225 metric tons/hr. The Extruder was recently deployed with the Orion ScanSort{sup SM} System to a site in northeast Ohio. The site has highly cohesive soils that result from the weathering of glacial sediments and consist of clay and clay loam, resulting in poor drainage and high moisture content. The project infrastructure only supported two bins of material (above and below criteria), which necessitated that soils with unfavorable geometries be discharged into the above criteria bin, along with material determined to have been contaminated. The Extruder was extremely successful in producing a stable soil column with favorable geometries for radiological assay. During production, less than 0.5% (by mass) of the soils processed were presented in unsatisfactory geometries requiring disposition as radioactive waste. The Extruder did not get clogged due to the cohesive soils (as is typical with conventional conveyors) and required very little maintenance, again minimizing impacts to sorting operations. The conveyor belt speed was operated at 14 cm/s with a typical belt loading of ∼109 kg/m resulting in a mean process rate of ∼49 metric tons/hr. The process rate was constrained by other project logistics, rather than by the capabilities of the Extruder, which was operated at less than 20% of the maximum design speed (75 cm/s). It is concluded that the Extruder is highly likely to produce a stable highly cohesive soil column suitable for efficient and confident radiological assay in production environments of 250 short tons/hr or more. The Extruder has the capability to drastically reduce the duration and costs of projects with radiologically impacted highly cohesive soils. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗