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,765 records · Page 98

Experimental observation of linear and nonlinear ion acoustic phenomena in a cylindrical geometry

Ion acoustic phenomena are studied in a cylindrical geometry for two distinct cases. A large amplitude compressive pulse is seen to evolve into solitons. The evolution of these solitons and their dependence on initial conditions show a similarity to previous work on one dimensional solitons. Dimensionless scaling arguments are used to distinguish the two cases. In the presence of a steady state uniform cylindrical beam, approximated by a ring in V sub r, V sub phi, an ion-ion beam instability is observed. This instability exists for a limited range of beam velocities and shows a marked similarity to the strictly one dimensional ion-ion beam instability. Solution of the appropriate dispersion relation shows agreement with the observed phenomenon.

Romesser, T. E.↗

Microwave Heating of a Liquid Stably Flowing in a Circular Channel Under the Conditions of Nonstationary Radiative-Convective Heat Transfer

A class of nonlinear problems of nonstationary radiative-convective heat transfer under the microwave action with a small penetration depth is considered in a stabilized coolant flow in a circular channel. The solutions to these problems are obtained, using asymptotic procedures at the stages of nonstationary and stationary convective heat transfer on the heat-radiating channel surface. The nonstationary and stationary stages of the solution are matched, using the 'longitudinal coordinate–time' characteristic. The approximate solutions constructed on such principles correlate reliably with the exact ones at the limiting values of the operation parameters, as well as with numerical and experimental data of other researchers. An important advantage of these solutions is that they allow the determination of the main regularities of the microwave and thermal radiation influence on convective heat transfer in a channel even before performing cumbersome calculations. It is shown that, irrespective of the heat exchange regime (nonstationary or stationary), the Nusselt number decreases and the rate of the surface temperature change increases with increase in the intensity of thermal action.

42 ENGINEERING↗

On the nonlinearity of spatial scales in extreme weather attribution statements

In the context of continuing climate change, extreme weather events are drawing increasing attention from the public and news media. A question often asked is how the likelihood of extremes might have changed by anthropogenic greenhouse-gas emissions. Answers to the question are strongly influenced by the model used, duration, spatial extent, and geographic location of the event—some of these factors often overlooked. Using output from four global climate models, we provide attribution statements characterised by a change in probability of occurrence due to anthropogenic greenhouse-gas emissions, for rainfall and temperature extremes occurring at seven discretised spatial scales and three temporal scales. An understanding of the sensitivity of attribution statements to a range of spatial and temporal scales of extremes allows for the scaling of attribution statements, rendering them relevant to other extremes having similar but non-identical characteristics. This is a procedure simple enough to approximate timely estimates of the anthropogenic contribution to the event probability. Furthermore, since real extremes do not have well-defined physical borders, scaling can help quantify uncertainty around attribution results due to uncertainty around the event definition. Results suggest that the sensitivity of attribution statements to spatial scale is similar across models and that the sensitivity of attribution statements to the model used is often greater than the sensitivity to a doubling or halving of the spatial scale of the event. The use of a range of spatial scales allows us to identify a nonlinear relationship between the spatial scale of the event studied and the attribution statement.

54 ENVIRONMENTAL SCIENCES↗

Construction of solutions for some nonlinear two-point boundary value problems

