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 775 records · Page 43

Model of a rotating magnetic cloud

The possibility that magnetic clouds rotate while they propagate antisunward was investigated. Magnetic clouds are modeled as magnetic flux ropes which rotate rigidly about the axis of symmetry. An ideal magnetohydrodynamic model, in which the evolution of the magnetic structure is related to the time evolution of the angular frequency, is developed. A class of 'separable' magnetic fields is employed to reduce the problem to a nonlinear ordinary differential equation for the evolution function, and it is solved numerically. The corresponding effective potential gives rise to two modes of evolution--expansion and oscillation--depending on the energy and on the value of a dimensionless parameter, k. Parameter k depends on the gas pressure, the ratio of the magnetic field components, and the frequency of rotation. There is a critical value of k, k(sub c), above which the oscillatory regime disappears and the flux rope invariably expands, regardless of the energy. Below k(sub c) the energy determines whether the configuration is confined or unbounded. Rotation always helps expansion by lowering the potential barrier. A data example was studied and features which are interpreted as signatures of rotation are presented. The angular speed is comparable to the Alfven speed, and the core of the rotating cloud completes on average one full revolution every three days at 1 AU. The parameter k is calculated from observations, and it is found to be close to, but below, critical. Only three out of the nine clouds examined showed signatures of rotation. Theoretical analysis suggests that close to the Sun rotation effects may play a more important role in the evolution of magnetic clouds than 1 AU.

Farrugia, C. J.↗

Broadband Achromatic Phase Shifter for a Nulling Interferometer

Nulling interferometry is a technique for imaging exoplanets in which light from the parent star is suppressed using destructive interference. Light from the star is divided into two beams and a phase shift of radians is introduced into one of the beams. When the beams are recombined, they destructively interfere to produce a deep null. For monochromatic light, this is implemented by introducing an optical path difference (OPD) between the two beams equal to lambda/2, where lambda is the wavelength of the light. For broadband light, however, a different phase shift will be introduced at each wavelength and the two beams will not effectively null when recombined. Various techniques have been devised to introduce an achromatic phase shift a phase shift that is uniform across a particular bandwidth. One popular technique is to use a series of dispersive elements to introduce a wavelength-dependent optical path in one or both of the arms of the interferometer. By intelligently choosing the number, material and thickness of a series of glass plates, a nearly uniform, arbitrary phase shift can be introduced between two arms of an interferometer. There are several constraints that make choosing the number, type, and thickness of materials a difficult problem, such as the size of the bandwidth to be nulled. Several solutions have been found for bandwidths on the order of 20 to 30 percent (Delta(lambda)/lambda(sub c)) in the mid-infrared region. However, uniform phase shifts over a larger bandwidth in the visible regime between 480 to 960 nm (67 percent) remain difficult to obtain at the tolerances necessary for exoplanet detection. A configuration of 10 dispersive glass plates was developed to be used as an achromatic phase shifter in nulling interferometry. Five glass plates were placed in each arm of the interferometer and an additional vacuum distance was also included in the second arm of the interferometer. This configuration creates a phase shift of pi radians with an average error of 5.97 x 10(exp -8) radians and standard deviation of 3.07 x 10(exp -4) radians. To reduce ghost reflections and interference effects from neighboring elements, the glass plates are tilted such that the beam does not strike each plate at normal incidence. Reflections will therefore walk out of the system and not contribute to the intensity when the beams are recombined. Tilting the glass plates, however, introduces several other problems that must be mitigated: (1) the polarization of a beam changes when refracted at an interface at non-normal incidence; (2) the beam experiences lateral chromatic spread as it traverses multiple glass plates; (3) at each surface, wavelength- dependent intensity losses will occur due to reflection. For a fixed angle of incidence, each of these effects must be balanced between each arm of the interferometer in order to ensure a deep null. The solution was found using a nonlinear optimization routine that minimized an objective function relating phase shift, intensity difference, chromatic beam spread, and polarization difference to the desired parameters: glass plate material and thickness. In addition to providing a uniform, broadband phase shift, the configuration achieves an average difference in intensity transmission between the two arms of the interferometer of 0.016 percent with a standard deviation of 3.64 x 10(exp -4) percent, an average difference in polarization between the two arms of the interferometer of 5.47 x 10(exp -5) percent with a standard deviation of 1.57 x 10(exp -6) percent, and an average chromatic beam shift between the two arms of the interferometer of -47.53 microns with a wavelength-by-wavelength spread of 0.389 microns.

