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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2012

Iterative Atmospheric Correction Scheme and the Polarization Color of Alpine Snow

Characterization of the Earth's surface is crucial to remote sensing, both to map geomorphological features and because subtracting this signal is essential during retrievals of the atmospheric constituents located between the surface and the sensor. Current operational algorithms model the surface total reflectance through a weighted linear combination of a few geometry-dependent kernels, each devised to describe a particular scattering mechanism. The information content of these measurements is overwhelmed by that of instruments with polarization capabilities: proposed models in this case are based on the Fresnel reflectance of an isotropic distribution of facets. Because of its remarkable lack of spectral contrast, the polarized reflectance of land surfaces in the shortwave infrared spectral region, where atmospheric scattering is minimal, can be used to model the surface also at shorter wavelengths, where aerosol retrievals are attempted based on well-established scattering theories. In radiative transfer simulations, straightforward separation of the surface and atmospheric contributions is not possible without approximations because of the coupling introduced by multiple reflections. Within a general inversion framework, the problem can be eliminated by linearizing the radiative transfer calculation, and making the Jacobian (i.e., the derivative expressing the sensitivity of the reflectance with respect to model parameters) available at output. We present a general methodology based on a Gauss-Newton iterative search, which automates this procedure and eliminates de facto the need of an ad hoc atmospheric correction. In this case study we analyze the color variations in the polarized reflectance measured by the NASA Goddard Institute of Space Studies Research Scanning Polarimeter during a survey of late-season snowfields in the High Sierra. This insofar unique dataset presents challenges linked to the rugged topography associated with the alpine environment and a likely high water content due to melting. The analysis benefits from ancillary information provided by the NASA Langley High Spectral Resolution Lidar deployed on the same aircraft. The results obtained from the iterative scheme are contrasted against the surface polarized reflectance obtained ignoring multiple reflections, via the simplistic subtraction of the atmospheric scattering contribution. Finally, the retrieved reflectance is modeled after the scattering properties of a dense collection of ice crystals at the surface. Confirming that the polarized reflectance of snow is spectrally flat would allow to extend the techniques already in use for polarimetric retrievals of aerosol properties over land to the large portion of snow-covered pixels plaguing orbital and suborbital observations.

Polarized BRDF↗

NASA NTRSDate not supplied

Mrk 421: 2 Weeks of Multi-Wavelength Monitoring with IXPE

At z=0.03, Mrk 421 is one of the closest Blazars and hence an exceptional target for X-ray polarization. In December of 2023, IXPE monitored Mrk 421 for approximately two weeks, in order to investigate the possible origins of X-ray polarization variability previously seen in this target, and whether it might originate in a shock propagating through the field lines of an energy-stratified jet. We present the results of the IXPE campaign, which is also complemented by NuSTAR, XMM-Newton, Swift, and ground-based monitoring.

IXPE↗

NASA NTRSDate not supplied

Mrk 421: Two Weeks of Multi-Wavelength Monitoring with IXPE

IXPE↗

DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2007

Quasi-Periodic Variability in NGC 5408 X-1

We report the discovery with XMM-Newton of quasiperiodic variability in the 0.2 - 10 keV X-ray flux from the ultraluminous X-ray source NGC 5408 X-1. The average power spectrum of all EPIC-pn data reveals a strong 20 mHz QPO with an average amplitude (rms) of 9%, and a coherence, Q identical with nu(sub 0)/sigma approximately equal to 6. In a 33 ksec time interval when the 20 mHz QPO is strongest we also find evidence for a 2nd QPO peak at 15 mHz, the first indication for a close pair of QPOs in a ULX source. Interestingly, the frequency ratio of this QPO pair is inconsistent with 3:2 at the 3 sigma level, but is consistent with a 4:3 ratio. A powerlaw noise component with slope near 1.5 is also present below 0.1 Hz with evidence for a break to a flatter slope at about 3 mHz. The source shows substantial broadband variability, with a total amplitude (rms) of about 30% in the 0.1 - 100 mHz frequency band, and there is strong energy dependence to the variability. The power spectrum of hard X-ray photons (greater than 2 keV) shows a "classic" flat-topped continuum breaking to a power law with index 1.5 - 2. Both the break and 20 mHz QPO are detected in the hard band, and the 20 mHz QPO is essentially at the break. The QPO is both strong and narrow in this band, having an amplitude (rms) of 15%, and Q approx. equal to 25. The energy spectrum is well fit by three components, a "cool" disk with kT = 0.15 keV, a steep power law with index 2.56, and a thermal plasma at kT = 0.87 keV. The disk, power law, and thermal plasma components contribute 35, 60, and 5% of the 0.3 - 10 keV flux, respectively. Both the timing and spectral properties of NGC 5408 X-1 are strikingly reminiscent of Galactic black hole systems at high inferred accretion rates, but with its characteristic frequencies (QPO and break frequencies) scaled down by a factor of 10 - 100. We discuss the implications of these findings in the context of models for ULXs, and their implications for the object's mass.

