Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear iteration”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Newton's method: A link between continuous and discrete solutions of nonlinear problems

Newton's method for nonlinear mechanics problems replaces the governing nonlinear equations by an iterative sequence of linear equations. When the linear equations are linear differential equations, the equations are usually solved by numerical methods. The iterative sequence in Newton's method can exhibit poor convergence properties when the nonlinear problem has multiple solutions for a fixed set of parameters, unless the iterative sequences are aimed at solving for each solution separately. The theory of the linear differential operators is often a better guide for solution strategies in applying Newton's method than the theory of linear algebra associated with the numerical analogs of the differential operators. In fact, the theory for the differential operators can suggest the choice of numerical linear operators. In this paper the method of variation of parameters from the theory of linear ordinary differential equations is examined in detail in the context of Newton's method to demonstrate how it might be used as a guide for numerical solutions.

Thurston, G. A.↗
DOE OSTI2022

Time-discretization of a plasma-neutral MHD model with a semi-implicit leapfrog algorithm

The semi-implicit leapfrog time-discretization is a workhorse algorithm for initial-value MHD codes to bridge between vastly separated time scales. Inclusion of atomic interactions with neutrals breaks the functional structure of the MHD equations that exploited by the leapfrog. In this work, we address how to best integrate atomic physics into the semi-implicit leapfrog. Following the Crank-Nicolson method, one approach is to time-center the atomic interactions in the linear solver and use a Newton method to include the nonlinear contributions. Alternatively, another family of methods are based on operator-splitting the terms associated with the atomic interactions using a Strang-splitting technique. These methods naturally break equations into constituent ODE and PDE parts and preserve the structure exploited by the semi-implicit leapfrog. We study the accuracy and efficiency of these methods through a battery of 0D and 1D cases and show that a second-order-in-time Douglas-Rachford inspired coupling between the ODE and PDE advances is effective in reducing the time-discretization error to be comparable to that of Crank-Nicolson with Newton iteration of the nonlinear terms. Splitting ODE and PDE parts results in independent matrix solves for each field which reduces the computational cost considerably and provides parallelization over species relative to Crank-Nicolson.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Analysis of pressure-broadened ozone spectra in the 3 micron region

The Molecular Spectroscopy Lab at NASA-Langley has been involved in a long term effort to carefully characterize the infrared spectra of small molecules of atmospheric interest, including methane, water vapor, ozone, and their isotopic counterparts. High resolution gas phase infrared spectra are obtained using both a tunable diode laser system, and the McMath Fourier transform spectrometer at the Kitt Peak Solar Observatory. Spectra are obtained at various pressures and temperatures for pure gas samples, and for samples containing mixtures of the species of interest in nitrogen, oxygen, or air. From these spectra, using a nonlinear least squares fitting technique, spectral parameters of position, intensity, and half-width were determined for varying laboratory conditions that approximate atmospheric conditions experienced in remote sensing situations. These parameters are of interest in theoretical studies of these species, as well as in allowing more accurate interpretation of remote sensing data. The current work involves the analysis of a series of McMath FTIR spectra of ozone broadened by mixing with air, nitrogen, or oxygen. Each spectrum covers the region from 2396 to 4057/cm. Each vibrational band is analyzed by first diving its region into small intervals containing a few well isolated absorption lines of reasonable intensity. Each of these small intervals is fit by multiple iterations of the nonlinear least square program until residuals are minimized to a reasonable value which corresponds to the noise level of the measured spectrum. Intervals for the 3 nu(sub 3) ozone band in the region from 3000 to 3060 wavenumbers are being examined.

Prochaska, Eleanor S.↗
DOE OSTI2021

Limited-angle computed tomography with deep image and physics priors

Computed tomography is a well-established x-ray imaging technique to reconstruct the three-dimensional structure of objects. It has been used extensively in a variety of fields, from diagnostic imaging to materials and biological sciences. One major challenge in some applications, such as in electron or x-ray tomography systems, is that the projections cannot be gathered over all the angles due to the sample holder setup or shape of the sample. This results in an ill-posed problem called the limited angle reconstruction problem. Typical image reconstruction in this setup leads to distortion and artifacts, thereby hindering a quantitative evaluation of the results. To address this challenge, we use a generative model to effectively constrain the solution of a physics-based approach. Our approach is self-training that can iteratively learn the nonlinear mapping from partial projections to the scanned object. Because our approach combines the data likelihood and image prior terms into a single deep network, it is computationally tractable and improves performance through an end-to-end training. We also complement our approach with total-variation regularization to handle high-frequency noise in reconstructions and implement a solver based on alternating direction method of multipliers. We present numerical results for various degrees of missing angle range and noise levels, which demonstrate the effectiveness of the proposed approach.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

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↗
DOE OSTI2021

Testing $f(R)$ gravity with scale dependent cosmic void velocity profiles