Bolcar, Matthew R.↗

Position Estimation Using Image Derivative

This paper describes an image processing algorithm to process Moon and/or Earth images. The theory presented is based on the fact that Moon hard edge points are characterized by the highest values of the image derivative. Outliers are eliminated by two sequential filters. Moon center and radius are then estimated by nonlinear least-squares using circular sigmoid functions. The proposed image processing has been applied and validated using real and synthetic Moon images.

Mortari, Daniele↗

A Formally Verified Conflict Detection Algorithm for Polynomial Trajectories

In air traffic management, conflict detection algorithms are used to determine whether or not aircraft are predicted to lose horizontal and vertical separation minima within a time interval assuming a trajectory model. In the case of linear trajectories, conflict detection algorithms have been proposed that are both sound, i.e., they detect all conflicts, and complete, i.e., they do not present false alarms. In general, for arbitrary nonlinear trajectory models, it is possible to define detection algorithms that are either sound or complete, but not both. This paper considers the case of nonlinear aircraft trajectory models based on polynomial functions. In particular, it proposes a conflict detection algorithm that precisely determines whether, given a lookahead time, two aircraft flying polynomial trajectories are in conflict. That is, it has been formally verified that, assuming that the aircraft trajectories are modeled as polynomial functions, the proposed algorithm is both sound and complete.

Narkawicz, Anthony↗

A Systematic Methodology for Modeling and Attitude Control of Multi-body Space Telescopes

This paper derives a symbolic multi-body rigid nonlinear model for a space telescope using Stoneking’s implementation of Kane’s method. This symbolic nonlinear model is linearized using Matlab symbolic functions diff and inv because the analytic linearization is intractable for manual derivation. The linearized system model is then used to design the controllers using both linear quadratic regulator (LQR) and robust pole assignment methods. The closed-loop systems for the two designs are simulated using both the rigid model as well as a second model containing flexible modes. The performances of the two designs are compared based on the simulation testing results. Our conclusion is that the robust pole assignment design offers better performance than that of the LQR system in terms of actuator usage and pointing accuracy. However, the LQR approach remains an effective first design step that can inform the selection of real eigenvalues for robust pole assignment. The proposed method may be used for the modeling and controller designs for various multi-body systems.

Modeling↗

Simulation studies for surfaces and materials strength

A realistic potential energy function comprising angle dependent terms was employed to describe the potential surface of the N+O2 system. The potential energy parameters were obtained from high level ab-initio results using a nonlinear fitting procedure. It was shown that the potential function is able to reproduce a large number of points on the potential surface with a small rms deviation. A literature survey was conducted to analyze exclusively the status of current small cluster research. This survey turned out to be quite useful in understanding and finding out the existing relationship between theoretical as well as experimental investigative techniques employed by different researchers. Additionally, the importance of the role played by computer simulation in small cluster research, was documented.

Halicioglu, Timur↗

A Method for Calculating Strain Energy Release Rates in Preliminary Design of Composite Skin/Stringer Debonding Under Multi-Axial Loading

