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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

The Effect of Viscosity on the Spherical Stability of Oscillating Gas Bubbles

Gas bubbles driven in radial oscillations are subject to an instability of the spherical shape that is opposed by surface tension and viscosity. An exact linear formulation for the study of the phenomenon has been available for many years, but its complexity has discouraged a detailed investigation. With the recent theory of sonoluminescence of Lohse and co-workers, there has been a renewed interest in the problem and new data have become available. This paper presents a numerical method for the solution of the pertinent equations and compares the theory with these new data. The coupling of the strong nonlinearity of the bubble radial oscillations with the parametric mechanism of the surface instability results in a very complex structure for the stability boundary. Nevertheless, a good agreement between theory and data is found. A comparison with earlier approximate models is also made.

Hao, Y.↗
NASA NTRS1957

A Fuselage Addition to Delay Drag-Rise Mach Number of Subsonic Airplanes at Lifting Conditions

The addition of fuselage volume, concentrated on top of the forward portion of the fuselage, for the purpose of delaying the drag-rise Mach number of subsonic airplanes at lifting conditions is investigated. The additions have been designed on the basis of the area rule and other important considerations to provide greater practicability of application compared with shapings previously investigated. The addition delayed the drag-rise Mach number by an increment of approximately 0.03 for a configuration having a wing with moderate thickness and 35 deg of sweepback at a lift coefficient of 0.3. A lesser delay was obtained for a configuration with a thicker wing. The additions increase the nonlinearities of the variations of pitching moment with lift,

Whitcomb, Richard T.↗
DOE OSTI2017

Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration

Highlights: • Two numerical methods are presented to accelerate DMD simulations. • The methods greatly reduce the computation time without much loss of accuracy. • The accelerated DMD model is capable of simulating long-term mass transport at atomistic length scales. • Demonstrated for the propagation of an atomically sharp phase boundary over 30 seconds. - Abstract: Diffusive Molecular Dynamics (DMD) is a class of recently developed computational models for the simulation of long-term diffusive mass transport at atomistic length scales. Compared to previous atomistic models, e.g., transition state theory based accelerated molecular dynamics, DMD allows the use of larger time-step sizes, but has a higher computational complexity at each time-step due to the need to solve a nonlinear optimization problem at every time-step. This paper presents two numerical methods to accelerate DMD based simulations. First, we show that when a many-body potential function, e.g., embedded atom method (EAM), is employed, the cost of DMD is dominated by the computation of the mean of the potential function and its derivatives, which are high-dimensional random variables. To reduce the cost, we explore both first- and second-order mean field approximations. Specifically, we show that the first-order approximation, which uses a point estimate to calculate the mean, can reduce the cost by two to three orders of magnitude, but may introduce relatively large error in the solution. We show that adding an approximate second-order correction term can significantly reduce the error without much increase in computational cost. Second, we show that DMD can be significantly accelerated through subcycling time integration, as the cost of integrating the empirical diffusion equation is much lower than that of the optimization solver. To assess the DMD model and the numerical approximation methods, we present two groups of numerical experiments that simulate the diffusion of hydrogen in palladium nanoparticles. In particular, we show that the computational framework is capable of capturing the propagation of an atomically sharp phase boundary over a time window of more than 30 seconds. The effects of the proposed numerical methods on solution accuracy and computation time are also assessed quantitatively.

36 MATERIALS SCIENCE↗
NASA NTRS1976

A nonlinear theory of cosmic-ray pitch-angle diffusion in homogeneous magnetostatic turbulence

A plasma strong turbulence, weak coupling, theory is applied to the problem of cosmic-ray pitch-angle scattering in magnetostatic turbulence. The theory used is a rigorous generalization of Weinstock's 'resonance broadening' theory and contains no ad hoc approximations. A detailed calculation is presented for a model of 'slab' turbulence with an exponential correlation function. The results agree well with numerical simulations. The rigidity dependence of the pitch-angle scattering coefficient differs from that found by previous researchers. The differences result from an inadequate treatment of particle trajectories near 90 deg pitch angle in earlier work.

Goldstein, M. L.↗
DOE OSTI2022

Physics-Informed Sparse Gaussian Process for Probabilistic Stability Analysis of Large-Scale Power System with Dynamic PVs and Loads

