Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Real-time reconstruction and control of pedestal-top electron density using RMP and gas puff at KSTAR

We report the experimental results of controlling the pedestal-top electron density by applying resonant magnetic perturbation (RMP) with in-vessel control coils and the main gas puff in the 2024-2025 KSTAR experimental campaign. The density is reconstructed using a parameterized $ψ_N$ grid and five channels of line-averaged density measured by the two-colored interferometer (TCI). The reconstruction procedure is accelerated by deploying a multi-layer perceptron to run in approximately 120 µ s and is sufficiently fast for real-time control. A proportional-integral controller was adopted, with the controller gains estimated from the system identification procedure. The experimental results demonstrate that the developed controller can follow a dynamic target while exclusively using both actuators. The absolute percentage errors between the electron density at $ψ_N$ = 0.89 and the target were approximately 1.5% median and a 2.5% average, respectively. The developed controller can even lower the density by using the pump-out mechanism under RMP, and it can follow a more dynamic range of density targets than a single actuator controller. The developed controller will enable experimental scenario exploration within a shot by dynamically setting the density target or maintaining a constant electron density within a discharge.

EFIT↗
DOE OSTI2021

Design and calibration of an optically segmented single volume scatter camera for neutron imaging

The Optically Segmented Single Volume Scatter Camera (OS-SVSC) aims to image neutron sources for non-proliferation applications using the kinematic reconstruction of elastic double-scatter events. Our prototype system consists of 64 EJ-204 organic plastic scintillator bars, each measuring 5mm × 5mm × 200mm and individually wrapped in Teflon tape. The scintillator array is optically coupled to two silicon photomultiplier ArrayJ-60035 64P-PCB arrays, each comprised of 64 individual 6x6mm J-Series sensors arranged in an 8 × 8 array. We report on the design details, including component selections, mechanical design and assembly, and the electronics system. The described design leveraged existing off-the-shelf solutions to support the rapid development of a phase 1 prototype. Several valuable lessons were learned from component and system testing, including those related to the detector's mechanical structure and electrical crosstalk that we conclude originates in the commercial photodetector arrays and the associated custom breakout cards. We detail our calibration efforts, beginning with calibrations for the electronics, based on the IRS3D application-specific integrated circuits, and their associated timing resolutions, ranging from 3090ps. With electronics calibrations applied, energy and position calibrations were performed for a set of edge bars using 22 Na and 90 Sr, respectively, reporting an average resolution of 12.07±0.03mm for energy depositions between 900keVee and 1000keVee. We further demonstrate a position calibration method for the internal bars of the matrix using cosmic-ray muons as an alternative to emission sources that cannot easily access these bars, with an average measured resolution of 14.86+-0.29mm for depositions between 900keVee and 1000keVee. The coincident time resolution reported between pairs of bars measured up to 400ps from muon acquisitions. Energy and position calibration values measured with muons are consistent with those obtained using particle emission sources.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Measurement of nuclear effects in neutrino-argon interactions using generalized kinematic imbalance variables with the MicroBooNE detector

We present a set of new generalized kinematic imbalance variables that can be measured in neutrino scattering. These variables extend previous measurements of kinematic imbalance on the transverse plane and are more sensitive to modeling of nuclear effects. We demonstrate the enhanced power of these variables using simulation and then use the MicroBooNE detector to measure them for the first time. We report flux-integrated single- and double-differential measurements of charged-current muon neutrino scattering on argon using a topology with one muon and one proton in the final state as a function of these novel kinematic imbalance variables. These measurements allow us to demonstrate that the treatment of charged current quasielastic interactions in genie version 2 is inadequate to describe data. As a result, they reveal tensions with more modern generator predictions particularly in regions of phase space where final state interactions are important.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Giant Nonreciprocity and Gyration through Modulation-Induced Hatano-Nelson Coupling in Integrated Photonics

