SEARCH · Engineering Papers

Results for “linearity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2023

An accurate and efficient Ehrenfest dynamics approach for calculating linear and nonlinear electronic spectra

Linear and nonlinear electronic spectra provide an important tool to probe the absorption and transfer of electronic energy. Here, we introduce a pure state Ehrenfest approach to obtain accurate linear and nonlinear spectra that is applicable to systems with large numbers of excited states and complex chemical environments. We achieve this by representing the initial conditions as sums of pure states and unfolding multi-time correlation functions into the Schrödinger picture. By doing this, we show that one can obtain significant improvements in accuracy over the previously used projected Ehrenfest approach and that these benefits are particularly pronounced in cases where the initial condition is a coherence between excited states. While such initial conditions do not arise when calculating linear electronic spectra, they play a vital role in capturing multidimensional spectroscopies. We demonstrate the performance of our method by showing that it is able to quantitatively capture the exact linear, 2D electronic spectroscopy, and pump–probe spectra for a Frenkel exciton model in slow bath regimes and is even able to reproduce the main spectral features in fast bath regimes.

Chemistry↗

DOE OSTI2023

Linear shaped-charge jet optimization using machine learning methods

Linear shaped charges are used to focus energy into rapidly creating a deep linear incision. The general design of a shaped charge involves detonating a confined mass of high explosive (HE) with a metal-lined concave cavity on one side to produce a high velocity jet for the purpose of striking and penetrating a given material target. This jetting effect occurs due to the interaction of the detonation wave with the cavity geometry, which produces an unstable fluid phenomenon known as the Richtmyer–Meshkov instability and results in the rapid growth of a long narrow jet. We apply machine learning and optimization methods to hydrodynamics simulations of linear shaped charges to improve the simulated jet characteristics. The designs that we propose and investigate in this work generally involve modifying the behavior of the detonation waves prior to interaction with the liner material. These designs include the placement of multiple detonators and the use of metal inclusions within the HE. In conclusion, we are able to produce a linear shaped-charge design with a higher penetration depth than the baseline case that we consider and accomplish this using the same amount of or less HE.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Natural formation of linear defect structures in MoS 2

Near surface defects can significantly impact the quality of metallic interconnects and other interfaces necessary to create device structures incorporating two-dimensional materials. Furthermore, the impact of such defects can strongly depend on their organization. In this study, we present scanning tunneling microscopy images and tunneling spectroscopy of point and linear defects near the surface of natural MoS 2 . The point defects share similar structural and electronic characteristics and occur with comparable frequency as subsurface sulfur vacancies observed previously on natural MoS 2 . The linear defects observed here occur less frequently than the point defects but share the same depth profile and electronic structure. These data indicate that the linear defects are actually a one-dimensional organization of subsurface sulfur vacancies. Our density functional calculations agree with this assessment in that, for sufficient local defect concentrations, it is energetically more favorable for the defects to be organized in a linear fashion rather than as clusters or even isolated single point defects. Given these measurements were taken from naturally formed MoS 2 , this organization likely occurs during crystal formation. Considering the impact of one-dimensional organization on the local properties of layered materials, and the potential for them to be introduced purposefully during crystal formation, research into the formation mechanism and properties of these defects could enable new paths for defect engineering in MoS 2 -based systems.

2D materials↗

DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Gauge-invariant gravitational waves in matter beyond linearized gravity

Modeling the propagation of gravitational waves (GWs) in media other than vacuum is complicated by the gauge freedom of linearized gravity in that, once nonlinearities are taken into consideration, gauge artifacts can cause spurious acceleration of the matter. To eliminate these artifacts, we propose how to keep the theory of dispersive GWs gauge-invariant beyond the linear approximation and, in particular, obtain an unambiguous gauge-invariant expression for the energy–momentum of a GW in a dispersive medium. Using analytic tools from plasma physics, we propose an exactly gauge-invariant 'quasilinear' theory, in which GWs are governed by linear equations and also affect the background metric on scales large compared to their wavelength. As a corollary, the gauge-invariant geometrical optics of linear dispersive GWs in a general background is formulated. As an example, we show how the well-known properties of vacuum GWs are naturally and concisely yielded by our theory in a manifestly gauge-invariant form. We also show how the gauge invariance can be maintained within a given accuracy to an arbitrary order in the GW amplitude. These results are intended to form a physically meaningful framework for studying dispersive GWs in matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

A method to test the coupling strength of the linear and nonlinear contributions to higher-order flow harmonics via event shape engineering

