Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximate Newton”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An aggregation-based nonlinear multigrid solver for two-phase flow and transport in porous media

A nonlinear multigrid solver for two-phase flow and transport in a mixed fractional-flow velocity-pressure-saturation formulation is proposed. The solver, which is under the framework of the full approximation scheme (FAS), extends our previous work on nonlinear multigrid for heterogeneous diffusion problems. The coarse spaces in the multigrid hierarchy are constructed by first aggregating degrees of freedom, and then solving some local flow problems. The mixed formulation and the choice of coarse spaces allow us to assemble the coarse problems without visiting finer levels during the solving phase, which is crucial for the scalability of multigrid methods. Specifically, a natural generalization of the upwind flux can be evaluated directly on coarse levels using the precomputed coarse flux basis vectors. We find that the resulting solver is applicable to problems discretized on general unstructured grids. The performance of the proposed nonlinear multigrid solver in comparison with the standard single level Newton's method is demonstrated through challenging numerical examples. It is observed that the proposed solver is robust for highly nonlinear problems and clearly outperforms Newton's method in the case of high Courant-Friedrichs-Lewy (CFL) numbers.

58 GEOSCIENCES↗
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
NASA NTRS1981

Monomethylhydrazine versus hydrazine fuels - Test results using a 100 pound thrust bipropellant rocket engine

A test program was performed to evaluate hydrazine (N2H4) as a fuel for a 445 Newton (100 lbf) thrust bipropellant rocket engine. Results of testing with an identical thruster utilizing monomethylhydrazine (MMH) are included for comparison. Engine performance with hydrazine fuel was essentially identical to that experienced with monomethylhydrazine although higher combustor wall temperatures (approximately 400 F) were obtained with hydrazine. Results are presented which indicate that hydrazine as a fuel is compatible with Marquardt bipropellant rocket engines which use monomethylhydrazine as a baseline fuel.

Smith, J. A.↗
NASA NTRS2010

Probing the Nature of Short Swift Bursts via Deep INTEGRAL Monitoring of GRB 050925

We present results from Swift, XMM-Newton, and deep INTEGRAL monitoring in the region of GRB 050925. This short Swift burst is a candidate for a newly discovered soft gamma-ray repeater (SGR) with the following observational burst properties: 1) galactic plane (b=-0.1 deg) localization, 2) 150 msec duration, and 3) a blackbody rather than a simple power-law spectral shape (with a significance level of 97%). We found two possible X-ray counterparts of GRB 050925 by comparing the X-ray images from Swift XRT and XMM-Newton. Both X-ray sources show the transient behavior with a power-law decay index shallower than -1. We found no hard X-ray emission nor any additional burst from the location of GRB 050925 in approximately 5 Ms of INTEGRAL data. We discuss about the three BATSE short bursts which might be associated with GRB 050925, based on their location and the duration. Assuming GRB 050925 is associated with the H(sub II), regions (W 58) at the galactic longitude of 1=70 deg, we also discuss the source frame properties of GRB 050925.

Sakamoto, T.↗
NASA NTRS2016

Nustar Reveals the Extreme Properties of the Super-Eddington Accreting Supermassive Black Hole in PG 1247+267

PG1247+267 is one of the most luminous known quasars at z approximately 2 and is a strongly super-Eddington accreting supermassive black hole (SMBH) candidate. We obtained NuSTAR data of this intriguing source in December 2014 with the aim of studying its high-energy emission, leveraging the broad band covered by the new NuSTAR and the archival XMM-Newton data. Several measurements are in agreement with the super-Eddington scenario for PG1247+267: the soft power law (gamma = 2.3 +/- 0.1); the weak ionized Fe emission line; and a hint of the presence of outflowing ionized gas surrounding the SMBH. The presence of an extreme reflection component is instead at odds with the high accretion rate proposed for this quasar. This can be explained with three different scenarios; all of them are in good agreement with the existing data, but imply very different conclusions: i) a variable primary power law observed in a low state, superimposed on a reflection component echoing a past, higher flux state; ii) a power law continuum obscured by an ionized, Compton thick, partial covering absorber; and iii) a relativistic disk reflector in a lamp-post geometry, with low coronal height and high BH spin. The first model is able to explain the high reflection component in terms of variability. The second does not require any reflection to reproduce the hard emission, while a rather low high-energy cutoff of approximately 100 keV is detected for the first time in such a high redshift source. The third model require a face-on geometry, which may affect the SMBH mass and Eddington ratio measurements. Deeper X-ray broad-band data are required in order to distinguish between these possibilities.

