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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions. Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz model as well as more realistic models of the means and atmosphere. A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter situations to allow for correct update of the ensemble members. The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to be quite puzzling in that results state estimates are worse than for their filter analogue. In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use the Lorenz model to test and compare the behavior of a variety of implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqui↗
NASA NTRS1999

The Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known Gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions (e.g., Miller 1994). Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz (1963) model as well as more realistic models of the oceans (Evensen and van Leeuwen 1996) and atmosphere (Houtekamer and Mitchell 1998). A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter equations to allow for correct update of the ensemble members (Burgers 1998). The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to quite puzzling in that results of state estimate are worse than for their filter analogue (Evensen 1997). In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use Lorenz (1963) model to test and compare the behavior of a variety implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqiu↗
NASA NTRS2002

The Sun's Impact on Climate

We provide an overview of the impact of the Sun on the Earth atmosphere and climate system, focused on heating of Earth's atmosphere and oceans. We emphasize the importance of the spectral measurements of SIM and SOLSTICE- that we must know how solar variations are distributed over ultraviolet, visible, and infrared wavelengths, since these have separate characteristic influences on Earth's ozone layer, clouds, and upper layers of the oceans. Emphasis is also given to understanding both direct and indirect influences of the Sun on the Earth, which involve feedbacks between Earth's stratosphere, troposphere, and oceans, each with unique time scales, dynamics, chemistry, and biology, interacting non-linearly. Especially crucial is the role of all three phases of water on Earth, water vapor being the primary greenhouse gas in the atmosphere, the importance of trace gases such as CO2 arising from their absorption in the "water vapor window" at 800 - 1250/cm (12.5 to 8 microns). Melting of polar ice is one major response to the post-industrial global warming, enhanced due to "ice-albedo" feedback. Finally, water in liquid form has a major influence due to cloud albedo feedback, and also due to the oceans' absorption of solar radiation, particularly at visible wavelengths, through the visible "liquid water window" that allows penetration of visible light deep into the mixed layer, while nearby ultraviolet and infrared wavelengths do not penetrate past the upper centimeter ocean surface skin layer. A large fraction of solar energy absorbed by the oceans goes into the latent heat of evaporation. Thus the solar heating of the atmosphere-ocean system is strongly coupled through the water cycle of evaporation, cloud formation, precipitation, surface runoff and ice formation, to Earth's energy budget and climate, each different climate component responding to variations in different solar spectral bands, at ultraviolet, visible and infrared wavelengths.

Cahalan, Robert↗
NASA NTRS2005

Single-Scale Retinex Using Digital Signal Processors

The Retinex is an image enhancement algorithm that improves the brightness, contrast and sharpness of an image. It performs a non-linear spatial/spectral transform that provides simultaneous dynamic range compression and color constancy. It has been used for a wide variety of applications ranging from aviation safety to general purpose photography. Many potential applications require the use of Retinex processing at video frame rates. This is difficult to achieve with general purpose processors because the algorithm contains a large number of complex computations and data transfers. In addition, many of these applications also constrain the potential architectures to embedded processors to save power, weight and cost. Thus we have focused on digital signal processors (DSPs) and field programmable gate arrays (FPGAs) as potential solutions for real-time Retinex processing. In previous efforts we attained a 21 (full) frame per second (fps) processing rate for the single-scale monochromatic Retinex with a TMS320C6711 DSP operating at 150 MHz. This was achieved after several significant code improvements and optimizations. Since then we have migrated our design to the slightly more powerful TMS320C6713 DSP and the fixed point TMS320DM642 DSP. In this paper we briefly discuss the Retinex algorithm, the performance of the algorithm executing on the TMS320C6713 and the TMS320DM642, and compare the results with the TMS320C6711.

Hines, Glenn↗
NASA NTRS2013

Trajectory Reconstruction and Uncertainty Analysis Using Mars Science Laboratory Pre-Flight Scale Model Aeroballistic Testing

