Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Exponential time differencing for the tracer equations appearing in primitive equation ocean models

The tracer equations are part of the primitive equations used in ocean modeling and describe the transport of tracers, such as temperature, salinity or chemicals, in the ocean. Depending on the number of tracers considered, several equations may be added to and coupled to the dynamics system. In many relevant situations, the time-step requirements of explicit methods imposed by the transport and mixing in the vertical direction are more restrictive than those for the horizontal, and this may cause the need to use very small time steps if a fully explicit method is employed. To overcome this issue, we propose an exponential time differencing (ETD) solver where the vertical terms (transport and diffusion) are treated with a matrix exponential, whereas the horizontal terms are dealt with in an explicit way. In this work, we investigate numerically the computational speed-ups that can be obtained over other semi-implicit methods, and we analyze the advantages of the method in the case of multiple tracers.

42 ENGINEERING↗
DOE OSTI2021

Parallel exponential time differencing methods for geophysical flow simulations

Two ocean models are considered for geophysical flow simulations: the multilayer shallow water equations and the multilayer primitive equations. For the former, we investigate the parallel performance of exponential time differencing (ETD) methods, including exponential Rosenbrock–Euler, ETD2wave, and B-ETD2wave. For the latter, we take advantage of the splitting of barotropic and baroclinic modes and propose a new two-level method in which an ETD method is applied to solve the fast barotropic mode. Furthermore, these methods could improve the computational efficiency of numerical simulations because ETD methods allow for much larger time step sizes than traditional explicit time-stepping techniques that are commonly used in existing computational ocean models. Several standard benchmark tests for ocean modeling are performed and comparison of the numerical results demonstrates a great potential of applying the parallel ETD methods for simulating real-world geophysical flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ignition delay time and speciation of dibutyl ether at high pressures

Dibutyl ether (DBE) is a promising biofuel due to its high cetane number (~ 100) and high volumetric energy density (31.6 MJ/L). It could either be used directly in compression ignition engines or blended with other conventional or renewable fuels. Oxidation and pyrolysis kinetics of DBE are not well known, particularly at high pressures. In this work, we have experimentally investigated the chemical kinetics of DBE in three domains: (a) ignition delay time measurements in a rapid compression machine over T = 550–650 K, P = 10, 20, 40 bar, φ = 0.5, 1; (b) ignition delay time measurements in a shock tube over T = 900–1300 K, P = 20, 40 bar, φ = 0.5, 1; (c) laser-based carbon monoxide speciation measurements in a shock tube during DBE pyrolysis and oxidation over T = 1100–1400 K, P = 20 bar. Pressure time-histories measured in RCM experiments exhibited unique 3-stage and 4-stage ignition behavior predominantly at fuel-lean conditions. Experimental data were compared with the predictions of two recent chemical kinetic models of DBE. Sensitivity analyses were carried out to identify key reactions which may have caused the discrepancy between experiments and simulations. We found that the rate of decomposition of DBE may need to be revisited to improve the oxidative and pyrolytic predictions of DBE kinetic model.

09 BIOMASS FUELS↗
DOE OSTI2021

Skeletal model reduction with forced optimally time dependent modes

Skeletal model reduction based on local sensitivity analysis of time dependent systems is presented in which sensitivities are modeled by forced optimally time dependent (f-OTD) modes. The f-OTD factorizes the sensitivity coefficient matrix into a compressed format as the product of two skinny matrices, i.e. f-OTD modes and f-OTD coefficients. The modes create a low-dimensional, time dependent, orthonormal basis which capture the directions of the phase space associated with most dominant sensitivities. These directions highlight the instantaneous active species, and reaction paths. Evolution equations for the f-OTD modes and coefficients are derived, and the implementation of f-OTD for skeletal reduction is described. For demonstration, skeletal reduction is conducted of the constant pressure ethylene-air burning in a zero-dimensional reactor, and new reduced models are generated. The laminar flame speed, the ignition delay, and the extinction curve as predicted by the models are compared against some existing skeletal models in literature for the same detailed model. Finally, the results demonstrate the capability of f-OTD to eliminate unimportant reactions and species in a systematic, efficient and accurate manner.