accretion disks↗
DOE OSTI2020

An investigation of Newton-Sketch and subsampled Newton methods

Sketching, a dimensionality reduction technique, has received much attention in the statistics community. In this paper, we study sketching in the context of Newton's method for solving finite-sum optimization problems in which the number of variables and data points are both large. In this work, we study two forms of sketching that perform dimensionality reduction in data space: Hessian subsampling and randomized Hadamard transformations. Each has its own advantages, and their relative tradeoffs have not been investigated in the optimization literature. Additionally, our study focuses on practical versions of the two methods in which the resulting linear systems of equations are solved approximately, at every iteration, using an iterative solver. The advantages of using the conjugate gradient method vs. a stochastic gradient iteration are revealed through a set of numerical experiments, and a complexity analysis of the Hessian subsampling method is presented.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

A Multi-Epoch Timing and Spectral Study of the ULX NGC 5408 X-1 with XMM-Newton

We report results from extensive new XMM- Newton observations of the ultraluminous X-ray source (ULX) NGC 5408 X-1, one of the few ULXs to show quasi-periodic X-ray variability. We detect quasi-periodic oscillations (QPOs) in each of four new (approximately equal 100 ks each) pointings, expanding the range of frequencies and rms amplitudes observed from the source to 10-40 mHz and 10-45 %, respectively. However, similarly significant variations in the power-law photon spectral index, Gamma, are not observed. We use the results of timing and energy spectral modeling to compare with the timing and spectral correlations seen in stellar-mass systems. We find that the qualitative nature of the timing and energy spectra of NGC 5408 X-1 are very similar to stellar-mass black holes in the steep power-law state exhibiting Type-C QPOs. However, in order for this analogy to quantitatively hold we must only be seeing the so-called saturated portion of the QPO frequency - photon index (or disk flux) relation. Assuming this to be the case, we place a lower limit on the mass of NGC 5408 X-1 of approx greater than 800 Solar Mass. Alternatively, the QPO centroid frequency is largely independent of the spectral parameters, in which case a close analogy of NGC 5408 X-1's mHz QPOs with Type-C QPOs in stellar systems is problematic. Measurement of the source's timing properties over a greater range of spectral parameters (in particular the spectral index) is needed in order to definitively resolve this ambiguity. We searched all the available data for both a broad Fe emission line as well as high frequency QPO analogs (0.1 - 1 Hz), but detected neither. We place upper limits on the equivalent width of any Fe emission feature in the 6 - 7 keY band, and of the amplitude (rms) of a high frequency QPO analog of approx equal 10 eV and approx equal 4%, respectively.

Strohmayer, Tod E.↗
NASA NTRS2017

Identification of the Hard X-Ray Source Dominating the E > 25 keV Emission of the Nearby Galaxy M31

We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 kiloelectronvolts from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-Sky Hard X-Ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 kiloelectronvolts as the counterpart for the first time. In the 0.5-10 kiloelectronvolt band, the source had been classified as an X-ray Binary candidate in various Chandra and XMM-Newton studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E is less than 10 kiloelectronvolts observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT approximately equal to 0.2 kiloelectronvolts) plus a hard spectrum with a power law of gamma approximately equal to 1 and a cutoff around 15-20 kiloelectronvolts, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST (Hubble Space Telescope) images exclude high-mass (greater than 3 times the radius of the moon) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (less than 3 times the radius of the moon M) companion or a symbiotic X-ray binary. We discuss other possibilities in more detail.

galaxies: bulges – galaxies: individual (M31)↗
DOE OSTI2025

Parallel-in-Time Solution of Scalar Nonlinear Conservation Laws

Here, we consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially nonoscillatory (WENO) reconstructions, and in time with high-order explicit Runge–Kutta methods. The solution of the global, discretized space-time problem is sought via a nonlinear iteration that uses a novel linearization strategy in cases of nondifferentiable equations. Under certain choices of discretization and algorithmic parameters, the nonlinear iteration coincides with Newton’s method, although, more generally, it is a preconditioned residual correction scheme. At each nonlinear iteration, the linearized problem takes the form of a certain discretization of a linear conservation law over the space-time domain in question. An approximate parallel-in-time solution of the linearized problem is computed with a single multigrid reduction-in-time (MGRIT) iteration; however, any other effective parallel-in-time method could be used in its place. The MGRIT iteration employs a novel coarse-grid operator that is a modified conservative semi-Lagrangian discretization and generalizes those we have developed previously for nonconservative scalar linear hyperbolic problems. Numerical tests are performed for the inviscid Burgers and Buckley–Leverett equations. For many test problems, the solver converges in just a handful of iterations with a convergence rate independent of mesh resolution, including problems with (interacting) shocks and rarefactions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

