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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Cascade Optimization for Aircraft Engines With Regression and Neural Network Analysis - Approximators

The NASA Engine Performance Program (NEPP) can configure and analyze almost any type of gas turbine engine that can be generated through the interconnection of a set of standard physical components. In addition, the code can optimize engine performance by changing adjustable variables under a set of constraints. However, for engine cycle problems at certain operating points, the NEPP code can encounter difficulties: nonconvergence in the currently implemented Powell's optimization algorithm and deficiencies in the Newton-Raphson solver during engine balancing. A project was undertaken to correct these deficiencies. Nonconvergence was avoided through a cascade optimization strategy, and deficiencies associated with engine balancing were eliminated through neural network and linear regression methods. An approximation-interspersed cascade strategy was used to optimize the engine's operation over its flight envelope. Replacement of Powell's algorithm by the cascade strategy improved the optimization segment of the NEPP code. The performance of the linear regression and neural network methods as alternative engine analyzers was found to be satisfactory. This report considers two examples-a supersonic mixed-flow turbofan engine and a subsonic waverotor-topped engine-to illustrate the results, and it discusses insights gained from the improved version of the NEPP code.

Patnaik, Surya N.↗
NASA NTRS1989

Recursive linearization of multibody dynamics equations of motion

The equations of motion of a multibody system are nonlinear in nature, and thus pose a difficult problem in linear control design. One approach is to have a first-order approximation through the numerical perturbations at a given configuration, and to design a control law based on the linearized model. Here, a linearized model is generated analytically by following the footsteps of the recursive derivation of the equations of motion. The equations of motion are first written in a Newton-Euler form, which is systematic and easy to construct; then, they are transformed into a relative coordinate representation, which is more efficient in computation. A new computational method for linearization is obtained by applying a series of first-order analytical approximations to the recursive kinematic relationships. The method has proved to be computationally more efficient because of its recursive nature. It has also turned out to be more accurate because of the fact that analytical perturbation circumvents numerical differentiation and other associated numerical operations that may accumulate computational error, thus requiring only analytical operations of matrices and vectors. The power of the proposed linearization algorithm is demonstrated, in comparison to a numerical perturbation method, with a two-link manipulator and a seven degrees of freedom robotic manipulator. Its application to control design is also demonstrated.

Lin, Tsung-Chieh↗
DOE OSTI2020

A least-squares finite element method based on the Helmholtz decomposition for hyperbolic balance laws

In this paper, a least-squares finite element method for scalar nonlinear hyperbolic balance laws is proposed and studied. The approach is based on a formulation that utilizes an appropriate Helmholtz decomposition of the flux vector and is related to the standard notion of a weak solution. This relationship, together with a corresponding connection to negative-norm least-squares, is described in detail. As a consequence, an important numerical conservation theorem is obtained, similar to the famous Lax–Wendroff theorem. The numerical conservation properties of the method in this paper do not fall precisely in the framework introduced by Lax and Wendroff, but they are similar in spirit as they guarantee that when L 2 convergence holds, the resulting approximations approach a weak solution to the hyperbolic problem. The least-squares functional is continuous and coercive in an H -1 -type norm, but not L 2 -coercive. Nevertheless, the L 2 convergence properties of the method are discussed. Convergence can be obtained either by an explicit regularization of the functional, that provides control of the L 2 norm, or by properly choosing the finite element spaces, providing implicit control of the L 2 norm. Numerical results for the inviscid Burgers equation with discontinuous source terms are shown, demonstrating the L 2 convergence of the obtained approximations to the physically admissible solution. The numerical method utilizes a least-squares functional, minimized on finite element spaces, and a Gauss–Newton technique with nested iteration. Finally, we believe that the linear systems encountered with this formulation are amenable to multigrid techniques and combining the method with adaptive mesh refinement would make this approach an efficient tool for solving balance laws (this is the focus of a future study).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 3rd, California State University, Long Beach, CA, January 21-24, 1985, Proceedings

