Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear response”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Optimal Storage Response to Utility Tariff Structures and Potential Use of Capacity Charges

Energy storage is increasingly being deployed in behind-the-meter use cases, partially in response to falling lithium-ion prices and new utility tariff structures. Utilities are grappling with new tariff design in the presence of distributed energy resources (DER), trying to motivate, fairly price the contribution of, and, in some cases, discourage, certain operation of DER. There are opportunities for energy storage under emerging tariff structures, but utility net load management objectives for storage remain unclear. Diverse tariff structures motivate the use of energy storage to provide one or a combination of: energy arbitrage, energy shifting, solar self-consumption, and import shaving. This paper formulates a linear optimization to demonstrate the optimal storage tariff response, examining customer net load metrics under diverse utility tariff structures, such as time-of-use, net metering, feed-in tariffs, zero export tariffs, and demand charges. A key contribution of this paper is the introduction and examination of a capacity charge as mechanism that can both motivate reductions in both peak import, and exports, along with greater load levelling.

behind-the-meter↗
DOE OSTI2022

Robust trajectory-constrained frequency control for microgrids considering model linearization error

Grid supportive modes integrated within inverter-based resources can improve the frequency response of renewable-rich microgrids. The synthesis of grid supportive modes to guarantee frequency trajectory constraints under a predefined disturbance set is challenging but essential. To tackle this challenge, a numerical optimal control (NOC)-based control synthesis methodology is proposed. Without loss of generality, a wind-diesel fed microgrid is studied, where we aim to design grid supportive functions in the wind turbine. In the control design, linearized models are used, and the linearization-induced errors are quantitatively analyzed by reachability and interval arithmetics and represented in the form of interval uncertainties. Then, the NOC problem can be formulated into a robust mixed-integer linear program. The control structure is strategically configured into two levels to realize online deployment. The proposed control is verified on the modified 33-node microgrid with a full-order three-phase nonlinear model in Simulink. In conclusion, the simulation results show the effectiveness of the proposed control paradigm and the necessity of considering linearization-induced uncertainty.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Unveiling the effect of composition on nuclear waste immobilization glasses’ durability by nonparametric machine learning

Abstract Ensuring the long-term chemical durability of glasses is critical for nuclear waste immobilization operations. Durable glasses usually undergo qualification for disposal based on their response to standardized tests such as the product consistency test or the vapor hydration test (VHT). The VHT uses elevated temperature and water vapor to accelerate glass alteration and the formation of secondary phases. Understanding the relationship between glass composition and VHT response is of fundamental and practical interest. However, this relationship is complex, non-linear, and sometimes fairly variable, posing challenges in identifying the distinct effect of individual oxides on VHT response. Here, we leverage a dataset comprising 654 Hanford low-activity waste (LAW) glasses across a wide compositional envelope and employ various machine learning techniques to explore this relationship. We find that Gaussian process regression (GPR), a nonparametric regression method, yields the highest predictive accuracy. By utilizing the trained model, we discern the influence of each oxide on the glasses’ VHT response. Moreover, we discuss the trade-off between underfitting and overfitting for extrapolating the material performance in the context of sparse and heterogeneous datasets.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quapprox: A Framework for Benchmarking the Approximability of Variational Quantum Circuit

Most of the existing quantum neural network models, such as variational quantum circuits (VQCs), are limited in their ability to explore the non-linear relationships in input data. This gradually becomes the main obstacle for it to tackle realistic applications, such as natural language processing, medical image processing, and wireless communications. Recently, there have emerged research efforts that enable VQCs to perform non-linear operations. However, it is still unclear on the approximability of a given VQC (i.e., the order of non-linearity that can be handled by a specified design). In response to this issue, we developed an automated tool designed to benchmark the approximation of a given VQC. The proposed tool will generate a set of synthetic datasets with different orders of non-linearity and train the given VQC on these datasets to estimate their approximability. Our experiments benchmark VQCs with different designs, where we know their theoretic approximability. We then show that the proposed tool can precisely estimate the approximability, which is consistent with the theoretic value, indicating that the proposed tool can be used for benchmarking the approximability of a given quantum circuit for learning tasks.

artificial intelligence↗
DOE OSTI2024

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

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

Ladder symmetries of black holes. Implications for love numbers and no-hair theorems

