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 1,801 records · Page 100

On the nonlinear versus linearized lower boundary conditions for topographically forced stationary long waves

For quasi-geostraphic stationary long waves forced by topography, the nonlinear lower boundary condition is derived in terms of the geopotential height and compared with the linearized version. The common practice of replacing terms describing the flow over and around a mountain by upstream zonal flow over the mountain and evaluating the resulting condition at sea level is found to be a good approximation for the cases considered and does not need to be modified as sometimes suggested. Specifically, it is found that this approximation does not affect, for most cases, the lower boundary condition expressed in terms of the geopotential height provided that the stationary wave is not near resonance. At resonance, the eddy advection terms may become important for large-amplitude waves when dissipation and surface diabatic heating are taken into account.

Tung, K. K.↗

A crystallographic model for nickel base single crystal alloys

The purpose of this research is to develop a tool for the mechanical analysis of nickel-base single-crystal superalloys, specifically Rene N4, used in gas turbine engine components. This objective is achieved by developing a rate-dependent anisotropic constitutive model and implementing it in a nonlinear three-dimensional finite-element code. The constitutive model is developed from metallurgical concepts utilizing a crystallographic approach. An extension of Schmid's law is combined with the Bodner-Partom equations to model the inelastic tension/compression asymmetry and orientation-dependence in octahedral slip. Schmid's law is used to approximate the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response and strain-rate sensitivity of the single-crystal superalloys. Methods for deriving the material constants from standard tests are also discussed. The model is implemented in a finite-element code, and the computed and experimental results are compared for several orientations and loading conditions.

Dame, L. T.↗

Crossflow-vortex instability and transition on a 45 deg swept wing

A crossflow vortex experiment on a 45 deg swept wing is currently being conducted in the Arizona State University Unsteady Wind Tunnel. The experimental apparatus is designed to produce crossflow-dominated transition by simulating infinite swept wing flow using contoured end liners in a closed throat wind tunnel. Stationary fixed-wavelength crossflow vortices are observed at several chord Reynolds numbers. The vortex wavelength which is fixed for a given Reynolds number varies with Reynolds number approximately as predicted by linear stability theory, but with the predicted wavelengths about 30 percent larger than the observed wavelengths. Travelling waves are observed both in the frequency range predicted by linear stability theory and at higher frequencies. These higher frequency waves may be harmonics of the primary crossflow waves generated by a nonlinear parametric resonance phenomena.

Dagenhart, J. Ray↗

Resonant decay of global MHD modes at 'thick' interfaces

A global surface mode can exist on a 'surface' which is not a true discontinuity, but without sufficient dissipation it is not a normal mode and it decays in time via resonance absorption. The decay rate can be calculated analytically when the 'surface' is thin. The goal of this paper is to numerically estimate the decay rate when the surface is not thin. A cold plasma in a uniform magnetic field is considered in which density varies linearly across the 'surface'. In this linearized calculation, the global surface mode is driven in steady state by an antenna located in one of the uniform regions external to the surface. The frequency is a free parameter, and resonance curves are computed numerically without approximations. The widths of the resonance curves are used to estimate the free decay times of undriven surface modes, via the uncertainty principle. When applied to the active solar corona, the decay rates are large enough to account for the coronal heating, but it should be kept in mind that the role of nonlinearity in resonance absorption is still undetermined.

Hollweg, Joseph V.↗

Effects of microgravity on the interaction of vestibular and optokinetic nystagmus in the vertical plane

The extent to which the slow phase velocity (SPV) of nystagmus elicited by a vertical optokinetic stimulation with constant velocity could be modulated by sinusoidal angular motion in the vertical plane was investigated under normal gravity condition and during the microgravity period of parabolic flight. In normal gravity, when the angular head motion and the optokinetic stimulation were in the same direction, the peak SPV was slower than the velocity of the optokinetic display. When the head motion and the optokinetic stimulation were in opposite directions, the peak SPV was equal to the velocity of the optokinetic display. In microgravity, the peak SPV was approximately equal to the velocity of the optokinetic display when head rotation and optokinetic stimulation were in the same direction, and was faster than the velocity of the optokinetic dispaly when head rotation and optokinetic stimulation were in opposite directions. In addition, the interaction of vestibular and optokinetic nystagmus was found to be nonlinear in microgravity, especially when the optokinetic stimulation was directed downward. These results suggest an interaction between the vestibular and the optokinetic systems modulated as a function of the gravitational state, and support the observation that visual input is more effective in reducing sensory conflict experienced in microgravity.

