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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations

In this work, we introduce a new residual distribution (RD) framework for the design of bound-preserving high-resolution finite element schemes. The continuous and discontinuous Galerkin discretizations of the linear advection equation are modified to construct local extremum diminishing (LED) approximations. To that end, we perform mass lumping and redistribute the element residuals in a manner which guarantees the LED property. The hierarchical correction procedure for high-order Bernstein finite element discretizations involves localization to subcells and definition of bound-preserving weights for subcell contributions. Using strong stability preserving (SSP) Runge–Kutta methods for time integration, we prove the validity of discrete maximum principles under CFL-like time step restrictions. The low-order version of our method has roughly the same accuracy as the one derived from a piecewise (multi)-linear approximation on a submesh with the same nodal points. In high-order extensions, we use an element-based flux-corrected transport (FCT) algorithm which can be interpreted as a nonlinear RD scheme. The proposed LED corrections are tailor-made for matrix-free implementations which avoid the rapidly growing cost of matrix assembly for high-order Bernstein elements. The results for 1D, 2D, and 3D test problems compare favorably to those obtained with the best matrix-based approaches.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Correlator convolutional neural networks as an interpretable architecture for image-like quantum matter data

Image-like data from quantum systems promises to offer greater insight into the physics of correlated quantum matter. However, the traditional framework of condensed matter physics lacks principled approaches for analyzing such data. Machine learning models are a powerful theoretical tool for analyzing image-like data including many-body snapshots from quantum simulators. Recently, they have successfully distinguished between simulated snapshots that are indistinguishable from one and two point correlation functions. Thus far, the complexity of these models has inhibited new physical insights from such approaches. Here, we develop a set of nonlinearities for use in a neural network architecture that discovers features in the data which are directly interpretable in terms of physical observables. Applied to simulated snapshots produced by two candidate theories approximating the doped Fermi-Hubbard model, we uncover that the key distinguishing features are fourth-order spin-charge correlators. Our approach lends itself well to the construction of simple, versatile, end-to-end interpretable architectures, thus paving the way for new physical insights from machine learning studies of experimental and numerical data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Activity-induced propulsion and separation of passive chiral particles in liquids

Recent literature shows that chiral particles driven by vorticity gradients of a base Newtonian liquid can separate only in nonlinear or nonstationary flows. Here, we show that a base Newtonian liquid endowed with rotational degrees of freedom, can give rise to chiral particle propulsion and separation even in the Stokes flow approximation. We consider active mechanisms generating spatially inhomogeneous vorticity in the base liquid in a channel such as a magnetic torque, a pressure gradient and a sliding wall. These mechanisms generate regions of spatially inhomogeneous vorticity maintained by the internal angular momentum, propelling the chiral particles according to their handedness. Similar conclusions can be reached when one of the channel walls is replaced by a free surface acted on by surface shear. This behavior is in contrast to the case of standard Newtonian liquids where chiral particle propulsion generated by vorticity gradients is not possible in Stokes flow. Furthermore, this effect can find applications in industry, biology and medicine.

42 ENGINEERING↗

DOE OSTI2017

Quantifying Fusion Born Ion Populations in Magnetically Confined Plasmas using Ion Cyclotron Emission

Ion cyclotron emission (ICE) offers unique promise as a diagnostic of the fusion born alpha-particle population in magnetically confined plasmas. Pioneering observations from JET and TFTR found that ICE intensity P ICE scales approximately linearly with the measured neutron flux from fusion reactions, and with the inferred concentration, n /n i , of fusion-born alpha-particles confined within the plasma. We present fully nonlinear self-consistent kinetic simulations that reproduce this scaling for the first time. This resolves a longstanding question in the physics of fusion alpha particle confinement and stability in MCF plasmas. It confirms the MCI as the likely emission mechanism and greatly strengthens the basis for diagnostic exploitation of ICE in future burning plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Connecting relativistic density functional theory to microscopic calculations

The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Linear Distribution System State Estimation Using Synchrophasor Data and Pseudo-Measurement