Three simple procedures were developed to determine strain energy release rates, G, in composite skin/stringer specimens for various combinations of unaxial and biaxial (in-plane/out-of-plane) loading conditions. These procedures may be used for parametric design studies in such a way that only a few finite element computations will be necessary for a study of many load combinations. The results were compared with mixed mode strain energy release rates calculated directly from nonlinear two-dimensional plane-strain finite element analyses using the virtual crack closure technique. The first procedure involved solving three unknown parameters needed to determine the energy release rates. Good agreement was obtained when the external loads were used in the expression derived. This superposition technique was only applicable if the structure exhibits a linear load/deflection behavior. Consequently, a second technique was derived which was applicable in the case of nonlinear load/deformation behavior. The technique involved calculating six unknown parameters from a set of six simultaneous linear equations with data from six nonlinear analyses to determine the energy release rates. This procedure was not time efficient, and hence, less appealing. A third procedure was developed to calculate mixed mode energy release rates as a function of delamination lengths. This procedure required only one nonlinear finite element analysis of the specimen with a single delamination length to obtain a reference solution for the energy release rates and the scale factors. The delamination was extended in three separate linear models of the local area in the vicinity of the delamination subjected to unit loads to obtain the distribution of G with delamination lengths. This set of sub-problems was Although additional modeling effort is required to create the sub- models, this local technique is efficient for parametric studies.

Krueger, Ronald↗

Nonlinear saturation of whistler modes driven by runaway electrons

Runaway electrons exhibit kinetic instabilities with potentially beneficial consequences. The anomalous Doppler resonance between the electrons and whistler modes is a primary underlying mechanism. These instabilities require a first-principle nonlinear theoretical analysis, which would ultimately enable assessment of their impact on disruption mitigation and diagnostics, especially in ITER-relevant conditions. This paper presents recent progress in developing the required theoretical framework for isolated nonlinear resonances. By employing the generic bump-on-tail model, we predict the wave saturation levels in both near-threshold and strongly driven regimes. Remarkably, for the waves of interest, the parallel component of the wave vector dictates the parallel momentum of the resonant particles. This feature, together with the expected wave saturation level, provides a complete description of the resonance impact on the runaway electron distribution function. We show that the parametric dependence of the nonlinearly saturated mode provides a way to recover certain features of the runaway momentum distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of the comparison principle to analysis of nonlinear systems

A comparison principle based on a Kamke theorem and Lipschitz conditions is presented along with its possible applications and modifications. It is shown that the comparison lemma can be used in the study of such areas as classical stability theory, higher order trajectory derivatives, Liapunov functions, boundary value problems, approximate dynamic systems, linear and nonlinear systems, and bifurcation analysis.

Gunderson, R. W.↗

Global transformations of nonlinear systems

Necessary and sufficient conditions for a nonlinear system of equations to be locally equivalent, in a neighborhood of the origin in the real number system, to a controllable linear system are combined with several versions of the global inverse function theorem to define sufficient conditions for transforming the nonlinear system into a linear system. Additionally, a technique is introduced for developing a transformation under the assumptions that the columns of a controllability matrix span an n-dimensional space. Finally, the n-l form of the controllability matrix columns is demonstrated to be involutive

Hunt, L. R.↗

Nonlinear nanoelectrodynamics of a Weyl metal

Chiral Weyl fermions with linear energy-momentum dispersion in the bulk accompanied by Fermi-arc states on the surfaces prompt a host of enticing optical effects. While new Weyl semimetal materials keep emerging, the available optical probes are limited. In particular, isolating bulk and surface electrodynamics in Weyl conductors remains a challenge. We devised an approach to the problem based on near-field photocurrent imaging at the nanoscale and applied this technique to a prototypical Weyl semimetal TaIrTe 4 . As a first step, we visualized nano-photocurrent patterns in real space and demonstrated their connection to bulk nonlinear conductivity tensors through extensive modeling augmented with density functional theory calculations. Notably, our nanoscale probe gives access to not only the in-plane but also the out-of-plane electric fields so that it is feasible to interrogate all allowed nonlinear tensors including those that remained dormant in conventional far-field optics. Surface- and bulk-related nonlinear contributions are distinguished through their “symmetry fingerprints” in the photocurrent maps. Robust photocurrents also appear at mirror-symmetry breaking edges of TaIrTe 4 single crystals that we assign to nonlinear conductivity tensors forbidden in the bulk. Here, nano-photocurrent spectroscopy at the boundary reveals a strong resonance structure absent in the interior of the sample, providing evidence for elusive surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Uncertainty propagation in feed-forward neural network models