Strohmayer, Tod E.↗

NASA NTRS1978

Efficiency trade-offs of steady-state methods using FEM and FDM

The efficiency characteristics of finite element and finite difference approximations for the steady-state solution of the Navier-Stokes equations are examined. The finite element method discussed is a standard Galerkin formulation of the incompressible, steady-state Navier-Stokes equations. The finite difference formulation uses simple centered differences that are O(delta x-squared). Operation counts indicate that a rapidly converging Newton-Raphson-Kantorovitch iteration scheme is generally preferable over a Picard method. A split NOS Picard iterative algorithm for the finite difference method was most efficient.

Gartling, D. K.↗

NASA NTRS2003

XMM-Newton Observation of Fe K(alpha) Emission from a BAL QSO: Mrk 231

We present results from a 20 ksec XMM-Newton observation of Mrk 231. EPIC spectral data reveal strong line emission due to Fe K alpha, which has rarely been detected in this class, as BAL QSOs are very faint in the X-ray band. The line energy is consistent with an origin in neutral Fe. The width of the line is equivalent to a velocity dispersion approximately 18,000 kilometers per second and thus the line may be attributed to transmission and/or reflection from a distribution of emitting clouds. If, instead, the line originates in the accretion disk then the line strength and flat X-ray continuum support some contribution from a reflected component, although the data disfavor a model where the hard X-ray band is purely reflected X-rays. The line parameters are similar to those obtained for the Fe Ka line detected in another BAL QSO, H1413 + 117.

Turner, T. J.↗

NASA NTRS2003

Elemental Abundances in NGC 3516

We present Reflection Grating Spectrometer data from an XMM-Newton observation of the Seyfert 1 galaxy NGC 3516, taken while the continuum source was in an extremely low flux state. This observation offers a rare opportunity for a detailed study of emission from a Seyfert 1 galaxy as these are usually dominated by high nuclear continuum levels and heavy absorption. The spectrum shows numerous narrow emission lines (FWHM approximately less than 1300 kilometers per second) in the 0.3 - 2 keV range, including the H-like lines of C, N, and O and the He-like lines of N, O and Ne. The emission-line ratios and the narrow width of the radiative recombination continuum of CVI indicate that the gas is photoionized and of fairly low temperature (kT approximately less than 0.01 keV). The availability of emission lines from different elements for two iso-electronic sequences allows us to constrain the element abundances. These data show that the N lines are far stronger than would be expected from gas of solar abundances. Based on our photoionization models we find that nitrogen is overabundant in the central regions of the galaxy, compared to carbon, oxygen and neon by at least a factor of 2.5. We suggest that this is the result of secondary production of nitrogen in intermediate mass stars, and indicative of the history of star formation in NGC 3516.

Turner, T. J.↗

NASA NTRS2004

Textbook Multigrid Efficiency for Leading Edge Stagnation

A multigrid solver is defined as having textbook multigrid efficiency (TME) if the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in evaluating the discrete residuals. TME in solving the incompressible inviscid fluid equations is demonstrated for leading- edge stagnation flows. The contributions of this paper include (1) a special formulation of the boundary conditions near stagnation allowing convergence of the Newton iterations on coarse grids, (2) the boundary relaxation technique to facilitate relaxation and residual restriction near the boundaries, (3) a modified relaxation scheme to prevent initial error amplification, and (4) new general analysis techniques for multigrid solvers. Convergence of algebraic errors below the level of discretization errors is attained by a full multigrid (FMG) solver with one full approximation scheme (F.4S) cycle per grid. Asymptotic convergence rates of the F.4S cycles for the full system of flow equations are very fast, approaching those for scalar elliptic equations.

Diskin, Boris↗

NASA NTRS2004

Textbook Multigrid Efficiency for Leading Edge Stagnation

A multigrid solver is defined as having textbook multigrid efficiency (TME) if the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in evaluating the discrete residuals. TME in solving the incompressible inviscid fluid equations is demonstrated for leading-edge stagnation flows. The contributions of this paper include (1) a special formulation of the boundary conditions near stagnation allowing convergence of the Newton iterations on coarse grids, (2) the boundary relaxation technique to facilitate relaxation and residual restriction near the boundaries, (3) a modified relaxation scheme to prevent initial error amplification, and (4) new general analysis techniques for multigrid solvers. Convergence of algebraic errors below the level of discretization errors is attained by a full multigrid (FMG) solver with one full approximation scheme (FAS) cycle per grid. Asymptotic convergence rates of the FAS cycles for the full system of flow equations are very fast, approaching those for scalar elliptic equations.

