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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
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↗
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.↗
NASA NTRS1994

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗
NASA NTRS1986

A numerical algorithm for optimal feedback gains in high dimensional LQR problems

A hybrid method for computing the feedback gains in linear quadratic regulator problems is proposed. The method, which combines the use of a Chandrasekhar type system with an iteration of the Newton-Kleinman form with variable acceleration parameter Smith schemes, is formulated so as 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 advantage of the proposed algorithm over the standard eigenvector (Potter, Laub-Schur) based techniques are discussed and numerical evidence of the efficacy of our ideas presented.

Banks, H. T.↗
NASA NTRS1974

Propagation of sound through a real jet flowfield

The sound propagation of harmonic disturbances through a real jet flowfield, which contains mean flow gradients in pressure and velocities, has been studied. A finite-difference approximation for the equations which govern the acoustic disturbance is obtained for a subsonic axisymmetric jet in terms of acoustic pressure and acoustic energy flux. Experimental results are used to describe the mean flowfield of a model jet. The equations are solved numerically using a Newton-type iterative scheme. The directivity at various distances from a point source located at two jet diameters downstream of the jet exit on its center line compared qualitatively well with the numerical result of Schubert (1972) and with the experimental result of Grande (1966).

Liu, C. H.↗
NASA NTRS1977

Optimization of flexible wing structures subject to strength and induced drag constraints

An optimization procedure for designing wing structures subject to stress, strain, and drag constraints is presented. The optimization method utilizes an extended penalty function formulation for converting the constrained problem into a series of unconstrained ones. Newton's method is used to solve the unconstrained problems. An iterative analysis procedure is used to obtain the displacements of the wing structure including the effects of load redistribution due to the flexibility of the structure. The induced drag is calculated from the lift distribution. Approximate expressions for the constraints used during major portions of the optimization process enhance the efficiency of the procedure. A typical fighter wing is used to demonstrate the procedure. Aluminum and composite material designs are obtained. The tradeoff between weight savings and drag reduction is investigated.

Haftka, R. T.↗
NASA NTRS1991

Numerical recovery of material parameters in Euler-Bernoulli beam models

A fully Sinc-Galerkin method for recovering the spatially varying stiffness parameter in fourth-order time-dependence problems with fixed and cantilever boundary conditions is presented. The forward problems are discretized with a sinc basis in both the spatial and temporal domains. This yields an approximation solution which converges exponentially and is valid on the infinite time interval. When the forward methods are applied to parameter recovery problems, the resulting inverse problems are ill-posed. Tikhonov regularization is applied and the resulting minimization problems are solved via a quasi-Newton/trust region algorithm. The L-curve method is used to determine an appropriate value of the regularization parameter. Numerical results which highlight the method are given for problems with both fixed and cantilever boundary conditions.

Smith, R. C.↗
NASA NTRS2007

X-Ray, UV, and Optical Observations of Supernova 2006bp with Swift: Detection of Early X-Ray Emission

