Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

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

Comparing complex impedance and bias step measurements of Simons Observatory transition edge sensors

The Simons Observatory (SO) will perform ground-based observations of the cosmic microwave background (CMB) with several small and large aperture telescopes, each outfitted with thousands to tens of thousands of superconducting aluminum manganese (AlMn) transition-edge sensor bolometers (TESs). In-situ characterization of TES responsivities and effective time constants will be required multiple times each observing-day for calibrating time-streams during CMB map-making. Effective time constants are typically estimated in the field by briefly applying small amplitude square-waves on top of the TES DC biases, and fitting exponential decays in the bolometer response. These so-called “bias step” measurements can be rapidly implemented across entire arrays and therefore are attractive because they take up little observing time. However, individual detector complex impedance measurements, while too slow to implement during observations, can provide a fuller picture of the TES model and a better understanding of its temporal response. Here, we present the results of dark TES characterization of many prototype SO bolometers and compare the effective thermal time constants measured via bias steps to those derived from complex impedance data.

Nicholas F. Cothard↗

BETO 2021 Peer Review - WBS 4.1.2.32: Bioeconomy Scenario Analysis

The Bioeconomy Scenario Analysis project uses systems thinking and analysis to assess current and/or prospective techno-economics, research and development, deployment strategies, policy, and market conditions and their impact on the potential development trajectories of the bioenergy industry over time. Results from this project include identification of opportunities and constraints to industrial development, quantification of multiple metrics (energy, economic, environmental) and informing researchers, decision makers, and industry of the steps needed for a sustainable, nationwide biofuels industry. Analyses from this project enable the creation of a bioenergy industry by (1) inciting policy-makers to explore scenarios for nationwide biofuels production, identifying policy actions consistent with pathways for growth; (2) improving industry’s understanding of industry growth potential under different technology and investment conditions, better targeting their development efforts; and (3) providing universities and other interested stakeholders with tools and analyses that can be adapted to meet research and teaching objectives, connecting students with careers that build the industry. One of the many modeling tools used in this project, the Biomass Scenario Model (BSM) is a publicly-available, unique, validated, state-of-the-art, award-winning, fourth-generation model of the domestic biofuels supply chain which explicitly focuses on how and under what conditions biofuel technologies might be deployed to contribute to the U.S. transportation energy sector. We use models like the BSM to examine the implications of policies and incentives as well as their potential side-effects. The BSM uses a system-dynamics simulation to model dynamic interactions and transitions across the supply chain; it tracks the deployment of biofuels given industrial learning and the reaction of the investment community in the context of land availability, projected oil markets, consumer demand for biofuels, and government policies over time. Under expected market conditions, analyses using the BSM suggest that the biofuels industry may require significant external actions in the early years to thrive. Interventions that accelerate the industrial learning process (e.g. operation of pre-commercial and commercial facilities) have been identified as having strong influence in starting the growth of a commercial biofuel industry. Policies which are coordinated across the whole supply chain in BSM foster the growth of the biofuels industry and production of tens of billions of gallons of biofuels may occur under sufficiently favorable conditions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

QuYBE - An Algebraic Compiler for Quantum Circuit Compression

QuYBE is an open-source algebraic compiler for the compression of quantum circuits. It has been applied for the efficient simulation of the Heisenberg Hamiltonian on quantum computers. Currently, it can simulate the time dynamics of one-dimensional chains. It includes modules to generate the quantum circuits for the above as well as produce the compressed circuits, which are independent of the time step. It utilizes the Yang-Baxter equation (YBE) to perform the compression. QuYBE enables users to seamlessly design, execute, and analyze the time dynamics of the Heisenberg Hamiltonian on quantum computers. QuYBE is the first step toward making the YBE technique available to a broader community of scientists from multiple domains. The QuYBE compiler is available at https://github.com/ZichangHe/QuYBE.

Gulania, Sahil↗

GCAM–GLORY v1.0: representing global reservoir water storage in a multi-sector human–Earth system model

