Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Emulator”

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 1,171 records · Page 65

Hofmann Stability Charts Revisited for PIP-II: From Classical Theory to Assumption-Free and ML-Driven Maps

The Hofmann stability chart remains a standard for visualizing parametric resonances in space-charge–dominated linacs, but its use typically relies on non-oscillatory Vlasov dispersion relations with simplifying assumptions (continuous focusing, KV phase space, linear optics, limited transverse–longitudinal coupling). We revisit the chart for the PIP-II linac along three tracks. (1) We reproduce the conventional maps in the (νz/νx, νx/ν0x) plane for relevant εz/εx, providing a validated reference. (2) We remove key assumptions by deriving stability surfaces directly from multi-particle tracking with realistic lattice discreteness, RF defocusing, solenoid/quad optics, and bunched-beam dynamics; local tunes and early-time growth rates are estimated from envelope oscillations and projected to the same coordinates. These assumption-reduced maps recover the canonical stopbands while revealing shifts and broadenings driven by tune modulation, non-KV distributions, and transverse–longitudinal coupling at PIP-II intensities. (3) We train a compact machine-learning surrogate that emulates the growth surface from zero-current optics, tune depression, emittance ratio, bunching factor, and selected lattice descriptors, enabling rapid scans and online working-point selection. We compare the three representations on representative PIP-II sections and discuss implications for commissioning guard bands, resonance avoidance, and routine operations.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

BIOSHOT

BioSHOT is an internally developed SCEPTRE emulation for bio lab environments, meant to train operators to better defend their environments. This project aims to take the existing base topology, add additional realisitic bio software (e.g. GNU Health) and add the final documentation and packaging for release to open source.

