Search NASA⌕ Search

SEARCH · Search NASA

Results for “flight computers”

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 163 records · Page 9

Nuclear Safety [Vol. 36, No. 2, July-December 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 195 The Chornobyl Accident Revisited, Part III: Chernobyl Source Term Release Dynamics and Reconstruction of Events During the Active Phase, A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 218 Second ANS Workshop on the Safety of Soviet-Designed Nuclear Power Plants, R. A. Bari; 234 Elements of a Nuclear Criticality Safety Program, C. M. Hopper; 243 Rickover, Excellence, and Criticality Safety Programs, R. E. Wilson; ACCIDENT ANALYSIS: 249 Transient Analysis of the PIUS Advanced Reactor Design with the TRAC-PF1/MOD2 Code, B. E. Boyack, J. L Steiner, S. C. Harmony, H. J. Stumpf, and J. F. Lime; 278 The Hierarchy-By-Interval Approach to Identifying Important Models that Need Improvement in Severe-Accident Simulation Codes, T. J. Heames, M. Khatib-Rahbar, J. E. Kelly, R. P. Jenks-Johnson, and Y.-S. Chen; 290 RELAP5/MOD3 Code Coupling Model, R. P. Martin; 299 Missiles Caused by Severe Pressurized-Water Reactor Accidents, R. Krieg; DESIGN FEATURES: 310 Validation of COMMIX with Westinghouse AP-600 PCCS Test Data, J. G. Sun, T. H. Chien, J. Ding, and W T. Sha; ENVIRONMENTAL EFFECTS: 321 Spent Nuclear Fuel Characterization for a Bounding Reference Assembly for the Receiving Basin for Off-Site Fuel, S. D. Kahook, R. L. Garrett, L. R. Canas, and M J. Beckum; OPERATING EXPERIENCES: 332 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 335 Reactor Coolant System Blowdown at Wolf Creek on September 17, 1994, J. V. Kauffman and S. L. Israel; RECENT DEVELOPMENTS: 344 Reports, Standards, and Safety Guides, D. S. Oueener; 349 Proposed Rule Changes as of June 30, 1995; ANNOUNCEMENTS: 320 Symposium on Acceptability of Risk From Radiation—Application to Manned Space Flight; 320 24th DOE/NRC Nuclear Air Cleaning and Treatment Conference; 361 Radiation Biology and Radiation Protection— Modern Developments and Tendencies in Radiation Biology; 362 1997 IEEE Sixth Conference on Human Factors and Power Plants; 354 The Authors; 360 Reviewers of Nuclear Safety, Vol. 36.

05 NUCLEAR FUELS↗

CHESS 2025: Leaf Area Index (LAI) for meadow, shrub, tree, and understory vegetation

This dataset contains Leaf Area Index (LAI) measurements made as part of the Colorado Headwaters Ecological Spectroscopy Study (CHESS) during June and July of 2025. Data were collected in the Upper Gunnison Basin, Colorado, across three study domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). Field observations of LAI were collected within 72 hours of airborne data collection by the National Ecological Observatory Network’s Aerial Observation Platform (NEON AOP). The NEON AOP collected waveform LiDAR (Light Detection and Ranging) and imaging spectrometer data in 426 spectral bands from the visible to shortwave infrared. LAI measurements were collected using the LICOR LAI-2200C Plant Canopy Analyzer following protocols outlined in the instrument manual (LI-COR 2019). Sampling targeted four distinct vegetation types: meadows, shrubs, trees, and aspen forest understory. We have archived data separately by site type because different field methods were used for each. At meadow sites, measurements were made at the four corners of 1m x 1m plots, with the instrument moving inward toward the center of the plot. At shrub sites, we measured the canopies of individual shrubs. At tree sites, we made measurements within a 10m x 10m subplot centered around a focal tree, with 30 observations taken on a regular grid. At aspen understory sites, we measured overstory trees following the tree protocol and understory herbaceous vegetation following the meadow protocol. All measurements included above-canopy (A) and below-canopy (B) readings, with specific protocols for scattering correction measurements in direct-sun conditions. Data were processed using the R package `rlai` (Worsham 2025). This package includes functions to calculate LAI, gap fraction, apparent clumping factor (Ω), scattering correction, and other canopy metrics. Package contents: Full file descriptions appear in ‘flmd.csv’. Files named according to the convention ‘lai_*_summary_data_cleaned.csv’ contain summary values of LAI, apparent clumping factor (Ωapp), and scattering correction factors for each site. These are the analysis-ready products that most data users will work with. Files named ‘lai_*_metadata_cleaned.csv’ contain additional site-level observations made during field collection. We have also archived intermediate and supplementary data for users who wish to check our processing approach or apply alternative methods. ‘raw_lai_2200C.zip’ contains the raw files as read from the LI-COR instrument, with no processing applied, in TXT format. The zip archive contains subdirectories by site type, which are further subdivided by sampling area. Filenames correspond to the sampling site number. ‘intermediate_results.zip’ contains detailed output from the processing routines, in JSON format. The zip archive contains subdirectories by site type; filenames correspond to the sampling site number. ‘scattering_correction_logs.zip’ contains logfiles from the implementation of Kobayashi et al.'s (2013) scattering correction algorithm. The logfiles report values of several parameters at each iteration of the algorithm, as the model converges toward a stable solution. They are intended for users who want to verify scattering correction performance. The zip archive contains subdirectories by site type; filenames correspond to the sampling site number. ‘spot_checks.csv’ reports LAI and other values for a small number of files processed with LI-COR FV2200 software (LI-COR 2013) using the same control parameters as in our R-based approach. Additional metadata are provided in a data dictionary describing column names and definitions (dd.csv), and in a file-level metadata file (flmd.csv). All zip files can be expanded with common archive utilities. TXT, CSV, and JSON files can be ingested into R or Python computing environments or read in common text editor utilities. Geospatial information: Geospatial data for mapping measurement site locations are in the files CHESS_polygons_lai_UTM.geojson, CHESS_polygons_shrub_UTM.geojson, and CHESS_polygons_meadow_UTM.geojson in the companion geospatial package for the 2025 CHESS campaign, ‘CHESS 2025: Location data for field observations and sampling’ (Henderson et al., 2026). CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. * Todorov and Worsham are co–first authors.

