Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear function”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Intelligent Control for Drag Reduction on the X-48B Vehicle

This paper focuses on the development of an intelligent control technology for in-flight drag reduction. The system is integrated with and demonstrated on the full X-48B nonlinear simulation. The intelligent control system utilizes a peak-seeking control method implemented with a time-varying Kalman filter. Performance functional coordinate and magnitude measurements, or independent and dependent parameters respectively, are used by the Kalman filter to provide the system with gradient estimates of the designed performance function which is used to drive the system toward a local minimum in a steepestdescent approach. To ensure ease of integration and algorithm performance, a single-input single-output approach was chosen. The framework, specific implementation considerations, simulation results, and flight feasibility issues related to this platform are discussed.

Griffin, Brian Joseph↗
NASA NTRS1988

Numerical solution methods for viscoelastic orthotropic materials

Numerical solution methods for viscoelastic orthotropic materials, specifically fiber reinforced composite materials, are examined. The methods include classical lamination theory using time increments, direction solution of the Volterra Integral, Zienkiewicz's linear Prony series method, and a new method called Nonlinear Differential Equation Method (NDEM) which uses a nonlinear Prony series. The criteria used for comparison of the various methods include the stability of the solution technique, time step size stability, computer solution time length, and computer memory storage. The Volterra Integral allowed the implementation of higher order solution techniques but had difficulties solving singular and weakly singular compliance function. The Zienkiewicz solution technique, which requires the viscoelastic response to be modeled by a Prony series, works well for linear viscoelastic isotropic materials and small time steps. The new method, NDEM, uses a modified Prony series which allows nonlinear stress effects to be included and can be used with orthotropic nonlinear viscoelastic materials. The NDEM technique is shown to be accurate and stable for both linear and nonlinear conditions with minimal computer time.

Gramoll, K. C.↗
DOE OSTI2026

Lyapunov-Based Iterative Learning of Regions of Attraction for Autonomous Systems

This paper proposes a novel algorithm for estimating the region of attraction of equilibrium points for nonlinear discrete-time autonomous systems. The method iteratively expands an initial estimate of the region of attraction by constructing unions of sublevel sets of learned functions parametrized as neural networks. Unlike conventional techniques that rely on a single global Lyapunov function, the proposed approach provides a collection of local Lyapunov-like functions, enabling richer representations and potentially larger region of attraction estimates. These functions are trained using sampled state-space data, and their Lipschitz continuity ensures that desirable properties extend beyond the training samples. The devised strategy is tested via numerical simulations, demonstrating the effectiveness of the proposed approach.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Unbalance Response Prediction for Rotors on Ball Bearings Using Speed and Load Dependent Nonlinear Bearing Stiffness

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus an accurate rotordynamic analysis requires that bearing forces corresponding to the actual bearing deflection be utilized. For this work bearing forces were calculated by COBRA-AHS, a recently developed rolling-element bearing analysis code. Bearing stiffness was found to be a strong function of bearing deflection, with higher deflection producing markedly higher stiffness. Curves fitted to the bearing data for a range of speeds and loads were supplied to a flexible rotor unbalance response analysis. The rotordynamic analysis showed that vibration response varied nonlinearly with the amount of rotor imbalance. Moreover, the increase in stiffness as critical speeds were approached caused a large increase in rotor and bearing vibration amplitude over part of the speed range compared to the case of constant bearing stiffness. Regions of bistable operation were possible, in which the amplitude at a given speed was much larger during rotor acceleration than during deceleration. A moderate amount of damping will eliminate the bistable region, but this damping is not inherent in ball bearings.

Fleming, David P.↗
DOE OSTI2020

Nonlinear time-domain spectroscopy near a band inversion

We develop a theory for the nonlinear time-domain response of a semimetal driven by an ultrafast optical pulse. At quadratic order in the driving field we find that the response near a band inversion transition contains a coherent oscillating component proportional to field intensity with a frequency that can be tuned over a wide range (from terahertz to the near IR) selected by the chemical potential. Finally, we illustrate the effect by calculating the induced current as a function of time for a two-band system with time-reversal symmetry but broken parity symmetry to study semimetals where large Berry curvature near a band inversion transition promotes this nonlinear response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1971

Dynamic responses of railroad car models to vertical and lateral rail inputs

