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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Optimal Energy Measurement in Nonlinear Systems: An Application of Differential Geometry

Design of TES microcalorimeters requires a tradeoff between resolution and dynamic range. Often, experimenters will require linearity for the highest energy signals, which requires additional heat capacity be added to the detector. This results in a reduction of low energy resolution in the detector. We derive and demonstrate an algorithm that allows operation far into the nonlinear regime with little loss in spectral resolution. We use a least squares optimal filter that varies with photon energy to accommodate the nonlinearity of the detector and the non-stationarity of the noise. The fitting process we use can be seen as an application of differential geometry. This recognition provides a set of well-developed tools to extend our work to more complex situations. The proper calibration of a nonlinear microcalorimeter requires a source with densely spaced narrow lines. A pulsed laser multi-photon source is used here, and is seen to be a powerful tool for allowing us to develop practical systems with significant detector nonlinearity. The combination of our analysis techniques and the multi-photon laser source create a powerful tool for increasing the performance of future TES microcalorimeters.

Nonlinear response↗
NASA NTRS2022

A Quantum Algorithm to Simulate Open Quantum Systems

Given the advent of quantum algorithms for a wide array of problems in linear algebra and machine learning, it is important to develop general methods for the simulation of arbitrary (ie non-unitary) operators on quantum hardware. In this talk, we present a novel quantum algorithm based on the quantum singular value transformation (QSVT) to apply an arbitrary operator K to some input state and subsequently estimate the expectation value of some observable. Our construction then immediately yields a route to estimating observables of states undergoing open quantum dynamics, whose effect is captured by a set of non-unitary Kraus operators. Our algorithm succeeds deterministically given the Sz-Nagy dilation, and we provide details on the algorithm's query and gate complexity, numerical verification, and comparisons with prior methods.

Quantum computing↗
NASA NTRSDate not supplied

Pterodactyl: 6-DOF Integration of Guidance and Control Algorithms in Genesis

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles. This paper builds upon the Pterodactyl architecture that employed eight individually articulating flaps with two options for guidance, bank angle modulation with the Fully Numerical Predictor Corrector Entry Guidance (FNPEG) technique and angle of attack and sideslip modulation with FNPEG uncoupled range control. These, with a linear quadratic regulator controller, have previously been presented in separate 3-DOF trajectory simulations, one for translational motion and another for rotational dynamics. This work will show results for integrated 6-DOF simulations, leveraging recent advancements in trajectory simulation tools, namely the Julia-based Genesis flight simulation tool. Results show similar performance for both schemes when compared to the previously presented 3-DOF results once the controllers were adequately tuned. Angle of attack and sideslip angle modulation necessitated a controller designed at high dynamic pressure conditions, and bank angle modulation necessitated a controller designed at lower dynamic pressure conditions. Notably, the angle of attack and sideslip modulation approach could not achieve the parachute deploy point target Mach number of 2, when only tuning the controller gains, so future work will be needed to improve guidance and control tracking.

EDL↗
NASA NTRS2001

Troposphere-Stratosphere Dynamic Coupling Under Strong and Weak Polar Vortex Conditions

The relationship between Northern Hemisphere (NH) tropospheric and stratospheric wave-like anomalies of spherical zonal wave number (ZWN) 1 is studied by applying Canonical Correlation Analysis (CCA). A lag-correlation technique is used with 10-day lowpass filtered daily time series of 50- and 500-hPa geopotential heights. Generally stratospheric circulation is determined by ultralong tropospheric planetary waves. During winter seasons characterized either by any anomalously strong or weak polar winter vortex different propagation characteristics for waves of ZWN 1 are observed. The non-linear perspective of the results have implications for medium range weather forecast and climate sensitivity experiments.

Perlwitz, Judith↗
NASA NTRSDate not supplied

Pterodactyl: 6-DOF Integration of Guidance and Control Algorithms in Genesis

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles. This paper builds upon the Pterodactyl architecture that employed eight individually articulating flaps with two options for guidance, bank angle modulation with the Fully Numerical Predictor Corrector Entry Guidance (FNPEG) technique and angle of attack and sideslip modulation with FNPEG uncoupled range control. These, with a linear quadratic regulator controller, have previously been presented in separate 3-DOF trajectory simulations, one for translational motion and another for rotational dynamics. This work will show results for integrated 6-DOF simulations, leveraging recent advancements in trajectory simulation tools, namely the Julia-based Genesis package. Results show similar performance for both schemes when compared to the previously presented 3-DOF results once the controllers were adequately tuned. Angle of attack and sideslip modulation necessitated a controller designed at high dynamic pressure conditions, and bank angle modulation necessitated a controller designed at lower dynamic pressure conditions. Notably, angle of attack and sideslip modulation could not achieve the parachute deploy point target Mach number of 2, when only tuning the controller gains, thus future work will be needed to tune guidance specific parameters to achieve acceptable range targeting and guidance command tracking.