As part of the Mars Science Laboratory (MSL) trajectory reconstruction effort at NASA Langley Research Center, free-flight aeroballistic experiments of instrumented MSL scale models was conducted at Aberdeen Proving Ground in Maryland. The models carried an inertial measurement unit (IMU) and a flush air data system (FADS) similar to the MSL Entry Atmospheric Data System (MEADS) that provided data types similar to those from the MSL entry. Multiple sources of redundant data were available, including tracking radar and on-board magnetometers. These experimental data enabled the testing and validation of the various tools and methodologies that will be used for MSL trajectory reconstruction. The aerodynamic parameters Mach number, angle of attack, and sideslip angle were estimated using minimum variance with a priori to combine the pressure data and pre-flight computational fluid dynamics (CFD) data. Both linear and non-linear pressure model terms were also estimated for each pressure transducer as a measure of the errors introduced by CFD and transducer calibration. Parameter uncertainties were estimated using a "consider parameters" approach.

Lugo, Rafael A.↗
NASA NTRS2015

Non-Linear Slosh Damping Model Development and Validation

Propellant tank slosh dynamics are typically represented by a mechanical model of spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control (GN&C) analysis. For a partially-filled smooth wall propellant tank, the critical damping based on classical empirical correlation is as low as 0.05%. Due to this low value of damping, propellant slosh is potential sources of disturbance critical to the stability of launch and space vehicles. It is postulated that the commonly quoted slosh damping is valid only under the linear regime where the slosh amplitude is small. With the increase of slosh amplitude, the critical damping value should also increase. If this nonlinearity can be verified and validated, the slosh stability margin can be significantly improved, and the level of conservatism maintained in the GN&C analysis can be lessened. The purpose of this study is to explore and to quantify the dependence of slosh damping with slosh amplitude. Accurately predicting the extremely low damping value of a smooth wall tank is very challenging for any Computational Fluid Dynamics (CFD) tool. One must resolve thin boundary layers near the wall and limit numerical damping to minimum. This computational study demonstrates that with proper grid resolution, CFD can indeed accurately predict the low damping physics from smooth walls under the linear regime. Comparisons of extracted damping values with experimental data for different tank sizes show very good agreements. Numerical simulations confirm that slosh damping is indeed a function of slosh amplitude. When slosh amplitude is low, the damping ratio is essentially constant, which is consistent with the empirical correlation. Once the amplitude reaches a critical value, the damping ratio becomes a linearly increasing function of the slosh amplitude. A follow-on experiment validated the developed nonlinear damping relationship. This discovery can lead to significant savings by reducing the number and size of slosh baffles in liquid propellant tanks.

Yang, H. Q.↗
DOE OSTI2025

Using "AI Poincare" to analyze non-linear integrable optics

This study dives into the applicability of using automated discovery of conserved quantities in dynamical systems relevant to accelerator physics. Specifically, we explore the performance of AI Poincaré in analyzing numerical trajectory data obtained using the McMillan system of non-linear integrable optics. A comprehensive evaluation of the algorithm's performance is conducted through diverse methodologies. These include the analysis of the estimated number of conserved quantities embedded in a dataset and the deviation of interpolated points on the inferred manifold with respect to points in actually in the dataset. the investigation identifies an optimal range of perturbation distances where the underlying manifold extraction algorithm inside AI Poincaré exhibits optimal performance. Additionally, an improved neural network architecture is proposed based on the observed results. Finally, we apply the algorithm to preliminary experimental data from the Integrable Optics Test Accelerator at Fermilab to successfully infer the number of conserved quantities even in the presence of fast decoherence of the measured signal.

Osmanov, Lazare [Free U. Tbilisi]↗
DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗
DOE OSTI2020

Theoretical Development of Controller Transfer applied to Dynamical Systems

