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,639 records · Page 91

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Detailed characterization of runaway electron driven whistler waves in low-density DIII-D discharges

RE-driven whistler waves during quiescent DIII-D shots have been investigated further. The waves are confirmed to be mostly perpendicularly propagating and are observed for the first time with frequencies up to 700 MHz. Phase-spectral analysis has been used to infer their toroidal mode numbers, n, which are expected to scale with the wavenumber, k, of the mode. Though we derive a theoretical scaling of k ≈ 4n, the measured mode numbers are found to exhibit a very weak dependence on k. In addition, increases in synchrotron emission have been found to consistently lag whistler wave bursts by roughly 3–5 ms, suggesting the waves are causing pitch-angle scattering, since the emitted synchrotron radiation is a strong function of the REs' perpendicular energy. The stronger the wave bursts, the greater the subsequent increase in synchrotron emission. A predator-prey model is used to describe these nonlinear wave-particle interactions, from which the wave damping rates and the loss parameter can be inferred. The damping rates are found to be of the order of (1.6 ± 0.8) × 10 4 /s, and the unitless loss parameter is found to be approximately 2, suggesting that the loss mechanism is diffusive. These observations will serve to validate models of RE-driven waves in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Limitations for Quantum Algorithms to Solve Turbulent and Chaotic Systems

We investigate the limitations of quantum computers for solving nonlinear dynamical systems. In particular, we tighten the worst-case bounds of the quantum Carleman linearisation (QCL) algorithm answering one of their open questions. We provide a further significant limitation for any quantum algorithm that aims to output a quantum state that approximates the normalized solution vector. Given a natural choice of coordinates for a dynamical system with one or more positive Lyapunov exponents and solutions that grow sub-exponentially, we prove that any such algorithm has complexity scaling at least exponentially in the integration time. As such, an efficient quantum algorithm for simulating chaotic systems or regimes is likely not possible.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1972

Application of the Galerkin Method in the Design of Stable Liquid Rocket Motors

A brief description of results obtained in the continuing investigations is presented. The following areas are covered: (1) the development of the third-order, multimode theory to study the behavior of large amplitude transverse instabilities, (2) nonlinear axial mode instability, and (3) the influence of the functional form of the unsteady combustion response function upon the stability characteristics of rocket motors. Also a study was begun to improve the second order potential theory in order to provide a better approximation to the nozzle boundary condition.

Zinn, B. T.↗

NASA NTRS1975

Research study on stabilization and control: Modern sampled-data control theory. Design of the large space telescope system

Conditions of self-sustained oscillations in a two-axis model of the nonlinear LST system are studied. The describing function of the CMG frictional nonlinearity of the LST system is used for the analysis, as well as continuous-data and discrete-data models of the simplified LST control system. A numerical-iterative method is described for the analysis of the two-axis system. Approximation methods and the direct plotting of the stability equation are implemented in the study. It is shown that although the dynamics of the two axes are identical, the amplitudes of self-sustained oscillations in the two axes may in principle be different. Analysis shows that the LST systems are of equal amplitudes but with 180-degree phase shift.

Kuo, B. C.↗

NASA NTRS1975

Performance analysis of wideband data and television channels

Several aspects are discussed of space shuttle communications, including the return link (shuttle-to-ground) relayed through a satellite repeater (TDRS). The repeater exhibits nonlinear amplification and an amplitude-dependent phase shift. Models were developed for various link configurations, and computer simulation programs based on these models are described. Certain analytical results on system performance were also obtained. For the system parameters assumed, the results indicate approximately 1 db degradation relative to a link employing a linear repeater. While this degradation is dependent upon the repeater, filter bandwidths, and modulation parameters used, the programs can accommodate changes to any of these quantities. Thus the programs can be applied to determine the performance with any given set of parameters, or used as an aid in link design.

Geist, J. M.↗

NASA NTRS1977

Variants and extensions of a fast direct numerical cauchy-riemann solver, with illustrative applications

Revised and extended versions of a fast, direct (noniterative) numerical Cauchy-Riemann solver are presented for solving finite difference approximations of first order systems of partial differential equations. Although the difference operators treated are linear and elliptic, one significant application of these extended direct Cauchy-Riemann solvers is in the fast, semidirect (iterative) solution of fluid dynamic problems governed by the nonlinear mixed elliptic-hyperbolic equations of transonic flow. Different versions of the algorithms are derived and the corresponding FORTRAN computer programs for a simple example problem are described and listed. The algorithms are demonstrated to be efficient and accurate.

Martin, E. D.↗

NASA NTRS1979

A statistical model for the initial stellar mass function

