Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Transient dynamics, damping, and mode coupling of nonlinear systems with internal resonances

The study of both linear and nonlinear structural vibrations routinely circles the concise yet complex problem of choosing a set of coordinates which yield simple equations of motion. In both experimental and mathematical methods, that choice is a difficult one because of measurement, computational, and interpretation difficulties. Often times, researchers choose to solve their problems in terms of linear, undamped mode shapes because they are easy to obtain; however, this is known to give rise to complicated phenomena such as mode coupling and internal resonance. This work considers the nature of mode coupling and internal resonance in systems containing non-proportional damping, linear detuning, and cubic nonlinearities through the method of multiple scales as well as instantaneous measures of effective damping. The energy decay observed in the structural modes is well approximated by the slow-flow equations in terms of the modal amplitudes, and it is shown how mode coupling enhances the damping observed in the system. Moreover, in the presence of a 3:1 internal resonance between two modes, the nonlinearities not only enhance the dissipation, but can allow for the exchange and transfer of energy between the resonant modes. However, this exchange depends on the resonant phase between the modes and is proportional to the energy in the lowest mode. The results of the analysis tie together interpretations used by both experimentalists and theoreticians to study such systems and provide a more concrete way to interpret these phenomena.

42 ENGINEERING↗
DOE OSTI2023

Impacts of Permeability Heterogeneity and Background Flow on Supercritical CO 2 Dissolution in the Deep Subsurface

Motivated by CO 2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO 2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi-Gaussian log-permeability fields: mean, variance, and correlation length. Hundreds of high-resolution simulations were performed using the PFLOTRAN finite volume software to model CO 2 dissolution in a kilometer-scale aquifer over 1,000 years. Predictive dimensionless scaling relationships relating CO 2 dissolution rate to heterogeneity statistics, Rayleigh (Ra) and Péclet (Pe) numbers were developed for both gravitationally dominated free convection to background flow-dominated forced convection regimes. An empirical criterion, Pe = Ra 3/4 , was discovered for regime transition. All simulations converged quickly to a quasi-steady, approximately linear dissolution rate. However, this rate displayed profound variability between permeability field realizations sharing the same heterogeneity statistics, even under mild permeability heterogeneity. In general, increased heterogeneity was associated with a lower mean and higher variance of dissolution rate, undesirable from a CCS design perspective. The relationship between dissolution rate and background flow was found to be complex and nonlinear. Dimensionless scaling relationships were uncovered for a number of special cases. Results call into question the validity of the Boussinesq approximation in the context of modest-to-high background flow rates and the general applicability of numerical simulations without background flow.

42 ENGINEERING↗
NASA NTRS1993

H(sub infinity)-type filter for spacecraft attitude estimation

A nonlinear filtering theory from a deterministic point of view is presented and an application to attitude determination is considered. The approach that is taken in this paper is motivated largely by the H(sun infinity) control and estimation theory for linear systems which has been evolved within the last decade. Rather than formulating the estimation problem as a game played by two adversaries, as has been done in the linear case, we employ in this work some notions from the theory of dissipative systems as a vehicle for arriving at certain Hamilton-Jacobi inequality, which in turn, provides a solution to the filtering problem, whenever it is satisfied. Application of this method to a linear estimation problem and to the problem of estimating a spacecraft attitude quaternion and gyro drift bias vector are presented. In limiting cases, these give the Kalman filter and the extended Kalman filter, respectively. The main advantages of this approach over its probabilistic counterpart are that this approach does not require a prior knowledge of any statistics, and that in general it is more amenable to a quantitative assessment regarding approximations such as linearization, and that in certain cases this approach yields an exact solution to the nonlinear filtering.

Markley, F. Landis↗
DOE OSTI2018

A two-grid penalty method for American options

In this paper we consider the pricing of American options, governed by a partial differential complementarity problem. The differential problem is first approximated by a semi-linear PDE using two distinct penalty approaches which are well known in computational finance. We then initiate the two-grid algorithm by solving the nonlinear problem on a coarse grid and further the linearized in the interpolated coarse-grid solution problem on a fine grid. By means of the maximum principle the algorithm is shown to be of fourth order convergence rate in space. Numerical experiments verify the presented two-grid approach where we draw some interesting conclusions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Numerical and experimental investigation of the bending response of thin-walled composite cylinders

