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,729 records · Page 96

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Aeroelastic modeling of swept tip rotor blades using finite elements

A special finite element for the modeling of a swept tip rotor blade is derived. The swept tip blade undergoes moderate deflections in flap, lag and torsion. The nonlinear, partial differential equations of motion are discretized using a Galerkin finite element method. Tip sweep introduces flap-torsion and lag-axial couplings, and may lead to aeroelastic instabilities associated with frequency coalescence. When frequency coalescence does not occur, sweep is usually stabilizing. A comparison of the approximate model of a swept tip blade, using a blade with straight elastic axis and offsets of aerodynamic centers of gravity, with the exact swept tip model developed in this paper indicates that the approximate model can produce inaccurate results for the case of hingeless rotor blades.

Celi, R.↗

NASA NTRS1989

Kinetic theory model for the flow of a simple gas from a two-dimensional nozzle

A system of nonlinear integral equations equivalent to the Krook kinetic equation for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust two-dimensional nozzles. The method of characteristics was used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it exhausts into a vacuum. Results predict backscatter and show the effect of the inside wall boundary layer on the flowfields external to the nozzle.

Riley, B. R.↗

NASA NTRS1993

Parameter estimation techniques based on optimizing goodness-of-fit statistics for structural reliability

New methods are presented that utilize the optimization of goodness-of-fit statistics in order to estimate Weibull parameters from failure data. It is assumed that the underlying population is characterized by a three-parameter Weibull distribution. Goodness-of-fit tests are based on the empirical distribution function (EDF). The EDF is a step function, calculated using failure data, and represents an approximation of the cumulative distribution function for the underlying population. Statistics (such as the Kolmogorov-Smirnov statistic and the Anderson-Darling statistic) measure the discrepancy between the EDF and the cumulative distribution function (CDF). These statistics are minimized with respect to the three Weibull parameters. Due to nonlinearities encountered in the minimization process, Powell's numerical optimization procedure is applied to obtain the optimum value of the EDF. Numerical examples show the applicability of these new estimation methods. The results are compared to the estimates obtained with Cooper's nonlinear regression algorithm.

Starlinger, Alois↗

DOE OSTI2023

Demonstration of acoustic monitoring for structural health of microreactors: Through use of neural networks and resonant ultrasound spectroscopy

Nuclear microreactors prioritize modularity and portability and are intended to be a cost-effective technology for non-conventional nuclear markets. As such, the development of microreactors into a safe and feasible solution for energy security applications will necessitate the development of non-destructive technologies to monitor the integrity of inaccessible reactors components during operation. This demonstration applies linear and nonlinear acoustic techniques, in combination with machine learning, to detect and classify mechanical changes (stress and damage) in a test article which are broadly representative of potential operating challenges within a functioning microreactor. All necessary data has been collected for this demonstration, with minor experimental issues identified that can be addressed in follow-on work. Motivated by the expected conditions within a functioning microreactor, we have demonstrated our monitoring techniques on a core-block-like test article using an unstructured excitation source that approximates the noisy acoustic environment expected during reactor operation. At all stress states mechanically applied to the test article, a machine learning model using an artificial neural network was able to classify with 100% accuracy whether a 3D laser vibrometry point measurement was made on an intact or artificially defective test article. Further, model predictions about whether the defect interface was rough or smooth were 95% accurate, indicating the ability of acoustic techniques to recover defect characteristics. Resonant ultrasound spectroscopy (RUS) was also applied to the dataset to provide further quantitative insights about material properties. RUS analysis was ultimately hampered by several minor experimental and data issues, limiting results to certain cases for this demonstration. Last, analysis using nonlinear RUS exhibited sensitivity to changing levels of applied stresses for each intact and defective state. As presented in this demonstration, acoustic monitoring exhibits sensitivity to stress changes, which are of concern due to high thermal gradients expected during startup and operation. Further, our techniques distinguish between measurements made on intact and damaged test articles.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1977

The influence of the directional energy distribution on the nonlinear dispersion relation in a random gravity wave field

