Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fast function extraction”

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 55 records · Page 3

Greenhouse Gas Concentration Data Recovery Algorithm for a Low Cost, Laser Heterodyne Radiometer

The goal of a coordinated effort between groups at GWU and NASA GSFC is the development of a low-cost, global, surface instrument network that continuously monitors three key carbon cycle gases in the atmospheric column: carbon dioxide (CO2), methane (CH4), carbon monoxide (CO), as well as oxygen (O2) for atmospheric pressure profiles. The network will implement a low-cost, miniaturized, laser heterodyne radiometer (mini-LHR) that has recently been developed at NASA Goddard Space Flight Center. This mini-LHR is designed to operate in tandem with the passive aerosol sensor currently used in AERONET (a well established network of more than 450 ground aerosol monitoring instruments worldwide), and could be rapidly deployed into this established global network. Laser heterodyne radiometry is a well-established technique for detecting weak signals that was adapted from radio receiver technology. Here, a weak light signal, that has undergone absorption by atmospheric components, is mixed with light from a distributed feedback (DFB) telecommunications laser on a single-mode optical fiber. The RF component of the signal is detected on a fast photoreceiver. Scanning the laser through an absorption feature in the infrared, results in a scanned heterodyne signal io the RF. Deconvolution of this signal through the retrieval algorithm allows for the extraction of altitude contributions to the column signal. The retrieval algorithm is based on a spectral simulation program, SpecSyn, developed at GWU for high-resolution infrared spectroscopies. Variations io pressure, temperature, composition, and refractive index through the atmosphere; that are all functions of latitude, longitude, time of day, altitude, etc.; are modeled using algorithms developed in the MODTRAN program developed in part by the US Air Force Research Laboratory. In these calculations the atmosphere is modeled as a series of spherically symmetric shells with boundaries specified at defined altitudes. Temperature, pressure, and species mixing ratios are defined at these boundaries. Between the boundaries, temperature is assumed to vary linearly with altitude while pressure (and thus gas density) vary exponentially. The observed spectrum at the LHR instrument will be the integration of the contributions along this light path. For any absorption measurement the signal at a particular spectral frequency is a linear combination of spectral line contributions from several species. For each species that might absorb in a spectral region, we have pre-calculated its contribution as a function of temperature and pressure. The integrated path absorption spectrum can then by calculated using the initial sun angle (from location, date, and time) and assumptions about pressure and temperature profiles from an atmospheric model. The modeled spectrum is iterated to match the experimental observation using standard multilinear regression techniques. In addition to the layer concentrations, the numerical technique also provides uncertainty estimates for these quantities as well as dependencies on assumptions inherent in the atmospheric models.

Miller, J. Houston↗

Genomic Signatures for Sedimentary Microbial Utilization of Phytoplankton Detritus in a Fast-Flowing Estuary

In fast-flowing, river-dominated estuaries, “hotspots” of microbial biogeochemical cycling can be found within areas of extended water retention. Lateral bays located off of the North and South channels of the Columbia River estuary are proposed to be such hotspots. Previous metagenomic studies on water samples indicated that these regions function both as sources and sinks of biogenic particles, with potential to impact organic matter fluxes in the estuary. To extend this work, we analyzed 11 sediment metagenomes from three disparate bays: the freshwater Cathlamet Bay, and the brackish Youngs Bay and more saline Baker Bay located nearer the mouth to the south and north of the main channel, respectively. Samples were collected from upper layers of sediments in August of 2011 and 2013 for DNA extraction and metagenome sequencing. All metagenomes were dominated by bacterial sequences, although diatom sequences as high as 26% of the total annotated sequences were observed in the higher salinity samples. Unsupervised 2D hierarchical clustering analysis resulted in the eleven metagenome samples clustered into four groups by microbial taxonomic composition, with Bacteroides, diatom, and phage levels driving most of the grouping. Results of functional gene clustering further indicated that diatom bloom degradation stage (early vs. late) was an important factor. While the Flavobacteriia and Cytophagia classes were well represented in metagenomes containing abundant diatoms, taxa from the Bacteroidia class, along with certain members of the Sphingobacteriia class, were particularly abundant in metagenomes representing later stages of diatom decomposition. In contrast, the sediment metagenomes with low relative abundance of diatom and Bacteroidetes sequences appeared to have a metabolic potential biased toward microbial growth under nutrient limitation. While differences in water salinity clearly also influenced the microbial community composition and metabolic potential, our results highlight a central role for allochthonous labile organic matter (i.e., diatom detritus), in shaping bacterial taxonomic and functional properties in the Columbia River estuary lateral bay sediments. These results suggest that in fast-flowing, river-dominated estuaries, sediment microbial communities in areas of extended water retention, such as the lateral bays, may contribute disproportionately to estuarine organic matter degradation and recycling.