2018 NEON and 2025 CHESS Campaigns↗

Design of a multi-detector, single line-of-sight, time-of-flight system to measure time-resolved neutron energy spectra

In the dynamic environment of burning, thermonuclear deuterium–tritium plasmas, diagnosing the time-resolved neutron energy spectrum is of critical importance. Strategies exist for this diagnosis in magnetic confinement fusion plasmas, which presently have a lifetime of ~10 12 longer than inertial confinement fusion (ICF) plasmas. In this work, we present a novel concept for a simple, precise, and scale-able diagnostic to measure time-resolved neutron spectra in ICF plasmas. The concept leverages general tomographic reconstruction techniques adapted to time-of-flight parameter space, and then employs an updated Monte Carlo algorithm and National Ignition Facility-relevant constraints to reconstruct the time-evolving neutron energy spectrum. Reconstructed spectra of the primary 14.028 MeV n DT peak are in good agreement with the exact synthetic spectra. The technique is also used to reconstruct the time-evolving downscattered spectrum, although the present implementation shows significantly more error.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Status and Ambitions of the US Heavy Element Program

The aim of this whitepaper is to highlight the current capabilities and priorities of the US Heavy Element community and to provide the framework for a coordinated advancement of nuclear science from these studies. This is an organized effort to reflect on what has been achieved in the field given the recommendations and initiatives of the 2015 NSAC Long Range Plan, and on what can be realized in the next decade given current and possibly expanded investments. Current investments have positioned the US community to be among the world leaders in studies of the nuclear and chemical properties of the heaviest elements. These include studies of reaction mechanisms, moving us ever closer to the “island of stability”, in spectroscopy, allowing us to better understand nuclear structure at these extreme proton numbers, in chemical behavior, looking to determine how these elements should be placed on the Periodic Table, in performing the first measurements where isotopes are directly identified by their mass numbers, and in laying the foundation for a potential US-led new element discovery experiment. At present, the highest priority of the US Heavy Element community is to capitalize on the current investments by supporting the operations of US facilities at optimal values. These facilities include the Argonne Tandem Linac Accelerator System at Argonne National Laboratory and other Department of Energy facilities such as the 88-Inch Cyclotron at Lawrence Berkeley National Laboratory, which has a dedicated superheavy element program, as well as university laboratories, including Texas A&M University. The High Flux Isotope Reactor at Oak Ridge National Laboratory is crucial to providing the radioactive isotopes required for heavy element science targets. This facility should be supported to provide the actinide materials that are essential for US-based science. Production of stable, rare isotopes for beam material, including 48 Ca, 50 Ti, 54 Cr and 58 Fe, at the Stable Isotope Production and Research Center is critical to continued research in heavy element science and should be a priority. The continued development of targets for heavy element science and retaining US-based expertise is critical for the heavy element community. This is an area that is currently under pressure. For example, the target laboratory at Argonne National Laboratory serves a broad community and is currently under threat due to loss of critical personnel. The skills needed to make targets for nuclear science and develop new targetry methods need to be supported long term at Argonne, Oak Ridge and Lawrence Livermore National Laboratories as well as maintaining the programs at Oregon State University and San José State University as vital pipelines for training students. Advances in theory are the foundation to understand how nuclei behave and to predict those behaviors in new circumstances. Progress in these studies will necessitate continued and new investment and access to high-performance computing. The future health of the heavy element field is dependent on the continuous support of talented early-career professionals at all levels. It is critical that opportunities continue to be created for the next generation to become established in heavy element research so that we can ensure the field is attracting and retaining the best minds for continued success. It is also clear that to ensure diversity of ideas, perspectives and techniques, we need to recruit diverse personnel that are trained at the best facilities. The heavy element community is in support of continued investment to programs with initiatives in diversity, equity, and inclusion. Looking to the next decade of research, support needs to maintain and grow US leadership in heavy element science. Specifically, new investments in state-of-the-art instrumentation will be essential to scientific development of the field and in expanding scientific knowledge. Advances in the next generation of electron cyclotron resonance ion sources, multi-reflection time-of-flight devices, laser spectroscopy, trapping methods and next generation alpha and gamma spectroscopy systems should be prioritized.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An investigation into the potential for wind turbines to cause barotrauma in bats

