Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear models”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Nonadiabatic transitions during a passage near a critical point

The passage through a critical point of a many-body quantum system leads to abundant nonadiabatic excitations. Here, we explore a regime, in which the critical point is not crossed although the system is passing slowly very close to it. Here, we show that the leading exponent for the excitation probability can then be obtained by standard arguments of the Dykhne formula, but the exponential prefactor is no longer simple and behaves as a power law on the characteristic transition rate. We derive this prefactor for the nonlinear Landau–Zener model by adjusting Dykhne’s approach. Then, we introduce an exactly solvable model of the transition near a critical point in the Stark ladder. We derive the number of excitations for it without approximations and find qualitatively similar results for the excitation scaling.

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

Landsat 8 monitoring of multi-depth suspended sediment concentrations in Lake Erie’s Maumee River using machine learning

Satellite remote sensing has been widely used to map suspended sediment concentration (SSC) in waterbodies. However, due to the complexity of sediment-water interactions, it has been difficult to derive linear and non-linear regression equations to reliably predict SSC, especially when trying to estimate depth of integrated sediment. Herein, this study uses Landsat 8 OLI (Operational Land Imager) sensor to map SSC within the Maumee River in Ohio, USA, at multiple depth intervals (15, 61, 91, and 182 cm). Simple linear least squares regression (LLSR), and three common machine learning models: random forest (RF), support vector regression (SVR), and model averaged neural network (MANN) were used to estimate SSC at the depth intervals. All machine learning models significantly outperformed LLSR while RF performed the best. In both RF and MANN, R2 (coefficient of determination) increases with depth with a maximum R2 of 0.89 and 0.83, respectively, at a depth of 0–182 cm. The results show that machine learning models can implement nonlinear relationships that produce better predictions than traditional linear regression methods in estimating depth integrated SSC, especially when samples are limited.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Neutrino winds on the sky

Abstract We develop a first-principles formalism to compute the distortion to the relic neutrino density field caused by the peculiar motions of large-scale structures. This distortion slows halos down due to dynamical friction, causes a local anisotropy in the neutrino-CDM cross-correlation, and reduces the global cross-correlation between neutrinos and CDM. The local anisotropy in the neutrino-CDM cross-spectrum is imprinted in the three point cross-correlations of matter and galaxies, or the bispectrum in Fourier space, producing a signal peaking at squeezed triangle configurations. This bispectrum signature of neutrino masses is not limited by cosmic variance or potential inaccuracies in the modeling of complicated nonlinear and galaxy formation physics, and it is not degenerate with the optical depth to reionization. We show that future surveys have the potential to detect the distortion bispectrum.

Astronomy & Astrophysics↗
DOE OSTI2021

Optimal Energy Scheduling and Sensitivity Analysis for Integrated Power-Water-Heat Systems

The conventionally independent power, water, and heating networks are becoming more tightly connected, which motivates their joint optimal energy scheduling to improve the overall efficiency of an integrated energy system. However, such a joint optimization is known as a challenging problem with complex network constraints and couplings of electric, hydraulic, and thermal models that are nonlinear and nonconvex. We formulate an optimal power-water-heat flow (OPWHF) problem and develop a computationally efficient heuristic to solve it. The proposed heuristic decomposes OPWHF into subproblems, which are iteratively solved via convex relaxation and convex-concave procedure. Simulation results validate that the proposed framework can improve operational flexibility and social welfare of the integrated system, wherein the water and heating networks respond as virtual energy storage to time-varying energy prices and solar photovoltaic generation. Moreover, we perform sensitivity analysis to compare two modes of heating network control: by flow rate and by temperature. Our results reveal that the latter is more effective for heating networks with a wider space of pipeline parameters.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Tractable Learning in Underexcited Power Grids

Estimating the structure of physical flow networks, such as power grids, is critical to secure delivery of energy. This article discusses statistical structure estimation in power grids in the “underexcited” regime, where a subset of internal nodes has zero injection fluctuations. Prior estimation algorithms based on nodal voltages fail for such grids as the voltage covariance matrix is not invertible. Here, we propose a novel topology learning algorithm for learning underexcited general networks. Our algorithm uses physics-informed conservation laws to first identify the zero-injection buses and their neighbors, and then estimates the remaining edges in the grid. We prove the asymptotic correctness of our algorithm for grids with nonadjacent internal zero-injection nodes. More important, we theoretically analyze our algorithm’s efficacy under noisy measurements, and determine bounds on maximum noise under which asymptotically correct recovery is guaranteed. Our approach is validated through simulations with voltage samples generated on test distribution grids with real injection data and nonlinear power flow models.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Stacked Low-Inertia Converter or Solid-State Transformer: Modeling and Model Predictive Priority-Shifting Control for Voltage Balance

