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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Properties of cerebellar fastigial neurons during translation, rotation, and eye movements

The most medial of the deep cerebellar nuclei, the fastigial nucleus (FN), receives sensory vestibular information and direct inhibition from the cerebellar vermis. We investigated the signal processing in the primate FN by recording single-unit activities during translational motion, rotational motion, and eye movements. Firing rate modulation during horizontal plane translation in the absence of eye movements was observed in all non-eye-movement-sensitive cells and 26% of the pursuit eye-movement-sensitive neurons in the caudal FN. Many non-eye-movement-sensitive cells recorded in the rostral FN of three fascicularis monkeys exhibited convergence of signals from both the otolith organs and the semicircular canals. At low frequencies of translation, the majority of these rostral FN cells changed their firing rates in phase with head velocity rather than linear acceleration. As frequency increased, FN vestibular neurons exhibited a wide range of response dynamics with most cells being characterized by increasing phase leads as a function of frequency. Unlike cells in the vestibular nuclei, none of the rostral FN cells responded to rotational motion alone, without simultaneously exhibiting sensitivity to translational motion. Modulation during earth-horizontal axis rotation was observed in more than half (77%) of the neurons, although with smaller gains than during translation. In contrast, only 47% of the cells changed their firing rates during earth-vertical axis rotations in the absence of a dynamic linear acceleration stimulus. These response properties suggest that the rostral FN represents a main processing center of otolith-driven information for inertial motion detection and spatial orientation.

Non-NASA Center↗
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2021

Retrieval Algorithm for the Column CO2 Mixing Ratio from Pulsed Multi-Wavelength Lidar Measurements

The retrieval algorithm for the column mixing ratio of CO2 from the measurements of a pulsed multi-wavelength integrated path differential absorption (IPDA) lidar is described. The lidar samples the shape of the 1572.33 nm CO2 absorption line at 15 or 30 wavelengths. The algorithm uses a least-squares fit between the CO2 line shape computed from a layered 10 atmosphere model to that sampled by the lidar. In addition to the column average CO2 dry air mole fraction (XCO2), several other parameters are also solved simultaneously from the fit. These include the Doppler shift in the received laser signal wavelengths, the product of the surface reflectivity and atmospheric transmission and a linear trend in the lidar receiver’s spectral response. The algorithm can also be used to solve for the average water vapor mixing ratio, which causes a secondary absorption in the wings of the CO2 absorption line under high humidity conditions. The least-squares fit is linearized about the 15 expected XCO2 value which allows the use of a standard linear least-squares fitting method and software tools. The standard deviation of the retrieved XCO2 is obtained from covariance matrix of the fit. An averaging kernel is defined similarly to that used for passive trace-gas sounding. Examples are presented of using the algorithm to retrieve XCO2 from the measurements from NASA Goddard’s airborne CO2 Sounder lidar made at a constant altitude and during spiral-down maneuvers.

Xiaoli Sun↗
NASA NTRS1985

Characteristics pertaining to a stiffness cross-coupled Jeffcott model

Rotordynamic studies of complex systems utilizing multiple degree-of-freedom analysis have been performed to understand response, loads, and stability. In order to understand the fundamental nature of rotordynamic response, the Jeffcott rotor model has received wide attention. The purpose of this paper is to provide a generic rotordynamic analysis of a stiffness cross-coupled Jeffcott rotor model to illustrate characteristics of a second order stiffness-coupled linear system. The particular characteristics investigated were forced response, force vector diagrams, response orbits, and stability. Numerical results were achieved through a fourth order Runge-Kutta method for solving differential equations and the Routh Hurwitz stability criterion. The numerical results were verified to an exact mathematical solution for the steady state response.

Spanyer, K. L.↗
DOE OSTI2020

Linking the Dynamic Response of the Carbon Dioxide-Concentrating Mechanism to Carbon Assimilation Behavior in Fremyella diplosiphon

