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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Acceleration of Power System Dynamic Simulations Using a Deep Equilibrium Layer and Neural ODE Surrogate

The dominant paradigm for power system dynamic simulation is to build system-level simulations by combining physics-based models of individual components. The sheer size of the system along with the rapid integration of inverter-based resources exacerbates the computational burden of running time domain simulations. Here, in this paper, we propose a data-driven surrogate model based on implicit machine learningspecifically deep equilibrium layers and neural ordinary differential equationsto learn a reduced order model of a portion of the full underlying system. The data-driven surrogate achieves similar accuracy and reduction in simulation time compared to a physics-based surrogate, without the constraint of requiring detailed knowledge of the underlying dynamic models. This work also establishes key requirements needed to integrate the surrogate into existing simulation workflows; the proposed surrogate is initialized to a steady state operating point that matches the power flow solution by design.

Neural ordinary differential equations↗
NASA NTRS2011

Dynamic Monitoring of Cleanroom Fallout Using an Air Particle Counter

The particle fallout limitations and periodic allocations for the James Webb Space Telescope are very stringent. Standard prediction methods are complicated by non-linearity and monitoring methods that are insufficiently responsive. A method for dynamically predicting the particle fallout in a cleanroom using air particle counter data was determined by numerical correlation. This method provides a simple linear correlation to both time and air quality, which can be monitored in real time. The summation of effects provides the program better understanding of the cleanliness and assists in the planning of future activities. Definition of fallout rates within a cleanroom during assembly and integration of contamination-sensitive hardware, such as the James Webb Space Telescope, is essential for budgeting purposes. Balancing the activity levels for assembly and test with the particle accumulation rate is paramount. The current approach to predicting particle fallout in a cleanroom assumes a constant air quality based on the rated class of a cleanroom, with adjustments for projected work or exposure times. Actual cleanroom class can also depend on the number of personnel present and the type of activities. A linear correlation of air quality and normalized particle fallout was determined numerically. An air particle counter (standard cleanroom equipment) can be used to monitor the air quality on a real-time basis and determine the "class" of the cleanroom (per FED-STD-209 or ISO-14644). The correlation function provides an area coverage coefficient per class-hour of exposure. The prediction of particle accumulations provides scheduling inputs for activity levels and cleanroom class requirements.

Perry, Radford↗
DOE OSTI2020

A State-Space Model of an Inverter-Based Microgrid for Multivariable Feedback Control Analysis and Design

In this work, a synchronous model for grid-connected and islanded microgrids is presented. The grid-connected model is based on the premise that the reference frame is synchronized with the AC bus. The quadrature component of the AC bus voltage can be cancelled, which allows to express output power as a linear equation for nominal values in the AC bus amplitude voltage. The model for the islanded microgrid is developed by integrating all the inverter dynamics using a state-space model for the load currents. This model is presented in a comprehensive way such that it could be scalable to any number of inverter-based generators using inductor–capacitor–inductor (LCL) output filters. The use of these models allows designers to assess microgrid stability and robustness using modern control methods such as eigenvalue analysis and singular value diagrams. Both models were tested and validated in an experimental setup to demonstrate their accuracy in describing microgrid dynamics. In addition, three scenarios are presented: non-controlled model, Linear-Quadratic Integrator (LQI) power control, and Power-Voltage (PQ/Vdq) droop–boost controller. Experimental results demonstrate the effectiveness of the control strategies and the accuracy of the models to describe microgrid dynamics.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2004

Continuation Methods for Qualitative Analysis of Aircraft Dynamics

A class of numerical methods for constructing bifurcation curves for systems of coupled, non-linear ordinary differential equations is presented. Foundations are discussed, and several variations are outlined along with their respective capabilities. Appropriate background material from dynamical systems theory is presented.

Cummings, Peter A.↗
DOE OSTI2025

The Water Table Model (WTM) (v2.0.1): coupled groundwater and dynamic lake modelling

