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 217 records · Page 12

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.↗

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.↗

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.↗

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.↗

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↗

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↗

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.↗

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↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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↗

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↗

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.↗

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.↗