Constructive existence and uniqueness results for boundary value problems associated with some simple special cases of the second order equation y'' = f(x,y,y') 0 or = x or = 1, are sought. The approach considered is to convert the differential equation and boundary conditions to an integral equation via Green's functions, and then to apply fixed point and contraction map principles to a sequence of successive approximations. The approach is tested on several applied problems. Difficulties in trying to prove general theorems are discussed.

Pennline, J. A.↗

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Designing stellarator coils by a modified Newton method using FOCUS

To find the optimal coils for stellarators, nonlinear optimization algorithms are applied in existing coil design codes. However, none of these codes have used the information from the second-order derivatives. In this paper, we present a modified Newton method in the recently developed code FOCUS. The Hessian matrix is calculated with analytically derived equations. Its inverse is approximated by a modified Cholesky factorization and applied in the iterative scheme of a classical Newton method. Using this method, FOCUS is able to recover the W7-X modular coils starting from a simple initial guess. Results demonstrate significant advantages.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uncovering extreme nonlinear dynamics in solids through time-domain field analysis

Time-domain analysis of harmonic fields with sub-cycle resolution is now experimentally viable due to the emergence of sensitive, on-chip techniques for petahertz-scale optical-field sampling. We demonstrate how such a time-domain, field-resolved analysis uncovers the extreme nonlinear electron dynamics responsible for high-harmonic generation within solids. Time-dependent density functional theory was used to simulate harmonic generation from a solid-state band-gap system driven by near- to mid-infrared waveforms. Particular attention was paid to regimes where both intraband and interband emission mechanisms play a critical role in shaping the nonlinear response. We show that a time-domain analysis of the harmonic radiation fields identifies the interplay between intra- and interband dynamical processes underlying the nonlinear light generation. With further analysis, we show that changes to the dominant emission regime can occur after only slight changes to the peak driving intensity and central driving wavelength. Furthermore, time-domain analysis of harmonic fields also reveals, for the first time, the possibility of rapid changes in the dominant emission mechanism within the temporal window of the driving pulse envelope. Finally, we examine the experimental viability of performing time-domain analysis of harmonic fields with sub-cycle resolution using realistic parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗

Average nonlinear dynamics of particles in gravitational pulses: Effective Hamiltonian, secular acceleration, and gravitational susceptibility

Particles interacting with a prescribed quasimonochromatic gravitational wave (GW) exhibit secular (average) nonlinear dynamics that can be described by Hamilton’s equations. We derive the Hamiltonian of this “ponderomotive” dynamics to the second order in the GW amplitude for a general background metric. For the special case of vacuum GWs, we show that our Hamiltonian is equivalent to that of a free particle in an effective metric, which we calculate explicitly. We also show that already a linear GW pulse displaces a particle from its unperturbed trajectory by a finite distance that is independent of the GW phase and proportional to the integral of the pulse intensity.This effect is independent from the nonlinear memory effects that has been known. We calculate the particle displacement analytically and show that our result is in agreement with numerical simulations. We also show how the Hamiltonian of the nonlinear averaged dynamics naturally leads to the concept of the linear gravitational susceptibility of a particle gas with an arbitrary phase-space distribution. We calculate this susceptibility explicitly to apply it, in a follow-up paper, toward studying self-consistent GWs in inhomogeneous media within the geometrical-optics approximation

79 ASTRONOMY AND ASTROPHYSICS↗

Uncertainty propagation and sensitivity analysis for constrained optimization of nuclear waste vitrification

Abstract The vitrification of high‐level waste (HLW) by heating a mixture of glass‐forming chemicals (GFCs) with the waste can be improved using a constrained optimization problem. This study explores how different uncertainty propagation (UP) methods implemented with the optimization process can affect the glass formulation of nuclear waste glasses. UP is the effort of propagating uncertain inputs through a system to understand and quantify output distributions. Uncertainty intervals are crafted from output distributions to inform the optimization algorithm. UP is often implemented with Monte Carlo (MC) sampling for large nonlinear systems, which can be difficult to implement within a constrained optimization algorithm that requires derivative information. Other UP methods often used for optimization under uncertainty (OUU) can be designed to work within an established constrained optimization framework. Methods of UP are evaluated in this study including iterative sampling approaches, first‐order approximations, and surrogate modeling with machine learning (ML). A method of dimensional reduction based on global sensitivity analysis is introduced to support the UP methods for the large dimensionality of the problem. Analytical UP methods able to achieve similar optimums 10 times faster than the baseline MC approach, and produce 93.9% similar output distributions are reported.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Implicit approximate-factorization schemes for the efficient solution of steady transonic flow problems

Implicit approximate-factorization algorithms (AF) are developed for the solution of steady-state transonic flow problems. The performance of the AF solution method is evaluated relative to that of the standard solution method for transonic flow problems, successive line over-relaxation (SLOR). Both methods are applied to the solution of the nonlinear, two-dimensional transonic small-disturbance equation. Results indicate that the AF method requires substantially less computer time than SLOR to solve the nonlinear finite-difference matrix equation for a transonic flow field. This increase in computational efficiency is achieved with no appreciable increase in computer storage or coding complexity.

Ballhaus, W. F.↗

Navigation simulator for the Space Tug vehicle

A general simulation program (GSP) for state estimation of a nonlinear space vehicle flight navigation system is developed and used as a basis for evaluating the performance of a Space Tug navigation system. An explanation of the iterative guidance mode (IGM) guidance law, derivation of the dynamics, coordinate frames and state estimation routines are given in order to clarify the assumptions and approximations made. A number of simulation and analytical studies are used to demonstrate the operation of the Tug system. Included in the simulation studies are (1) initial offset vector parameter study; (2) propagation time vs accuracy; (3) measurement noise parametric study and (4) reduction in computational burden of an on-board implementable scheme. From the results of these studies, conclusions and recommendations concerning future areas of practical and theoretical work are presented.

Colburn, B. K.↗

Flight through thunderstorm outflows

Computer simulation of aircraft landing through thunderstorm gust fronts is carried out. The two-dimensional, nonlinear equations or aircraft motion containing all wind shear terms are solved numerically. The gust front spatial wind field inputs are provided in the form of tabulated experimental data which are coupled with a computer table lookup routine to provide the required wind components and shear at any given position within an approximate 500 m by 1 km vertical plane. The aircraft is considered to enter the wind field at a specified position under trimmed conditions. Both fixed control and automatic control landings are simulated. Flight paths, as well as control inputs necessary to maintain specified trajectories, are presented and discussed for aircraft having characteristics of a DC-8, B-747, augmentor-wing STOL, and a DHC-6.

Frost, W.↗

Internal gravity waves in the solar atmosphere. II - Effects of radiative damping

In the solar photosphere, temperature fluctuations associated with acoustic-gravity waves may be rapidly smoothed by the transfer mechanism of radiation between hotter and cooler regions. The present investigation of the radiative effects on internal gravity waves takes into account the parameterization of the radiative energy, employing the Newtonian cooling approximation. A linear analysis of the propagation of internal gravity waves is carried out in a model of the solar atmosphere which is taken to be homogeneous in the horizontal coordinates. Linear wave properties both with and without radiative cooling are summarized, and the variation with height of energy fluxes and of nonlinearities in the waves is discussed. Attention is given to the significance of the obtained results in terms of energy balance in the chromosphere and in relation to spectral line observations.

Mihalas, B. W.↗

Development of an electronically-scanned pressure module for operation at cryogenic temperatures

Pressure and temperature characteristics were measured for a number of multichannel electronically scanned pressure sensors. The tests were made on commercially available units designed to operate in a controlled temperature environment. Measurements of zero shift, sensitivity, and nonlinearity for each transducer were taken over a temperature range from 100 K to 340 K using a computer controlled data acquisition system. The units tested failed to meet accuracy specifications over the complete temperature range, which was expected. However, the sensors showed acceptable and predictable behavior over the temperature range from approximately -40 C (233 K) to 70 C (343 K). It was determined that a combination of local heating and accurate temperature monitoring can result in a device that can be compensated for temperature as well as its other physical properties. The design of a prototype for operation in a cryogenic environment is proposed, and a method for compensation is developed.

Johnson, Preston B.↗

Estimating short-period dynamics using an extended Kalman filter

An extended Kalman filter (EKF) is used to estimate the parameters of a low-order model from aircraft transient response data. The low-order model is a state space model derived from the short-period approximation of the longitudinal aircraft dynamics. The model corresponds to the pitch rate to stick force transfer function currently used in flying qualities analysis. Because of the model chosen, handling qualities information is also obtained. The parameters are estimated from flight data as well as from a six-degree-of-freedom, nonlinear simulation of the aircraft. These two estimates are then compared and the discrepancies noted. The low-order model is able to satisfactorily match both flight data and simulation data from a high-order computer simulation. The parameters obtained from the EKF analysis of flight data are compared to those obtained using frequency response analysis of the flight data. Time delays and damping ratios are compared and are in agreement. This technique demonstrates the potential to determine, in near real time, the extent of differences between computer models and the actual aircraft. Precise knowledge of these differences can help to determine the flying qualities of a test aircraft and lead to more efficient envelope expansion.

Bauer, Jeffrey E.↗

A locally refined rectangular grid finite element method - Application to computational fluid dynamics and computational physics

The present FEM technique addresses both linear and nonlinear boundary value problems encountered in computational physics by handling general three-dimensional regions, boundary conditions, and material properties. The box finite elements used are defined by a Cartesian grid independent of the boundary definition, and local refinements proceed by dividing a given box element into eight subelements. Discretization employs trilinear approximations on the box elements; special element stiffness matrices are included for boxes cut by any boundary surface. Illustrative results are presented for representative aerodynamics problems involving up to 400,000 elements.

Young, David P.↗

Reduced basis technique for calculating sensitivity coefficients of nonlinear structural response

An efficient reduced basis technique is presented for calculating the sensitivity of nonlinear structural response to variations in the design variables. The structure is discretized by using two-field mixed finite element models. The vector of structural response and its sensitivity coefficients (derivatives with respect to design variables) are each expressed as a linear combination of a small number of basis (or global approximation) vectors. The Bubnov-Galerkin technique is then used to approximate each of the finite element equations governing the response and the sensitivity coefficients by a small number of algebraic equations in the amplitudes of these vectors. The path derivatives (derivatives of the response vector with respect to path parameters, e.g., load parameters) are used as basis vectors for approximating the response. A combination of the path derivatives and their derivatives with respect to the design variables is used for approximating the sensitivity coefficients. The potential of the proposed technique is discussed and its effectiveness is demonstrated by means of numerical examples of laminated composite plates subjected to mechanical and thermal loads.

Noor, Ahmed K.↗

Acoustic Analysis of Sensor Ports

Acquiring data in rocket engines that is representative of the actual dynamic environment can often be difficult due to a multitude of influences. One source of contamination that is often not considered entirely is the response associated with the acoustic cavity created by a sensor offset. It is common to offset a sensor due to various reasons such as for mounting accessibility, thermal isolation, shock reduction, or prevention of debris impingement. While estimating the natural frequency of the acoustic cavity is straightforward, limited analysis has been described on the determination of the overall frequency response. The sensor port design approach usually attempts to ensure the port is short enough such that the acoustic response is negligible near the frequency-of-interest, but this requires knowledge of the frequency response and simple rules-of-thumb are not always guaranteed. Data correction and/or data interpretation is also often desired for an unsatisfactory response. The limited response analysis in the literature only offers approximations or neglects important contributions. A new approach is devised theoretically and computationally that captures the true acoustic response of a sensor port. This paper summarizes the acoustics background, the port response theoretical development, and provides comparisons of a port acoustic response using an analytical model and computational acoustics. The effects of nonlinear acoustics are also examined. Additionally, the paper summarizes the design of a specialized filter using the predicted sensor port response that can be applied to data for correction.

Casiano, M. J.↗