We present results on the X-ray and optical/UV emission from the Type IIP supernova (SN) 2006bp and the interaction of the SW shock with its environment, obtained with the X-Ray Telescope (XRT) and UV/Optical Telescope (UVOT) on-board the Swift observatory. SN 2006bp is detected in X-rays at a 4.5 sigmalevel of significance in the merged XRT data from days 1 to 12 after the explosion. If the (0.2-10 keV band) X-ray luminosity of L(sub 0.2-10) = (1.8 plus or minus 0.4) x l0(exp 39 ergs s(exp -1) is caused by interaction of the SN shock with circumstellar material (CSM), deposited by a stellar wind from the progenitor's companion star, a mass-loss rate of M is approximately 2x10(exp -6) solar mass yr(exp -1) (v(sub w)/10 km s(exp -l) is inferred. The mass-loss rate is one of the lowest ever recorded for a core-collapse SN and consistent with the non-detection in the radio with the VLA on days 2, 9, and 11 after the explosion. The Swift data further show a fading of the X-ray emission starting around day 12 after the explosion. In combination with a follow-up XMM-Newton observation obtained on day 21 after the explosion, an X-ray rate of decline Lx, varies as t(exp -n) with index n = 1.2 plus or minus 0.6 is inferred. Since no other SN has been detected in X-rays prior to the optical peak and since Type IIP SNe have an extended 'plateau' phase in the optical, we discuss the scenario that the X-rays might be due to inverse Compton scattering of photospheric optical photons off relativistic electrons produced in circumstellar shocks. However, due to the high required value of the Lorentz factor (approximately 10-100), inconsistent with the ejecta velocity inferred from optical line widths, we conclude that Inverse Compton scattering is an unlikely explanation for the observed X-ray emission. The fast evolution of the optical/ultraviolet (1900-5500A) spectral energy distribution and the spectral changes observed with Swift reveal the onset of metal line-blanketing and cooling of the expanding photosphere during the first few weeks after the outburst.

Immler, S.↗
NASA NTRS2006

Highly Ionized Iron Absorption Lines from Outflowing Gases in the X-ray Spectrum of NGC 1365

We present the discovery of four absorption lines in the X-ray spectrum of the Seyfert galaxy NGC 1365, at energies between 6.7 and 8.3 keV. The lines are detected with high statistical confidence (from >20 sigma for the strongest to -4 sigma for the weakest) in two XMM-Newton observations 60 ks long. We also detect the same lines, with a lower signal-to-noise ratio (but still >2 sigma for each line), in two previous shorter (-10 ks) XMM-Newton observations. The spectral analysis identifies these features as Fe XXV and Fe XXVI Kalpha and Kbeta lines, outflowing with velocities varying between -1000 and -5000 km/s among the observations. These are the highest quality detections of such lines so far. The high equivalent widths [EW (Kalpha) approximately 100 eV] and the Kalpha/Kbeta ratios imply that the lines are due to absorption of the AGN continuum by a highly ionized gas with column density NH-5?1023 cm(exp -2) at a distance of -(50-100)RS from the continuum source.

Risaliti, G.↗
NASA NTRS2012

The Outer Limits of Galaxy Clusters: Observations to the Virial Radius with Suzaku, XMM,and Chandra

The outskirts of galaxy clusters, near the virial radius, remain relatively unexplored territory and yet are vital to our understanding of cluster growth, structure, and mass. In this presentation, we show the first results from a program to constrain the sate of the outer intra-cluster medium (ICM) in a large sample of galaxy clusters, exploiting the strengths of three complementary X-ray observatories: Suzaku (low, stable background), XMM-Newton (high sensitivity),and Chandra (good spatial resolution). By carefully combining observations from the cluster core to beyond r200, we are able to identify and reduce systematic uncertainties that would impede our spatial and spectral analysis using a single telescope. Our sample comprises nine clusters at z is approximately 0.1-0.2 fully covered in azimuth to beyond r200, and our analysis indicates that the ICM is not in hydrostatic equilibrium in the cluster outskirts, where we see clear azimuthal variations in temperature and surface brightness. In one of the clusters, we are able to measure the diffuse X-ray emission well beyond r200, and we find that the entropy profile and the gas fraction are consistent with expectations from theory and numerical simulations. These results stand in contrast to recent studies which point to gas clumping in the outskirts; the extent to which differences of cluster environment or instrumental effects factor in this difference remains unclear. From a broader perspective, this project will produce a sizeable fiducial data set for detailed comparison with high-resolution numerical simulations.

Miller, Eric D.↗
NASA NTRS2013

X-ray View of Four High-Luminosity Swift-BAT AGN: Unveiling Obscuration and Reflection with Suzaku

Aims. A complete census of obscured Active Galactic Nuclei (AGN) is necessary to reveal the history of the super massive black hole (SMBH) growth and galaxy evolution in the Universe given the complex feedback processes and the fact that much of this growth occurs in an obscured phase. In this context, hard X-ray surveys and dedicated follow-up observations represent a unique tool for selecting highly absorbed AGN and for characterizing the obscuring matter surrounding the SMBH. Here we focus on the absorption and reflection occurring in highly luminous, quasar-like AGN, to study the relation between the geometry of the absorbing matter and the AGN nature (e.g. X-ray, optical, and radio properties), and to help to determine the column density dependency on the AGN luminosity. Methods. The Swift/BAT nine-month survey observed 153 AGN, all with ultra-hard X-ray BAT fluxes in excess of 10(exp −11) erg per square centimeter and an average redshift of 0.03. Among them, four of the most luminous BAT AGN (44.73 less than LogLBAT less than 45.31) were selected as targets of Suzaku follow-up observations: J2246.0+3941 (3C 452), J0407.4+0339 (3C 105), J0318.7+6828, and J0918.5+0425. The column density, scattered/reflected emission, the properties of the Fe K line, and a possible variability are fully analyzed. For the latter, the spectral properties from Chandra, XMM-Newton and Swift/XRT public observations were compared with the present Suzaku analysis, adding an original spectral analysis when non was available from the literature. Results. Of our sample, 3C 452 is the only certain Compton-thick AGN candidate because of i) the high absorption (N(sub H) approximately 4 × 10(exp 23) per square centimeter) and strong Compton reflection; ii) the lack of variability; iii) the "buried" nature, i.e. the low scattering fraction (less than 0.5%) and the extremely low relative [OIII] luminosity. In contrast 3C 105 is not reflection-dominated, despite the comparable column density, X-ray luminosity and radio morphology, but shows a strong long-term variability in flux and scattering fraction, consistent with the soft emission being scattered from a distant region (e.g., the narrow emission line region). The sample presents high (greater than 100) X-to- [OIII] luminosity ratios, with an extreme value of R(sup X)(sub [OIII]) approximately 800 for 3C 452, confirming the [OIII] luminosity to be affected by residual extinction in presence of mild absorption, especially for "buried" AGN such as 3C 452. Three of our targets are powerful FRII radio galaxies, which is shown by their high luminosity and absorption; this makes them the most luminous and absorbed AGN of the BAT Seyfert survey despite the inversely proportional N(sub H) - L(sub X) relation.

