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 577 records · Page 32

A componential model of human interaction with graphs: 1. Linear regression modeling

Task analyses served as the basis for developing the Mixed Arithmetic-Perceptual (MA-P) model, which proposes (1) that people interacting with common graphs to answer common questions apply a set of component processes-searching for indicators, encoding the value of indicators, performing arithmetic operations on the values, making spatial comparisons among indicators, and repsonding; and (2) that the type of graph and user's task determine the combination and order of the components applied (i.e., the processing steps). Two experiments investigated the prediction that response time will be linearly related to the number of processing steps according to the MA-P model. Subjects used line graphs, scatter plots, and stacked bar graphs to answer comparison questions and questions requiring arithmetic calculations. A one-parameter version of the model (with equal weights for all components) and a two-parameter version (with different weights for arithmetic and nonarithmetic processes) accounted for 76%-85% of individual subjects' variance in response time and 61%-68% of the variance taken across all subjects. The discussion addresses possible modifications in the MA-P model, alternative models, and design implications from the MA-P model.

Gillan, Douglas J.↗

Aeroelastic Indicial Response Reduced-Order Modeling for Flexible Flight Vehicles

A reduced-order modeling method capable of providing computationally efficient predictions of the nonlinear, unsteady aerodynamics encountered by flexible flight vehicles under- going forced oscillations is presented. Models are developed using indicial response theory, which characterizes a vehicle’s dynamics through identification of time-accurate aerodynamic responses due to step changes in the vehicle-state parameters, e.g., angle-of-attack, pitch rate. A coupled computational fluid dynamics aeroelastic analysis is proposed for identifying step (indicial) responses of flexible vehicles. In this approach, aeroelastic indicial responses are simulated via prescribed rigid body motions, while fluid-structure interactions are captured at the subiterative level through coupling to a linear modal structural solver. A nonlinear extension of indicial response theory is applied through time-dependent linear interpolation of a database of locally linear aeroelastic step responses. Reduced-order models are then created using the mathematical principle of convolution applied to the interpolated aeroelastic indicial responses to predict the time-dependent aerodynamic response of a flexible vehicle to any arbitrary prescribed flight maneuver. The NASA FUN3D computational fluid dynamics solver is utilized for simulating full-order trajectories and indicial response functions. Aerodynamic predictions were generated for the X-56A aircraft undergoing a series of harmonic forced oscillations. The s are shown to provide a practical option for evaluating the unsteady aerodynamics of flexible vehicles using high-fidelity simulations.

Hiller, Brett↗

Sensitivity of the human circadian pacemaker to nocturnal light: melatonin phase resetting and suppression

Ocular exposure to early morning room light can significantly advance the timing of the human circadian pacemaker. The resetting response to such light has a non-linear relationship to illuminance. The dose-response relationship of the human circadian pacemaker to late evening light of dim to moderate intensity has not been well established. Twenty-three healthy young male and female volunteers took part in a 9 day protocol in which a single experimental light exposure6.5 h in duration was given in the early biological night. The effects of the light exposure on the endogenous circadian phase of the melatonin rhythm and the acute effects of the light exposure on plasma melatonin concentration were calculated. We demonstrate that humans are highly responsive to the phase-delaying effects of light during the early biological night and that both the phase resetting response to light and the acute suppressive effects of light on plasma melatonin follow a logistic dose-response curve, as do many circadian responses to light in mammals. Contrary to expectations, we found that half of the maximal phase-delaying response achieved in response to a single episode of evening bright light ( approximately 9000 lux (lx)) can be obtained with just over 1 % of this light (dim room light of approximately 100 lx). The same held true for the acute suppressive effects of light on plasma melatonin concentrations. This indicates that even small changes in ordinary light exposure during the late evening hours can significantly affect both plasma melatonin concentrations and the entrained phase of the human circadian pacemaker.

NASA Discipline Regulatory Physiology↗

Gender-related Changes in Dorsal Hand and Foot Vein Function Following 60 Days of Head Down Bed Rest