This paper presents control challenges of stacked low-inertia converter (SLIC) or cascaded reduced dc-link solid-state transformer (SST) and proposes a novel model predictive priority-shifting (MPPS) control with implicit modulator and a discrete-time large-signal model for voltage balancing and dc-link regulation. Low-inertia converters, featuring small electrolytic capacitor-less dc links, dramatically reduce cost, size, and weight compared to conventional solutions. However, without a large dc-link buffer, the input and output are tightly coupled, leading to significant control challenges. The control becomes even more challenging with these converters stacked input-series output-parallel (ISOP) for medium-voltage (MV) grid, which causes coupling between the modules besides the coupling within each module. This paper analyzes the multi-objective, multi-degree of freedom control problem, using the modular soft-switching solid-state transformer (M-S4T) as an example of the SLIC. First, distribution of control efforts under controller saturation is critical because multiple control objectives can be conflicting, especially when the module voltages are unbalanced and are being restored. The MPPS can shift the priorities to address this issue. Second, due to the low inertia and high dc-link ripple, classic space vector pulse-width modulation (SVPWM), average model with small-ripple assumption, and control design based on small-signal model cannot accurately modulate, model, and control the nonlinear reduced dc link. Therefore, a discrete-time large-signal model of the M-S4T is established to derive the predictive control in the MPPS. The MPPS and the PI control are compared in MV simulations to show the issue of applying the PI to the SLIC and the effectiveness of the MPPS for voltage balancing and dc-link regulation in a deadbeat manner. Finally, the proposed control is tested on a 5 kV ISOP SiC SST prototype to verify priority shifting to address controller saturation issue and fast and robust voltage balancing.

14 SOLAR ENERGY↗
DOE OSTI2021

An Adaptive-Importance-Sampling-Enhanced Bayesian Approach for Topology Estimation in an Unbalanced Power Distribution System

The reliable operation of a power distribution system relies on a good prior knowledge of its topology and its system state. Although crucial, due to the lack of direct monitoring devices on the switch statuses, the topology information is often unavailable or outdated for the distribution system operators for real-time applications. Apart from the limited observability of the power distribution system, other challenges are the nonlinearity of the model, the complicated, unbalanced structure of the distribution system, and the scale of the system. To overcome the above challenges, we, in this paper, propose a Bayesian-inference framework that allows us to simultaneously estimate the topology and the state of a three-phase, unbalanced power distribution system. Specifically, by using the very limited number of measurements available that are associated with the forecast load data, we efficiently recover the full Bayesian posterior distributions of the system topology under both normal and outage operation conditions. This is performed through an adaptive importance sampling procedure that greatly alleviates the computational burden of the traditional Monte-Carlo (MC)-sampling-based approach while maintaining a good estimation accuracy. The simulations conducted on the IEEE 123-bus test system and an unbalanced 1282-bus system reveal the excellent performances of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Investigating the Influence of Ni, ZrO 2 , and Y 2 O 3 from SOFC Anodes on Siloxane Deposition

Siloxanes, as a type of impurity in biogas, can poison the Ni-YSZ anode of SOFCs. However, the influence of individual components of the anode, such as Ni, ZrO 2 , and Y 2 O 3 , on the siloxane deposition process has not been investigated extensively. In this study, Ni, ZrO 2 , and Y 2 O 3 pellets were exposed to H 2 + N 2 + H 2 O + D4 (octamethylcyclotetrasiloxane, 2.5 ppmv) and H 2 + N 2 + D4 (2.5 ppmv siloxane) gas mixtures at 750 °C to investigate their affinity and tolerance for siloxane degradation. Surface morphology analysis and electrochemical analysis including electrochemical impedance spectroscopy (EIS), related distribution of relaxation times (DRT) analysis and equivalent circuit modeling with complex nonlinear least square (CNLS) fitting were conducted. Here, a microstructure parameter—tortuosity factor to porosity ratio $\tau /\varepsilon $ calculated by diffusion polarization resistance was utilized for siloxane deposition evaluation. After comparing pellets surface morphology changes before and after experiments and $\tau /\varepsilon $ change following the contamination, Ni is considered as a major factor in siloxane deposition reactions in Ni-YSZ anode.