Simplified dynamic models were applied in a study of vibration in a high-speed railroad car. The mathematical models used were a four-degree-of-freedom model for vertical responses to vertical rail inputs and a ten-degree-of-freedom model for lateral response to lateral or rolling (cross-level) inputs from the rails. Elastic properties of the passenger car body were represented by bending and torsion of a uniform beam. Rail-to-car (truck) suspensions were modeled as spring-mass-dashpot oscillators. Lateral spring nonlinearities approximating certain complicated truck mechanisms were introduced. The models were excited by displacement and, in some cases, velocity inputs from the rails by both deterministic (including sinusoidal) and random input functions. Results were obtained both in the frequency and time domains. Solutions in the time domain for the lateral model were obtained for a wide variety of transient and random inputs generated on-line by an analog computer. Variations in one of the damping properties of the lateral car suspension gave large fluctuations in response over a range of car speeds for a given input. This damping coefficient was significant in reducing lateral car responses that were higher for nonlinear springs for three different inputs.

Sewall, J. L.↗
DOE OSTI2025

Pressure-induced structural and dielectric changes in liquid water at room temperature

Understanding the pressure-dependent dielectric properties of water is crucial for a wide range of scientific and practical applications. In this study, we employ a deep neural network trained on density functional theory data to investigate the dielectric properties of liquid water at room temperature across a pressure range of 0.1–1000 MPa. We observe a nonlinear increase in the static dielectric constant ɛ 0 with increasing pressure, a trend that is qualitatively consistent with experimental observations. This increase in ɛ 0 is primarily attributed to the increase in water density under compression, which enhances collective dipole fluctuations within the hydrogen-bonding network as well as the dielectric response. Despite the increase in ɛ 0 , our results reveal a decrease in the Kirkwood correlation factor G K with increasing pressure. Furthermore, this decrease in G K is attributed to pressure-induced structural distortions in the hydrogen-bonding network, which weaken dipolar correlations by disrupting the ideal tetrahedral arrangement of water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Nonlinear Elastic Effects on the Energy Flux Deviation of Ultrasonic Waves in GR/EP Composites

In isotropic materials, the direction of the energy flux (energy per unit time per unit area) of an ultrasonic plane wave is always along the same direction as the normal to the wave front. In anisotropic materials, however, this is true only along symmetry directions. Along other directions, the energy flux of the wave deviates from the intended direction of propagation. This phenomenon is known as energy flux deviation and is illustrated. The direction of the energy flux is dependent on the elastic coefficients of the material. This effect has been demonstrated in many anisotropic crystalline materials. In transparent quartz crystals, Schlieren photographs have been obtained which allow visualization of the ultrasonic waves and the energy flux deviation. The energy flux deviation in graphite/epoxy (gr/ep) composite materials can be quite large because of their high anisotropy. The flux deviation angle has been calculated for unidirectional gr/ep composites as a function of both fiber orientation and fiber volume content. Experimental measurements have also been made in unidirectional composites. It has been further demonstrated that changes in composite materials which alter the elastic properties such as moisture absorption by the matrix or fiber degradation, can be detected nondestructively by measurements of the energy flux shift. In this research, the effects of nonlinear elasticity on energy flux deviation in unidirectional gr/ep composites were studied. Because of elastic nonlinearity, the angle of the energy flux deviation was shown to be a function of applied stress. This shift in flux deviation was modeled using acoustoelastic theory and the previously measured second and third order elastic stiffness coefficients for T300/5208 gr/ep. Two conditions of applied uniaxial stress were considered. In the first case, the direction of applied uniaxial stress was along the fiber axis (x3) while in the second case it was perpendicular to the fiber axis along the laminate stacking direction (x1).

Prosser, William H.↗
NASA NTRS1976

Gradient optimization and nonlinear control

The book represents an introduction to computation in control by an iterative, gradient, numerical method, where linearity is not assumed. The general language and approach used are those of elementary functional analysis. The particular gradient method that is emphasized and used is conjugate gradient descent, a well known method exhibiting quadratic convergence while requiring very little more computation than simple steepest descent. Constraints are not dealt with directly, but rather the approach is to introduce them as penalty terms in the criterion. General conjugate gradient descent methods are developed and applied to problems in control.

Hasdorff, L.↗
NASA NTRS1977

Nonlinear viscoelasticity and relaxation phenomena of polymer solids

In the light of a three-chain model of statistical network theories of rubberlike elastic models, it is assumed that the free energy function of incompressible viscoelastic polymer solids is a separable, symmetric function of the principle stretch ratios and the hidden thermodynamic coordinates along the same directions. This assumption leads to a characterization of those viscoelastic polymer solids which exhibit the property of factorizability between the time and strain functions.

Peng, S. T. J.↗
DOE OSTI2023

Optimization of transformer ratio and beam loading in a plasma wakefield accelerator with a structure-exploiting algorithm