Diskin, Boris↗

NASA NTRS1999

Atomic Data Needs for X-Ray Astronomy

This publication contains written versions of most of the invited talks presented at the workshop on Atomic Data Needs for X-ray Astronomy which was held at NASA's Goddard Space Flight Center on December 16-1 7 1999. The idea of hosting such a workshop emerged from an imminent need to update and complete current atomic datasets in anticipation of a new era of high quality X-ray spectra starting with the launching of Chandra and XMM-Newton observatories. At first, our vision of the workshop was of a short and limited attendance event, given the specialization of the topic. But it, was soon realized, from the response to the first workshop announcement, that the topic was of much interest, to researchers working in X-ray spectra (physicists and astronomers). As a result, the workshop grew to approximately 120 participants from several countries. The kind of atomic data that interests us are those parameters needed for analysis and modeling of spectra shortward of about about 100 A and relevant to ionic species of astronomical interest. The physical mechanisms of interest in the formation of spectra include photoionization. collisional ionization, recombination (radiative and dielectronic). collisional excitation (by electrons and protons). and radiative deexcitation. Unique to X-ray spectroscopy are the ionization and excitation processes from inner-closed shells. in addition to the challenges in interpret,ing the medium resolution (epsilon/delta epsilon is about 0.05 - 0.1) data obtained by current X-ray astronomy experiments. Line wavelengths are of interest too, particularly owing to the high resolution spectra from the new experiments. The workshop was divided into five major areas: Observational Spectroscopy, Theoretical Calculations of Atomic Data, Laboratory Measurements of Atomic Parameters. Spectra Modeling, and Atomic Databases. One comforting finding from the work shop is that the enthusiasm felt by X-ray astronomers about the new observational missions seems to be shared by theoretical and experimental physicists. Talks were presented about several exciting new projects and experimental and theoretical techniques devoted to X-ray spectroscopy. Simultaneously, several new tools for spectral analysis and modeling have recently been developed, together with improved atomic databases. These proceeding are expected to be of interests to producers and users of atomic data. Moreover. the contributions presented here have been written in a way that can be used by a general audience of scientists and graduate students in X-ray astronomy, modelling, and in computational and experimental atomic physics.

Kallman, Timothy↗

NASA NTRS2017

Diverse Long-Term Variability of Five Candidate High-Mass X-Ray Binaries from Swift Burst Alert Telescope Observations

We present an investigation of long-term modulation in the X-ray light curves of five little-studied candidate high-mass X-ray binaries using the Swift Burst Alert Telescope (SWIFT-BAT). IGR J14488-5942 and AX J1700.2-4220 show strong modulation at periods of 49.6 and 44 days, respectively, which are interpreted as orbital periods of Be star systems. For IGR J14488-5942, observations with the Swift X-ray Telescope show a hint of pulsations at 33.4 seconds. For AX J1700.2-4220, 54 second-pulsations were previously found with XMM-Newton. Swift J1816.7-1613 exhibits complicated behavior. The strongest peak in the power spectrum is at a period near 150 days, but this conflicts with a determination of a period of 118.5 days by La Parola et al. AX J1820.5-1434 has been proposed to exhibit modulation near 54 days, but the extended BAT observations suggest modulation at slightly longer than double this at approximately 111 days. There appears to be a long-term change in the shape of the modulation near 111 days, which may explain the apparent discrepancy. The X-ray pulsar XTE J1906+090,which was previously proposed to be a Be star system with an orbital period of approximately 30 days from pulse timing, shows peaks in the power spectrum at 81 and 173 days. The origins of these periods are unclear, although theymight be the orbital period and a superorbital period respectively. For all five sources, the long-term variability, together with the combination of orbital and proposed pulse periods, suggests that the sources contain Be starmass donors.

Corbet, Robin H. D.↗

NASA NTRS1995

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗

NASA NTRS2016

The NuSTAR View of the Non-Thermal Emission from PSR J0437-4715

We present a hard X-ray Nuclear Spectroscopic Telescope Array (NuSTAR) observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period approximately 5.76 ms are observed with a significance of 3.7sigma, at energies up to 20 keV above which the NuSTAR background dominates. We measure a photon index gamma = 1.50 +/- 0.25(90 per cent confidence) for the power-law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton spectrum of PSR J0437-4715, depending on the slope of the power-law component, and the amount of absorption of soft X-rays. The new constraint on the high-energy emission provided by NuSTAR removes ambiguities regarding the thermal components of the emission below 3 keV. We performed a simultaneous spectral analysis of the XMM-Newton and NuSTAR data to confirm that three thermal components and a power law are required to fit the 0.3-20 keV emission of PSR J0437-4715. Adding a ROSAT-PSPC spectrum further confirmed this result and allowed us to better constrain the temperatures of the three thermal components. A phase resolved analysis of the NuSTAR data revealed no significant change in the photon index of the high-energy emission. This NuSTAR observation provides further impetus for future observations with the NICER mission (Neutron Star Interior Composition Explorer) whose sensitivity will provide much stricter constraints on the equation of state of nuclear matter by combining model fits to the pulsars phase-folded light curve with the pulsars well-defined mass and distance from radio timing observations.

