Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIGITAL SIMULATION”

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

RFI in hybrid loops - Simulation and experimental results.

A digital simulation of an imperfect second-order hybrid phase-locked loop (HPLL) operating in radio frequency interference (RFI) is described. Its performance is characterized in terms of phase error variance and phase error probability density function (PDF). Monte-Carlo simulation is used to show that the HPLL can be superior to the conventional phase-locked loops in RFI backgrounds when minimum phase error variance is the goodness criterion. Similar experimentally obtained data are given in support of the simulation data.

Ziemer, R. E.↗

Non-linear dynamics of a spur gear pair

The backlash nonlinearity excited primarily by transmission error between spur gear pairs is studied for both external and internal excitations. The digital simulation technique and the method of harmonic balance are used to develop steady state solutions for the internal sinuosidal excitations. The analytic predictions agreed well with available experimental data. Digital simulation is used to observe that at the chaotic and subharmonic resonances may exist in a gear pair depending on the mean or design load, mean to alternating force ratio, damping, and backlash.

Kahraman, A.↗

Thermalization and criticality on an analogue–digital quantum simulator

Abstract Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators 1,2 . Unlocking the full potential of such systems towards this goal requires flexible initial state preparation, precise time evolution and extensive probes for final state characterization. Here we present a quantum simulator comprising 69 superconducting qubits that supports both universal quantum gates and high-fidelity analogue evolution, with performance beyond the reach of classical simulation in cross-entropy benchmarking experiments. This hybrid platform features more versatile measurement capabilities compared with analogue-only simulators, which we leverage here to reveal a coarsening-induced breakdown of Kibble–Zurek scaling predictions 3 in theXYmodel, as well as signatures of the classical Kosterlitz–Thouless phase transition 4 . Moreover, the digital gates enable precise energy control, allowing us to study the effects of the eigenstate thermalization hypothesis 5–7 in targeted parts of the eigenspectrum. We also demonstrate digital preparation of pairwise-entangled dimer states, and image the transport of energy and vorticity during subsequent thermalization in analogue evolution. These results establish the efficacy of superconducting analogue–digital quantum processors for preparing states across many-body spectra and unveiling their thermalization dynamics.

Science & Technology - Other Topics↗

Frequency Domain Reflectometry (FDR) Simulation Techniques for Digital Twin Representation of an Electrical Cable

Simulation of cable system response to frequency domain reflectometry (FDR) tests can be instrumental to understanding these tests and the nature and influence of various cable anomalies on test signatures. Reflectometry simulations are based upon a finite element representation of cable conductors and insulation to produce an S-parameter at each evaluated frequency. The aggregate collection of cable model S-parameters can simulate the influence on a test signal injected into a physical cable. Such an approach was undertaken in this work to produce a digital twin simulation of a low-voltage electrical cable. The electrical cable digital twin examined the influence of test simulation parameters and the relative influence of cable anomalies, including thermal aging, water or moisture exposure, water or moisture ingress, and other anomalies. The digital twin in this work included modeling of the conductors, insulation, jacket, and surrounding environment (air, water, etc.). The digital twin could be expanded to include cable bends, junctions and splices, branch or T systems, and termination impedances of motors or instruments. Observations and conclusions of this work include: 1. Fully 3D digital twin simulation of an electrical cable using an FDR approach is possible. However, there are tradeoffs between simulation fidelity and solution time, which must be balanced to ensure the simulation solves in an adequate amount of time (e.g., less than 20 minutes). Simulation parameters to balance include frequency bandwidth, number of frequencies, mesh density, connection impedance, and permittivity tolerance. 2. The digital twin simulation can explain FDR sensitivity to various cable anomalies, including entry and exit from an oven or water bath. 3. The digital twin simulation FDR response attenuates with distance along the cable and is further affected by the frequency bandwidth, which is similar to that observed with physical measurements. 4. The resolution of the digital twin FDR peaks increased with increasing bandwidth and with increasing number of frequencies, again similar to physical measurements. 5. The presence of multiple anomalies in the digital twin does not substantially attenuate the FDR response to anomalies located beyond the first encountered anomaly and impedance mismatch. 6. Spectral variation of the permittivity did not have a significant effect on the FDR response compared to a fixed nominal value. 7. Extension of the digital twin to 1000 ft still allowed for detection of distal anomalies near the far end of the electrical cable from the instrument connection point. 8. The ARENA test bed facilitates efficient NDE evaluations of well understood cable anomalies with various NDE methods without risking actual plant damage.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Floating-point system quantization errors in digital control systems