Among the topics discussed are hydrodynamic stability and turbulent transition, separated flow simulations employing finite difference methods, a quasi-simultaneous finite difference approach for strongly interacting flow, the Newton solution of coupled Euler and boundary layer equations, multielement airfoil aerodynamics, vorticity dynamics, the significance of the thin layer Navier-Stokes approximation, a comparison of interactive and Navier-Stokes calculations of separating boundary layer flows, and a Navier-Stokes analysis for the mechanism of shock wave/boundary layer interactions. Also discussed are transonic interactions of unsteady vortical flows, boundary layers on oscillating airfoils, a laser-interferometer skin friction meter, flow over surface-mounted semibluff bodies, computations for spinning and spin-stabilized supersonic projectiles, transition calculations in three-dimensional flows, the effect of projectile boattailing in transonic flow, flow characteristics of a body of revolution at incidence, the influence of the boundary layer on the leading edge flow of a swept wing, and an interactive scheme for three-dimensional transonic flows.

Source record↗
NASA NTRS1987

Analysis of a parallelized nonlinear elliptic boundary value problem solver with application to reacting flows

A parallelized finite difference code based on the Newton method for systems of nonlinear elliptic boundary value problems in two dimensions is analyzed in terms of computational complexity and parallel efficiency. An approximate cost function depending on 15 dimensionless parameters is derived for algorithms based on stripwise and boxwise decompositions of the domain and a one-to-one assignment of the strip or box subdomains to processors. The sensitivity of the cost functions to the parameters is explored in regions of parameter space corresponding to model small-order systems with inexpensive function evaluations and also a coupled system of nineteen equations with very expensive function evaluations. The algorithm was implemented on the Intel Hypercube, and some experimental results for the model problems with stripwise decompositions are presented and compared with the theory. In the context of computational combustion problems, multiprocessors of either message-passing or shared-memory type may be employed with stripwise decompositions to realize speedup of O(n), where n is mesh resolution in one direction, for reasonable n.

Keyes, David E.↗
DOE OSTI2023

Discovery of the luminous X-ray ignition eRASSt J234402.9-352640: I. Tidal disruption event or a rapid increase in accretion in an active galactic nucleus?

In November 2020, a new, bright object, eRASSt J234402.9-352640, was discovered in the second all-sky survey of SRG/eROSITA. The object brightened by a factor of at least 150 in 0.2–2.0 keV flux compared to an upper limit found six months previous, reaching an observed peak of 1.76 -0.24 +0.03 × 10 -11 erg cm -2 s -1 . The X-ray ignition is associated with a galaxy at z = 0.10, making the peak luminosity log 10 (L 0.2-2 keV /[erg s -1 ]) = 44.7 ± 0.1. Around the time of the rise in X-ray flux, the nucleus of the galaxy brightened by approximately 3 mag. in optical photometry, after correcting for the host contribution. We present X-ray follow-up data from Swift, XMM-Newton, and NICER, which reveal a very soft spectrum as well as strong 0.2–2.0 keV flux variability on multiple timescales. Optical spectra taken in the weeks after the ignition event show a blue continuum with broad, asymmetric Balmer emission lines, and high-ionisation ([OIII]λλ4959,5007) and low-ionisation ([NII]λ6585, [SII]λλ6716,6731) narrow emission lines. Following the peak in the optical light curve, the X-ray, UV, and optical photometry all show a rapid decline. The X-ray light curve shows a decrease in luminosity of ~0.45 over 33 days and the UV shows a drop of ~0.35 over the same period. eRASSt J234402.9-352640 also shows a brightening in the mid-infrared, likely powered by a dust echo of the luminous ignition. We find no evidence in Fermi-LAT γ-ray data for jet-like emission. The event displays characteristics of a tidal disruption event (TDE) as well as of an active galactic nucleus (AGN), complicating the classification of this transient. Based on the softness of the X-ray spectrum, the presence of high-ionisation optical emission lines, and the likely infrared echo, we find that a TDE within a turned-off AGN best matches our observations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1996

