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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Experiences on p-Version Time-Discontinuous Galerkin's Method for Nonlinear Heat Transfer Analysis and Sensitivity Analysis

The focus of this research is on the development of analysis and sensitivity analysis equations for nonlinear, transient heat transfer problems modeled by p-version, time discontinuous finite element approximation. The resulting matrix equation of the state equation is simply in the form ofA(x)x = c, representing a single step, time marching scheme. The Newton-Raphson's method is used to solve the nonlinear equation. Examples are first provided to demonstrate the accuracy characteristics of the resultant finite element approximation. A direct differentiation approach is then used to compute the thermal sensitivities of a nonlinear heat transfer problem. The report shows that only minimal coding effort is required to enhance the analysis code with the sensitivity analysis capability.

Hou, Gene↗
NASA NTRS1988

Implicit methods in CFD

A class of implicit approximate factorization schemes is examined for stability and convergence characteristics. These schemes include Newton's method, factorization, and flux-vector splitting. Examples are used to show that all practical methods suffer from some limited stability or asymptotic convergence restriction. A three-dimensional factored scheme is shown which suffers from unconditional instability which can only be ameliorated by added artificial dissipation. An F3D + or - flux split scheme is described which avoids unconditional instability, but in the end has similar convergence characteristics.

Pulliam, Thomas H.↗
NASA NTRS1991

Numerical solutions of 2-D multi-stage rotor/stator unsteady flow interactions

The Rai method of single-stage rotor/stator flow interaction is extended to handle multistage configurations. In this study, a two-dimensional Navier-Stokes multi-zone approach was used to investigate unsteady flow interactions within two multistage axial turbines. The governing equations are solved by an iterative, factored, implicit finite-difference, upwind algorithm. Numerical accuracy is checked by investigating the effect of time step size, the effect of subiteration in the Newton-Raphson technique, and the effect of full viscous versus thin-layer approximation. Computer results compared well with experimental data. Unsteady flow interactions, wake cutting, and the associated evolution of vortical entities are discussed.

Yang, R.-J.↗
NASA NTRS2016

NASA's Space Launch System Marks Critical Design Review

With completion of its Critical Design Review (CDR) in 2015, NASA is deep into the manufacturing and testing phases of its new Space Launch System (SLS) for beyond-Earth exploration. This CDR was the first in almost 40 years for a NASA human launch vehicle and marked another successful milestone on the road to the launch of a new era of deep space exploration. The review marked the 90-percent design-complete, a final look at the design and development plan of the integrated vehicle before full-scale fabrications begins and the prelude to the next milestone, design certification. Specifically, the review looked at the first of three increasingly capable configurations planned for SLS. This "Block I" design will stand 98.2 meters (m) (322 feet) tall and provide 39.1 million Newtons (8.8 million pounds) of thrust at liftoff to lift a payload of approximately 70 metric tons (154,000 pounds). This payload is more than double that of the retired space shuttle program or other current launch vehicles. It dramatically increases the mass and volume of human and robotic exploration. Additionally, it will decrease overall mission risk, increase safety, and simplify ground and mission operations - all significant considerations for crewed missions and unique, high-value national payloads. The Block 1 SLS will launch NASA's Orion Multi-Purpose Crew Vehicle (MPCV) on an uncrewed flight beyond the moon and back and the first crewed flight around the moon. The current design has a direct evolutionary path to a vehicle with a 130t lift capability that offers even more flexibility to reduce planetary trip times, simplify payload design cycles, and provide new capabilities such as planetary sample returns. Every major element of SLS has hardware in production or testing, including flight hardware for the Exploration 1 (EM-1) test flight. In fact, the SLS MPCV-to-Stage-Adapter (MSA) flew successfully on the Exploration Flight Test (EFT) 1 launch of a Delta IV and Orion spacecraft in December 2014. This paper will discuss these and other technical and programmatic successes and challenges over the past year and provide a preview of work ahead before the first flight of this new capability.

Singer, Chris↗
NASA NTRS1992

Embedding impedance approximations in the analysis of SIS mixers