State estimation is often a challenging task in distribution systems due to deploying a limited number of measurement devices. Moreover, the integration of precise distributionlevel phasor measurement unites, a.k.a micro-PMUs, along with inaccurate pseudo-measurements in state estimation introduces another challenge. These issues might decrease the efficiency of traditional standard weighted least squares (WLS) for distribution system state estimation. This paper proposes a novel linear distribution system state estimation with taken different types of measurements, including micro-PMU measurements and line current measurements, into account. To involve pseudomeasurements into the linear state estimation, a linearization method based on the Taylor’s approximation is adopted to reformulate pseudo-measurement functions in a linear form. The results obtained via numerical simulations show that the proposed linear state estimation method performs similarly to the standard nonlinear WLS estimator. In addition, the sensitivity analyses results show that our method has a better performance compared to WLS once a limited number of highly precise micro-PMUs are accompanied with inaccurate pseudo-measurements in state estimation.

State estimation Distribution system Micro-PMU Lin↗

DOE OSTI2018

Global Solution Strategies for the Network-Constrained Unit Commitment Problem with AC Transmission Constraints

Here, we propose a novel global solution algorithm for the network-constrained unit commitment problem that incorporates a nonlinear alternating current (AC) model of the transmission network, which is a nonconvex mixed-integer nonlinear programming (MINLP) problem. Our algorithm is based on the multi-tree global optimization methodology, which iterates between a mixed-integer lower-bounding problem and a nonlinear upper-bounding problem. We exploit the mathematical structure of the unit commitment problem with AC power flow constraints (UC-AC) and leverage second-order cone relaxations, piecewise outer approximations, and optimization-based bounds tightening to guarantee a globally optimal solution at convergence. Numerical results on four benchmark problems illustrate the effectiveness of our algorithm, both in terms of convergence rate and solution quality.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1975

Investigation B: Laminar confined coaxial entrance flow with heat generation

The results of a parametric study on the entrance flow region in a gas core nuclear reactor are presented. The physical system is modeled as laminar confined, coaxial flow with heat generation in the inner fluid. The governing equations include the boundary layer approximations and the assumptions of only radial radiative transport of energy represented as an energy diffusion term. The Von Mises transformation and a zeta transformation are used to transform the equations into nonlinear nonhomogeneous convective-diffusion equations. A unique combination of forward and backward difference equations which yields accurate results at moderate computational times, is used in the numerical method. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Bobba, G. K. M.↗

NASA NTRS1976

On the response of superpressure balloons to displacements from equilibrium density level

The response of a superpressure balloon to an initial displacement from its constant-density floating level is examined. An approximate solution is obtained to the governing vertical equation of motion for constant-density superpressure balloons. This solution is used to filter out neutrally buoyant oscillations in balloon records despite the nonlinear behavior of the balloon. The graph depicting the pressure data after deconvolution between the raw pressure data and the normalized balloon wavelet shows clearly the strong filtering-out of the neutral buoyancy oscillations.

Levanon, N.↗

NASA NTRS1976

Effects of stream-associated fluctuations upon the radial variation of average solar-wind parameters

The effects of nonlinear fluctuations due to solar wind streams upon radial gradients of average solar wind parameters are computed, using a numerical MHD model for both spherically symmetric time dependent and corotating equatorial flow approximations. Significant effects of correlations are found between fluctuations upon the gradients of azimuthal magnetic fields, radial velocity, density and azimuthal velocity. Between 400 to 900 solar radii stream interactions have transferred the major portion of the angular momentum flux to the magnetic field; at even greater distances the plasma again carries the bulk of the angular momentum flux. The average azimuthal component of the magnetic field may decrease as much as 10% faster than the Archimedean spiral out to 6 AU due to stream interactions, but this result is dependent upon inner boundary conditions.

Goldstein, B. E.↗

NASA NTRS1980

The free-electron laser - Maxwell's equations driven by single-particle currents

It is shown that if single particle currents are coupled to Maxwell's equations, the resulting set of self-consistent nonlinear equations describes the evolution of the electron beam and the amplitude and phase of the free-electron-laser field. The formulation is based on the slowly varying amplitude and phase approximation, and the distinction between microscopic and macroscopic scales, which distinguishes the microscopic bunching from the macroscopic pulse propagation. The capabilities of this new theoretical approach become apparent when its predictions for the ultrashort pulse free-electron laser are compared to experimental data; the optical pulse evolution, determined simply and accurately, agrees well with observations.

Colson, W. B.↗

NASA NTRS1987

Acoustic gravity waves: A computational approach

This paper discusses numerical solutions of a hyperbolic initial boundary value problem that arises from acoustic wave propagation in the atmosphere. Field equations are derived from the atmospheric fluid flow governed by the Euler equations. The resulting original problem is nonlinear. A first order linearized version of the problem is used for computational purposes. The main difficulty in the problem as with any open boundary problem is in obtaining stable boundary conditions. Approximate boundary conditions are derived and shown to be stable. Numerical results are presented to verify the effectiveness of these boundary conditions.