Clement, Gilles↗

An intercomparison of plasma turbulence at three comets: Grigg-Skjellerup, Giacobini-Zinner, and Halley

We examine and intercompare the LF plasma wave turbulence at three comets: Grigg-Skjellerup (GS), Giacobini-Zinner (GZ), and Halley (H). All three have power spectral peaks at the local ion cyclotron frequency (the pump wave) at approx. 10(exp -2) Hz, and a power-law fall-off at higher frequencies that suggest the development of turbulent cascades. The power laws for the three comets are approximately f(exp -1.9), f(exp -1.9) and f(exp -2.1), respectively. However, other than the similarities in the power spectra, we find the magnetic field turbulence is considerably different at the three comets. Phase steepening is demonstrated to occur at the trailing edges of the GS waves. This is probably due to nonlinear steepening plus dispersion of the left-hand mode components. A coherency analysis of GZ turbulence indicates that it is primarily composed of righthanded mode components, i.e., the turbulence is 'whistlermode.' This too can be explained by nonlinear steepening plus dispersion of the magnetosonic waves. At the level of GS and GZ turbulence development when the spacecraft measurements were made, classical three-wave processes, such as the decay or modulation instabilities do not appear to play important roles. It is most likely that the nonlinear steepening and dispersive time scales are more rapid than three-wave processes, and the latter had not had time to develop for the relatively 'new' turbulence. The wave turbulence at Halley is linearly polarized. The exact nature of this turbulence is still not well understood at this time. Several possibilities are suggested, based on our preliminary analyses.

Tsurutani, Bruce T.↗

Kinetic fluid moments closure for a magnetized plasma with collisions

A novel method aimed at a kinetic moments closure for a magnetized plasma with collisions is proposed. The velocity distribution function for each species is expanded in 8 Gaussian Radial Basis Functions (GRBFs) which are essentially shifted Maxwellians at eight representative 3D-velocity points of drift. The vector of 8 fluid moments (for particle density, 3 particle fluxes, total energy density, and 3 energy fluxes) has an 8x8 analytic linear matrix relation to the vector of 8 GRBF density weighs in 3D-real space. The 8 fluid moments with sources for each species are advanced in time while the 8 GRBF weighs are determined from the 8x8 inverse matrix. The two closure moments (for the stress tensor and the energy weighted stress tensor) are linearly determined from the GRBF weighs. Most importantly the velocity moments of the nonlinear Coulomb Fokker-Planck collision operator [Rosenbluth et al, Phys. Rev 107, 1957] are evaluated from the GRBF weights. Generalization from 8 to 12, 16, 20 .., in an energy weighted moment hierarchy is straightforward. The electric field follows from a generalized vorticity (quasi-neutral charge conservation) equation. A strong drift ordering approximation can be applied to eliminate any spuriously unstable high frequency cyclotron motions. A novel weak drift ordering two-time step scheme avoids the vorticity equation by following ion cyclotron motion in time to get the electric field with ion gyroaveraging. Inclusion of low-beta magnetic perturbations is straightforward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Elastic Scattering of Ultrashort Laser Pulses on Noble Gas Atoms

The work deals with a theoretical investigation of elastic scattering of ultrashort electromagnetic pulses (USP) on the noble gas atoms in the far ultra-violet and soft X-ray regions. The investigation is performed within the framework of a dipole approximation and the first-order quantum-mechanics perturbation theory. Major attention is given to the study of the dependence of the total probability of scattering within the pulse time on the USP duration at different carrying frequency values. It is shown that at certain values of the problem parameters this dependence is nonlinear. The calculations are performed for the corrected Gaussian pulses and wavelet pulses without the carrying frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enhancement of the sensitivity of the light-induced drift effect to interatomic interaction potentials using a mixture of two buffer gases