It is well known that asymptotically flat black holes in general relativity have a vanishing static, conservative tidal response. Here, we show that this is a result of linearly realized symmetries governing static (spin 0,1,2) perturbations around black holes. The symmetries have a geometric origin: in the scalar case, they arise from the (E)AdS isometries of a dimensionally reduced black hole spacetime. Underlying the symmetries is a ladder structure which can be used to construct the full tower of solutions, and derive their general properties: (1) solutions that decay with radius spontaneously break the symmetries, and must diverge at the horizon; (2) solutions regular at the horizon respect the symmetries, and take the form of a finite polynomial that grows with radius. Taken together, these two properties imply that static response coefficients — and in particular Love numbers — vanish. Moreover, property (1) is consistent with the absence of black holes with linear (perturbative) hair. We also discuss the manifestation of these symmetries in the effective point particle description of a black hole, showing explicitly that for scalar probes the worldline couplings associated with a non-trivial tidal response and scalar hair must vanish in order for the symmetries to be preserved.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Growth-controlled twinning and magnetic anisotropy in CeSb 2

Cerium diantimonide (CeSb 2 ) is a layered heavy-fermion Kondo lattice material that hosts complex magnetism and pressure-induced superconductivity. The interpretation of its in-plane anisotropy has remained unsettled due to structural twinning, which superimposes orthogonal magnetic responses. Here, in this study, we combine controlled crystal growth with magnetization and rotational magnetometry to disentangle the effects of twinning. Nearly untwinned high-quality single crystals reveal the intrinsic in-plane anisotropy: The in-plane easy axis saturates at 𝑀 easy ⁡(4 T) ≈ 1.8 𝜇 B /Ce, while the in-plane hard axis magnetization is strongly suppressed, nearly linear, and comparable to the out-of-plane response. These results resolve long-standing discrepancies in reported magnetic measurements, in which in-plane metamagnetic transition fields and saturation magnetization varied significantly across previous studies. Growth experiments demonstrate that avoiding the proposed 𝛼 𝛽 structural transition—through Sb-rich flux and slower cooling—systematically reduces twinning. However, powder x-ray diffraction and differential thermal analysis measurements show no clear evidence of a distinct 𝛽 phase. Our results establish a consistent magnetic phase diagram and provide essential constraints for crystal-electric field models, enabling a clearer understanding of the interplay between anisotropic magnetism and unconventional superconductivity in CeSb 2 .