Abstract. Reservoirs play a significant role in modifying the spatiotemporal availability of surface water to meet multi-sector human demands, despite representing a relatively small fraction of the global water budget. Yet the integrated modeling frameworks that explore the interactions among climate, land, energy, water, and socioeconomic systems at a global scale often contain limited representations of water storage dynamics that incorporate feedbacks from other systems. In this study, we implement a representation of water storage in the Global Change Analysis Model (GCAM) to enable the exploration of the future role (e.g., expansion) of reservoir water storage globally in meeting demands for, and evolving in response to interactions with, the climate, land, and energy systems. GCAM represents 235 global water basins, operates at 5-year time steps, and uses supply curves to capture economic competition among renewable water (now including reservoirs), non-renewable groundwater, and desalination. Our approach consists of developing the GLObal Reservoir Yield (GLORY) model, which uses a linear programming (LP)-based optimization algorithm and dynamically linking GLORY with GCAM. The new coupled GCAM–GLORY approach improves the representation of reservoir water storage in GCAM in several ways. First, the GLORY model identifies the cost of supplying increasing levels of water supply from reservoir storage by considering regional physical and economic factors, such as evolving monthly reservoir inflows and demands, and the leveled cost of constructing additional reservoir storage capacity. Second, by passing those costs to GCAM, GLORY enables the exploration of future regional reservoir expansion pathways and their response to climate and socioeconomic drivers. To guide the model toward reasonable reservoir expansion pathways, GLORY applies a diverse array of feasibility constraints related to protected land, population, water sources, and cropland. Finally, the GLORY–GCAM feedback loop allows evolving water demands from GCAM to inform GLORY, resulting in an updated supply curve at each time step, thus enabling GCAM to establish a more meaningful economic value of water. This study improves our understanding of the sensitivity of reservoir water supply to multiple physical and economic dimensions, such as sub-annual variations in climate conditions and human water demands, especially for basins experiencing socioeconomic droughts.

54 ENVIRONMENTAL SCIENCES↗

Hamiltonian learning using machine-learning models trained with continuous measurements

Here, we build upon recent work on the use of machine-learning models to estimate Hamiltonian parameters using continuous weak measurement of qubits as input. We consider two settings for the training of our model: (1) supervised learning, where the weak-measurement training record can be labeled with known Hamiltonian parameters, and (2) unsupervised learning, where no labels are available. The first has the advantage of not requiring an explicit representation of the quantum state, thus potentially scaling very favorably to a larger number of qubits. The second requires the implementation of a physical model to map the Hamiltonian parameters to a measurement record, which we implement using an integrator of the physical model with a recurrent neural network to provide a model-free correction at every time step to account for small effects not captured by the physical model. We test our construction on a system of two qubits and demonstrate accurate prediction of multiple physical parameters in both the supervised context and the unsupervised context. We demonstrate that the model benefits from larger training sets, establishing that it is “learning,” and we show robustness regarding errors in the assumed physical model by achieving accurate parameter estimation in the presence of unanticipated single-particle relaxation.

97 MATHEMATICS AND COMPUTING↗

Energy yields in the prebiotic synthesis of hydrogen cyanide and formaldehyde

Prebiotic experiments are usually reported in terms of carbon yields, i.e., the yield of product based on the total carbon in the system. These experiments usually involve a large input of energy and are designed to maximize the yields of product. However, large inputs of energy result in multiple activation of the reactants and products. A more realistic prebiotic experiment is to remove the products of the activation step so they are not exposed a second time to the energy source. This is equivalent to transporting the products synthesized in the primitive atmosphere to the ocean, and thereby protecting them from destruction by atmospheric energy sources. Experiments of this type, using lower inputs of energy, give energy yields (moles of products/joule) which can be used to estimate the relative importance of the different energy sources on the primitive earth. Simulated prebiotic atmospheres containing either CH4, CO or CO2 with N2, H2O and variable amounts of H2 were subjected to a high frequency Tesla coil. Samples of the aqueous phase were taken at various time intervals from 1 hr to 7 days, and the energy yields were obtained by extrapolation to zero time. The samples were analyzed for HCN with the cyanide electrode and for H2CO by chromotropic acid. The spark energy was estimated by calorimetry. The temperature rise in an insulated discharge flask was compared with the temperature rise from a resistance heater in the same flask. These results will be compared with calculated production rates of HCN and H2CO from lightning and a number of photochemical processes on the primitive Earth.

Stribling, R.↗

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗

Poplar lignin structural changes during extraction in γ-valerolactone (GVL)

