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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Microwave remote sensing and its application to soil moisture detection

The author has identified the following significant results. Experimental measurements were utilized to demonstrate a procedure for estimating soil moisture, using a passive microwave sensor. The investigation showed that 1.4 GHz and 10.6 GHz can be used to estimate the average soil moisture within two depths; however, it appeared that a frequency less than 10.6 GHz would be preferable for the surface measurement. Average soil moisture within two depths would provide information on the slope of the soil moisture gradient near the surface. Measurements showed that a uniform surface roughness similar to flat tilled fields reduced the sensitivity of the microwave emission to soil moisture changes. Assuming that the surface roughness was known, the approximate soil moisture estimation accuracy at 1.4 GHz calculated for a 25% average soil moisture and an 80% degree of confidence, was +3% and -6% for a smooth bare surface, +4% and -5% for a medium rough surface, and +5.5% and -6% for a rough surface.

Newton, R. W.↗
NASA NTRS2003

XMM-Newton X-ray Observatory Guest Observer program (AO-1) at CASA

In this research program, we obtained and analyzed X-ray observations of the Wolf-Rayet (WR) star WR 110 (HD 165688) using the XMM-Newton space-based observatory. Radio observations were also obtained using the Very Large Array (VLA) radio telescope located in New Mexico and operated by the Natl. Radio Astronomy Observatory (NRAO). This star was targeted for observations primarily because it is believed to be a single WR star without a companion. Single WR stars are thought to emit X-rays from cool plasma in shocks distributed throughout their powerful stellar winds. However, there has been little observational work done to test this idea since single WR stars are relatively weak X-ray sources and have been difficult to detect with previous generation telescopes. The launch of XMM-Newton provides a new telescope that is much more sensitive than its predecessors, allowing single WR stars to be studied in detail for the first time. X-ray emission was clearly detected from WR 110. Analysis of its spectrum yields a surprising result. Its X-ray emitting plasma is distributed over a range of temperatures and is dominated by relatively cool plasma with a characteristic temperature T is approximately 6 million K. Such plasma can be explained by existing theoretical wind shock models. However, the spectrum also shows hotter plasma whose temperature is uncertain but is thought to be in excess of T approximately 30 million K. The origin of this hotter plasma is yet unknown, but possible mechanisms are identified

Skinner, Stephen L.↗
DOE OSTI2024

Revisiting Artifacts of Kohn–Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Computation of the shock-wave boundary layer interaction with flow separation

The boundary layer concept is used to describe the flow near the wall. The external flow is approximated by a pressure displacement relationship (tangent wedge in linearized supersonic flow). The boundary layer equations are solved in finite difference form and the question of the presence and unicity of the solution is considered for the direct problem (assumed pressure) or converse problem (assumed displacement thickness, friction ratio). The coupling algorithm presented implicitly processes the downstream boundary condition necessary to correctly define the interacting boundary layer problem. The algorithm uses a Newton linearization technique to provide a fast convergence.

Ardonceau, P.↗
NASA NTRS2004

The X-Ray Weakness of GPS Radio Galaxies: A Volume-Limited Complete Sample

The XMM observations of Mkn 668 have been analyzed. We found soft X-ray signatures of a hot plasma (kT approximately 10^7 approximately K) and a hard X-ray emission from the nucleus. The X-ray spectrum above 2.5 approximately keV is characterized by a very flat (observed photon index, Gamma approximately 0.5) power-law continuum, alongside with a strong Fe-K-alpha neutral iron fluorescent line (EW approximately 600 approximately eV). The best explanation for the origin of this high energy X-ray emission is in terms of the Compton-reflection of the nuclear emission. The primary X-ray emission is obscured by a Compton-thick (N_H approximately 10^24 approximately cm-2) matter which becomes transparent at higher energies. The observed above 2.5-keV X-rays are mostly due to reflection which is indicated by a strong Fe-K-alpha line. This represents the second hard X-ray detection of the GPS galaxy ever (the first one being 1345+125; O Dea et al. 2000). Interestingly, the both such trend is confirmed by our on going XMM-Newton observations of a larger GPS sample, it would lead us to looking into the question on how the dense nuclear environment impacts the nature and evolution of a GPS source, and more generally, on the history of radio power in the universe. The paper summarizing the results has been submitted to Astronomy and Astrophysics in December 2003.

Mushotzky, Richard F.↗
DOE OSTI2023

Solution and sensitivity analysis of nonlinear equations using a hypercomplex-variable Newton-Raphson method