Jobe, Tyler Lee [Sandia National Laboratories (SNL↗

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence N - Raw Data

Sequence N: Sin AOA, Rotating (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the presence of rotational influences by varying blade pitch angle. Test sequence N used an upwind, rigid turbine with a 0° cone angle. The wind speed was 15 m/s, and the yaw angle was 0° throughout the sequence. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the Ohio State University (OSU) wind tunnel (Reuss Ramsey et al. 1995) and the 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation N, followed by two digits for the span location (30, 47, 63, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, most of the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence O - Raw Data

Sequence O: Sin AOA, Parked (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the absence of rotational influences by varying blade pitch angle. This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speed ranged from 18.9 m/s to 39.3 m/s, and the yaw angle was 0° throughout the sequence. The rotor was parked with the instrumented blade fixed at 0° azimuth, and the rotor lock was installed (see Appendix A). Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to -99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the OSU wind tunnel (Reuss Ramsey et al. 1995) as well as 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation O, followed by two digits for the span location (30, 47, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗

OC6 Phase III: Validation of Wind Turbine Aerodynamic Loading During Surge/Pitch Motion

The objective of Phase III of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project was to validate the accuracy of aerodynamic load predictions by offshore wind modeling tools for a floating offshore wind turbine as it experiences large surge-translational and pitch-rotational motion, as would occur during normal operation. The test data considered were generated at Politecnico di Milano – a wind tunnel with a robotic excitation system to emulate wave loading on the wind turbine. Testing was performed using a scaled version of the DTU 10-MW reference model, and motion was prescribed as harmonic oscillations in the surge and pitch directions, independently. A variety of models were examined in the project, performing a three-way validation between engineering-level tools, higher-fidelity tools, and measurement data from two wind tunnel experimental campaigns. The Load Cases (LC) considered in this testing are as follows: LC 1.X - Steady Wind LC 2.X - Unsteady Wind - Surge Motion LC 3.X - Unsteady Wind - Pitch Motion Details on the results of the OC6 Phase III campaign can be found in the following reference: Bergua R., et al. OC6 project Phase III: validation of the aerodynamic loading on a wind turbine rotor undergoing large motion caused by a floating support structure, Wind Energ. Sci., 8, 465–485, https://doi.org/10.5194/wes-8-465-2023, 2023.

17 WIND ENERGY↗

Latency considerations for stochastic optimizers in variational quantum algorithms

Variational quantum algorithms, which have risen to prominence in the noisy intermediate-scale quantum setting, require the implementation of a stochastic optimizer on classical hardware. To date, most research has employed algorithms based on the stochastic gradient iteration as the stochastic classical optimizer. In this work we propose instead using stochastic optimization algorithms that yield stochastic processes emulating the dynamics of classical deterministic algorithms. This approach results in methods with theoretically superior worst-case iteration complexities, at the expense of greater per-iteration sample (shot) complexities. We investigate this trade-off both theoretically and empirically and conclude that preferences for a choice of stochastic optimizer should explicitly depend on a function of both latency and shot execution times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ground state energy and magnetization curve of a frustrated magnetic system from real-time evolution on a digital quantum processor

Models of interacting many-body quantum systems that may realize new exotic phases of matter, notably quantum spin liquids, are challenging to study using even state-of-the-art classical methods such as tensor network simulations. Quantum computing provides a promising route for overcoming these difficulties to find ground states, dynamics, and more. In this paper, we argue that recently developed hybrid quantum-classical algorithms based on real-time evolution are promising methods for solving a particularly important model in the search for spin liquids, the antiferromagnetic Heisenberg model on the two-dimensional kagome lattice. We show how to construct efficient quantum circuits to implement time evolution for the model and to evaluate key observables on the quantum computer, and we argue that the method has favorable scaling with increasing system size. We then restrict to a 12-spin star plaquette from the kagome lattice and a related 8-spin system, and we give an empirical demonstration on these small systems that the hybrid algorithms can efficiently find the ground state energy and the magnetization curve. For these demonstrations, we use four levels of approximation: exact state vectors, exact state vectors with statistical noise from sampling, noisy classical emulators, and (for the 8-spin system only) real quantum hardware, specifically the Quantinuum H1-1 processor; for the noisy simulations and hardware demonstration, we also employ error mitigation strategies based on the symmetries of the Hamiltonian. Our results strongly suggest that these hybrid algorithms present a promising direction for studying quantum spin liquids and more generally for resolving important unsolved problems in condensed matter theory and beyond.

97 MATHEMATICS AND COMPUTING↗

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

Nominal and adversarial synthetic PMU data for standard IEEE test systems

GridSTAGE (Spatio-Temporal Adversarial scenario GEneration) is a framework for the simulation of adversarial scenarios and the generation of multivariate spatio-temporal data in cyber-physical systems. GridSTAGE is developed based on Matlab and leverages Power System Toolbox (PST) where the evolution of the power network is governed by nonlinear differential equations. Using GridSTAGE, one can create several event scenarios that correspond to several operating states of the power network by enabling or disabling any of the following: faults, AGC control, PSS control, exciter control, load changes, generation changes, and different types of cyber-attacks. Standard IEEE bus system data is used to define the power system environment. GridSTAGE emulates the data from PMU and SCADA sensors. The rate of frequency and location of the sensors can be adjusted as well. Detailed instructions on generating data scenarios with different system topologies, attack characteristics, load characteristics, sensor configuration, control parameters are available in the Github repository - https://github.com/pnnl/GridSTAGE. There is no existing adversarial data-generation framework that can incorporate several attack characteristics and yield adversarial PMU data. The GridSTAGE framework currently supports simulation of False Data Injection attacks (such as a ramp, step, random, trapezoidal, multiplicative, replay, freezing) and Denial of Service attacks (such as time-delay, packet-loss) on PMU data. Furthermore, it supports generating spatio-temporal time-series data corresponding to several random load changes across the network or corresponding to several generation changes. A Koopman mode decomposition (KMD) based algorithm to detect and identify the false data attacks in real-time is proposed in https://ieeexplore.ieee.org/document/9303022. Machine learning-based predictive models are developed to capture the dynamics of the underlying power system with a high level of accuracy under various operating conditions for IEEE 68 bus system. The corresponding machine learning models are available at https://github.com/pnnl/grid_prediction.

99 GENERAL AND MISCELLANEOUS↗

Intrinsic value of the international benchmark projects, ICSBEP and IRPhEP, for advanced reactor development

The International Criticality Safety Benchmark Evaluation Project, a sanctioned program under the auspices of the Nuclear Energy Agency of the Organisation for Economic Co-operation and Development, has been a highly successful and productive collaboration, now encompassing over 5,000 evaluated experimental benchmarks trusted and relied upon throughout the international nuclear communities. The success of this project led to the development of the International Reactor Physics Experiment Evaluation Project, which is dedicated to the evaluation of benchmark experiment data to sustain current and future reactor physics validation needs. These exemplary programs, and their widely utilized handbooks, serve as gold standards to which other databases strive to emulate. The purpose of the two projects is to preserve modern and legacy experimental data and evaluate it in a standardized handbook format to provide quality benchmarks to support modern and future criticality safety and reactor physics validation. These two projects have often served as the mechanism through which historic and modern neutronics experiments are evaluated and shared across international borders, to best provide unique, high-quality peer-reviewed, and often otherwise unavailable, benchmark data. The contents of these handbooks are utilized not only in validating criticality safety, reactor physics, and advanced reactor calculations, but are used to validate neutronics calculations and nuclear data for most other nuclear applications. This paper discusses both international programs and available content to enable advanced reactor design validation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Synthesis of Si-Fe Chondrule-like Dust Analogues in RF Discharge Plasmas

Chondrules are tiny particles that occur in stony meteorites and are considered as the building blocks of early asteroids and planets. It is believed that they were formed by the fast heating of the dust in the solar nebula. To date, there is no lab-scale experimental study of the formation of chondrules from the initial gas phase precursors following fast heating and crystallisation. The motivation of this work is a pre-trial study of the formation of chnodrule-like particles. The formation of meteorites in the space environment is associated with the aggregation of small particles or molecular clouds under the influence of shock waves or high-energy gas discharges in the solar nebula. In this work, the properties of product formation at the nanoscale-level were investigated using different feedstock materials which are the dominant elements in the meteorite. The structural and morphological properties of the synthesised Si-Fe nanomaterials were analysed by scanning/transmission electron microscopy (SEM/TEM), and chemical composition was analysed by X-ray energy-dispersive spectroscopy (EDS). The identification of crystalline phases was carried out by X-ray diffraction (XRD), whereas the presence of an Fe-Si system in the synthesised particles was demonstrated by Mössbauer spectroscopy. The obtained materials were exposed to the relatively high-energy pulsed plasma beam on the substrate with the aim to emulate the possible fast heating and melting of the formed nanoparticles. The formation steps of growing synthetic (engineered) chondro-like particles and nanostructures in laboratory conditions is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of Envelope- and Machine Learning-Based Electrical Fault Type Detection Algorithms for Electrical Distribution Grids

This study introduces envelope- and machine learning (ML)-based electrical fault type detection algorithms for electrical distribution grids, advancing beyond traditional logic-based methods. The proposed detection model involves three stages: anomaly area detection, ML-based fault presence detection, and ML-based fault type detection. Initially, an envelope-based detector identifying the anomaly region was improved to handle noisier power grid signals from meters. The second stage acts as a switch, detecting the presence of a fault among four classes: normal, motor, switching, and fault. Finally, if a fault is detected, the third stage identifies specific fault types. This study explored various feature extraction methods and evaluated different ML algorithms to maximize prediction accuracy. The performance of the proposed algorithms is tested in an emulated software–hardware electrical grid testbed using different sample rate meters/relays, such as SEL735, SEL421, SEL734, SEL700GT, and SEL351S near and far from an inverter-based photovoltaic array farm. The performance outcomes demonstrate the proposed model’s robustness and accuracy under realistic conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Review of Design Elements within Power Infrastructure Cyber–Physical Test Beds as Threat Analysis Environments

Cyber–physical systems (CPSs) are an integral part of modern society; thus, enhancing these systems’ reliability and resilience is paramount. Cyber–physical testbeds (CPTs) are a safe way to test and explore the interplay between the cyber and physical domains and to cost-effectively enhance the reliability and resilience of CPSs. Here a review of CPT elements, broken down into physical components (simulators, emulators, and physical hardware), soft components (communication protocols, network timing protocols), and user interfaces (visualization-dashboard design considerations) is presented. Various methods used to validate CPS performance are reviewed and evaluated for potential applications in CPT performance validation. Last, initial simulated results for a CPT design, based on the IEEE 33 bus system, are presented, along with a brief discussion on how model-based testing and fault–injection-based testing (using scaling and ramp-type attacks) may be used to help validate CPT performance.

14 SOLAR ENERGY↗

A Methodology for Robust Load Reduction in Wind Turbine Blades Using Flow Control Devices

Decades of wind turbine research, development and installation have demonstrated reductions in levelized cost of energy (LCOE) resulting from turbines with larger rotor diameters and increased hub heights. Further reductions in LCOE by up-scaling turbine size can be challenged by practical limitations such as the square-cube law: where the power scales with the square of the blade length and the added mass scales with the volume (the cube). Active blade load control can disrupt this trend, allowing longer blades with less mass. This paper presents the details of the development of a robust load control system to reduce blade fatigue loads. The control system, which we coined sectional lift control or SLC, uses a lift actuator model to emulate an active flow control device. The main contributions of this paper are: (1) Methodology for SLC design to reduce dynamic blade root moments in a neighborhood of the rotor angular frequency (1P). (2) Analysis and numerical evidence supporting the use of a single robust SLC for all wind speeds, without the need for scheduling on wind speed or readily available measurements such as collective pitch or generator angular speed. (3) Intuition and numerical evidence to demonstrate that the SLC and the turbine controller do not interact. (4) Evaluation of the SLC using a full suite of fatigue and turbine performance metrics.

17 WIND ENERGY↗

Dynamical System Scaling of a Thermocline Thermal Storage System in the Thermal Energy Distribution System (TEDS) Facility

The purpose of this study was to develop a process to convert input signals from one facility into another by reflecting geometric and environmental settings. The Dynamic Energy Transport and Integration Laboratory (DETAIL) is one facility in development that aims to emulate the daily interactions among power production industry systems and be capable of receiving real-time data from those systems as inputs. To convert signals and ensure that the temporal sequences and magnitudes reflect the laboratory settings, the ability to scale and project data is essential. To demonstrate this ability, Dynamical System Scaling (a methodology that enables systems to scale and project or extrapolate datasets to desired environments while conserving the observed transient behavior based on first principles) was applied to DETAIL’s thermocline thermal storage system in the Thermal Energy Distribution System. The thermocline system was successfully scaled and a test case was conducted to generate a doubly accelerated energy charge and discharge in reference to past experimental data from the facility. The accelerated data were determined as able to conserve the amount of energy stored and the associated test boundary conditions were charge line maximum temperature, charge line velocity, and thermocline maximum temperature at 354 °C, 0.458 m/s, and 418 °C, respectively. The research results represented a case that required signals to be accelerated without altering the stored energy.

Yoshiura, Ramon (ORCID:0000000319781313)↗

Simulation of the Multi-Wake Evolution of Two Sandia National Labs/National Rotor Testbed Turbines Operating in a Tandem Layout

The future of wind power systems deployment is in the form of wind farms comprised of scores of such large turbines, most likely at offshore locations. Individual turbines have grown in span from a few tens of meters to today’s large turbines with rotor diameters that dwarf even the largest commercial aircraft. These massive dynamical systems present unique challenges at scales unparalleled in prior applications of wind science research. Fundamental to this effort is the understanding of the wind turbine wake and its evolution. Furthermore, the optimization of the entire wind farm depends on the evolution of the wakes of different turbines and their interactions within the wind farm. In this article, we use the capabilities of the Common ODE Framework (CODEF) model for the analysis of the effects of wake–rotor and wake-to-wake interactions between two turbines situated in a tandem layout fully and partially aligned with the incoming wind. These experiments were conducted in the context of a research project supported by the National Rotor Testbed (NRT) program of Sandia National Labs (SNL). Results are presented for a layout which emulates the turbine interspace and relative turbine emplacement found at SNL’s Scaled Wind Technologies Facility (SWiFT), located in Lubbock, Texas. The evolution of the twin-wake interaction generates a very rich series of secondary transitions in the vortex structure of the combined wake. These ultimately affect the wake’s axial velocity patterns, altering the position, number, intensity, and shape of localized velocity-deficit zones in the wake’s cross-section. This complex distribution of axial velocity patterns has the capacity to substantially affect the power output, peak loads, fatigue damage, and aeroelastic stability of turbines located in subsequent rows downstream on the farm.

Baruah, Apurva (ORCID:0000000252354068)↗

Power Production, Inter- and Intra-Array Wake Losses from the U.S. East Coast Offshore Wind Energy Lease Areas

There is an urgent need to develop accurate predictions of power production, wake losses and array–array interactions from multi-GW offshore wind farms in order to enable developments that maximize power benefits, minimize levelized cost of energy and reduce investment uncertainty. New, climatologically representative simulations with the Weather Research and Forecasting (WRF) model are presented and analyzed to address these research needs with a specific focus on offshore wind energy lease areas along the U.S. east coast. These, uniquely detailed, simulations are designed to quantify important sources of wake-loss projection uncertainty. They sample across different wind turbine deployment scenarios and thus span the range of plausible installed capacity densities (ICDs) and also include two wind farm parameterizations (WFPs; Fitch and explicit wake parameterization (EWP)) and consider the precise WRF model release used. System-wide mean capacity factors for ICDs of 3.5 to 6.0 MWkm−2 range from 39 to 45% based on output from Fitch and 50 to 55% from EWP. Wake losses are 27–37% (Fitch) and 11–19% (EWP). The discrepancy in CF and wake losses from the two WFPs derives from two linked effects. First, EWP generates a weaker ‘deep array effect’ within the largest wind farm cluster (area of 3675 km2), though both parameterizations indicate substantial within-array wake losses. If 15 MW wind turbines are deployed at an ICD of 6 MWkm−2 the most heavily waked wind turbines generate an average of only 32–35% of the power of those that experience the freestream (undisturbed) flow. Nevertheless, there is no evidence for saturation of the resource. The wind power density (electrical power generation per unit of surface area) increases with ICD and lies between 2 and 3 Wm−2. Second, EWP also systematically generates smaller whole wind farm wakes. Sampling across all offshore wind energy lease areas and the range of ICD considered, the whole wind farm wake extent for a velocity deficit of 5% is 1.18 to 1.38 times larger in simulations with Fitch. Over three-quarters of the variability in normalized wake extents is attributable to variations in freestream wind speeds, turbulent kinetic energy and boundary layer depth. These dependencies on meteorological parameters allow for the development of computationally efficient emulators of wake extents from Fitch and EWP.

17 WIND ENERGY↗

Exponential Backoff and Its Security Implications for Safety-Critical OT Protocols over TCP/IP Networks

The convergence of Operational Technology (OT) and Information Technology (IT) networks has become increasingly prevalent with the growth of Industrial Internet of Things (IIoT) applications. This shift, while enabling enhanced automation, remote monitoring, and data sharing, also introduces new challenges related to communication latency and cybersecurity. Oftentimes, legacy OT protocols were adapted to the TCP/IP stack without an extensive review of the ramifications to their robustness, performance, or safety objectives. To further accommodate the IT/OT convergence, protocol gateways were introduced to facilitate the migration from serial protocols to TCP/IP protocol stacks within modern IT/OT infrastructure. However, they often introduce additional vulnerabilities by exposing traditionally isolated protocols to external threats. This study investigates the security and reliability implications of migrating serial protocols to TCP/IP stacks and the impact of protocol gateways, utilizing two widely used OT protocols: Modbus TCP and DNP3. Our protocol analysis finds a significant safety-critical vulnerability resulting from this migration, and our subsequent tests clearly demonstrate its presence and impact. A multi-tiered testbed, consisting of both physical and emulated components, is used to evaluate protocol performance and the effects of device-specific implementation flaws. Through this analysis of specifications and behaviors during communication interruptions, we identify critical differences in fault handling and the impact on time-sensitive data delivery. The findings highlight how reliance on lower-level IT protocols can undermine OT system resilience, and they inform the development of mitigation strategies to enhance the robustness of industrial communication networks.

DNP3↗