DEV↗
NASA NTRSDate not supplied

Efficient Parametric Uncertainty Analysis of an Earth Entry Vehicle Concept Using Least Angle Regression

The objective of this work was to outline and apply an efficient and accurate parametric un-certainty propagation approach to the analysis of convective heating on an Earth entry vehicle concept. The described approach was based on Least Angle Regression used to solve a sparse and underdetermined linear system in the point-collocation non-intrusive polynomial chaos surrogate method. This approach involved an iterative process to computing the non-zero terms of the underlying polynomial chaos model using only enough samples to converge uncertainty interval predictions and Sobol index values based global nonlinear sensitivity estimates. The Earth entry vehicle was analyzed at three points along a representative trajectory for a Mars return mission. 329 sources of uncertainty were identified in the computational fluid dynamics model used to predict the forebody convective heating. These included uncertainty in flow field chemical rates, collision integrals, heats of formation, surface finite rate char model reaction rates, wall roughness height, and the turbulent Schmidt number. Results from this study showed that convective heating uncertainty as high as 50% of the nominal was predicted with only about 50 evaluations of the computational model. This was far fewer than would be required for a sampling-based approach or a full basis polynomial chaos model, which would have required over 50,000 samples. Additionally, results showed that over 90% of the total convective heating uncertainty was due to uncertainty in the N2catalytic rate on the surface, while the remainder of the uncertainty was attributed to the turbulent Schmidt number and the wall roughness uncertainties.

Thomas K West IV↗
NASA NTRSDate not supplied

Quantification of Numerical Uncertainty via Nonlinear Dynamical Approach

Motivations (Ensure a Higher Level of Confidence in the Predictability & Reliability of Numerical Simulation for Multiscale Complex Nonlinear Fluid Problems) - The last two decades have been an era when computation is ahead of analysis & when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems & non-traditional fields (Especially when computations offer the ONLY way of generating this type of data limited simulations). - At present some of the numerical uncertainties can be explained and minimized by traditional numerical analysis and standard CFD practices. However, such practices, usually based on linearized analysis, MIGHT NOT be sufficient for strongly nonlinear and/or stiff problems. - We need a good understanding of the nonlinear behavior of numerical schemes being used as an integral part of code verification, validation and certification.

HEC↗
NASA NTRS2011

Atomic Simulation of Complex DNA DSBs and the Interactions with the Ku70/80 Heterodimer

DNA double strand breaks (DSBs) induced by ionizing radiation (IR) usually contain modified bases such as 8-oxo-7,8-dihydroguanine (8-oxoG) and thymine glycol, apurinic/apyrimidinic (AP) sites, 2-deoxyribonolactone, or single-strand breaks (SSBs). The presence of such lesions in close proximity to the DSB terminus makes the DNA nicks more difficult to repair and rejoin than endogenously induced simple DSBs, and as such a major determinant of the biological effects of high linear energy transfer (LET) radiation as encountered in space travel. In this study we conducted molecular dynamics simulations on a series of DNA duplexes with various complex lesions of 8-oxoG and AP sites, in an effort to investigate the effects of such lesions to the structural integrity and stability of DNA after insulted by IR. We also simulated the interaction of such complex DSBs with the Ku70/80 heterodimer, the first protein in mammalian cells to embark the non-homologous end joining (NHEJ) DNA repair pathway. The results indicate, compared to DNA with simple DSBs, the complex lesions can enhance the hydrogen bonds opening rate at the DNA terminus, and increase the mobility of the whole duplex, thus they present more deleterious effects to the genome integrity if not captured and repaired promptly in cells. Simulations also demonstrate the binding of Ku drastically reduces structural disruption and flexibility caused by the complex lesions, and the interactions of Ku with complex DSBs have a different potential energy landscape from the bound structure with simple DSB. In all complex DSBs systems, the binding of DSB terminus with Ku70 is softened while the binding of the middle duplex with Ku80 is tightened. This energy shift may help the Ku protein to secure at the DSB terminus for a longer time, so that other end processing factors or repair pathways can proceed at the lesions before NHEJ repair process starts. These atomic simulations may provide valuable new insight into the selective action of repair proteins on damaged DNA.