The high rates of bat mortality caused by operating wind turbines is a concern for wind energy and wildlife stakeholders. One theory that explains the mortality is that bats are not only killed by impact trauma, but also by barotrauma that results from exposure to the pressure variations caused by rotating turbine blades. To date, no published research has calculated the pressure changes that bats may be exposed to when flying near wind turbines and then used these data to estimate the likelihood that turbines cause barotrauma in bats. To address this shortcoming, we performed computational fluid dynamics simulations of a wind turbine and analytical calculations of blade-tip vortices to estimate the characteristics of the sudden pressure changes bats may experience when flying near a utility-scale wind turbine. Because there are no data available that characterize the pressure changes that cause barotrauma in bats, we compared our results to changes in pressure levels that cause barotrauma and mortality in other mammals of similar size. This comparison shows that the magnitude of the low-pressures bats experience when flying near wind turbines is approximately 8 times smaller than the pressure that causes mortality in rats, the smallest mammal for which data are available. The magnitude of the high-pressures that bats may experience are approximately 80 times smaller than the exposure level that causes 50% mortality in mice, which have a body mass similar to several bat species that are killed by wind turbines. Further, our results show that for a bat to experience the largest possible magnitude of low- and high-pressures, they must take very specific and improbable flight paths that skim the surface of the blades. Even a small change in the flight path results in the bat being hit by the blade or experiencing a much smaller pressure change. Accordingly, if bats have a physiological response to rapid low- and high-pressure exposure that is similar to other mammals, we conclude that it is unlikely that barotrauma is responsible for a significant number of turbine-related bat fatalities, and that impact trauma is the likely cause of the majority of wind-turbine-related bat fatalities.

