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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1960

Landing Characteristics of a Reentry Capsule with a Torus-Shaped Air Bag for Load Alleviation

An experimental investigation has been made to determine the landing characteristics of a conical-shaped reentry capsule by using torus-shaped air bags for impact-load alleviation. An impact bag was attached below the large end of the capsule to absorb initial impact loads and a second bag was attached around the canister to absorb loads resulting from impact on the canister when the capsule overturned. A 1/6-scale dynamic model of the configuration was tested for nominal flight paths of 60 deg. and 90 deg. (vertical), a range of contact attitudes from -25 deg. to 30 deg., and a vertical contact velocity of 12.25 feet per second. Accelerations were measured along the X-axis (roll) and Z-axis (yaw) by accelerometers rigidly installed at the center of gravity of the model. Actual flight path, contact attitudes, and motions were determined from high-speed motion pictures. Landings were made on concrete and on water. The peak accelerations along the X-axis for landings on concrete were in the order of 3Og for a 0 deg. contact attitude. A horizontal velocity of 7 feet per second, corresponding to a flight path of 60 deg., had very little effect upon the peak accelerations obtained for landings on concrete. For contact attitudes of -25 deg. and 30 deg. the peak accelerations along the Z-axis were about +/- l5g, respectively. The peak accelerations measured for the water landings were about one-third lower than the peak accelerations measured for the landings on concrete. Assuming a rigid body, computations were made by using Newton's second law of motion and the force-stroke characteristics of the air bag to determine accelerations for a flight path of 90 deg. (vertical) and a contact attitude of 0 deg. The computed and experimental peak accelerations and strokes at peak acceleration were in good agreement for the model. The special scaling appears to be applicable for predicting full-scale time and stroke at peak acceleration for a landing on concrete from a 90 deg. flight path at a 0 deg. It appears that the full-scale approximately the same as those obtained from the model for the range of attitudes and flight paths investigated.

McGehee, John R.↗
NASA NTRS2004

Detecting the Sun in 2-Gyr Old Open Clusters

The aim of this observation was to take advantage of the large collecting area and sensitivity of the European Space Agency's XMM-Newton telescope to survey the X ray brightness of solar-mass stars in the 2 Gyr old NGC 752 cluster in order to chart out the dependence of coronal X ray emission of Sun like stars as a function of age. Dr. Simon's stated role as one of the co-investigators in the project was limited to assisting the project's European Principal Investigator with the analysis and interpretation of the X ray measurements, and to assist in summarizing the results for publication. The observation was executed successfully, and approximately 130 X ray sources were detected in the XMM field of view and subsequently measured to obtain X ray count rates. To date neither the measurements nor the raw data have been distributed by the European P.I. to the co-investigators, It is expected that the measurements will be distributed to the co-investigators by the end of 2004, and once that is done, the scientific analysis by Dr. Simon can then begin. At present, since the data and measurements have not been shared with the co-investigators, there are no results to report.

Simon, Theodore↗
NASA NTRS2012

Radio and X-Ray Observations of SN 2006jd: Another Strongly Interacting Type IIn Supernova

We report four years of radio and X-ray monitoring of the Type IIn supernova SN 2006jd at radio wavelengths with the Very Large Array, Giant Metrewave Radio Telescope and Expanded Very Large Array at X-ray wavelengths with Chandra, XMM-Newton and Swift-XRT. We assume that the radio and X-ray emitting particles are produced by shock interaction with a dense circumstellar medium. The radio emission shows an initial rise that can be attributed to free-free absorption by cool gas mixed into the nonthermal emitting region external free-free absorption is disfavored because of the shape of the rising light curves and the low gas column density inferred along the line of sight to the emission region. The X-ray luminosity implies a preshock circumstellar density approximately 10(exp 6) per cubic meter at a radius r approximately 2 x 10(exp 16) centimeter, but the column density inferred from the photoabsorption of X-rays along the line of sight suggests a significantly lower density. The implication may be an asymmetry in the interaction. The X-ray spectrum shows Fe line emission at 6.9 keV that is stronger than is expected for the conditions in the X-ray emitting gas. We suggest that cool gas mixed into the hot gas plays a role in the line emission. Our radio and X-ray data both suggest the density profile is flatter than r2 because of the slow evolution of the unabsorbed emission.