An efficient solution technique for shockwave-boundary layer interactions with flow separation and slot suction effects

An efficient method for computing two-dimensional compressible Navier-Stokes flow fields is presented. The solution algorithm is a fully-implicit approximate factorization technique based on an unsymmetric line Gauss-Seidel splitting of the equation system Jacobian matrix. Convergence characteristics are improved by the addition of acceleration techniques based on Shamanskii's method for nonlinear equations and Broyden's quasi-Newton update. Characteristic-based differencing of the equations is provided by means of Van Leer's flux vector splitting. In this investigation, emphasis is placed on the fast and accurate computation of shock-wave-boundary layer interactions with and without slot suction effects. In the latter context, a set of numerical boundary conditions for simulating the transpiration flow in an open slot is devised. Both laminar and turbulent cases are considered, with turbulent closure provided by a modified Cebeci-Smith algebraic model. Comparisons with computational and experimental data sets are presented for a variety of interactions, and a fully-coupled simulation of a plenum chamber/inlet flowfield with shock interaction and suction is also shown and discussed.

Edwards, Jack R.↗
NASA NTRS2002

The Case for General Relativistic Effects in the Fe K(alpha) Profile of an Active Galaxy

We present results from a simultaneous Chandra HETG (High Energy Transmission Grating) and XMM (X-ray Multi-mirror Mission)-Newton observation of NGC 3516. We find evidence for several narrow components of Fe K(alpha) along with a broad line. We consider the possibility that the lines arise in a blob of material ejected from the nucleus with velocity of approximately 0.25c. We also consider an origin in a neutral accretion disk, suffering enhanced illumination at 35 and 175 R(sub g), perhaps due to magnetic reconnection. The presence of these narrow features indicates there is no Comptonizing region along the line-of-sight to the nucleus. This in turn is compelling support for the hypothesis that broad Fe K(alpha) components are, in general, produced by strong gravity.

Turner, T. J.↗
NASA NTRS1978

Effect of filtration on rolling-element-bearing life in contaminated lubricant environment

Fatigue tests were conducted on groups of 65 millimeter-bore ball bearings under four levels of filtration with and without a contaminated MIL-L-23699 lubricant. The baseline series used noncontaminated oil with 49 micron absolute filtration. In the remaining tests contaminants of the composition found in aircraft engine filters were injected into the filter's supply line at a constant rate of 125 milligrams per bearing-hour. The test filters had absolute particle removal ratings of 3, 30, 49, and 105 microns (0.45, 10, 30, and 70 microns nominal), respectively. Bearings were tested at 15,000 rpm under 4580 newtons radial load. Bearing life and running tract condition generally improved with finer filtration. The 3 and 30 micron filter bearings in a contaminated lubricant had statistically equivalent lives, approaching those from the baseline tests. The experimental lives of 49 micron bearings were approximately half the baseline bearing's lives. Bearings tested with the 105 micron filter experienced wear failures. The degree of surface distress, weight loss, and probable failure mode were found to be dependent on filtration level, with finer filtration being clearly beneficial.

Loewenthal, S. H.↗
NASA NTRS1991

A numerical algorithm for optimal feedback gains in high dimensional linear quadratic regulator problems

A hybrid method for computing the feedback gains in linear quadratic regulator problem is proposed. The method, which combines use of a Chandrasekhar type system with an iteration of the Newton-Kleinman form with variable acceleration parameter Smith schemes, is formulated to efficiently compute directly the feedback gains rather than solutions of an associated Riccati equation. The hybrid method is particularly appropriate when used with large dimensional systems such as those arising in approximating infinite-dimensional (distributed parameter) control systems (e.g., those governed by delay-differential and partial differential equations). Computational advantages of the proposed algorithm over the standard eigenvector (Potter, Laub-Schur) based techniques are discussed, and numerical evidence of the efficacy of these ideas is presented.

Banks, H. T.↗
NASA NTRS2012

A Chandra Observation of the Ultraluminous Infrared Galaxy IRAS 19254-7245 (The Superantennae): X-Ray Emission from the Compton-Thick Active Galactic Nucleus and the Diffuse Starburst