Future millimeter-wave radio astronomy instruments will use arrays of many SIS receivers, either as focal plane arrays on individual radio telescopes, or as individual receivers on the many antennas of radio interferometers. Such applications will require broadband integrated mixers without mechanical tuners. To produce such mixers, it will be necessary to improve present mixer design techniques, most of which use the three-frequency approximation to Tucker's quantum mixer theory. This paper examines the adequacy of three approximations to Tucker's theory: (1) the usual three-frequency approximation which assumes a sinusoidal LO voltage at the junction, and a short-circuit at all frequencies above the upper sideband; (2) a five-frequency approximation which allows two LO voltage harmonics and five small-signal sidebands; and (3) a quasi five-frequency approximation in which five small-signal sidebands are allowed, but the LO voltage is assumed sinusoidal. These are compared with a full harmonic-Newton solution of Tucker's equations, including eight LO harmonics and their corresponding sidebands, for realistic SIS mixer circuits. It is shown that the accuracy of the three approximations depends strongly on the value of omega R(sub N)C for the SIS junctions used. For large omega R(sub N)C, all three approximations approach the eight-harmonic solution. For omega R(sub N)C values in the range 0.5 to 10, the range of most practical interest, the quasi five-frequency approximation is a considerable improvement over the three-frequency approximation, and should be suitable for much design work. For the realistic SIS mixers considered here, the five-frequency approximation gives results very close to those of the eight-harmonic solution. Use of these approximations, where appropriate, considerably reduces the computational effort needed to analyze an SIS mixer, and allows the design and optimization of mixers using a personal computer.

Kerr, A. R.↗
NASA NTRS1984

Maximum likelihood algorithm using an efficient scheme for computing sensitivities and parameter confidence intervals

Improved techniques for estimating airplane stability and control derivatives and their standard errors are presented. A maximum likelihood estimation algorithm is developed which relies on an optimization scheme referred to as a modified Newton-Raphson scheme with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort compared to integrating the analytically-determined sensitivity equations or using a finite difference scheme. An aircraft estimation problem is solved using real flight data to compare MNRES with the commonly used modified Newton-Raphson technique; MNRES is found to be faster and more generally applicable. Parameter standard errors are determined using a random search technique. The confidence intervals obtained are compared with Cramer-Rao lower bounds at the same confidence level. It is observed that the nonlinearity of the cost function is an important factor in the relationship between Cramer-Rao bounds and the error bounds determined by the search technique.

Murphy, P. C.↗
NASA NTRS1998

Congruence Approximations for Entrophy Endowed Hyperbolic Systems

Building upon the standard symmetrization theory for hyperbolic systems of conservation laws, congruence properties of the symmetrized system are explored. These congruence properties suggest variants of several stabilized numerical discretization procedures for hyperbolic equations (upwind finite-volume, Galerkin least-squares, discontinuous Galerkin) that benefit computationally from congruence approximation. Specifically, it becomes straightforward to construct the spatial discretization and Jacobian linearization for these schemes (given a small amount of derivative information) for possible use in Newton's method, discrete optimization, homotopy algorithms, etc. Some examples will be given for the compressible Euler equations and the nonrelativistic MHD equations using linear and quadratic spatial approximation.

Barth, Timothy J.↗
NASA NTRS1993

An efficient nonlinear relaxation technique for the three-dimensional, Reynolds-averaged Navier-Stokes equations

An efficient implicit method for the computation of steady, three-dimensional, compressible Navier-Stokes flowfields is presented. A nonlinear iteration strategy based on planar Gauss-Seidel sweeps is used to drive the solution toward a steady state, with approximate factorization errors within a crossflow plane reduced by the application of a quasi-Newton technique. A hybrid discretization approach is employed, with flux-vector splitting utilized in the streamwise direction and central differences with artificial dissipation used for the transverse fluxes. Convergence histories and comparisons with experimental data are presented for several 3-D shock-boundary layer interactions. Both laminar and turbulent cases are considered, with turbulent closure provided by a modification of the Baldwin-Barth one-equation model. For the problems considered (175,000-325,000 mesh points), the algorithm provides steady-state convergence in 900-2000 CPU seconds on a single processor of a Cray Y-MP.

Edwards, Jack R.↗
NASA NTRS2016

NuSTAR and XMM-Newton Observations of 1E1743.1-2843: Indications of a Neutron Star LMXB Nature of the Compact Object

