Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “non linear dynamics”

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 433 records · Page 24

Digital optical recorder-reproducer system

A mass archival optical recording and reproduction system includes a recording light source such as a laser beam focussed and directed upon an acousto-optic linear modulator array (or page composer) that receives parallel blocks of data converted from a serial stream of digital data to be stored. The page composer imparts to the laser beam modulation representative of a plurality of parallel channels of data and through focussing optics downstream of the page composer parallel arrays of optical spots are recorded upon a suitable recording medium such as a photographic film floppy disc. The recording medium may be substantially frictionlessly and stably positioned for recording at a record/read station by an air-bearing platen arrangement which is preferably thermodynamically non-throttling so that the recording film may be positioned in the path of the information-carrying light beam in a static or dynamic mode. During readout, the page composer is bypassed and a readout light beam is focussed directly upon the recording medium containing an array of previously recorded digital spots, a sync bit, data positioning bits, and a tracking band. The readout beam which has been directed through the recording medium is then imaged upon a photodetector array, the output of which may be coupled to suitable electronic processing circuitry, such as a digital multiplexer, whereby the parallel spot array is converted back into the original serial data stream.

Reddersen, Brad R.↗

Sensitivity analysis and approximation methods for general eigenvalue problems

Optimization of dynamic systems involving complex non-hermitian matrices is often computationally expensive. Major contributors to the computational expense are the sensitivity analysis and reanalysis of a modified design. The present work seeks to alleviate this computational burden by identifying efficient sensitivity analysis and approximate reanalysis methods. For the algebraic eigenvalue problem involving non-hermitian matrices, algorithms for sensitivity analysis and approximate reanalysis are classified, compared and evaluated for efficiency and accuracy. Proper eigenvector normalization is discussed. An improved method for calculating derivatives of eigenvectors is proposed based on a more rational normalization condition and taking advantage of matrix sparsity. Important numerical aspects of this method are also discussed. To alleviate the problem of reanalysis, various approximation methods for eigenvalues are proposed and evaluated. Linear and quadratic approximations are based directly on the Taylor series. Several approximation methods are developed based on the generalized Rayleigh quotient for the eigenvalue problem. Approximation methods based on trace theorem give high accuracy without needing any derivatives. Operation counts for the computation of the approximations are given. General recommendations are made for the selection of appropriate approximation technique as a function of the matrix size, number of design variables, number of eigenvalues of interest and the number of design points at which approximation is sought.

Murthy, D. V.↗

Response properties of pigeon otolith afferents to linear acceleration

In the present study, the sensitivity to sinusoidal linear accelerations in the plane of the utricular macula was tested in afferents. The head orientation relative to the translation axis was varied in order to determine the head position that elicited the maximal and minimal responses for each afferent. The response gain and phase values obtained to 0.5-Hz and 2-Hz linear acceleration stimuli were then plotted as a function of head orientation and a modified cosine function was fit to the data. From the best-fit cosine function, the predicted head orientations that would produce the maximal and minimal response gains were estimated. The estimated maximum response gains to linear acceleration in the utricular plane for the afferents varied between 75 and 1420 spikes s-1 g-1. The mean maximal gains for all afferents to 0.5-Hz and 2-Hz sinusoidal linear acceleration stimuli were 282 and 367 spikes s-1 g-1, respectively. The minimal response gains were essentially zero for most units. The response phases always led linear acceleration and remained constant for each afferent, regardless of head orientation. These response characteristics indicate that otolith afferents are cosine tuned and behave as one-dimensional linear accelerometers. The directions of maximal sensitivity to linear acceleration for the afferents varied throughout the plane of the utricle; however, most vectors were directed out of the opposite ear near the interaural axis. The response dynamics of the afferents were tested using stimulus frequencies ranging between 0.25 Hz and 10 Hz (0.1 g peak acceleration). Across stimulus frequencies, most afferents had increasing gains and constant phase values. These dynamic properties for individual afferents were fit with a simple transfer function that included three parameters: a mechanical time constant, a gain constant, and a fractional order distributed adaptation operator.