Chandra, Poonam↗
NASA NTRS1990

Linear iterative solvers for implicit ODE methods

The numerical solution of stiff initial value problems, which lead to the problem of solving large systems of mildly nonlinear equations are considered. For many problems derived from engineering and science, a solution is possible only with methods derived from iterative linear equation solvers. A common approach to solving the nonlinear equations is to employ an approximate solution obtained from an explicit method. The error is examined to determine how it is distributed among the stiff and non-stiff components, which bears on the choice of an iterative method. The conclusion is that error is (roughly) uniformly distributed, a fact that suggests the Chebyshev method (and the accompanying Manteuffel adaptive parameter algorithm). This method is described, also commenting on Richardson's method and its advantages for large problems. Richardson's method and the Chebyshev method with the Mantueffel algorithm are applied to the solution of the nonlinear equations by Newton's method.

Saylor, Paul E.↗
DOE OSTI2025

A randomized sketching trust-region secant method for low-memory dynamic optimization

The numerical solution of dynamic optimization problems is often limited by the memory required to store the state trajectory, which is used to evaluate the objective function and its derivatives. Recently, [R. Muthukumar et al., SIAM Journal on Optimization 31(2), pp. 1242–1275 (2021)] introduced a trust-region method for dynamic optimization that employs randomized sketching to compress the state trajectory, resulting in inexact derivative computations. By adaptively learning the sketch rank, the trust-region algorithm achieves rigorous convergence guarantees. Here, we extend this approach to use secant Hessian approximations. Due to the randomness introduced by the sketch, the traditional secant update formulae can produce poor Hessian approximations. In particular, the difference of two gradients, computed from two different sketches, may be inconsistent. To overcome this, we employ a sketched approximation of the Hessian application, in lieu of computing the gradient difference. We numerically demonstrate the improved stability of this approach on an example from PDE-constrained optimization.

dynamic optimization↗
NASA NTRS1973

Neutral thermosphere temperatures from density scale height measurements.

During each 4-min satellite interrogation period the Explorer 32 density gages measured the atmospheric density approximately every 2 sec. Over certain segments of the satellite orbit these measurements determined the neutral atmospheric density scale height. The scale heights measured at an altitude of 400 (plus or minus 50) km have been analyzed to infer thermospheric temperatures. The results confirm an earlier conclusion from the density data of the same experiment that the diurnal temperature variation is latitude dependent.

Newton, G. P.↗
DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

GLAD-M35: a joint P and S global tomographic model with uncertainty quantification

We present our third and final generation joint P and S global adjoint tomography (GLAD) model, GLAD-M35, and quantify its uncertainty based on a low-rank approximation of the inverse Hessian. Starting from our second-generation model, GLAD-M25, we added 680 new earthquakes to the database for a total of 2160 events. New P-wave categories are included to compensate for the imbalance between P- and S-wave measurements, and we enhanced the window selection algorithm to include more major-arc phases, providing better constraints on the structure of the deep mantle and more than doubling the number of measurement windows to 40 million. Two stages of a Broyden–Fletcher–Goldfarb–Shanno (BFGS) quasi-Newton inversion were performed, each comprising five iterations. With this BFGS update history, we determine the model’s standard deviation and resolution length through randomized singular value decomposition.

58 GEOSCIENCES↗
NASA NTRS2001

FE K EMISSION AND ABSORPTION FEATURES IN THE XMM-EPIC SPECTRUM OF THE SEYFERT GALAXY IC 4329A