This work proposes a physics-informed sparse Gaussian process (SGP) for probabilistic stability assessment of large-scale power systems in the presence of uncertain dynamic PVs and loads. The differential and algebraic equations considering uncertainties from dynamic PVs and loads are reformulated to a nonlinear mapping relationship that allows the application of SGP. Thanks to the nonparametric characteristic of Gaussian process, the proposed framework does not require distributions of uncertain inputs and this distinguishes it from existing approaches. As the original Gaussian process is not scalable to large-scale systems with high dimensional uncertain inputs, this paper develops the SGP with a stochastic variational inference technique. It leads to approximately two orders of complex reduction. A data pre-processing step is also introduced to tackle the coexistence of stable and unstable cases by sample clustering and constructing separate SGPs. The probabilistic transient stability index is analyzed to assess system stability under different uncertain dynamics loads and PVs. Comparisons are performed with the sampling-based, the polynomial chaos expansion-based, and traditional Gaussian process-based methods on the modified IEEE 118-bus and Texas 2000-bus systems under various scenarios, including different levels of uncertainties and the existence of nonlinear correlations among dynamic PVs. The impacts of data quality and quantity issues are also investigated. It is shown that the proposed SGP achieves significantly improved computational efficiency while maintaining high accuracy with a limited number of data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗
DOE OSTI2024

A simulation study of the ability to detect power distribution perturbations in the texas A&M TRIGA reactor with self-powered neutron detectors

Given the variety of ways that nuclear reactor core power may be perturbed, reactor operators and developers are keen on understanding the accuracy and convergence time during which perturbations in reactor power distribution may be synthesized (i.e., inferred) from an array of in-core radiation detectors. A simulation study was conducted as described herein using a highly detailed model of the Texas A&M Training, Research, Isotopes, General Atomics Reactor, in which an array of self-powered neutron detectors (SPNDs) was considered for input to the power synthesis methodology. The core power synthesis is conducted using a point-based iterative method with an iterative loop built in to ensure working equation consistency. The forward problem of SPND response to simulated perturbations in reactor power was solved for Gaussian peak-type perturbations in the reactor power distribution. These perturbations varied in variance, amplitude, and core location to assess their impact on synthesis error and to determine the number of iterations required for convergence. A relation between the unique resolvability limit and perturbation width was identified such that the maximum synthesis error increased rapidly when the peak width went beneath this limit (a width approximating half the reactor’s fuel pin-to-pin pitch); this resolvability limit is specific to the SPND configuration and fuel segmentation considered herein. The synthesis error increased linearly with perturbation peak amplitude, whereas the convergence time increased nonlinearly. Perturbations located closer to the center of the core were synthesized more accurately, albeit with a higher number of required iterations. These findings provide a qualitative and quantitative understanding of the accuracy and speed at which different types of spatial power perturbations can be resolved in light-water reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Cr 2.37 Ga 3 Se 8 : A Quasi-Two-Dimensional Magnetic Semiconductor

We present a novel magnetic semiconductor, Cr 2.37 Ga 3 Se 8 , synthesized by partially replacing magnetic Cr 3+ in antiferromagnetic Cr 5+δ Se 8 with nonmagnetic Ga 3+ . The crystal structure of Cr 2.37 Ga 3 Se 8 was determined by both powder and single-crystal X-ray diffraction. The title compound crystallizes in a monoclinic structure with space group C2/m (No. 12). In Cr 2.37 Ga 3 Se 8 , the Cr atoms are surrounded by 6 Se atoms and form filled octahedral clusters, while Ga atoms are centered in the Se 4 tetrahedral clusters. The two kinds of clusters pack alternatingly along the c-axis, which results in a quasi-two-dimensional layered structure. The magnetization (M) measurement shows the development of short-range ferromagnetic coupling below the Curie-Weiss (CW) temperature θ CW ~ 92 K, evidenced by the nonlinear field dependence of M. However, the magnetic susceptibility exhibits a peak at low fields at ~18 K, indicating the existence of an antiferromagnetic interaction as well. Electronic structure calculations using the WIEN2k program in the local spin density approximation indicate that the magnetism arises exclusively from local moments of the Cr atoms. Here, the electrical resistivity measurement of the Cr 2.37 Ga 3 Se 8 sample confirms that this material is a semiconductor with the band gap ~0.26 eV. Meanwhile, the experimental band gap (~0.26 eV) is close to the theoretical prediction using the WIEN2k program (~0.35 eV).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

The roles of kinematic constraint and diffusion in non-equilibrium solid state phase transformations of Ti-6Al-4V

