Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “classic nonlinearity”

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

Mid-Field Sonic Boom Extrapolation Methodology

In the design cycle of low-boom airplanes, sonic boom prediction must be accurate and efficient. The classical linear method, Whitham's F-function theory, has been widely applied to predict sonic boom signatures. However, linear theory fails to capture the nonlinear effects created by large civil transport. Computational fluid dynamics (CFD) has been used successfully to predict sonic boom signals at the near and mid fields. Nevertheless, it is computationally expansive in airplane design runs. In the present study, the method of characteristics is used to predict sonic boom signals in an efficient fashion. The governing equations are the axisymmetric Euler's equations with constant enthalpy. Since the method solves Euler's equations, it captures more nonlinear effects than the classical Whitham's F-function technique. Furthermore, the method of characteristics is an efficient marching scheme for initial value problems. In this study, we will first review the current CFD extrapolation technique and the work previously done in sonic boom extrapolation. Then, we will introduce the governing equations and the method of characteristics. Finally, we will show that the present method yields the same accurate results as previous CFD techniques, but with higher efficiency.

Cheung, Samson↗
NASA NTRS1993

Control of nonlinear systems with applications to constrained robots and spacecraft attitude stabilization

This thesis is organized in two parts. In Part 1, control systems described by a class of nonlinear differential and algebraic equations are introduced. A procedure for local stabilization based on a local state realization is developed. An alternative approach to local stabilization is developed based on a classical linearization of the nonlinear differential-algebraic equations. A theoretical framework is established for solving a tracking problem associated with the differential-algebraic system. First, a simple procedure is developed for the design of a feedback control law which ensures, at least locally, that the tracking error in the closed loop system lies within any given bound if the reference inputs are sufficiently slowly varying. Next, by imposing additional assumptions, a procedure is developed for the design of a feedback control law which ensures that the tracking error in the closed loop system approaches zero exponentially for reference inputs which are not necessarily slowly varying. The control design methodologies are used for simultaneous force and position control in constrained robot systems. The differential-algebraic equations are shown to characterize the slow dynamics of a certain nonlinear control system in nonstandard singularly perturbed form. In Part 2, the attitude stabilization (reorientation) of a rigid spacecraft using only two control torques is considered. First, the case of momentum wheel actuators is considered. The complete spacecraft dynamics are not controllable. However, the spacecraft dynamics are small time locally controllable in a reduced sense. The reduced spacecraft dynamics cannot be asymptotically stabilized using continuous feedback, but a discontinuous feedback control strategy is constructed. Next, the case of gas jet actuators is considered. If the uncontrolled principal axis is not an axis of symmetry, the complete spacecraft dynamics are small time locally controllable. However, the spacecraft attitude cannot be asymptotically stabilized using continuous feedback, but a discontinuous stabilizing feedback control strategy is constructed. If the uncontrolled principal axis is an axis of symmetry, the complete spacecraft dynamics cannot be stabilized. However, the spacecraft dynamics are small time locally controllable in a reduced sense. The reduced spacecraft dynamics cannot be asymptotically stabilized using continuous feedback, but again a discontinuous feedback control strategy is constructed.

Krishnan, Hariharan↗
NASA NTRS1988

Effects of nonlinear damping in flexible space structures

The classical Krylov-Bogoliubov "averaging" technique is used to study a class of nonlinear damping models, for which the damping force is proportional to the product of positive integer or fractional power of absolute values of displacement and that of velocity. The results are compared with linear viscous damping models. The amplitude decrement of free vibration for a single mode system with nonlinear models depends not only on damping ratio, but also on the initial amplitude, the time to measure the response, frequency of the system, and the powers of displacement and velocity. For the distributed system, the action of nonlinear damping is found to reduce energy of the system as well as to pass energy to higher modes. Experimental evidence such as in Spacecraft Control Laboratory Experiment seems to support the need for nonlinear models.

Hu, Anren↗
NASA NTRS2013

Planetary Protection Bioburden Analysis Program