We present a re-analysis of the XMM-Newton long-look of the X-ray bright Seyfert galaxy IC 4329a. The Fe K bandpass is dominated by two peaks, consistent with emission from neutral or near-neutral Fe Ka and KP. A relativistic diskline model whereby both peaks are the result of one doubly-peaked diskline profile is found to be a poor description of the data. Models using two relativistic disklines are found to describe the emission profile well. A low-inclination, moderately-relativistic dual-diskline model is possible if the contribution from narrow components, due to distant material, is small or absent. A high-inclination, moderately relativistic profile for each peak is possible if there are roughly equal contributions from both the broad and narrow components. Upper limits on Fe XXV and Fe XXVI emission and absorption at the systemic velocity of IC 4329a are obtained. We also present the results of RXTE monitoring of this source obtained so far; the combined XMM-Newton and RXTE data sets allow us to explore the time-resolved spectral behavior of this source on time scales ranging from hours to 2 years. We find no strong evidence for variability of the Fe Ka emission line on any time scale probed, likely due to the minimal level of continuum variability. We detect a narrow absorption line, at a energy of 7.68 keV in the rest frame of the source; its significance has been confirmed using Monte Carlo simulations. This feature is most likely due to absorption from Fe XXVI blueshifted to approximately 0.1c relative to the systemic velocity, making IC 4329a the lowest-redshift AGN known with a high-velocity, highly-ionized outflow component. As is often the case with similar outflows seen in high-luminosity quasars, the estimated mass outflow rate is larger than the inflow accretion rate, signaling that the outflow represents a substantial portion of the total energy budget of the AGN. The outflow could arise from a radiatively-driven disk wind, or it may be in the form of a discrete, transient blob of ejected material.

Markowitz, A.↗
DOE OSTI2020

Deep XMM–Newton observations of the most distant SPT-SZ galaxy cluster

ABSTRACT We present results from a 577 ks XMM–Newton observation of SPT-CL J0459–4947, the most distant cluster detected in the South Pole Telescope 2500 square degree (SPT-SZ) survey, and currently the most distant cluster discovered through its Sunyaev–Zel’dovich effect. The data confirm the cluster’s high redshift, z = 1.71 ± 0.02, in agreement with earlier, less precise optical/IR photometric estimates. From the gas density profile, we estimate a characteristic mass of $M_{500}=(1.8\pm 0.2)\times 10^{14}\, {\rm M}_{\odot }$; cluster emission is detected above the background to a radius of $\sim \!2.2\, r_{500}$, or approximately the virial radius. The intracluster gas is characterized by an emission-weighted average temperature of 7.2 ± 0.3 keV and metallicity with respect to Solar of $Z/\, Z_{\odot }=0.37\pm 0.08$. For the first time at such high redshift, this deep data set provides a measurement of metallicity outside the cluster centre; at radii $r\gt 0.3\, r_{500}$, we find $Z/\, Z_{\odot }=0.33\pm 0.17$ in good agreement with precise measurements at similar radii in the most nearby clusters, supporting an early enrichment scenario in which the bulk of the cluster gas is enriched to a universal metallicity prior to cluster formation, with little to no evolution thereafter. The leverage provided by the high redshift of this cluster tightens by a factor of 2 constraints on evolving metallicity models, when combined with previous measurements at lower redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2005

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