25 ENERGY STORAGE↗
DOE OSTI2024

Thermodynamics-informed latent space dynamics identification

This software showcases a latent space dynamics identification method, namely tLaSDI, that embeds the first and second principles of thermodynamics. The latent variables are learned through an autoencoder as a nonlinear dimension reduction model. The latent dynamics are constructed by a neural network-based model that precisely preserves certain structures for the thermodynamic laws through the GENERIC formalism. An abstract error estimate is established, which provides a new loss formulation involving the Jacobian computation of autoencoder. The autoencoder and the latent dynamics are simultaneously trained to minimize the new loss. Computational examples demonstrate the effectiveness of tLaSDI, which exhibits robust generalization ability, even in extrapolation. In addition, an intriguing correlation is empirically observed between a quantity from tLaSDI in the latent space and the behaviors of the full-state solution.

Cheung, Siu Wun↗
DOE OSTI2023

Physics of the thermal quench in MST tokamak plasmas (Final Technical Report)

This project at the University of Wisconsin-Madison targeted measurements in tokamak plasmas in the Madison Symmetric Torus (MST) to help elucidate the physics of the tokamak disruption. We intended to focus on the thermal quench, during which energy is lost from the plasma in a very short time. Understanding and ultimately mitigating the thermal quench are among the highest priority research thrusts for ITER, a large international fusion experiment presently under construction. The measurements in MST exploited an advanced set of diagnostics for internal measurements with high time resolution, targeting linear and nonlinear magnetohydrodynamic (MHD) activity as well as the electron temperature and density profiles. We also exploited other features of the MST device, including very flexible power supplies for the toroidal magnetic field and toroidal plasma current that are capable of producing almost arbitrary waveforms. While the work funded by this project was purely experimental, we continued and expanded an existing collaboration with UW-Madison experts in nonlinear MHD computational modeling of tokamak plasmas. That modeling incorporated realistic, experimentally measured parameters. The results of the modeling were anticipated to deepen considerably our understanding of the experimental measurements, and in turn, the experimental measurements were anticipated to help validate the modeling results. The science of transients such as disruptions is one of the highest priority research areas in the U.S. and international fusion programs. The work carried out in this project was targeted to help advance the science of tokamak disruptions which will in turn contribute to the success of ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Utilizing Reinforcement Learning to Continuously Improve a Primitive-Based Motion Planner

We report in this paper describes how the performance of motion primitive-based planning algorithms can be improved using reinforcement learning. Specifically, we describe and evaluate a framework that autonomously improves the performance of a primitive-based motion planner. The improvement process consists of three phases: exploration, extraction, and reward updates. This process can be iterated continuously to provide successive improvement. The exploration step generates new trajectories, and the extraction step identifies new primitives from these trajectories. These primitives are then used to update rewards for continued exploration. This framework required novel shaping rewards, development of a primitive extraction algorithm, and modification of the Hybrid A* algorithm. The framework is tested on a navigation task using a nonlinear F-16 model. The framework autonomously added 91 motion primitives to the primitive library and reduced average path cost by 21.6 seconds, or 35.75% of the original cost. The learned primitives are applied to an obstacle field navigation task, which was not used in training, and reduced path cost by 16.3 seconds, or 24.1%. Additionally, two heuristics for the modified Hybrid A* algorithm are designed to improve effective branching factor.

42 ENGINEERING↗
DOE OSTI2022

Deep Generative Modeling of Periodic Variable Stars Using Physical Parameters