In this paper, we describe an approach for producing both high quality and high quantity of lignin through studying the structural change of lignin during treatment of poplar wood in γ-valerolactone (GVL) for a range of temperatures (from 80 to 120 °C) and reaction time at temperature (from 1 to 24 h). Throughout the study, various techniques, including nuclear magnetic resonance (NMR) spectroscopies (solution- and gel-state 1 H –13 C 2D HSQC and 31 P) and gel-permeation chromatography (GPC) were applied to characterize the lignin structures. As the GVL-extracted lignin yield increases, the level of β-ether units decreases and the level of condensation products increases. The β-ether content, the aliphatic hydroxyl group content, and the molecular weight of the GVL-extracted lignin fractions were close to the poplar lignin from other preparation methods (e.g., enzyme lignin). A two-step hydrolytic process (120 °C, 2 × 15 min) gave a higher lignin yield (56.5% vs. 54.8%) with three times higher β-ether content (31.9% vs. 10.6%) than lignin extracted from a single-step process at 120 °C for 1 h. The results demonstrate that multiple-step cycling of cosolvent-assisted hydrolysis can help preserve more of the virgin ether-bond structures of GVL-extracted poplar lignin. Such a strategy can also be applied to a fully continuous-flow reactor system in future research to further improve both the productivity and quality of GVL-extracted lignin.

2D HSQC NMR↗

Proto-Examples of Data Access and Visualization Components of a Potential Cloud-Based GEOSS-AI System

Once a research or application problem has been identified, one logical next step is to search for available relevant data products. Thus, an early component of a potential GEOSS-AI system, in the continuum between observations and end point research, applications, and decision making, would be one that enables transparent data discovery and access by users. Such a component might be effected via the systems data agents. Presumably, some kind of data cataloging has already been implemented, e.g., in the GEOSS Common Infrastructure (GCI). Both the agents and cataloging could also leverage existing resources external to the system. The system would have some means to accept and integrate user-contributed agents. The need or desirability for some data format internal to the system should be evaluated. Another early component would be one that facilitates browsing visualization of the data, as well as some basic analyses.Three ongoing projects at the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) provide possible proto-examples of potential data access and visualization components of a cloud-based GEOSS-AI system. 1. Reorganizing data archived as time-step arrays to point-time series (data rods), as well as leveraging the NASA Simple Subset Wizard (SSW), to significantly increase the number of data products available, at multiple NASA data centers, for production as on-the-fly (virtual) data rods. SSWs data discovery is based on OpenSearch. Both pre-generated and virtual data rods are accessible via Web services. 2. Developing Web Feature Services to publish the metadata, and expose the locations, of pre-generated and virtual data rods in the GEOSS Portal and enable direct access of the data via Web services. SSW is also leveraged to increase the availability of both NASA and non-NASA data.3.Federating NASA Giovanni (Geospatial Interactive Online Visualization and Analysis Interface), for multi-sensor data exploration, that would allow each cooperating data center, currently the NASA Distributed Active Archive Centers (DAACs), to configure its own Giovanni deployment, while also allowing all the deployments to incorporate each others data. A federated Giovanni comprises Giovanni Virtual Machines, which can be run on local servers or in the cloud.

access↗

Small Body Hopper Mobility Concepts

A propellant-saving hopper mobility system was studied that could help facilitate the exploration of small bodies such as Phobos for long-duration human missions. The NASA Evolvable Mars Campaign (EMC) has proposed a mission to the moons of Mars as a transitional step for eventual Mars surface exploration. While a Mars transit habitat would be parked in High-Mars Orbit (HMO), crew members would visit the surface of Phobos multiple times for up to 14 days duration (up to 50 days at a time with logistics support). This paper describes a small body surface mobility concept that is capable of transporting a small, two-person Pressurized Exploration Vehicle (PEV) cabin to various sites of interest in the low-gravity environment. Using stored kinetic energy between bounces, a propellant-saving hopper mobility system can release the energy to vector the vehicle away from the surface in a specified direction. Alternatively, the stored energy can be retained for later use while the vehicle is stationary in respect to the surface. The hopper actuation was modeled using a variety of launch velocities, and the hopper mobility was evaluated using NASA Exploration Systems Simulations (NExSyS) for transit between surface sites of interest. A hopper system with linear electromagnetic motors and mechanical spring actuators coupled with Control Moment Gyroscope (CMG) for attitude control will use renewable electrical power, resulting in a significant propellant savings.

Phobos↗

Refinement and application of acoustic impulse technique to study nozzle transmission characteristics