Here, we study the impact of cosmological scale modifications to general relativity on the dynamics of halos within voids by comparing N-body simulations incorporating Hu-Sawicki f(R) gravity, with |f R0 | = 10 -6 and 10 -5 , to those of Λ CDM. By examining the radial velocity statistics within voids classified based on their size and density-profile, as “rising” (R-type) or “shell” (S-type), we find that halo motions in small R-type voids, with effective radius <15 Mpc/h, reveal distinctive differences between f(R) and Λ CDM cosmologies. To understand this observed effect, we study the linear and nonlinear fifth forces, and develop an iterative algorithm to accurately solve the nonlinear fifth force equation. We use this to characterize the Chameleon screening mechanism in voids and contrast the behavior with that observed in gravitationally collapsed objects. The force analysis underscores how smaller R-type voids exhibit the highest ratios of fifth force to Newtonian force, which source distinguishable differences in the velocity profiles and thereby provide rich environments in which to constrain gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Impact of K-factor on Short Circuit Program Convergence for Inverter-Based Resources during Faults

Classical short-circuit programs that linearize the power network are no longer applicable for inverter based resources (IBRs), necessitating an iterative approach. Phasor domain programs can model the IBRs using an iterative approach considering nonlinear fault responses. In phasor domain models, the IBR can be modeled as a voltage controlled current source (VCCS) in tabular form with positive, negative, and zero sequence information for balanced and unbalance short-circuit faults. In the VCCS modeling of the IBR, positive and negative incremental reactive current, also known as the k-factor, plays an important role in short circuit program convergence. In this work, a few approaches: conventional VCCS modeling with a k-factor of 2, conventional VCCS modeling with a k-factor of 2 with modified pre-fault voltages, VCCS characteristics based on the power flow solution with a k-factor of 2, and VCCS characteristics based on the power flow solution with a k-factor of 1 are investigated for short circuit program convergence under higher IBR penetration. The IEEE 39 bus New England Test System is taken as the test system, and simulations are carried out in PSS®CAPE 15.0.26 simulation software. Simulation results demonstrate that IBR penetration is higher for the VCCS model, which corresponds to the power flow solution with k-factor 1, compared to other approaches.

42 ENGINEERING↗
NASA NTRS1983

Continuation in a parameter - Experience with viscous and free surface flows

The results of modifications in continuation methods applied to obtain solutions to the Navier-Stokes systems of equations for incompressible, two-dimensional, steady flows are reported. It is shown that parameter continuation permits prediction of accurate, initial estimates for iterative processing of nonlinear finite difference and finite element equations of motions. The new parameter steps are derived from values of the preceding parameter steps. The accuracy of the estimates is ensured with appropriate choices of the step size. The continuation predictor/iterative corrector is demonstrated to trace the branches of parameter space along which steady flow states are found, and techniques are available for tracing multiply branching paths. The techniques are applied to solving the Navier-Stokes equations for flow through a rotating square channel, the formation of a falling liquid curtain, and gyrostatic equilibria of rotating cylindrical drops.

Kheshgi, H. S.↗
NASA NTRS1986

Self-adaptive difference method for the effective solution of computationally complex problems of boundary layer theory

An implicit difference procedure for the solution of equations for a chemically reacting hypersonic boundary layer is described. Difference forms of arbitrary error order in the x and y coordinate plane were used to derive estimates for discretization error. Computational complexity and time were minimized by the use of this difference method and the iteration of the nonlinear boundary layer equations was regulated by discretization error. Velocity and temperature profiles are presented for Mach 20.14 and Mach 18.5; variables are velocity profiles, temperature profiles, mass flow factor, Stanton number, and friction drag coefficient; three figures include numeric data.

Schoenauer, W.↗
NASA NTRS1987

Program VSAERO theory document: A computer program for calculating nonlinear aerodynamic characteristics of arbitrary configurations

The VSAERO low order panel method formulation is described for the calculation of subsonic aerodynamic characteristics of general configurations. The method is based on piecewise constant doublet and source singularities. Two forms of the internal Dirichlet boundary condition are discussed and the source distribution is determined by the external Neumann boundary condition. A number of basic test cases are examined. Calculations are compared with higher order solutions for a number of cases. It is demonstrated that for comparable density of control points where the boundary conditions are satisfied, the low order method gives comparable accuracy to the higher order solutions. It is also shown that problems associated with some earlier low order panel methods, e.g., leakage in internal flows and junctions and also poor trailing edge solutions, do not appear for the present method. Further, the application of the Kutta conditions is extremely simple; no extra equation or trailing edge velocity point is required. The method has very low computing costs and this has made it practical for application to nonlinear problems requiring iterative solutions for wake shape and surface boundary layer effects.

Maskew, Brian↗
NASA NTRS1991

Analysis of pressure-broadened ozone spectra in the 3 micron region