The influence of the directional distribution of wave energy on the dispersion relation is calculated numerically using various directional wave spectrum models. The results indicate that the dispersion relation varies both as a function of the directional energy distribution and the direction of propagation of the wave component under consideration. Furthermore, both the mean deviation and the random scatter from the linear approximation increase as the energy spreading decreases. Limited observational data are compared with the theoretical results. The agreement is favorable.

Huang, N. E.↗

DOE OSTI2021

B-PINNs: Bayesian Physics-informal Neural Networks for Forward and Inverse PDE Problems with Noisy Data

We propose a Bayesian physics-informed neural network (B-PINN) to solve both forward and inverse nonlinear problems described by partial differential equations (PDEs) and noisy data. In this Bayesian framework, the Bayesian neural network (BNN) combined with a PINN for PDEs serves as the prior while the Hamiltonian Monte Carlo (HMC) or the variational inference (VI) could serve as an estimator of the posterior. B-PINNs make use of both physical laws and scattered noisy measurements to provide predictions and quantify the aleatoric uncertainty arising from the noisy data in the Bayesian framework. Compared with PINNs, in addition to uncertainty quantification, B-PINNs obtain more accurate predictions in scenarios with large noise due to their capability of avoiding overfitting. We conduct a systematic comparison between the two different approaches for the B-PINNs posterior estimation (i.e., HMC or VI), along with dropout used for quantifying uncertainty in deep neural networks. Our experiments show that HMC is more suitable than VI with mean field Gaussian approximation for the B-PINNs posterior estimation, while dropout employed in PINNs can hardly provide accurate predictions with reasonable uncertainty. Finally, we replace the BNN in the prior with a truncated Karhunen-Loève (KL) expansion combined with HMC or a deep normalizing flow (DNF) model as posterior estimators. The KL is as accurate as BNN and much faster but this framework cannot be easily extended to high-dimensional problems unlike the BNN based framework.

Non linear PDEs, Noisy data, Bysian physics inform↗

DOE OSTI2018

An improved iterative HDG approach for partial differential equations

Here we propose and analyze an iterative high-order hybridized discontinuous Galerkin (iHDG) discretization for linear partial differential equations. We improve our previous work [45] in several directions: 1) the improved iHDG approach converges in a finite number of iterations for the scalar transport equation; 2) it is unconditionally convergent for both the linearized shallow water system and the convection–diffusion equation; 3) it has improved stability and convergence rates; 4) we uncover a relationship between the number of iterations and time stepsize, solution order, meshsize and the equation parameters. This allows us to choose the time stepsize such that the number of iterations is approximately independent of the solution order and the meshsize; and 5) we provide both strong and weak scalings of the improved iHDG approach up to 16,384 cores. A connection between iHDG and time integration methods such as parareal and implicit/explicit methods are discussed. Extensive numerical results for linear (and nonlinear) PDEs are presented to verify the theoretical findings.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Sources of error and methods to improve accuracy in interface state density analysis using quasi-static capacitance–voltage measurements in wide bandgap semiconductors

Characterizing interface trap states in commercial wide bandgap devices using frequency-based measurements requires unconventionally high probing frequencies to account for both fast and slow traps associated with wide bandgap materials. The C - ψs technique has been suggested as a viable quasi-static method for determining the interface trap state densities in wide bandgap systems, but the results are shown to be susceptible to errors in the analysis procedure. Here this work explores the primary sources of errors present in the C - ψs technique using an analytical model that describes the apparent response for wide bandgap MOS capacitor devices. Measurement noise is shown to greatly impact the linear fitting routine of the 1/C$^{*2}_s$ vs Ψ s plot to calibrate the additive constant in the surface potential/gate voltage relationship, and an inexact knowledge of the oxide capacitance is also shown to impede interface trap state analysis near the band edge. In addition, a slight nonlinearity that is typically present throughout the 1/C$^{*2}_s$ vs Ψ s plot hinders the accurate estimation of interface trap densities, which is demonstrated for a fabricated n-SiC MOS capacitor device. Methods are suggested to improve quasi-static analysis, including a novel method to determine an approximate integration constant without relying on a linear fitting routine.

42 ENGINEERING↗

DOE OSTI2025