A solid state phase transformation of Ti-6Al-4V was studied using high speed in situ x-ray diffraction measurements made during rapid cooling of a cold metal transfer arc weld bead deposited onto a water cooled substrate. Analysis of body centered cubic (BCC) and hexagonal close packed (HCP) lattices revealed an abrupt, nonlinear shift in the lattice parameters of both phases just after the HCP phase had nucleated. Postmortem transmission electron microscopy confirmed that V diffusion was mostly suppressed during cooling. Together, these results indicate that at this cooling rate of approximately 10 4 K/s, which is representative of cooling rates of many additive manufacturing and welding processes, kinematic coherency of the BCC–HCP interfaces gives rise to the anomalous lattice expansion and contraction behaviors of both phases during the initial nucleation and growth stages of (mostly) martensitic transformation from BCC to HCP; the role of diffusion in such lattice anomalies is shown to be minimal.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Stability of new exact solutions of the nonlinear Schrödinger equation in a Pöschl–Teller external potential

Here, we discuss the stability properties of the solutions of the general nonlinear Schrödinger equation in 1 + 1 dimensions in an external potential derivable from a parity-time ($ \newcommand{\PT}{\mathcal{PT}} \PT$ ) symmetric superpotential $W(x)$ that we considered earlier in Kevrekidis et al (2015 Phys. Rev. E 92 042901). In particular we consider the nonlinear partial differential equation $ \{i \, \partial_t + \partial_x^2 - V(x) + g \vert \psi(x, t) \vert ^{2\kappa} \} \, \psi(x, t) = 0 \>, $ for arbitrary nonlinearity parameter κ, where $g= \pm1$ and V is the well known Pöschl–Teller potential which we allow to be repulsive as well as attractive. Using energy landscape methods, linear stability analysis as well as a time dependent variational approximation, we derive consistent analytic results for the domains of instability of these new exact solutions as a function of the strength of the external potential and κ. For the repulsive potential we show that there is a translational instability which can be understood in terms of the energy landscape as a function of a stretching parameter and a translation parameter being a saddle near the exact solution. In this case, numerical simulations show that if we start with the exact solution, the initial wave function breaks into two pieces traveling in opposite directions. If we explore the slightly perturbed solution situations, a 1% change in initial conditions can change significantly the details of how the wave function breaks into two separate pieces. For the attractive potential, changing the initial conditions by 1% modifies the domain of stability only slightly. For the case of the attractive potential and negative g perturbed solutions merely oscillate with the oscillation frequencies predicted by the variational approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Thermal instability of a reconnecting current layer as a trigger for solar flares

The stability of small perturbations of a reconnecting current layer (CL) in a plasma with a strong magnetic field has been investigated in the approximation of dissipative magnetohydrodynamics. The case where the wavevector of the perturbations is parallel to the electric current in the CL has been considered. The suppression of plasma heat conduction by a magnetic field perturbation inside the CL is shown to be responsible for the instability. At the linear stage of instability development, the perturbations grow with the characteristic radiative plasma cooling time calculated in the approximation of an optically thin plasma with cosmic abundances of elements. The formation of a periodic structure of cold and hot magnetic flux tubes, viz., filaments, located across the direction of the electric current, should be expected at the nonlinear stage of the instability in the CL. The proposed mechanism of the thermal CL instability can explain the sequential brightening (ignition) in the arcades of magnetic loops in solar flares.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

A model of the effect of the satellite photosheath on a double floating probe system

Description of a model developed to predict the interactions between multiple electrostatic probes separated by distances that are shorter than or comparable to the characteristic plasma shielding length. The computation involves crude approximations accounting for space charge effects, currents due to the emission and exchange of photoelectrons, and currents from the ambient plasma. Floating potentials are calculated by requiring current balance for each probe and solving by iterating the resulting coupled nonlinear equations. Although only qualitative accuracy is claimed, predictions using the model are found to be in good agreement with observations from the Explorer 45 satellite. The Explorer 45 satellite measures the potential difference between a pair of symmetrical spherical floating probes. The Explorer 45 data show that with increasing altitude, first, the sunward sphere will typically assume a more negative potential than the antisunward sphere; second, the polarities of the probes in the region L = 2 to 4 may switch, sometimes several times; and, finally, near the plasmapause the detector output may saturate, the sunward sphere being more negative.

Cauffman, D. P.↗
NASA NTRS1974