A Multi-Phase Transport (AMPT) model is used to study the efficacy of shape-engineered events to delineate the degree of coupling between the linear and nonlinear contributions to the higher-order flow harmonics v 4 and v 5 . Here, the study shows that the nonlinear contributions are strongly shape-dependent while the linear contributions are shape-independent, indicating little if any, coupling between the linear and nonlinear flow coefficients. The experimental verification of such patterns could be an invaluable tool for robust extraction of the linear and mode-coupled flow coefficients, especially for beam energies where the charged particle multiplicity and the event statistics precludes the use of current methods to establish the coupling strength.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Non-linearities in the Lyman-α forest and in its cross-correlation with dark matter halos

Here three-dimensional correlations of the Lyman-α (Lyα) forest and cross correlations between the Lyα forest and quasars have been measured on large scales, allowing a precise measurement of the baryon acoustic oscillation (BAO) feature at redshifts z > 2. These 3D correlations are often modelled using linear perturbation theory, but full-shape analyses to extract cosmological information beyond BAO will require more realistic models capable of describing non-linearities present at smaller scales. We present a measurement of the Lyα forest flux power spectrum from large hydrodynamic simulations — the Sherwood simulations — and compare it to different models describing the small-scale deviations from linear theory. We confirm that the model presented in Arinyo-i-Prats et al. (2015) fits the measured 3D power up to k = 10 h Mpc -1 with an accuracy better than 5%, and show that the same model can also describe the 1D correlations with similar precision. We also present, for the first time, an equivalent study for the cross-power spectrum of halos with the Lyα forest, and we discuss different challenges we face when modelling the cross-power spectrum beyond linear scales. We make all our measured power spectra public in https://github.com/andreufont/sherwoo_p3d. This study is a step towards joint analyses of 1D and 3D flux correlations, and towards using the quasar-Lyα cross-correlation beyond BAO analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

A direct detection method of galaxy intrinsic ellipticity-gravitational shear correlation in non-linear regimes using self-calibration

Intrinsic alignment (IA) of galaxies is a challenging source of contamination in the Cosmic shear (GG) signals. The galaxy intrinsic ellipticity-gravitational shear (IG) correlation is generally the most dominant component of such contamination for cross-correlating redshift bins. One of the most effective techniques to mitigate such contamination is the self-calibration (SC) method which extracts the IG correlation and allows for its removal from the GG signal. In a photometric survey, the SC method first extracts the galaxy number density-galaxy intrinsic ellipticity (gI) correlation from the observed galaxy-galaxy lensing correlation using the redshift dependence of lens-source pairs. The IG correlation is computed through a scaling relation using the gI correlation and other lensing observables. The applicability of the SC method has so far been focused on the linear IA scales and the linear galaxy bias. We extend the SC method beyond the linear regime by modifying its scaling relation which can account for the non-linear galaxy bias model and various IA models. In this study, we provide a framework to detect the IG correlation for the redshift bins for source galaxies for the proposed year 1 survey of the Rubin Legacy Survey of Space and Time (LSST Y1). We tested the method for the tidal alignment and tidal torquing (TATT) model of IA and we found that the scaling relation is accurate within 10% and 20% for cross-correlating and auto-correlating redshift bins, respectively. Hence the suppression of IG contamination in observed GG correlation can be accomplished with a factor of 10 and 5, for cross-correlating and auto-correlating redshift bins, respectively. We tested the method's robustness and found that the suppression of IG contamination by a factor of 5 is still achievable for all combinations of cross-correlating bins even with the inclusion of a moderate amount of uncertainties on IA and bias parameters, respectively. We also make available, a branch of the code FAST-PT to provide gI correlations up to 1-loop order term used by the new SC method.

gravitational lensing↗

DOE OSTI2020

A digital twin based on OpenFAST linearizations for real-time load and fatigue estimation of land-based turbines

Monitoring a wind turbine requires intensive instrumentation, which would be too cost-prohibitive to deploy to an entire wind plant. This work presents a technique that uses readily available measurements to estimate signals that would otherwise require additional instrumentation. This study presents a digital twin concept with a focus on estimating wind speed, thrust, torque, tower-top position, and loads in the tower using supervisory control and data acquisition (SCADA) measurements. The model combines a linear state-space model obtained using OpenFAST linearizations, a wind speed estimator, and a Kalman filter algorithm that integrates measurements with the state model to perform state estimations. The measurements are: top acceleration, generator torque, pitch, and rotational speed. The article extends previous work that derived the linear state-space model using a different method. The new implementation, based on OpenFAST linearization capability, allows for a systematic extension of the method to more states, inputs, outputs, and to the offshore environment. Results from the two methods are compared, and the validation is made with additional measurements using the GE 1.5-MW turbine located at the National Renewable Energy Laboratory test site. Real-time damage equivalent loads of the tower bottom moment are estimated with an average accuracy of approximately 10%. Overall, the results from this proof of concept are encouraging, and further application of the model will be considered.

17 WIND ENERGY↗

DOE OSTI2023

A non-linear solution to the S 8 tension – II. Analysis of DES Year 3 cosmic shear