59 BASIC BIOLOGICAL SCIENCES↗

A deep learning-based workflow for fast prediction of 3D state variables in geological carbon storage: A dimension reduction approach

Deep learning (DL) models are extensively used as surrogate models for high-fidelity simulations of multiphase fluid flow in porous media at large scales, enabling fast forecasts of the spatial–temporal evolution of three-dimensional (3D) state variables in geological carbon storage (GCS). However, training these models in high-dimensional space remains computationally demanding and prone to overfitting because of limited training data. This paper presents a novel workflow to address these challenges by integrating dimension reduction (DR) methods. Here, the proposed workflow employed pre-trained DR models to extract the latent variables of geological models and state variables and utilized the multi-layer perceptron (MLP) for constructing mapping functions between the input and output variables in latent spaces. Subsequently, the pre-trained reconstruction models converted the MLP-predicted latent state variables to their original high-dimensional form. Furthermore, we proposed a novel strategy for the DR and reconstruction of 3D saturation fields to account for the unique data characteristics of sparsity, nonuniformity, and discontinuity. The proposed strategy applied PCA and inverse PCA for 2D average saturation fields and developed a DL-based 3D reconstruction model, leveraging three 2D average saturation fields as input to produce a 3D saturation field as output. The pre-training of DR and reconstruction models and training of MLP models were conducted on 84 Gulf of Mexico (GoM) simulations and evaluated on 12 testing simulations. Each simulation contained 720 monthly time steps, with the first 360 months as the injection period and the rest as the post-injection period. The proposed workflow, incorporating DR and DL models, accurately predicts the normalized 3D pressure fields, achieving mean square error (MSE) of 2.92 × 10 -7 compared to the ground truth obtained from a full-physics simulator. Furthermore, the proposed strategy outperformed PCA and convolutional autoencoder (CAE) models on 3D saturation fields, resulting in minor workflow prediction errors with an MSE of 2.93 × 10 -5 . The results suggest the proposed workflow provides sufficient predictive fidelity across temporal and spatial scales, and enables a speedup of 160 times compared to the full-physics simulator, facilitating improved decision-making and risk assessment for large-scale GCS management in real-time scenarios.

3D reconstruction model↗

Efficient machine-learning model for fast assessment of elastic properties of high-entropy alloys

We combined descriptor-based analytical models for stiffness-matrix and elastic-moduli with mean-field methods to accelerate assessment of technologically useful properties of high-entropy alloys, such as strength and ductility. Model training for elastic properties uses Sure-Independence Screening (SIS) and Sparsifying Operator (SO) method yielding an optimal analytical model, constructed with meaningful atomic features to predict target properties. Computationally inexpensive analytical descriptors were trained using a database of elastic properties determined from density functional theory for binary and ternary subsets of Nb-Mo-Ta-W-V refractory alloys. The optimal Elastic-SISSO models, extracted from an exponentially large feature space, give an extremely accurate prediction of target properties, similar to or better than other models, with some verified from existing experiments. Here we also show that electronegativity variance and elastic-moduli can directly predict trends in ductility and yield strength of refractory HEAs, and reveals promising alloy concentration regions.

36 MATERIALS SCIENCE↗

Modeling the co-assembly of binary nanoparticles