Hariharan, S. I.↗

NASA NTRS1986

Electrostatic instabilities of velocity-space-shell distributions in magnetized plasmas

Electron instabilities of magnetized spherical shell distributions in velocity space with a colder Maxwellian background are investigated analytically with simulations using electrostatic particle codes. The resonant and nonresonant instabilities observed in the particle simulations are in agreement with zeros of the dielectric function, as found from the resonant approximations for waves with an electric field component along the magnetic field or by computation with a root solver code in the case of perpendicular propagation. Saturation of the instabilities is by nonlinear cyclotron resonance with the cold background in the resonant case or by nonstochastic cyclotron harmonic damping by the cold background in the nonresonant case. Instabilities invariably lead to perpendicular acceleration and heating of the cold background to velocities sometimes exceeding the shell velocity.

Sentman, D. D.↗

NASA NTRS1988

Acoustic gravity waves - A computational approach

Hariharan, S. I.↗

NASA NTRS1989

Model and simulation development for turbopump health monitoring by parameter estimation

The nonlinear equations describing the flow, pressure and temperature dynamics of the High Pressure Fuel Turbopump of the Space Shuttle Main Engine are summarized and implemented in a digital simulation. Some closed form approximations to the operating curves of the turbopump are made so that closed form derivatives exist. These equations and approximations are then used to develop an extended Kalman filter that estimates the dynamic states of the turbopump and a parameter that can reflect turbopump degradation. The filter is implemented as a digital computer program and applied to data generated by the digital simulation of the turbopump with substantial noise added to the output signals. The results indicate that the parameter estimate converges quickly and accurately. This suggests that the parameter estimate can be used to monitor for turbopump degradations affecting its value.

Walker, Bruce K.↗

NASA NTRS1990

Krylov subspace methods - Theory, algorithms, and applications

Projection methods based on Krylov subspaces for solving various types of scientific problems are reviewed. The main idea of this class of methods when applied to a linear system Ax = b, is to generate in some manner an approximate solution to the original problem from the so-called Krylov subspace span. Thus, the original problem of size N is approximated by one of dimension m, typically much smaller than N. Krylov subspace methods have been very successful in solving linear systems and eigenvalue problems and are now becoming popular for solving nonlinear equations. The main ideas in Krylov subspace methods are shown and their use in solving linear systems, eigenvalue problems, parabolic partial differential equations, Liapunov matrix equations, and nonlinear system of equations are discussed.

Sad, Youcef↗

NASA NTRS2010

Improving Computational Efficiency of Prediction in Model-Based Prognostics Using the Unscented Transform

Model-based prognostics captures system knowledge in the form of physics-based models of components, and how they fail, in order to obtain accurate predictions of end of life (EOL). EOL is predicted based on the estimated current state distribution of a component and expected profiles of future usage. In general, this requires simulations of the component using the underlying models. In this paper, we develop a simulation-based prediction methodology that achieves computational efficiency by performing only the minimal number of simulations needed in order to accurately approximate the mean and variance of the complete EOL distribution. This is performed through the use of the unscented transform, which predicts the means and covariances of a distribution passed through a nonlinear transformation. In this case, the EOL simulation acts as that nonlinear transformation. In this paper, we review the unscented transform, and describe how this concept is applied to efficient EOL prediction. As a case study, we develop a physics-based model of a solenoid valve, and perform simulation experiments to demonstrate improved computational efficiency without sacrificing prediction accuracy.

Daigle, Matthew John↗

NASA NTRS2011

Uncertainty Analysis via Failure Domain Characterization: Unrestricted Requirement Functions

This paper proposes an uncertainty analysis framework based on the characterization of the uncertain parameter space. This characterization enables the identification of worst-case uncertainty combinations and the approximation of the failure and safe domains with a high level of accuracy. Because these approximations are comprised of subsets of readily computable probability, they enable the calculation of arbitrarily tight upper and lower bounds to the failure probability. The methods developed herein, which are based on nonlinear constrained optimization, are applicable to requirement functions whose functional dependency on the uncertainty is arbitrary and whose explicit form may even be unknown. Some of the most prominent features of the methodology are the substantial desensitization of the calculations from the assumed uncertainty model (i.e., the probability distribution describing the uncertainty) as well as the accommodation for changes in such a model with a practically insignificant amount of computational effort.

Crespo, Luis G.↗