Here, the classical Newton-Raphson (NR) method for solving nonlinear equations is enhanced in two ways through the use of hypercomplex variables and algebra. In particular, i) the Jacobian is computed in a highly accurate and automated way, and ii) the derivative of the solution to the nonlinear equations is computed with respect to any parameter contained within the system of equations. These advances provide two significant enhancements in that it is straightforward to provide an accurate Jacobian and to construct a reduced order model (ROM) of arbitrary order with respect to any parameter of the system. The ROM can then be used to approximate the solution for other parameter values without requiring additional solutions of the nonlinear equations. Several case studies are presented including 1D and 2D academic examples with fully functioning Python code provided. Additionally, a case of study of the catenary of an elastic cable subject to its own weight and a vertical point load. Derivatives up to 10th order were computed with respect to material, loading, and geometrical parameters. The derivatives were used to generate reduced order models of the cable deformation and reaction forces at its ends with respect to multiple input parameters. Results show that from a single hypercomplex evaluation of the cable under a single vertical point load, it is possible to generate an accurate reduced order model capable of predicting the cable deformation with 1.5 times the load in the opposite direction and with 3.5 times the load in the same direction without resolving the system of equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Solution of free-boundary problems using finite-element/Newton methods and locally refined grids - Application to analysis of solidification microstructure

A new method is presented for the solution of free-boundary problems using Lagrangian finite element approximations defined on locally refined grids. The formulation allows for direct transition from coarse to fine grids without introducing non-conforming basis functions. The calculation of elemental stiffness matrices and residual vectors are unaffected by changes in the refinement level, which are accounted for in the loading of elemental data to the global stiffness matrix and residual vector. This technique for local mesh refinement is combined with recently developed mapping methods and Newton's method to form an efficient algorithm for the solution of free-boundary problems, as demonstrated here by sample calculations of cellular interfacial microstructure during directional solidification of a binary alloy.

Tsiveriotis, K.↗
NASA NTRS2004

Complex X-ray Absorption and the Fe K(alpha) Profile in NGC 3516

We present data from simultaneous Chandra, XMM-Newton and BeppoSAX observations of the Seyfert 1 galaxy NGC 3516, taken during 2001 April and November. We have investigated the nature of the very flat observed X-ray spectrum. Chandra grating data show the presence of X-ray absorption lines, revealing two distinct components of the absorbing gas, one which is consistent with our previous model of the UV/X-ray absorber while the other, which is outflowing at a velocity of approximately 1100 kilometers per second, has a larger column density and is much more highly ionized. The broad-band spectral characteristics of the X-ray continuum observed with XMM during 2001 April, reveal the presence of a third layer of absorption consisting of a very large column (approximately 2.5 x 10(exp 23) per square centimeter) of highly ionized gas with a covering fraction approximately 50%. This low covering fraction suggests that the absorber lies within a few 1t-days of the X-ray source and/or is filamentary in structure. Interestingly, these absorbers are not in thermal equilibrium with one another. The two new components are too highly ionized to be radiatively accelerated, which we suggest is evidence for a hydromagnetic origin for the outflow. Applying our model to the November dataset, we can account for the spectral variability primarily by a drop in the ionization states of the absorbers, as expected by the change in the continuum flux. When this complex absorption is accounted for we find the underlying continuum to be typical of Seyfert 1 galaxies. The spectral curvature attributed to the high column absorber, in turn, reduces estimates of the flux and extent of any broad Fe emission line from the accretion disk.

Turner, T. J.↗
NASA NTRS2012

Peak Seeking Control for Reduced Fuel Consumption with Preliminary Flight Test Results

The Environmentally Responsible Aviation project seeks to accomplish the simultaneous reduction of fuel burn, noise, and emissions. A project at NASA Dryden Flight Research Center is contributing to ERAs goals by exploring the practical application of real-time trim configuration optimization for enhanced performance and reduced fuel consumption. This peak-seeking control approach is based on Newton-Raphson algorithm using a time-varying Kalman filter to estimate the gradient of the performance function. In real-time operation, deflection of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of a modified F-18 are directly optimized, and the horizontal stabilators and angle of attack are indirectly optimized. Preliminary results from three research flights are presented herein. The optimization system found a trim configuration that required approximately 3.5% less fuel flow than the baseline trim at the given flight condition. The algorithm consistently rediscovered the solution from several initial conditions. These preliminary results show the algorithm has good performance and is expected to show similar results at other flight conditions and aircraft configurations.

Brown, Nelson↗
NASA NTRS2002

High Resolution X-ray Imaging

