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 91 records · Page 5

Dynamic analysis of a flexible spacecraft with rotating components. Volume 1: Analytical developments

Analytical procedures and digital computer code are presented for the dynamic analysis of a flexible spacecraft with rotating components. Topics, considered include: (1) nonlinear response in the time domain, and (2) linear response in the frequency domain. The spacecraft is assumed to consist of an assembly of connected rigid or flexible subassemblies. The total system is not restricted to a topological connection arrangement and may be acting under the influence of passive or active control systems and external environments. The analytics and associated digital code provide the user with the capability to establish spacecraft system nonlinear total response for specified initial conditions, linear perturbation response about a calculated or specified nominal motion, general frequency response and graphical display, and spacecraft system stability analysis.

Bodley, C. S.↗

Time-dependent density functional theory applied to average atom opacity

We focus on studying the opacity of iron, chromium, and nickel plasmas at conditions relevant to experiments carried out at Sandia National Laboratories. In this work, we calculate the photoabsorption cross sections and subsequent opacity for plasmas using linear-response time-dependent density functional theory (TD-DFT). Our results indicate that the physics of channel mixing accounted for in linear-response TD-DFT leads to an increase in the opacity in the bound-free quasicontinuum, where the Sandia experiments indicate that models underpredict iron opacity. However, the increase seen in our calculations is only in the range of 5%–10%. Further, we do not see any change in this trend for chromium and nickel. This behavior indicates that channel mixing effects do not explain the trends in opacity observed in the Sandia experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗

Toward witnessing molecular exciton entanglement from spectroscopy

Entanglement is a defining feature of quantum mechanics that can be a resource in engineered and natural systems, but measuring entanglement in experiment remains elusive especially for large chemical systems. Most practical approaches require determining and measuring a suitable entanglement witness which provides some level of information about the entanglement structure of the probed state. A fundamental quantity of quantum metrology is the quantum Fisher information (QFI), which is a rigorous witness of multipartite entanglement that can be evaluated from linear response functions for certain states. Here, in this paper, we explore measuring the QFI of molecular exciton states of the first-excitation subspace from spectroscopy. In particular, we utilize the fact that the linear response of a pure state subject to a weak electric field over all possible driving frequencies encodes the variance of the collective dipole moment in the probed state, which is a valid measure for QFI. The systems that are investigated include the molecular dimer, N-site linear aggregate with nearest-neighbor coupling, and N-site circular aggregate, all modeled as a collection of interacting qubits. Our theoretical analysis shows that the variance of the collective dipole moment in the brightest dipole-allowed eigenstate is the maximum QFI. The optical response of a thermally equilibrated state in the first-excitation subspace is also a valid QFI. Theoretical predictions of the measured QFI for realistic linear dye aggregates as a function of temperature and energetic disorder due to static variations of the host matrix show that two- to three-partite entanglement is realizable. This paper lays some groundwork and inspires measurement of multipartite entanglement of molecular excitons with ultrafast pump-probe experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Vortex shedding flowmeters for liquids at high flow velocities

A number of vortex shedding flowmeter designs for flow measurements in liquid oxygen ducts on the space shuttle main engines have been tested in a high head water flow test facility. The results have shown that a vortex shedding element or vane spanning the duct can give a linear response to an average flow velocity of 46 m/s (150 ft/s) in a 1 1/2 inch nominal (41 mm actual) diameter duct while a vane partially spanning the duct can give a linear response to velocities exceeding 55 m/s (180 ft/s). The maximum pressure drops across the flow sensing elements extrapolate to less than 0.7 MPa (100 psi) at 56 m/s (184 ft/s) for liquid oxygen. The test results indicate that the vanes probably cannot be scaled up with pipe size, at least not linearly.

Siegwarth, J. D.↗

Nonequilibrium steady-state thermoelectrics of Kondo-correlated quantum dots