Asymmetric energy exchange interactions, also known as Hatano-Nelson type couplings, enable the study of non-Hermitian physics and associated phenomena like the non-Hermitian skin effect and exceptional points (EPs). Since these interactions are by definition nonreciprocal, there have been very few options for implementations in integrated photonics. Here, in this work, we show that asymmetric couplings are readily achievable in integrated photonic systems through time-domain dynamic modulation. We experimentally study this concept using a two-resonator photonic molecule produced in a lithium niobate on insulator platform that is electro-optically modulated by rf stimuli. We demonstrate the dynamic tuning of the Hatano-Nelson coupling between the resonators, surpassing the asymmetry that has been achieved in previous work, to reach an EP for the first time. We are additionally able to flip the relative sign of the couplings for opposite directions by going past the EP. Using this capability, we show that the through-chain transport can be configured to exhibit both giant (~60 dB) optical contrast as well as photonic gyration or nonreciprocal π phase contrast.

42 ENGINEERING↗
DOE OSTI2024

Intelligent Hierarchical Resilient Operation of Distribution Systems: Implementation and Validation in a Power Hardware-in-the-Loop Simulation Testbed

This paper reports on the structure of a power hardware-in-the-loop (PHIL) simulation testbed that implements, tests, and validates a novel AI-based hierarchical resilient operation model for distribution systems. The testbed implements the central and distributed controllers of the hierarchical resilient operation model and integrates a Digital Real-Time Simulator (DRTS), protective relays, a Real-Time Automation Controller (RTAC), a Software Defined Network (SDN) switch, and a battery energy storage (BES) system. The testbed provides comprehensive real-time visualization and monitoring capability as an advanced situational awareness and operator interface solution. The IEEE 33-node system is used as a test case to test and validate the operation of the model in normal operation and recovery operation after major outages in a fully automated fashion.

Ganjkhani, Mehdi↗
DOE OSTI2023

Automated Controller Hardware-In-The-Loop Testbed for EV Charger Resilience Analysis

This paper focuses on the development of a tool that includes an automated testbed with controls, protection, and communications integrated into a real-time system to provide a platform to generate data sets for failure modes and effects analysis. This tool establishes a value for automation of data generation for different scenarios and addresses the gap of nonexistent field data for different applications and use cases. The features of this tool can further be expanded to include multiple power electronics models, communication protocols, and scaled system architectures. This general framework was evaluated for a DC fast charger system use case to provide quantitative solution for resiliency.

Starke, Michael↗
DOE OSTI2024

Anomaly Detection and Mitigation for Dynamic Frequency Regulation in Hydropower-Battery Systems

Hydropower operators and energy storage providers are increasingly interested in participating in frequency regulation services, driven by the incentives offered by independent system operators, such as the PJM Interconnection. This transition, however, unfolds against the backdrop of a modernizing and rapidly digitizing power grid, exposing the integrated legacy infrastructure to a multitude of cybersecurity threats. This work presents an approach for developing an anomaly detection and mitigation system to address cybersecurity challenges during the participation of a hydropower-integrated battery energy storage system (BESS) in a frequency regulation market. The applied anomaly detector utilizes machine learning algorithms to provide detailed classification of cyber-physical events. Later, the applied mitigation system triggers predefined corrective actions to minimize the impact of data integrity attacks on the regulation market and system stability. We evaluated the proposed approach on a hydropower-integrated BESS topology, specifically analyzing the slow regulation signal (Reg A) coming from the PJM market. Our simulation results demonstrate that the proposed approach performs well in detecting data integrity attacks within the allocated time frame and also minimizes the system's transient instability during the participation of hydropower and BESS in the regulation market.

battery energy storage system↗
DOE OSTI2024

Modeling and ZVS Operation of the Isolated Modular Multilevel DC–DC Converter With a Unified Trapezoidal Wave Modulation

Here, this article introduces a unified trapezoidal wave (UTW) modulation scheme for the isolated modular multilevel dc–dc (IM2dc) converter, which consolidates multiple existing modulation strategies for the IM2dc converter, including the quasi-square wave (QSW) modulation, trapezoidal wave modulation, and sinusoidal wave modulation, into a unified framework. Furthermore, this article introduces the harmonic state-space (HSS) equations to model the IM2dc converter based on the UTW modulation method. The HSS model operates in the frequency domain, enabling it to circumvent the complexities associated with time-domain analysis and seamlessly integrate with the UTW modulation. This article proceeds to analyze the real and reactive power transfer characteristics of the IM2dc converter, as well as the power factor, considering the influence of multiple modulation parameters. Subsequently, it delves into the examination of zero voltage switching (ZVS) conditions for the IM2dc converter based on the UTW modulation and the HSS model. The complete ZVS boundaries of the IM2dc converter, taking various voltage ratios into account, are derived. This article also illustrates the effects of including harmonic orders in the modeling process, modulation parameters, and internal harmonic ripples on the ZVS boundaries. Finally, experimental validation of the analyses is conducted on a down-scaled prototype.