Abstract In this work, we present a binary assembly model that can predict the co-assembly structure and spatial frequency spectra of monodispersed nanoparticles with two different particle sizes. The approach relies on an iterative algorithm based on geometric constraints, which can simulate the assembly patterns of particles with two distinct diameters, size distributions, and at various mixture ratios on a planar surface. The two-dimensional spatial-frequency spectra of the modeled assembles can be analyzed using fast Fourier transform analysis to examine their frequency content. The simulated co-assembly structures and spectra are compared with assembled nanoparticles fabricated using transfer coating method are in qualitative agreement with the experimental results. The co-assembly model can also be used to predict the peak spatial frequency and the full-width at half-maximum bandwidth, which can lead to the design of the structure spectra by selection of different monodispersed particles. This work can find applications in fabrication of non-periodic nanostructures for functional surfaces, light extraction structures, and broadband nanophotonics.

Mohanty, Saurav (ORCID:0009000302592855)↗

Ionic Control over Ferroelectricity in 2D Layered van der Waals Capacitors

The van der Waals layered material CuInP 2 S 6 features interesting functional behavior, including the existence of four uniaxial polarization states, polarization reversal against the electric field through Cu ion migration, a negative-capacitance regime, and reversible extraction of Cu ions. At the heart of these characteristics lies the high mobility of Cu ions, which also determines the spontaneous polarization. Therefore, Cu migration across the lattice results in unusual ferroelectric behavior. Here, we demonstrate how the interplay of polar and ionic properties provides a path to ionically controlled ferroelectric behavior, achieved by applying selected DC voltage pulses and subsequently probing ferroelectric switching during fast triangular voltage sweeps. Using current measurements and theoretical calculations, we observe that increasing DC pulse duration results in higher ionic currents, the buildup of an internal electric field that shifts polarization loops, and an increase in total switchable polarization by ~50% due to the existence of a high polarization phase which is stabilized by the internal electric field. Apart from tuning ferroelectric behavior by selected square pulses, hysteretic polarization switching can even be entirely deactivated and reactivated, resulting in three-state systems where polarization switching is either inhibited or can be performed in two different directions.

36 MATERIALS SCIENCE↗

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗

Graph Contractions for Calculating Correlation Functions in Lattice QCD

Computing correlation functions for many-particle systems in Lattice QCD is vital to extract nuclear physics observables like the energy spectrum of hadrons such as protons. However, this type of calculation has long been considered to be very challenging and computing-resource intensive because of the complex nature of a hadron composed of quarks with many degrees of freedom. In particular, a correlation function can be calculated through a sum of all possible pairs of quark contractions, each of which is a batched tensor contraction, dictated by Wick's theorem. Because the number of terms of this sum can be very large for any hadronic system of interest, fast evaluation of the sum faces several challenges: an extremely large number of contractions, a huge memory footprint at runtime, and the speed of tensor contractions. In this paper, we present a Lattice QCD analysis software suite, Redstar, which addresses these challenges by utilizing novel algorithmic and software engineering methods targeting modern computing platforms such as many-core CPUs and GPUs. In particular, Redstar represents every term in the sum of a correlation function by a graph, applies efficient graph algorithms to reduce the number of contractions to lower the cost of computations, and minimizes the total memory footprint. Moreover, Redstar carries out the contractions on either CPUs or GPUs utilizing an internal and highly efficient Hadron contraction library. Specifically, we illustrate some important algorithmic optimizations of Redstar, show various key design features of Hadron library, and present the speedup values due to the optimizations along with performance figures for calculating six correlations functions on four computing platforms.

Chen, Jie↗

Problem areas in the use of the firefly luciferase assay for bacterial detection

By purifying the firefly luciferase extract and adding all necessary chemicals but ATP in excess, an assay for ATP was performed by measuring the amount of light produced when a sample containing soluble ATP is added to the luciferase reaction mixture. Instrumentation, applications, and basic characteristics of the luciferase assay are presented. Effect of the growth medium and length of time grown in this medium on ATP per viable E. coli values is shown in graphic form, along with an ATP concentration curve showing relative light units versus ATP injected. Reagent functions and concentration methods are explored. Efforts to develop a fast automatable system to detect the presence of bacteria in biological fluids, especially urine, resulted in the optimization of procedures for use with different types of samples.

Picciolo, G. L.↗

Laboratory modeling and analysis of aircraft-lightning interactions