Practical Aspects of Krylov Subspace Iterative Methods in CFD

Implementation issues associated with the application of Krylov subspace iterative methods, such as Newton-GMRES, are presented within the framework of practical computational fluid dynamic (CFD) applications. This paper categorizes, evaluates, and contrasts the major ingredients (function evaluations, matrix-vector products, and preconditioners) of Newton-GMRES Krylov subspace methods in terms of their effect on the local linear and global nonlinear convergence, memory requirements, and accuracy. The discussion focuses on Newton-GMRES in both a structured multi-zone incompressible Navier-Stokes solver and an unstructured mesh finite-volume Navier-Stokes solver. Approximate versus exact matrix-vector products, effective preconditioners, and other pertinent issues are addressed.

Pulliam, Thomas H.↗
NASA NTRS2005

XMM-Newton Spectroscopy of the X-ray Detected Broad Absorption Line QSO CSO 755

We present the results from XMM-Newton observations of the highly optically polarized broad absorption line quasar (BALQSO) CSO 755. By analyzing its X-ray spectrum with a total of approximately 3000 photons we find that this source has an X-ray continuum of "typical" radio-quiet quasars, with a photon index of Gamma=1.83, and a rather flat (X-ray bright) intrinsic optical-to-X-ray spectral slope of alpha_ox=- 1.51. The source shows evidence for intrinsic absorption, and fitting the spectrum with a neutral-absorption model gives a column density of N_H approximately 1.2x10^22 cm^{-2}; this is among the lowest X-ray columns measured for BALQSOs. We do not detect, with high significance, any other absorption features in the X-ray spectrum. Upper limits we place on the rest-frame equivalent width of a neutral (ionized) Fe K-alpha line, less than =180 eV (less than =120 eV), and on the Compton-reflection component parameter, R less than =0.2, suggest that most of the X-rays from the source are directly observed rather than being scattered or reflected; this is also supported by the relatively flat intrinsic alpha ox we measure. The possibility that most of the X-ray flux is scattered due to the high level of UV-optical polarization is ruled out. Considering data for 46 BALQSOs from the literature, including CSO 755, we have found that the UV-optical continuum polarization level of BALQSOs is not correlated with any of their X-ray properties. A lack of significant short-term and long-term X-ray flux variations in the source may be attributed to a large black-hole mass in CSO 755. We note that another luminous BALQSO, PG 2112+059, has both similar shallow C IV BALs and moderate X-ray absorption.

Brandt, Niel↗
NASA NTRS2012

Evidence for a Circum-Nuclear and Ionised Absorber in the X-ray Obscured Broad Line Radio Galaxy 3C 445

Here we present the results of a Suzaku observation of the Broad Line Radio Galaxy 3C 445. We confirm the results obtained with the previous X-ray observations which unveiled the presence of several soft X-ray emission lines and an overall X-ray emission which strongly resembles a typical Seyfert 2 despite of the optical classification as an unobscured AGN. The broad band spectrum allowed us to measure for the first time the amount of reflection (R approximately 0.9) which together with the relatively strong neutral Fe Ka emission line (EW approximately 100 eV) strongly supports a scenario where a Compton-thick mirror is present. The primary X ray continuum is strongly obscured by an absorber with a column density of NH = 2 - 3 x 10(exp 23) per square centimeter. Two possible scenarios are proposed for the absorber: a neutral partial covering or a mildly ionised absorber with an ionisation parameter log xi approximately 1.0 erg centimeter per second. A comparison with the past and more recent X-ray observations of 3C 445 performed with XMM-Newton and Chandra is presented, which provided tentative evidence that the ionised and outflowing absorber varied. We argue that the absorber is probably associated with an equatorial diskwind located within the parsec scale molecular torus.

Braito, V.↗
DOE OSTI2024

Local convergence analysis of an inexact trust-region method for nonsmooth optimization