A numerical and experimental investigation of the bending behavior of six eight-ply graphite-epoxy circular cylinders is presented. Bending is induced by applying a known end-rotation to each end of the cylinders, analogous to a beam in bending. The cylinders have a nominal radius of 6 inches, a length-to-radius ratio of 2 and 5, and a radius-to-thickness ratio of approximately 160. A (+/- 45/0/90)S quasi-isotropic layup and two orthotropic layups, (+/- 45/0 sub 2)S and (+/- 45/90 sub 2)S, are studied. A geometrically nonlinear special-purpose analysis, based on Donnell's nonlinear shell equations, is developed to study the prebuckling responses and gain insight into the effects of non-ideal boundary conditions and initial geometric imperfections. A geometrically nonlinear finite element analysis is utilized to compare with the prebuckling solutions of the special-purpose analysis and to study the buckling and post buckling responses of both geometrically perfect and imperfect cylinders. The imperfect cylinder geometries are represented by an analytical approximation of the measured shape imperfections. Extensive experimental data are obtained from quasi-static tests of the cylinders using a test fixture specifically designed for the present investigation. A description of the test fixture is included. The experimental data are compared to predictions for both perfect and imperfect cylinder geometries. Prebuckling results are presented in the form of displacement and strain profiles. Buckling end-rotations, moments, and strains are reported, and predicted mode shapes are presented. Observed and predicted moment vs. end-rotation relations, deflection patterns, and strain profiles are illustrated for the post buckling responses. It is found that a geometrically nonlinear boundary layer behavior characterizes the prebuckling responses. The boundary layer behavior is sensitive to laminate orthotropy, cylinder geometry, initial geometric imperfections, applied end-rotation, and non-ideal boundary conditions. Buckling end-rotations, strains, and moments are influenced by laminate orthotropy and initial geometric imperfections. Measured buckling results correlate well with predictions for the geometrically imperfect specimens. The postbuckling analyses predict equilibrium paths with a number of scallop-shaped branches that correspond to unique deflection patterns. The observed postbuckling deflection patterns and measured strain profiles show striking similarities to the predictions in some cases. Ultimate failure of the cylinders is attributed to an interlaminar shear failure mode along the nodal lines of the postbuckling deflection patterns.

Fuchs, J. P.↗
DOE OSTI2024

Second harmonic generation for estimating state of charge of lithium-ion batteries

This study applied the nonlinear ultrasonic method, second harmonic generation, to precisely estimate the state of charge (SoC) in lithium-ion batteries. The second harmonic of the longitudinal wave is generated on a pouch cell battery at 5 MHz with a through-transmission setup. The relative nonlinear parameter β' is determined by analyzing the amplitudes at the fundamental and second harmonic frequencies. To enhance the nonlinear parameter's measurement accuracy, multiple excitation amplitudes are employed. Two separate charge/discharge tests (four-cycle and eight-cycle) are conducted on the battery at a rate of C/10. The nonlinear parameter is measured periodically during the charge/discharge process, and temperature compensation is applied to the measurement. The correlation curves between the nonlinear parameter and the actual SoC align well for the four-cycle and eight-cycle tests, and a robust linear relationship is observed for both correlation curves. A linear model and a second-order polynomial model are applied to fit the correlation using all data points from both tests. The two models are employed to validate the SoC prediction on a second battery by using a four-cycle test. Furthermore, the results indicate that both models can predict the SoC with an accuracy of approximately 3%, whereas the polynomial model demonstrates smaller errors in the regions near 0% and 100% SoC. Therefore, the nonlinear parameter β', measured through the second harmonic generation, can effectively predict lithium-ion battery SoC with an accuracy of less than 3%.

25 ENERGY STORAGE↗
NASA NTRS2021

Expected Collision Rates for Tracked Satellites

This analysis estimates collision risks between tracked satellites based on the statistically expected collision rate _Nc and number of collisions Nc, which are closely related to the collision probability Pc. For isolated encounters, Nc equals Pc. For multi-encounter interactions, Nc can exceed Pc and is significantly easier to approximate semi-analytically. Estimating Nc requires integrating the collision rate over time throughout an interaction, with each required _Nc value calculated as a statistical expectation value based on the uncertainty distributions of initial orbital states. The uncertainty distributions are estimated from orbit determination analyses of satellite trajectory measurements. The formulation accounts for the nonlinear orbital motion of the satellites, and derives Nc and _Nc expressions for Monte Carlo simulations as well as semi-analytical approximations for single- and multi-encounter interactions.

Collision Rates↗
NASA NTRS1991

Investigation of the structure of the electromagnetic field and related phenomena, generated by the Active Satellite