This work involves the analysis of a series of McMath Fourier Transform Infrared (FTIR) spectra of ozone broadened by mixing with air (four different pressures), nitrogen (three pressures), or oxygen (three pressures). Each spectrum covers the region from 2396 to 4057 cm(-1). This study focused on the 3 sub nu sub 3 band in t 3000 to 3060 cm(-1). The band is analyzed by first dividing its region into small intervals containing a few well isolated absorption lines of reasonable intensity. Each of these small intervals is fit by multiple iterations of the nonlinear least squares program until residuals (the difference between calculated and observed spectrum, as a percent of the strongest intensity in the interval) are minimized to a reasonable value which corresponds to the noise level of the measured spectrum. Position, intensity, and half-width are recorded for later analysis. From the measured half-widths, a pressure broadening coefficient was determined for each absorption line. Pressure shifts were determined by comparing observed line positions in the spectra of the diluted ozone samples to tabulated line positions determined from spectra of pure gas samples. Comparisons to other work on ozone indicate that the broadening and shift coefficients determined in this study are consistent with those determined in other spectral regions.

Prochaska, Eleanor S.↗
DOE OSTI2022

Nonlinear Model Predictive Control Based on Real-Time Iteration Scheme for Wave Energy Converters Using WEC-Sim

One of several challenges that wave energy technologies face is their inability to generate electricity cost-competitively with other grid-scale energy generation sources. Several studies have identified two approaches to lower the levelised cost of electricity: reduce the cost over the device's lifetime or increase its overall electrical energy production. Several advanced control strategies have been developed to address the latter. However, only a few take into account the overall efficiency of the power take-off (PTO) system, and none of them solve the optimisation problem that arises at each sampling time on real-time. In this paper, a detailed Nonlinear model predictive control (NMPC) approach based on the real-time iteration (RTI) scheme is presented, and the controller performance is evaluated using a time-domain hydrodynamics model (WEC-Sim). The proposed control law incorporates the PTO system's efficiency in a control law to maximise the energy extracted. The study also revealed that RTI-NMPC clearly outperforms a simple resistive controller.

model predictive control↗
DOE OSTI2022

Nonlinear Model Predictive Control Based on Real-Time Iteration Scheme for Wave Energy Converters Using WEC-Sim: Preprint

One of several challenges that wave energy technologies face is their inability to generate electricity cost-competitively with other grid-scale energy generation sources. Several studies have identified two approaches to lower the levelised cost of electricity: reduce the cost over the device's lifetime or increase its overall electrical energy production. Several advanced control strategies have been developed to address the latter. However, only a few take into account the overall efficiency of the power take-off (PTO) system, and none of them solve the optimisation problem that arises at each sampling time on real-time. In this paper, a detailed Nonlinear model predictive control (NMPC) approach based on the real-time iteration (RTI) scheme is presented, and the controller performance is evaluated using a time-domain hydrodynamics model (WEC-Sim). The proposed control law incorporates the PTO system's efficiency in a control law to maximise the energy extracted. The study also revealed that RTI-NMPC clearly outperforms a simple resistive controller.

model predictive control↗
NASA NTRS1993

Efficient Minimum-Polynomial And Reduced-Rank Extrapolation

MPERRE computer program accelerates convergence of sequence of vectors by use of minimum-polynomial extrapolation (MPE) and reduced-rank extrapolation (RRE). Effective in accelerating convergences of such sequences of vectors as those obtained from iterative solution of systems of linear and nonlinear equations. In conjunction with various iterative techniques, successfully employed in finite-difference solution of large-scale elliptic boundary-value problems and problems in computational-fluid-dynamics. Only input required is sequence of vectors, convergence of which is accelerated. Program economical and easy to use. Written in FORTRAN 77.

Sidi, Avram↗
DOE OSTI2024

Iterative methods in GPU-resident linear solvers for nonlinear constrained optimization

Linear solvers are major computational bottlenecks in a wide range of decision support and optimization computations. The challenges become even more pronounced on heterogeneous hardware, where traditional sparse numerical linear algebra methods are often inefficient. For example, methods for solving ill-conditioned linear systems have relied on conditional branching, which degrades performance on hardware accelerators such as graphical processing units (GPUs). To improve the efficiency of solving ill-conditioned systems, our computational strategy separates computations that are efficient on GPUs from those that need to run on traditional central processing units (CPUs). Our strategy maximizes the reuse of expensive CPU computations. Iterative methods, which thus far have not been broadly used for ill-conditioned linear systems, play an important role in our approach. In particular, we extend ideas from Arioli et al., (2007) to implement iterative refinement using inexact LU factors and flexible generalized minimal residual (FGMRES), with the aim of efficient performance on GPUs. In conclusion, we focus on solutions that are effective within broader application contexts, and discuss how early performance tests could be improved to be more predictive of the performance in a realistic environment.

97 MATHEMATICS AND COMPUTING↗