Plasma-based acceleration has emerged as a promising candidate as an accelerator technology for a future linear collider or a next-generation light source. We consider the plasma wakefield accelerator (PWFA) concept where a plasma wave wake is excited by a particle beam and a trailing beam surfs on the wake. For a linear collider, the energy transfer efficiency from the drive beam to the wake and from the wake to the trailing beam must be large, while the emittance and energy spread of the trailing bunch must be preserved. One way to simultaneously achieve this when accelerating electrons is to use longitudinally shaped bunches and nonlinear wakes. In the linear regime, there is an analytical formalism to obtain the optimal shapes. In the nonlinear regime, however, the optimal shape of the driver to maximize the energy transfer efficiency cannot be precisely obtained because currently no theory describes the wake structure and excitation process for all degrees of nonlinearity. In addition, the ion channel radius is not well defined at the front of the wake where the plasma electrons are not fully blown out by the drive beam. We present results using a novel optimization method to effectively determine a current profile for the drive and trailing beam in PWFA that provides low energy spread, low emittance, and high acceleration efficiency. We parameterize the longitudinal beam current profile as a piecewise-linear function and define optimization objectives. For the trailing beam, the algorithm converges quickly to a nearly inverse trapezoidal trailing beam current profile similar to that predicted by the ultrarelativistic limit of the nonlinear wakefield theory. For the drive beam, the beam profile found by the optimization in the nonlinear regime that maximizes the transformer ratio also resembles that predicted by linear theory. Furthermore, the current profiles found from the optimization method provide higher transformer ratios compared with the linear ramp predicted by the relativistic limit of the nonlinear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

Low-Velocity Impact Response of Sandwich Beams with Functionally Graded Core

The problem of low-speed impact of a one-dimensional sandwich panel by a rigid cylindrical projectile is considered. The core of the sandwich panel is functionally graded such that the density, and hence its stiffness, vary through the thickness. The problem is a combination of static contact problem and dynamic response of the sandwich panel obtained via a simple nonlinear spring-mass model (quasi-static approximation). The variation of core Young's modulus is represented by a polynomial in the thickness coordinate, but the Poisson's ratio is kept constant. The two-dimensional elasticity equations for the plane sandwich structure are solved using a combination of Fourier series and Galerkin method. The contact problem is solved using the assumed contact stress distribution method. For the impact problem we used a simple dynamic model based on quasi-static behavior of the panel - the sandwich beam was modeled as a combination of two springs, a linear spring to account for the global deflection and a nonlinear spring to represent the local indentation effects. Results indicate that the contact stiffness of thc beam with graded core Increases causing the contact stresses and other stress components in the vicinity of contact to increase. However, the values of maximum strains corresponding to the maximum impact load arc reduced considerably due to grading of thc core properties. For a better comparison, the thickness of the functionally graded cores was chosen such that the flexural stiffness was equal to that of a beam with homogeneous core. The results indicate that functionally graded cores can be used effectively to mitigate or completely prevent impact damage in sandwich composites.

Apetre, N. A.↗
NASA NTRS2001

Damage Tolerance of Sandwich Plates With Debonded Face Sheets

A nonlinear finite element analysis was performed to simulate axial compression of sandwich beams with debonded face sheets. The load - end-shortening diagrams were generated for a variety of specimens used in a previous experimental study. The energy release rate at the crack tip was computed using the J-integral, and plotted as a function of the load. A detailed stress analysis was performed and the critical stresses in the face sheet and the core were computed. The core was also modeled as an isotropic elastic-perfectly plastic material and a nonlinear post buckling analysis was performed. A Graeco-Latin factorial plan was used to study the effects of debond length, face sheet and core thicknesses, and core density on the load carrying capacity of the sandwich composite. It has been found that a linear buckling analysis is inadequate in determining the maximum load a debonded sandwich beam can carry. A nonlinear post-buckling analysis combined with an elastoplastic model of the core is required to predict the compression behavior of debonded sandwich beams.

Sankar, Bhavani V.↗
NASA NTRS2016

On the Connection Between Microbursts and Nonlinear Electronic Structures in Planetary Radiation Belts

Using a dynamical-system approach, we have investigated the efficiency of large-amplitude whistler waves for causing microburst precipitation in planetary radiation belts by modeling the microburst energy and particle fluxes produced as a result of nonlinear wave-particle interactions. We show that wave parameters, consistent with large amplitude oblique whistlers, can commonly generate microbursts of electrons with hundreds of keV-energies as a result of Landau trapping. Relativistic microbursts (greater than 1 MeV) can also be generated by a similar mechanism, but require waves with large propagation angles Theta (sub k)B greater than 50 degrees and phase-speeds v(sub phi) greater than or equal to c/9. Using our result for precipitating density and energy fluxes, we argue that holes in the distribution function of electrons near the magnetic mirror point can result in the generation of double layers and electron solitary holes consistent in scales (of the order of Debye lengths) to nonlinear structures observed in the radiation belts by the Van Allen Probes. Our results indicate a relationship between nonlinear electrostatic and electromagnetic structures in the dynamics of planetary radiation belts and their role in the cyclical production of energetic electrons (E greater than or equal to 100 keV) on kinetic timescales, which is much faster than previously inferred.