We report on the results of NuSTAR and XMM-Newton observations of the persistent X-ray source 1E1743.1-2843, located in the Galactic Center region. The source was observed between 2012 September and October by NuSTAR and XMM-Newton, providing almost simultaneous observations in the hard and soft X-ray bands. The high X-ray luminosity points to the presence of an accreting compact object. We analyze the possibilities of this accreting compact object being either a neutron star (NS) or a black hole, and conclude that the joint XMM-Newton and NuSTAR spectrum from 0.3 to 40 keV fits a blackbody spectrum with kT approximately 1.8 keV emitted from a hot spot or an equatorial strip on an NS surface. This spectrum is thermally Comptonized by electrons with kTe approximately 4.6 keV. Accepting this NS hypothesis, we probe the low-mass X-ray binary (LMXB) or high-mass X-ray binary (HMXB) nature of the source. While the lack of Type-I bursts can be explained in the LMXB scenario, the absence of pulsations in the 2 MHz-49 Hz frequency range, the lack of eclipses and of an IR companion, and the lack of a Kaline from neutral or moderately ionized iron strongly disfavor interpreting this source as a HMXB. We therefore conclude that 1E1743.1-2843 is most likely an NS-LMXB located beyond the Galactic Center. There is weak statistical evidence for a soft X-ray excess which may indicate thermal emission from an accretion disk. However, the disk normalization remains unconstrained due to the high hydrogen column density (N(sub H) approximately 1.6 x 10(exp 23) cm(exp -2)).

accretion↗
NASA NTRSDate not supplied

Implicit Preconditioning for Explicit Multigrid Solvers on Cut-Cell Cartesian Meshes

This work assesses the effectiveness of linearized implicit Euler preconditioning for multigrid solvers using an unpreconditioned, Jacobian-free Newton Krylov method to converge the linear system of equations. Multigrid convergence rates improve to approximately 0.75 across the cases tested including a Mach 2 supersonic wedge, transonic NACA 0012 airfoil, and ONERA M6 wing. While larger Krylov subspaces increase the convergence rate, they also increase the computational cost, such that 4-8 Krylov vectors often offers the fastest turnaround. Further reductions in computational cost are achieved with a sequential hybrid preconditioner that begins with the explicit multigrid solver before transitioning to the preconditioned algorithm later on. In addition, a novel implementation of dual time stepping is extended to include both common BDF methods as well as high-order implicit Runge-Kutta schemes. This particular formulation, which uses A −1 preconditioning, is amenable to matrix-free solvers, and the L-stable methods are especially suited for meshes with arbitrarily small cut-cells. Asymptotic order of convergence is demonstrated for BDF1, BDF2, SDIRK2, and 3rd-order Radau IIA time integration with unsteady 2D vortex simulations.

ARMD↗
NASA NTRS1982

Comparison of starting values for iterative solutions to a universal Kepler's equation

General starting values for the iterative numerical solution of a universal Kepler's equation for position in a conic orbit at a specified time are investigated. Three starting values based on recent refinements of previously obtained bounds on the solution are derived and tested numerically. Of these, a simple starting value based on a cubic approximation to Kepler's equation provides the most rapid convergence using both first and second order Newton algorithms. The performance of the starting values are compared with similar studies which used the restricted case of elliptical orbits with the initial epoch at periapse.

Bergam, M. J.↗
NASA NTRS2004

Mass and Dynamical Structures of the Lensing Clusters CL0024+17 and CL2244+02

We present a detailed gravitational mass measurement based on the XMM-Newton imaging spectroscopy analysis of the lensing cluster of galaxies CL0024+17 at $z = 0.395$. The emission appears approximately symmetric. However, on the scale of $r\sim3.3'$, some indication of elongation is visible in the northwest-southeast direction from the hardness ratio map. Within $3'$, we measure a global gas temperature of $3.52\pm0.17$ keV, metallicity of $0.22\pm0.07$, and a bolometric luminosity of $2.9\pm0. l\times10(exp 44)$ erg/s. We derive a temperature distribution with an isothermal temperature of 3.9 keV up to a radius of $1.5'$ and a strong temperature gradient in the outskirts ($1.3' less than r less than 3.3'$). Under the assumption of hydrostatic equilibrium, we measure the gravitational mass and gas mass fraction to be $M-{200} = 2.0\pm0.3\times 10(exp 14)$ solar masses and $f-{gas} = 0.20\pm0.03$ at $r-{200} = 1.05$ Mpc (all for a Hubble constant of 70 km/sec/Mpc) using the observed gas temperature profile. The complex core structure is the key to explaining the discrepancy between the gravitational mass determined from the XMM-Newton observations and HST optical lensing measurements.