The nonlinear stages of the fragmentation of a collapsing molecular cloud are modeled by coagulation theory. Several distinct physical processes are discussed, including protostellar fragment coalescence, gas accretion, and binary formation. This work confirms and extends the earlier result of Nakano (1966) that an approximately self-similar limiting form of the mass spectrum develops after several mean initial collision times. An approximate solution to the velocity-averaged coagulation equation is given for an arbitrary power-law dependence of the coalescence rate on mass, with dimensional dependence proportional to m to the power lambda; i.e., the asymptotic mass spectrum varies as m to the -3 lambda/2 power at small masses and cuts off exponentially at large masses, the characteristic mass depending on the number of collision times elapsed. Simple physical arguments suggest that lambda may increase with increasing mass, but is restricted to the range from 2/3 to 4/3. A large fraction of collisions could result in binary formation.

Silk, J.↗

NASA NTRS1981

Bifurcation and post-buckling analysis of laminated composite plates via reduced basis technique

A reduced basis technique and a problem-adaptive computational algorithm are presented for the bifurcation and post-buckling analysis of laminated anisotropic plates. The computational algorithm can be conveniently divided into three distinct stages. The first stage is that of determining the bifurcation point. The plate is discretized by using displacement finite element (or finite difference) models. The special symmetries exhibited by the response of the anisotropic plate are used to reduce the size of the analysis region. The vector of unknown nodal parameters is expressed as a linear combination of a small number of basis vectors, and a Rayleigh-Ritz technique is used to approximate the finite element equations by a small system of algebraic equations. The reduced equations are used to determine the bifurcation point and the associated eigenmode of the panel. In the second stage of the bifurcation buckling mode is used to obtain a nonlinear solution in the vicinity of the bifurcation point and new (updated) sets of basis vectors and reduced equations are generated. In the third stage the reduced equations are used to trace the post-buckling paths. The effectiveness of the proposed technique for predicting the bifurcation and post-buckling behavior of plates is demonstrated by means of numerical examples for plates loaded by means of prescribed edge displacements.

Noor, A. K.↗

NASA NTRS1986

Anisotropic constitutive model for nickel base single crystal alloys: Development and finite element implementation

A tool for the mechanical analysis of nickel base single crystal superalloys, specifically Rene N4, used in gas turbine engine components is developed. This is achieved by a rate dependent anisotropic constitutive model implemented in a nonlinear three dimensional finite element code. The constitutive model is developed from metallurigical concepts utilizing a crystallographic approach. A non Schmid's law formulation is used to model the tension/compression asymmetry and orientation dependence in octahedral slip. Schmid's law is a good approximation to the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response, and strain rate sensitivity of these alloys. Methods for deriving the material constants from standard tests are presented. The finite element implementation utilizes an initial strain method and twenty noded isoparametric solid elements. The ability to model piecewise linear load histories is included in the finite element code. The constitutive equations are accurately and economically integrated using a second order Adams-Moulton predictor-corrector method with a dynamic time incrementing procedure. Computed results from the finite element code are compared with experimental data for tensile, creep and cyclic tests at 760 deg C. The strain rate sensitivity and stress relaxation capabilities of the model are evaluated.

Dame, L. T.↗

NASA NTRS1987

Model-size reduction for the analysis of symmetric structures with asymmetric boundary conditions

