Comparison of several methods for obtaining the time response of linear systems to either a unit impulse or arbitrary input from frequency-response data
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Methods are presented that use general correlative time-response input and output data for a linear system to determine the frequency-response function of that system. These methods give an exact description of any linear system for which such transient data are available. Examples are shown of application of a method to both an underdamped and a critically damped exact second-order system, and to an exact first-order system with and without dead time. Experimental data for a turbine-propeller engine showing the response of engine speed to change in propeller-blade angle are presented and analyzed.
Data-reduction methods using general-purpose analog computer equipment and compatible testing techniques for determining the frequency response of linear physical systems are examined. The techniques considered may be classed as steady state or transient depending on the method of excitation. The relative merits of periodic, slow sweep, and transient (rapid sweep) forcing functions are discussed and applications are given that relate to dynamic-response tests of aeroelastic systems. Two frequency-sweep-input methods are considered in detail. one case the sweep rate is sufficiently slow that the response is approximately the same as that for steady-state conditions. With this input the frequency response can be evaluated and displayed in real time while the test is in progress. Errors due to treating sweep data as steady state can be eliminated, when desired, by reanalyzing tape-recorded time histories of the input and output as transient rather than as periodic data. In the second method the frequency-response function is deter- mined from the system's transient response to a very rapid sweep input. The purpose of frequency sweep in this case is to provide sufficient harmonic content in the input to overcome noise while keeping the test time as short as possible. tests and limited flight-test data presented herein, it appears that a transient-type rapid-sweep forcing function offers a considerable saving in test time while preserving the accuracy possible with steady-state sinusoidal inputs.
Accurate approximations of the change of a system's output and its statistics with respect to the input are highly desired in computational dynamics. Ruelle's linear response theory provides breakthrough mathematical machinery for computing the linear response of chaotic dynamical systems. In this paper, we propose an algorithm for sensitivity analysis of discrete chaos with an arbitrary number of positive Lyapunov exponents. We combine the concept of perturbation space-splitting, which regularizes Ruelle's original expression, together with measure-based parameterization of the expanding subspace. We use these tools to rigorously derive trajectory-following recursive relations that converge exponentially fast, and construct a memory-efficient Monte Carlo scheme for derivatives of the output statistics. Thanks to the regularization and lack of simplifying assumptions on the system's behavior, our method is immune to the common problems of other popular methods such as the exploding tangent solutions and unphysical shadowing directions. Here, we provide a ready-to-use algorithm, analyze its complexity, and demonstrate several numerical examples of sensitivity computation using physically-inspired low-dimensional systems.
Recent Enhanced Low Dose Rate Sensitivity (ELDRS) investigations carried out by RLP Research, Arizona State University (ASU) and Jet Propulsion Laboratory (JPL) have shown significant differences in the degradation of bipolar micro-circuits with total dose in the presence of molecular hydrogen (H2) in packages. This has a significant impact on radiation hardness assurance and opens up opportunities to improve device performance. The general objectives of this program are to: 1. Determine the extent to which hydrogen contamination affects the total dose and dose rate response of linear bipolar circuits; 2. Develop a model that will enable the prediction of high dose rate (HDR) and low dose rate (LDR) response asymptotes and transition dose rates as a function of total dose, temperature, pressurized hydrogen, defect precursors, and other process dependent variables; 3. Explore the possibility of an accelerated hardness assurance method and possible hardening approaches; and 4. Extend the work to other technologies that have total dose response affected by hydrogen contamination.In this document, we xperimentally demonstrate with test transistors and circuits that hydrogen is correlated with ELDRS in bipolar linear circuits. We show that the amount of hydrogen determines: 1) the total dose response versus dose rate and 2) the transition dose rates between the high and low dose rate responses. The experimental results are supported with a steady state drift/diffusion analytical model as well as modeling calculations using COMSOL Multiphysics.
Vestibular evoked potentials (VsEPs) were measured in normal mice and in mice homozygous for the head tilt mutation (het/het, abbr. het). The het mice lack otoconia, the inertial mass critical for natural stimulation of inner ear gravity receptors. Our findings demonstrate that vestibular neural responses to pulsed linear acceleration are absent in het mice. The results: (1) confirm that adequate sensory stimuli fail to activate gravity receptors in the het model; and (2) serve as definitive evidence that far-field vestibular responses to pulsed linear acceleration depend critically on otolith end organs. The C57BL/6JEi-het mouse may be an excellent model of gravity receptor sensory deprivation.
To support the reliable and resilient operation of modular reactors and microreactors, anticipatory control strategies have been proposed for achieving faster-than-real-time predictions and decision-making capabilities in anticipation of potential anomalies, including setpoint changes and cyber incidents. Here this work presents how anticipatory control strategies can be implemented via model predictive control (MPC) of a single heat pipe’s temperature. Considering the uncertainty in developing and applying MPC, this work evaluates MPC performance given three different model forms: a linear response surface model, an artificial neural network (ANN), and an autoregressive model with exogenous input (ARX). This work also evaluates the impacts of different input biases and variance on MPC performance in order to account for potential sensor reading variations due to cyber incidents. We observe that nonparametric models such as the ANN and ARX result in more fluctuated control actions compared to the MPC applied to the linear response surface model. However, when the cyber incidents are of a large magnitude, the linear response surface model produces smaller feasible regions than the nonparametric models under identical constraints.
Recent enhanced low dose rate sensitivity (ELDRS) investigations carried out by RLP Research, Crane, Arizona State University (ASU), and Jet Propulsion Laboratory (JPL) have shown significant differences in the degradation of bipolar microcircuits with total dose in the presence of molecular hydrogen (H2) in packages. This has a significant impact on radiation hardness assurance and opens up opportunity to improve device performance. The objectives of this program are 1) to investigate and confirm the causal relationship between packaging recipes, hydrogen contamination, and total dose response of linear bipolar microcircuits; and 2) to develop a guideline that will take into account these effects at the radiation level for future NASA space missions. This program is geared to benefit all future NASA space missions using bipolar or BiCMOS linear devices.
This paper is a follow-on to recent work by the authors in which the response and high-cycle fatigue of a nonlinear structure subject to non-Gaussian loadings was found to vary markedly depending on the nature of the loading. There it was found that a non-Gaussian loading having a steady rate of short-duration, high-excursion peaks produced essentially the same response as would have been incurred by a Gaussian loading. In contrast, a non-Gaussian loading having the same kurtosis, but with bursts of high-excursion peaks was found to elicit a much greater response. This work is meant to answer the question of when consideration of a loading probability distribution other than Gaussian is important. The approach entailed nonlinear numerical simulation of a beam structure under Gaussian and non-Gaussian random excitations. Whether the structure responded in a Gaussian or non-Gaussian manner was determined by adherence to, or violations of, the Central Limit Theorem. Over a practical range of damping, it was found that the linear response to a non-Gaussian loading was Gaussian when the period of the system impulse response is much greater than the rate of peaks in the loading. Lower damping reduced the kurtosis, but only when the linear response was non-Gaussian. In the nonlinear regime, the response was found to be non-Gaussian for all loadings. The effect of a spring-hardening type of nonlinearity was found to limit extreme values and thereby lower the kurtosis relative to the linear response regime. In this case, lower damping gave rise to greater nonlinearity, resulting in lower kurtosis than a higher level of damping.
Van der Waals interactions are ubiquitous and they play an important role for the stability of materials. Current understanding of this type of coupling is based on linear response theory, while optical nonlinearities are rarely considered in this context. Many materials, however, exhibit strong optical nonlinear response, which prompts further evaluation of dispersive forces beyond linear response. Here we present a approach that takes into account linear and nonlinear properties of all dipolar nanoparticles in a given system. This method is based on a Hamiltonian for nonlinear dipoles, which we apply in different systems uncovering a complex interplay of distance, anisotropy, polarizabilities, and hyperpolarizabilities in the vdW energy. This investigation broadens our basic understanding of dispersive interactions, especially in the context of nonlinear materials. Published by the American Physical Society 2024
1. The dynamic properties of otolith-ocular reflexes elicited by sinusoidal linear acceleration along the three cardinal head axes were studied during off-vertical axis rotations in rhesus monkeys. As the head rotates in space at constant velocity about an off-vertical axis, otolith-ocular reflexes are elicited in response to the sinusoidally varying linear acceleration (gravity) components along the interaural, nasooccipital, or vertical head axis. Because the frequency of these sinusoidal stimuli is proportional to the velocity of rotation, rotation at low and moderately fast speeds allows the study of the mid-and low-frequency dynamics of these otolith-ocular reflexes. 2. Animals were rotated in complete darkness in the yaw, pitch, and roll planes at velocities ranging between 7.4 and 184 degrees/s. Accordingly, otolith-ocular reflexes (manifested as sinusoidal modulations in eye position and/or slow-phase eye velocity) were quantitatively studied for stimulus frequencies ranging between 0.02 and 0.51 Hz. During yaw and roll rotation, torsional, vertical, and horizontal slow-phase eye velocity was sinusoidally modulated as a function of head position. The amplitudes of these responses were symmetric for rotations in opposite directions. In contrast, mainly vertical slow-phase eye velocity was modulated during pitch rotation. This modulation was asymmetric for rotations in opposite direction. 3. Each of these response components in a given rotation plane could be associated with an otolith-ocular response vector whose sensitivity, temporal phase, and spatial orientation were estimated on the basis of the amplitude and phase of sinusoidal modulations during both directions of rotation. Based on this analysis, which was performed either for slow-phase eye velocity alone or for total eye excursion (including both slow and fast eye movements), two distinct response patterns were observed: 1) response vectors with pronounced dynamics and spatial/temporal properties that could be characterized as the low-frequency range of "translational" otolith-ocular reflexes; and 2) response vectors associated with an eye position modulation in phase with head position ("tilt" otolith-ocular reflexes). 4. The responses associated with two otolith-ocular vectors with pronounced dynamics consisted of horizontal eye movements evoked as a function of gravity along the interaural axis and vertical eye movements elicited as a function of gravity along the vertical head axis. Both responses were characterized by a slow-phase eye velocity sensitivity that increased three- to five-fold and large phase changes of approximately 100-180 degrees between 0.02 and 0.51 Hz. These dynamic properties could suggest nontraditional temporal processing in utriculoocular and sacculoocular pathways, possibly involving spatiotemporal otolith-ocular interactions. 5. The two otolith-ocular vectors associated with eye position responses in phase with head position (tilt otolith-ocular reflexes) consisted of torsional eye movements in response to gravity along the interaural axis, and vertical eye movements in response to gravity along the nasooccipital head axis. These otolith-ocular responses did not result from an otolithic effect on slow eye movements alone. Particularly at high frequencies (i.e., high speed rotations), saccades were responsible for most of the modulation of torsional and vertical eye position, which was relatively large (on average +/- 8-10 degrees/g) and remained independent of frequency. Such reflex dynamics can be simulated by a direct coupling of primary otolith afferent inputs to the oculomotor plant. (ABSTRACT TRUNCATED).
Imaging detecting systems based on microchannel plates (MCPs) are the most common for low energy plasma measurements for both space borne and ground applications. One of the key parameters of these detection systems is the dynamic range of the MCP's response to the input fluxes of charged particles. For most applications the dynamic range of the linear response should be as wide as possible. This is especially true for the Dual Electron Spectrometers (DESs) of the Fast Plasma Investigation (FPI) on NASA's Magnetospheric MultiScale (MMS) mission because a wide range of input fluxes are expected. To make use of the full available dynamic range, it is important to understand the MCP response behavior beyond the linear regime where the MCPs start to saturate. We have performed extensive studies of this during the characterization and calibration of the DES instruments and have identified several saturation effects of the detection system. The MCP itself exhibits saturation when the channels lack the ability to replenish charge sufficiently rapidly. It is found and will be shown that the ground system can significantly impact the correct measurement of this effect. As the MCP starts to saturate, the resulting pulse height distribution (PHD) changes shape and location (with less pulse height values), which leads to truncation of the PHD by the threshold set on the detection system discriminator. Finally, the detection system pulse amplifier exhibits saturation as the input flux drives pulse rates greater than its linear response speed. All of these effects effectively change the dead time of the overall detection system and as a result can affect the quality and interpretation of the flight data. We present results of detection system saturation effects and their interaction with special emphasis on the MCP related effects.
Vestibular responses of single otoliths of frog to linear acceleration and short periods of weightlessness during parabolic flight
Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Delta, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. Here, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Delta. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Delta variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfven wave resonances as rotation slows down. These are important features that can alter the nonaxisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.
Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Δ, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. In this work, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Δ. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Δ variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfvén wave resonances as rotation slows down. These are important features that can alter the non-axisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.
Transient linear acceleration stimuli have been shown to elicit eighth nerve vestibular compound action potentials in birds and mammals. The present study was undertaken to better define the nature of the adequate stimulus for neurons generating the response in the chicken (Gallus domesticus). In particular, the study evaluated the question of whether the neurons studied are most sensitive to the maximum level of linear acceleration achieved or to the rate of change in acceleration (da/dt, or jerk). To do this, vestibular response thresholds were measured as a function of stimulus onset slope. Traditional computer signal averaging was used to record responses to pulsed linear acceleration stimuli. Stimulus onset slope was systematically varied. Acceleration thresholds decreased with increasing stimulus onset slope (decreasing stimulus rise time). When stimuli were expressed in units of jerk (g/ms), thresholds were virtually constant for all stimulus rise times. Moreover, stimuli having identical jerk magnitudes but widely varying peak acceleration levels produced virtually identical responses. Vestibular response thresholds, latencies and amplitudes appear to be determined strictly by stimulus jerk magnitudes. Stimulus attributes such as peak acceleration or rise time alone do not provide sufficient information to predict response parameter quantities. Indeed, the major response parameters were shown to be virtually independent of peak acceleration levels or rise time when these stimulus features were isolated and considered separately. It is concluded that the neurons generating short latency vestibular evoked potentials do so as "jerk encoders" in the chicken. Primary afferents classified as "irregular", and which traditionally fall into the broad category of "dynamic" or "phasic" neurons, would seem to be the most likely candidates for the neural generators of short latency vestibular compound action potentials.
The spontaneous activity of 40 otolith afferents and 44 canal afferents was recorded in 4 alert, intact squirrel monkeys. Polarization vectors and response properties of otolith afferents were determined during static re-orientations relative to gravity and during Earth-horizontal, sinusoidal, linear oscillations. Canal afferents were tested for sensitivity to linear accelerations. For regular otolith afferents, a significant correlation between upright discharge rate and sensitivity to dynamic acceleration in the horizontal plane was observed. This correlation was not present in irregular units. The sensitivity of otolith afferents to both static tilts and dynamic linear acceleration was much greater in irregularly discharging units than in regularly discharging units. The spontaneous activity and static and dynamic response properties of regularly discharging otolith afferents were similar to those reported in barbiturate-anesthetized squirrel monkeys. Irregular afferents also had similar dynamic response properties when compared to anesthetized monkeys. However, this sample of irregular afferents in alert animals had higher resting discharge rates and greater sensitivity to static tilts. The majority of otolith polarization vectors were oriented near the horizontal in the plane of the utricular maculae; however, directions of maximum sensitivity were different during dynamic and static testing. Canal afferents were not sensitive to static tilts or linear oscillations of the head.
Linear analysis of damped flexibly mounted rolling element bearings influence on dynamic rotor unbalance response, including optimum support characteristics