ABSTRACT Cyanobacteria use a carbon dioxide (CO 2 )-concentrating mechanism (CCM) that enhances their carbon fixation efficiency and is regulated by many environmental factors that impact photosynthesis, including carbon availability, light levels, and nutrient access. Efforts to connect the regulation of the CCM by these factors to functional effects on carbon assimilation rates have been complicated by the aqueous nature of cyanobacteria. Here, we describe the use of cyanobacteria in a semiwet state on glass fiber filtration discs—cyanobacterial discs—to establish dynamic carbon assimilation behavior using gas exchange analysis. In combination with quantitative PCR (qPCR) and transmission electron microscopy (TEM) analyses, we linked the regulation of CCM components to corresponding carbon assimilation behavior in the freshwater, filamentous cyanobacterium Fremyella diplosiphon . Inorganic carbon (C i ) levels, light quantity, and light quality have all been shown to influence carbon assimilation behavior in F. diplosiphon . Our results suggest a biphasic model of cyanobacterial carbon fixation. While behavior at low levels of CO 2 is driven mainly by the C i uptake ability of the cyanobacterium, at higher CO 2 levels, carbon assimilation behavior is multifaceted and depends on C i availability, carboxysome morphology, linear electron flow, and cell shape. Carbon response curves (CRCs) generated via gas exchange analysis enable rapid examination of CO 2 assimilation behavior in cyanobacteria and can be used for cells grown under distinct conditions to provide insight into how CO 2 assimilation correlates with the regulation of critical cellular functions, such as the environmental control of the CCM and downstream photosynthetic capacity. IMPORTANCE Environmental regulation of photosynthesis in cyanobacteria enhances organismal fitness, light capture, and associated carbon fixation under dynamic conditions. Concentration of carbon dioxide (CO 2 ) near the carbon-fixing enzyme RubisCO occurs via the CO 2 -concentrating mechanism (CCM). The CCM is also tuned in response to carbon availability, light quality or levels, or nutrient access—cues that also impact photosynthesis. We adapted dynamic gas exchange methods generally used with plants to investigate environmental regulation of the CCM and carbon fixation capacity using glass fiber-filtered cells of the cyanobacterium Fremyella diplosiphon . We describe a breakthrough in measuring real-time carbon uptake and associated assimilation capacity for cells grown in distinct conditions (i.e., light quality, light quantity, or carbon status). These measurements demonstrate that the CCM modulates carbon uptake and assimilation under low-C i conditions and that light-dependent regulation of pigmentation, cell shape, and downstream stages of carbon fixation are critical for tuning carbon uptake and assimilation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The effect of non-linear mutual friction on pulsar glitch sizes and rise times

ABSTRACT Observations of pulsar glitches have the potential to provide constraints on the dynamics of the high density interior of neutron stars. However, to do so, realistic glitch models must be constructed and compared to the data. We take a step towards this goal by testing non-linear models for the mutual friction force, which is responsible for the exchange of angular momentum between the neutron superfluid and the observable normal component in a glitch. In particular, we consider a non-linear dependence of the drag force on the relative velocity between superfluid vortices and the normal component, in which the contributions of both kelvin and phonon excitations are included. This non-linear model produces qualitatively new features, and is able to reproduce the observed bimodal distribution of glitch sizes in the pulsar population. The model also suggests that the differences in size distributions in individual pulsars may be due to the glitches being triggered in regions with different pinning strengths, as stronger pinning leads to higher vortex velocities and a qualitatively different mutual friction coupling with respect to the weak pinning case. Glitches in pulsars that appear to glitch quasi-periodically with similar sizes may thus be due to the same mechanisms as smaller events in pulsars that have no preferred glitch size, but simply originate in stronger pinning regions, possibly in the core of the star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Peak Power Minimization for Commercial Thermostatically Controlled Loads in Multi-Unit Grid-Interactive Efficient Buildings

The load profiles of most commercial and industrial consumers are characterized by brief periods of very high power consumption followed by intervals of lower demand. To encourage such consumers to flatten their load profiles, power utilities in and around the world often levy a monthly demand charge (DC) on the peak demand measured over brief intervals. In this work, we consider the joint optimization of energy costs (EC) and the instantaneous peak power of a multi-unit building which uses a hydronic heating, ventilation and cooling (HVAC) system and responds to a demand response (DR) program. Despite the non-linear structure of the problem, we show how optimal solutions can be obtained efficiently using linear programming. Next, we study the power demand patterns resulting from our proposed strategy for thermostatically controlled loads (TCLs), and evaluate the strategy’s performance for various climate zones in the US, under both typical and atypical weather conditions. Finally, the results show that depending on the ambient conditions and the tariff structure, our strategy can result in utility bill savings of up to nearly 19% compared to the baseline. The results also indicate that our power control strategy can significantly reduce the instantaneous peak power consumption in commercial TCLs.