An improved acoustic impulse technique was developed and was used to study the transmission characteristics of duct/nozzle systems. To accomplish the above objective, various problems associated with the existing spark-discharge impulse technique were first studied. These included (1) the nonlinear behavior of high intensity pulses, (2) the contamination of the signal with flow noise, (3) low signal-to-noise ratio at high exhaust velocities, and (4) the inability to control or shape the signal generated by the source, specially when multiple spark points were used as the source. The first step to resolve these problems was the replacement of the spark-discharge source with electroacoustic driver(s). These included (1) synthesizing on acoustic impulse with acoustic driver(s) to control and shape the output signal, (2) time domain signal averaging to remove flow noise from the contaminated signal, (3) signal editing to remove unwanted portions of the time history, (4) spectral averaging, and (5) numerical smoothing. The acoustic power measurement technique was improved by taking multiple induct measurements and by a modal decomposition process to account for the contribution of higher order modes in the power computation. The improved acoustic impulse technique was then validated by comparing the results derived by an impedance tube method. The mechanism of acoustic power loss, that occurs when sound is transmitted through nozzle terminations, was investigated. Finally, the refined impulse technique was applied to obtain more accurate results for the acoustic transmission characteristics of a conical nozzle and a multi-lobe multi-tube supressor nozzle.

Salikuddin, M.↗

Feature-Based PMU Event Classification under Variable PMU Participation and Overlapping Events

Danovo Energy Solution's presented its paper named: Feature-Based PMU Event Classification under Variable PMU Participation and Overlapping Events at the 2026 Georgia Tech Fault & Disturbance Analysis Conference. The full paper can be found at OSTI ID# 3169150 Paper Abstract—Phasor Measurement Units (PMUs) stream time synchronized, high-resolution measurements from the grid, enabling data-driven techniques for event detection and classification. Accurate event classification improves grid reliability and stability. Events can be detected by varying numbers of PMUs and exhibit different durations depending on the event type. This variability challenges standard classifiers that require uniform input sizes. Moreover, multiple events may coincide, which increases classification complexity. Standard classifiers assign each instance to the class with the highest predicted probability, whereas overlapping events may exhibit comparable probabilities across multiple classes. In this study, to handle data size variability, we extract a wide range of time–frequency domain features from all available PMUs for each event into a fixed-length vector, facilitating the application of standard machine learning classifiers, including Random Forest, XGBoost, LightGBM, Support Vector Machine, and Multilayer Perceptron. To account for overlapping events, a probabilistic post-processing step is applied. For a given data instance, if multiple predicted class probabilities exceed 30% and the differences between them are less than 10%, the event is assigned to multiple classes. Experiments using real-world PMU data demonstrate that the Random Forest and XGBoost models achieve the highest accuracy, while the proposed post-processing method yields perfect classification performance on external unseen test sets.

Nematirad, Reza [Danova Energy Solutions]↗

Feature-Based PMU Event Classification under Variable PMU Participation and Overlapping Events

This paper is the basis for a presentation help at the 2026 Georgia Tech Fault & Disturbance Analysis Conference, which can be found at OSTI # 3168287 Paper Abstract—Phasor Measurement Units (PMUs) stream time synchronized, high-resolution measurements from the grid, enabling data-driven techniques for event detection and classification. Accurate event classification improves grid reliability and stability. Events can be detected by varying numbers of PMUs and exhibit different durations depending on the event type. This variability challenges standard classifiers that require uniform input sizes. Moreover, multiple events may coincide, which increases classification complexity. Standard classifiers assign each instance to the class with the highest predicted probability, whereas overlapping events may exhibit comparable probabilities across multiple classes. In this study, to handle data size variability, we extract a wide range of time–frequency domain features from all available PMUs for each event into a fixed-length vector, facilitating the application of standard machine learning classifiers, including Random Forest, XGBoost, LightGBM, Support Vector Machine, and Multilayer Perceptron. To account for overlapping events, a probabilistic post-processing step is applied. For a given data instance, if multiple predicted class probabilities exceed 30% and the differences between them are less than 10%, the event is assigned to multiple classes. Experiments using real-world PMU data demonstrate that the Random Forest and XGBoost models achieve the highest accuracy, while the proposed post-processing method yields perfect classification performance on external unseen test sets.

Nematirad, Reza [Danovo Energy Solutions]↗

Full event interpretation with machine-learning-based particle-flow reconstruction in the CMS detector