Non-NASA Center↗

Response of the QBO to Abrupt Increases in CO2 Using Three Atmospheric Chemistry Configurations

Long-term projections of the Quasi-Biennial Oscillation (QBO) remain highly uncertain. This is partly due to the paucity of models which are able to properly simulate that phenomenon. Only 5 of the 47 CMIP5 models are capable of spontaneously generating a realistic QBO (Butchart et al., 2018), and even those models exhibit large biases in key QBO characteristics (e.g. amplitude, period, vertical extent) when compared with observations. Furthermore, only 1 of these 5 employed interactive atmospheric chemistry, which is known to modulate QBO dynamics. We here investigate the QBO response to increased greenhouse gases using the NASA Goddard Institute for Space Studies Middle Atmosphere Model E2.2. Compared to lower vertical resolution versions of Model E, version 2.2 has a higher model top (0.002 hPa), and additional interactive non-orographic gravity wave drag sources from convection and shear, which produce a sufficiently realistic QBO, thus rendering it suitable for use in climate change studies. Steady-state responses to doubled and quadrupled CO2 from a pre-industrial control are analyzed, as well as the transient response to a 1% per year CO2 increase. In addition, we systematically explore the impact of interactive chemistry in modulating the QBO response to increased CO2 by contrasting interactive, prescribed, and linearized ozone chemistry configurations of the model. Overall, in response to increase CO2 concentrations the QBO is seen to increase in frequency and weaken in amplitude, consistent with previous results, but the memory of the tropical stratosphere may complicate assessments of trends in chemistry and surface impacts. We also discuss implications for the trade-off between ensemble size and the complexity of the chemistry scheme in the model.

Quasi-Biennial Oscillation (QBO)↗

Measurement of the linear viscoelastic behavior of antimisting kerosene

Measurements of dynamic viscoelastic properties in very small oscillating shear deformations was made on solutions of a jet fuel, Jet A, containing an antimisting polymeric additive, FM-9. A few measurements were also made on solutions of FM-9 in a mixed solvent of mineral oil, Tetralin, and 0-terphenyl. Two samples of FM-9 had approximate number-average molecular weights of 12,000,000 and 8,100,000 as deduced from analysis of the measurements. The ranges of variables were 2.42 to 4.03 g/1 in concentration (0.3 to 0.5% by weight), 1 to 35 in temperature, 1.3 to 9.4 cp in solvent viscosity, and 103 to 6100 Hz in frequency. Measurements in the Jet A solvent were made both with and without a modifying carrier. The results were compared with the Zimm theory and the viscoelastic behavior was found to resemble rather closely that of ordinary non-polar polymers in theta solvents. The relation of the results to the antithixotropic behavior of such solutions a high shear rates is discussed in terms of intramolecular and intermolecular interactions.

Ferry, J. D.↗

Probing transmembrane mechanical coupling and cytomechanics using magnetic twisting cytometry

We recently developed a magnetic twisting cytometry technique that allows us to apply controlled mechanical stresses to specific cell surface receptors using ligand-coated ferromagnetic microbeads and to simultaneously measure the mechanical response in living cells. Using this technique, we have previously shown the following: (i) beta 1 integrin receptors mediate mechanical force transfer across the cell surface and to the cytoskeleton, whereas other transmembrane receptors (e.g., scavenger receptors) do not; (ii) cytoskeletal stiffness increases in direct proportion to the level of stress applied to integrins; and (iii) the slope of this linear stiffening response differs depending on the shape of the cell. We now show that different integrins (beta 1, alpha V beta 3, alpha V, alpha 5, alpha 2) and other transmembrane receptors (scavenger receptor, platelet endothelial cell adhesion molecule) differ in their ability to mediate force transfer across the cell surface. In addition, the linear stiffening behavior previously observed in endothelial cells was found to be shared by other cell types. Finally, we demonstrate that dynamic changes in cell shape that occur during both cell spreading and retraction are accompanied by coordinate changes in cytoskeletal stiffness. Taken together, these results suggest that the magnetic twisting cytometry technique may be a powerful and versatile tool for studies analyzing the molecular basis of transmembrane mechanical coupling to the cytoskeleton as well as dynamic relations between changes in cytoskeletal structure and alterations in cell form and function.