The ability to generate physically plausible ensembles of variable sources is critical to the optimization of time domain survey cadences and the training of classification models on data sets with few to no labels. Traditional data augmentation techniques expand training sets by reenvisioning observed exemplars, seeking to simulate observations of specific training sources under different (exogenous) conditions. Unlike fully theory-driven models, these approaches do not typically allow principled interpolation nor extrapolation. Moreover, the principal drawback of theory-driven models lies in the prohibitive computational cost of simulating source observables from ab initio parameters. In this work, we propose a computationally tractable machine learning approach to generate realistic light curves of periodic variables capable of integrating physical parameters and variability classes as inputs. Our deep generative model, inspired by the transparent latent space generative adversarial networks, uses a variational autoencoder (VAE) architecture with temporal convolutional network layers, trained using the OGLE-III optical light curves and physical characteristics (e.g., effective temperature and absolute magnitude) from Gaia DR2. A test using the temperature–shape relationship of RR Lyrae demonstrates the efficacy of our generative “physics-enhanced latent space VAE” (PELS-VAE) model. Such deep generative models, serving as nonlinear nonparametric emulators, present a novel tool for astronomers to create synthetic time series over arbitrary cadences.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment

We present a detection of 21 cm emission from large-scale structure (LSS) between redshift 0.78 and 1.43 made with the Canadian Hydrogen Intensity Mapping Experiment. Radio observations acquired over 102 nights are used to construct maps that are foreground filtered and stacked on the angular and spectral locations of luminous red galaxies (LRGs), emission-line galaxies (ELGs), and quasars (QSOs) from the eBOSS clustering catalogs. We find decisive evidence for a detection when stacking on all three tracers of LSS, with the logarithm of the Bayes factor equal to 18.9 (LRG), 10.8 (ELG), and 56.3 (QSO). An alternative frequentist interpretation, based on the likelihood ratio test, yields a detection significance of 7.1σ (LRG), 5.7σ (ELG), and 11.1σ (QSO). These are the first 21 cm intensity mapping measurements made with an interferometer. We constrain the effective clustering amplitude of neutral hydrogen (H I ), defined as ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}\equiv {10}^{3}\,{{\rm{\Omega }}}_{{\rm{H}}\,{\rm\small{I}}}\left({b}_{{\rm{H}}\,{\rm\small{I}}}+\langle \,f{\mu }^{2}\rangle \right)$, where Ω H I is the cosmic abundance of H I , b H I is the linear bias of H I , and $\langle$fμ 2 $\rangle$ = 0.552 encodes the effect of redshift-space distortions at linear order. We find ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={1.51}_{-0.97}^{+3.60}$ for LRGs (z = 0.84), ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={6.76}_{-3.79}^{+9.04}$ for ELGs (z = 0.96), and ${{ \mathcal A }}_{{\rm{H}}\,{\rm\small{I}}}={1.68}_{-0.67}^{+1.10}$ for QSOs (z = 1.20), with constraints limited by modeling uncertainties at nonlinear scales. We are also sensitive to bias in the spectroscopic redshifts of each tracer, and we find a nonzero bias Δ v = - 66 ± 20 km s -1 for the QSOs. We split the QSO catalog into three redshift bins and have a decisive detection in each, with the upper bin at z = 1.30 producing the highest-redshift 21 cm intensity mapping measurement thus far.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Vertical dependence of horizontal variation of cloud microphysics: observations from the ACE-ENA field campaign and implications for warm-rain simulation in climate models

Abstract. In the current global climate models (GCMs), the nonlinearity effect of subgrid cloud variations on the parameterization of warm-rain process, e.g., the autoconversion rate, is often treated by multiplying the resolved-scale warm-rain process rates by a so-called enhancement factor (EF). In this study, we investigate the subgrid-scale horizontal variations and covariation of cloud water content (qc) and cloud droplet number concentration (Nc) in marine boundary layer (MBL) clouds based on the in situ measurements from a recent field campaign and study the implications for the autoconversion rate EF in GCMs. Based on a few carefully selected cases from the field campaign, we found that in contrast to the enhancing effect of qc and Nc variations that tends to make EF > 1, the strong positive correlation between qc and Nc results in a suppressing effect that tends to make EF < 1. This effect is especially strong at cloud top, where the qc and Nc correlation can be as high as 0.95. We also found that the physically complete EF that accounts for the covariation of qc and Nc is significantly smaller than its counterpart that accounts only for the subgrid variation of qc, especially at cloud top. Although this study is based on limited cases, it suggests that the subgrid variations of Nc and its correlation with qc both need to be considered for an accurate simulation of the autoconversion process in GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

DROP DURABILITY ASSESSMENT OF ELECTRONIC ASSEMBLIES UNDER OFF-AXIS LOADING WITH SKEWED FIXTURES