A simple computational procedure is presented for reducing the size of the analysis model for a symmetric structure with asymmetric boundary conditions to that of the corresponding structure with symmetric boundary conditions. The procedure is based on approximating the asymmetric response of the structure by a linear combination of symmetric and antisymmetric global approximation vectors (or modes). The key elements of the procedure are (1) restructuring the governing finite-element equations to delineate the contributions to the symmetric and antisymmetric components of the asymmetric response, (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The effectiveness of the computational procedure is demonstrated by means of numerical examples of linear static problems of shells, and its potential for solving nonlinear problems is discussed.

Noor, Ahmed K.↗

NASA NTRS1998

Nonlinear Analysis of the Space Shuttle Superlightweight LO2 Tank: Behavior Under 3g End-of-Flight Loads - Part 2

Results of linear bifurcation and nonlinear analyses of the Space Shuttle super lightweight (SLWT) external liquid-oxygen (LO2) tank are presented for an important end-of-flight loading condition. These results illustrate an important type of response mode for thin-walled shells, that are subjected to combined mechanical and thermal loads, that may be encountered in the design of other liquid-fuel launch vehicles. Linear bifurcation analyses are presented that predict several nearly equal eigenvalues that correspond to local buckling modes in the aft dome of the LO2 tank. In contrast, the nonlinear response phenomenon is shown to consist of a short-wavelength bending deformation in the aft elliptical dome of the LO2 tank that grows in amplitude in a stable manner with increasing load. Imperfection sensitivity analyses are presented that show that the presence of several nearly equal eigenvalues does not lead to a premature general instability mode for the aft dome. For the linear bifurcation and nonlinear analyses, the results show that accurate predictions of the response of the shell generally require a large-scale, high fidelity finite-element model. Results are also presented that show that the SLWT LO2 tank can support loads in excess of approximately 1.9 times the values of the operational loads considered.

Nemeth, Michael P.↗

NASA NTRS2014

Fuzzy Model-based Pitch Stabilization and Wing Vibration Suppression of Flexible Wing Aircraft.

This paper presents a fuzzy nonlinear controller to regulate the longitudinal dynamics of an aircraft and suppress the bending and torsional vibrations of its flexible wings. The fuzzy controller utilizes full-state feedback with input constraint. First, the Takagi-Sugeno fuzzy linear model is developed which approximates the coupled aeroelastic aircraft model. Then, based on the fuzzy linear model, a fuzzy controller is developed to utilize a full-state feedback and stabilize the system while it satisfies the control input constraint. Linear matrix inequality (LMI) techniques are employed to solve the fuzzy control problem. Finally, the performance of the proposed controller is demonstrated on the NASA Generic Transport Model (GTM).

Stability and Control↗

NASA NTRS2013

Characterization of Swirl-Venturi Lean Direct Injection Designs for Aviation Gas-Turbine Combustion

Injector geometry, physical mixing, chemical processes, and engine cycle conditions together govern performance, operability and emission characteristics of aviation gas-turbine combustion systems. The present investigation explores swirl-venturi lean direct injection combustor fundamentals, characterizing the influence of key geometric injector parameters on reacting flow physics and emission production trends. In this computational study, a design space exploration was performed using a parameterized swirl-venturi lean direct injector model. From the parametric geometry, 20 three-element lean direct injection combustor sectors were produced and simulated using steady-state, Reynolds-averaged Navier-Stokes reacting computations. Species concentrations were solved directly using a reduced 18-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear κ-ε model. Results demonstrate sensitivities of the geometric perturbations on axially averaged flow field responses. Output variables include axial velocity, turbulent kinetic energy, static temperature, fuel patternation and minor species mass fractions. Significant trends have been reduced to surrogate model approximations, intended to guide future injector design trade studies and advance aviation gas-turbine combustion research.

Combustion Chamber↗

DOE OSTI2022

Determination of the nonlinear thermo-optic coefficient of silicon nitride and oxide using an effective index method

There is little literature characterizing the temperature-dependent thermo-optic coefficient (TOC) for low pressure chemical vapor deposition (LPCVD) silicon nitride or plasma enhanced chemical vapor deposition (PECVD) silicon dioxide at temperatures above 300 K. In this study, we characterize these material TOC’s from approximately 300-460 K, yielding values of (2.51 ± 0.08) · 10 −5 K −1 for silicon nitride and (5.67 ± 0.53) · 10 −6 K −1 for silicon oxide at room temperature (300 K). We use a simplified experimental setup and apply an analytical technique to account for thermal expansion during the extraction process. We also show that the waveguide geometry and method used to determine the resonant wavelength have a substantial impact on the precision of our results, a fact which can be used to improve the precision of numerous ring resonator index sensing experiments.

36 MATERIALS SCIENCE↗

NASA NTRS1975

Multidimensional Hermite-Gaussian quadrature formulae and their application to nonlinear estimation

A simplified technique is proposed for calculating multidimensional Hermite-Gaussian quadratures that involves taking the square root of a matrix by the Cholesky algorithm rather than computation of the eigenvectors of the matrix. Ways of reducing the dimension, number, and order of the quadratures are set forth. If the function f(x) under the integral sign is not well approximated by a low-order algebraic expression, the order of the quadrature may be reduced by factoring f(x) into an expression that is nearly algebraic and one that is Gaussian.

Mcreynolds, S. R.↗

DOE OSTI2024

VpROM: a novel variational autoencoder-boosted reduced order model for the treatment of parametric dependencies in nonlinear systems

Reduced Order Models (ROMs) are of considerable importance in many areas of engineering in which computational time presents difficulties. Established approaches employ projection-based reduction, such as Proper Orthogonal Decomposition. The limitation of the linear nature of such operators is typically tackled via a library of local reduction subspaces, which requires the assembly of numerous local ROMs to address parametric dependencies. Our work attempts to define a more generalisable mapping between parametric inputs and reduced bases for the purpose of generative modeling. We propose the use of Variational Autoencoders (VAEs) in place of the typically utilised clustering or interpolation operations, for inferring the fundamental vectors, termed as modes, which approximate the manifold of the model response for any and each parametric input state. The derived ROM still relies on projection bases, built on the basis of full-order model simulations, thus retaining the imprinted physical connotation. However, it additionally exploits a matrix of coefficients that relates each local sample response and dynamics to the global phenomena across the parametric input domain. The VAE scheme is utilised for approximating these coefficients for any input state. This coupling leads to a high-precision low-order representation, which is particularly suited for problems where model dependencies or excitation traits cause the dynamic behavior to span multiple response regimes. Moreover, the probabilistic treatment of the VAE representation allows for uncertainty quantification on the reduction bases, which may then be propagated to the ROM response. The performance of the proposed approach is validated on an open-source simulation benchmark featuring hysteresis and multi-parametric dependencies, and on a large-scale wind turbine tower characterised by nonlinear material behavior and model uncertainty.

Conditional VAEs↗