NASA Discipline Number 40-20↗

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Morozov, S. K.↗

Transonic Flutter Suppression Control Law Design Using Classical and Optimal Techniques with Wind-Tunnel Results

The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf. The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Mukhopadhyay, Vivek↗

Transonic Flutter Suppression Control Law Design, Analysis and Wind Tunnel Results

The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Mukhopadhyay, Vivek↗

Transonic Flutter Suppression Control Law Design, Analysis and Wind-Tunnel Results

The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf. The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Mukhopadhyay, Vivek↗

Transonic Flutter Suppression Control Law Design Using Classical and Optimal Techniques with Wind-Tunnel Results

The benchmark active controls technology and wind tunnel test program at NASA Langley Research Center was started with the objective to investigate the nonlinear, unsteady aerodynamics and active flutter suppression of wings in transonic flow. The paper will present the flutter suppression control law design process, numerical nonlinear simulation and wind tunnel test results for the NACA 0012 benchmark active control wing model. The flutter suppression control law design processes using (1) classical, (2) linear quadratic Gaussian (LQG), and (3) minimax techniques are described. A unified general formulation and solution for the LQG and minimax approaches, based on the steady state differential game theory is presented. Design considerations for improving the control law robustness and digital implementation are outlined. It was shown that simple control laws when properly designed based on physical principles, can suppress flutter with limited control power even in the presence of transonic shocks and flow separation. In wind tunnel tests in air and heavy gas medium, the closed-loop flutter dynamic pressure was increased to the tunnel upper limit of 200 psf. The control law robustness and performance predictions were verified in highly nonlinear flow conditions, gain and phase perturbations, and spoiler deployment. A non-design plunge instability condition was also successfully suppressed.

Mukhopadhyay, Vivek↗

Validation of High-Fidelity CFD/CAA Framework for Launch Vehicle Acoustic Environment Simulation against Scale Model Test Data

A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed for launch vehicle liftoff acoustic environment predictions. The framework couples the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin solver developed in the same production framework, Loci/THRUST, to accurately resolve and propagate acoustic physics across the entire launch environment. Time-accurate, Hybrid RANS/LES CFD modeling is applied for predicting the acoustic generation physics at the plume source, and a high-order accurate unstructured discontinuous Galerkin (DG) method is employed to propagate acoustic waves away from the source across large distances using high-order accurate schemes. The DG solver is capable of solving 2nd, 3rd, and 4th order Euler solutions for non-linear, conservative acoustic field propagation. Initial application testing and validation has been carried out against high resolution acoustic data from the Ares Scale Model Acoustic Test (ASMAT) series to evaluate the capabilities and production readiness of the CFD/CAA system to resolve the observed spectrum of acoustic frequency content. This paper presents results from this validation and outlines efforts to mature and improve the computational simulation framework.

Liever, Peter A.↗

Validation of High-Fidelity CFD/CAA Framework for Launch Vehicle Acoustic Environment Simulation against Scale Model Test Data

A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed for launch vehicle liftoff acoustic environment predictions. The framework couples the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate Discontinuous Galerkin solver developed in the same production framework, Loci/THRUST, to accurately resolve and propagate acoustic physics across the entire launch environment. Time-accurate, Hybrid RANS/LES CFD modeling is applied for predicting the acoustic generation physics at the plume source, and a high-order accurate unstructured mesh Discontinuous Galerkin (DG) method is employed to propagate acoustic waves away from the source across large distances using high-order accurate schemes. The DG solver is capable of solving 2nd, 3rd, and 4th order Euler solutions for non-linear, conservative acoustic field propagation. Initial application testing and validation has been carried out against high resolution acoustic data from the Ares Scale Model Acoustic Test (ASMAT) series to evaluate the capabilities and production readiness of the CFD/CAA system to resolve the observed spectrum of acoustic frequency content. This paper presents results from this validation and outlines efforts to mature and improve the computational simulation framework.