We develop new uncertainty propagation methods for feed-forward neural network architectures with leaky ReLU activation functions subject to random perturbations in the input vectors. In particular, we derive analytical expressions for the probability density function (PDF) of the neural network output and its statistical moments as a function of the input uncertainty and the parameters of the network, i.e., weights and biases. A key finding is that an appropriate linearization of the leaky ReLU activation function yields accurate statistical results even for large perturbations in the input vectors. This can be attributed to the way information propagates through the network. We also propose new analytically tractable Gaussian copula surrogate models to approximate the full joint PDF of the neural network output. To validate our theoretical results, we conduct Monte Carlo simulations and a thorough error analysis on a multi-layer neural network representing a nonlinear integro-differential operator between two polynomial function spaces. Our findings demonstrate excellent agreement between the theoretical predictions and Monte Carlo simulations.

MLP networks↗

Paraxial, Thin-Lens Analysis of Fixed-Tune, Non-Scaling FFAs with Two Magnets per Cell

Fixed-field accelerators (FFAs) have closed orbits that change as a function of beam momentum. It is sometimes useful to avoid various resonances by keeping the tunes constant as these orbits change. A well-known example is the scaling FFA where the entire beam orbit and optics are geometrically scaled as a function of momentum. However, this is stricter than necessary: cells with three or more lenses can have fixed tunes even when the focussing strengths in the lenses change. Recently, Dejan Trbojevic has found a pair of nonlinear magnets that produce fixed tunes and fixed beta functions while not being a scaling FFA. This improves on a more approximate solution in. Notably a scaling FFA would require one magnet to be entirely reverse-bending, whereas the Trbojevic solution does not do this and instead resembles an intermediate point between the traditional scaling field profiles and the nonscaling profiles centred around a momentum with equal and positive fields. This suggests there are at least three levels of stringency that one can apply to a fixed-tune FFA design: 1. Fixed cell tunes (in both planes) as a function of momentum; 2. Fixed optics (beta functions) as a function of momentum; 3. Similarity of all orbits via a scaling symmetry law ↔ traditional scaling FFA. This note studies the interesting case #2 above (#3 being fully characterised by the orbit at a single energy) in the simplest possible example: a cell of two thin lenses in the small angle (paraxial) approximation.

43 PARTICLE ACCELERATORS↗

Equilibrium thermodynamics and stochastic nonlinear acoustic fields

A crystalline solid is considered to consist of a large number of incoherent nonlinear acoustic radiation sources identified with the vibrating particles of the crystalline lattice. Randomization of the field, together with the assumption of a stochastically independent, fluctuating, radiation field at the absolue zero of temperature, leads to an expression of the temperature-dependent radiation field in terms of the zero-point field. The equation is identified with the Planck distribution formula of quantum mechanics in the linear field limit. The thermodynamic state functions are also obtained in terms of the nonlinear acoustic modal energies per unit mass and reduce to the results of the Debye-Einstein stochastic quantum oscillator model in the linear field limit.

Cantrell, J. H.↗

The 1996 Mid-Atlantic Winter Flood: Exploring Climate Risk through a Storyline Approach