In Baraldi, we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function in Hilbert space—a class of problems that is ubiquitous in data science, learning, optimal control, and inverse problems. Furthermore, this algorithm has demonstrated excellent performance and scalability with problem size. In this paper, we enrich the convergence analysis for this algorithm, proving strong convergence of the iterates with guaranteed rates. In particular, we demonstrate that the trust-region algorithm recovers superlinear, even quadratic, convergence rates when using a second-order Taylor approximation of the smooth objective function term.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Arakelian 564: An XMM-Newton View

We report on two XMM-Newton observations of the bright narrow-line Seyfert 1 galaxy Ark 564 taken one year apart (2000 June and 2001 June). The 0.6-10 keV continuum is well described by a soft blackbody component (kTau approximately equal 140-150 eV) plus a steep power law (Tau approximately equal to 2.50-2.55). No significant spectral changes are observed between the two observations, although the X-ray flux in the second observation is approximately equal to 40-50 per cent lower. In both observations we detect a significant absorption edge at a rest-frame energy of approximately equal to 0.73 keV, corresponding to O VII. The presence of the absorption feature is confirmed by a simultaneous Chandra grating observation in 2000 June, although the best-fitting edge threshold is at a slightly lower energy in the Chandra data, possibly because of a different parameterization of the underlying X-ray continuum. We find tentative evidence for a broad iron emission line in the 2000 June observation. The results from an analysis of the power spectral density (PSD) function are also presented. The present XMM-Newton data support the idea that the PSD shown two breads, although the location of the high-frequency break requires further constraints.

Vignali, Cristian↗
NASA NTRS1976

Design and evaluation of a 3 million DN series-hybrid thrust bearing

The bearing, consisting of a 150-mm ball bearing and a centrifugally actuated, conical, fluid-film bearing, was fatigue tested. Test conditions were representative of a mainshaft ball bearing in a gas turbine engine operating at maximum thrust load to simulate aircraft takeoff conditions. Tests were conducted up to 16000 rpm and at this speed an axial load of 15568 newtons (3500 lb) was safely supported by the hybrid bearing system. Through the series-hybrid bearing principle, the effective ball bearing speed was reduced to approximately one-half of the shaft speed. It was concluded that a speed reduction of this magnitude results in a ten-fold increase in the ball bearing fatigue life. A successful evaluation of fluid-film bearing lubricant supply failure was performed repeatedly at an operating speed of 10,000 rpm. A complete and smooth changeover to full-scale ball bearing operation was effected when the oil supply to the fluid-film bearing was cut off. Reactivation of the fluid-film oil supply system resulted in a flawless return to the original mode of hybrid operation.

Scibbe, H. W.↗
NASA NTRS1929

Viscosity of diesel engine fuel oil under pressure

In the development of Diesel engine fuel injection systems it is necessary to have an approximate knowledge of the absolute viscosity of the fuel oil under high hydrostatic pressures. This report presents the results of experimental tests conducted by Mr. Jackson Newton Shore, utilizing the A.S.M.E. high pressure equipment.

Hersey, Mayo D↗
DOE OSTI2024

Projective Integral Updates for High-Dimensional Variational Inference

Variational inference is an approximation framework for Bayesian inference that seeks to improve quantified uncertainty in predictions by optimizing a simplified distribution over parameters to stand in for the full posterior. Capturing model variations that remain consistent with training data enables more robust predictions by reducing parameter sensitivity. This work introduces a fixed-point optimization for variational inference that is applicable when every feasible log density can be expressed as a linear combination of functions from a given basis. In such cases, the optimizer becomes a fixed-point of projective integral updates. When the basis spans univariate quadratics in each parameter, the feasible distributions are Gaussian mean-fields and the projective integral updates yield quasi-Newton variational Bayes (QNVB). Other bases and updates are also possible. Since these updates require high-dimensional integration, this work begins by proposing an efficient quasirandom sequence of quadratures for mean-field distributions. Each iterate of the sequence contains two evaluation points that combine to correctly integrate all univariate quadratic functions and, if the mean-field factors are symmetric, all univariate cubics. More importantly, averaging results over short subsequences achieves periodic exactness on a much larger space of multivariate polynomials of quadratic total degree. The corresponding variational updates require four loss evaluations with standard (not second-order) backpropagation to eliminate error terms from over half of all multivariate quadratic basis functions. Furthermore, this integration technique is motivated by first proposing stochastic blocked mean-field quadratures, which may be useful in other contexts. A PyTorch implementation of QNVB allows for better control over model uncertainty during training than competing methods. Experiments demonstrate superior generalizability for multiple learning problems and architectures.