17 WIND ENERGY↗

A super-resolution technique to analyze single-crystal inelastic neutron scattering measurements using direct-geometry chopper spectrometers

Direct-geometry time-of-flight chopper neutron spectroscopy is instrumental in studying dynamics in liquid, powder, and single crystal systems. We report here that real-space techniques in optical imagery can be adapted to obtain reciprocal-space super resolution dispersion for phonon or magnetic excitations from single-crystal neutron spectroscopy measurements. The procedure to reconstruct super-resolution energy dispersion of excitations relies on an accurate determination of the momentum and energy-dependent point spread function and a dispersion correction technique inspired by an image disparity calculation technique commonly used in stereo imaging. Here by, applying these methods to spinwave dispersion data from a virtual neutron experiment demonstrates ~5-fold improvement over nominal energy resolution.

47 OTHER INSTRUMENTATION↗

Derivation and verification of the direct-sampling method for simulating Monte Carlo flight paths in tetrahedral meshes with linear finite-element cross sections

This paper provides a derivation of a direct-sampling approach for modeling continuously varying cross sections in tetrahedral-mesh-based Monte Carlo codes. Specifically, cross sections are spatially approximated using linear nodal finite elements. A linearization strategy is provided for non-linearly varying cross sections. The method is verified against seven analytical pure-absorber test problems. These test problems also highlight the benefit of using linear finite elements over element-wise-constant cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA’s Pandora SmallSat Mission: Simulating the Impact of Stellar Photospheric Heterogeneity and Its Correction

Stellar photospheric heterogeneity is a dominant astrophysical systematic impacting exoplanet transmission spectroscopy. NASA’s Pandora SmallSat Mission is designed to address this challenge through contemporaneous visible-band photometry and near-infrared spectroscopy of exoplanet host stars. Here, we present an end-to-end simulation study quantifying Pandora’s ability to infer stellar photospheric properties and correct stellar contamination using out-of-transit observations. We construct eight representative stellar activity scenarios and generate 160 simulated Pandora datasets, incorporating time-dependent stellar spectra, instrument response, and noise. Given accurate models, Bayesian retrievals of joint visible photometry (0.4–0.7 μm) and near-infrared spectroscopy (0.9–1.6 μm, R ≈ 120) recover photospheric temperatures with typical uncertainties of ≈30 K, with no significant bias. Models with two spectral components (i.e., a quiescent photosphere and spots) are strongly favored in 95% of cases; one-component models are preferred when true spot filling factors fall below a detection threshold of ≈0.3%. We propagate the true and inferred stellar parameters to compute true, inferred, and residual contamination signals under physically motivated spot geometries. For simple spot distributions, contamination signals of 10 2 –10 3 ppm are reduced to ≲10 ppm—well below Pandora’s expected transmission spectroscopy precision (30–100 ppm). For more complex spot distributions, geometric degeneracies limit deterministic corrections, leaving residual contamination at the 10 3 ppm level that must be mitigated using additional constraints, such as spot-crossing events and joint stellar–planetary retrievals of transmission spectra. These results define regimes in which stellar contamination can be corrected from stellar observations alone and show how Pandora stellar observations can identify cases where additional information is required.

Astronomy and AstroPhysics↗

Rapid Simulation of Unconventional Reservoirs by Multidomain Multiresolution Modeling Based on the Diffusive Time of Flight