The results are reported of research into the effects on system operation of signal quantization in a digital control system. The investigation considered digital controllers (filters) operating in floating-point arithmetic in either open-loop or closed-loop systems. An error analysis technique is developed, and is implemented by a digital computer program that is based on a digital simulation of the system. As an output the program gives the programing form required for minimum system quantization errors (either maximum of rms errors), and the maximum and rms errors that appear in the system output for a given bit configuration. The program can be integrated into existing digital simulations of a system.

Phillips, C. L.↗

Augmented Kalman filter with a reduced mechanical model to estimate tower loads on a land-based wind turbine: a step towards digital-twin simulations

This article presents an application of the Kalman filtering technique to estimate loads on a wind turbine. The approach combines a mechanical model and a set of measurements to estimate signals that are not available in the measurements, such as wind speed, thrust, tower position, and tower loads. The model is severalfold faster than real time and is intended to be run online, for instance, to evaluate real-time fatigue life consumption of a field turbine using a digital twin, perform condition monitoring, or assess loads for dedicated control strategies. The mechanical model is built using a Rayleigh–Ritz approach and a set of joint coordinates. We present a general method and illustrate it using a 2-degrees-of-freedom (DOF) model of a wind turbine and using rotor speed, generator torque, pitch, and tower-top acceleration as measurement signals. The different components of the model are tested individually. The overall method is evaluated by computing the errors in estimated tower-bottom-equivalent moment from a set of simulations. From this preliminary study, it appears that the tower-bottom-equivalent moment is obtained with about 10 % accuracy. The limitation of the model and the required steps forward are discussed.

17 WIND ENERGY↗

Design and Cosimulation of Hierarchical Architecture for Demand Response Control and Coordination

Demand response (DR) plays a key role for optimum asset utilization and to avoid or delay the need of new infrastructure investment. However, coordinated execution of multiple DRs is desired to maximize the DR benefits. In this paper, we propose a hierarchical DR architecture (HDRA) to control and coordinate the performance of various DR categories such that the operation of every DR category is backed-up by time delayed action of the others. A reliable, cost-effective communication infrastructure based on ZigBee, WiMAX, and fibers is designed to facilitate the HDRA execution. The performance of the proposed HDRA is demonstrated from the power system and communication perspectives in a cosimulation environment applied to a 0.4 kV/400 kVA real distribution network considering electric vehicles as a potential DR resource (DRR). The power simulation is performed employing a real time digital simulator whereas the communication simulation is performed using OMNeT++. Finally, the HDRA performance demonstrated the maximum utilization of available DR potential by facilitating simultaneous execution of multiple DRs and enabling participation of single DRR for multiple grid applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Facile Thermal W–W Bond Homolysis in the N-Heterocyclic Carbene Containing Tungsten Dimer [CpW(CO) 2 (IMe)] 2

The thermal W-W bond homolysis in [CpW(CO)2(IMe)]2 (IMe = 1,3-dimethylimidazol-2-ylidene) was investigated and was found to occur to a large extent compared to other tungsten dimers such as [CpW(CO)3]2. CpW(CO)2(IMe)H was prepared by heating a solution of [IMeH]+[CpW(CO)2(PMe3)]-, and exists in solution as a mixture of interconverting cis and trans isomers. The carbene rotation in CpW(CO)2(IMe)H was explored by DFT calculations, and low enthalpic barriers (< 3.5 kcal mol-1) are predicted. CpW(CO)2(IMe)H has pKaMeCN = 31.5(3) and deprotonation with KH gives K+[CpW(CO)2(IMe)]- (• MeCN). Hydride abstraction from CpW(CO)2(IMe)H with Ph3C+PF6- in the presence of a coordinating ligand L (MeCN or THF) gives [CpW(CO)2(IMe)(L)]+PF6-. Electrochemical measurements on the anion [CpW(CO)2(IMe)]- in MeCN, together with digital simulations, give an E1/2 of -1.54(2) V vs Cp2Fe+/0 for the [CpW(CO)2(IMe)]•/- couple. A thermochemical cycle provides the solution bond dissociation free energy of the W-H bond of CpW(CO)2(IMe)H as 61.3(6) kcal mol-1. In the electrochemical oxidation of [CpW(CO)2(IMe)]-, reversible dimerization of the electrogenerated radical CpW(CO)2(IMe)• occurs, and digital simulation provides kinetic and thermodynamic parameters for the monomer-dimer equilibrium: kdimerization ~ 2.5 x 104 M-1 s-1, khomolysis ~ 0.5 s-1 (i.e., Kdim ~ 5 x 104 M-1). Reduction of [CpW(CO)2(IMe)(MeCN)]+PF6- with cobaltocene gives the dimer [CpW(CO)2(IMe)]2, which in solution exists as a mixture of anti and gauche rotomers. As expected from the electrochemical experiments, the dimer is in equilibrium with detectable amounts of CpW(CO)2(IMe)•. This species was observed by IR spectroscopy, and its presence in solution is also in accordance with the observed reactivity toward 2,6-di-tert-butyl-1,4-benzoquinone, chloroform and dihydrogen. This work was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Biosciences and Geosciences. Pacific Northwest National Laboratory is a multiprogram national laboratory operated for DOE by Battelle. The EPR studies were performed at EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research located at Pacific Northwest National Laboratory.