The transport across a Kondo-correlated quantum dot coupled to two leads with independent temperatures and chemical potentials is studied using a controlled nonperturbative, and in this sense numerically exact, treatment based on a hybrid numerical renormalization group combined with time-dependent density matrix renormalization group (NRG-tDMRG). In the Kondo regime, for sufficiently large fixed voltage bias V ≳ T K , with T K the Kondo temperature, we find a peak in the conductance vs the temperature gradient Δ⁢T = T R - T L across left and right lead. Focusing then on zero voltage bias but finite ΔT far beyond linear response, we reveal the dependence of the characteristic zero-bias conductance on the individual lead temperatures. Here, we find that the finite-Δ⁢T data behaves quantitatively similar to linear response with an effective equilibrium temperature derived from the different lead temperatures. The regime of sign changes in the Seebeck coefficient, signaling the presence of Kondo correlations, and its dependence on the individual lead temperatures provide a complete picture of the Kondo regime in the presence of finite-temperature gradients. The results from the zero-bias conductance and Seebeck coefficient studies unveil an approximate “Kondo circle” in the T L /T R plane as the regime within which the Kondo correlations dominate. We also study the heat current and the corresponding heat conductance vs finite Δ⁢T. We provide a polynomial fit for our numerical results for the thermocurrent as a function of the individual lead temperatures, which may be used to fit experimental data in the Kondo regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Perspectives on the Calibration of CNN Energy Reconstruction in Highly Granular Calorimeters

We present a study which shows encouraging stability of the response linearity for a simulated high granularity calorimeter module reconstructed by a CNN model to miscalibration, bias, and noise effects. Our results also show an intuitive, quantifiable relationship between these factors and the calibration parameters. We trained a CNN model to reconstruct energy in the calorimeter module using simulated single-pion events; we then observed the response of the model under various miscalibration, bias, and noise conditions that affected the model input. From these data, we estimated linear response models to calibrate the CNN. We also quantified the relationship between these factors and the calibration parameters by regression analysis.

Akchurin, Nural↗

Investigations of the linear and non-linear flow harmonics using the a multi-phase transport model

The multi-phase transport model (AMPT) is used to study the effects of the parton-scattering cross-sections (σ pp ) and hadronic re-scattering on the linear contributions to the flow harmonic v 4 , the non-linear response coefficients, and the correlations between different order flow symmetry planes in Au + Au collisions at 200 GeV. The model results, which agree with current experimental measurements, indicate that the higher-order flow harmonics are sensitive to the σ pp variations. However, the non-linear response coefficients and the correlations between different order flow symmetry planes are σ pp independent. Furthermore, these results suggest that further detailed experimental measurements which span a broad range of collision systems and beam energies could serve as an additional constraint for the theoretical models' calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Linear Time-Invariant Models of a Large Cumulus Ensemble

Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.

Meteorology & Atmospheric Sciences↗

Unified, Geometric Framework for Nonequilibrium Protocol Optimization

Controlling thermodynamic cycles to minimize the dissipated heat is a long-standing goal in thermodynamics, and more recently, a central challenge in stochastic thermodynamics for nanoscale systems. Here, we introduce a theoretical and computational framework for optimizing nonequilibrium control protocols that can transform a system between two distributions in a minimally dissipative fashion. These protocols optimally transport a system along paths through the space of probability distributions that minimize the dissipative cost of a transformation. Furthermore, we show that the thermodynamic metric—determined via a linear response approach—can be directly derived from the same objective function that is optimized in the optimal transport problem, thus providing a unified perspective on thermodynamic geometries. As a result, we investigate this unified geometric framework in two model systems and observe that our procedure for optimizing control protocols is robust beyond linear response.

36 MATERIALS SCIENCE↗

The Leading Modes of Asian Summer Monsoon Variability as Pulses of Atmospheric Energy Flow

Monsoon rainfall anomalies are often organized into dynamical modes that are more predictable than rainfall at specific locations. Here a linear response function (LRF) is constructed to extract the neutral modes (NM) of the dynamical atmospheric system and their optimal forcing that govern precipitation variability, providing a dynamical underpinning for the empirically derived, recurrent rainfall patterns. The leading mode (NM1) features a shift in the Atlantic ITCZ, a drying over the Indian Ocean and Indian subcontinent, and a “South-Flood-North Drought (SFND)” pattern over eastern China, a pattern inferable from changes in the atmospheric energy flow forced by the corresponding optimal forcing and further enhanced by the water vapor and cloud feedbacks. Moreover, the slow evolving component of the NM1 in response to CO2 forcing can be skillfully predicted by running the ocean heat uptake through the LRF, implicating a deterministic ocean dynamic source for the NM1-like trend under climate warming.

