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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Deep Koopman Neural Network for Analyzing High-Energy-Density Simulations of Electrical Wire Explosions

Megaampere-scale electrical wire experiments (EWEs) provide a platform for studying magnetohydrodynamic (MHD) instability growth in magneto-inertial fusion (MIF) devices. Even when nonlinear simulations of these experiments can digitally reproduce much of the experimentally observed instability growth, interpreting the results and understanding mode growth and evolution can be non-trivial. As a first step toward providing better interpretation of these simulation features, this work investigates the use of a deep neural network that uses Koopman operator theory to analyze the dynamics of pulsed-power-driven explosions of EWEs. This deep neural network is trained on 1-D resistive MHD simulations of EWEs. This neural network learns to transform the nonlinear data into a lower-dimensional representation where the time dynamics are linear. Layers of this neural network are shown to learn features of the simulations, including the locations of shock waves and different physical regimes of the simulation. Using the learned features, the network can compress a time state of the simulation consisting of 5120 data point into a 36-parameter lower-dimensional latent space embedding. Furthermore, these embeddings are shown to be clustered in the latent space by initial radius and time state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Development of A Multi-Spectral Pyrometry Sensor for High-Speed Transient Surface-Temperature Measurements in Combustion-Relevant Harsh Environments

Accurate and high-speed transient surface-temperature measurements of combustion devices including internal combustion (IC) engines, gas turbines, etc., provide validation targets and boundary conditions for computational fluid dynamics models, and are broadly relevant to technology advancements such as performance improvement and emissions reduction. Development and demonstration of a multi-infrared-channel pyrometry-based optical instrument for high-speed surface-temperature measurement is described. The measurement principle is based on multi-spectral radiation thermometry (MRT) and uses surface thermal radiation at four discrete spectral regions and a corresponding emissivity model to obtain surface temperature via non-linear least squares (NLLS) optimization. Rules of thumb for specifying the spectral regions and considerations to avoid interference with common combustion products are developed; the impact of these along with linear and non-linear MRT analysis are assessed as a function of temperature and signal-to-noise ratio. A multi-start method to determine the MRT-solution global optimum is described and demonstrated. The resulting multi-channel transient pyrometry instrument is described along with practical considerations including optical-alignment drift, matching intra-channel transient response, and solution-confidence indicators. The instrument demonstrated excellent >97% accuracy and >99% 2-sigma precision over the 400–800 °C range, with ~20 µs (50 kHz, equivalent to 0.2 cad at 2000 RPM IC-engine operation) transient response in the bench validation.

47 OTHER INSTRUMENTATION↗
NASA NTRS1998

A Magnetic Suspension and Excitation System for Spin Vibration Testing of Turbomachinery Blades

The Dynamic Spin Rig (DSR) is used to perform vibration tests of turbomachinery blades and components under spinning conditions in a vacuum. A heteropolar radial active magnetic bearing was integrated into the DSR to provide non-contact magnetic suspension and mechanical excitation of the rotor to induce turbomachinery blade vibrations. The magnetic bearing replaces one of the two existing conventional radial ball bearings. Prior operation of the DSR used two voice-coil type linear electromagnetic shakers which provided axial excitation of the rotor. The new magnetic suspension and excitation system has provided enhanced testing capabilities. Tests were performed at high rotational speeds for longer duration and higher vibration amplitudes. Some characteristics of the system include magnetic bearing stiffness values up to 60,000 lb./in., closed loop control bandwidth around 500 Hz, and multi-directional radial excitation of the rotor. This paper reports on the implementation and operation of this system and presents some test results using this system.

Johnson, Dexter↗
NASA NTRS2012

Nonlinear Reduced-Order Simulation Using An Experimentally Guided Modal Basis

A procedure is developed for using nonlinear experimental response data to guide the modal basis selection in a nonlinear reduced-order simulation. The procedure entails using nonlinear acceleration response data to first identify proper orthogonal modes. Special consideration is given to cases in which some of the desired response data is unavailable. Bases consisting of linear normal modes are then selected to best represent the experimentally determined transverse proper orthogonal modes and either experimentally determined inplane proper orthogonal modes or the special case of numerically computed in-plane companions. The bases are subsequently used in nonlinear modal reduction and dynamic response simulations. The experimental data used in this work is simulated to allow some practical considerations, such as the availability of in-plane response data and non-idealized test conditions, to be explored. Comparisons of the nonlinear reduced-order simulations are made with the surrogate experimental data to demonstrate the effectiveness of the approach.

Rizzi, Stephen A.↗
NASA NTRS1993

Evolution in vitro: analysis of a lineage of ribozymes