The polarized-signal density matrix: A practical way to recover molecular frame information from isotropic samples

We present a novel approach to model ultrafast time-dependent nonlinear optical polarization sensitive signals emitted from randomly oriented molecules. By projecting the laboratory-frame analyzer polarization axis into the molecular-frame and linking that axis with the density matrix through a tensor product, we demonstrate an approach to find a specific molecular orientation that yields a good approximation to simulated four-wave mixing signals produced by the same model but with averaging over molecular orientation.

Thurston, Richard L↗

DOE OSTI2018

Final technical report

This final report covers the final year of the project. Reports from previous years were included in the group report for that year. The objective of this part of the project was to compare results from two large-scale, nonlinear plasma simulation codes based on magnetohydrodynamics, M3D and M3D-C1, for plasma edge instabilities. This is a difficult undertaking and the comparison was only partially successful because, while the MHD equations are standard, there are a number of valid MHD models of the plasma dissipation that make different approximations to their physics and the numerical implementation of a given model can be quite different depending on the way that the rest of the equations are solved. Large differences arise, and have major effects on the predicted plasma behavior, in the case of the greatest interest to this project, the plasma edge region under conditions where the plasma edge instabilities are strongest and the background plasma changes abruptly from interior hot, dense conditions to an unconfined, cool, diffuse plasma outside the magnetic boundary surface. Experimental observations only partially measure the background plasma and there are many non-MHD processes that can have important effects (studied in the rest of the project). Comparisons for edge instabilities at very small amplitude, the so-called linear modes, had showed good agreement for a few interior plasma instabilities, but strong differences between the default models for the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1980

Investigation of multiple scattering effects in aerosols