Modeling studies of the interaction of a delta wing aircraft with direct lightning strikes were carried out using an approximate scale model of an F-106B. The model, which is three feet in length, is subjected to direct injection of fast current pulses supplied by wires, which simulate the lightning channel and are attached at various locations on the model. Measurements are made of the resulting transient electromagnetic fields using time derivative sensors. The sensor outputs are sampled and digitized by computer. The noise level is reduced by averaging the sensor output from ten input pulses at each sample time. Computer analysis of the measured fields includes Fourier transformation and the computation of transfer functions for the model. Prony analysis is also used to determine the natural frequencies of the model. Comparisons of model natural frequencies extracted by Prony analysis with those for in flight direct strike data usually show lower damping in the in flight case. This is indicative of either a lightning channel with a higher impedance than the wires on the model, only one attachment point, or short streamers instead of a long channel.

Turner, C. D.↗

Dimensionally Aligned Signal Projection Algorithms Library

Dimensionally aligned signal projection (DASP) algorithms are used to analyze fast Fourier transforms (FFTs) and generate visualizations that help focus on the harmonics for specific signals. At a high level, these algorithms extract the FFT segments around each harmonic frequency center, and then align them in equally sized arrays ordered by increasing distance from the base frequency. This allows for a focused view of the harmonic frequencies, which, among other use cases, can enable machine learning algorithms to more easily identify salient patterns. This work seeks to provide an effective open-source implementation of the DASP algorithms proposed by Vann et al. (2018) as well as functionality to help explore and test how these algorithms work with an interactive dashboard and signal-generation tool. The DASP library is implemented in Python and contains four types of algorithms for implementing these feature engineering techniques: fixed harmonically aligned signal projection (HASP), decimating HASP, interpolating HASP, and frequency aligned signal projection (FASP). Each algorithm returns a numerical array, which can be visualized as an image. The HASP algorithms are variations of the algorithms originally presented by Vann et al. (2018). For consistency, FASP, which is the terminology used for the short-time Fourier transform (STFT), has been implemented as part of the library to provide a similar interface to the STFT of the raw signal. Additionally, the library contains an algorithm to generate artificial signals with basic customizations such as the base frequency, sample rate, duration, number of harmonics, noise, and number of signals. Finally, the library provides multiple interactive visualizations, each of which is implemented using IPyWidgets and works in a Jupyter environment. A dashboard-style visualization is provided, which contains some common signal-processing visual components (signal, FFT, spectogram) updating in unison with the HASP functions (see Figure 1 below). Separate from the dashboard, an independent visualization is provided for each of the DASP algorithms as well as for the artifical signal generator. These visualizations are included in the library to aid in developing an intuitive understanding how the algorithms are affected by different input signals and parameter selections.

harmonics↗

Determination of Molecular Structure and Dynamics of Molten Salts by Advanced Neutron and X-ray Scattering Measurements and Computer Modeling

The design and development of fully functional Molten-Salt Reactors (MSR) require detailed knowledge of the molten salt properties in order to understand and predict the salt’s behavior. Fundamental properties of interest include molecular structure, speciation, and dynamics (such as diffusion coefficients) of salt components and dissolved corrosion and fission products. Computer modeling is necessary to predict changes in physical and chemical properties due to irradiation, burning of dissolved fuel, and corrosion. The modeling requires experimental data, and advanced neutron and x-ray scattering and spectroscopy provide the most reliable and direct determination of the structure (Pair-Distribution Functions, PDF), and dynamics of ions in the melt. This project dealt with both fluoride and chloride salts. The PDFs have been measured by a combination of neutron and x-ray diffraction. We utilized the techniques of isotope substitutions, a very powerful tool available for neutron-scattering, to extract the details of the liquid structure. Although similar measurements have been done before, modern advanced neutron and x-ray-scattering techniques allow collecting the data at much higher resolution and in a wider range of temperatures. Importantly, we were among the first to study fluoride salts by neutron scattering. The importance of impurities and their effects on salt properties have become apparent recently and so new methods of salt purification were developed. We took advantage of these developments to produce reliable data, which have been used for computer simulations of both clean salts and those with added fission and corrosion products most relevant for MSRs. Ab initio molecular dynamics simulations have been performed to understand the multi-component liquid solution, in particular solubility of impurities and thermodynamic interactions in relation to the ionic-cluster structure of the fluid. We applied machine learning to regress from the simulation and experimental data in order to develop a fast-acting model that can handle molten salt with an arbitrary (≥ 10) number of chemical elements and be able to predict chemical potential as a function of composition and temperature. This project resulted in a number of experimental and computer-simulation publications, a patent application, and numerous conference presentations (American Physical Society, American Chemical Society, and The Electrochemical Society among others). Multiple students and postdocs participated and collaborated on aspects of this project. This project seeded new collaborations between MIT and other institutions, such as the University of Massachusetts Lowell, the University of Illinois Urbana-Champaign, the University of California Berkeley, and Oak Ridge and Los Alamos National Labs. As such, this project has had a broad and lasting impact beyond its original scientific scope.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fast HARDI Uncertainty Quantification and Visualization with Spherical Sampling