stars; neutro↗

DOE OSTI2021

Train Like a (Var)Pro: Efficient Training of Neural Networks with Variable Projection

Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Combustor exhaust-emissions and blowout-limits with diesel number 2 and jet A fuels utilizing air-atomizing and pressure atomizing nozzles

Experimental tests with diesel number 2 and Jet A fuels were conducted in a combustor segment to obtain comparative data on exhaust emissions and blowout limits. An air-atomizing nozzle was used to inject the fuels. Tests were also made with diesel number 2 fuel using a pressure-atomizing nozzle to determine the effectiveness of the air-atomizing nozzle in reducing exhaust emissions. Test conditions included fuel-air ratios of 0.008 to 0.018, inlet-air total pressures and temperatures of 41 to 203 newtons per square centimeter and 477 to 811 K, respectively, and a reference velocity of 21.3 meters per second. Smoke number and unburned hydrocarbons were twice as high with diesel number 2 as with Jet A fuel. This was attributed to diesel number 2 having a higher concentration of aromatics and lower volatility than Jet A fuel. Oxides of nitrogen, carbon monoxide, and blowout limits were approximately the same for the two fuels. The air-atomizing nozzle, as compared with the pressure-atomizing nozzle, reduced oxides-of-nitrogen by 20 percent, smoke number by 30 percent, carbon monoxide by 70 percent, and unburned hydrocarbons by 50 percent when used with diesel number 2 fuel.

Ingebo, R. D.↗

NASA NTRS1984

Comet rendezvous

A Mariner Mark II spacecraft rendezvous mission with comet Kopff has been recommended by NASA's Solar System Exploration Committee, and is scheduled for a Shuttle launch in 1990. The spacecraft, which is scheduled to encounter the comet in 1994, will conduct a series of experiments including the study of the cometary nucleus, coma, and tail during perihelion passage; the determination of the nucleus and coma chemical/isotropic composition; the description of the nucleus's size, mass, rotation period, and pole occultation; distinctions of gas and dust hydrodynamics; and distinctions of solar wind interactions with the coma. The three-axis-stabilized spacecraft will be modular in design, and will incorporate a 490-Newton Insat engine, Viking Orbiter propulsion tanks, an Integrated Platform Pointing and Attitude Control Subsystem, and several scientific instruments. The spacecraft's Radio Frequency Subsystem will only use X-band telemetry with a downlink frequency of 8415 MHz, and an uplink frequency of 7161 MHz. The power sources in the present design include one radioisotope thermoelectric generator, three 3 Ah batteries, and a solar panel of approximately 7 sq m.

Draper, R. F.↗

DOE OSTI2021

Magnetars and axion-like particles: probes with the hard X-ray spectrum

Quiescent hard X-ray and soft gamma-ray emission from neutron stars constitute a promising frontier to explore axion-like-particles (ALPs). ALP production in the core peaks at energies of a few keV to a few hundreds of keV; subsequently, the ALPs escape and convert to photons in the magnetosphere. The emissivity goes as ~T 6 while the conversion probability is enhanced for large magnetic fields, making magnetars, with their high core temperatures and strong magnetic fields, ideal targets for probing ALPs. We compute the energy spectrum of photons resulting from conversion of ALPs in the magnetosphere and then compare it against hard X-ray data from NuSTAR, INTEGRAL, and XMM-Newton, for a set of eight magnetars for which such data exists. Upper limits are placed on the product of the ALP-nucleon and ALP-photon couplings. For the production in the core, we perform a calculation of the ALP emissivity in degenerate nuclear matter modeled by a relativistic mean field theory. The reduction of the emissivity due to improvements to the one-pion exchange approximation is incorporated, as is the suppression of the emissivity due to proton superfluidity in the neutron star core. A range of core temperatures is considered, corresponding to different models of the steady heat transfer from the core to the stellar surface. For the subsequent conversion, we solve the coupled differential equations mixing ALPs and photons in the magnetosphere. The conversion occurs due to a competition between the dipolar magnetic field and the photon refractive index induced by the external magnetic field. Semi-analytic expressions are provided alongside the full numerical results. Furthermore, we also present an analysis of the uncertainty on the axion limits we derive due to the uncertainties in the magnetar masses, nuclear matter equation of state, and the proton superfluid critical temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