The XMM-Newton observations of the GPS sample was completed last summer. We are in process of finalizing the paper describing the data and the results. The main goal of the project was to determine the X-ray spectra of the GPS galaxies in comparison to regular radio galaxies. Our XMM observations show evidence that the GPS galaxies are heavily obscured with the large absorbing columns exceeding N(H)greater than le22 cm^-2. Taking into account the obscuration we determined that the intrinsic X-ray luminosities of GPS galaxies are of order le43-le44 erg/s, comparable to low luminosity radio loud quasars. The large GPS samples can confirm the result, as at this moment our evidence is based only on 7 GPS galaxies observed with good S/N/ in X-rays. The first paper summarizing the results of the XMM observation of Mkn 668 has been published Astronomy & Astrophysics. We found soft X-ray signatures of a hot plasma (kT approximately 10^7~K) and a hard X-ray emission from the nucleus. The X-ray spectrum above 2.5-keV is characterized by a very flat (observed photon index, Gamma-0.5) power-law continuum, alongside with a strong Fe-K-alpha neutral iron fluorescent line (EW~600-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 ~ to 10^24 ~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 source represented the second hard X-ray detection of the GPS galaxy ever (the first one being 1345+125; O'Dea et al. 2000). The observations of the other GPS galaxies in our sample confirmed the trend of the large obscuration present in the spectra. However, we do not have a compelling evidence for a hot gas in the nucleus. The two other GPS galaxies observed with Chandra were added to the total of 7 GPS galaxies. This GPS sample was compared against the sample of radio loud galaxies and quasar. The paper is about to be submitted. Our preliminarily results based on the XMM-Newton data for the entire sample were presented at the meeting in July 2004.

Mushotsky, Richard↗
DOE OSTI2021

Lorentz symmetry in ghost-free massive gravity

The role of Lorentz symmetry in ghost-free massive gravity is studied, emphasizing features emerging in approximately Minkowski spacetime. The static extrema and saddle points of the potential are determined and their Lorentz properties identified. Solutions preserving Lorentz invariance and ones breaking four of the six Lorentz generators are constructed. Locally, globally, and absolutely stable Lorentz-invariant extrema are found to exist for certain parameter ranges of the potential. Gravitational waves in the linearized theory are investigated. Deviations of the fiducial metric from the Minkowski metric are shown to lead to pentarefringence of the five wave polarizations, which can include superluminal modes and subluminal modes with negative energies in certain observer frames. The Newton limit of ghost-free massive gravity is explored. The propagator is constructed and used to obtain the gravitational potential energy between two point masses. Here, the result extends the Fierz-Pauli limit to include corrections generically breaking both rotation and boost invariance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

A simple suboptimal least-squares algorithm for attitude determination with multiple sensors

Three-axis attitude determination is equivalent to finding a coordinate transformation matrix which transforms a set of reference vectors fixed in inertial space to a set of measurement vectors fixed in the spacecraft. The attitude determination problem can be expressed as a constrained optimization problem. The constraint is that a coordinate transformation matrix must be proper, real, and orthogonal. A transformation matrix can be thought of as optimal in the least-squares sense if it maps the measurement vectors to the reference vectors with minimal 2-norm errors and meets the above constraint. This constrained optimization problem is known as Wahba's problem. Several algorithms which solve Wahba's problem exactly have been developed and used. These algorithms, while steadily improving, are all rather complicated. Furthermore, they involve such numerically unstable or sensitive operations as matrix determinant, matrix adjoint, and Newton-Raphson iterations. This paper describes an algorithm which minimizes Wahba's loss function, but without the constraint. When the constraint is ignored, the problem can be solved by a straightforward, numerically stable least-squares algorithm such as QR decomposition. Even though the algorithm does not explicitly take the constraint into account, it still yields a nearly orthogonal matrix for most practical cases; orthogonality only becomes corrupted when the sensor measurements are very noisy, on the same order of magnitude as the attitude rotations. The algorithm can be simplified if the attitude rotations are small enough so that the approximation sin(theta) approximately equals theta holds. We then compare the computational requirements for several well-known algorithms. For the general large-angle case, the QR least-squares algorithm is competitive with all other know algorithms and faster than most. If attitude rotations are small, the least-squares algorithm can be modified to run faster, and this modified algorithm is faster than all but a similarly specialized version of the QUEST algorithm. We also introduce a novel measurement averaging technique which reduces the n-measurement case to the two measurement case for our particular application, a star tracker and earth sensor mounted on an earth-pointed geosynchronous communications satellite. Using this technique, many n-measurement problems reduce to less than or equal to 3 measurements; this reduces the amount of required calculation without significant degradation in accuracy. Finally, we present the results of some tests which compare the least-squares algorithm with the QUEST and FOAM algorithms in the two-measurement case. For our example case, all three algorithms performed with similar accuracy.

Brozenec, Thomas F.↗
NASA NTRS1985

Minimal norm constrained interpolation

In computational fluid dynamics and in CAD/CAM, a physical boundary is usually known only discreetly and most often must be approximated. An acceptable approximation preserves the salient features of the data such as convexity and concavity. In this dissertation, a smooth interpolant which is locally concave where the data are concave and is locally convex where the data are convex is described. The interpolant is found by posing and solving a minimization problem whose solution is a piecewise cubic polynomial. The problem is solved indirectly by using the Peano Kernal theorem to recast it into an equivalent minimization problem having the second derivative of the interpolant as the solution. This approach leads to the solution of a nonlinear system of equations. It is shown that Newton's method is an exceptionally attractive and efficient method for solving the nonlinear system of equations. Examples of shape-preserving interpolants, as well as convergence results obtained by using Newton's method are also shown. A FORTRAN program to compute these interpolants is listed. The problem of computing the interpolant of minimal norm from a convex cone in a normal dual space is also discussed. An extension of de Boor's work on minimal norm unconstrained interpolation is presented.

Irvine, L. D.↗
NASA NTRS1990

On the numerical solution of the dynamically loaded hydrodynamic lubrication of the point contact problem

The transient analysis of hydrodynamic lubrication of a point-contact is presented. A body-fitted coordinate system is introduced to transform the physical domain to a rectangular computational domain, enabling the use of the Newton-Raphson method for determining pressures and locating the cavitation boundary, where the Reynolds boundary condition is specified. In order to obtain the transient solution, an explicit Euler method is used to effect a time march. The transient dynamic load is a sinusoidal function of time with frequency, fractional loading, and mean load as parameters. Results include the variation of the minimum film thickness and phase-lag with time as functions of excitation frequency. The results are compared with the analytic solution to the transient step bearing problem with the same dynamic loading function. The similarities of the results suggest an approximate model of the point contact minimum film thickness solution.

Lim, Sang G.↗
NASA NTRS2010

The Atacama Cosmology Telescope: Physical Properties and Purity of a Galaxy Cluster Sample Selected Via the Sunyaev-Zel'Dovich Effect

We present optical and X-ray properties for the first confirmed galaxy cluster sample selected by the Sunyaev-Zel'dovich Effect from 148 GHz maps over 455 square degrees of sky made with the Atacama Cosmology Telescope. These maps. coupled with multi-band imaging on 4-meter-class optical telescopes, have yielded a sample of 23 galaxy clusters with redshifts between 0.118 and 1.066. Of these 23 clusters, 10 are newly discovered. The selection of this sample is approximately mass limited and essentially independent of redshift. We provide optical positions, images, redshifts and X-ray fluxes and luminosities for the full sample, and X-ray temperatures of an important subset. The mass limit of the full sample is around 8.0 x 10(exp 14) Stellar Mass. with a number distribution that peaks around a redshift of 0.4. For the 10 highest significance SZE-selected cluster candidates, all of which are optically confirmed, the mass threshold is 1 x 10(exp 15) Stellar Mass and the redshift range is 0.167 to 1.066. Archival observations from Chandra, XMM-Newton. and ROSAT provide X-ray luminosities and temperatures that are broadly consistent with this mass threshold. Our optical follow-up procedure also allowed us to assess the purity of the ACT cluster sample. Eighty (one hundred) percent of the 148 GHz candidates with signal-to-noise ratios greater than 5.1 (5.7) are confirmed as massive clusters. The reported sample represents one of the largest SZE-selected sample of massive clusters over all redshifts within a cosmologically-significant survey volume, which will enable cosmological studies as well as future studies on the evolution, morphology, and stellar populations in the most massive clusters in the Universe.

Menanteau, Felipe↗
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