Summary The fast marching method (FMM)-based rapid flow simulation has been shown to accelerate simulation efficiency by orders of magnitude by transforming 3D simulation to equivalent 1D simulation using the concept of the “diffusive time of flight” (DTOF). However, the 1D transformation does not directly apply to multiwell problems. In this paper, we propose a novel DTOF-based multidomain multiresolution discretization scheme to accelerate multiwell simulation of unconventional reservoirs. Our method formulates multiwell simulation problems based on the DTOF which displays the pressure front propagation in unconventional reservoirs. The DTOF contours are used to partition the reservoir into local and shared domains. A local domain is where the flow is dominated by a single well, and the shared domain is where the fluid flow is influenced by multiple wells. The DTOF contours expand independently in local domains and interfere in the shared domain. After the partitioning, each domain is discretized using a multiresolution scheme whereby the original 3D fine mesh is preserved near the wells to account for detailed physics including gravity, and the rest of the domain is discretized into 1D mesh based on the DTOF contours to alleviate the simulation workload. The power and efficacy of our approach are demonstrated using synthetic and field-scale simulation models with different degrees of geologic and well-completion complexity. The simulation results, number of active cells, and computation time for the proposed discretization scheme are compared with the original high-fidelity 3D model for each case. The results show that the proposed method is suitable for multiwell simulation problems in unconventional reservoirs and can accelerate flow simulations by orders of magnitude with minimal loss of accuracy. The novelty of this work is the creation of DTOF-derived multiresolution discretization with local and shared domains to simplify and accelerate the calculation of subsurface flow problems, especially in unconventional reservoirs. Our workflow can be easily interfaced with commercial simulators, making it suitable for large-scale field applications.

Engineering↗

Disjunctive linear separation conditions and mixed-integer formulations for aircraft conflict resolution

In this paper, we address the aircraft conflict resolution problem in air traffic control. We introduce new mixed-integer programming formulations for aircraft conflict resolution with speed, heading and altitude control which are based on disjunctive linear separation conditions. We first examine the two-dimensional aircraft conflict resolution problem with speed and heading control represented as continuous decision variables. We show that the proposed disjunctive linear separation conditions are equivalent to the classical nonlinear conditions for aircraft separation. Further, we characterise conflict-free trajectories based on aircraft velocity bounds and propose a simple pre-processing algorithm to identify aircraft pairs which are either always conflict-free, or which cannot be separated using speed and heading control only. We then incorporate altitude control and propose a lexicographic optimisation formulation that aims to minimise the number of flight level changes before resolving outstanding conflicts via two-dimensional velocity control. The proposed mixed-integer programming formulations are nonconvex, and we propose convex relaxations, decomposition methods and constraint generation algorithms to solve the two-dimensional and lexicographic optimisation formulations to guaranteed optimality. Numerical experiments on four types of conflict resolution benchmarking instances are conducted to test the performance of the proposed mixed-integer formulations. Further, the proposed method is compared against two benchmarks based on state-of-the-art approaches for the aircraft conflict resolution problem. Our numerical results show that the proposed method largely outperforms both benchmarks in terms of runtime and is able to solve significantly more instances to global optimality.

97 MATHEMATICS AND COMPUTING↗

Understanding Peelle’s Pertinent Puzzle bias in generalized least squares regression through eigenspectrum analysis

Certain correlation structures in the data covariance matrix (DCM) used for generalized least squares (GLS) regression can result in biased estimates, commonly known in the field of nuclear data evaluation as Peele’s Pertinent Puzzle (PPP). This article introduces a generative, forward modeling framework within which the PPP bias is characterized through an eigenspectrum analysis of the DCM. This analysis highlights the root cause of the bias, generalizes the problem beyond the nuclear data field, and provides insight to the problem regimes where it can occur. What follows is an understanding that the bias can show up for any experimental neutron time-of-flight data for which systematic uncertainties have been quantified. Lastly, a discussion of the adaptation of cross validation approaches that require pre-whitening to incorporate the known ‘fix’ to the PPP bias in the GLS estimator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Toward nanoscale molecular mass spectrometry imaging via physically constrained machine learning on co-registered multimodal data