Liever, Peter A.↗

Challenges to Aboveground Biomass Prediction from Waveform Lidar

Accurate accounting of aboveground biomass density (AGBD) is crucial for carbon cycle, biodiversity, and climate change science. The Global Ecosystem Dynamics Investigation (GEDI), which maps global AGBD from waveform lidar, is the first of a new generation of Earth observation missions designed to improve carbon accounting. This paper explores the possibility that lidar waveforms may not be unique to AGBD—that forest stands with different AGBD may produce highly similar waveforms—and we hypothesize that non-uniqueness may contribute to the large uncertainties in AGBD predictions. Our analysis integrates simulated GEDI waveforms from 428 in situ stem maps with output from an individual-based forest gap model, which we use to generate a database of potential forest stands and simulate GEDI waveforms from those stands. We use this database to predict the AGBD of the 428 in situ stem maps via two different methods: a linear regression from waveform metrics, and a waveform-matching approach that accounts for waveform-AGBD non-uniqueness. We find that some in situ waveforms are more unique to AGBD than others, which notably impacts AGBD prediction uncertainty (7–411 Mg ha−1, average of 167 Mg ha−1). We also find that forest structure complexity may influence the non-uniqueness effect; stands with low structural complexity are more unique to AGBD than more mature stands with multiple cohorts and canopy layers. These findings suggest that the non-uniqueness phenomena may be introduced by the measuring characteristics of waveform lidar in combination with how forest structure manifests at small scales, and we discuss how this complexity may complicate uncertainty estimation in AGBD prediction. This analysis suggests a limit to the accuracy and precision of AGBD predictions from lidar waveforms seen in empirical studies, and underscores the need for further exploration of the relationships between lidar remote sensing measurements, forest structure, and AGBD.

GEDI↗

Mesospheric Non-Migrating Tides Generated With Planetary Waves: Characteristics - 1

We discuss results from a modeling study with our Numerical Spectral Model (NSM) that specifically deals with the non-migrating tides generated in the mesosphere. The NSM extends from the ground to the thermosphere, incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GWs), and it describes the major dynamical features of the atmosphere including the wave driven equatorial oscillations (QBO and SAO), and the seasonal variations of tides and planetary waves. Accounting solely for the excitation sources of the solar migrating tides, the NSM generates through dynamical interactions also non-migrating tides in the mesosphere that are comparable in magnitude to those observed. Large non-migrating tides are produced in the diurnal and semi-diurnal oscillations for the zonal mean (m = 0) and in the semidiurnal oscillation for m = 1. In general, significant eastward and westward propagating tides are generated for all the zonal wave numbers m = 1 to 4. To identify the cause, the NSM is run without the solar heating for the zonal mean (m = 0), and the amplitudes of the resulting non-migrating tides are then negligibly small. In this case, the planetary waves are artificially suppressed, which are generated in the NSM through instabilities. This leads to the conclusion that the non-migrating tides are generated through non-linear interactions between planetary waves and migrating tides, as Forbes et al. and Talaat and Liberman had proposed. In an accompanying paper, we present results from numerical experiments, which indicate that gravity wave filtering contributes significantly to produce the non-linear coupling that is involved.

Mayr, H. G.↗

CLASP2: The Chromospheric LAyer Spectro-Polarimeter