Background: Catalytic RNAs, or ribozymes, possessing both a genotype and a phenotype, are ideal molecules for evolution experiments in vitro. A large, heterogeneous pool of RNAs can be subjected to multiple rounds of selection, amplification and mutation, leading to the development of variants that have some desired phenotype. Such experiments allow the investigator to correlate specific genetic changes with quantifiable alterations of the catalytic properties of the RNA. In addition, patterns of evolutionary change can be discerned through a detailed examination of the genotypic composition of the evolving RNA population. Results: Beginning with a pool of 10(13) variants of the Tetrahymena ribozyme, we carried out in vitro evolution experiments that led to the generation of ribozymes with the ability to cleave an RNA substrate in the presence of Ca2+ ions, an activity that does not exist for the wild-type molecule. Over the course of 12 generations, a seven-error variant emerged that has substantial Ca(2+)-dependent RNA-cleavage activity. Advantageous mutations increased in frequency in the population according to three distinct dynamics--logarithmic, linear and transient. Through a comparative analysis of 31 individual variants, we infer how certain mutations influence the catalytic properties of the ribozyme. Conclusions: In vitro evolution experiments make it possible to elucidate important aspects of both evolutionary biology and structural biochemistry on a reasonable short time scale.

Non-NASA Center↗
NASA NTRS1996

Computational Analysis of Static and Dynamic Behaviour of Magnetic Suspensions and Magnetic Bearings

Static modelling of magnetic bearings is often carried out using magnetic circuit theory. This theory cannot easily include nonlinear effects such as magnetic saturation or the fringing of flux in air-gaps. Modern computational tools are able to accurately model complex magnetic bearing geometries, provided some care is exercised. In magnetic suspension applications, the magnetic fields are highly three-dimensional and require computational tools for the solution of most problems of interest. The dynamics of a magnetic bearing or magnetic suspension system can be strongly affected by eddy currents. Eddy currents are present whenever a time-varying magnetic flux penetrates a conducting medium. The direction of flow of the eddy current is such as to reduce the rate-of-change of flux. Analytic solutions for eddy currents are available for some simplified geometries, but complex geometries must be solved by computation. It is only in recent years that such computations have been considered truly practical. At NASA Langley Research Center, state-of-the-art finite-element computer codes, 'OPERA', 'TOSCA' and 'ELEKTRA' have recently been installed and applied to the magnetostatic and eddy current problems. This paper reviews results of theoretical analyses which suggest general forms of mathematical models for eddy currents, together with computational results. A simplified circuit-based eddy current model proposed appears to predict the observed trends in the case of large eddy current circuits in conducting non-magnetic material. A much more difficult case is seen to be that of eddy currents in magnetic material, or in non-magnetic material at higher frequencies, due to the lower skin depths. Even here, the dissipative behavior has been shown to yield at least somewhat to linear modelling. Magnetostatic and eddy current computations have been carried out relating to the Annular Suspension and Pointing System, a prototype for a space payload pointing and vibration isolation system, where the magnetic actuator geometry resembles a conventional magnetic bearing. Magnetostatic computations provide estimates of flux density within airgaps and the iron core material, fringing at the pole faces and the net force generated. Eddy current computations provide coil inductance, power dissipation and the phase lag in the magnetic field, all as functions of excitation frequency. Here, the dynamics of the magnetic bearings, notably the rise time of forces with changing currents, are found to be very strongly affected by eddy currents, even at quite low frequencies. Results are also compared to experimental measurements of the performance of a large-gap magnetic suspension system, the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Eddy current effects are again shown to significantly affect the dynamics of the system. Some consideration is given to the ease and accuracy of computation, specifically relating to OPERA/TOSCA/ELEKTRA.

Britcher, Colin P.↗
NASA NTRS2012

Low-Dissipation Advection Schemes Designed for Large Eddy Simulations of Hypersonic Propulsion Systems

The 2nd-order upwind inviscid flux scheme implemented in the multi-block, structured grid, cell centered, finite volume, high-speed reacting flow code VULCAN has been modified to reduce numerical dissipation. This modification was motivated by the desire to improve the codes ability to perform large eddy simulations. The reduction in dissipation was accomplished through a hybridization of non-dissipative and dissipative discontinuity-capturing advection schemes that reduces numerical dissipation while maintaining the ability to capture shocks. A methodology for constructing hybrid-advection schemes that blends nondissipative fluxes consisting of linear combinations of divergence and product rule forms discretized using 4th-order symmetric operators, with dissipative, 3rd or 4th-order reconstruction based upwind flux schemes was developed and implemented. A series of benchmark problems with increasing spatial and fluid dynamical complexity were utilized to examine the ability of the candidate schemes to resolve and propagate structures typical of turbulent flow, their discontinuity capturing capability and their robustness. A realistic geometry typical of a high-speed propulsion system flowpath was computed using the most promising of the examined schemes and was compared with available experimental data to demonstrate simulation fidelity.

