Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient problem”

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 127 records · Page 7

Convex Optimization for Nonequilibrium Steady States on a Hybrid Quantum Processor

Finding the transient and steady state properties of open quantum systems is a central problem in various fields of quantum technologies. Here, in this work, we present a quantum-assisted algorithm to determine the steady states of open system dynamics. By reformulating the problem of finding the fixed point of Lindblad dynamics as a feasibility semidefinite program, we bypass several well-known issues with variational quantum approaches to solving for steady states. We demonstrate that our hybrid approach allows us to estimate the steady states of higher dimensional open quantum systems and discuss how our method can find multiple steady states for systems with symmetries.

97 MATHEMATICS AND COMPUTING↗

Dating submarine landslides using the transient response of gas hydrate stability

Abstract Submarine landslides are prevalent on the modern-day seafloor, yet an elusive problem is constraining the timing of past slope failure. We present a novel age-dating technique based on perturbations to underlying gas hydrate stability caused by slide-impacted seafloor changes. Using three-dimensional (3-D) seismic data, we mapped an irregular bottom simulating reflection (BSR) underneath a submarine landslide in the Orca Basin, Gulf of Mexico. The irregular BSR mimics the pre-slide seafloor geometry rather than the modern bathymetry. Therefore, we suggest that the gas hydrate stability zone (GHSZ) is still adjusting to the post-slide sediment temperature. We applied transient conductive heat-flow modeling to constrain the response of the GHSZ to the slope failure, which yielded a most likely age of ca. 8 ka, demonstrating that gas hydrate can respond to landslides even on multimillennial time scales. We further provide a generalized analytical solution that can be used to remotely date submarine slides in the absence of traditional dating techniques.

Geology↗

Multichannel Deconvolution of Vibrational Shock Signals: An Inverse Filtering Approach

When transporting critical systems of national security interest, out-of-the ordinary, impulsive events that can potentially be undetected and affect overall system performance are of great concern. Impulsive events that can occur are essentially pulse-like, transient signals of short duration that evolve from various phenomena. Here the event can be created by either the dropping of a test object, the system, subjecting it to a compact high-energy blow or being struck unintentionally during transit resulting in potential damage. The intensity and location of the strike can cause an inoperability condition that is unacceptable in a national security environment. Therefore, it is essential to detect, classify and localize damage of any test object subjected to an impulsive-event. This effort was targeted to evaluate the vibrational response of test objects that are subjected to “transport” shocks and roadway vibrations during shipping and handling. Any potential damage that could be inflicted during transportation must not only be detected, but also be evaluated to determine the operational readiness of a test object before and after transport. This event is a critical task that must be addressed as part of the Lawrence Livermore National Laboratory (LLNL) national security mission. The estimation of excitation signals from noisy data is termed the deconvolution problem in the signal processing literature. The deconvolution problem is based on recovering the input excitation signal from a system characterized by its impulse response sequence. Using this model of the system, an “inverse” representation or filter is developed to remove the system from the measured data and recover the input. Deconvolution techniques have existed for a long-time; however, transient deconvolution presents a few uncommon problems, since the signal has a finite-length time duration resulting in a limited amount of data containing information about the excitation process. The transient is wideband in the frequency domain relative to any measurement sensor implying that the smaller bandwidth sensor system “filters” the excitation eliminating some of its essential information for recovery. This fact, coupled with the filtering effect of the test object itself makes this excitation recovery (deconvolution) problem a challenge for signal processing.

42 ENGINEERING↗

MSR Transient Simulation and MSRE Transient Benchmark with SAM and SPECTRA

In recent years, there has been renewed interest in molten salt reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. In response to such interest, many methods and codes have been developed to capture the unique features of MSRs. Among them, SPECTRA and SAM are two system analysis codes that have been enhanced to include MSR-specific modeling capabilities, including delayed neutron precursor drift and modified point kinetics equations.This paper discusses the efforts taken to verify and validate these features. A standard MSR system test problem was developed to verify and demonstrate the capability of SPECTRA and SAM on the MSR transient simulation. Sixteen transients were simulated. Further, the results obtained from SPECTRA and SAM show good agreement. The Molten Salt Reactor Experiment transient experiments were reviewed and selected to validate the SPECTRA and SAM codes. The experiments included pump startup and coastdown tests at zero power, reactivity insertion tests at different power levels, frequency tests, and a natural convection test. The simulation results from SPECTRA and SAM show good agreement with the experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Understanding Transients: Needs from the Laboratory Astrophysics Community