Abstract. Ice-free land comprises 26 % of the Earth's surface and holds liquid water that delineates ecosystems, affects global geochemical cycling, and modulates sea levels. However, we currently lack the capacity to simulate and predict these terrestrial water changes across the full range of relevant spatial (watershed to global) and temporal (monthly to millennial) scales. To address this knowledge gap, we present the Water Table Model (WTM), which integrates coupled components to compute dynamic lake and groundwater levels. The groundwater component solves the 2D horizontal groundwater flow equation using non-linear equation solvers from the C++ PETSc (Portable, Extensible Toolkit for Scientific Computation) library. The dynamic lake component makes use of the Fill–Spill–Merge (FSM) algorithm to move surface water into lakes, where it may evaporate or affect groundwater flow. In a proof-of-concept application, we demonstrate the continental-scale capabilities of the WTM by simulating the steady-state climate-driven water table for the present day and the Last Glacial Maximum (LGM; 21 000 calendar years before present) across the North American continent. During the LGM, North America stored an additional 14.98 cm of sea-level equivalent (SLE) in lakes and groundwater compared to the climate-driven present-day scenario. We compare the present-day result to other simulations and real-world data. Open-source code for the WTM is available on GitHub and Zenodo.

Callaghan, Kerry L. (ORCID:0000000226740838)↗
NASA NTRS1998

Radio to Gamma-Ray Emission from Shell-Type Supernova Remnants: Predictions from Non-Linear Shock Acceleration Models

Supernova remnants (SNRs) are widely believed to be the principal source of galactic cosmic rays, produced by diffusive shock acceleration in the environs of the remnant's expanding blast wave. Such energetic particles can produce gamma-rays and lower energy photons via interactions with the ambient plasma. The recently reported observation of TeV gamma-rays from SN1006 by the CANGAROO Collaboration, combined with the fact that several unidentified EGRET sources have been associated with known radio/optical/X-ray-emitting remnants, provides powerful motivation for studying gamma-ray emission from SNRs. In this paper, we present results from a Monte Carlo simulation of non-linear shock structure and acceleration coupled with photon emission in shell-like SNRs. These non-linearities are a by-product of the dynamical influence of the accelerated cosmic rays on the shocked plasma and result in distributions of cosmic rays which deviate from pure power-laws. Such deviations are crucial to acceleration efficiency considerations and impact photon intensities and spectral shapes at all energies, producing GeV/TeV intensity ratios that are quite different from test particle predictions.

Baring, Matthew G.↗
NASA NTRS2004

Formations Near the Libration Points: Design Strategies Using Natural and Non-Natural Arcs

Space based observatory and interferometry missions, such as Terrestrial Planet Finder (TPF), Stellar Imager, and MAXIM, have sparked great interest in multi-spacecraft formation flight in the vicinity of the Sun-Earth/Moon (SEM) libration points. The initial phase of this research considered the formation keeping problem from the perspective of continuous control as applied to non-natural formations. In the present study, closer inspection of the flow, corresponding to the stable and center manifolds near the reference orbit, reveals some interesting natural relative motions as well as some discrete control strategies for deployment. A hybrid control strategy is also employed that combines both the natural formation dynamics with non-natural motions via input feedback linearization techniques.

Howell, K. C.↗
NASA NTRS2013

A Fiber-Optic Aircraft Lightning Current Measurement Sensor

A fiber-optic current sensor based on the Faraday Effect is developed for aircraft installations. It can measure total lightning current amplitudes and waveforms, including continuing current. Additional benefits include being small, lightweight, non-conducting, safe from electromagnetic interference, and free of hysteresis and saturation. The Faraday Effect causes light polarization to rotate in presence of magnetic field in the direction of light propagation. Measuring the total induced light polarization change yields the total current enclosed. The system operates at 1310nm laser wavelength and can measure approximately 300 A - 300 kA, a 60 dB range. A reflective polarimetric scheme is used, where the light polarization change is measured after a round-trip propagation through the fiber. A two-detector setup measures the two orthogonal polarizations for noise subtraction and improved dynamic range. The current response curve is non-linear and requires a simple spline-fit correction. Effects of high current were achieved in laboratory using combinations of multiple fiber and wire loops. Good result comparisons against reference sensors were achieved up to 300 kA. Accurate measurements on a simulated aircraft fuselage and an internal structure illustrate capabilities that maybe difficult with traditional sensors. Also tested at a commercial lightning test facility from 20 kA to 200 kA, accuracy within 3-10% was achieved even with non-optimum setups.