In this paper, we study uncertainty quantification and visualization of orientation distribution functions (ODF), which corresponds to the diffusion profile of high angular resolution diffusion imaging (HARDI) data. The shape inclusion probability (SIP) function is the state‐of‐the‐art method for capturing the uncertainty of ODF ensembles. The current method of computing the SIP function with a volumetric basis exhibits high computational and memory costs, which can be a bottleneck to integrating uncertainty into HARDI visualization techniques and tools. We propose a novel spherical sampling framework for faster computation of the SIP function with lower memory usage and increased accuracy. In particular, we propose direct extraction of SIP isosurfaces, which represent confidence intervals indicating spatial uncertainty of HARDI glyphs, by performing spherical sampling of ODFs. Our spherical sampling approach requires much less sampling than the state‐of‐the‐art volume sampling method, thus providing significantly enhanced performance, scalability, and the ability to perform implicit ray tracing. Our experiments demonstrate that the SIP isosurfaces extracted with our spherical sampling approach can achieve up to 8164× speedup, 37282× memory reduction, and 50.2% less SIP isosurface error compared to the classical volume sampling approach. We demonstrate the efficacy of our methods through experiments on synthetic and human‐brain HARDI datasets.

97 MATHEMATICS AND COMPUTING↗

Direct observation of C 3 S particle dissolution using fast nano X-ray computed tomography