Based on five different nonempirical (calculated ab initio) interaction potentials for pairs of colliding Li – Ne particles and three different interaction potentials for pairs of colliding Li – Ar particles, we have theoretically investigated the spectral features of the light-induced drift (LID) rate of Li atoms in the buffer Ne gas in a mixture of buffer Ne + Ar gases. The calculations of LID of Li atoms in the buffer Ne gas for two interaction potentials predict anomalous LID and, as a result, strong sensitivity of the spectral shape of the LID lines to the differences in these potentials. For three other potentials (out of the five in question), the shape of the LID line of Li atoms in Ne is insensitive to the shape of the potential, since calculations with these potentials predict a normal LID effect. In this case, as it turned out, by adding a small fraction (approximately 10 %) of Ar to Ne, one can go from normal LID to anomalous LID and thereby radically increase the sensitivity of the LID line shape of Li atoms to the difference in these interaction potentials. The results obtained enable high-precision testing of the interatomic interaction potentials in experiments on anomalous LID. (nonlinear optical phenomena)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generative Vulnerability Assessment for Cyber-Physical Systems

Cyber-physical systems (CPS) are highly susceptible to malicious attacks due to their complex dynamics and interconnectivity. A comprehensive understanding of their vulnerabilities is essential for designing effective resilience measures. This paper presents a data-driven attack generative system for evaluating the vulnerability of CPS. The proposed approach formulates the vulnerability assessment problem as determining the feasibility of a specific attack set based on two boundary functions that represent the effectiveness and stealthiness of attacks. The attack generative model is trained using a custom loss function, with two universal approximators designed to learn the effectiveness and stealthiness functions simultaneously. Theoretical results for successful generation and asymptotic convergence of the resulting training algorithm are given. As a result, the proposed approach is evaluated via numerical simulation of an IEEE 14-bus system and gas pipeline systems, demonstrating its viability in learning how to attack nonlinear CPS and identify potential vulnerabilities.

Computer systems organization↗

Asymptotic Behavior of Memristive Circuits

The interest in memristors has risen due to their possible application both as memory units and as computational devices in combination with CMOS. This is in part due to their nonlinear dynamics, and a strong dependence on the circuit topology. We provide evidence that also purely memristive circuits can be employed for computational purposes. In the present paper we show that a polynomial Lyapunov function in the memory parameters exists for the case of DC controlled memristors. Such a Lyapunov function can be asymptotically approximated with binary variables, and mapped to quadratic combinatorial optimization problems. This also shows a direct parallel between memristive circuits and the Hopfield-Little model. In the case of Erdos-Renyi random circuits, we show numerically that the distribution of the matrix elements of the projectors can be roughly approximated with a Gaussian distribution, and that it scales with the inverse square root of the number of elements. This provides an approximated but direct connection with the physics of disordered system and, in particular, of mean field spin glasses. Using this and the fact that the interaction is controlled by a projector operator on the loop space of the circuit. We estimate the number of stationary points of the approximate Lyapunov function and provide a scaling formula as an upper bound in terms of the circuit topology only.

97 MATHEMATICS AND COMPUTING↗

Radio observations of the planets - The importance of laboratory measurements

Laboratory data on the line broadening parameters of H2S gas under Uranian/Neptunian conditions, on the far wings of the H2S and NH3 line profiles, and on the dielectric properties of CH4-, NH3-, H2S-, and NH4SH-ice are needed to constrain elemental abundances and understand the dynamics and cloud physics in the atmospheres of the giant planets. Measurements of the absorption coefficient of gaseous H2SO4 at millimeter wavelengths are needed in order to obtain a better understanding of Venus' atmosphere. To determine wind velocity fields in Venus' and Mars' atmospheres, accurate measurements of the center frequencies of the CO lines are necessary. The absorption and scattering properties of lunar soils and/or terrestrial rock powders at frequencies from approximately 1 to 200 GHz, determined in laboratory experiments, would provide a valuable addition to existing data at 450 MHz, 35 GHz, and far infrared frequencies. These data would be used to analyze the microwave spectra of planetary surfaces. Such studies may be helpful in distinguishing the effects of radiative transfer from those of nonlinear heat conduction and internal heat sources.

De Pater, Imke↗

Dynamic analysis of Space Shuttle/RMS configuration using continuum approach

The initial assembly of Space Station Freedom involves the Space Shuttle, its Remote Manipulation System (RMS) and the evolving Space Station Freedom. The dynamics of this coupled system involves both the structural and the control system dynamics of each of these components. The modeling and analysis of such an assembly is made even more formidable by kinematic and joint nonlinearities. The current practice of modeling such flexible structures is to use finite element modeling in which the mass and interior dynamics is ignored between thousands of nodes, for each major component. The model characteristics of only tens of modes are kept out of thousands which are calculated. The components are then connected by approximating the boundary conditions and inserting the control system dynamics. In this paper continuum models are used instead of finite element models because of the improved accuracy, reduced number of model parameters, the avoidance of model order reduction, and the ability to represent the structural and control system dynamics in the same system of equations. Dynamic analysis of linear versions of the model is performed and compared with finite element model results. Additionally, the transfer matrix to continuum modeling is presented.