42 ENGINEERING↗
DOE OSTI2026

An adaptive model-free robotic force control strategy for hydrodynamic real-time hybrid simulation of floating offshore wind turbines

Real-time hybrid simulation (RTHS) - a cyber-physical testing approach - promises to enhance the simulation fidelity of the model-scale experiments used to prototype floating offshore wind turbines (FOWTs). In hydrodynamic RTHS (hydro-RTHS), actuators emulate aerodynamic forces on model-scale FOWT specimens subjected to physical waves in a hydrodynamic laboratory. Robotic arms are promising candidates for actuation in hydro-RTHS due to their compact multi-degree-of-freedom (DOF) capabilities. Unlike classical RTHS for seismic applications, which typically relies on displacement control, hydro-RTHS requires 6-DOF force control on newly designed floating prototypes in a model-scale setting, which presents significant challenges, including modeling uncertainties, directional asymmetry, configuration drift, bandwidth limitations, and time-varying delays. To mitigate these constraints without extensive pre-test calibration, this study proposes an adaptive model-free robotic force control strategy that combines task-space explicit force control with a secondary joint-space pose-keeping task. The Adaptive Feedforward Compensator (AFC) is integrated into the force control loop to compensate for time-varying delay. Experimental testing was conducted using a Franka Emika Panda robotic arm with a 1:50 scale FOWT specimen under operational wind and wave conditions. Results demonstrate stable and consistent 6-DOF force tracking. Effective delay compensation was observed, with low-frequency delay reductions ranging from 71.4% to 91.8% and improvements in low-frequency surge force tracking of 25.0% to 52.1%. This study enhances robotic actuation performance in hydro-RTHS and introduces a force control strategy that supports reliable robotic operation in uncertain floating environments. Future work will explore disturbance-observer mechanisms to further enhance wave rejection capabilities under extreme wind and wave conditions.

17 WIND ENERGY↗
DOE OSTI2021

Keras2c: A library for converting Keras neural networks to real-time compatible C

With the growth of machine learning models and neural networks in measurement and control systems comes the need to deploy these models in a way that is compatible with existing systems. Existing options for deploying neural networks either introduce very high latency, require expensive and time consuming work to integrate into existing code bases, or only support a very limited subset of model types. We have therefore developed a new method called Keras2c, which is a simple library for converting Keras/TensorFlow neural network models into real-time compatible C code. It supports a wide range of Keras layers and model types including multidimensional convolutions, recurrent layers, multi-input/output models, and shared layers. Keras2c re-implements the core components of Keras/TensorFlow required for predictive forward passes through neural networks in pure C, relying only on standard library functions considered safe for real-time use. The core functionality consists of ~ 1500 lines of code, making it lightweight and easy to integrate into existing codebases. Keras2c has been successfully tested in experiments and is currently in use on the plasma control system at the DIII-D National Fusion Facility at General Atomics in San Diego.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Coupled aero-hydro-geotech real-time hybrid simulation of offshore wind turbine monopile structures

Real-time hybrid simulation (RTHS) divides a structural system into an analytical and experimental substructure. The former is based on a well-established analytical model while the latter consists of a physical model in the laboratory, for which there is not a well-established analytical model. This paper extends real-time hybrid simulation to monopile-type Offshore Wind Turbines (OWTs) to enable the investigation of their behavior considering the response of pile foundations under operational and more severe conditions. The embedded foundation and surrounding soil of the OWT are modeled physically in a soil box in the laboratory while the remaining parts of the system and loading are modeled analytically. The program OpenFAST, developed by the National Renewable Energy Laboratory (NREL), is linked to the RTHS coordinator to determine the hydrodynamic and aerodynamic loads acting on the OWT, along with modeling the dynamics of the electric power generation equipment and associated controller for the OWT. The RTHS framework along with its initial implementation and validation are described in this paper. RTHSs of a 5 MW OWT subjected to operational and more severe conditions are performed to experimentally validate the framework. The framework offers a realistic approach to investigate the behavior of OWT structures supported on monopiles. Furthermore, this approach accounts for the coupled response of the OWT structure with its foundation, while experimentally capturing the nonlinearities of the soil-foundation interaction in real-time.