Tricalcium silicate (C 3 S) occupies 50 % to 70 % of ordinary portland cement (OPC) by mass and it is an important component affecting the hydration of OPC [1], [2], [3], [4], [5], [6], [7]. Generally, the hydration of C 3 S is described by two processes: the dissolution of C 3 S particles and the precipitation of hydration products. While it is understood that the dissolution rates of C 3 S vary with time, more precise measurements are needed to understand this process. Many mechanisms have been proposed to explain the time-evolving dissolution rates of C 3 S [2]. The metastable barrier hypothesis suggests that a thin metastable layer of hydrates forms around the C 3 S particle surface and prohibits the access of grains to the aqueous solution [8], [9], [10], [11], [12]. The slow dissolution step hypothesis suggests that the increased ion concentration from the initial reaction delays the C 3 S dissolution [2], [13], [14], [15], [16], [17]. More recent publications suggest that C 3 S may react differently depending on the existence of crystallographic defects [18], [19], [20]. Etch pits are thought to open on the particle surface during the initial reaction and this contributes to the C 3 S dissolution [21], [22]. As hydrates precipitate and cover these highly reactive surfaces, hydration slows down and the induction period starts [18], [23], [24], [25], [26]. Many experiments have been conducted to test the aforementioned mechanisms. Some hydration studies utilize bulk measurements, such as isothermal calorimetry [27], [28], [29], pore solution analysis [30], quasi-elastic neutron scattering [31], and nuclear magnetic resonance spectroscopy [32], [33]. One limitation of these measurements is that they do not provide direct and detailed information on the individual C 3 S particles. Some other studies utilize imaging techniques, such as scanning electron microscopy (SEM) [34], [35], [36] and transmission electron microscopy (TEM) [37]. However, SEM/TEM cannot track the evolution of individual particles throughout hydration [34], [35], [38], [39], [40] and they do not give insights into the microstructure of materials before hydration [34], [35], [39]. This makes it challenging to draw strong conclusions from only SEM or TEM observations. Synchrotron X-ray tomography techniques have been used more broadly in recent years to study cement hydration. They are not only non-destructive but also able to image a sample in full 3D with resolutions that can reach from micron to nanoscale. Nano computed tomography (nCT) is one technique that has been applied to study cement hydration at the nanoscale [26], [41]. A typical nCT can reach a pixel size from 15 to 65 nm, providing enough detail for observing features <1 μm. However, nCT often takes >0.5 h to finish one scan. This makes the application of this technique on continuous scans for in-situ observations challenging. Fast X-ray computerized tomography (fCT) is another technique that has shown success in studying the time-evolving cement microstructures [20], [42], [43], [44], [45], [46], [47], [48], [49]. Due to the high flux of the X-ray beam from the synchrotron ring, fCT allows a scan to be captured within 1 min at a pixel size of 1 μm. This allows a paste sample to be continuously scanned during the hydration process. However, the micron-sized resolutions limit does not provide detailed insights for particles <5 μm [20], [49]. Fortunately, the combination of nCT and fCT has allowed the development of fast nano X-ray computed tomography (fnCT). fnCT can capture a 3D data set in <2 min at a pixel size of 50 nm. This makes this procedure an exciting method to evaluate hydrating pastes. fnCT collects multiple X-ray radiographs at various rotation angles and generates a 3D model of the scanned sample, which is also referred to as a 3D tomography [50], [51]. In one tomography, the X-ray absorptions of different components (e.g., C 3 S and hydrates) differ as functions of density and chemistry [52], [53]. These X-ray absorption contrasts can be used to extract detailed information about the 3D microstructure [26], [54], [55]. In this paper, fnCT is used to collect time-lapse tomographs of hydrating C 3 S paste from 18 min after mixing to 7 h of hydration. The bulk measurements of anhydrous C 3 S, as well as the microstructural changes of individual C 3 S particles, are directly observed, quantified, and discussed. The dissolution behavior of C 3 S particles at various size scales is systematically analyzed and compared. This work aims to find the relationship between the size of C 3 S particle sizes and their dissolution rates. This provides significant insights into the early-age hydration of C 3 S on length and time scales not previously possible. Because of the magnitude of the data and the substantial amount of observations, this work will solely focus on the change in the anhydrous particles. Changes in the hydration products will be reported in future work.

42 ENGINEERING↗

On the analysis of signal peaks in pulse-height spectra

The estimation of the signal location and intensity of a peak in a pulse height spectrum is important for X- ray and γ-ray spectroscopy, charged-particle spectrometry, liquid chromatography, and many other subfields. However, both the ‘‘centroid’’ and ‘‘signal intensity’’ of a peak in a pulse-height spectrum are ill-defined quantities and different methods of analysis will yield different numerical results. Here, we apply three methods of analysis. Method A is based on simple count summation and is likely the technique most frequently applied in practice. The analysis is straightforward and fast, and does not involve any statistical modeling. We find that it provides reliable results only for high signal-to-noise data, but has severe limitations in all other cases. Method B employs a Bayesian model to extract signal counts and centroid from the measured total and background counts. The resulting values are derived from the respective posteriors and, therefore, have a rigorous statistical meaning. The method makes no assumptions about the peak shape. It yields reliable and relatively small centroid uncertainties. However, it provides relatively large signal count uncertainties. Method C makes a strong assumption regarding the peak shape by fitting a Gaussian function to the data. Furthermore, the fit is based again on a Bayesian model. Although Method C requires careful consideration of the Gaussian width (usually given by the detector resolution) used in the fitting, it provides reliable values and relatively small uncertainties both for the signal counts and the centroid.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Scattering phase shift in quantum mechanics on quantum computers