Ramakrishnan, Jayant↗

Hybrid Upwind Splitting (HUS) by a Field-by-Field Decomposition

We introduce and develop a new approach for upwind biasing: the hybrid upwind splitting (HUS) method. This original procedure is based on a suitable hybridization of current prominent flux vector splitting (FVS) and flux difference splitting (FDS) methods. The HUS method is designed to naturally combine the respective strengths of the above methods while excluding their main deficiencies. Specifically, the HUS strategy yields a family of upwind methods that exhibit the robustness of FVS schemes in the capture of nonlinear waves and the accuracy of some FDS schemes in the resolution of linear waves. We give a detailed construction of the HUS methods following a general and systematic procedure directly performed at the basic level of the field by field (i.e. waves) decomposition involved in FDS methods. For such a given decomposition, each field is endowed either with FVS or FDS numerical fluxes, depending on the nonlinear nature of the field under consideration. Such a design principle is made possible thanks to the introduction of a convenient formalism that provides us with a unified framework for upwind methods. The HUS methods we propose bring significant improvements over current methods in terms of accuracy and robustness. They yield entropy-satisfying approximate solutions as they are strongly supported in numerical experiments. Field by field hybrid numerical fluxes also achieve fairly simple and explicit expressions and hence require a computational effort between that of the FVS and FDS. Several numerical experiments ranging from stiff 1D shock-tube to high speed viscous flows problems are displayed, intending to illustrate the benefits of the present approach. We assess in particular the relevance of our HUS schemes to viscous flow calculations.

Coquel, Frederic↗

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE↗

Train Like a (Var)Pro: Efficient Training of Neural Networks with Variable Projection

Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.

97 MATHEMATICS AND COMPUTING↗

Halo Profiles and the Concentration–Mass Relation for a ΛCDM Universe

Profiles of dark matter-dominated halos at the group and cluster scales play an important role in modern cosmology. Using results from two very large cosmological N-body simulations, which increase the available volume at their mass resolution by roughly two orders of magnitude, we robustly determine the halo concentration-mass (c-M) relation over a wide range of masses, employing multiple methods of concentration measurement. We then characterize individual halo profiles, as well as stacked profiles, relevant for galaxy-galaxy lensing and next-generation cluster surveys; the redshift range covered is 0 ≤ z ≤ 4, with a minimum halo mass of M 200c ~ 2 x 10 11 M ⊙ . Despite the complexity of a proper description of a halo (environmental effects, merger history, nonsphericity, relaxation state), when the mass is scaled by the nonlinear mass scale M * (z), we find that a simple non-power-law form for the c-M/M * relation provides an excellent description of our simulation results across eight decades in M/M * and for 0 ≤ z ≤ 4. Over the mass range covered, the c-M relation has two asymptotic forms: an approximate power law below a mass threshold M/M * ~500-1000, transitioning to a constant value, c 0 ~ 3 at higher masses. The relaxed halo fraction decreases with mass, transitioning to a constant value of ~0.5 above the same mass threshold. We compare Navarro-Frenk-White (NFW) and Einasto fits to stacked profiles in narrow mass bins at different redshifts; as expected, the Einasto profile provides a better description of the simulation results. At cluster scales at low redshift, however, both NFW and Einasto profiles are in very good agreement with the simulation results, consistent with recent weak lensing observations.

79 ASTRONOMY AND ASTROPHYSICS↗

A method for digital image registration using a mathematical programming technique

A new algorithm based on a nonlinear programming technique to correct the geometrical distortions of one digital image with respect to another is discussed. This algorithm promises to be superior to existing ones in that it is capable of treating localized differential scaling, translational and rotational errors over the whole image plane. A series of piece-wise 'rubber-sheet' approximations are used, constrained in such a manner that a smooth approximation over the entire image can be obtained. The theoretical derivation is included. The result of using the algorithm to register four channel S065 Apollo IX digitized photography over Imperial Valley, California, is discussed in detail.

Yao, S. S.↗

Refine your search to explore more results.