Nguyen, Truong X.↗
NASA NTRS2005

A Reduced-frequency Approach for Calculating Dynamic Derivatives

Computational Fluid Dynamics (CFD) is increasingly being used to both augment and create an aerodynamic performance database for aircraft configurations. This aerodynamic database contains the response of the aircraft to varying flight conditions and control surface deflections. The current work presents a novel method for calculating dynamic stability derivatives which reduces the computational cost over traditional unsteady CFD approaches by an order of magnitude, while still being applicable to arbitrarily complex geometries over a wide range of flow regimes. The primary thesis of this work is that the response to a forced motion can often be represented with a small, predictable number of frequency components without loss of accuracy. By resolving only those frequencies of interest, the computational effort is significantly reduced so that the routine calculation of dynamic derivatives becomes practical. The current implementation uses this same non-linear, frequency-domain approach and extends the application to the 3-D Euler equations. The current work uses a Cartesian, embedded-boundary method to automate the generation of dynamic stability derivatives.

Murman, Scott M.↗
DOE OSTI2020

The effect of non-linear mutual friction on pulsar glitch sizes and rise times

ABSTRACT Observations of pulsar glitches have the potential to provide constraints on the dynamics of the high density interior of neutron stars. However, to do so, realistic glitch models must be constructed and compared to the data. We take a step towards this goal by testing non-linear models for the mutual friction force, which is responsible for the exchange of angular momentum between the neutron superfluid and the observable normal component in a glitch. In particular, we consider a non-linear dependence of the drag force on the relative velocity between superfluid vortices and the normal component, in which the contributions of both kelvin and phonon excitations are included. This non-linear model produces qualitatively new features, and is able to reproduce the observed bimodal distribution of glitch sizes in the pulsar population. The model also suggests that the differences in size distributions in individual pulsars may be due to the glitches being triggered in regions with different pinning strengths, as stronger pinning leads to higher vortex velocities and a qualitatively different mutual friction coupling with respect to the weak pinning case. Glitches in pulsars that appear to glitch quasi-periodically with similar sizes may thus be due to the same mechanisms as smaller events in pulsars that have no preferred glitch size, but simply originate in stronger pinning regions, possibly in the core of the star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Uncertainty Error Modeling for Non-Linear State Estimation With Unsynchronized SCADA and µPMU Measurements

Distribution systems of the future smart grid require enhancements to the reliability of distribution system state estimation (DSSE) in the face of low measurement redundancy, unsynchronized measurements, and dynamic load profiles. Micro phasor measurement units (µPMUs) facilitate co-synchronized measurements with high granularity, albeit at an often prohibitively expensive installation cost. Supervisory control and data acquisition (SCADA) measurements can supplement µPMU data, although they are received at a slower sampling rate. Further complicating matters is the uncertainty associated with load dynamics and unsynchronized measurements–not only are the SCADA and µPMU measurements not synchronized with each other, but the SCADA measurements themselves are received at different time intervals with respect to one another. This paper proposes a non-linear state estimation framework which models dynamic load uncertainty error by updating the variances of the unsynchronized measurements, leading to a time-varying system of weights in the weighted least squares state estimator. Case studies are performed on the 33-Bus Distribution System in MATPOWER, using Ornstein–Uhlenbeck stochastic processes to simulate dynamic load conditions.

Cooper, Austin↗
DOE OSTI2020