Abstract This article explores the application of thermodynamic perturbations to a historical midlatitude, wintertime, rain-on-snow flood event to evaluate how similar events may evolve under different climate forcings. In particular, we generate a hindcast of the 1996 Mid-Atlantic flood using an ensemble of 14-km variable-resolution simulations completed with the U.S. Department of Energy’s global Energy Exascale Earth System Model (E3SM). We show that the event is skillfully reproduced over the Susquehanna River Basin (SRB) by E3SM when benchmarked against in situ observational data and high-resolution reanalyses. In addition, we perform five counterfactual experiments to simulate the flood under preindustrial conditions and four different levels of warming as projected by the Community Earth System Model Large Ensemble. We find a nonlinear response in simulated surface runoff and streamflow as a function of atmospheric warming. This is attributed to changing contributions of liquid water input from a shallower initial snowpack (decreased snowmelt), increased surface temperatures and rainfall rates, and increased soil water storage. Flooding associated with this event peaks from around +1 to +2 K of global average surface warming and decreases with additional warming beyond this. There are noticeable timing shifts in peak runoff and streamflow associated with changes in the flashiness of the event. This work highlights the utility of using storyline approaches for communicating climate risk and demonstrates the potential nonlinearities associated with hydrologic extremes in areas that experience ephemeral snowpack, such as the SRB.

Meteorology & Atmospheric Sciences↗

Real-Time Onboard Global Nonlinear Aerodynamic Modeling from Flight Data

Flight test and modeling techniques were developed to accurately identify global nonlinear aerodynamic models onboard an aircraft. The techniques were developed and demonstrated during piloted flight testing of an Aermacchi MB-326M Impala jet aircraft. Advanced piloting techniques and nonlinear modeling techniques based on fuzzy logic and multivariate orthogonal function methods were implemented with efficient onboard calculations and flight operations to achieve real-time maneuver monitoring and analysis, and near-real-time global nonlinear aerodynamic modeling and prediction validation testing in flight. Results demonstrated that global nonlinear aerodynamic models for a large portion of the flight envelope were identified rapidly and accurately using piloted flight test maneuvers during a single flight, with the final identified and validated models available before the aircraft landed.

Brandon, Jay M.↗

Compositional dependence of linear and nonlinear optical response in crystalline hafnium zirconium oxide thin films

Composition dependence of second harmonic generation, refractive index, extinction coefficient, and optical bandgap in 20 nm thick crystalline Hf 1– x Zr x O 2 (0 ≤ x ≤ 1) thin films is reported here. The refractive index exhibits a general increase with increasing ZrO 2 content with all values within the range of 1.98–2.14 from 880 nm to 400 nm wavelengths. A composition dependence of the indirect optical bandgap is observed, decreasing from 5.81 eV for HfO 2 to 5.17 eV for Hf 0.4 Zr 0.6 O 2 . The bandgap increases for compositions with x > 0.6, reaching 5.31 eV for Hf 0.1 Zr 0.9 O 2 . Second harmonic signals are measured for 880 nm incident light. The magnitude of the second harmonic signal scales with the magnitude of the remanant polarization in the composition series. Film compositions that display near zero remanent polarizations exhibit minimal second harmonic generation while those with maximum remanent polarization also display the largest second harmonic signal. The results are discussed in the context of ferroelectric phase assemblage in the hafnium zirconium oxide films and demonstrate a path toward a silicon-compatible integrated nonlinear optical material.

36 MATERIALS SCIENCE↗

Linear and nonlinear ARMA model parameter estimation using an artificial neural network

This paper addresses parametric system identification of linear and nonlinear dynamic systems by analysis of the input and output signals. Specifically, we investigate the relationship between estimation of the system using a feedforward neural network model and estimation of the system by use of linear and nonlinear autoregressive moving-average (ARMA) models. By utilizing a neural network model incorporating a polynomial activation function, we show the equivalence of the artificial neural network to the linear and nonlinear ARMA models. We compare the parameterization of the estimated system using the neural network and ARMA approaches by utilizing data generated by means of computer simulations. Specifically, we show that the parameters of a simulated ARMA system can be obtained from the neural network analysis of the simulated data or by conventional least squares ARMA analysis. The feasibility of applying neural networks with polynomial activation functions to the analysis of experimental data is explored by application to measurements of heart rate (HR) and instantaneous lung volume (ILV) fluctuations.

NASA Discipline Regulatory Physiology↗