A growing number of astrophysical transients are pushing astronomers to develop increasingly complex computational tools to model both radiation hydrodynamics and electron transport. Here we review the physics needs and simulation uncertainties associated with modeling these transients. Although this review will focus on theory and simulation aspects of this problem, we also review some experiments designed to study this physics.

79 ASTRONOMY AND ASTROPHYSICS↗

The Transient Multi-Level method for Monte Carlo reactor statics calculations

The Transient Multi-Level (TML) method is applied to a time-dependent Monte Carlo transport solver to offload some of the computational burden of the expensive Monte Carlo solve to lower-order Coarse Mesh Finite Difference (CMFD) and Exact Point Kinetics Equations (EPKE) solvers via factorization of the neutron flux at the transport and CMFD levels using the Predictor Corrector Quasi-Static Method (PCQM). The Monte Carlo transient is solved by a modified fission source iteration scheme that introduces a single transient source bank. The method is implemented in the production-level Monte Carlo code, Shift, and verified with prescribed reactivity ramps from the two-dimensional version of the C5G7-TD reactor benchmark. The results show that, as compared to other quasi-static methods, the TML reduces the stochastic noise inherent to the transient Monte Carlo solver by factors of ~2 to 6 for various norm comparisons of the reactor power amplitude. Finally, the TML additionally reduces the number of Monte Carlo evaluations needed to simulate the transient, leading to roughly an order of magnitude improvement in CPU time relative to the standard PCQM for the problems tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Transient Stability Preventive Control via Tuning the Parameters of Virtual Synchronous Generators

This paper presents an optimal preventive control (OPC) method to improve the power system transient stability via tuning the transient parameter of virtual synchronous generators. The novelty of this work is that we formulate the preventive control as an optimization problem so that inverter parameters can be adjusted at the pre-contingency stage. A reinforcement learning (RL)-driven method is proposed to solve the OPC problem with the fault energy-based reward function. An ANDES-based RL environment is also developed. Versatile functions included in the proposed environment have been presented in this paper. The proposed OPC formulation, the RLdriven method, and the fault energy-based reward function are verified on several standard test systems.

Inverter-based resources, virtual synchronous gene↗

Serrated Flow in Alloy Systems

This chapter presents the problem of the complexity of plastic flow in alloys, which is manifested by serrated deformation curves and transient plastic strain localizations. This phenomenon, which uncovers an inherently collective nature of the dynamics of crystal defects, is well known for conventional alloys. Surprisingly, despite the particular microstructure of HEAs, which strongly impacts the microscopic dynamics of defects, not only the HEAs are prone to the macroscopically jerky flow, but the serration patterns display many similar features, thus bearing evidence to similar laws for the collective dynamics. To provide a comprehensive approach to this problem, the chapter is subdivided into two parts. The first part presents the state-of-the art of this research under the angle of distinct regimes of collective dislocation dynamics found in conventional alloys, e.g., self-organized criticality, deterministic chaos, or synchronization. The second part summarizes recent investigations into serrated flows in HEAs, uncovering common and specific types of behaviors. Both parts conclude with the formulation of diverse open questions concerning these strongly coupled fields of research.

Lebyodkin, Mikhail↗

EMT Model Validation of a 2 MVA PV Inverter Using Transient and Frequency-Domain Hardware Testing: Preprint