Extrapolation of thermal conductivity in non-equilibrium molecular dynamics simulations to bulk scale

Predictions of the bulk scale thermal conductivity of solids using non-equilibrium molecular dynamics (MD) simulations have relied on the linear extrapolation of the thermal resistivity versus the reciprocal of the system length in the simulations. Several studies have reported deviation of the extrapolation from linearity near the micro-scale, raising a concern of its applicability to large systems. To investigate this issue, present work conducted extensive MD simulations of silicon with two different potentials (EDIP and Tersoff-II) for unprecedented length scales up to 10.3 μm and simulation times up to 530 ns. For large systems ≥0.35 μm in size the non-linearity of the extrapolation of the reciprocal of the thermal conductivity is mostly due to ignoring the dependence of the thermal conductivity on temperature. To account for such dependence, the present analysis fixes the temperature range for determining the gradient for calculating the thermal conductivity values. However, short systems ≤0.23 μm in size show significant non-linearity in the calculated thermal conductivity values using a temperature window of 500 ± 10 K from the simulations results with the EDIP potential. Since these system sizes are shorter than the mean phonon free path in EDIP (~0.22 μm), the nonlinearity may be attributed to phonon transport. For the MD simulations with the Tersoff-II potential there is no significant non-linearity in the calculated thermal conductivity values for systems ranging in size from 0.05 to 5.4 μm.

42 ENGINEERING↗
DOE OSTI2021

Effects of interaction strength of associating groups on linear and star polymer dynamics

A small number of associating groups incorporated onto a polymer backbone have dramatic effects on the mobility and viscoelastic response of the macromolecules in melts. These associating groups assemble, driving the formation of clusters, whose lifetime affects the properties of the polymers. In this study, we probe the effects of the interaction strength on the structure and dynamics of two topologies, linear and star polymer melts, and further investigate blends of associative and non-associating polymers using molecular dynamics simulations. Polymer chains of approximately one entanglement length are described by a bead–spring model, and the associating groups are incorporated in the form of interacting beads with an interaction strength between them that is varied from 1 to 20 k B T. We find that, for all melts and blends, interaction of a few k B T between the associating groups drives cluster formation, where the size of the clusters increases with increasing interaction strength. These clusters act as physical crosslinkers, which slow the chain mobility. Blends of chains with and without associating groups macroscopically phase separate for interaction strength between the associating groups of a few k B T and above. For weakly interacting associating groups, the static structure function S(q) is well fit by functional form predicted by the random phase approximation where a clear deviation occurs as phase segregation takes place, providing a quantitative assessment of phase segregation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamic Aperture of the PIP-II Accumulator Ring (PAR)

The Proton Improvement Plan-II (PIP-II) will include a superconducting linear accelerator which will accelerate H- ions to 800MeV, at an intensity never seen before. A proposed PIP-II Accumulator Ring (PAR) can reduce the injection losses to the booster, enhance the ramp up of the Long Baseline Neutrino Facility (LBNF) and support experiments in Dark Sector Physics. When the beam is circulating in PAR at its full intensity, it will then be transferred to the Booster via single turn injection. Strong PAR focusing results in larger phase advances, shorter focal lengths and higher machine tunes, which creates robust beam envelope control. Particles will be able to survive modest excursions from the elliptical phase-space trajectories of linear optics, and any loss of particles would be due to reduced dynamic aperture. The values for tunes and chromaticities are flexible and will continue to evolve as the lattice matures. Since there will be non-linear terms impacting the particles’ trajectories, tracking will include realistic magnetic field errors and installation alignment errors to determine the dynamic aperture as a function of the machine tunes. This will allow us to determine the optimum operating point in phase-space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2000

Determination of anionic surfactants during wastewater recycling process by ion pair chromatography with suppressed conductivity detection