This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware. This work was done by David Shores and Sharon P. Goza of Johnson Space Center; Cheyenne McKeegan, Rick Easley, Janet Way, and Shonn Everett of MEI Technologies; Mark Manning of PTI; and Mark Guerra, Ray Kraesig, and William Leu of Tietronix Software, Inc. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24211-1 Spitzer Telemetry Processing System NASA's Jet Propulsion Laboratory, Pasadena, California The Spitzer Telemetry Processing System (SirtfTlmProc) was designed to address objectives of JPL's Multi-mission Image Processing Lab (MIPL) in processing spacecraft telemetry and distributing the resulting data to the science community. To minimize costs and maximize operability, the software design focused on automated error recovery, performance, and information management. The system processes telemetry from the Spitzer spacecraft and delivers Level 0 products to the Spitzer Science Center. SirtfTlmProc is a unique system with automated error notification and recovery, with a real-time continuous service that can go quiescent after periods of inactivity. The software can process 2 GB of telemetry and deliver Level 0 science products to the end user in four hours. It provides analysis tools so the operator can manage the system and troubleshoot problems. It automates telemetry processing in order to reduce staffing costs. This work was done by Alice Stanboli, Elmain M. Martinez, and James M. McAuley of Caltech for NASA's Jet Propulsion Laboratory. For more information, contact iaoffice @jpl.nasa.gov. This software is available for commercial licensing. Please contact Dan Broderick at Daniel.F. Broderick@jpl.nasa.gov. Refer to NPO-47803. NASA Tech Briefs, September 2013 29 This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known. This work was done by Robert Clark, Jr., Paul Cotter, and Constantine Michalopoulos of The Boeing Company for Johnson Space Center. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24988-1 Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2) Lyndon B. Johnson Space Center, Houston, Texas (3) the MSL statistics for only the accountable assays. Other options on the main menu include a data editing form and utility programs that produce various reports requested by the microbiologists and the project, and tools to generate the groupings for the final analyses. The analyses can be carried out in three ways: Each assay can be treated separately, the assays can be collectively treated for the whole zone as a group, or the assays can be collected in groups designated by the JPL Planetary Protection Manager. The latter approach was used to generate the final report because assays on the same equipment or similar equipment can be assumed to have been exposed to the same environment and cleaning. Thus, the statistics are improved by having a larger population, thereby reducing the standard deviation by the square root of N. For each method mentioned above, three reports are available. The first is a detailed report including all the data. This version was very useful in verifying the calculations. The second is a brief report that is similar to the full detailed report, but does not print out the data. The third is a grand total and summary report in which each assay requires only one line. For the first and second reports, most of the calculations are performed in the report section itself. For the third, all the calculations are performed directly in the query bound to the report. All the numerical results were verified by comparing them with Excel templates, then exporting the data from the Planetary Protection Analysis program to Excel.

Beaudet, Robert A.↗
NASA NTRS2013

Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2)

This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known.

Clark, Robert↗
DOE OSTI2021

Ergodic and nonergodic many-body dynamics in strongly nonlinear lattices

The study of nonlinear oscillator chains in classical many-body dynamics has a storied history going back to the seminal work of Fermi et al. [Los Alamos Scientific Laboratory Report No. LA-1940, 1955 (unpublished)]. Here, we introduce a family of such systems which consist of chains of N harmonically coupled particles with the nonlinearity introduced by confining the motion of each individual particle to a box or stadium with hard walls. The stadia are arranged on a one-dimensional lattice but they individually do not have to be one dimensional, thus permitting the introduction of chaos already at the lattice scale. For the most part we study the case where the motion is entirely one dimensional. We find that the system exhibits a mixed phase space for any finite value of N . Computations of Lyapunov spectra at randomly picked phase space locations and a direct comparison between Hamiltonian evolution and phase space averages indicate that the regular regions of phase space are not significant at large system sizes. While the continuum limit of our model is itself a singular limit of the integrable sinh Gordon theory, we do not see any evidence for the kind of nonergodicity famously seen in the work of Fermi et al. Finally, we examine the chain with particles confined to two-dimensional stadia where the individual stadium is already chaotic and find a much more chaotic phase space at small system sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Simulating plasma wave propagation on a superconducting quantum chip

Quantum computers may one day enable the efficient simulation of strongly coupled plasmas that lie beyond the reach of classical computation in regimes where quantum effects are important and the scale separation is large. Here, in this article, we take a first step toward efficient simulation of quantum plasmas by demonstrating linear plasma wave propagation on a superconducting quantum chip. Using high-fidelity and highly expressive device-native gates, combined with an error-mitigation technique, we simulate the scattering of laser pulses from inhomogeneous plasmas. Our approach is made feasible by the identification of a suitable local spin model whose excitations mimic plasma waves, and whose circuit implementation requires a lower gate count than other proposed approaches that would require a future fault-tolerant quantum computer. This work opens avenues to study more complicated phenomena that cannot be simulated efficiently on classical computers, such as nonlinear quantum dynamics when strongly coupled plasmas are driven out of equilibrium.

general physics↗
DOE OSTI2022

Hong–Ou–Mandel sensing via superradiant coupling of discrete fluorescent emitters

The Hong–Ou–Mandel (HOM) effect is a fascinating quantum phenomenon that defies classical explanation. Traditionally, remote nonlinear sources have been used to achieve coincident photons at the HOM beam splitter. Here, we suggest that the coincident emission source required for HOM interference can be created locally using superradiant near field coupled emitters positioned across the beam splitter gap. We show that sensitivity to permittivity changes in the beam splitter gap, and corresponding Fisher information can be substantially enhanced with HOM photon detection. Subsequently, we outline several strategies for integration of superradiant emitters with practical sensor systems. Taken together, these findings should pave a way for a wide array of near field HOM quantum sensors and novel quantum devices.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Extracting dynamical frequencies from invariants of motion in finite-dimensional nonlinear integrable systems

Integrable dynamical systems play an important role in many areas of science, including accelerator and plasma physics. An integrable dynamical system with n degrees of freedom possesses n nontrivial integrals of motion, and can be solved, in principle, by covering the phase space with one or more charts in which the dynamics can be described using action-angle coordinates. To obtain the frequencies of motion, both the transformation to action-angle coordinates and its inverse must be known in explicit form. However, no general algorithm exists for constructing this transformation explicitly from a set of n known (and generally coupled) integrals of motion. In this study we describe how one can determine the dynamical frequencies of the motion as functions of these n integrals in the absence of explicitly known action-angle variables, and we provide several examples.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗
DOE OSTI2024