17 WIND ENERGY↗
DOE OSTI2020

High-pressure ignition delay time measurements of a four-component gasoline surrogate and its high-level blends with ethanol and methyl acetate

Ignition delay times of a gasoline surrogate (iso-octane/toluene/n-heptane/1-hexene at 55/25/15/5% by liquid volume, developed in McCormick et al. 2017) with high-level bioblendstock-gasoline surrogate blends (50% and 85% biofuel by liquid volume of each ethanol and methyl acetate) were collected behind reflected shock waves. Post-reflected-shock temperatures ranged from 968 to 1361 K at pressures of about 4, 10, 25, and 50 atm. Data were collected for real fuel–air mixtures at fuel lean and stoichiometric conditions (φ = 0.5, 1.0) with focus on the more-practical, fuel-lean conditions. Ignition delay times were measured from OH* chemiluminescence around 307 nm using an endwall diagnostic. The RON and MON of the surrogate were 90.3 and 84.7, respectively. Ethanol was chosen due to its wide use in flex-fuel vehicles and availability, while also providing an increased octane rating. Methyl acetate was chosen for its especially high octane rating to investigate the effect of an extreme case. To validate the gasoline surrogate, the data are compared to real gasoline (RD387) and several gasoline surrogate experiments from the literature. Using the wide range of pressures studied, the gasoline surrogate’s pressure dependence was quantified to account for test-to-test variations using regression analysis. Similarly, a global correlation for gasoline and its surrogates was developed using all available data from the literature. Two modern chemical kinetics models targeting gasoline and its surrogates are compared to the ignition delay time measurements. These new high-pressure, high-bioblendstock concentration tests provide required chemical kinetic data for optimizing fuel and engine design.

42 ENGINEERING↗
DOE OSTI2024

Utilization of an Advanced Sensor network to determine fuel heating value and Real-Time net unit heat rate during transient operation

Coal-fired utility boilers are being increasingly used as variable electricity generation to resolve the imbalance in the energy market from the expansion of intermittent renewable energy. The frequent transient operation required to meet residual energy demand has created a challenge for coal-fired units to operate efficiently. This work utilizes an advanced sensor network (ASN) to calculate net unit heat rate (NUHR) of a coal-fired boiler in real time through combustion calculations and statistical correlations to provide the tools for optimizing dynamic operation. Real-time heating values that were necessary to determine fuel input energy to calculate accurate NUHR were found using both fundamental and data-driven methods. Real-time NUHR shows distinct shifts that reflect changes in process conditions that will improve the ability to optimize transient operation. Data-driven heating value correlations had 24% lower root mean square error (RMSE) than the fundamental combustion calculation approach when compared to daily retrospective proximate analysis. Furthermore, the data-driven method RMSE improved by 7% with the inclusion of ASN data. Future work is to validate by comparing unit performance with and without the inclusion of NUHR as a control parameter for the dynamic neural network.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Extended real-time voltage instability identification method based on synchronized phasor measurements