ABSTRACT Weak galaxy lensing surveys have consistently reported low values of the S8 parameter compared to the Planck lambda cold dark matter (ΛCDM) cosmology. Amon & Efstathiou used KiDS-1000 cosmic shear measurements to propose that this tension can be reconciled if the matter fluctuation spectrum is suppressed more strongly on non-linear scales than assumed in state-of-the-art hydrodynamical simulations. In this paper, we investigate cosmic shear data from the Dark Energy Survey (DES) Year 3. The non-linear suppression of the matter power spectrum required to resolve the S8 tension between DES and the Planck ΛCDM model is not as strong as inferred using KiDS data, but is still more extreme than predictions from recent numerical simulations. An alternative possibility is that non-standard dark matter contributes to the required suppression. We investigate the redshift and scale dependence of the suppression of the matter power spectrum. If our proposed explanation of the S8 tension is correct, the required suppression must extend into the mildly non-linear regime to wavenumbers $k\sim 0.2 \, h\, {\rm Mpc}^{-1}$. In addition, all measures of S8 using linear scales should agree with the Planck ΛCDM cosmology, an expectation that will be testable to high precision in the near future.

Preston, Calvin (ORCID:0000000281136916)↗

DOE OSTI2025

Variational quantum and neural quantum states algorithms for the linear complementarity problem

Variational quantum algorithms (VQAs) are promising hybrid quantum-classical methods designed to leverage the computational advantages of quantum computing while mitigating the limitations of current noisy intermediate-scale quantum (NISQ) hardware. Although VQAs have been demonstrated as proofs of concept, their practical utility in solving real-world problems—and whether quantum-inspired classical algorithms can match their performance—remains an open question. We present a novel application of the variational quantum linear solver (VQLS) and its classical neural quantum states-based counterpart, the variational neural linear solver (VNLS), as key components within a minimum map Newton solver for a complementarity-based rigid-body contact model. We demonstrate using the VNLS that our solver accurately simulates the dynamics of rigid spherical bodies during collision events. These results suggest that quantum and quantum-inspired linear algebra algorithms can serve as viable alternatives to standard linear algebra solvers for modelling certain physical systems.

neural quantum states↗

DOE OSTI2022

Topological alloy engineering and locally linearized gap dependence on concentration

Alloy engineering is a well-established approach to tune various materials’ properties, but its application to topological alloys remains rudimentary. Of special interest is the band gap, the most defining property of topological materials; however, the concentration dependence of energy gaps in topological alloys remains unknown. Here we systematically investigate the band gap evolution of a topological alloy as a function of alloy concentration, using KZnSb 1-x Bi x as a prototype, based on first-principles calculations. In contrast to the well-established smooth bowing curve for a trivial gap in semiconductor alloys, we found that the topological gap evolves generally with a complex fragmented pattern due to topological phase transitions, and most strikingly a linear dependence on concentration locally in each distinct phase. Such gap linearization is fundamentally rooted in the linear dependence on alloy concentration of spin-orbit coupling (SOC) that predominantly determines a topological gap. Furthermore, we demonstrate topological alloy engineering as a general approach to tune the topological order by modulating the band edge composition and degeneracy through the alloying-induced interplay of SOC and atomic orbital on-site energy, while the linear gap dependence on alloy concentration remains independent of the degree of topological order.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Controlled skyrmion ratchet in linear protrusion defects

Using atomistic simulations, we investigate the dynamical behavior of a single skyrmion interacting with an asymmetric linear protrusion array under external ac driving. The linear protrusion array is placed in the xy plane and is composed of magnetic walls near the edge of the nanotrack that have angled linear protrusions into the nanotrack. The structure forms a hard axis for –x direction motion and an easy axis for +x direction motion. When the ac drive is applied along the x direction, the skyrmion moves along the hard direction of the substrate asymmetry in three phases: a pinned phase with localized skyrmion orbits, a constant velocity phase where the orbits become delocalized, and a reentrant pinned phase with larger localized orbits. We measure the dependence of the skyrmion velocity on the frequency and amplitude of the ac drive. All three phases appear for all frequency values simulated here, and in the constant velocity phase the skyrmion velocity depends only on the frequency and not on the amplitude of the ac drive. When ac driving is applied in the y direction, the skyrmion moves along the easy direction of the substrate asymmetry and exhibits the same three phases as for x-direction driving along with a fourth phase which, at high driving frequencies, consists of a series of constant velocity phases, each with different average skyrmion velocities. For low frequencies, the constant velocity phase is lost and the skyrmion speed increases linearly with increasing ac drive amplitude due to a Magnus boost effect. Finally, our findings suggest new ways to create reliable data transport for spintronic devices using skyrmions as information carriers, where the skyrmion direction and speed can be controlled by varying only the ac drive amplitude and frequency.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Temperature scaling behavior of the linear magnetoresistance observed in high-temperature superconductors