This paper presents results and new insights gained from a hardware test campaign on a 2 MVA PV inverter for validating its vendor-supplied EMT model. The test campaign was conducted using a 7 MVA grid simulator and a 2 MW PV emulator. It considered both time-domain transient tests and frequency-domain impedance scan tests. The paper highlights the inadequacy of transient tests in capturing all critical resonance modes of the inverter, and the effectiveness of the frequency scan testing in addressing this problem. The paper shows the frequency scan testing as an effective tool for EMT model validation of IBR units which can highlight inaccuracies in the EMT models that are easy to overlook when model validation is performed using only the time-domain transient tests such as voltage ride-through and phase jump tests.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Unitarity of quantum tunneling decay for an analytical exact non-Hermitian resonant-state approach

Highlights: • Foundations of quantum mechanics • Quantum decay • Non-Hermitian Hamiltonian • Unitarity • Resonant state By using an analytical exact non-Hermitian formalism for quantum tunneling decay that involves the expansion of the decaying wave function as a linear combination of resonant states and transient functions associated with the complex poles of the outgoing Green’s function to the problem, it is shown that the integrated decaying probability density in the whole space satisfies unitarity at each value of time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗

Predicting distribution of malaria vector larval habitats in Ethiopia by integrating distributed hydrologic modeling with remotely sensed data

Abstract Larval source management has gained renewed interest as a malaria control strategy in Africa but the widespread and transient nature of larval breeding sites poses a challenge to its implementation. To address this problem, we propose combining an integrated high resolution (50 m) distributed hydrological model and remotely sensed data to simulate potential malaria vector aquatic habitats. The novelty of our approach lies in its consideration of irrigation practices and its ability to resolve complex ponding processes that contribute to potential larval habitats. The simulation was performed for the year of 2018 using ParFlow-Common Land Model (CLM) in a sugarcane plantation in the Oromia region, Ethiopia to examine the effects of rainfall and irrigation. The model was calibrated using field observations of larval habitats to successfully predict ponding at all surveyed locations from the validation dataset. Results show that without irrigation, at least half of the area inside the farms had a 40% probability of potential larval habitat occurrence. With irrigation, the probability increased to 56%. Irrigation dampened the seasonality of the potential larval habitats such that the peak larval habitat occurrence window during the rainy season was extended into the dry season. Furthermore, the stability of the habitats was prolonged, with a significant shift from semi-permanent to permanent habitats. Our study provides a hydrological perspective on the impact of environmental modification on malaria vector ecology, which can potentially inform malaria control strategies through better water management.

59 BASIC BIOLOGICAL SCIENCES↗

Flux-based three-dimensional electrodynamic modeling approach to superconducting circuits and materials

Modeling the behavior of superconducting electronic circuits containing Josephson junctions is crucial for the design of superconducting information processors and devices. In this paper, we introduce DEC-QED, a computational approach for modeling the electrodynamics of superconducting electronic circuits containing Josephson junctions in arbitrary three-dimensional electromagnetic environments. DEC-QED captures the nonlinear response and induced currents in BCS superconductors and accurately captures phenomena such as the Meissner effect, flux quantization, and Josephson effects. Using a spatial coarse-graining formulation based on discrete exterior calculus (DEC), DEC-QED can accurately simulate transient and long-time dynamics in superconductors. Finally, the expression of the entire electrodynamic problem in terms of the gauge-invariant flux field and charges makes the resulting classical field theory suitable for second quantization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Developing a PSCAD Model of the Reduced 240-Bus WECC Test System: Preprint

Fast-acting controls of inverter-based resources (IBR) introduce unprecedented dynamics to bulk power systems and island power grids, especially the high-frequency dynamics ranging from $>5$ Hz to hundreds of Hz. These oscillatory stability problems cannot be captured by the phasor domain simulations, e.g., in PSS/E, Powerworld and TSAT, due to insufficient modeling. Instead, they should be studied by electromagnetic transient (EMT) simulations by using proper EMT dynamic models. This paper develops and validates the PSCAD model of the reduced 240-bus Western Electricity Coordinating Council (WECC) test system. The PSCAD model will be released to the public, which can be used by both academia and power industry practitioners for studying the fast dynamic problems without any Critical Energy/Electric Infrastructure Information (CEII).

240-bus WECC test system↗

Development of an Improved RELAP5-3D Model for the High Temperature Test Facility