Divide and conquer: Learning chaotic dynamical systems with multistep penalty neural ordinary differential equations

Forecasting high-dimensional dynamical systems is a fundamental challenge in various fields, such as geosciences and engineering. Neural Ordinary Differential Equations (NODEs), which combine the power of neural networks and numerical solvers, have emerged as a promising algorithm for forecasting complex nonlinear dynamical systems. However, classical techniques used for NODE training are ineffective for learning chaotic dynamical systems. In this work, we propose a novel NODE-training approach that allows for robust learning of chaotic dynamical systems. Here, our method addresses the challenges of non-convexity and exploding gradients associated with underlying chaotic dynamics. Training data trajectories from such systems are split into multiple, non-overlapping time windows. In addition to the deviation from the training data, the optimization loss term further penalizes the discontinuities of the predicted trajectory between the time windows. The window size is selected based on the fastest Lyapunov time scale of the system. Multi-step penalty(MP) method is first demonstrated on Lorenz equation, to illustrate how it improves the loss landscape and thereby accelerates the optimization convergence. MP method can optimize chaotic systems in a manner similar to least-squares shadowing with significantly lower computational costs. Our proposed algorithm, denoted the Multistep Penalty NODE, is applied to chaotic systems such as the Kuramoto-Sivashinsky equation, the two-dimensional Kolmogorov flow, and ERA5 reanalysis data for the atmosphere. It is observed that MP-NODE provide viable performance for such chaotic systems, not only for short-term trajectory predictions but also for invariant statistics that are hallmarks of the chaotic nature of these dynamics.

Chaotic dynamical systems↗
DOE OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Control of spinning flexible spacecraft by modal synthesis

A procedure is presented for the active control of a spinning flexible spacecraft. Such a system exhibits gyroscopic effects. The design of the controller is based on modal decomposition of the gyroscopic system. This modal decoupling procedure leads to a control mechanism implemented in modular form, which represents a distinct computational advantage over the control of the coupled system. Design procedures are demonstrated for two types of control algorithms, linear and nonlinear. The first represents classical linear feedback approach, and the second represents an application of on-off control, both types made feasible by the modal decomposition scheme.

Meirovitch, L.↗
NASA NTRS1976

Control of spinning flexible spacecraft by modal synthesis

A procedure is presented for the active control of a spinning flexible spacecraft. Such a system exhibits gyroscopic effects. The design of the controller is based on modal decomposition of the gyroscopic system. This modal decoupling procedure leads to a control mechanism implemented in modular form, which represents a distinct computational advantage over the control of the coupled system. Design procedures are demonstrated for two types of control algorithms, linear and nonlinear. The first represents classical linear feedback approach and the second represents an application of on-off control, both types made feasible by the modal decomposition scheme.

Meirovitch, L.↗
NASA NTRS1987

Extension of a semi-implicit shock-capturing algorithm for 3-D fully coupled, chemically reacting flows in generalized coordinates

A semiimplicit high-resolution shock-capturing method for multidimensional systems of hyperbolic conservation laws with stiff source terms has been developed by Yee and Shinn (1987). The goal of this work is to extend this method to solve the three-dimensional fully coupled Navier-Stokes equations for a hypersonic chemically reacting flow in generalized coordinates. In this formulation, the global continuity equation was replaced by all the species continuity equations. The shock-capturing technique is a second-order-accurate, symmetric total-variation-diminishing method which accounts fully and directly for the coupling among the fluid and all the species. To verify the current approach, it was implemented into an existing computer code which contained the MacCormack method. Test results for a five-species reacting flow are shown to be oscillation-free around the shock, and the time spent per iteration only doubles when compared to the result using classical way of supplying numerical dissipation. The extra computation is more than justified by the elimination of spurious oscillation and nonlinear instability associated with the classical shock-capturing schemes in computing hypersonic reacting flows.

Shinn, J. L.↗
DOE OSTI2025

Noise Constraints on Sensitivity Scaling in Quantum Nonlinear Metrology

Quantum-enhanced metrology surpasses classical metrology by improving estimation precision scaling with a resource 𝑁 (e.g., particle number or energy) from 1/√𝑁 to 1/𝑁. Through the use of nonlinear effects, Roy and Braunstein [Exponentially enhanced quantum metrology, Phys. Rev. Lett. 100, 220501 (2008)] derived a 1/𝑁 2 scaling. However, later works argued this exponential improvement is unphysical and that even modest gains, like 1/𝑁 2 , may vanish under noise. We show that, in the presence of small errors, the nonlinear interactions enabling metrological enhancement induce emergent errors. Here, the errors propagate through the sensing protocol and are magnified proportional to any intended nonlinear enhancement. We identify a critical value of the parameter to be estimated, for a fixed error, below which the emergent errors can be avoided.

Quantum Fisher information↗
DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