monsoon mode of variability, EBM, climate change, ↗

Thermal control across a chain of electronic nanocavities

Here, we study a chain of alternating hot and cold electronic nanocavities—connected to one another via resonant-tunneling quantum dots—with the intent of achieving precise thermal control across the chain. This is accomplished by positioning the dots' energy levels such that a predetermined distribution of heat currents is realized across the chain in the steady state. The number of electrons in each cavity is conserved in the steady state which constrains the cavities' chemical potentials. We determine these chemical potentials analytically in the linear response regime where the energy differences between the dots' resonant levels and the neighboring chemical potentials are much smaller than the thermal energy. In this regime, the thermal control problem can be solved exactly, while, in the general case, thermal control can only be achieved in a relative sense, that is, when one only preassigns the ratios between different heat currents. We apply our results to two different cases: We first demonstrate that a “heat switch” can be easily realized with three coupled cavities, and then we show that our linear response results can provide accurate results in situations with a large number of cavities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effects of Tunable Data Compression on Geophysical Products Retrieved from Surface Radar Observations with Applications to Spaceborne Meteorological Radars

This paper presents results and analyses of applying an international space data compression standard to weather radar measurements that can easily span 8 orders of magnitude and typically require a large storage capacity as well as significant bandwidth for transmission. By varying the degree of the data compression, we analyzed the non-linear response of models that relate measured radar reflectivity and/or Doppler spectra to the moments and properties of the particle size distribution characterizing clouds and precipitation. Preliminary results for the meteorologically important phenomena of clouds and light rain indicate that for a 0.5 dB calibration uncertainty, typical for the ground-based pulsed-Doppler 94 GHz (or 3.2 mm, W-band) weather radar used as a proxy for spaceborne radar in this study, a lossless compression ratio of only 1.2 is achievable. However, further analyses of the non-linear response of various models of rainfall rate, liquid water content and median volume diameter show that a lossy data compression ratio exceeding 15 is realizable. The exploratory analyses presented are relevant to future satellite missions, where the transmission bandwidth is premium and storage requirements of vast volumes of data, potentially problematic.

cloud properties↗

Thermomechanical Postbuckling of Shells

A number of issues involving the postbuckling of doubly-curved panels exposed to a non-uniform temperature field and subjected to a system of mechanical loads are considered in this paper. Effects played by transverse shear, initial geometric imperfections, membrane and thickness-wise temperature gradient, tangential edge constraints, linear/non-linear elastic foundation are incorporated and their implications on the non-linear response behavior of shells are discussed. Issues related with the delay of the occurrence of the snap-through buckling, alleviation of its intensity and even with its removal are addressed and remarks on conditions rendering a curved panel insensitive to initial geometric imperfections are included. Results on small vibration of shells about sub- and post-critical ranges are provided and structural modeling issues playing a non-negligible role towards a reliable prediction of the postbuckling behavior of shells under complex loading conditions are discussed. Results on non-linear response of flat reinforced panels subjected to thermomechanical loads are also presented. Finally, problems not covered by this article are briefly indicated.

Librescu, Liviu↗

Machine learning the Hubbard U parameter in DFT+U using Bayesian optimization

Abstract Within density functional theory (DFT), adding a Hubbard U correction can mitigate some of the deficiencies of local and semi-local exchange-correlation functionals, while maintaining computational efficiency. However, the accuracy of DFT+U largely depends on the chosen Hubbard U values. We propose an approach to determining the optimal U parameters for a given material by machine learning. The Bayesian optimization (BO) algorithm is used with an objective function formulated to reproduce the band structures produced by more accurate hybrid functionals. This approach is demonstrated for transition metal oxides, europium chalcogenides, and narrow-gap semiconductors. The band structures obtained using the BO U values are in agreement with hybrid functional results. Additionally, comparison to the linear response (LR) approach to determining U demonstrates that the BO method is superior.