It is well known that female astronauts are more likely to experience post-flight orthostatic hypotension and presyncope compared to male astronauts. It has been suggested that the disproportionally higher incidence of presyncope (83% of female vs. 20% male crewmembers) may be due to sex-related differences in vascular function between the upper and lower limbs. However, much of this evidence is specific to changes in resistance vessels. Given that more than 70% of the circulating blood volume resides in compliance vessels, it is conceivable that even small changes in venous function may contribute to post-flight orthostatic hypotension. In spite of this, little is currently known regarding the influence of microgravity exposure on venous function between males and females. PURPOSE: To determine the influence of 60 days of HDBR on dorsal foot and hand vein function between healthy males (M) and females (F). METHODS: Using 2-D ultrasound, dorsal hand and foot vein diameter responses to intravenous infusions phenylephrine (PE), acetylcholine (ACh), and nitroglycerine (NTG) were determined in 26 adults; 10 females (age:37 +/- 2 yr ) and 16 males (age:34 +/- 2 yr ). Changes in venous function were calculated as the difference between diameter at baseline and following each venoactive drug. Differences in venous function between limb and sexes across HDBR were determined using mixed-effects linear regression. RESULTS: In response to 60 days of HDBR, the change in venousconstrictor response to PE in the dorsal hand veins was not significantly different between M and F. Interestingly, the change in constrictor response in the dorsal foot veins (compared to pre HDBR) was approximately 30% greater in the F, whereas the constrictor response was approximately 45% less in the M (p=0.026). HDBR had no influence on the change in dilator response to ACh, or NTG between M and F and between vascular beds. CONCLUSION: These results demonstrate that 60 days of HDBR contributes to sex-related differences in venous function. Importantly, the changes are specific to dorsal foot veins and do not appear to be related to changes in endothelial or smooth muscle function. Previous data from our laboratory indicate that females have an impaired sympathetic response following HDBR. However, the results from the current investigation imply that F maintain the ability to respond to sympathetic stimulation but may not be able to adequately activate a proper response following HDBR.

Westby, Christian M.↗

Nonlinear dynamic analysis of sandwich panels

Two analytical techniques applicable to large deflection dynamic response calculations for pressure loaded composite sandwich panels are demonstrated. One technique utilizes finite element modeling with a single equivalent layer representing the face sheets and core. The other technique utilizes the modal analysis computer code DEPROP which was recently modified to include transverse shear deformation in a core layer. The example problem consists of a simply supported rectangular sandwich panel. Included are comparisons of linear and nonlinear static response calculations, in addition to dynamic response calculations.

Lush, A. M.↗

Three-dimensional organization of otolith-ocular reflexes in rhesus monkeys. III. Responses To translation

The three-dimensional (3-D) properties of the translational vestibulo-ocular reflexes (translational VORs) during lateral and fore-aft oscillations in complete darkness were studied in rhesus monkeys at frequencies between 0.16 and 25 Hz. In addition, constant velocity off-vertical axis rotations extended the frequency range to 0.02 Hz. During lateral motion, horizontal responses were in phase with linear velocity in the frequency range of 2-10 Hz. At both lower and higher frequencies, phase lags were introduced. Torsional response phase changed more than 180 degrees in the tested frequency range such that torsional eye movements, which could be regarded as compensatory to "an apparent roll tilt" at the lowest frequencies, became anticompensatory at all frequencies above approximately 1 Hz. These results suggest two functionally different frequency bandwidths for the translational VORs. In the low-frequency spectrum (<<0.5 Hz), horizontal responses compensatory to translation are small and high-pass-filtered whereas torsional response sensitivity is relatively frequency independent. At higher frequencies however, both horizontal and torsional response sensitivity and phase exhibit a similar frequency dependence, suggesting a common role during head translation. During up-down motion, vertical responses were in phase with translational velocity at 3-5 Hz but phase leads progressively increased for lower frequencies (>90 degrees at frequencies <0.2 Hz). No consistent dependence on static head orientation was observed for the vertical response components during up-down motion and the horizontal and torsional response components during lateral translation. The frequency response characteristics of the translational VORs were fitted by "periphery/brain stem" functions that related the linear acceleration input, transduced by primary otolith afferents, to the velocity signals providing the input to the velocity-to-position neural integrator and the oculomotor plant. The lowest-order, best-fit periphery/brain stem model that approximated the frequency dependence of the data consisted of a second order transfer function with two alternating poles (at 0.4 and 7.2 Hz) and zeros (at 0.035 and 3.4 Hz). In addition to clearly differentiator dynamics at low frequencies (less than approximately 0.5 Hz), there was no frequency bandwidth where the periphery/brain stem function could be approximated by an integrator, as previously suggested. In this scheme, the oculomotor plant dynamics are assumed to perform the necessary high-frequency integration as required by the reflex. The detailed frequency dependence of the data could only be precisely described by higher order functions with nonminimum phase characteristics that preclude simple filtering of afferent inputs and might be suggestive of distributed spatiotemporal processing of otolith signals in the translational VORs.