An analytical model invoking variations in the charge-carrier density is used to generate magnetoresistance curves that are almost indistinguishable from those produced by sophisticated numerical models. This demonstrates that, though disorder is pivotal in causing linear magnetoresistance, the form of the magnetoresistance thus generated is insensitive to details of the disorder. Additionally, taken in conjunction with the temperature (T) dependence of the zero-field resistivity, realistic levels of disorder are shown to be sufficient to explain the linear magnetoresistance and field-T resistance scaling observed in high-temperature pnictide and cuprate superconductors. Hence, though the T-linear zero-field resistance is a definite signature of the “strange metal” state of high-temperature superconductors, their linear magnetoresistance and its scaling is unlikely to be so.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Studying the Peak Force of Magnetic Linear Motion Devices Using Dimensional Neutrality

The merits of different sizing equations for linear motion magnetic force devices are assessed. A peak force scaling analysis for a linear magnetic coupling with respect to the volumetric force density, mass force density, magnetic shear stress and energy density is compared. By utilizing a 3-D analytic based model of a Halbach cylinder magnetic coupling it is shown that the maximizing of volumetric and mass force density is not an appropriate sizing metric for a linear magnetic actuator. Only energy density and magnetic shear stress are dimensionally neutral sizing metrics for linear actuators.

16 TIDAL AND WAVE POWER↗

DOE OSTI2024

Linear Solvers for Collector Systems of Generalized Large-scale Inverter-Based Resources

Collector systems for inverter-based resources (IBRs) are typically represented by equivalent circuits for electromagnetic transient (EMT) simulations. Recent studies have revealed that modeling a detailed collector system is essential to accurately represent the behavior of IBRs, especially when dealing with partial tripping during external disturbances. However, there are several challenges in simulating a detailed EMT model of a collector system due to the time required to simulate such systems. Thus, this paper investigates the modeling of a detailed collector system, taking into account its configuration and components as defined in IEEE standard 2800. The configurations include the collector systems of generalized large-scale IBR plants. The components include the main IBR transformer, collector bus, and feeders with lines and/or cables. The EMT model of the collector system is represented by differential algebraic equations (DAEs) that are discretized to form linear equations that are solved using linear solvers. In this paper, linear solvers are proposed based on the Schur complement method, which are utilized for simulation of the EMT model of collector systems of generalized large-scale IBRs to accelerate simulation speed while maintaining the accuracy of the results. The proposed solvers are verified by comparing the performance to that of linear solvers provided in MATLAB.

Choi, Jongchan↗

DOE OSTI2021

Learning Optimal Power Flow Solutions using Linearized Models in Power Distribution Systems

Solving nonlinear optimal power flow (OPF) problem is computationally expensive, and poses scalability challenges for power distribution networks. An alternative to solving the original nonlinear OPF is the linear approximated OPF models. Although, these linear approximated OPF models are fast, the resulting solutions may result in significant optimality gap. Lately, the application of machine learning (ML) methods in successfully solving the nonlinear OPF has been reported. These methods learn and estimate the nonlinear control policies using a purely data-driven approach. In this paper, we propose an approach to complements the ML based approach to solving OPF using solutions from known linearized OPF model. Specifically, we use supervised learning to map the solutions of linear OPF to nonlinear control variables. Unlike, the traditional ML based methods for OPF that approximate the full distribution feeder model using function approximation, our approach uses a two-node approximation of radial networks. The proposed approach is validated using IEEE 123 bus test system for OPF solutions obtained using the nonlinear OPF models.

optimal power flow, power distribution systems, su↗

DOE OSTI2024

Estimation and Visualization of Isosurface Uncertainty from Linear and High-Order Interpolation Methods

Isosurface visualization is fundamental for exploring and analyzing 3D volumetric data. Marching cubes (MC) algorithms with linear interpolation are commonly used for isosurface extraction and visualization. Although linear interpolation is easy to implement, it has limitations when the underlying data is complex and high-order, which is the case for most real-world data. Linear interpolation can output vertices at the wrong location. Its inability to deal with sharp features and features smaller than grid cells can lead to an incorrect isosurface with holes and broken pieces. Despite these limitations, isosurface visualizations typically do not include insight into the spatial location and the magnitude of these errors. We utilize high-order interpolation methods with MC algorithms and interactive visualization to highlight these uncertainties. Our visualization tool helps identify the regions of high interpolation errors. It also allows users to query local areas for details and compare the differences between isosurfaces from different interpolation methods. In addition, we employ high-order methods to identify and reconstruct possible features that linear methods cannot detect. We showcase how our visualization tool helps explore and understand the extracted isosurface errors through synthetic and real-world data.

Ouermi, Timbwaoga↗