The altitude dependencies of the moduli of the electric field E in the VLF and LF frequency bands (f sub B much less than F less than f sub B) and in the altitude range of the ionosphere Z equals (400 to 2500) km up to Z equals 6000 km (the bottom of the magnetosphere) were calculated by the linear theory. The amplitudes of the field have large maxima in four regions: the axis field (E sub o) close to the direction of the Earth's magnetic field line B sub o, beta approximately 0 degrees, the fields (E sub St), (E sub RevSt) and (E sub Res) in the Storey, Reversed Story and Resonance cones, beta approximately (0 approaches 20) degrees. Their maxima are very pronounced close to the low hybrid frequency F sub L. The nonlinear heating of the magnetoplasma under the action of an electric field Ee (sup iwt) was recently expanded by the macroscopic theory by the author. The velocities, collision frequencies and temperatures of all the constituents of a magnetoplasma - electrons, ions, and neutral particles - are taken into account. Formulae and numerical results are presented for the ionosphere in the frequency band F equals (1 to 10 exp 4) kHz and altitude range Z approximately (100 - 1000) km. Some results of calculations by the self consistent solution of the basis system of equations are also discussed.

Alpert, Yakov L.↗
DOE OSTI2017

On non-local energy transfer via zonal flow in the Dimits shift

The two-dimensional Terry–Horton equation is shown to exhibit the Dimits shift when suitably modified to capture both the nonlinear enhancement of zonal/drift-wave interactions and the existence of residual Rosenbluth–Hinton states. This phenomenon persists through numerous simplifications of the equation, including a quasilinear approximation as well as a four-mode truncation. It is shown that the use of an appropriate adiabatic electron response, for which the electrons are not affected by the flux-averaged potential, results in an $\boldsymbol{E}\times \boldsymbol{B}$ nonlinearity that can efficiently transfer energy non-locally to length scales of the order of the sound radius. The size of the shift for the nonlinear system is heuristically calculated and found to be in excellent agreement with numerical solutions. The existence of the Dimits shift for this system is then understood as an ability of the unstable primary modes to efficiently couple to stable modes at smaller scales, and the shift ends when these stable modes eventually destabilize as the density gradient is increased. This non-local mechanism of energy transfer is argued to be generically important even for more physically complete systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Nonlinear thermodynamic computing out of equilibrium

We present the design for a thermodynamic computer that can perform arbitrary nonlinear calculations in or out of equilibrium. Simple thermodynamic circuits, fluctuating degrees of freedom in contact with a thermal bath and confined by a quartic potential, display an activity that is a nonlinear function of their input. Such circuits can therefore be regarded as thermodynamic neurons, and can serve as the building blocks of networked structures that act as thermodynamic neural networks, universal function approximators whose operation is powered by thermal fluctuations. We simulate a digital model of a thermodynamic neural network, and show that its parameters can be adjusted by genetic algorithm to perform nonlinear calculations at specified observation times, regardless of whether the system has attained thermal equilibrium. This work expands the field of thermodynamic computing beyond the regime of thermal equilibrium, enabling fully nonlinear computations, analogous to those performed by classical neural networks, at specified observation times.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
DOE OSTI2023

An Ensemble Score Filter for Tracking High-Dimensional Nonlinear Dynamical Systems

We propose an ensemble score filter (EnSF) for solving high-dimensional nonlinear filtering problems with superior accuracy. A major drawback of existing filtering methods, e.g., particle filters or ensemble Kalman filters, is the low accuracy in handling high-dimensional and highly nonlinear problems. EnSF attacks this challenge by exploiting the score-based diffusion model, defined in a pseudo-temporal domain, to characterizing the evolution of the filtering density. EnSF stores the information of the recursively updated filtering density function in the score function, in stead of storing the information in a set of finite Monte Carlo samples (used in particle filters and ensemble Kalman filters). Unlike existing diffusion models that train neural networks to approximate the score function, we develop a training-free score estimation that uses mini-batch-based Monte Carlo estimator to directly approximate the score function at any pseudo-spatial-temporal location, which provides sufficient accuracy in solving high-dimensional nonlinear problems as well as saves tremendous amount of time spent on training neural networks. Another essential aspect of EnSF is its analytical update step, gradually incorporating data information into the score function, which is crucial in mitigating the degeneracy issue faced when dealing with very high-dimensional nonlinear filtering problems. High-dimensional Lorenz systems are used to demonstrate the performance of our method. EnSF provides surprisingly impressive performance in reliably tracking extremely high-dimensional Lorenz systems (up to 1,000,000 dimension) with highly nonlinear observation processes, which is a well-known challenging problem for existing filtering methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Numerical Simulation of the Instability of an Inviscid Liquid Jet in a Coflowing Compressible Airstream