Non-NASA Center↗

Estimating Single-Trial Responses in EEG

Accurate characterization of single-trial field potential responses is critical from a number of perspectives. For example, it allows differentiation of an evoked response from ongoing EEG. We previously developed the multiple component Event Related Potential (mcERP) algorithm to improve resolution of the single-trial evoked response. The mcERP model states that multiple components, each specified by a stereotypic waveform varying in latency and amplitude from trial to trial, comprise the evoked response. Application of the mcERP algorithm to simulated data with three independent, synthetic components has shown that the model is capable of separating these components and estimating their variability. Application of the model to single trial, visual evoked potentials recorded simultaneously from all V1 laminae in an awake, fixating macaque yielded local and far-field components. Certain local components estimated by the model were distributed in both granular and supragranular laminae. This suggests a linear coupling between the responses of thalamo-recipient neuronal ensembles and subsequent responses of supragranular neuronal ensembles, as predicted by the feedforward anatomy of V1. Our results indicate that the mcERP algorithm provides a valid estimation of single-trial responses. This will enable analyses that depend on trial-to-trial variations and those that require separation of the evoked response from background EEG rhythms

Shah, A. S.↗

Recursive Inversion Of Externally Defined Linear Systems

Technical memorandum discusses mathematical technique described in "Recursive Inversion by Finite-Impulse-Response Filters" (ARC-12247). Technique is recursive algorithm yielding finite-impulse-response approximation of unknown single-input/single-output, causal, time-invariant, linear, real system, response of which is sequence of impulses. Useful in such diverse applications as medical diagnoses, identification of military targets, geophysical exploration, and nondestructive testing.

Bach, Ralph E., Jr.↗

ASTROP2-LE: A Mistuned Aeroelastic Analysis System Based on a Two Dimensional Linearized Euler Solver

An aeroelastic analysis system for flutter and forced response analysis of turbomachines based on a two-dimensional linearized unsteady Euler solver has been developed. The ASTROP2 code, an aeroelastic stability analysis program for turbomachinery, was used as a basis for this development. The ASTROP2 code uses strip theory to couple a two dimensional aerodynamic model with a three dimensional structural model. The code was modified to include forced response capability. The formulation was also modified to include aeroelastic analysis with mistuning. A linearized unsteady Euler solver, LINFLX2D is added to model the unsteady aerodynamics in ASTROP2. By calculating the unsteady aerodynamic loads using LINFLX2D, it is possible to include the effects of transonic flow on flutter and forced response in the analysis. The stability is inferred from an eigenvalue analysis. The revised code, ASTROP2-LE for ASTROP2 code using Linearized Euler aerodynamics, is validated by comparing the predictions with those obtained using linear unsteady aerodynamic solutions.

Reddy, T. S. R.↗

Non-Additivity of the Midlatitude Circulation Response to Regional Arctic Temperature Anomalies: The Role of the Stratosphere