SwiftBAG AGN↗
NASA NTRS1985

Structural parameter identification of distributed systems using finite element approximation

A system identification technique is developed for classes of distributed systems using finite element approximations. Vibrating systems represented by partial differential equations have physical parameters associated with mass, stiffness, and damping distributions which need to be known in order to properly control and design mathematical models of the system. In order to identify these parameters a weighted least-squares algorithm and modified Newton-Raphson method is used for the identification process. The theory and technique is demonstrated by estimating the system parameters of a vibrating cantilever beam made up of several different structural properties.

Lee, K. Y.↗
NASA NTRS1989

Gradiometry and gravitomagnetic field detection

Gravitomagnetism was apparently first introduced into physics about 120 years ago when major developments in electrodynamics and the strong similarity between Coulomb's law of electricity and Newton's law of gravity led to the hypothesis that mass current generates a fundamental force of gravitational origin analogous to the magnetic force caused by charge current. According to general relativity, the rotation of a body leads to the dragging of the local inertial frames. In the weak-field approximation, the dragging frequency can be interpreted, up to a constant proportionality factor, as a gravitational magnetic field. There is, as yet, no direct evidence regarding the existence of such a field. The possibility is examined of detecting the gravitomagnetic field of the Earth by gravity gradiometry.

Mashhoon, Bahram↗
NASA NTRS1995

A Hybrid Parallel Preconditioning Algorithm For CFD