A direct approach utilizing ion pairing reversed-phase chromatography coupled with suppressed conductivity detection was developed to monitor biodegradation of anionic surfactants during wastewater recycling through hydroponic plant growth systems and fixed-film bioreactors. Samples of hydroponic nutrient solution and bioreactor effluent with high concentrations (up to 120 mS electrical conductance) of inorganic ions can be analyzed without pretreatment or interference. The presence of non-ionic surfactants did not significantly affect the analysis. Dynamic linear ranges for tested surfactants [Igepon TC-42, ammonium lauryl sulfate, sodium laureth sulfate and sodium alkyl (C10-C16) ether sulfate] were 2 to approximately 500, 1 to approximately 500, 2.5 to approximately 550 and 3.0 to approximately 630 microg/ml, respectively.

NASA Program Advanced Human Support Technology↗
NASA NTRS2005

Dynamic Weakening (Extinction) of Simple Hydrocarbon-air Counterflow Diffusion Flames by Oscillatory Inflows

This study of laminar non-premixed HC-air flames used an Oscillatory-input Opposed Jet Burner (OOJB) system developed from a previously well-characterized 7.2-mm Pyrex-nozzle OJB system. Over 600 dynamic Flame Strength (FS) measurements were obtained on unanchored (free-floating) laminar Counterflow Diffusion Flames (CFDF's). Flames were stabilized using plug inflows having steady-plus-sinusoidal axial velocities of varied magnitude, frequency, f, up to 1600 Hz, and phase angle from 0 (most data) to 360 degrees. Dynamic FS is defined as the maximum average air input velocity (U(sub air), at nozzle exit) a CFDF can sustain before strain-induced extinction occurs due to prescribed oscillatory peak-to-peak velocity inputs superimposed on steady inputs. Initially, dynamic flame extinction data were obtained at low f, and were supported by 25-120 Hz Hot-Wire cold-flow velocity data at nozzle exits. Later, expanded extinction data were supported by 4-1600 Hz Probe Microphone (PM) pk/pk P data at nozzle exits. The PM data were first obtained without flows, and later with cold stagnating flows, which better represent speaker-diaphragm dynamics during runs. The PM approach enabled characterizations of Dynamic Flame Weakening (DFW) of CFDF's from 8 to 1600 Hz. DFW was defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = - 100 d(U(sub air) / U(sub air),0Hz) / d(pkpk P). The linear normalization with respect to acoustic pressure magnitude (and steady state (SS) FS) led to a DFW unaffected by strong internal resonances. For the C2H4/N2-air system, from 8 to 20 Hz, DFW is constant at 8.52 plus or minus 0.20 (% weakening)/Pa. This reflects a quasi-steady flame response to an acoustically induced dU(sub air)/dP. Also, it is surprisingly independent of C2H4/N2 mole fraction due to normalization by SS FS. From 20 to approximately 150 Hz, the C2H4/N2 air-flames weakened progressively less, with an inflection at approximately 70 Hz, and became asymptotically insensitive (DFW approximately 0) at approximately 300 Hz, which continued to 1600 Hz. The DFW of CH4-air flames followed a similar pattern, but showed much greater weakening than C2H4/N2-air flames; i.e., the quasi-steady DFW (8 to approximately 15 Hz) was 44.3 %/Pa, or approximately 5x larger, even though the 0 Hz (SS) FS was only 3.0 x smaller. The quasi-steady DFW's of C3H8-air and C2H6-air were intermediate at 34.8 and 20.9 %Pa, respectively. The DFW profiles of all four fuels, at various frequencies, correlated well but non-linearly with respective SS FS's. Notably, the DFW profile for C3H8 air fell more rapidly in the range greater than 15 to 60 Hz, compared with the 1- and 2-carbon fuels. This may indicate a shift in chemical kinetics, and/or O2 transport to a flame that moved closer to the fuel-side. In conclusion, Dynamic Flame Weakening limits appear significant and unique for each fuel, and correlate closely, but non-linearly, with Steady-State Flame Strengths at any given frequency. For reasons unknown, the dynamic flames didn't weaken more at intermediate frequencies (e.g., at 20-50 Hz) than they did at low frequencies (less than 15 Hz), where quasi-steady weakening appears to dominate. Quasi-steady flame weakening ostensibly represents a transient input strain rate maximum that just exceeds the steady-state strain-rate-limited extinction limit for a few cycles. Clearly, further detailed mechanistic understanding is needed in the fall-off region.

Pellett, G.↗
NASA NTRS2012

Active Combustion Control for Aircraft Gas-Turbine Engines-Experimental Results for an Advanced, Low-Emissions Combustor Prototype

Lean combustion concepts for aircraft engine combustors are prone to combustion instabilities. Mitigation of instabilities is an enabling technology for these low-emissions combustors. NASA Glenn Research Center s prior activity has demonstrated active control to suppress a high-frequency combustion instability in a combustor rig designed to emulate an actual aircraft engine instability experience with a conventional, rich-front-end combustor. The current effort is developing further understanding of the problem specifically as applied to future lean-burning, very low-emissions combustors. A prototype advanced, low-emissions aircraft engine combustor with a combustion instability has been identified and previous work has characterized the dynamic behavior of that combustor prototype. The combustor exhibits thermoacoustic instabilities that are related to increasing fuel flow and that potentially prevent full-power operation. A simplified, non-linear oscillator model and a more physics-based sectored 1-D dynamic model have been developed to capture the combustor prototype s instability behavior. Utilizing these models, the NASA Adaptive Sliding Phasor Average Control (ASPAC) instability control method has been updated for the low-emissions combustor prototype. Active combustion instability suppression using the ASPAC control method has been demonstrated experimentally with this combustor prototype in a NASA combustion test cell operating at engine pressures, temperatures, and flows. A high-frequency fuel valve was utilized to perturb the combustor fuel flow. Successful instability suppression was shown using a dynamic pressure sensor in the combustor for controller feedback. Instability control was also shown with a pressure feedback sensor in the lower temperature region upstream of the combustor. It was also demonstrated that the controller can prevent the instability from occurring while combustor operation was transitioning from a stable, low-power condition to a normally unstable high-power condition, thus enabling the high-power condition.

DeLaat, John C.↗
NASA NTRS2021

Truncated ARQ Statistical Link Analysis for Dynamic Links

The future deep space links are migrating towards higher frequency bands such as Ka band and optical. These links are susceptible to non Gaussian and non linear effects such as atmospheric turbulence, scintillation, antenna mis-pointing, jitter, etc. These dynamic links thus will experience various degrees of fading loss, and some of these link disruptions cannot be effectively mitigated by forward error correction coding and/or interleaving. One effective way to ensure reliable communication is by using Automatic Repeat Request (ARQ) protocol, where the receiver acknowledges to the transmitter whether or not a data unit is successfully received. If a data unit is not successfully received (such as after a pre-set time-out), the transmitter would then re-transmit the lost data unit to the receiver. In a previous paper, we derived a statistical link analysis method of finding the optimal operating Signal-to-Noise Ratio (SNR) and estimating the latency of an ARQ scheme. In a more recent paper, we demonstrated the above method using the SNR distribution constructed from the Ka-band (32 GHz) flight data. To simplify the discussion, we considered the academic approach that the ARQ scheme allows for an infinite number of retransmissions. In this paper, we consider the more practical case of a truncated ARQ scheme, where there is a limit on the number of retransmissions. We derive the error probability, the optimal SNR setting, and the latency statistics of the correctly received frames of the truncated ARQ schemes. We first discuss the truncated ARQ link analysis principles using the Gaussian assumption for SNR distribution with a large variance. Next, we demonstrate the statistical truncated ARQ link analysis using the SNR distribution constructed from the Ka-band flight data. The results in this paper can be applied in the design of reliable communication systems such as the Consultative Committee for Space Data System (CCSDS) File Transfer Protocol (CFTP) and the Delay Tolerant Network (DTN).

Morabito, David↗