Mass spectrometry imaging (MSI) plays a pivotal role in investigating the chemical nature of complex systems that underly our understanding in biology and medicine. Multiple fields of life science such as proteomics, lipidomics and metabolomics benefit from the ability to simultaneously identify molecules and pinpoint their distribution across a sample. However, achieving the necessary submicron spatial resolution to distinguish chemical differences between individual cells and generating intact molecular spectra is still a challenge with any single imaging approach. Here, we developed an approach that combines two MSI techniques, matrix-assisted laser desorption/ionization (MALDI) and time-of-flight secondary ion mass spectrometry (ToF-SIMS), one with low spatial resolution but intact molecular spectra and the other with nanometer spatial resolution but fragmented molecular signatures, to predict molecular MSI spectra with submicron spatial resolution. The known relationships between the two MSI channels of information are enforced via a physically constrained machine-learning approach and directly incorporated in the data processing. We demonstrate the robustness of this method by generating intact molecular MALDI-type spectra and chemical maps at ToF-SIMS resolution when imaging mouse brain thin tissue sections. This approach can be readily adopted for other types of bioimaging where physical relationships between methods have to be considered to boost the confidence in the reconstruction product.

97 MATHEMATICS AND COMPUTING↗

Implementation of a High-Mach Integral Boundary Layer Method for Arbitrary Streamlined Body Geometry

The Momentum-Energy Integral Technique (MEIT) is an integral boundary layer method for the high-Mach flow regime used to approximate heat transfer and viscous force quantities of interest along streamlines of an inviscid flow solution on the surface of a flight vehicle. This method allows rapid mid-fidelity estimation of these quantities which would otherwise require a much more expensive viscous flow solution to produce. Integral boundary layer methods like MEIT have been around for decades, though usually only formulated for simple geometries such as 2-dimensional wing shapes or axi-symmetric nose shapes. The implementation discussed herein has been generalized to apply to any 3-dimensional streamlined body geometry through correct treatment of the curvilinear axes (streamline attached) momentum and energy entrainment terms, and handling of arbitrary stagnation region geometry. This implementation is provided as a software package for the Python environment, along with readers for common inviscid flow solution providers such as NASA’s CART3D flow solver.

97 MATHEMATICS AND COMPUTING↗

Trajectory Optimization via Unsupervised Probabilistic Learning On Manifolds

This report investigates the use of unsupervised probabilistic learning techniques for the analysis of hypersonic trajectories. The algorithm first extracts the intrinsic structure in the data via a diffusion map approach. Using the diffusion coordinates on the graph of training samples, the probabilistic framework augments the original data with samples that are statistically consistent with the original set. The augmented samples are then used to construct conditional statistics that are ultimately assembled in a path-planing algorithm. In this framework the controls are determined stage by stage during the flight to adapt to changing mission objectives in real-time. A 3DOF model was employed to generate optimal hypersonic trajectories that comprise the training datasets. The diffusion map algorithm identfied that data resides on manifolds of much lower dimensionality compared to the high-dimensional state space that describes each trajectory. In addition to the path-planing worflow we also propose an algorithm that utilizes the diffusion map coordinates along the manifold to label and possibly remove outlier samples from the training data. This algorithm can be used to both identify edge cases for further analysis as well as to remove them from the training set to create a more robust set of samples to be used for the path-planing process.

42 ENGINEERING↗

Three-dimensional magnetohydrodynamic modeling of auto-magnetizing liner implosions on the Z accelerator