This paper presents an extended adaptive approach designed to accurately estimate the Thévenin equivalent parameters using phasor measurements at a given bus for measurements lying in any of four quadrants of the PQ-plane. The improvement is achieved by using a new condition to properly update the estimated parameters after an initial guess. Based on an adaptive philosophy, the proposed approach can correctly account for the intrinsic nonlinearities of a power system, can provide a real-time estimation of Thévenin parameters, and does not require network topology knowledge. The method is validated using the Kundur 2-area system, showing estimation improvements compared to the current adaptive approach and the classical recursive least-squares method. The proposed approach is able to estimate both sides of the system with respect to the measurement bus. In addition, a data-driven voltage stability index is developed. To illustrate the performance of the proposed approach in a larger power system, a voltage stability assessment is carried out on the IEEE 39-bus system, considering the action of overexcitation limiters of generators and nonlinear loads. The proposed approach is suitable for applications that require an accurate Thévenin equivalent estimation in real-time, such as for voltage stability assessment. The new approach provides a reliable tool for the system operators to make proper and timely decisions.

42 ENGINEERING↗
DOE OSTI2025

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Computer-controlled finishing via dynamically constraint position-velocity-time scheduler

In a Computer Numerical Controlled (CNC) finishing process, the target material removal from an optical surface is guided by the convolution between the influence function of a machine tool and its dwell time at certain points over the surface. To reduce dynamics stressing and increase machining efficiency, the dwell time must be converted to varying velocities, which are the actual inputs to the machine tool controller. Conventionally, the conversion assumed constant acceleration and relied on linear motion interpolation, which caused discontinuities in velocities. This unsmooth motion affects the material removal distribution, and, thus, the accuracy of the finished surface shape. Many modern CNC machines support the smoother, cubic-polynomial interpolated Position-Velocity-Time (PVT) motion mode; however, the conventional scheduler may fail to provide suitable velocities for the PVT. Here in this study answers this challenge by proposing a novel PVT-based velocity scheduler that achieves smooth motion while considering CNC dynamic limits. Firstly, the principle of the PVT is explained, and the PVT-based velocity scheduler is formulated. Secondly, a quadratic programming is used to optimize the velocities by imposing the CNC dynamic constraints and the C 1 continuities (zeroth and first derivatives are continuous) simultaneously. Thirdly, the smoothness and accuracy of the scheduled velocities are studied on different kinds of tool paths via simulation. Finally, a sub-0.3 nm level surface finishing experiment using ion beam figuring is demonstrated to verify the feasibility of the proposed method. The PVT-based scheduler and simulator code is open-sourced.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Machine learning-based real-time monitoring system for smart connected worker to improve energy efficiency

Recent advances in machine learning and computer vision brought to light technologies and algorithms that serve as new opportunities for creating intelligent and efficient manufacturing systems. In this study, the real-time monitoring system of manufacturing workflow for the Smart Connected Worker (SCW) is developed for the small and medium-sized manufacturers (SMMs), which integrates state-of-the-art machine learning techniques with the workplace scenarios of advanced manufacturing systems. Specifically, object detection and text recognition models are investigated and adopted to ameliorate the labor-intensive machine state monitoring process, while artificial neural networks are introduced to enable real-time energy disaggregation for further optimization. The developed system achieved efficient supervision and accurate information analysis in real-time for prolonged working conditions, which could effectively reduce the cost related to human labor, as well as provide an affordable solution for SMMs. The competent experiment results also demonstrated the feasibility and effectiveness of integrating machine learning technologies into the realm of advanced manufacturing systems.

42 ENGINEERING↗
DOE OSTI2022

Optimal checkpointing for adjoint multistage time-stepping schemes

Here, we consider checkpointing strategies that minimize the number of recomputations needed when performing discrete adjoint computations using multistage time-stepping schemes that require computing several substeps within one complete time step. Specifically, we propose two algorithms that can generate optimal checkpoint-ing schedules under weak assumptions. The first is an extension of the seminal Revolve algorithm adapted to multistage schemes. The second algorithm, named CAMS, is developed based on dynamic programming, and it requires the least number of recomputations when compared with other algorithms. The CAMS algorithm is made publicly available in a library with bindings to C and Python. Numerical results show that the proposed algorithms can deliver up to two times the speedup compared with that of classical Revolve. Moreover, we discuss the utilization of the CAMS library in mature scientific computing libraries and demonstrate the ease of using it in an adjoint workflow. The proposed algorithms have been adopted by the PETSc TSAdjoint library. Their performance has been demonstrated with a large-scale PDE-constrained optimization problem on a leadership-class supercomputer. This work is a significant extension of the authors' conference paper.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Identifying electrochemical processes by distribution of relaxation times in proton exchange membrane electrolyzers