Forman, William↗
DOE OSTI2020

BISON Robustness and Performance Improvements

BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2002

Higher Order Time Integration Schemes for the Unsteady Navier-Stokes Equations on Unstructured Meshes

The efficiency gains obtained using higher-order implicit Runge-Kutta schemes as compared with the second-order accurate backward difference schemes for the unsteady Navier-Stokes equations are investigated. Three different algorithms for solving the nonlinear system of equations arising at each timestep are presented. The first algorithm (NMG) is a pseudo-time-stepping scheme which employs a non-linear full approximation storage (FAS) agglomeration multigrid method to accelerate convergence. The other two algorithms are based on Inexact Newton's methods. The linear system arising at each Newton step is solved using iterative/Krylov techniques and left preconditioning is used to accelerate convergence of the linear solvers. One of the methods (LMG) uses Richardson's iterative scheme for solving the linear system at each Newton step while the other (PGMRES) uses the Generalized Minimal Residual method. Results demonstrating the relative superiority of these Newton's methods based schemes are presented. Efficiency gains as high as 10 are obtained by combining the higher-order time integration schemes with the more efficient nonlinear solvers.

Jothiprasad, Giridhar↗
NASA NTRS1957

Bodies of Revolution Having Minimum Drag at High Supersonic Airspeeds

Approximate shapes of nonlifting bodies having minimum pressure foredrag at high supersonic airspeeds are calculated. With the aid of Newton's law of resistance, the investigation is carried out for various combinations of the conditions of given body length, base diameter, surface area, and volume. In general, it is found that when body length is fixed, the body has a blunt nose; whereas, when the length is not fixed, the body has a sharp nose. The additional effect of curvature of the flow over the surface is investigated to determine its influence on the shapes for minimum drag. The effect is to increase the bluntness of the shapes in the region of the nose and the curvature in the region downstream of the nose. These shape modifications have, according to calculation, only a slight tendency to reduce drag. Several bodies of revolution of fineness ratios 3 and 5, including the calculated shapes of minimum drag for given length and base diameter and for given base diameter and surface area, were tested at Mach numbers from 2.73 to 6.28. A comparison of theoretical and experimental foredrag coefficients indicates that the calculated minimum-drag bodies are reasonable approximations to the correct shape.

Eggers, A J , Jr↗
NASA NTRS2014

The Soft-X-Ray Emission of Ark 120. XMM-Newton, NuSTAR, and the Importance of Taking the Broad View

We present simultaneous XMM-Newton and NuSTAR observations of the 'bare' Seyfert 1 galaxy, Ark 120, a system in which ionized absorption is absent. The NuSTAR hard-X-ray spectral coverage allows us to constrain different models for the excess soft-X-ray emission. Among phenomenological models, a cutoff power law best explains the soft-X-ray emission. This model likely corresponds to Comptonization of the accretion disc seed UV photons by a population of warm electrons: using Comptonization models, a temperature of approximately 0.3 kiloelectronvolts and an optical depth of approximately 13 are found. If the UV-to-X-ray OPTXAGNF model is applied, the UV fluxes from the XMM-Newton Optical Monitor suggest an intermediate black hole spin. Contrary to several other sources observed by NuSTAR, no high-energy cutoff is detected with a lower limit of 190 kiloelectronvolts.

accretion↗
NASA NTRS1988

Advances in dual algorithms and convex approximation methods

A new algorithm for solving the duals of separable convex optimization problems is presented. The algorithm is based on an active set strategy in conjunction with a variable metric method. This first order algorithm is more reliable than Newton's method used in DUAL-2 because it does not break down when the Hessian matrix becomes singular or nearly singular. A perturbation technique is introduced in order to remove the nondifferentiability of the dual function which arises when linear constraints are present in the approximate problem.

Smaoui, H.↗
DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