36 MATERIALS SCIENCE↗

Horizons of modern molecular dynamics simulation in digitalized solid freeform fabrication with advanced materials

Our ability to shape and finish a component by combined methods of fabrication including (but not limited to) subtractive, additive and/or no theoretical mass-loss/addition during the fabrication is now popularly known as solid free form fabrication. Fabrication of a telescope mirror is a typical example where grinding and polishing processes are first applied to shape the mirror and thereafter an optical coating is usually applied to enhance its optical performance. The area of nanomanufacturing cannot grow without a deep knowledge of the fundamentals of materials and consequently, the use of computer simulations is becoming ubiquitous. This article is intended to introduce the most recent advances in the computation benefit specific to the area of solid free form fabrication as these systems are traversing through the journey of digitalisation and Industry-4.0. Specifically, this article demonstrates that the application of the latest materials modelling approaches, based on techniques such as molecular dynamics, are enabling breakthroughs in applied precision manufacturing techniques.

36 MATERIALS SCIENCE↗

Improved reconstruction stability for chemical shift encoded hyperpolarized 13 C magnetic resonance spectroscopic imaging using k‐t spiral acquisitions

Purpose A multiecho, field of view (FOV)‐oversampled k‐t spiral acquisition and direct iterative decomposition of water and fat with echo asymmetry and least‐squares estimation reconstruction is demonstrated to improve the stability of hyperpolarized 13 C magnetic resonance spectroscopic imaging (MRSI) in the presence of signal ambiguities attributed to low‐SNR (signal‐to‐noise‐ratio) species, local uncertainties in metabolite peaks, and echo‐to‐echo signal inconsistencies. Theory k‐t spiral acquisitions redistribute readout points to be more densely spaced radially in k‐space by acquiring an FOV and matrix that are oversampled by η . These more densely spaced spiral turns constitute effective intraspiral echoes and can supplement conventional interspiral echoes to improve spectral separation and reduce spectral cross‐talk to better resolve 13 C‐labeled species for spectroscopic imaging. Methods Digital simulations and imaging phantom experiments were performed for a range of interspiral echo spacings and η using multiecho, k‐t spiral acquisitions. Image spectral cross‐talk artifacts were evaluated both qualitatively and quantitatively as the percent error in measured metabolite ratios. In vivo murine experiments evaluated the feasibility of multiecho, k‐t spiral [1‐ 13 C]pyruvate MRSI to reduce spectral cross‐talk for 3 scenarios of different expected reconstruction stability. Results Digital simulations and imaging phantom experiments both demonstrated reduced or comparable image spectral cross‐talk and percent errors in measured metabolite ratios with increasing η and better choices of echo spacings. In vivo images displayed markedly reduced spectral cross‐talk in lactate images acquired with η = 7 versus η = 1. Conclusion The precision of hyperpolarized 13 C metabolic imaging and quantification in the presence of low‐SNR species, local uncertainties in metabolite resonances, and echo‐to‐echo signal inconsistencies can be improved with the use of FOV‐oversampled k‐t spiral acquisitions.

Macdonald, Erin B.↗

Investigation of aircraft landing in variable wind fields

A digital simulation study is reported of the effects of gusts and wind shear on the approach and landing of aircraft. The gusts and wind shear are primarily those associated with wind fields created by surface wind passing around bluff geometries characteristic of buildings. Also, flight through a simple model of a thunderstorm is investigated. A two-dimensional model of aircraft motion was represented by a set of nonlinear equations which accounted for both spatial and temporal variations of winds. The landings of aircraft with the characteristics of a DC-8 and a DHC-6 were digitally simulated under different wind conditions with fixed and automatic controls. The resulting deviations in touchdown points and the controls that are required to maintain the desired flight path are presented. The presence of large bluff objects, such as buildings in the flight path is shown to have considerable effect on aircraft landings.

