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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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.↗
DOE OSTI2020

Gravitational Fluctuations as an Alternative to Inflation III. Numerical Results

Power spectra play an important role in the theory of inflation, and their ability to reproduce current observational data to high accuracy is often considered a triumph of inflation, largely because of a lack of credible alternatives. In previous work we introduced an alternative picture for the cosmological power spectra based on the nonperturbative features of the quantum version of Einstein’s gravity, instead of currently popular inflation models based on scalar fields. The key ingredients in this new picture are the appearance of a nontrivial gravitational vacuum condensate (directly related to the observed cosmological constant), and a calculable renormalization group running of Newton’s G on cosmological scales. More importantly, one notes the absence of any fundamental scalar fields in this approach. Results obtained previously were largely based on a semi-analytical treatment, and thus, while generally transparent in their implementation, often suffered from the limitations of various approximations and simplifying assumptions. In this work, we extend and refine our previous calculations by laying out an updated and extended analysis, which now utilizes a set of suitably modified state-of-the-art numerical programs (ISiTGR, MGCAMB and MGCLASS) developed for observational cosmology. As a result, we are able to remove some of the approximations employed in our previous studies, leading to a number of novel and detailed physical predictions. These should help in potentially distinguishing the vacuum condensate picture of quantum gravity from that of other models such as scalar field inflation. Here, besides the matter power spectrum P m (k) , we work out, in detail, predictions for what are referred to as the TT, TE, EE, BB angular spectra, as well as their closely related lensing spectra. However, the current limited precision of observational data today (especially on large angular scales) does not allow us yet to clearly prove or disprove either set of ideas. Nevertheless, by exploring in more details the relationship between gravity and cosmological matter and radiation both analytically and numerically, together with an expected future influx of increasingly accurate observational data, one can hope that the new quantum gravitational picture can be subjected to further stringent tests in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2025

First‐Order Empirical Interpolation Method for Real‐Time Solution of Parametric Time‐Dependent Nonlinear PDEs

ABSTRACT We present a model reduction approach for the real‐time solution of time‐dependent nonlinear partial differential equations (PDEs) with parametric dependencies. A major challenge in constructing efficient and accurate reduced‐order models for nonlinear PDEs is the efficient treatment of nonlinear terms. We address this by unifying the implementation of hyperreduction methods to deal with nonlinear terms. Furthermore, we introduce a first‐order empirical interpolation method (EIM) to provide an efficient approximation of the nonlinear terms in time‐dependent PDEs. We demonstrate the effectiveness of our approach on the Allen–Cahn equation, which models phase separation, and the Buckley–Leverett equation, which describes two‐phase fluid flow in porous media. Numerical results highlight the accuracy, efficiency, and stability of the proposed method compared with both the Galerkin–Newton approach and hyper‐reduced models using the standard EIM.

Nguyen, Ngoc Cuong [Center for Computational Engin↗
NASA NTRS2012

Black Holes in Bulgeless Galaxies: An XMM-Newton Investigation of NGC 3367 AND NGC 4536

The vast majority of optically identified active galactic nuclei (AGNs) in the local Universe reside in host galaxies with prominent bulges, supporting the hypothesis that black hole formation and growth is fundamentally connected to the build-up of galaxy bulges. However, recent mid-infrared spectroscopic studies with Spitzer of a sample of optically "normal" late-type galaxies reveal remarkably the presence of high-ionization [NeV] lines in several sources, providing strong evidence for AGNs in these galaxies. We present follow-up X-ray observations recently obtained with XMM-Newton of two such sources, the late-type optically normal galaxies NGC 3367 and NGC 4536. Both sources are detected in our observations. Detailed spectral analysis reveals that for both galaxies, the 2-10 keV emission is dominated by a power law with an X-ray luminosity in the L(sub 2- 10 keV) approximates 10(exp 39) - 10(exp 40) ergs/s range, consistent with low luminosity AGNs. While there is a possibility that X-ray binaries account for some fraction of the observed X-ray luminosity, we argue that this fraction is negligible. These observations therefore add to the growing evidence that the fraction of late-type galaxies hosting AGNs is significantly underestimated using optical observations alone. A comparison of the midinfrared [NeV] luminosity and the X-ray luminosities suggests the presence of an additional highly absorbed X-ray source in both galaxies, and that the black hole masses are in the range of 10(exp 5) - 10(exp 7) solar M for NGC 3367 and 10(exp 4) - (exp 10) solar M for NGC 4536

McAlpine, W.↗