Hu, Shaowen↗
NASA NTRS2005

Simultaneous quantification of poly-dispersed anionic, amphoteric and nonionic surfactants in simulated wastewater samples using C18 high-performance liquid chromatography-quadrupole ion-trap mass spectrometry

This paper describes the development of a guantitative method for direct and simultaneous determination of three frequently encountered surfactants, amphoteric (cocoamphoacetate, CAA), anionic (sodium laureth sulfate, SLES), and nonionic (alcohol ethoxylate, AE) using a reversed-phase C18 HPLC coupled with an ESI ion-trap mass spectrometer (MS). Chemical composition, ionization characteristics and fragmentation pathways of the surfactants are presented. Positive ESI was effective for all three surfactants in agueous methanol buffered with ammonium acetate. The method enables rapid determinations in small sample volumes containing inorganic salts (up to 3.5 g L(-1)) and multiple classes of surfactants with high specificity by applying surfactant specific tandem mass spectrometric strategies. It has dynamic linear ranges of 2-60, 1.5-40, 0.8-56 mg L(-1) with R2 egual or greater than 0.999, 0.98 and 0.999 (10 microL injection) for CAA, SLES, and AE, respectively.

Chromatography, High Pressure Liquid/methods↗
NASA NTRS1988

On-orbit parametric identification methodology

On-orbit system identification (ID) of large space systems is essential for various reasons. For example, the complex composite structure of such systems cannot be ground-tested; their structural dynamic characteristics must be known accurately in order to accomplish active control. Furthermore, such capability can be used to characterize/identify various disturbances. The identification process is consisted of four principal elements: (1) modeling, (2) the estimation algorithm, (3) input system, and (4) measurement system. These elements are highly correlated and all togerher determine the success of the identification problem. Accurate modeling of large space systems is the most important element of the identification process. Large flexible structures are non-linear and infinite dimensional systems with highly coupled parameters and low frequency packed modes. In addition, these systems are subject to stochastic and time-varying disturbances, they have structural parameters which can vary due to on-orbit assembly deployment, and operations. These systems are generally; however, represented by constant coefficient, finite order differential equations. The non-linearities, coupling and noise effects are also often neglected. Moreover, identification experiment designs which lead to highly complex optimization problems usually require the simultaneous choice of ID algorithm, sensor, and actuator type and placement. On-orbit bandwidth and power restrictions on excitation, limited data window, and restrictions on sensor/actuator type, placement and number, has led to practical questions of implementations.

Hadaegh, Fred Y.↗
NASA NTRS1995

Linear operating region in the ozone dial photon counting system

Ozone is a relatively unstable molecule found in Earth's atmosphere. An ozone molecule is made up of three atoms of oxygen. Depending on where ozone resides, it can protect or harm life on Earth. High in the atmosphere, about 15 miles up, ozone acts as a shield to protect Earth's surface from the sun's harmful ultraviolet radiation. Without this shield, we would be more susceptible to skin cancer, cataracts, and impaired immune systems. Closer to Earth, in the air we breathe, ozone is a harmful pollutant that causes damage to lung tissue and plants. Since the early 1980's, airborne lidar systems have been used for making measurements of ozone. The differential absorption lidar (DIAL) technique is used in the remote measurement of O3. This system allows the O3 to be measured as function of the range in the atmosphere. Two frequency-doubled Nd:YAG lasers are used to pump tunable dye lasers. The lasers are operating at 289 nm for the DIAL on-line wavelength of O3, and the other one is operated at 300 nm for the off-line wavelength. The DIAL wavelengths are produced in sequential laser pulses with a time separation of 300 micro s. The backscattered laser energy is collected by telescopes and measured using photon counting systems. The photon counting system measures the light signal by making use of the photon nature of light. The output pulse from the Photo-Multiplier Tube (PE), caused by a photon striking the PMT photo-cathode, is amplified and passed to a pulse height discriminator. The peak value of the pulse is compared to a reference voltage (discrimination level). If the pulse amplitude exceeds the discrimination level, the discriminator generates a standard pulse which is counted by the digital counter. Non-linearity in the system is caused by the overlapping of pulses and the finite response time of the electronics. At low count rates one expects the system to register one event for each output pulse from the PMT corresponding to a photon incident upon the photocathode, however, at higher rates the limitations of the discrimination/counting system will cause the observed count rate to be non-linear with respect to the true count rate. Depending on the pulse height distribution and the discriminator level, the overlapping of pulses (pulse pile-up) can cause count loss or even an additional apparent count gain as the signal levels increase. Characterization of the system, including the pulse height distribution, the signal to noise ratio, and the effect of the discriminator threshold level, is critical in maximizing the linear operating region of the system, thus greatly increasing the useful dynamic range of the system.