NAG5-5020 covered a period of 7.5 years during which a great deal of progress was made in x-ray optical techniques under this grant. We survived peer review numerous times during the effort to keep the grant going. In 1994, when the grant started we were actively pursuing the application of spherical mirrors to improving x-ray telescopes. We had found that x-ray detectors were becoming rapidly more sophisticated and affordable, but that x-ray telescopes were only being improved through the intense application of money within the AXAF program. Clearly new techniques for the future were needed. We were successful in developing and testing at the HELSTF facility in New Mexico a four reflection coma-corrected telescope made from spheres. We were able to demonstrate 0.3 arcsecond resolution, almost to the diffraction limit of the system. The community as a whole was, at that time, not particularly interested in looking past AXAF (Chandra) and the effort needed to evolve. Since we had reached the diffraction limit using non-Wolter optics we then decided to see if we could build an x-ray interferometer in the laboratory. In the lab the potential for improved resolution was substantial. If synthetic aperture telescopes could be built in space, then orders of magnitude improvement would become feasible. In 1998 NASA, under the direction of Dr Nick White of Goddard, started a study to assess the potential and feasibility of x-ray interferometry in space. My work became of central interest to the committee because it indicated that such was possible. In early 1999 we had the breakthrough that allowed us build a practical interferometer. By using flats and hooking up with the Marshall Space Flight Center facilities we were able to demonstrate fringes at 1.25keV on a one millimeter baseline. This actual laboratory demonstration provided the solid proof of concept that NASA needed. As the year progressed the future of x-ray astronomy jelled around the Maxim program. Maxim is a set of two major x-ray astronomy missions based on the concepts I developed and demonstrated under this SR&T grant. The first Maxim is to image the sky at 100 micro-arcsecond resolution. That is one thousand times higher resolution than Hubble. The full Maxim has the ultimate goal of imaging the event horizon of a black hole in an active galactic nucleus (ALAN). This will require 0.1 micro-arcsecond resolution - one million times higher than Hubble! Nonetheless, using the techniques developed under this grant, it has become possible. Maxim Pathfinder is now in the NASA planning for a new start in approximately 20 10. The full Maxim is carried as a vision mission for the post 2015 timeframe. Finally, this grant is the evolved version of the SR&T grant we carried during the 1980s and up to 1994. At that point in time this grant was also working on x-ray optics, but concentrating on x-ray spectroscopy. The techniques developed by 1990 were not chosen for use on Chandra or XMM-Newton because they were too new. During the last year, however, the Constellation-X mission recognized the need for better spectroscopy techniques and tapped our expertise. We were able to support the initial work on Con-X through this program. It now appears that the off-plane mount will be used in Con-X, increasing performance and decreasing cost and risk.

Cash, Webster↗
NASA NTRS1988

Domain decomposition methods for the parallel computation of reacting flows

Domain decomposition is a natural route to parallel computing for partial differential equation solvers. Subdomains of which the original domain of definition is comprised are assigned to independent processors at the price of periodic coordination between processors to compute global parameters and maintain the requisite degree of continuity of the solution at the subdomain interfaces. In the domain-decomposed solution of steady multidimensional systems of PDEs by finite difference methods using a pseudo-transient version of Newton iteration, the only portion of the computation which generally stands in the way of efficient parallelization is the solution of the large, sparse linear systems arising at each Newton step. For some Jacobian matrices drawn from an actual two-dimensional reacting flow problem, comparisons are made between relaxation-based linear solvers and also preconditioned iterative methods of Conjugate Gradient and Chebyshev type, focusing attention on both iteration count and global inner product count. The generalized minimum residual method with block-ILU preconditioning is judged the best serial method among those considered, and parallel numerical experiments on the Encore Multimax demonstrate for it approximately 10-fold speedup on 16 processors.

Keyes, David E.↗
NASA NTRS1991

Development of a steady potential solver for use with linearized, unsteady aerodynamic analyses

A full potential steady flow solver (SFLOW) developed explicitly for use with an inviscid unsteady aerodynamic analysis (LINFLO) is described. The steady solver uses the nonconservative form of the nonlinear potential flow equations together with an implicit, least squares, finite difference approximation to solve for the steady flow field. The difference equations were developed on a composite mesh which consists of a C grid embedded in a rectilinear (H grid) cascade mesh. The composite mesh is capable of resolving blade to blade and far field phenomena on the H grid, while accurately resolving local phenomena on the C grid. The resulting system of algebraic equations is arranged in matrix form using a sparse matrix package and solved by Newton's method. Steady and unsteady results are presented for two cascade configurations: a high speed compressor and a turbine with high exit Mach number.

Hoyniak, Daniel↗
NASA NTRS1996

Comments on SEE: Comparative Advantages and and Experimental Consequences