Data-driven control strategies suffer from the necessity of training and its associated computational time and cost. This motivates the continuous increase of popularity of transfer learning based methods. In this paper, we have investigated transfer learning in the context of controller for dynamical systems. First, we have derived conditions under which controller can be transferred between two linear dynamical systems. Furthermore, we have introduced a novel algorithm to design the transferable and non-transferable control components for linear systems. Secondly, we have derived conditions under which a stabilizing controller can be transferred between two different nonlinear dynamical systems. Furthermore, we have used the notion of zero dynamics to transform the nonlinear dynamics into global normal form, to derive the transferability conditions on the controller. Finally, we have numerically evaluated the performance of our transfer learning based methods for two sets of example linear and nonlinear dynamical systems.

Chakraborty, Indrasis↗
NASA NTRS1995

Heart rate and respiratory rhythm dynamics on ascent to high altitude

OBJECTIVE--To investigate the alterations in autonomic control of heart rate at high altitude and to test the hypothesis that hypoxaemic stress during exposure to high altitude induces non-linear, periodic heart rate oscillations, similar to those seen in heart failure and the sleep apnoea syndrome. SUBJECTS--11 healthy subjects aged 24-64. MAIN OUTCOME MEASURES--24 hour ambulatory electrocardiogram records obtained at baseline (1524 m) and at 4700 m. Simultaneous heart rate and respiratory dynamics during 2.5 hours of sleep by fast Fourier transform analysis of beat to beat heart rate and of an electrocardiographically derived respiration signal. RESULTS--All subjects had resting hypoxaemia at high altitude, with an average oxyhaemoglobin saturation of 81% (5%). There was no significant change in mean heart rate, but low frequency (0.01-0.05 Hz) spectral power was increased (P < 0.01) at high altitude. Time series analysis showed a complex range of non-linear sinus rhythm dynamics. Striking low frequency (0.04-0.06 Hz) heart rate oscillations were observed during sleep in eight subjects at high altitude. Analysis of the electrocardiographically derived respiration signal indicated that these heart rate oscillations correlated with low frequency respiratory oscillations. CONCLUSIONS--These data suggest (a) that increased low frequency power during high altitude exposure is not simply attributable to increased sympathetic modulation of heart rate, but relates to distinctive cardiopulmonary oscillations at approximately 0.05 Hz and (b) that the emergence of periodic heart rate oscillations at high altitude is consistent with an unstable cardiopulmonary control system that may develop on acute exposure to hypoxaemic stress.

Non-NASA Center↗
NASA NTRS1994

A fundamental aeroservoelastic study combining unsteady CFD with adaptive control

This paper describes a two-dimensional aeroservoelastic study in the time domain. The model, which is based on exact inviscid aerodynamics, correctly represents the large amplitude motions and the associated strong shock dynamics in the transonic regime. The aeroservoelastic system consists of a two degree-of-freedom airfoil with a trailing edge control surface. Using first-order actuator dynamics, a digital adaptive controller is applied to provide active flutter suppression. Comparisons between time-responses of the open-loop and closed loop systems show the ability of the trailing edge control surface to suppress non-linear transonic aeroelastic phenomena. A relation between actuator dynamics, sampling time-step and limits on the flap deflection angle to guarantee the effectiveness of the adaptive controller was demonstrated by the results generated.

Friedmann, P.↗
NASA NTRS1992

Detector noise statistics in the non-linear regime

The statistical behavior of an idealized linear detector in the presence of threshold and saturation levels is examined. It is assumed that the noise is governed by the statistical fluctuations in the number of photons emitted by the source during an exposure. Since physical detectors cannot have infinite dynamic range, our model illustrates that all devices have non-linear regimes, particularly at high count rates. The primary effect is a decrease in the statistical variance about the mean signal due to a portion of the expected noise distribution being removed via clipping. Higher order statistical moments are also examined, in particular, skewness and kurtosis. In principle, the expected distortion in the detector noise characteristics can be calibrated using flatfield observations with count rates matched to the observations. For this purpose, some basic statistical methods that utilize Fourier analysis techniques are described.

Shopbell, P. L.↗
DOE OSTI2020