Andrawis, Madeleine↗
NASA NTRS2012

Orion MPCV Service Module Avionics Ring Pallet Testing, Correlation, and Analysis

The NASA Orion Multi-Purpose Crew Vehicle (MPCV) is being designed to replace the Space Shuttle as the main manned spacecraft for the agency. Based on the predicted environments in the Service Module avionics ring, an isolation system was deemed necessary to protect the avionics packages carried by the spacecraft. Impact, sinusoidal, and random vibration testing were conducted on a prototype Orion Service Module avionics pallet in March 2010 at the NASA Glenn Research Center Structural Dynamics Laboratory (SDL). The pallet design utilized wire rope isolators to reduce the vibration levels seen by the avionics packages. The current pallet design utilizes the same wire rope isolators (M6-120-10) that were tested in March 2010. In an effort to save cost and schedule, the Finite Element Models of the prototype pallet tested in March 2010 were correlated. Frequency Response Function (FRF) comparisons, mode shape and frequency were all part of the correlation process. The non-linear behavior and the modeling the wire rope isolators proved to be the most difficult part of the correlation process. The correlated models of the wire rope isolators were taken from the prototype design and integrated into the current design for future frequency response analysis and component environment specification.

Staab, Lucas↗
NASA NTRS2018

STEREO Observations of Interplanetary Coronal Mass Ejections in 2007-2016

We have conducted a survey of 341 interplanetary coronal mass ejections (ICMEs) using STEREO A/B data, analyzing their properties while extending a Level 3 product through 2016. Among the 192 ICMEs with distinguishable sheath region and magnetic obstacle, the magnetic field maxima in the two regions are comparable, and the dynamic pressure peaks mostly in the sheath. The north/south direction of the magnetic field does not present any clear relationship between the sheath region and the magnetic obstacle. About 71% of ICMEs are expanding at 1 au, and their expansion speed varies roughly linearly with their maximum speed except for ICMEs faster than 700 km/s. The total pressure generally peaks near the middle of the well-defined magnetic cloud (MC) passage, while it often declines along with the non-MC ICME passage, consistent with our previous interpretation concerning the effects of sampling geometry on what is observed. The hourly average iron charge state reaches above 12+ approx. 31% of the time for MCs, approx. 16% of the time for non-MC ICMEs, and approx. 1% of the time for non- ICME solar wind. In four ICMEs abrupt deviations of the magnetic field from the nominal field rotations occur in the magnetic obstacles, coincident with a brief drop or increase in field strength-features could be related to the interaction with dust. In comparison with the similar phases of solar cycle 23, the STEREO ICMEs in this cycle occur less often and are generally weaker and slower, although their field and pressure compressions weaken less than the background solar wind.

Jian, L. K.↗
NASA NTRS1990

A nonrecursive order N preconditioned conjugate gradient: Range space formulation of MDOF dynamics

While excellent progress has been made in deriving algorithms that are efficient for certain combinations of system topologies and concurrent multiprocessing hardware, several issues must be resolved to incorporate transient simulation in the control design process for large space structures. Specifically, strategies must be developed that are applicable to systems with numerous degrees of freedom. In addition, the algorithms must have a growth potential in that they must also be amenable to implementation on forthcoming parallel system architectures. For mechanical system simulation, this fact implies that algorithms are required that induce parallelism on a fine scale, suitable for the emerging class of highly parallel processors; and transient simulation methods must be automatically load balancing for a wider collection of system topologies and hardware configurations. These problems are addressed by employing a combination range space/preconditioned conjugate gradient formulation of multi-degree-of-freedom dynamics. The method described has several advantages. In a sequential computing environment, the method has the features that: by employing regular ordering of the system connectivity graph, an extremely efficient preconditioner can be derived from the 'range space metric', as opposed to the system coefficient matrix; because of the effectiveness of the preconditioner, preliminary studies indicate that the method can achieve performance rates that depend linearly upon the number of substructures, hence the title 'Order N'; and the method is non-assembling. Furthermore, the approach is promising as a potential parallel processing algorithm in that the method exhibits a fine parallel granularity suitable for a wide collection of combinations of physical system topologies/computer architectures; and the method is easily load balanced among processors, and does not rely upon system topology to induce parallelism.