42 ENGINEERING↗
DOE OSTI2021

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA’s diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Fully Homomorphic Encryption

This code implements a Fully Homomorphic Encryption (FHE) system, enabling secure computation on encrypted data without requiring decryption. It supports encryption, decryption, and homomorphic operations like matrix multiplication and addition. This code is adaptable for integrating FHE into linear-time invariant (LTI) systems, including digital control and filtering. With proper configuration from subject matter expertise, encrypted system parameters and signals can be manipulated to perform tasks like state updates, output calculations, and convolution in the encrypted domain. By preserving the structure of LTI systems while ensuring privacy, the framework facilitates secure applications in areas such as autonomous systems, signal processing, and industrial automation. The code initializes the encryption system using parameters provided in the env dictionary. These parameters include the ciphertext modulus, key dimension, plaintext fixed-point scaling factor, and noise bound. During initialization, a secret key is generated, which is essential for encrypting and decrypting data securely. The modular design allows users to tailor these parameters to specific use cases or security requirements. The code implements multiple cryptographic schemes. The learning with errors (LWE) encryption method encodes cleartext message to their plaintext fixed-point representation then encrypted into ciphertext space with additive noise. This noise ensures the security of the scheme, relying on the computational hardness of the LWE problem. The code also includes the Gentry-Sahai-Waters (GSW) scheme based off the LWE problem. Homomorphic matrix multiplication is performed between the LWE and GSW to encrypted data. This is achieved using a decomposition function on the LWE ciphertext during the multiplication operation. For higher-dimensional data, the code includes a method to encrypt entire matrices (GSWMat) using GSW encryption. These encrypted matrices can then be used for homomorphic matrix multiplications (MatMult). The decryption function uses the secret key to recover the original plaintext, removing the added noise and scaling that was originally applied during encryption.