This thesis studies drop durability of electronic assemblies when the acceleration vector is oriented at 45° to the out-of-plane direction of the circuit card. The off-axis drop tests are accomplished with a skewed fixture and are conducted as a proxy for multiaxial drop testing. Advanced shock testing and vibration test methods have been developed over the last few decades to better represent real-world field environments during ground-based laboratory testing. However, many of these test methods require expensive and specialized equipment not available in most laboratories. An alternative approach for approximating simultaneous loading along multiple axes on conventional equipment utilizes skewed fixtures which have seen use in off-axis random vibration and drop impact testing. These methods generally rely on the conversion of a uniaxial input load from the test equipment (using a uniaxial drop tower or shaker) into a multiaxial load when resolved in the reference frame of the test article (mounted on a skewed fixture). Skewed fixture design is presented and recommendations for conducting skewed angle drop testing are introduced based on local measurements along the skewed face of the fixture to accurately monitor the impact event. Characterization tests were performed with a skewed fixture, at simultaneous acceleration loads from 500 to 3,000 g in two (in-plane and out-of-plane) directions, while meeting standard time domain tolerances. Upon experimental characterization, drop shock durability tests were conducted on a printed circuit assembly (PCA). Mean drops-to-failure were measured and quantified with Weibull statistics. Dominant solder joint failure modes were identified via failure analysis. Prior work on inclined angle impact testing is limited, and the majority of solder joint interconnect level fatigue studies are conducted considering perpendicular loading normal the circuit card. Low-cycle fatigue curves are generated based on plastic strain and plastic work density within the solder joint. A multiscale nonlinear finite element model is used to relate board-level flexure to solder joint interconnect level plastic strain. A high strain rate solder constitutive model allows for accurate modeling of solder plasticity resulting from high-impact drop shock. Fatigue parameters are computed from the Coffin-Manson relation and Palmgren-Miner damage accumulation. This work serves to apply established low-cycle fatigue methods for conventional drop shock loading (impact normal to circuit card) to non-perpendicular loading with a skewed fixture.

Hower, Jonathan [Kansas City National Security Cam↗
DOE OSTI2025

Explicit Dynamic Impact Analysis of the Building 3525 Packages

The safe transport and handling of radioactive material containers is critical to ensure the protection of personnel, facilities, and the environment. Accidental drops during handling pose a significant hazard, necessitating robust design and analysis to mitigate potential consequences. This study presents an explicit dynamic impact analysis of two container designs, the Baby Sugarman and the Core Conduction Cooldown Test Facility (CCCTF) Tall Boy, used in Building 3525 at Oak Ridge National Laboratory. The analysis employed nonlinear finite element modeling in LS-DYNA to simulate free-fall impacts from a height of 16.5 ft across eight critical orientations. The study evaluated deformation, strain failure, and shielding integrity under these extreme conditions. Results demonstrated that both designs maintain structural integrity, with no breaches or loss of shielding under specified administrative controls. Notably, a maximum shielding loss of 10 lbs was observed during a top-corner drop, highlighting the need to limit drop heights to less than 1 ft for certain orientations. These findings underscore the effectiveness of engineering controls and design measures in preventing containment breaches and ensuring compliance with safety requirements. This work provides a comprehensive methodology for impact analysis, contributing valuable insights to the field of radioactive materials packaging and transportation safety.

Martinez, Oscar [ORNL] (ORCID:0000000181814046)↗
DOE OSTI2020

Fast Approximate Methods for Modified Gravity Cosmological Simulations

The accelerated expansion of the Universe poses a major theoretical puzzle. Although the assumption of a non-zero cosmological constant provides a minimal extension of general relativity that is consistent with observational data, many theories of modified gravity have been suggested as possible alternatives. Predictions of structure formation for these models in the fully nonlinear regime are very expensive and it is difficult, if not impossible, to explore such a huge space of models and parameters using high-resolution N-body simulations. Even in the mildly nonlinear regime, perturbative methods can become extremely complex. We explore whether simplified dynamical approximations, applicable for a certain set of cosmological probes, can be used to investigate models of modified gravity with acceptable accuracy in the latter instance. For the case of chameleon gravity, we find that these methods can indeed be used to explore the region around the baryon acoustic oscillation scale, k ~ 0:1 hMpc-1 but not much further.

cosmology↗