Weber, Jan T. [Goethe University, Frankfurt (Germa↗
DOE OSTI2022

Broadband Asymmetric Transmission of Linearly Polarized Mid‐Infrared Light Based on Quasi‐3D Metamaterials

Abstract Metamaterials consisting of subwavelength resonators offer an exciting opportunity for realizing asymmetric transmission (AT) of linearly polarized light. However, to date, only moderate/narrow‐band AT responses have been obtained in metadevices based on stacked planar nanostructures. Here, leveraging a combination of a genetic algorithm (GA) based optimization method and a membrane projection lithography (MPL) fabrication approach, a quasi‐3D metamaterial for broadband AT of linearly polarized mid‐infrared light is demonstrated. Facilitated by the customized GA, an efficient exploration of 3D plasmonic meta‐atoms with broken mirror symmetry in the light propagation direction allows the satisfaction of the rigorous conditions for AT of linearly polarized waves over a broad wavelength range. Confirmed by surface current analysis, the observed AT behavior is attributed to the resonant coupling between the plasmonic nanostructures located on the two orthogonal walls of the MPL cavities. Incorporating an advanced inverse‐design method and a state‐of‐art fabrication technique, the methodology used in the present study provides a promising route for exploiting 3D metamaterials with sophisticated functionalities via effectively exploring the high‐dimensional parametric space offered by true 3D meta‐atoms.

Whiting, Eric B.↗
DOE OSTI2020

DNA Reaction–Diffusion Attractor Patterns

Abstract Living systems can form and recover complex chemical patterns with precisely sized features in the ranges of tens or hundreds of microns. We show how designed reaction–diffusion processes can likewise produce precise patterns, termed attractor patterns, that reform their precise shape after being perturbed. We use oligonucleotide reaction networks, photolithography, and microfluidic delivery to form precisely controlled attractor patterns and study the responses of these patterns to different localized perturbations. Linear and “hill”‐shaped patterns formed and stabilized into shapes and at time scales consistent with reaction–diffusion models. When patterns were perturbed in particular locations with UV light, they reliably reformed their steady‐state profiles. Recovery also occurred after repeated perturbations. By designing the far‐from‐equilibrium dynamics of a chemical system, this study shows how it is possible to design spatial patterns of molecules that are sustained and regenerated by continually evolving towards a specific steady state configuration.

Dorsey, Phillip James↗
DOE OSTI2020

DNA Reaction–Diffusion Attractor Patterns

Living systems can form and recover complex chemical patterns with precisely sized features in the ranges of tens or hundreds of microns. We show how designed reaction-diffusion processes can likewise produce precise patterns, termed attractor patterns, that reform their precise shape after being perturbed. We use oligonucleotide reaction networks, photolithography and microfluidic delivery to form precisely controlled attractor patterns and study the responses of these patterns to different localized perturbations. Linear and ‘hill’-shaped patterns formed and stabilized into shapes and at time scales consistent with reaction-diffusion models. When patterns were perturbed in particular locations with UV light, they reliably reformed their steady state profiles. Recovery also occurred after repeated perturbations. As a result, by designing the far-from-equilibrium dynamics of a chemical system, this study shows how it is possible to design spatial patterns of molecules that are sustained and regenerated by continually evolving towards a specific steady state configuration

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

A variational phase-field model For ductile fracture with coalescence dissipation

A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. Additionally, the model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.

42 ENGINEERING↗
DOE OSTI2024

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Autonomous Polycrystalline Material Decomposition For Hyperspectral Neutron Tomography

Hyperspectral neutron tomography is an effective method for analyzing crystalline material samples with complex compositions in a non-destructive manner. Since the counts in the hyperspectral neutron radiographs directly depend on the neutron cross-sections, materials may exhibit contrasting neutron responses across wavelengths. Therefore, it is possible to extract the unique signatures associated with each material and use them to separate the crystalline phases simultaneously.We introduce an autonomous material decomposition (AMD) algorithm to automatically characterize and localize polycrystalline structures using Bragg edges with contrasting neutron responses from hyperspectral data. The algorithm estimates the linear attenuation coefficient spectra from the measured radiographs and then uses these spectra to perform polycrystalline material decomposition and reconstructs 3D material volumes to localize materials in the spatial domain. Our results demonstrate that the method can accurately estimate both the linear attenuation coefficient spectra and associated reconstructions on both simulated and experimental neutron data.

Samin nur chowdhury, Mohammad↗
DOE OSTI2021

Multilane Automated Driving With Optimal Control and Mixed-Integer Programming

Road vehicle lane changes often initiate traffic disturbances and can therefore impact road networks’ energy and time efficiency. Furthermore, unexpected changes in traffic conditions may also render lane changes counterproductive for the lane-changing vehicle. Vehicle-to-vehicle connectivity combined with anticipative control could address these challenges via improved lane change decisions by automated vehicles. In a move toward this objective, receding horizon control cast as a mixed-integer quadratic program is used to plan lane changing and acceleration in a coupled optimization. A long-term pacing module, based on Pontryagin’s minimum principle from optimal control theory, sets terminal and input references for receding horizon control to target a user’s expected travel time. To remove nonlinear vehicle dynamics from the receding horizon controller, lane change commands are passed to a pure pursuit steering module whose response is approximated by a second-order linear model. Here, comparison against a rule-based reactive algorithm in arterial and highway scenarios shows an 8.9%–13.7% reduction in energy consumption and a 5.2%–10.3% reduction in the travel time, along with navigational improvements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

Work supported by the grant, "Stability and Control of Burning Tokamak Plasmas with Resistive Walls" in the period Aug 2017 to 2023 can be cast in three general categories: (1) linear tearing regimes and the effects on the response to error fields and resistive walls: (2) nonlinear tearing modes in the presence of error fields, with locking; and (3) energetic particle effects. In the first category, we have completed a study of tearing mode regimes with extended MHD effects, including parallel dynamics and toroidal effects, notably favorable average curvature. This work has been presented at the 2018 Sherwood Conference. This work has been published in the ArXiv and is being prepared for publication in Physics of Plasmas. A systematic study of shaping in a tokamak plasma with resistive tearing layers and a resistive wall has been presented in the PhD dissertation by Dov Rhodes, a student of Andrew Cole at Columbia University. John Finn was on Rhodes's PhD committee, and was a co-author of a paper published in Physics of Plasmas on that work. A detailed study of nonlinear tearing modes in the presence of finite beta, favorable average curvature, and an error field has been nearly completed. In this study it was found that the response to an error field with real frequency a tearing layer agreed with analytic results, and that tearing mode saturation by flattening of the current profile as well as destabilization due to flattening of the pressure profile were important factors in the locking process. This work was presented at the 2017 APS/DPP Conference and at the 2018 Sherwood Conference. The effects of energetic ions on tearing modes was presented at the 2017 APS/DPP Conference and in the PhD dissertation of Michael Halfmoon. John Finn served as a committee member on Halfmoon's PhD committee. Extensions to this work including two fluid effects were performed and presented at the APS/DPP Conference in 2017.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