The nonlinear interfacial instability of a liquid jet in a coflowing compressible airstream is studied numerically. A high-resolution scheme which has second-order accuracy in space and time is coupled with a Lagrangian marker particle algorithm to visualize the large-scale motion of the interfaces in compressible flow. A numerical algorithm based on an approximate equation of state of a compressible liquid is developed to allow this two-fluid system to be governed by the nonlinear unsteady Euler equations in conservative form. The initial growth of small disturbances given by the simulations agrees well with linear theory. The process of jet disruption in compressible flow is demonstrated to consist of the formation of liquid spikes, interweaving of the gas and liquid and stretching and detachment of the liquid main center core.

Li, H. -S.↗
DOE OSTI2017

A simple method of measuring the effective SRS coefficient in single-mode optical fibres and its applicability limits

This paper describes techniques for measuring the SRS coefficient in a wide spectral range, including the region of small Stokes shifts. A simple, approximate method is proposed for evaluating the SRS coefficient near a gain peak. Spectral dependences of the SRS coefficient are presented for various telecom fibres. (nonlinear optical phenomena)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRSDate not supplied

Improvements in Iterative Convergence of FUN3D Solutions

This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency of FUN3Dsolutions. In contrast to the baseline iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear CFL ramping, HANIM is based upon a hierarchy of modules including pre conditioner, generalized conjugate residual, realizability check, nonlinear control,and CFL adaption modules. HANIM performance is systematically compared with the performance of the baseline solver. The iterative solutions are compared for three aerodynamic benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-body juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmar as model and its nonlinear extension based on quadratic constitutive relations.

Li Wang↗
DOE OSTI2026

Analytic Chromaticity Formulas for the Muon g -2 experiment at fermilab

The Muon g-2 Experiment (E989) at Fermilab measured the muon anomalous magnetic moment aμ with unprecedented precision of 127 ppb using a storage ring. The experiment pushed the limits in terms of accuracy and systematic errors. Achieving the 78-ppb total systematic uncertainty required precise beam dynamics modeling and corrections for effects on the measured muon precession frequency. Here, we derived the first analytic aberration formulas up to the second order for the muon g-2 storage ring’s combined-function electrostatic quadrupoles (superimposed magnetic dipole and electric quadrupole fields) using an order-by-order perturbation method. From these, we obtained the exact chromaticity formulas for three ring models of different granularity and validated them against numerical calculations using COSY INFINITY, achieving analytic-numerical agreement to ⁠$\mathscr{O}(10^{-10})$. This work resolved discrepancies between previous approximate derivations and provided essential beam dynamics results for Runs 4-6 analyses. We also calculated nonlinear chromaticities up to ninth order. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous muon g-2 experiment at Brookhaven National Laboratory, with the final result announced in June 2025.

aberrations↗
DOE OSTI2022

Modeling pre-Exascale AMR Parallel I/O Workloads via Proxy Applications

The present work investigates the modeling of preexascale input/output (I/O) workloads of Adaptive Mesh Refinement (AMR) simulations through a simple proxy application. We collect data from the AMReX Castro framework running on the Summit supercomputer for a wide range of scales and mesh partitions for the hydrodynamic Sedov case as a baseline to provide sufficient coverage to the formulated proxy model. The non-linear analysis data production rates are quantified as a function of a set of input parameters such as output frequency, grid size, number of levels, and the Courant-Friedrichs-Lewy (CFL) condition number for each rank, mesh level and simulation time step. Linear regression is then applied to formulate a simple analytical model which allows to translate AMReX inputs into MACSio proxy I/O application parameters, resulting in a simple “kernel” approximation for data production at each time step. Results show that MACSio can simulate actual AMReX nonlinear “static” I/O workloads to a certain degree of confidence on the Summit supercomputer using the present methodology. The goal is to provide an initial level of understanding of AMR I/O workloads via lightweight proxy applications models to facilitate autotune data management strategies in anticipation of exascale systems.

Godoy, William↗
DOE OSTI2024

A Data-Driven Approach for High-Impedance Fault Localization in Distribution Systems [SWR-24-13]

This software provides a data-driven approach for efficiently identifying high impedance faults (HIFs) in distribution systems. To tackle the nonlinearity of the voltage current trajectory of HIFs, we first formulate linear least squares problems to approximate the trajectory with piecewise functions. Then we collect the function features of all segments as inputs and use the support vector machine approach to efficiently identify HIFs at different locations.

Zhou, Yuqi↗