Previous studies have documented the impact of the Arctic sea ice loss and associated warming on the midlatitude weather and climate, especially the influence of sea ice retreat over the Barents-Kara Sea on the North Atlantic and Europe regions. However, less attention has been given to other geographical locations over the Arctic, and to the linear additivity of the circulation response to regional Arctic sea ice loss and temperature anomalies. Using a simplified dry dynamical core model, we demonstrate that responses to regional Arctic temperature anomalies over the Barents-Kara Sea, Baffin Bay-Davis Strait-Labrador Sea, and East Siberia-Chukchi Sea, separately, cause similar equatorward shift of the tropospheric jet, but different stratospheric polar vortex responses. Furthermore, responses to regional Arctic temperature anomalies are not linearly additive, and the residual resembles a positive Northern Annular Mode-like structure. Additional targeted experiments highlight the stratospheric influence in the non-additivity of the midlatitude tropospheric response.

Bithi De↗

The Response of the Ozone Layer Under Abrupt 4xCO2 in CMIP6

Previous studies indicate a possible role of stratospheric ozone chemistry feedbacks in the climate response to 4xCO2, either via a reduction in equilibrium climate sensitivity (ECS) (Nowack et al., 2015), or via changes in the tropospheric circulation (Chiodo and Polvani, 2017). However, these effects are subject to uncertainty. Part of the uncertainty may stem from the dependency of the feedback on the pattern of the ozone response, as the radiative efficiency of ozone largely depends on its vertical distribution (Lacis et al., 1990). Here, an analysis is presented of the ozone layer response to 4xCO2 in chemistry–climate models (CCMs) which participated to CMIP intercomparisons. In a previous study using CMIP5 models, it has been shown that under 4xCO2, ozone decreases in the tropical lower stratosphere, and increases over the high latitudes and throughout the upper stratosphere. It was also found that a substantial portion of the spread in the tropical column ozone is tied to inter-model spread in tropical upwelling, which is in turn partly tied to ECS (Chiodo et al., 2018). Here, we revisit this connection using4xCO2 data from CMIP6, thereby exploiting the larger number of CCMs than in CMIP5. In addition, we explore the linearity of the ozone response, by complementing the analysis with simulations using transient CO2 forcing (1pc/year). We show that the pattern of the ozone response in is similar to CMIP5. In some models (e.g. WACCM), we find larger ozone responses in CMIP6 than in CMIP5, partly because of the larger ECS and thus larger upwelling response in the tropical pipe. In this presentation, we will discuss the relationship between radiative forcing, transport and ozone, as well as further implications for CMIP6 models.

Ozone layer↗

The Response of the Ozone Layer Under Abrupt 4xCO2 in CMIP6

Previous studies indicate a possible role of stratospheric ozone chemistry feedbacks in the climate response to 4xCO2, either via a reduction in equilibrium climate sensitivity (ECS) (Nowack et al.,2015), or via changes in the tropospheric circulation (Chiodo and Polvani, 2017). However, these effects are subject to uncertainty. Part of the uncertainty may stem from the dependency of the feedback on the pattern of the ozone response, as the radiative efficiency of ozone largely depends on its vertical distribution (Lacis et al., 1990). Here, an analysis is presented of the ozone layer response to 4xCO2 in chemistry–climate models (CCMs) which participated to CMIP intercomparisons. In a previous study using CMIP5 models, it has been shown that under 4xCO2, ozone decreases in the tropical lower stratosphere, and increases over the high latitudes and throughout the upper stratosphere. It was also found that a substantial portion of the spread in the tropical column ozone is tied to inter-model spread in tropical upwelling, which is in turn tied to ECS (Chiodo et al., 2018). Here, we revisit this connection using 4xCO2 data from CMIP6, thereby exploiting the larger number of CCMs than in CMIP5. In addition, we explore the linearity of the ozone response, by complementing the analysis with simulations using lower CO 2 forcing levels (2xCO2). We show that the pattern of the ozone response in is similar to CMIP5. In some models (e.g. WACCM), we find larger ozone responses in CMIP6 than in CMIP5, partly because of the larger ECS and thus larger upwelling response in the tropical pipe. In this presentation, we will discuss the relationship between radiative forcing, transport and ozone, as well as further implications for CMIP6 models.

Ozone layer↗

Development of a GCR Event-based Risk Model