HVAC↗
NASA NTRS2002

Context-specific adaptation of the gain of the oculomotor response to lateral translation using roll and pitch head tilts as contexts

Previous studies established that vestibular and oculomotor behaviors can have two adapted states (e.g., gain) simultaneously, and that a context cue (e.g., vertical eye position) can switch between the two states. The present study examined this phenomenon of context-specific adaptation for the oculomotor response to interaural translation (which we term "linear vestibulo-ocular reflex" or LVOR even though it may have extravestibular components). Subjects sat upright on a linear sled and were translated at 0.7 Hz and 0.3 gpeak acceleration while a visual-vestibular mismatch paradigm was used to adaptively increase (x2) or decrease (x0) the gain of the LVOR. In each experimental session, gain increase was asked for in one context, and gain decrease in another context. Testing in darkness with steps and sines before and after adaptation, in each context, assessed the extent to which the context itself could recall the gain state that was imposed in that context during adaptation. Two different contexts were used: head pitch (26 degrees forward and backward) and head roll (26 degrees or 45 degrees, right and left). Head roll tilt worked well as a context cue: with the head rolled to the right the LVOR could be made to have a higher gain than with the head rolled to the left. Head pitch tilt was less effective as a context cue. This suggests that the more closely related a context cue is to the response being adapted, the more effective it is.

Non-NASA Center↗
DOE OSTI2025

Physics augmented machine learning discovery of composition-dependent constitutive laws for 3D printed digital materials

Multi-material 3D printing, particularly through polymer jetting, enables the fabrication of digital materials by mixing distinct photopolymers at the micron scale within a single build to create a composite with tunable mechanical properties. Here, this work presents an integrated experimental and computational investigation into the composition-dependent mechanical behavior of 3D printed digital materials. We experimentally characterize five formulations, combining soft and rigid UV-cured polymers under uniaxial tension and torsion across three strain and twist rates. The results reveal nonlinear and rate-dependent responses that strongly depend on composition. To model this behavior, we develop a physics-augmented neural network (PANN) that combines a partially input convex neural network (pICNN) for learning the composition-dependent hyperelastic strain energy function with a quasi-linear viscoelastic (QLV) formulation for time-dependent response. The pICNN ensures convexity with respect to strain invariants while allowing non-convex dependence on composition. To enhance interpretability, we apply $L_0$ sparsification. For the time-dependent response, we introduce a multilayer perceptron (MLP) to predict viscoelastic relaxation parameters from composition. The proposed model accurately captures the nonlinear, rate-dependent behavior of 3D printed digital materials in both uniaxial tension and torsion, achieving high predictive accuracy for interpolated material compositions. This approach provides a scalable framework for automated, composition-aware constitutive model discovery for multi-material 3D printing.

Constitutive modeling↗
NASA NTRS1977

Effects of random phase changes on the formation of synthetic aperture radar imagery

The effects of Gaussian random fluctuation and linear change of the echo phase on the response of the azimuth matched processor of a synthetic aperture radar (airborne or spaceborne for ocean-surface sensing) is analyzed numerically. It is found that the random fluctuation would not alter appreciably the processor response if its standard deviation is less than pi, while its normalized correlation time is larger than 1. If these two bounding conditions are not satisfied, then the sidelobe level becomes relatively large and overshadows the central peak. In the case of linear displacement it is found that the main effect is the spatial displacement of the point scatterer in the image plane.

Elachi, C.↗
NASA NTRS1978

Investigation of ODE integrators using interactive graphics

Two FORTRAN programs using an interactive graphic terminal to generate accuracy and stability plots for given multistep ordinary differential equation (ODE) integrators are described. The first treats the fixed stepsize linear case with complex variable solutions, and generates plots to show accuracy and error response to step driving function of a numerical solution, as well as the linear stability region. The second generates an analog to the stability region for classes of non-linear ODE's as well as accuracy plots. Both systems can compute method coefficients from a simple specification of the method. Example plots are given.

Brown, R. L.↗
NASA NTRS1978

Investigation of ODE integrators using interactive graphics

In the present paper, two FORTRAN programs are described which use an interactive graphic terminal to generate accuracy and stability plots for a given multistep integrator of ordinary differential equations (ODE). One program treats the fixed-stepsize linear case with complex variable solutions, and generates plots to show accuracy and error response to step driving function of a numerical solution. It also generates the linear stability region. The other program generates an analog with an accuracy plot. Both systems are capable of computing method coefficients from a simple specification of the method. Example plots are given.