Chemistry↗

Geoscience Laser Altimeter System: Characteristics and Performance of the Altimeter Receiver

The Geoscience Laser Altimeter System (GLAS) on board ICESat spacecraft measures the surface height (altimetry) via the time of flight of its 1064 nm laser pulse. The GLAS laser transmitter produces 6 ns wide pulses with 70 mJ energy at 1064 nm at a 40 Hz rate. The altimeter receiver consists of a telescope, aft optics, a silicon avalanche photodiode, and electronic amplifiers. The transmitted and echo pulse waveforms are digitized at 1 GHz rate. The laser pulse time of flight is determined on the ground from the two digitized pulse waveforms and their positions in the full waveform record (about 5.4 ms ong) by computing the pulse centroids or by curve fitting. The GLAS receiver algorithms in on board software selects the two waveform segments containing the transmitted and the echo pulses and sends them to ground. The probability of echo pulse detection and the accuracy of time of flight measurement depend on the received signal level, the background light within the receiver field of view, the inherent detector and amplifier noise, the quantization of the digitizer, and some times by cloud obscurations. A receiver model has been developed to calculate the probability of detection and accuracy of the altimeter measurements with these noise sources. From prelaunch testing, the minimum detectable echo pulse energy for 90% detection probability was about 0.1 fj/pulse onto the detector. Such a receiver sensitivity allows GLAS to measure the surface height through clouds with optical density less than 2. The echo pulse energy required to achieve 10 cm ranging accuracy was found to be about 3 times higher than the minimum detectable signal level. The smallest single shot range measurement error, which was determined by ranging to a fixed target with strong echo pulses and no background light, was 2 to 3cm. The maximum linear response echo pulse energy was 10 fJ/pulse for the strongest echo signals, assuming a Lambertian scattering snow surface, clear sky atmosphere transmission, and no pulse width spreading due to surface slopes. Results from the initial 38 days of GLAS operation on-orbit with Laser 1, showed the echo pulse energies varied, as expected, over a wide dynamic range due to the variability of atmosphere transmission and Earth surface characteristics. The receiver was able to reliably detect the ground surface echoes and measure time of flight under such dynamic conditions, except when there were dense clouds. When over Antarctica, the echo pulse energies were several times stronger than predicted. This suggests that the 1064 nm backscattering from the snow pack and ice-sheet &ce are not completely Lambertian but is somewhat peaked at opposit.ion. With Laser 1, the peak amplitude of the echo pulses from flat ice sheets under clear sky conditions exceeded the linear response range of the receiver, causing some pulse waveform distortion due to saturation. We have characterized the effects of receiver saturation on the time of flight, pulse width, and pulse energy measurements for flat surface by testing a night spare detector in the lab with simulated echo pulses. A data processing algorithm that minimizes the errors due to saturation for these measurements are described.

Sun, Xiao-Li↗

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗

QuDPy: A Python-based tool for computing ultrafast non-linear optical responses

Nonlinear Optical Spectroscopy is a well-developed field with theoretical and experimental advances that have benefited multiple disciplines, including chemistry, biology, and physics. However, for the accurate interpretation of the corresponding multi-dimensional spectra, there is a need for precise quantum dynamical simulations based on model Hamiltonians. In this article, we present the initial release of our code, QuDPy (Quantum Dynamics in Python), which provides a robust numerical platform for performing quantum dynamics simulations based on model systems, including open quantum systems. Furthermore, a distinguishing feature of our approach is the ability to specify various high-order optical response pathways in the form of double-sided Feynman diagrams through a straightforward input syntax. This syntax outlines the time-ordering of ket-sided or bra-sided optical interactions acting on the time-evolving density matrix of the system. We utilize the quantum dynamics capabilities of QuTip to simulate the spectral response of complex systems, allowing us to compute virtually any $n$-th order optical response of the model system. To illustrate the utility of our approach, we provide a series of example calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