Here, we investigate the feasibility of extracting infinite volume scattering phase shift on quantum computers in a simple one-dimensional quantum mechanical model, using the formalism established in the work by Guo and Gasparian [Phys. Rev. D 108, 074504 (2023)] that relates the integrated correlation functions for a trapped system to the infinite volume scattering phase shifts through a weighted integral. The system is first discretized in a finite box with periodic boundary conditions, and the formalism in real time is verified by employing a contact interaction potential with exact solutions. Quantum circuits are then designed and constructed to implement the formalism on current quantum computing architectures. To overcome the fast oscillatory behavior of the integrated correlation functions in real-time simulation, different methods of postdata analysis are proposed and discussed. Test results on IBM hardware show that good agreement can be achieved with two qubits, but complete failure ensues with three qubits due to two-qubit gate operation errors and thermal relaxation errors.

Guo, Peng [Dakota State Univ., Madison, SD (United↗

On Fast Post-Processing of Global Positioning System Simulator Truth Data and Receiver Measurements and Solutions Data

Post-processing of data, related to a GPS receiver test in a GPS simulator and test facility, is an important step towards qualifying a receiver for space flight. Although the GPS simulator provides all the parameters needed to analyze a simulation, as well as excellent analysis tools on the simulator workstation, post-processing is not a GPS simulator or receiver function alone, and it must be planned as a separate pre-flight test program requirement. A GPS simulator is a critical resource, and it is desirable to move off the pertinent test data from the simulator as soon as a test is completed. The receiver and simulator databases are used to extract the test data files for postprocessing. These files are then usually moved from the simulator and receiver systems to a personal computer (PC) platform, where post-processing is done typically using PC-based commercial software languages and tools. Because of commercial software systems generality their functions are notoriously slow and more than often are the bottleneck even for short duration simulator-based tests. There is a need to do post-processing faster and within an hour after test completion, including all required operations on the simulator and receiver to prepare and move off the post-processing files. This is especially significant in order to use the previous test feedback for the next simulation setup or to run near back-to-back simulation scenarios. Solving the post-processing timing problem is critical for a pre-flight test program success. Towards this goal an approach was developed that allows to speed-up post-processing by an order of a magnitude. It is based on improving the post-processing bottleneck function algorithm using a priory information that is specific to a GPS simulation application and using only the necessary volume of truth data. The presented postprocessing scheme was used in support of a few successful space flight missions carrying GPS receivers.

Kizhner, Semion↗

Multilevel Analysis, Design, and Modeling of Coupling Advanced Nuclear Reactors and Thermal Energy Storage in an Integrated Energy System

This report discusses the different options for coupling thermal energy storage (TES) systems to advanced nuclear power plants (A-NPPs) in order to enable flexible and hybrid plant operation. An advanced light-water reactor (A LWR), a high-temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating a thermally balanced energy storage coupling design for thermal power extraction. The models presented herein showcase several design considerations, focusing on optimal deployment methodologies for achieving steady-state and transient-state operation with minimum disruption to the nuclear power cycle. This first part of the study presents steady-state models developed using Aspen HYSYS®, with the thermal energy bypass for NPP-TES coupling being varied at up to 50%. The various components were sized using the Aspen Process Economic Analyzer (APEA) and Aspen Exchanger Design and Rating (EDR), when applicable. Cost functions from these models were developed using the latest publicly available data obtained from APEA V11. The TES-coupled A-NPP steady-state models and cost functions then provided a baseline for dynamic operation and process optimization by using Idaho National Laboratory (INL)’s Framework for Optimization of Resources and Economics (FORCE) tools. A stochastic optimization of the various energy storage systems coupled to the A-NPPs was then performed using the Risk Analysis Virtual Environment (RAVEN) and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON). The signal processing and synthetic history capabilities of RAVEN were used to account for the unpredictable behavior of electricity markets. An autoregressive moving average (ARMA) model was used to analyze price signals from the Pennsylvania-New Jersey-Maryland (PJM) market and were applied to the HERON analysis in order to optimize a system with the best economics. Transient modeling evaluation was then performed using Modelica models within the HYBRID repository, which was developed at INL for the Department of Energy Integrated Energy Systems program for the characterization of dynamic integrated system behavior and feedback. This includes evaluation of the TES-coupled A-LWR systems’ impact on physical and thermal system response during imposed system demands. Additional TES-coupled reactor types, coupling approaches, markets, and TES technologies will be evaluated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