The results are presented of investigations on the various aspects of multiple scattering effects on visible and infrared laser beams transversing dense fog oil aerosols contained in a chamber (4' x 4' x 9'). The report briefly describes: (1) the experimental details and measurements; (2) analytical representation of the aerosol size distribution data by two analytical models (the regularized power law distribution and the inverse modified gamma distribution); (3) retrieval of aerosol size distributions from multispectral optical depth measurements by two methods (the two and three parameter fast table search methods and the nonlinear least squares method); (4) modeling of the effects of aerosol microphysical (coagulation and evaporation) and dynamical processes (gravitational settling) on the temporal behavior of aerosol size distribution, and hence on the extinction of four laser beams with wavelengths 0.44, 0.6328, 1.15, and 3.39 micrometers; and (5) the exact and approximate formulations for four methods for computing the effects of multiple scattering on the transmittance of laser beams in dense aerosols, all of which are based on the solution of the radiative transfer equation under the small angle approximation.

Deepak, A.↗

NASA NTRS1991

Kinetic theory model for the flow of a simple gas from a three-dimensional axisymmetric nozzle

A system of nonlinear integral equations equivalent to the Krook kinetic equations for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust three-dimensional axisymmetric nozzles. The method of characteristics is used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it expands into a vacuum. Results predict backscatter and show the effect of the nozzle wall boundary layer on the external flowfields.

Riley, B. R.↗

NASA NTRS1993

Near grazing scattering from non-Gaussian ocean surfaces

We investigate the behavior of the scattered electromagnetic waves from non-Gaussian ocean surfaces at near grazing incidence. Even though the scattering mechanisms at moderate incidence angles are relatively well understood, the same is not true for near grazing rough surface scattering. However, from the experimental ocean scattering data, it has been observed that the backscattering cross section of a horizontally polarized wave can be as large as the vertical counterpart at near grazing incidence. In addition, these returns are highly intermittent in time. There have been some suggestions that these unexpected effects may come from shadowing or feature scattering. Using numerical scattering simulations, it can be shown that the horizontal backscattering cannot be larger than the vertical one for the Gaussian surfaces. Our main objective of this study is to gain a clear understanding of scattering mechanisms underlying the near grazing ocean scattering. In order to evaluate the backscattering cross section from ocean surfaces at near grazing incidence, both the hydrodynamic modeling of ocean surfaces and an accurate near grazing scattering theory are required. For the surface modeling, we generate Gaussian surfaces from the ocean surface power spectrum which is derived using several experimental data. Then, weakly nonlinear large scale ocean surfaces are generated following Longuet-Higgins. In addition, the modulation of small waves by large waves is included using the conservation of wave action. For surface scattering, we use MOM (Method of Moments) to calculate the backscattering from scattering patches with the two scale shadowing approximation. The differences between Gaussian and non-Gaussian surface scattering at near grazing incidence are presented.

Kim, Yunjin↗

NASA NTRS2012

Monitoring Location and Angular Orientation of a Pill

A mobile pill transmitter system moves through, or adjacent to, one or more organs in an animal or human body, while transmitting signals from its present location and/or present angular orientation. The system also provides signals from which the present roll angle of the pill, about a selected axis, can be determined. When the location coordinates angular orientation and the roll angle of the pill are within selected ranges, an aperture on the pill container releases a selected chemical into, or onto, the body. Optionally, the pill, as it moves, provides a sequence of visually perceptible images. The times for image formation may correspond to times at which the pill transmitter system location or image satisfies one of at least four criteria. This invention provides and supplies an algorithm for exact determination of location coordinates and angular orientation coordinates for a mobile pill transmitter (PT), or other similar device that is introduced into, and moves within, a GI tract of a human or animal body. A set of as many as eight nonlinear equations has been developed and applied, relating propagation of a wireless signal between either two, three, or more transmitting antennas located on the PT, to four or more non-coplanar receiving antennas located on a signal receiver appliance worn by the user. The equations are solved exactly, without approximations or iterations, and are applied in several environments: (1) association of a visual image, transmitted by the PT at each of a second sequence of times, with a PT location and PT angular orientation at that time; (2) determination of a position within the body at which a drug or chemical substance or other treatment is to be delivered to a selected portion of the body; (3) monitoring, after delivery, of the effect(s) of administration of the treatment; and (4) determination of one or more positions within the body where provision and examination of a finer-scale image is warranted.

Schipper, John F.↗

NASA NTRS2019

High-Fidelity Multiple-Flyby Trajectory Optimization Using Multiple-Shooting

Rendering a complex spacecraft trajectory in high fidelity can be an expensive endeavor, both computationally and from a human time/cost standpoint. However, in many cases, a low-fidelity trajectory that reasonably approximates a high-fidelity counterpart is much easier to obtain. Thus, it is important to have an efficient process for converting a trajectory from lower-fidelity model to high fidelity. We present a method for converting low-fidelity trajectories into high fidelity that relies on multiple shooting, nonlinear programming, and numerical integration. The procedure converts any zero-radius sphere-of-influence gravity-assist events to fully integrated flyby events. Several numerical examples are presented that showcase the flexibility of the high-fidelity rendering process across multiple mission types and flight regimes.

high-fidelity↗

DOE OSTI2022

Improved CdTe PLR Estimates: Self-Shading and Spectral Mismatch

The RdTools year-on-year method of estimating performance loss rate (PLR) employs a simple normalization to remove the confounding effect of irradiance and temperature variation. However, the normalization's assumption that PV production scales linearly with in-plane broadband irradiance is a worse approximation for CdTe and other technologies with larger spectral sensitivities than it is for the more common c-Si technology. Additionally, CdTe systems using single-axis trackers (-20% of installed US utility-scale capacity) are subject to self-shading in the morning and afternoon, introducing another nonlinearity between PV output and broadband irradiance. Ignoring these effects may subject the estimated PLR to increased uncertainty, and perhaps bias, depending on the character of their short- and long-term variability. In this work we show that including self-shading and spectral mismatch models in the normalization for tracking CdTe systems can result not only in tighter PLR confidence intervals but different median PLRs as well. The shading and spectral models are kept simple to maintain consistency with the RdTools ethos of not requiring detailed system metadata or unusual measurements.

CdTe↗

NASA NTRS2022

Nonlinear Evolution of Instability in a Laminar Separation Bubble at Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

High speed flow↗

NASA NTRS2022

Nonlinear Evolution of Instabilities in a Laminar Separation Bubble at a Hypersonic Mach Number

Hypersonic flow↗