Distribution of relaxation time (DRT) is used to interpret electrochemical impedance spectroscopy (EIS) for proton exchange membrane (PEM) water electrolyzers, with an attempt to separate overlapped relaxation processes in Nyquist plots. By varying operating conditions and catalyst loadings, four main relaxation peaks arising from EIS can be identified and successfully separated from low to high frequencies as (P1) mass transport, (P2) oxygen evolution reaction kinetics, (P3) reaction kinetics (with faster time constant than P2), and (P4) ionic transport. Here, the shape, height, and frequency of the DRT peaks change with different membrane electrode assembly (MEA) configurations. Electron microscopy reveals distinct features from the cross-sectioned MEAs which verify critical DRT results in that increasing the iridium (Ir)-anode loading from 0.2 mgIr/cm 2 to 1.5 mgIr/cm 2 reduces kinetic losses due to higher site-access; a thick and compacted anode, however, also triggers higher ohmic resistances from membrane/catalyst layer hydration and increases transport losses due to longer ionomer pathways. DRT provides higher resolution to EIS for deconvoluting processes with different relaxation times and the quantification of DRT peaks improves the accounting of total losses from each process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Creating an Advanced Sensor Network to calculate real-time, mass-weighted flue gas composition and air heater leakage of a coal-fired utility boiler under dynamic operating conditions

Utilization of renewable energy sources to minimize the environmental impact of energy production has changed the way utility boilers operate, requiring frequent load cycling between full load and partial loads as low as 30%. Dynamic operation of coal-fired utility boilers significantly reduces boiler efficiency when compared to steady state at full load. Data-driven plant optimization has shown success with coal-fired utility boilers under dynamic operating conditions. The purpose of this work was to create an Advanced Sensor Network (ASN) to provide more extensive real-time data to inform dynamic plant optimization of Net Unit Heat Rate (NUHR). The ASN consists of gas sampling grids in the convective pass of the boiler and downstream of the air heater. These sampling grids allow for quantification of spatial variation of flue gas within the boiler and calculation of mass-weighted composition of flue gas through the combination of composition, velocity, and temperature measurements. The comparison of O 2 between the inlet and outlet of the air heater is used to calculate air leakage in real time. Flue gas composition and air heater leakage are both important factors in boiler efficiency and NUHR. Further, the results of this work support the value of mass-weighted averages for determining flue gas composition accurately. The measurements from the ASN show increased composition stratification during dynamic operation, with an average standard deviation 38% higher than observed during steady-state operation. Air heater leakage was also observed to increase from 2.8% to 5.1% following a load change. Prior to the installation of the ASN, these data would not have been available for dynamic control. These real-time data will be leveraged to calculate and optimize for NUHR during dynamic operation in future work.

42 ENGINEERING↗
DOE OSTI2021

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗
DOE OSTI2022

BoxScore - A real-time beam-diagnosis program for the CAEN digitizer x730 series

BoxScore is a real-time beam diagnosis and monitoring program for the CAEN x730 series digitizer, developed for the ATLAS in-flight system at Argonne National Laboratory. The CAEN x730 series digitizer, with built-in Digital Pulse Processing for the Pulse-Height-Analysis, digitizes the input signal in real-time and processes it using a trapezoidal filter. BoxScore reads the digitizer’s buffer directly, builds and saves events to local files, plots histograms for particle identification, and outputs the rates of selected isotopes every second. Implementation of BoxScore has shortened the time needed for in-flight beam-tuning and has potential applications for other nuclear physics experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