A major remaining challenge for heliophysicsis to decipher the magnetic structure of the chromosphere, due to its "large role in defining how energy is transported into the corona and solar wind" (NASA's Heliophysics Roadmap). Recent observational advances enabled by the Interface Region Imaging Spectrometer (IRIS) have revolutionized our view of the critical role this highly dynamic interface between the photosphere and corona plays in energizing and structuring the outer solar atmosphere. Despite these advances, a major impediment to better understanding the solar atmosphere is our lack of empirical knowledge regarding the direction and strength of the magnetic field in the upper chromosphere. Such measurements are crucial to address several major unresolved issues in solar physics: for example, to constrain the energy flux carried by the Alfven waves propagating through the chromosphere (De Pontieuet al., 2014), and to determine the height at which the plasma Beta = 1 transition occurs, which has important consequences for the braiding of magnetic fields (Cirtainet al., 2013; Guerreiroet al., 2014), for propagation and mode conversion of waves (Tian et al., 2014a; Straus et al., 2008) and for non-linear force-free extrapolation methods that are key to determining what drives instabilities such as flares or coronal mass ejections (e.g.,De Rosa et al., 2009). The most reliable method used to determine the solar magnetic field vector is the observation and interpretation of polarization signals in spectral lines, associated with the Zeeman and Hanle effects. Magnetically sensitive ultraviolet spectral lines formed in the upper chromosphere and transition region provide a powerful tool with which to probe this key boundary region (e.g., Trujillo Bueno, 2014). Probing the magnetic nature of the chromosphere requires measurement of the Stokes I, Q, U and V profiles of the relevant spectral lines (of which Q, U and V encode the magnetic field information).

Chromosphere↗

CLASP2: The Chromospheric LAyer Spectro-Polarimeter

A major remaining challenge for heliophysicsis to decipher the magnetic structure of the chromosphere, due to its 'large role in defining how energy is transported into the corona and solar wind' (NASA's Heliophysics Roadmap). Recent observational advances enabled by the Interface Region Imaging Spectrometer (IRIS) have revolutionized our view of the critical role this highly dynamic interface between the photosphere and corona plays in energizing and structuring the outer solar atmosphere. Despite these advances, a major impediment to better understanding the solar atmosphere is our lack of empirical knowledge regarding the direction and strength of the magnetic field in the upper chromosphere. Such measurements are crucial to address several major unresolved issues in solar physics: for example, to constrain the energy flux carried by the Alfven waves propagating through the chromosphere (De Pontieuet al., 2014), and to determine the height at which the plasma β = 1 transition occurs, which has important consequences for the braiding of magnetic fields (Cirtainet al., 2013; Guerreiroet al., 2014), for propagation and mode conversion of waves (Tian et al., 2014a; Straus et al., 2008) and for non-linear force-free extrapolation methods that are key to determining what drives instabilities such as flares or coronal mass ejections (e.g., De Rosa et al., 2009). The most reliable method used to determine the solar magnetic field vector is the observation and interpretation of polarization signals in spectral lines, associated with the Zeeman and Hanle effects. Magnetically sensitive ultraviolet spectral lines formed in the upper chromosphere and transition region provide a powerful tool with which to probe this key boundary region (e.g., Trujillo Bueno, 2014). Probing the magnetic nature of the chromosphere requires measurement of the Stokes I, Q, U and V profiles of the relevant spectral lines (of which Q, U and V encode the magnetic field information).

Solar; Chromosphere; Polarization↗

High-Order Shock-Capturing Methods for Modeling Dynamics of the Solar Atmosphere

We use one-dimensional high-order central shock capturing numerical methods to study the response of various model solar atmospheres to forcing at the solar surface. The dynamics of the atmosphere is modeled with the Euler equations in a variable-sized flux tube in the presence of gravity. We study dynamics of the atmosphere suggestive of spicule formation and coronal oscillations. These studies are performed on observationally-derived model atmospheres above the quiet sun and above sunspots. To perform these simulations, we provide a new extension of existing second- and third- order shock-capturing methods to irregular grids. We also solve the problem of numerically maintaining initial hydrostatic balance via the introduction of new variables in the model equations and a careful initialization mechanism. We find several striking results: all model atmospheres respond to a single impulsive perturbation with several strong shock waves consistent with the rebound-shock model. These shock waves lift material and the transition region well into the initial corona, and the sensitivity of this lift to the initial impulse depends non-linearly on the details of the atmosphere model. We also reproduce an observed 3-minute coronal oscillation above sunspots compared to 5-minute oscillations above the quiet sun.

Bryson, Steve↗