plasmas↗
DOE OSTI2025

Nonlinear analog processing with anisotropic nonlinear films

Digital signal processing is the cornerstone of several modern-day technologies, yet in multiple applications it faces critical bottlenecks related to memory and speed constraints. Thanks to recent advances in metasurface design and fabrication, light-based analog computing has emerged as a viable option to partially replace or augment digital approaches. Several light-based analog computing functionalities have been demonstrated using patterned flat optical elements, with great opportunities for integration in compact nanophotonic systems. So far, however, the available operations have been restricted to the linear regime, limiting the impact of this technology to a compactification of Fourier optics systems. In this paper, we introduce nonlinear operations to the field of metasurface-based analog optical processing, demonstrating that nonlinear optical phenomena, combined with nonlocality in flat optics, can be leveraged to synthesize kernels beyond linear Fourier optics, paving the way to a broad range of new opportunities. As a practical demonstration, we report the experimental synthesis of a class of nonlinear operations that can be used to realize broadband, polarization-selective analog-domain edge detection.

analog image processing↗
DOE OSTI2022

Measuring nonlinearities of a cantilever beam using a low-cost efficient wireless intelligent sensor for strain (LEWIS-S)

In the context of experimental vibration data, strain gauges can obtain linear and nonlinear dynamic measurements. However, measuring strain can be disincentivizing and expensive due to the complexity of data acquisition systems, lack of portability, and high costs. This research introduces the use of a low-cost efficient wireless intelligent sensor for strain (LEWIS-S) that is based on a portable-sensor-design platform that streamlines strain sensing. Additionally, the softening behavior of a cantilever beam with geometric and inertial nonlinearities is characterized by the LEWIS-S based on high force level inputs. Two experiments were performed on a nonlinear cantilever beam with measurements obtained by the LEWIS-S sensor and an accelerometer. First, a sine sweep test was performed through the fundamental resonance of the system, then a ring-down test was performed from a large initial static deformation. Good agreement was revealed in quantities of interest such as frequency response functions, the continuous wavelet transforms, and softening behavior in the backbone curves.

42 ENGINEERING↗
DOE OSTI2021

Pyomo.GDP: an ecosystem for logic based modeling and optimization development

We present three core principles for engineering-oriented integrated modeling and optimization tool sets—intuitive modeling contexts, systematic computer-aided reformulations, and flexible solution strategies—and describe how new developments in Pyomo.GDP for Generalized Disjunctive Programming (GDP) advance this vision. We describe a new logical expression system implementation for Pyomo.GDP allowing for a more intuitive description of logical propositions. The logical expression system supports automated reformulation of these logical constraints to linear constraints. We also describe two new logic-based global optimization solver implementations built on Pyomo.GDP that exploit logical structure to avoid “zero-flow” numerical difficulties that arise in nonlinear network design problems when nodes or streams disappear. These new solvers also demonstrate the capability to link to external libraries for expanded functionality within an integrated implementation. We present these new solvers in the context of a flexible array of solution paths available to GDP models. Finally, we present results on a new library of GDP models demonstrating the value of multiple solution approaches.

42 ENGINEERING↗
DOE OSTI2021

Minimal two-body quantum absorption refrigerator

Here, we study the phenomenon of absorption refrigeration, where refrigeration is achieved by heating instead of work, in two different setups: a minimal set up based on coupled qubits and two nonlinearly coupled resonators. Considering ZZ interaction between the two qubits, we outline the basic ingredients required to achieve cooling. Using local as well as global master equations, we observe that inclusion of an XX type term in the qubit-qubit coupling is detrimental to cooling. We compare the cooling effect obtained in the qubit case with that of nonlinearly coupled resonators (multilevel system) where the ZZ interaction translates to a Kerr-type nonlinearity. For small to intermediate strengths of nonlinearity, we observe that multilevel quantum systems, for example qutrits, give better cooling effect compared to the qubits. Using Keldysh nonequilibrium Green's function formalism, we go beyond first order sequential tunneling processes and study the effect of higher order processes on refrigeration. We find a reduced cooling effect compared to the master equation calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