The Satellite Energy Exchange experiment measures the periodic, near-miss encounters between a sheppard satellite and a small test body (satellite) in approximately the same orbit about a primary. Several important experimental requirements have been chosen to enhance capabilities: (a) The satellite be flown in a sun-synchronous orbit at an altitude of about 1350 Km, (b) Passive temperature system stabilized by spacecraft axial rotation with sunshade baffles at the end of the spacecraft, (c) Test bodies with different material composition be available for experiments, (d) The containment spacecraft fly about the sheppard mass in a zero-g environment whereas the test bodies, experience average zero-g environment over an orbital period, (e) Primary attitude and station-keeping uses magnetic field alignment plus micro-Newton thrusters such as Field Emission Electric Propulsion, and (f) Very low power (nW) laser tracking systems minimize impulse delivered to test bodies. With the above conditions, SEE has the capabilities: (1) Long duration (several years life-time) flight experiment (2) Long-term, active (with historical time record), self-calibration of satellite mass distribution (capsule geodesy) over lifetime of the spacecraft. (3) Novel passive thermal stabilization systems designed to attain cryogenic temperatures around 78K. (4) Novel spacecraft stabilization systems. (5) Ability to measure G to 1 part in 10(exp 6-7) depending on ultimate duration of experiment. (6) Ability to place limits on both temporal and spacial variations on G. (7) Ability to set experimental limits on the Post Newtonian parameters (PPN) alpha(2) and zeta(2). (8) Ability to measure (or place limits on) the non Einsteinian eccentricity of the Earth-Sun system (and the parameter alpha(1)) for long duration flight. (9) Ability to measure Delta((dot)-G)/G to 1 part in 10(exp 12-13). The MiniSTEP, competes in a limited way with Project SEE. It is designed to improve the measurement of the equivalence principle by seven orders of magnitude using active, low temperature (1.8 K) cooling for SQUID based, differential superconducting circuits. The experiment consists of a small cylinder concentrically located within a larger cylinder at its null gravitational point. The satellite is operated in zero-g mode using four differential accelerometers consisting to two test bodies of different material composition. The SQUIDS are needed to measure test body motion to precisions of 10(exp -18) over a four orbit period. The entire satellite moves in a very precise zero-g mode since the accelerometers are rigidly attached to the satellite. This limits the experiment to an approximately six month due to limitations on helium storage used in cryogenic cooling and thrust control to maintain the zero-g operation.

Smalley, Larry L.↗
DOE OSTI2025

Methods for Measuring and Computing the Reference Temperature in Newton’s Law of Cooling for External Flows

Newton’s law of cooling requires a reference temperature (𝑇 𝑟𝑒𝑓 ) to define the heat-transfer coefficient (ℎ). For external flows with multiple temperatures in the freestream, obtaining 𝑇 𝑟𝑒𝑓 is a challenge. One widely used method, referred to as the adiabatic-wall (AW) method, obtains 𝑇 𝑟𝑒𝑓 by requiring the surface of the solid exposed to convective heat transfer to be adiabatic. Another widely used method, referred to as the linear-extrapolation (LE) method, obtains 𝑇 𝑟𝑒𝑓 by measuring/computing the heat flux ($𝑞^{′′}_𝑠$) on the solid surface at two different surface temperatures (𝑇 𝑠 ) and then linearly extrapolating to $𝑞^{′′}_𝑠$ = 0. A third recently developed method, referred to as the state-space (SS) method, obtains 𝑇 𝑟𝑒𝑓 by probing the temperature space between the highest and lowest in the flow to account for the effects of 𝑇 𝑠 or $𝑞^{′′}_𝑠$ on 𝑇 𝑟𝑒𝑓 . This study examines the foundation and accuracy of these methods via a test problem involving film cooling of a flat plate where $𝑞^{′′}_𝑠$ switches signs on the plate’s surface. Results obtained show that only the SS method could guarantee a unique and physically meaningful 𝑇 𝑟𝑒𝑓 where 𝑇 𝑠 =𝑇 𝑟𝑒𝑓 on a nonadiabatic surface $𝑞^{′′}_𝑠$ = 0. The AW and LE methods both assume 𝑇 𝑟𝑒𝑓 to be independent of 𝑇 𝑠 , which the SS method shows to be incorrect. Though this study also showed the adiabatic-wall temperature, 𝑇 𝐴𝑊 , to be a good approximation of 𝑇 𝑟𝑒𝑓 (<10% relative error), huge errors can occur in ℎ about the solid surface where |𝑇 𝑠 −𝑇 𝐴𝑊 | is near zero because where 𝑇 𝑠 =𝑇 𝐴𝑊 , $𝑞^{′′}_𝑠$ ≠ 0.