A goal at NASA is to develop event-based systems biology models of space radiation risks that will replace the current dose-based empirical models. Complex and varied biochemical signaling processes transmit the initial DNA and oxidative damage from space radiation into cellular and tissue responses. Mis-repaired damage or aberrant signals can lead to genomic instability, persistent oxidative stress or inflammation, which are causative of cancer and CNS risks. Protective signaling through adaptive responses or cell repopulation is also possible. We are developing a computational simulation approach to galactic cosmic ray (GCR) effects that is based on biological events rather than average quantities such as dose, fluence, or dose equivalent. The goal of the GCR Event-based Risk Model (GERMcode) is to provide a simulation tool to describe and integrate physical and biological events into stochastic models of space radiation risks. We used the quantum multiple scattering model of heavy ion fragmentation (QMSFRG) and well known energy loss processes to develop a stochastic Monte-Carlo based model of GCR transport in spacecraft shielding and tissue. We validated the accuracy of the model by comparing to physical data from the NASA Space Radiation Laboratory (NSRL). Our simulation approach allows us to time-tag each GCR proton or heavy ion interaction in tissue including correlated secondary ions often of high multiplicity. Conventional space radiation risk assessment employs average quantities, and assumes linearity and additivity of responses over the complete range of GCR charge and energies. To investigate possible deviations from these assumptions, we studied several biological response pathway models of varying induction and relaxation times including the ATM, TGF -Smad, and WNT signaling pathways. We then considered small volumes of interacting cells and the time-dependent biophysical events that the GCR would produce within these tissue volumes to estimate how GCR event rates mapped to biological signaling induction and relaxation times. We considered several hypotheses related to signaling and cancer risk, and then performed simulations for conditions where aberrant or adaptive signaling would occur on long-duration space mission. Our results do not support the conventional assumptions of dose, linearity and additivity. A discussion on how event-based systems biology models, which focus on biological signaling as the mechanism to propagate damage or adaptation, can be further developed for cancer and CNS space radiation risk projections is given.

Cucinotta, Francis A.↗

A geometrically nonlinear analysis of interlaminar stresses in unsymmetrically laminated plates subjected to uniform thermal loading

An analytical study of interlaminar stresses in unsymmetrically laminated plates is presented. The study examines the linear elastic, large deflection response of square laminated composite plates subjected to uniform thermal loading. Both cross-ply and angle-ply, symmetric and unsymmetric, laminates are examined to evaluate the effects of mismatch between adjacent layers in elastic constants and coefficients of thermal expansion. A geometrically nonlinear kinematic description is used to predict the large out-of-plane (transverse) deflections. The nonlinear, three-dimensional boundary value problems are formulated from elasticity theory and approximate solutions are determined using the finite element method. A global/local analysis procedure is used to obtain improved free edge stress predictions. For the laminates and loading conditions considered, the results indicate that the out-of-plane deflections of the unsymmetric laminates reduce interlaminar shear stresses. These deflections also reduce interlaminar normal stresses in some laminates and increase these stresses for others. The results indicate that as the out-of-plane deflections become large, the differences in interlaminar stress predictions between linear and nonlinear theory can become quite large.

Norwood, D. S.↗

On identifying transfer functions and state equations for linear systems.

Two methods are established for identifying constant-coefficient, C to the 2n power type of noise-free linear systems if the time response data of the input-output or of all states are known. 2n response data are required to identify an nth-order transfer function or state equation for an unknown linear system. The order of the unknown system can be identified by checking a sequence of determinants. The Z transform and its inversion are mainly used.

Shieh, L. S.↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗

Linearized flexibility models in multibody dynamics and control

Simulation of structural response of multi-flexible-body systems by linearized flexible motion combined with nonlinear rigid motion is discussed. Advantages and applicability of such an approach for accurate simulation with greatly reduced computational costs and turnaround times are described, restricting attention to the control design environment. Requirements for updating the linearized flexibility model to track large angular motions are discussed. Validation of such an approach by comparison with other existing codes is included. Application to a flexible robot manipulator system is described.

Cimino, William W.↗