Kurdila, Andrew J.↗
NASA NTRS2002

Solar System Dynamics

In these 18 years, the research has touched every major dynamical problem in the solar system, including: the effect of chaotic zones on the distribution of asteroids, the delivery of meteorites along chaotic pathways, the chaotic motion of Pluto, the chaotic motion of the outer planets and that of the whole solar system, the delivery of short period comets from the Kuiper belt, the tidal evolution of the Uranian arid Galilean satellites, the chaotic tumbling of Hyperion and other irregular satellites, the large chaotic variations of the obliquity of Mars, the evolution of the Earth-Moon system, and the resonant core- mantle dynamics of Earth and Venus. It has introduced new analytical and numerical tools that are in widespread use. Today, nearly every long-term integration of our solar system, its subsystems, and other solar systems uses algorithms that was invented. This research has all been primarily Supported by this sequence of PGG NASA grants. During this period published major investigations of tidal evolution of the Earth-Moon system and of the passage of the Earth and Venus through non-linear core-mantle resonances were completed. It has published a major innovation in symplectic algorithms: the symplectic corrector. A paper was completed on non-perturbative hydrostatic equilibrium.

Wisdom, Jack↗
NASA NTRS1994

Gender- and age-related differences in heart rate dynamics: are women more complex than men?

OBJECTIVES. This study aimed to quantify the complex dynamics of beat-to-beat sinus rhythm heart rate fluctuations and to determine their differences as a function of gender and age. BACKGROUND. Recently, measures of heart rate variability and the nonlinear "complexity" of heart rate dynamics have been used as indicators of cardiovascular health. Because women have lower cardiovascular risk and greater longevity than men, we postulated that there are important gender-related differences in beat-to-beat heart rate dynamics. METHODS. We analyzed heart rate dynamics during 8-min segments of continuous electrocardiographic recording in healthy young (20 to 39 years old), middle-aged (40 to 64 years old) and elderly (65 to 90 years old) men (n = 40) and women (n = 27) while they performed spontaneous and metronomic (15 breaths/min) breathing. Relatively high (0.15 to 0.40 Hz) and low (0.01 to 0.15 Hz) frequency components of heart rate variability were computed using spectral analysis. The overall "complexity" of each heart rate time series was quantified by its approximate entropy, a measure of regularity derived from nonlinear dynamics ("chaos" theory). RESULTS. Mean heart rate did not differ between the age groups or genders. High frequency heart rate power and the high/low frequency power ratio decreased with age in both men and women (p < 0.05). The high/low frequency power ratio during spontaneous and metronomic breathing was greater in women than men (p < 0.05). Heart rate approximate entropy decreased with age and was higher in women than men (p < 0.05). CONCLUSIONS. High frequency heart rate spectral power (associated with parasympathetic activity) and the overall complexity of heart rate dynamics are higher in women than men. These complementary findings indicate the need to account for gender-as well as age-related differences in heart rate dynamics. Whether these gender differences are related to lower cardiovascular disease risk and greater longevity in women requires further study.

NASA Discipline Cardiopulmonary↗
NASA NTRS2012

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗
NASA NTRS1987

The electrostatic effect for the collisionless tearing mode

Electron dynamics has not been self-consistently considered in collisionless tearing mode theories to date because of the mathematical complexity of the Vlasov-Maxwell equations. It has been found, using computer simulations, that electrostatic fields play an important role in the tearing mode. Vlasov theory, including the electrostatic field, is investigated for topologies with both antiparallel and nonantiparallel magnetic field lines. The electrostatic field influences the resonant current in the neutral sheet which is a non-MHD effect, and modifies the linear growth rate. At the magnetopause, where the field lines are not antiparallel, the electrostatic effect acts to raise the linear growth rate of the tearing mode. On the other hand, in the magnetotail, where magnetic field lines are antiparallel, the electrostatic effect reduces the tearing mode growth rate.

Hoshino, Masahiro↗