SEARCH · Engineering Papers
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.
This work was focused on developing novel methods to facilitate rapid data to knowledge transformation for subsurface applications, particularly those related to unconventional oil and gas.
We present an ongoing work in which a surrogate model is being developed to reproduce the response dynamics of the third-integer resonant extraction process in the Delivery Ring (DR) at Fermilab. This is in pursuit of smoothly extracting circulating beam to the Mu2e Experiment’s production target, whereby the goal is to extract a uniform slice of the circulating 1e12 protons in the DR over 25,000 turns (43 ms). The DR contains 3 harmonic sextupoles that excite a third-integer resonance and three fast, tune-ramping quadrupole magnets that drive the horizontal tune towards the 29/3 resonance. In our initial work, the surrogate model trains on a semi-analytical simulation provided in the same format as live data. Using Reinforcement Learning (and other potential ML methods), the trained surrogate acts as the “environment” in which a simple ML control agent could learn to dynamically adjust the quadrupole ramp at 430 break points within the 43 microsecond spill window. The controller will be hosted on a dedicated Arria 10 FPGA. In this work, we report the accuracy and fidelity of the surrogate model in comparison to the response dynamics of the physics simulator.
Plant cell wall polysaccharide analysis encompasses the utilization of a variety of analytical tools, including gas and liquid chromatography, mass spectrometry (MS), and nuclear magnetic resonance (NMR) spectroscopy. These methods provide complementary data, which enable confident structural proposals of the many complex polysaccharide structures that exist in the complex matrices of plant cell walls. However, cell walls contain fractions of varying solubilities, and a few techniques are available that can analyze all fractions simultaneously. We have discovered that permethylation affords the complete dissolution of both soluble and insoluble polysaccharide fractions of plant cell walls in organic solvents such as chloroform or acetonitrile, which can then be analyzed by a number of analytical techniques including MS and NMR. In this work, NMR structure analysis of 10 permethylated polysaccharide standards was undertaken to generate chemical shift data providing insights into spectral changes that result from permethylation of polysaccharide residues. This information is of especial relevance to the structure analysis of insoluble polysaccharide materials that otherwise are not easily investigated by solution-state NMR methodologies. The preassigned NMR chemical shift data is shown to be vital for NMR structure analysis of minor polysaccharide components of plant cell walls that are particularly difficult to assign by NMR correlation data alone. With the assigned chemical shift data, we analyzed the permethylated samples of destarched, alcohol-insoluble residues of switchgrass and poplar by two-dimensional NMR spectral profiling. Thus, we identified, in addition to the major polysaccharide components, two minor polysaccharides, namely, <5% 3-linked arabinoxylan (switchgrass) and <2% glucomannan (poplar). In particular, the position of the arabinose residue in the arabinoxylan of the switchgrass sample was confidently assigned based on chemical shift values, which are highly sensitive to local chemical environments. Furthermore, the high resolution afforded by the 1H NMR spectra of the permethylated switchgrass and poplar samples allowed facile relative quantitative analysis of their polysaccharide composition, utilizing only a few milligrams of the cell wall material. As a result, the concepts herein developed will thus facilitate NMR structure analysis of insoluble plant cell wall polysaccharides, more so of minor cell wall components that are especially challenging to analyze with current methods.
Here, in this paper, we describe a simple method for performing multifrequency eddy current characterization of free-standing uniform-thickness metallic foils using a forked inductive coil arrangement. The method involves measuring the mutual inductance between two coils when a foil is present between the coils, and when it is not present; the ratio of these mutual inductances is compared with an analytical solution, and foil conductivity, thickness, and sheet resistance are simultaneously estimated using numerical inversion and least-squares fitting. This method was used to characterize 34 non-ferrous metallic samples with thicknesses between 50 and 640 μm and with conductivities between 0.8 × 10 7 and 5.8 × 10 7 S/m. The estimated thicknesses from eddy current characterization agreed well with those measured using confocal optical techniques; the two approaches agreed to within 1 μm for samples that were thinner than 200 μm, and to within 0.5% for samples that had a thickness of 200 μm or greater. The estimated conductivities from eddy current characterization were in close agreement with expected values, given knowledge of the materials used. A particular strength of this approach is that the instrumentation needed is broadly available in research and development laboratories and the associated fixturing is easy to manufacture and assemble. A calibration procedure is described that can be used to reduce errors from geometric uncertainties. This calibration requires a sample that has only a known conductivity or thickness; both do not need to be known. The method described herein is likely extensible to conductivities and thickness well outside the ranges measured as part of this work.
In this study, a phase fault detection algorithm is developed by employing the envelope detector method. The proposed method diagnoses faults among phases and defines fault areas in the incoming signal. The designed algorithm consists of three steps: analytical signal conversion, complex magnitude, and fault detection. Initially, an analytical signal is obtained from the incoming power signal to determine the instantaneous amplitude and phase of the signal. A complex magnitude operation is applied to analytical signals to display changes in amplitude. On the basis of the threshold values specified by the user, the last step identifies the distortion signal in terms of the type of error and size. The proposed method is tested with realistically simulated substation power signal data and real power system data from the Grid Event Signature Library. The obtained results revealed that the proposed method detects distortions accurately.
Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab-initio wavefunction calculations at the multi-reference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron – nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on non-relativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538–549] and therefore had only limited applicability. Furthermore, results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex, demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.
The microgrid (MG) plays a crucial role in the energy transition, but its nonlinearity presents a significant challenge for large-signal power systems studies in the electromagnetic transient (EMT) time scale. In this paper, we develop a large-signal linear MG model that considers the detailed dynamics of the primary and zero-control levels based on the Koopman operator (KO) theory. Firstly, a set of observable functions is carefully designed to capture the nonlinear dynamics of the MG. The corresponding linear KO is then analytically derived based on these observables, resulting in the linear representation of the original nonlinear MG with observables as the new coordinate. The influence of external input on the system dynamics is also considered during the derivation, enabling control of the MG. We solve the voltage control problem using the traditional linear quadratic integrator (LQI) method to demonstrate that textbook linear control techniques can accurately control the original nonlinear MG via the developed KO-linearized MG model. Our proposed KO linearization method is generic and can be easily extended for different control objectives and MG structures using our analytical derivation procedure. Here, we validate the effectiveness of our methodology through various case studies.
Abstract The constant molar overflow (CMO) framework, while useful for shortcut distillation models, assumes that all components have the same latent heats of vaporization. A simple transformation, from molar flows to latent‐heat flows, allows shortcut models to retain the mathematical simplicity of the CMO framework while accounting for different latent heats, resulting in the constant heat transport (CHT) framework for adiabatic distillation columns. Although several past works have already proposed this transformation in the literature, it has not been well utilized in recent times. In this article, we show the utility of this transformation in upgrading various applications such as identifying energy‐efficient multicomponent distillation configurations based on heat duty rather than surrogate vapor flow. The method transforms the diagram to a diagram. Furthermore, we derive new and insightful analytical results in distillation, such as cumulative latent‐heat stage fractions having monotonic profiles within a distillation column under the CHT framework.
Site specific microstructure control is a critical research area within the field of additive manufacturing due to its potential to revolutionize part performance. One way to achieve site specific microstructure control is through control of the solidification conditions via the construction of intricate scan paths; however, the search space for such a problem is large. Previous attempts only considered the solidification conditions at the top surface while also requiring either lots of manual-fine tuning or large amounts of computational resources. This paper introduces a general method for scan path optimization which considers the solidification conditions in the bulk of the material without an increase in computational expense. This method consists of three core components:1. A heat transfer model for simulating the temperature field at a given time.2. A surrogate model which takes scan pattern information and temperature data and predicts the solidification conditions of the bulk as well as the meltpool depths for a spot melt.3. A decision algorithm to decide which spot melt should be printed next based on the outputs of the surrogate model.Each of these components can be changed without changing the overall method. Within this work, this method is applied in the creation of an algorithm containing a semi-analytic heat transfer model to simulate the temperature field, a fully convolutional neural network (FCNN) as the surrogate model, and a greedy decision algorithm. The resulting algorithm produced complex scan patterns which gave strong results for simulated microstructure control.
Many physics problems involve integration in multi-dimensional space whose analytic solution is not available. The integrals can be evaluated using numerical integration methods, but it requires a large computational cost in some cases, so an efficient algorithm plays an important role in solving the physics problems. We propose a novel numerical multi-dimensional integration algorithm using machine learning (ML). After training a ML regression model to mimic a target integrand, the regression model is used to evaluate an approximation of the integral. Then, the difference between the approximation and the true answer is calculated to correct the bias in the approximation of the integral induced by ML prediction errors. Because of the bias correction, the final estimate of the integral is unbiased and has a statistically correct error estimation. Three ML models of multi-layer perceptron, gradient boosting decision tree, and Gaussian process regression algorithms are investigated. The performance of the proposed algorithm is demonstrated on six different families of integrands that typically appear in physics problems at various dimensions and integrand difficulties. The results show that, for the same total number of integrand evaluations, the new algorithm provides integral estimates with more than an order of magnitude smaller uncertainties than those of the VEGAS algorithm in most of the test cases.
Summary Fiber-optic cables cemented outside of the casing of an unconventional well measure crosswell strain changes during fracturing of neighboring wells with low-frequency distributed acoustic sensing (LF-DAS). As a hydraulic fracture intersects an observation well instrumented with fiber-optic cables, the fracture fluid injected at ambient temperatures can cool a section of the sensing fiber. Often, LF-DAS and distributed temperature sensing (DTS) cables are run in tandem, enabling the detection of such cooling events. The increasing use of LF-DAS for characterizing unconventional hydraulic fracture completions demands an investigation of the effects of temperature on the measured strain response by LF-DAS. Researchers have demonstrated that LF-DAS can be used to extract the temporal derivative of temperature for use as a differential-temperature-gradient sensor. However, differential-temperature-gradient sensing is predicated on the ability to filter strain components out of the optical signal. In this work, beginning with an equation for optical phase shift of LF-DAS signals, a model relating strain, temperature, and optical phase shift is explicitly developed. The formula provides insights into the relative strength of strain and temperature effects on the phase shift. The uncertainty in the strain-rate measurements due to thermal effects is estimated. The relationship can also be used to quantify uncertainties in differential-temperature-gradient sensors due to strain perturbations. Additionally, a workflow is presented to simulate the LF-DAS response accounting for both strain and temperature effects. Hydraulic fracture geometries are generated with a 3D fracture simulator for a multistage unconventional completion. The fracture width distributions are imported by a displacement discontinuity method (DDM) program to compute the strain rates along an observation well. An analytic model is used to approximate the temperature in the fracture. Using the derived formulae for optical phase shift, the model outputs are then used to compute the LF-DAS response at a fiber-optic cable, enabling the generation of waterfall plots including both strain and thermal effects. The model results suggest that before, during, and immediately following a fracture intersecting a well instrumented with fiber, the strain on the fiber drives the LF-DAS signal. However, at later times, as completion fluid cools the observation well, the temperature component of the LF-DAS signal can be equal to or exceed the strain component. The modeled results are compared to a published field case in an attempt to enhance the interpretation of LF-DAS waterfall plots. Finally, we propose a sensing configuration to identify the events when “wet fractures” (fractures with fluids) intersect the observation well.
Sensors are increasingly deployed for process monitoring and control. These produce on-line measurements at a high frequency, in parallel with low-frequency laboratory measurements. Compared to laboratory practices, sensor data quality assessment and control practices are far less structured at most utilities. This leads to inaccurate sensor data with unknown uncertainty factors.This chapter shows how to establish standard operating procedures (SOPs) to support sensor data quality assessment and control and subsequent maintenance actions by producing relevant sensor metadata. Furthermore, SOPs are provided for the most commonly used wastewater quality sensors, inspired by utility and academic best practices. This chapter builds on definitions provided in Chapter 3 and provides additional definitions specifically related to sensors maintenance. Chapter 6 complements the methods in this chapter, which are based on reference measurements, with data-analytical techniques.
The paper discusses the HOM interaction effect in a 5-cell 650 MHz HB cavity prototype for the PIP-II project. This effect may lead to flat HOM passbands. In the presence of misalignments, this passband disintegrates into individual cell resonances that are not coupled to the beam pipes and, therefore, have a high-quality factor and R/Q ratio. These resonances can lead to excessive cryogenic loading, increased emittance, and beam instabilities. An analytical theory has been developed to explain this effect, as well as special methods for eliminating it.
This presentation provides an update on the OPTIMA GOAAT project, with emphasis on the cloud-native data platform developed in-house to ingest, manage, and operationalize high-resolution power system data. Since our last NASPI presentation, accessible via OSTI ID #2671437, the project team advanced the design and deployment of a scalable architecture capable of handling both synchronized and non-synchronized streams, including PMU, point-on-wave (POW), COMTRADE, and SCADA data. These materials review the project status, recent progress, and key lessons learned. The core of the presentation examines the architecture and engineering of our cloud-native ingestion and data management platform. We then explain how pipelines were designed to collect, normalize, time-align, store, and serve heterogeneous data at scale. We will discuss design choices such as data models, streaming versus batch ingestion, storage tiers, and interoperability with analytics applications. Practical experiences with cloud-native technologies were shared during the event, including benefits, limitations, and integration challenges in a utility environment, along with methods used to improve performance, reduce latency, and optimize resource usage. The presentation also showcases user interface designs and visualization tools that convert raw measurements and analytics results into intuitive, actionable insights for operators and engineers. During the presentation examples were provided demonstrating how visualization, event views, and summarized analytics enhance situational awareness and support operational decision-making. These use cases illustrate how a well-designed data infrastructure can bridge the gap between high-volume measurements and practical grid operations.
The lipidome comprises the total content of molecular species of each lipid class, and is measured using the analytical techniques of lipidomics. Many liquid chromatography-mass spectrometry (LC-MS) methods have previously been described to characterize the lipidome. However, many lipidomic approaches may not fully uncover the subtleties of lipid molecular species, such as the full fatty acid (FA) composition of certain lipid classes. Here, we describe a stepwise targeted lipidomics approach to characterize the polar and non-polar lipid classes using complementary LC-MS methods. Our “polar” method measures 260 molecular species across 12 polar lipid classes, and is performed using hydrophilic interaction chromatography (HILIC) on a NH2 column to separate lipid classes by their headgroup. Our “non-polar” method measures 254 molecular species across three non-polar lipid classes, separating molecular species on their FA characteristics by reverse phase (RP) chromatography on a C30 column. Five different extraction methods were compared, with an MTBE-based extraction chosen for the final lipidomics workflow. A state-of-the-art strategy to determine and relatively quantify the FA composition of triacylglycerols is also described. This lipidomics workflow was applied to developing, mature, and germinated pennycress seeds/seedlings and found unexpected changes among several lipid molecular species. During development, diacylglycerols predominantly contained long chain length FAs, which contrasted with the very long chain FAs of triacylglycerols in mature seeds. Potential metabolic explanations are discussed. The lack of very long chain fatty acids in diacylglycerols of germinating seeds may indicate very long chain FAs, such as erucic acid, are preferentially channeled into beta-oxidation for energy production.
Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.
Maintaining the superior performance of new analytical instrumentation and pushing the boundaries of ever-smaller particle analysis by methods such as secondary ion mass spectrometry (SIMS) has become limited by the availability of certified particle standards for calibration, quality control, and validation. To meet this growing demand for a reliable and universal approach to the generation of uranium particle reference material, the Pacific Northwest National Laboratory (PNNL) and Savannah River National Laboratory (SRNL) are collaborating on a joint venture that establishes a number of chemical pathways to the fabrication, purification, and stabilization of uranium particle material to within a fixed particle size range, and well-characterized isotopics. PNNL particle standards are designed and tailored to meet the criteria essential to the calibration and benchmarking of instruments used in both the non-destructive and destructive assay of particulate material, typically collected from environmental swipe sampling. Previous testing and optimization of this colloidal approach to particulate material has centred on tight size distributions, singular composition and density, uniform morphology, and tailored isotopic abundances. Extensive validation of this material has since been performed, both at the National Institute of Standards and Technology (NIST) and the International Atomic Energy Agency (IAEA), demonstrating its suitability for quality control needs of large geometry (LG)-SIMS analysis. Building on this momentum, the production of low-enriched uranium particle reference material, with specially tailored isotopics, has moved from development to full operation. Herein, various aspects of the production cycle will be discussed, including detailed accounts of the technical methodologies being employed, as well as insights into the sample characterization and acceptance testing requirements.
We report on several efforts to search for physics beyond the standard model of particle physics at broad energy scales. The Mu2e experiment at Fermilab will search for charged lepton flavor violation via the muon to electron conversion process, which is suppressed in the Standard Model. Mu2e will be operated at a low energy, yet can probe New Physics at very high mass scales (O(1e3 - 1e4 ) TeV). At high energies, the CMS experiment at the CERN LHC continues to deliver an impressive suite of Standard Model measurements and limits on a variety of New Physics signatures. Mu2e is under construction and slated to collect its first physics data in the coming years. This thesis describes work done during the construction phase of Mu2e and focuses on two critical areas: magnetic field modeling and statistical analysis. We describe a novel method for field modeling which we validate using a simulated dataset representing the expected magnetic field in the Detector Solenoid. This method blends a standard least-squares fitting technique that utilizes physically motivated analytical model functions with a novel physics informed network that is constructed to obey Maxwell’s equations. We show the technique can model the field with an accuracy of 10−7 despite the presence of injected noise in the pseudo-measurements at the 10−5 level. We then present preliminary results of the calibration of 3D Hall probes at the sub-10−4 level. These probes will be used to directly measure the Mu2e Detector Solenoid magnetic field on a sparse grid; these measurements serve as the input to the field model fitting. Finally, we describe the first implementation of both an unbinned shape analysis and a Bayesian interpretation applied to Mu2e pseudo-data. Up to 20% tighter limits can be set by the shape analysis compared to a standard cut & count analysis. The AlCap experiment collected data at PSI in 2015 to measure several important quantities related to nuclear muon capture on an aluminum target, which is a significant background process for Mu2e. The neutron emission from muon capture can introduce background hits in the Mu2e detectors and can increase radiation damage in various elements of the apparatus. We present measurements of the neutron group fluence and mean neutron multiplicity for muon capture on aluminum nuclei. Finally, we discuss an analysis of triboson production at CMS using an Effective Field Theory framework. Standard Model triboson production, which was first observed at CMS in 2020, has a relatively small cross section and provides direct access to both anomalous triple gauge couplings and quartic gauge couplings. These couplings, interpreted in the Standard Model Effective Field Theory, are studied in the present work. We target the boosted regime where the background rate is low and yields are enhanced when dimension-6 and dimension-8 Wilson coefficients are non-zero. We do not observe an excess in the data and therefore set bounds on the Wilson coefficients. For dimension-6 coefficients the tightest observed (expected) bounds are set on cW /Λ2 where Λ is the mass scale of new physics; the bounds are [−0.13, 0.12] TeV−2 ([−0.12, 0.12] TeV−2 ) at 95% CL. The tightest bounds in dimension-8 are set on fT,0 / Λ4 ; the observed (expected) bounds at 95% CL are [−0.63, 0.69] TeV−4 ([−0.54, 0.62] TeV−4 ). Additional results are presented which include scenarios where multiple Wilson coefficients are non-zero, the application of signal model clipping to address unitarity violation in Effective Field Theories, and a novel template fit developed for easier reinterpretation of our results.