The particle-flow (PF) algorithm constructs a global description of each particle collision by producing a comprehensive list of final-state particles, and is central to event reconstruction in the CMS experiment at the CERN LHC. The existing PF implementation relies on physics-motivated heuristics and assumptions that can be replaced by machine-learning (ML) models trained directly on simulated data and naturally suited to modern graphics processing units (GPUs). A state-of-the-art ML-based PF (MLPF) reconstruction algorithm, implemented within the CMS software framework, is presented. The MLPF algorithm performs a learnable full-event reconstruction on GPUs, generalizes across detector conditions and collision energies, and replaces multiple modular reconstruction steps with a single unified model. Physics performance comparable to standard PF reconstruction is achieved in both simulation and data, with improved jet energy resolution and inference time. In simulated top quark-antiquark events under LHC Run-3 (2023-2024) conditions, the jet energy resolution improves by 10-20% for jets with transverse momentum between 30-100 GeV. Inference time is evaluated using simulated multijet events, with a median of $20\,\hbox {ms}$ per event on an Nvidia L4 GPU, compared to approximately $110\,\hbox {ms}$ for the standard CMS PF reconstruction.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

Integrated circuit reliability testing

A technique is described for use in determining the reliability of microscopic conductors deposited on an uneven surface of an integrated circuit device. A wafer containing integrated circuit chips is formed with a test area having regions of different heights. At the time the conductors are formed on the chip areas of the wafer, an elongated serpentine assay conductor is deposited on the test area so the assay conductor extends over multiple steps between regions of different heights. Also, a first test conductor is deposited in the test area upon a uniform region of first height, and a second test conductor is deposited in the test area upon a uniform region of second height. The occurrence of high resistances at the steps between regions of different height is indicated by deriving the measured length of the serpentine conductor using the resistance measured between the ends of the serpentine conductor, and comparing that to the design length of the serpentine conductor. The percentage by which the measured length exceeds the design length, at which the integrated circuit will be discarded, depends on the required reliability of the integrated circuit.

Buehler, Martin G.↗

Direct Adaptive Aircraft Control Using Dynamic Cell Structure Neural Networks

A Dynamic Cell Structure (DCS) Neural Network was developed which learns topology representing networks (TRNS) of F-15 aircraft aerodynamic stability and control derivatives. The network is integrated into a direct adaptive tracking controller. The combination produces a robust adaptive architecture capable of handling multiple accident and off- nominal flight scenarios. This paper describes the DCS network and modifications to the parameter estimation procedure. The work represents one step towards an integrated real-time reconfiguration control architecture for rapid prototyping of new aircraft designs. Performance was evaluated using three off-line benchmarks and on-line nonlinear Virtual Reality simulation. Flight control was evaluated under scenarios including differential stabilator lock, soft sensor failure, control and stability derivative variations, and air turbulence.

Jorgensen, Charles C.↗

Integrated circuit reliability testing

A technique is described for use in determining the reliability of microscopic conductors deposited on an uneven surface of an integrated circuit device. A wafer containing integrated circuit chips is formed with a test area having regions of different heights. At the time the conductors are formed on the chip areas of the wafer, an elongated serpentine assay conductor is deposited on the test area so the assay conductor extends over multiple steps between regions of different heights. Also, a first test conductor is deposited in the test area upon a uniform region of first height, and a second test conductor is deposited in the test area upon a uniform region of second height. The occurrence of high resistances at the steps between regions of different height is indicated by deriving the measured length of the serpentine conductor using the resistance measured between the ends of the serpentine conductor, and comparing that to the design length of the serpentine conductor. The percentage by which the measured length exceeds the design length, at which the integrated circuit will be discarded, depends on the required reliability of the integrated circuit.

Buehler, Martin G.↗

Microgravity Foam Structure and Rheology

To exploit rheological and multiple-light scattering techniques, and ultimately microgravity conditions, in order to quantify and elucidate the unusual elastic character of foams in terms of their underlying microscopic structure and dynamics. Special interest is in determining how this elastic character vanishes, i.e. how the foam melts into a simple viscous liquid, as a function of both increasing liquid content and shear strain rate. The unusual elastic character of foams will be quantified macroscopically by measurement of the shear stress as a function of static shear strain, shear strain rate, and time following a step strain; such data will be analyzed in terms of a yield stress, a static shear modulus, and dynamical time scales. Microscopic information about bubble packing and rearrangement dynamics, from which these macroscopic non-Newtonian properties presumably arise, will be obtained non-invasively by novel multiple-light scattering diagnostics such as Diffusing-Wave Spectroscopy (DWS). Quantitative trends with materials parameters, such as average bubble size, and liquid content, will be sought in order to elucidate the fundamental connection between the microscopic structure and dynamics and the macroscopic rheology.

Durian, Douglas J.↗