Operator inference for non-intrusive model reduction of systems with non-polynomial nonlinear terms

Here in this work we present a non-intrusive model reduction method to learn low-dimensional models of dynamical systems with non-polynomial nonlinear terms that are spatially local and that are given in analytic form. In contrast to state-of-the-art model reduction methods that are intrusive and thus require full knowledge of the governing equations and the operators of a full model of the discretized dynamical system, the proposed approach requires only the non-polynomial terms in analytic form and learns the rest of the dynamics from snapshots computed with a potentially black-box full-model solver. The proposed method learns operators for the linear and polynomially nonlinear dynamics via a least-squares problem, where the given non-polynomial terms are incorporated on the right-hand side. The least-squares problem is linear and thus can be solved efficiently in practice. The proposed method is demonstrated on three problems governed by partial differential equations, namely the diffusion–reaction Chafee–Infante model, a tubular reactor model for reactive flows, and a batch-chromatography model that describes a chemical separation process. The numerical results provide evidence that the proposed approach learns reduced models that achieve comparable accuracy as models constructed with state-of-the-art intrusive model reduction methods that require full knowledge of the governing equations.

42 ENGINEERING↗
NASA NTRS2003

Non-Migrating Tides, with Zonally Symmetric Component, Generated in the Mesosphere

For comparison with measurements from the TIMED satellite and coordinated ground based observations, we discuss results from our Numerical Spectral Model (NSM) that incorporates the Doppler Spread Parameterization (Hines, 1997) for small-scale gravity waves (GWs). The NSM extends from the ground into the thermosphere and 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. With emphasis on the non-migrating tides, having periods of 24 and 12 hours, we discuss our modeling results that account for the classical migrating solar excitation sources only. As reported earlier, the NSM reproduces the observed seasonal variations and in particular the large equinoctial maxima in the amplitude of the migrating diurnal tide at altitudes around 90 km. Filtering of the tide by the zonal circulation and GW momentum deposition was identified as the cause. The GWs were also shown to produce a strong non-linear interaction between the diurnal and semi-diurnal tides. Confined largely to the mesosphere, the NSM produces through dynamical interactions a relatively large contribution of non-migrating tides. A striking feature is seen in the diurnal and semi-diurnal oscillations of the zonal mean (m = 0). Eastward propagating tides are also generated for zonal wave numbers m = 1 to 4. When the NSM is run without GWs, the amplitudes for the non-migrating tides, including m = 0, are generally small. Planetary wave interaction and non-linear coupling that involves the filtering of GWs and related height integration of dynamical features are discussed as possible mechanisms for generating these non-migrating tides in the NSM. As is the case for the solar migrating tides, the non-migrating tides reveal persistent seasonal variations. Under the influence of the QBO and SAO, interannual variations are produced.

Mayr, H. G.↗
DOE OSTI2025

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Experimental verification of integrability in a Danilov-Nagaitsev lattice using machine learning

In non-linear optics, achieving integrability can enhance the dynamic aperture in storage rings. We analyze turn-by-turn phase-space data from our Danilov-Nagaitsev lattice implementation at Fermilab's Integrable Optics Test Accelerator using machine learning. AI Poincaré estimates conserved quantities from experimental data without prior knowledge of the invariant structure, showing qualitative agreement with theoretical predictions. Additionally, one of the two learned invariants exhibits comparable or better conservation compared to known theoretical expressions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Experimental Verification of Integrability in a Danilov-Nagaitsev Lattice using Machine Learning

In non-linear optics, achieving integrability can enhance the dynamic aperture in storage rings. We analyze turn-by-turn phase-space data from our Danilov-Nagaitsev lattice implementation at Fermilab's Integrable Optics Test Accelerator using machine learning. \textit{AI Poincar\'e} estimates conserved quantities from experimental data without prior knowledge of the invariant structure, showing qualitative agreement with theoretical predictions. Additionally, one of the two learned invariants exhibits comparable or better conservation compared to known theoretical expressions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