White, Jeffrey A.↗
DOE OSTI2023

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Extended nonergodic regime and spin subdiffusion in disordered SU(2)-symmetric Floquet systems

In this work, we explore thermalization and quantum dynamics in a one-dimensional disordered SU(2)-symmetric Floquet model, where a many-body localized phase is prohibited by the non-Abelian symmetry. Despite the absence of localization, we find an extended nonergodic regime at strong disorder where the system exhibits nonthermal behaviors. In the strong disorder regime, the level spacing statistics exhibit neither a Wigner-Dyson nor a Poisson distribution, and the spectral form factor does not show a linear-in-time growth at early times characteristic of random matrix theory. The average entanglement entropy of the Floquet eigenstates is subthermal, although violating an area-law scaling with system sizes. We further compute the expectation value of local observables and find strong deviations from the eigenstate thermalization hypothesis. The infinite-temperature spin autocorrelation function decays at long times as $t^{–β}$ with $\textit{β}$ < 0.5, indicating subdiffusive transport at strong disorders.

36 MATERIALS SCIENCE↗
NASA NTRS1989

A simple nonlinear model for the return to isotropy in turbulence

A quadratic nonlinear generalization of the linear Rotta model for the slow pressure-strain correlation of turbulence is developed. The model is shown to satisfy realizability and to give rise to no stable non-trivial equilibrium solutions for the anisotropy tensor in the case of vanishing mean velocity gradients. The absence of stable non-trivial equilibrium solutions is a necessary condition to ensure that the model predicts a return to isotropy for all relaxational turbulent flows. Both the phase space dynamics and the temporal behavior of the model are examined and compared against experimental data for the return to isotropy problem. It is demonstrated that the quadratic model successfully captures the experimental trends which clearly exhibit nonlinear behavior. Direct comparisons are also made with the predictions of the Rotta model and the Lumley model.

Sarkar, Sutanu↗
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↗
DOE OSTI2025

Analysis of soot formation from aviation fuels in laminar counterflow flames

Combustion emissions from aviation contribute to the formation of condensation trail (contrail) that can lead to the formation of anthropogenic cirrus clouds. Ice particles that form contrails are observed to have a linear correlation with soot particle number density. Synthetic aviation fuels (SAFs) offer a promising route to mitigate the production of soot particles while also increasing energy security. Although studies have focused on combustion and spray behavior, the detailed investigation of soot formation processes for different jet fuels and their impact on models for computational fluid dynamics (CFD) applications is not well understood. Moreover, experimental measurements of soot for canonical flames using Synthetic aviation fuels (SAF) for model validation remain scarce. To address this, we use employed the Lawrence Livermore National Laboratory (LLNL) detailed soot model based on the discrete sectional method. Additionally, we develop two reduced chemical mechanisms for Jet-A and Alcohol-to-Jet (C1) that are suitable for turbulent flame simulations and couple them with the Hybrid Method of Moments (HMOM). The detailed and reduced model frameworks are validated against experimental measurements of soot volume fraction (ƒ ν ) from a counterflow burner experiment previously reported in the literature. Given the good agreement between modeling results and experimental measurements for the (1) spatial distribution of ƒ ν and (2) the non-linear variation of peak ƒ ν with strain rate, we further investigate the modeled sub-processes (nucleation, condensation, surface growth, and oxidation) using the LLNL model to analyze the assumptions in the reduced model framework. Furthermore, the results indicate a significant contribution from resonant radicals to the surface growth of soot particles, which are not accounted for in the current implementation of HMOM and could help reconcile soot predictions by the reduced model with observations.

Counterflow↗
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↗
DOE OSTI2024

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2021

On applications of quantum computing to plasma simulations

Quantum computing is gaining increased attention as a potential way to speed up simulations of physical systems, and it is also of interest to apply it to simulations of classical plasmas. However, quantum information science is traditionally aimed at modeling linear Hamiltonian systems of a particular form that is found in quantum mechanics, so extending the existing results to plasma applications remains a challenge. Here, we report a preliminary exploration of the long-term opportunities and likely obstacles in this area. First, we show that many plasma-wave problems are naturally representable in a quantumlike form and thus are naturally fit for quantum computers. Second, we consider more general plasma problems that include non-Hermitian dynamics (instabilities, irreversible dissipation) and nonlinearities. We show that by extending the configuration space, such systems can also be represented in a quantumlike form and thus can be simulated with quantum computers too, albeit that requires more computational resources compared to the first case. Third, we outline potential applications of hybrid quantum–classical computers, which include analysis of global eigenmodes and also an alternative approach to nonlinear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