We present a Chandra observation of IRAS 19254-7245, a nearby ultraluminous infrared galaxy also known as the Superantennae. The high spatial resolution of Chandra allows us to disentangle for the first time the diffuse starburst (SB) emission from the embedded Compton-thick active galactic nucleus (AGN) in the southern nucleus. No AGN activity is detected in the northern nucleus. The 2-10 keV spectrum of the AGN emission is fitted by a flat power law (TAU = 1.3) and an He-like Fe Kalpha line with equivalent width 1.5 keV, consistent with previous observations. The Fe K line profile could be resolved as a blend of a neutral 6.4 keV line and an ionized 6.7 keV (He-like) or 6.9 keV (H-like) line. Variability of the neutral line is detected compared with the previous XMM-Newton and Suzaku observations, demonstrating the compact size of the iron line emission. The spectrum of the galaxy-scale extended emission excluding the AGN and other bright point sources is fitted with a thermal component with a best-fit kT of approximately 0.8 keV. The 2-10 keV luminosity of the extended emission is about one order of magnitude lower than that of the AGN. The basic physical and structural properties of the extended emission are fully consistent with a galactic wind being driven by the SB. A candidate ultraluminous X-ray source is detected 8 south of the southern nucleus. The 0.3 - 10 keV luminosity of this off-nuclear point source is approximately 6 x 10(exp 40) erg per second if the emission is isotropic and the source is associated with the Superantennae.

Jia, Jianjun↗
DOE OSTI2022

Inexact Newton-CG algorithms with complexity guarantees

Abstract We consider variants of a recently developed Newton-CG algorithm for nonconvex problems (Royer, C. W. & Wright, S. J. (2018) Complexity analysis of second-order line-search algorithms for smooth nonconvex optimization. SIAM J. Optim., 28, 1448–1477) in which inexact estimates of the gradient and the Hessian information are used for various steps. Under certain conditions on the inexactness measures, we derive iteration complexity bounds for achieving $\epsilon $-approximate second-order optimality that match best-known lower bounds. Our inexactness condition on the gradient is adaptive, allowing for crude accuracy in regions with large gradients. We describe two variants of our approach, one in which the step size along the computed search direction is chosen adaptively, and another in which the step size is pre-defined. To obtain second-order optimality, our algorithms will make use of a negative curvature direction on some steps. These directions can be obtained, with high probability, using the randomized Lanczos algorithm. In this sense, all of our results hold with high probability over the run of the algorithm. We evaluate the performance of our proposed algorithms empirically on several machine learning models. Our approach is a first attempt to introduce inexact Hessian and/or gradient information into the Newton-CG algorithm of Royer & Wright (2018, Complexity analysis of second-order line-search algorithms for smooth nonconvex optimization. SIAM J. Optim., 28, 1448–1477).

Mathematics↗
NASA NTRS1987

Recursive inverse kinematics for robot arms via Kalman filtering and Bryson-Frazier smoothing

This paper applies linear filtering and smoothing theory to solve recursively the inverse kinematics problem for serial multilink manipulators. This problem is to find a set of joint angles that achieve a prescribed tip position and/or orientation. A widely applicable numerical search solution is presented. The approach finds the minimum of a generalized distance between the desired and the actual manipulator tip position and/or orientation. Both a first-order steepest-descent gradient search and a second-order Newton-Raphson search are developed. The optimal relaxation factor required for the steepest descent method is computed recursively using an outward/inward procedure similar to those used typically for recursive inverse dynamics calculations. The second-order search requires evaluation of a gradient and an approximate Hessian. A Gauss-Markov approach is used to approximate the Hessian matrix in terms of products of first-order derivatives. This matrix is inverted recursively using a two-stage process of inward Kalman filtering followed by outward smoothing. This two-stage process is analogous to that recently developed by the author to solve by means of spatial filtering and smoothing the forward dynamics problem for serial manipulators.

Rodriguez, G.↗
NASA NTRS1994

Numerical solution of the two-dimensional time-dependent incompressible Euler equations

A numerical method is presented for solving the artificial compressibility form of the 2D time-dependent incompressible Euler equations. The approach is based on using an approximate Riemann solver for the cell face numerical flux of a finite volume discretization. Characteristic variable boundary conditions are developed and presented for all boundaries and in-flow out-flow situations. The system of algebraic equations is solved using the discretized Newton-relaxation (DNR) implicit method. Numerical results are presented for both steady and unsteady flow.

Whitfield, David L.↗