Gaussian mean-field↗
NASA NTRS1982

Numerical methods and computers used in elastohydrodynamic lubrication

Some of the methods of obtaining approximate numerical solutions to boundary value problems that arise in elastohydrodynamic lubrication are reviewed. The highlights of four general approaches (direct, inverse, quasi-inverse, and Newton-Raphson) are sketched. Advantages and disadvantages of these approaches are presented along with a flow chart showing some of the details of each. The basic question of numerical stability of the elastohydrodynamic lubrication solutions, especially in the pressure spike region, is considered. Computers used to solve this important class of lubrication problems are briefly described, with emphasis on supercomputers.

Hamrock, B. J.↗
NASA NTRS1985

A Residuals Approach to Filtering, Smoothing and Identification for Static Distributed Systems

An approach for state estimation and identification of spatially distributed parameters embedded in static distributed (elliptic) system models is advanced. The method of maximum likelihood is used to find parameter values that maximize a likelihood functional for the system model, or equivalently, that minimize the negative logarithm of this functional. To find the minimum, a Newton-Raphson search is conducted that from an initial estimate generates a convergent sequence of parameter estimates. For simplicity, a Gauss-Markov approach is used to approximate the Hessian in terms of products of first derivatives. The gradient and approximate Hessian are computed by first arranging the negative log likelihood functional into a form based on the square root factorization of the predicted covariance of the measurement process. The resulting data processing approach, referred to here by the new term of predicted data covariance square root filtering, makes the gradient and approximate Hessian calculations very simple. A closely related set of state estimates is also produced by the maximum likelihood method: smoothed estimates that are optimal in a conditional mean sense and filtered estimates that emerge from the predicted data covariance square root filter.

Rodriguez, G.↗
NASA NTRS2002

Groups and the Entropy Floor: XMM-Newton Observations of Two Groups

Using XMM-Newton spatially resolved X-ray imaging spectroscopy we obtain the temperature, density, entropy, gas mass, and total mass profiles for two groups of galaxies out to approximately 0.3 R(sub vir)(R(sub vir), the virial radius). Our density profiles agree well with those derived previously, and the temperature data are broadly consistent with previous results but are considerably more precise. Both of these groups are at the mass scale of 2x10(exp 13) M(solar mass), but have rather different properties. Both have considerably lower gas mass fractions at r < 0.3 R(sub vir), than the rich clusters. NGC2563, one of the least luminous groups for its X-ray temperature, has a very low gas mass fraction of approximately 0.004 inside 0.1 R(sub vir), which increases with radius. NGC4325, one of the most luminous groups at the same average temperature, has a higher gas mass fraction of 0.02. The entropy profiles and the absolute values of the entropy as a function of virial radius also differ, with NGC4325 having a value of approximately 100 keV cm(exp -2) and NGC2563 a value of approximately 300 keV cm(exp -2) at r approximately 0.1 R(sub vir). For both groups the profiles rise monotonically with radius and there is no sign of an entropy 'floor'. These results are inconsistent with pre-heating scenarios that have been developed to explain a possible entropy floor in groups, but are broadly consistent with models of structure formation that include the effects of heating and/or the cooling of the gas. The total entropy in these systems provides a strong constraint on all models of galaxy and group formation, and on the poorly defined feedback process that controls the transformation of gas into stars and thus the formation of structure in the universe.

Mushotzky, R. F.↗