Lois, Roberts [Idaho National Laboratory (INL), Id↗
DOE OSTI2026

Packets Design System (Packets) v2.0

A Design System is a comprehensive framework of reusable components, standards, and guidelines. Packets is the design system created by developers at ESnet to give our software suite a common look and feel, as well as speed up production time. We would like to integrate it across all our software platforms, including some that are already or are in the process of being open sourced. In order to do this, we would also need to open source a version of Packets.

Kafader, James [Lawrence Berkeley National Laborat↗
DOE OSTI2020

Accelerate M-TIP on GPUs and deploy to Summit and NERSC-9 (against simulated data) WBS 2.2.4.05 ExaFEL, Milestone ADSE13-199

We evaluated the orientation matching step in the M-TIP SPI workflow for potential offloading to accelerators. We ported the code to GPUs, benchmarked it, optimized and down-selected the best versions. The accelerated version of the orientation matching code that was developed at LANL (LANL GPU v3) is 34-55X faster than sequential, 2.4-4.9X faster than the fastest OpenMP open source version we found (FAISS OpenMP) and 1.5-4X faster than the fastest GPU open source version we found (FAISS GPU). Summit single-node GPU versions were somewhat faster than Cori GPU. Image size plays a role; mid-range image sizes take more time. The LANL CUDA multi-node, multi-GPU implementation shows mostly linear strong scaling. I/O also plays a large role; splitting data into parts improves read time and burst buffers dramatically improve read times. This work will be integrated into the M-TIP workflow as part of the next milestone ADSE13-193.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Demonstration of high-rate manufacturing of π-joint using additive manufacturing

Bonded composite π-joints are critical to support wing geometries in aerospace/aircraft structures. Its lightweight and robust load-bearing capabilities/design guarantee high structural integrity and reliability. Current manufacturing techniques of π-joints are complex and laborious, and do not always offer high-quality joining. This study presents a high-rate manufacturing approach to fabricate π-joints using an integrated Additive Manufacturing and Compression Molding (AMCM) process. The joints were manufactured using unidirectional aerospace-grade LMPAEK/CF at 60% fiber loading for the stiffeners/skin and LM-PAEK/GF as the adhesive and characterized to evaluate their microstructure and strength. While the results indicate void formation at the joint caused by inadequate pressure, the peak tensile load of 600N recorded is akin to injection-autoclave molded π-joints, which validates the potential of this high-rate integrated approach to eliminate time-consuming and labour-intensive processes of making π-joints to reduce cost and produce high-quality joints.

Marathe, Umesh [ORNL]↗
DOE OSTI2021

Valuation of Wind Energy Turbines Using Volatility of Wind and Price

The limitedness of the nonrenewable local energy resources in Israel, even in the background of the later gas fields’ findings, continues to force the state to devote various efforts towards ‘green’ energy development. These efforts include installations, both for the solar and for wind energy, thus improving the diversity of energy sources. While the standard discounted cash flow (DCF) method using the net present value (NPV) criterion is extensively adopted to evaluate investments, the standard DCF method is inappropriate for the rapidly changing investment climate and for the managerial flexibility in investment decisions. In recent years, the real options analysis (ROA) technique has been widely applied in many studies for the valuation of renewable energy investment projects. Taking into account the above background, we apply, in this study, the real options analysis approach for the valuation of wind energy turbines and apply it to the analysis of wind energy economic potential in Israel, which is the context of our work. We hypothesize that due to nature of wind energy production uncertainties, the ROA method is better than the alternative. The novelty of this paper includes the following: real world wind statistics of the Merom Golan site in Israel (velocity 3.73 m/s, with a standard deviation of 2.03 m/s), a realistic power generation estimation (power generation of 1205.84 kW with a standard deviation of about 0.5% in annual value which is worth about 1.3 M$ per annum), and an economic model to evaluate the profitability of such a project. We thus discuss the existing challenges of diversifying renewable energy sources in Israel by adding wind installations. Our motivation is to introduce a method which will allow investors and officials to take into account uncertainties when deciding in investing in such wind installations. The outcomes of the paper, which are obtained using the method of Weibull statistics and the Black–Scholes ROA technique, include the result that market price volatility adds to the uncertainties much more than any wind fluctuations, provided that the analysis is integrated over a long enough time.

17 WIND ENERGY↗
DOE OSTI2022

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗
DOE OSTI2023

Prompt Detection of Fast Optical Bursts with the Vera C. Rubin Observatory

Abstract The transient optical sky has remained largely unexplored on very short timescales. While there have been some experiments searching for optical transients from minutes to years, none have had the capability to distinguish millisecond fast optical bursts (FOBs). Such very fast transients could be the optical counterparts of fast radio bursts, the prompt emission from γ -ray bursts, or other previously unknown phenomena. Here, we investigate a novel approach to the serendipitous detection of FOBs, which relies on searching for anomalous spatial images. In particular, due to their short duration, the seeing-distorted images of FOBs should look characteristically different than those of steady sources in a standard optical exposure of finite duration. We apply this idea to simulated observations with the Vera C. Rubin Observatory, produced by tracing individual photons through a turbulent atmosphere, and down through the optics and camera of the Rubin telescope. We compare these simulated images to steady-source star simulations in 15 s integrations, the nominal Rubin exposure time. We report the classification accuracy results of a neural network classifier for distinguishing FOBs from steady sources. From this classifier, we derive constraints in duration–intensity parameter space for unambiguously identifying FOBs in Rubin observations. We conclude with estimates of the total number of detections of FOB counterparts to FRBs expected during the 10 yr Rubin Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Microwave diagnostics and passive sensors for pipeline, well-bore, and boiler-tube monitoring

Materials, methods of making, and methods of using an integrated wireless detector for real time interrogating metallic tubular structures comprising: an RF patch antenna; a passive surface acoustic wave (SAW) sensor; and data analytic methodologies. An embodiment relates to interrogating a metallic structure having a uniform cross section using an antenna which launches electromagnetic radiation. A sensor may be located within the structure is configured to re-emit electromagnetic radiation modified depending on parameters for which the sensor has been functionalized. An antenna may receive radiation as modified by the sensor, or may receive the transmitted or backscattered radiation directly, without use of a sensor. The antenna then communicates wirelessly with an interrogator providing data which may be used to understand the operational status of the structure in real-time.

Ohodnicki, Paul R.↗