Miniature load cell instrumentation for finite deformation biaxial testing of elastomers

Accuracy of biaxial test equipment to examine the nonlinear mechanical behavior of thin sheet specimens of elastomeric materials has been hampered by lack of precise determination of the force distribution along the sides of the specimen. It has been necessary to use an effective width established by approximate means to obtain the stress from the overall force applied along the edges. To avoid this experimental difficulty, individual miniature proof-ring load cells utilizing semiconductor strain gages have been developed and applied to the support hooks for the thin sheet specimens. Typical results are shown for time-dependent stress distributions for all degrees of biaxiality. An accurate evaluation of the effective specimen width is now possible.

Arenz, R. J.↗
NASA NTRS1975

A General Simulator Using State Estimation for a Space Tug Navigation System

A general simulation program is presented (GSP) involving nonlinear state estimation for space vehicle flight navigation systems. A complete explanation of the iterative guidance mode guidance law, derivation of the dynamics, coordinate frames, and state estimation routines are given so as to fully clarify the assumptions and approximations involved so that simulation results can be placed in their proper perspective. A complete set of computer acronyms and their definitions as well as explanations of the subroutines used in the GSP simulator are included. To facilitate input/output, a complete set of compatable numbers, with units, are included to aid in data development. Format specifications, output data phrase meanings and purposes, and computer card data input are clearly spelled out. A large number of simulation and analytical studies were used to determine the validity of the simulator itself as well as various data runs.

Boland, J. S., III↗
NASA NTRS1977

Inversion of passive microwave remote sensing data from satellites

Global passive microwave observations from earth-orbiting satellites have mapped humidity and liquid water over ocean, temperature profiles, ice and snow, and other geophysical parameters. In most applications, the inversion problem is adequately approximated as linear with jointly Gaussian statistics, and, thus, a linear retrieval performs well. In some cases, the problem is typically factored into a decision process followed by appropriate linear or quasilinear processes. Certain problems, however, require more powerful nonlinear or nonstationary procedures, such as Kalman filtering.

Staelin, D. H.↗
NASA NTRS1985

Reduction method for thermal analysis of complex aerospace structures

A reduction method which combines classical Rayleigh-Ritz modal superposition techniques with contemporary finite-element methods is applied to transient nonlinear thermal analysis of aerospace structures. The essence of the method is the use of a few thermal modes from eigenvalue analyses as basis vectors to represent the temperature response in the structure. The method is used to obtain approximate temperature histories for a portion of the Shuttle orbiter wing subject to reentry heating and for a large space antenna reflector subject to heating associated with a low Earth orbit. The reduction method has excellent potential for significant size reduction for radiation-dominated problems such as the antenna reflector. However, for conduction-dominated problems such as the Shuttle wing, especially those with complex spatial and temporal variations in the applied heating, additional work appears necessary to find alternate sources of basis vectors which will permit significant problem size reductions.

Shore, C. P.↗
NASA NTRS1988

The centring dynamics of a thin liquid shell in capillary oscillations

The physical mechanism governing the centring of a hollow liquid shell in capillary oscillations, which has been observed in experiments, is investigated theoretically. First, the shell is assumed to be inviscid and to have a thickness that is much less than its spherical radius. A system of one-dimensional nonlinear equations of motion is derived using a thin-sheet model. From a numerical study the nonlinear effects of the wave are found to cause the core to oscillate slowly relative to the shell while the centre of mass of the whole system remains stationary. The effects of small viscosity are then considered in an approximation. Finally the strength of the centring mechanism is compared with that of the decentring effect due to buoyancy. The findings are consistent with the limited experimental information available.

Lee, Chun P.↗
NASA NTRS1992

Quantum spatial propagation of squeezed light in a degenerate parametric amplifier

Differential equations which describe the steady state spatial evolution of nonclassical light are established using standard quantum field theoretic techniques. A Schroedinger equation for the state vector of the optical field is derived using the quantum analog of the slowly varying envelope approximation (SVEA). The steady state solutions are those that satisfy the time independent Schroedinger equation. The resulting eigenvalue problem then leads to the spatial propagation equations. For the degenerate parametric amplifier this method shows that the squeezing parameter obey nonlinear differential equations coupled by the amplifier gain and phase mismatch. The solution to these differential equations is equivalent to one obtained from the classical three wave mixing steady state solution to the parametric amplifier with a nondepleted pump.

Deutsch, Ivan H.↗