A new hybrid preconditioning algorithm will be presented which combines the favorable attributes of incomplete lower-upper (ILU) factorization with the favorable attributes of the approximate inverse method recently advocated by numerous researchers. The quality of the preconditioner is adjustable and can be increased at the cost of additional computation while at the same time the storage required is roughly constant and approximately equal to the storage required for the original matrix. In addition, the preconditioning algorithm suggests an efficient and natural parallel implementation with reduced communication. Sample calculations will be presented for the numerical solution of multi-dimensional advection-diffusion equations. The matrix solver has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The full paper will show numerous examples in CFD to demonstrate the efficiency and robustness of the method.

Barth,Timothy J.↗
NASA NTRS2012

Investigating SLIM Disk Solutions FOR HLX-1 IN ESO 243-49

The hyperluminous X-ray source HLX-1 in the galaxy ESO 243-49, currently the best intermediate-mass blackhole (BH) candidate, displays spectral transitions similar to those observed in Galactic BH binaries, but with aluminosity 100-1000 times higher. We investigated the X-ray properties of this unique source by fitting multiepochdata collected by Swift, XMM-Newton, and Chandra with a disk model computing spectra for a wide rangeof sub- and super-Eddington accretion rates assuming a non-spinning BH and a face-on disk (i=0 deg.). Under theseassumptions we find that the BH in HLX-1 is in the intermediate-mass range (approximately 2 x 10(exp 4) solar mass) and the accretionflow is in the sub-Eddington regime. The disk radiation efficiency is eta = 0.11 plus or minus 0.03. We also show that the source does follow the LX is proportional to T(exp 4) relation for our mass estimate. At the outburst peaks, the source radiates near the Eddington limit. The accretion rate then stays constant around 4 x 10(exp 4) solar mass yr (sup -1) for several days and then decreases exponentially. Such plateaus in the accretion rate could be evidence that enhanced mass-transfer rateis the driving outburst mechanism in HLX-1. We also report on the new outburst observed in 2011 August by theSwift X-Ray Telescope. The time of this new outburst further strengthens the approximately 1 year recurrence timescale.

solutions↗
NASA NTRS2006

A Possible Magnetar Nature for IGR J16358-4726

We present detailed spectral and timing analysis of the hard x-ray transient IGR J16358-4726 using multi-satellite archival observations. A study of the source flux time history over 6 years, suggests that this transient outbursts can be occurring in intervals of at most 1 year. Joint spectral fits using simultaneous Chandra/ACIS and INTEGRAL/ISGRI data reveal a spectrum well described by an absorbed cut-off power law model plus an Fe line. We detected the pulsations initially reported using Chandra/ACIS also in the INTEGRAL/ISGRI light curve and in subsequent XMM-Newton observations. Using the INTEGRAL data we identified a pulse spin up of 94 s (P = 1.6 x 10(exp -4), which strongly points to a neutron star nature for IGR J16358-4726. Assuming that the spin up is due to disc accretion, we estimate that the source magnetic field ranges between 10(sup 13) approximately 10(sup 15) depending on its distance, possibly supporting a magnetar nature for IGR J16358-4726.

Patel, S.↗
NASA NTRS2023

A Finite Element Approach for Simplified 2D Nonlinear Dynamic Contact/Impact Analysis

In this paper, a simplified numerical approach for finite element dynamic analysis of an inelastic solid structure subjected to solid object impact is presented. The approach approximates the impacting solid as the selected multiple nodes, for which mass of the impactor is distributed. The node-to-segment contact formulation with the penalty constraint technique incorporated is employed to impose contact conditions between the nodes and the surface of the receiver structure. The node-to-segment algorithm is integrated into Newton–Raphson time integration scheme and the Lagrange multiplier technique is applied to enforce the identical displacements for the selected nodes throughout the analysis process. The approach is verified using two-dimensional plane strain models considering elastic-perfectly-plastic material behavior. The results obtained using the proposed approach are in a good agreement with those simulated using a commercial finite element code, ABAQUS dynamic/implicit, in terms of displacements and stress distribution fields. The proposed approach is shown to be computationally superior to general finite element method-based contact/impact analysis without significantly sacrificing the accuracy.

Seungwook Seok↗