Auto-magnetizing (AutoMag) liners are cylindrical tubes that employ helical current flow to produce strong internal axial magnetic fields prior to radial implosion on ~100 ns timescales. AutoMag liners have demonstrated strong uncompressed axial magnetic field production (>100 T) and remarkable implosion uniformity during experiments on the 20 MA Z accelerator. However, both axial field production and implosion morphology require further optimization to support the use of AutoMag targets in magnetized liner inertial fusion (MagLIF) experiments. Data from experiments studying the initiation and evolution of dielectric flashover in AutoMag targets on the Mykonos accelerator have enabled the advancement of magnetohydrodynamic (MHD) modeling protocols used to simulate AutoMag liner implosions. Implementing these protocols using ALEGRA has improved the comparison of simulations to radiographic data. Specifically, both the liner in-flight aspect ratio and the observed width of the encapsulant-filled helical gaps during implosion in ALEGRA simulations agree more closely with radiography data compared to previous GORGON simulations. Although simulations fail to precisely reproduce the measured internal axial magnetic field production, improved agreement with radiography data inspired the evaluation of potential design improvements with newly developed modeling protocols. Three-dimensional MHD simulation studies focused on improving AutoMag target designs, specifically seeking to optimize the axial magnetic field production and enhance the cylindrical implosion uniformity for MagLIF. Importantly, by eliminating the driver current prepulse and reducing the initial inter-helix gap widths in AutoMag liners, simulations indicate that the optimal 30–50 T range of precompressed axial magnetic field for MagLIF on Z can be accomplished concurrently with improved cylindrical implosion uniformity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic simulations underestimate the effects of waves during magnetic reconnection

Collisionless plasma systems are often studied using fully kinetic simulations, where protons and electrons are treated as particles. Due to their computational expense, it is necessary to reduce the ion-to-electron mass ratio $m_i/m_e$ or the ratio between plasma and cyclotron frequencies in simulations of large systems. In this Letter we show that when electron-scale waves are present in larger-scale systems, numerical parameters affect their amplitudes and effects on the larger system. Using lower-hybrid drift waves during magnetic reconnection as an example, we find that the ratio between the wave electric field and the reconnection electric field scales as $\sqrt{m_i/m_e}$, while the phase relationship is also affected. The combination of these effects means that the anomalous drag that contributes to momentum balance in the reconnection region can be underestimated by an order of magnitude. The results are relevant to the coupling of electron-scale waves to ion-scale reconnection regions, and other systems such as collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantification of regional net CO 2 flux errors in the Orbiting Carbon Observatory-2 (OCO-2) v10 model intercomparison project (MIP) ensemble using airborne measurements

Inverse model intercomparison projects (MIPs) provide a chance to assess the uncertainties in inversion estimates arising from various sources. However, accurately quantifying ensemble CO 2 flux errors remains challenging and often relies on the ensemble spread. This study proposes a method for quantifying the errors in regional net surface–atmosphere CO 2 flux estimates from models taken from the Orbiting Carbon Observatory-2 (OCO-2) v10 MIP by using independent airborne CO 2 measurements for the period 2015–2017. We first calculate the root mean square error (RMSE) between the ensemble mean of posterior CO 2 concentrations and airborne observations and then isolate the CO 2 concentration errors caused solely by the ensemble mean of posterior net fluxes by subtracting the observation, representation, and transport errors from seven regions. Our analysis reveals that the flux errors projected onto CO 2 space account for 55 %–85 % of the regional average RMSE over the 3 years, ranging from 0.88 to 1.91 ppm. In five regions, the error estimates based on observations exceed those computed from the ensemble spread of posterior fluxes by a factor of 1.33–1.93, implying an underestimation of the actual flux errors, while their magnitudes are comparable in two regions. The adjoint sensitivity analysis identifies that the underestimation of flux errors is prominent where the magnitudes of fossil fuel emissions exceed those of terrestrial-biosphere fluxes by a factor of 3–31 over the 3 years. This suggests the presence of systematic biases in the inversion estimates associated with errors in the prescribed fossil fuel emissions common to all models. Our study emphasizes the value of airborne measurements for quantifying regional errors in ensemble net CO 2 flux estimates.

54 ENVIRONMENTAL SCIENCES↗