Newton's law of cooling↗
NASA NTRS1984

Numerical optimization design of advanced transonic wing configurations

A computationally efficient and versatile technique for use in the design of advanced transonic wing configurations has been developed. A reliable and fast transonic wing flow-field analysis program, TWING, has been coupled with a modified quasi-Newton method, unconstrained optimization algorithm, QNMDIF, to create a new design tool. Fully three-dimensional wing designs utilizing both specified wing pressure distributions and drag-to-lift ration minimization as design objectives are demonstrated. Because of the high computational efficiency of each of the components of the design code, in particular the vectorization of TWING and the high speed of the Cray X-MP vector computer, the computer time required for a typical wing design is reduced by approximately an order of magnitude over previous methods. In the results presented here, this computed wave drag has been used as the quantity to be optimized (minimized) with great success, yielding wing designs with nearly shock-free (zero wave drag) pressure distributions and very reasonable wing section shapes.

Cosentino, G. B.↗
NASA NTRS1985

Numerical optimization design of advanced transonic wing configurations

A computationally efficient and versatile technique for use in the design of advanced transonic wing configurations has been developed. A reliable and fast transonic wing flow-field analysis program, TWING, has been coupled with a modified quasi-Newton method, unconstrained optimization algorithm, QNMDIF, to create a new design tool. Fully three-dimensional wing designs utilizing both specified wing pressure disributions and drag-to-lift ratio minimization as design objectives are demonstrated. Because of the high computational efficiency of each of the components of the design code, in particular the vectorization of TWING and the high speed of the Cray X-MP vector computer, the computer time required for a typical wing design is reduced by approximately an order of magnitude over previous methods. In the results presented here, this computed wave drag has been used as the quantity to be optimized (minimized) with great success, yielding wing designs with nearly shock-free (zero wave drag) pressure distributions and very reasonable wing section shapes.

Cosentino, G. B.↗
NASA NTRS2004

The Sombrero: An Edge-On Spiral: Halo Dynamics, A Massive Dark Object, and XRBS

We present an analysis of the XMM-Newton and Chandra ACIS-S observations of the LINER nucleus of the Sombrero galaxy, and we discuss possible explanations for its very Sub-Eddington luminosity by complementing the X-ray results with high angular resolution observations in other bands. The X-ray investigation shows a hard (gamma=1.89) and moderately absorbed (NH=1.8x1021 per square centimeter) nuclear source of 1.5x1040 ergs per second in the 2-10 keV band, surrounded by hot gas at a temperature of approximately 0.6 keV. The bolometric nuclear luminosity is at least approximately 200 times lower than expected if mass accreted on the supermassive black hole, which Hubble Space Telescope shows to reside at the center of this galaxy, at the rate predicted by the spherical and adiabatic Bondi accretion theory and with the high radiative efficiency of a standard accretion disk. The low luminosity, coupled to the observed absence of Fe K emission in the nuclear spectrum, indicates that such a disk is not present. This nucleus also differs from bright unobscured active galactic nuclei in the lack of high flux variability and prominent broad H emission. However, it is also too faint for the predictions of simple radiatively inefficient accretion taking place at the Bondi rate; it could instead be too radio bright for radiatively inefficient accretion that includes strong mass outflows or convection. This discrepancy could be solved by the possible presence of nuclear radio jets. An alternative explanation of the low luminosity, in place of radiative inefficiency, could be unsteady accretion.

Fabbiano, G.↗
NASA NTRS1995

A General Method for Solving Systems of Non-Linear Equations

The method of steepest descent is modified so that accelerated convergence is achieved near a root. It is assumed that the function of interest can be approximated near a root by a quadratic form. An eigenvector of the quadratic form is found by evaluating the function and its gradient at an arbitrary point and another suitably selected point. The terminal point of the eigenvector is chosen to lie on the line segment joining the two points. The terminal point found lies on an axis of the quadratic form. The selection of a suitable step size at this point leads directly to the root in the direction of steepest descent in a single step. Newton's root finding method not infrequently diverges if the starting point is far from the root. However, the current method in these regions merely reverts to the method of steepest descent with an adaptive step size. The current method's performance should match that of the Levenberg-Marquardt root finding method since they both share the ability to converge from a starting point far from the root and both exhibit quadratic convergence near a root. The Levenberg-Marquardt method requires storage for coefficients of linear equations. The current method which does not require the solution of linear equations requires more time for additional function and gradient evaluations. The classic trade off of time for space separates the two methods.

Nachtsheim, Philip R.↗