Brown, R. L.↗
NASA NTRS1984

Coupling of the solar wind to measures of magnetic activity

Linear prediction filtering was used to generate empirical response functions relating the solar wind electric field to the magnetic indices, AL, AU, Dst and ASYM. The empirical response functions were convolved with solar wind observations obtained during the International Magnetospheric Study to predict the indices. The predictions are compared with the observed indices during two, 3-day intervals. Differences between the observed and predicted indices are discussed in terms of the linear assumption and in terms of physical processes other than direct solar wind-magnetosphere interaction.

Mcpherron, R. L.↗
DOE OSTI2021

Adaptive Critical Balance and Firehose Instability in an Expanding, Turbulent, Collisionless Plasma

Abstract Using a hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvénic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expansion (and consequent decrease in the mean magnetic-field strength) gradually reduces the effective elasticity of the field lines, causing reductions in the linear frequency and residual energy of the Alfvénic fluctuations. In response, these fluctuations modify their interactions and spatial anisotropy to maintain a scale-by-scale “critical balance” between their characteristic linear and nonlinear frequencies. Eventually the plasma becomes unstable to kinetic firehose instabilities, which excite rapidly growing magnetic fluctuations at ion-Larmor scales. The consequent pitch-angle scattering of particles maintains the temperature anisotropy near marginal stability, even as the turbulent plasma continues to expand. The resulting evolution of parallel and perpendicular temperatures does not satisfy double-adiabatic conservation laws, but is described accurately by a simple model that includes anomalous scattering. Our results have implications for understanding the complex interplay between macro- and microscale physics in various hot, dilute, astrophysical plasmas, and offer predictions concerning power spectra, residual energy, ion-Larmor-scale spectral breaks, and non-Maxwellian features in ion distribution functions that may be tested by measurements taken in high-beta regions of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Inverse modeling of circular lattices via orbit response measurements in the presence of degeneracy

The number and location of beam position monitors (BPMs) and steerers with respect to the quadrupoles in a circular lattice can lead to degeneracy in the context of fitting linear optics and extracting lattice information from measured closed orbits. Furthermore, the measurement uncertainties due to the imperfection of BPMs and steerers can be propagated by the fitting process in ways that prohibit the successful extraction of discrepancies between lattice elements in the real machine and their description in the corresponding model. We systematically studied the influence of the placement of BPMs and steerers on the reconstruction of linear optics and corresponding lattice information. The derivative of orbit response coefficients with respect to the quadrupole strengths, the Jacobian, is derived as an analytical formula. This analytical version of the Jacobian is used to further derive the theoretical limitations of fitting linear optics from closed orbits in terms of the placement of BPMs and steerers. It is further demonstrated that when evaluating the Jacobian during the fitting procedure, the analytical version can be used in place of the conventional finite-difference computation. This allows for greatly improved efficiency when computing the Jacobian during each iteration of the fitting procedure. The approach is tested with large-scale simulations and the findings are verified by measurement data taken on SIS18 synchrotron at GSI Helmholtz Centre for Heavy Ion Research. The presented methods are of general nature and can be applied to other accelerator lattices as well. The fitting procedure by using the analytical Jacobian is tested in conjunction with various methods for mitigating quasidegeneracy and the results agree with those obtained by using the conventional Jacobian via finite-difference approximation.

47 OTHER INSTRUMENTATION↗
NASA NTRS1995

On the uniqueness of linear moving-average filters for the solar wind-auroral geomagnetic activity coupling

The relation between the solar wind input to the magetosphere, VB(sub South), and the auroral geomagnetic index AL is modeled with two linear moving-average filtering methods: linear prediction filters and a driven harmonic oscillator in the form of an electric circuit. Although the response of the three-parameter oscillator is simpler than the filter's, the methods yield similar linear timescales and values of the prediction-observation correlation and the prediction Chi(exp 2). Further the filter responses obtained by the two methods are similar in their long-term features. In these aspects the circuit model is equivalent to linear prediction filtering. This poses the question of uniqueness and proper interpretation of detailed features of the filters such as response peaks. Finally, the variation of timescales and filter responses with the AL activity level is discussed.

Vassiliadis, D.↗