Frost, W.↗

Time of Flight Electrochemistry: Diffusion Coefficient Measurements Using Interdigitated Array (IDA) Electrodes

A simple and straightforward method for measuring diffusion coefficients using interdigitated array (IDA) electrodes is reported. The method does not require that the exact electrode area be known but depends only the size of the gap between the IDA electrode pairs. Electroactive molecules produced at the generator electrode of the IDA by a voltage step or scan can diffuse to the collector electrode and the time delay before the current for the reverse electrochemical reaction is detected at the collector is used to calculate the diffusion coefficient. The measurement of the diffusion rate of Ru(NH3)6+2 in aqueous solution has been used as an example measuring diffusion coefficients using this method. Additionally, a digital simulation of the electrochemical response of the IDA electrodes was used to simulate the entire current/voltage/time behavior of the system and verify the experimentally measured diffusion coefficients. This work was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the Department of Energy, Office of Science, Office of Basic Energy Sciences.

Liu, Fei↗

Kinetic Analysis of Competitive Electrocatalytic Pathways: New Insights into Hydrogen Production with Nickel Electrocatalysts

The hydrogen production electrocatalyst Ni(PPh2NPh2)22+ (1) is capable of traversing multiple electrocatalytic pathways. When using dimethylformamidium, DMF(H)+, the mechanism of formation of H2 catalyzed by 1 changes from an ECEC to an EECC mechanism as the potential approaches the Ni(I/0) couple. Two recent electrochemical methods, current-potential analysis and foot-of-the-wave analysis (FOWA), were performed on 1 to measure the detailed chemical kinetics of the competing ECEC and EECC pathways. A sensitivity analysis was performed on the electrochemical methods using digital simulations to gain a better understanding of their strengths and limitations. Notably, chemical rate constants were significantly underestimated when not accounting for electron transfer kinetics, even when electron transfer was fast enough to afford a reversible non-catalytic wave. The EECC pathway of 1 was found to be faster than the ECEC pathway under all conditions studied. Using buffered DMF: DMF(H)+ mixtures led to an increase in the catalytic rate constant (kobs) of the EECC pathway, but kobs for the ECEC pathway did not change when using buffered acid. Further kinetic analysis of the ECEC path revealed that added base increases the rate of isomerization of the exo-protonated Ni(0) isomers to the catalytically active endo-isomers, but decreases the net rate of protonation of Ni(I). FOWA on 1 did not provide accurate rate constants due to incomplete reduction of the exo-protonated Ni(I) intermediate at the foot of the wave, but FOWA could be used to estimate the reduction potential of this previously undetected intermediate. This research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences. Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy.

electrocatalysis↗

Digital quantum simulation of cavity quantum electrodynamics: insights from superconducting and trapped ion quantum testbeds

We explore the potential for hybrid development of quantum hardware where currently available quantum computers simulate open cavity quantum electrodynamical (CQED) systems for applications in optical quantum communication, simulation and computing. Our simulations make use of a recent quantum algorithm that maps the dynamics of a singly excited open Tavis–Cummings model containing N atoms coupled to a lossy cavity. We report the results of executing this algorithm on two noisy intermediate-scale quantum computers: a superconducting processor and a trapped ion processor, to simulate the population dynamics of an open CQED system featuring N = 3 atoms. By applying technology-specific transpilation and error mitigation techniques, we minimize the impact of gate errors, noise, and decoherence in each hardware platform, obtaining results which agree closely with the exact solution of the system. These results can be used as a recipe for efficient and platform-specific quantum simulation of cavity–emitter systems on contemporary and future quantum computers.

cavity QED↗

Quantum Computer-Aided Design: Digital Quantum Simulation of Quantum Processors

With the increasing size of quantum processors, submodules that constitute the processor hardware will become too large to accurately simulate on a classical computer. Therefore, one would soon have to fabricate and test each new design primitive and parameter choice in time-consuming coordination between design, fabrication, and experimental validation. Here we show how one can design and test the performance of next-generation quantum hardware—by using existing quantum computers. Focusing on superconducting transmon processors as a prominent hardware platform, we compute the static and dynamic properties of individual and coupled transmons. We show how the energy spectra of transmons can be obtained by variational hybrid quantum-classical algorithms that are well suited for near-term noisy quantum computers. In addition, single- and two-qubit gate simulations are demonstrated via Suzuki-Trotter decomposition. Our methods pave a promising way towards designing candidate quantum processors when the demands of calculating submodule properties exceed the capabilities of classical computing resources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