High-temperature gas-cooled reactors (HTGRs) are rapidly approaching deployment. Confidence in transient analysis of these systems requires modeling and simulation tools that have been validated against data relevant to HTGR conditions. The High Temperature Test Facility (HTTF) is an integral effects thermal hydraulics test facility for prismatic HTGRs. In spring and summer of 2019, HTTF was used for a series of experiments that now serve as the basis for the Organization of Economic Cooperation and Development / Nuclear Energy Agency Thermal Hydraulic Code Validation Benchmark for High Temperature Gas-Cooled Reactors using HTTF Data (HTGR T/H Benchmark). Previous analyses as part of the HTGR T/H benchmark used a RELAP5-3D model developed at Idaho National Laboratory (INL) and demonstrated an ability to reproduce trends in the measured data but difficulties reproducing experimental values within their uncertainty. These difficulties were largely attributed to assumptions made during the development of the initial RELAP5-3D model, which predated the HTTF experiments. A significant cause of difficulty reproducing the measured temperatures may be the radial nodalization of the previous RELAP5-3D model. In this paper, we present a new RELAP5-3D model of HTTF with finer radial nodalization built to assess the impact of radial heat transfer. We describe the new model and compare it against the old one at full-power steady state and for the pressurized conduction cooldown (PCC) transient. These analyses are based on the code-to-code comparison exercise for the PCC problem of the HTGR T/H benchmark. We compare maximum block temperature as the primary figure of merit and include discussion on intracore natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Transient Studies of Nonequilibrium Transport in Two-Dimensional Semiconductors

The overarching objective of research under DOE support has been to investigate emergent carrier phenomena in different two-dimensional (2D) systems, under conditions where they are driven far from equilibrium. The description of nonequilibrium transport has long represented one of the most enduring challenges in condensed-matter physics, in large part due to the need to invoke the role of complex many-body phenomena (electron-phonon and electron-electron coupling) in this regime. In this project, we have explored aspects of this problem in a variety of different two-dimensional systems. These include the modulation-doped two-dimensional electrons gas, or 2DEG, in GaAs/AlGaAs heterostructures; the electron and hole systems in graphene and its related hetero-structures, and; the strongly-confined electron layers in transition-metal dichalcogenide (MoS 2 and WS 2 ) and trichalcogenide (TiS 3 ) materials. Significant understanding has been achieved in a number of different areas, as we elaborate upon in subsequent sections of this proposal: Nanoscale spectroscopy of defect-related barriers in MoS2 Negative differential conductance due to inter-valley scattering in WS 2 under high fields Incipient localization due to scattering from frozen-phonon disorder in quantum point contacts High-field transient transport in graphene/hexagonal BN heterostructures Universal voltage scaling of quantum corrections in graphene under nonequilibrium A defining feature of much of this research has been an approach to problem solving in which novel experimental techniques (such as transient electrical pulsing and scanning-probe microscopy) are combined with the development (in collaboration with colleagues from condensed-matter theory) of sophisticated models, a strategy that we propose to continue in our future research. In recent years, we have developed an increasing capacity for handling atomically-thin materials, and their van der Waals systems, a capacity that is reflected in the balance of topics listed above.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Impact of Integration Scheme on Performance of Anisotropic Plasticity Models

Given the prevalent role of metals in a variety of industries, schemes to integrate corresponding constitutive models in finite element applications have long been studied. A number of formulations have been developed to accomplish this task; each with their own advantages and costs. Often the focus has been on ensuring the accuracy and numerical stability of these algorithms to enable robust integration. While important, emphasis on these performance metrics may often come at the cost of computational expense potentially neglecting the needs of individual problems. In the current work, the performance of two of the most common integration methods for anisotropic plasticity -- the convex cutting plane (CCP) and closest point projection (CPP) -- across a variety of metrics is assessed; including accuracy and cost. A variety of problems are considered ranging from single elements to large representative simulations including both implicit quasistatic and explicit transient dynamic type responses. The relative performance of each scheme in the different instances is presented with an eye towards guidance on when the different algorithms may be beneficial.

36 MATERIALS SCIENCE↗