Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient methods”

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
DOE OSTI2015

SU-C-BRA-02: Gradient Based Method of Target Delineation On PET/MR Image of Head and Neck Cancer Patients

Purpose: Validate the consistency of a gradient-based segmentation tool to facilitate accurate delineation of PET/CT-based GTVs in head and neck cancers by comparing against hybrid PET/MR-derived GTV contours. Materials and Methods: A total of 18 head and neck target volumes (10 primary and 8 nodal) were retrospectively contoured using a gradient-based segmentation tool by two observers. Each observer independently contoured each target five times. Inter-observer variability was evaluated via absolute percent differences. Intra-observer variability was examined by percentage uncertainty. All target volumes were also contoured using the SUV percent threshold method. The thresholds were explored case by case so its derived volume matched with the gradient-based volume. Dice similarity coefficients (DSC) were calculated to determine overlap of PET/CT GTVs and PET/MR GTVs. Results: The Levene’s test showed there was no statistically significant difference of the variances between the observer’s gradient-derived contours. However, the absolute difference between the observer’s volumes was 10.83%, with a range from 0.39% up to 42.89%. PET-avid regions with qualitatively non-uniform shapes and intensity levels had a higher absolute percent difference near 25%, while regions with uniform shapes and intensity levels had an absolute percent difference of 2% between observers. The average percentage uncertainty between observers was 4.83% and 7%. As the volume of the gradient-derived contours increased, the SUV threshold percent needed to match the volume decreased. Dice coefficients showed good agreement of the PET/CT and PET/MR GTVs with an average DSC value across all volumes at 0.69. Conclusion: Gradient-based segmentation of PET volume showed good consistency in general but can vary considerably for non-uniform target shapes and intensity levels. PET/CT-derived GTV contours stemming from the gradient-based tool show good agreement with the anatomically and metabolically more accurate PET/MR-derived GTV contours, but tumor delineation accuracy can be further improved with the use PET/MR.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

A posteriori superlinear convergence bounds for block conjugate gradient

In this paper, we extend to the block case the a posteriori bound showing superlinear convergence of the conjugate gradient method developed by van der Vorst and Vuik in [J. Comput. Applied Math., 48 (1993), pp. 327–341]. That is, we obtain similar bounds but now for the block conjugate gradient method. We also present a series of computational experiments, illustrating the validity of the bound developed here as well as the bound by Simoncini and Szyld from [SIAM Review, 47 (2005), pp. 247–272] using angles between subspaces. Using these bounds, we make some observations on the onset of superlinearity and how this onset depends on the eigenvalue distribution and the block size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Force field preconditioned ab initio structure relaxation method

We present a general method to accelerate ab initio atomic structural relaxations. In this method, the conjugate gradient ab initio relaxation is preconditioned by a force field relaxation. This force field is constructed on-the-flight with the information of the ab initio forces and the current atomic configuration. At each ab initio relaxation step, the so-constructed force field is used to prerelax the system, with the relaxation direction as the preconditioned direction for the ab initio conjugated gradient method. The force field model and its parameters are rather general making this method applicable for a wide range of systems commonly used in material simulations. We report more than 80 different systems have been tested representing different cases in material simulations. Across the board, we find accelerations mostly over a factor of 2 and for some large systems with a factor of 6 ~ 9 . There is no case where the relaxation becomes worse. The code and lib for this method are provided, which can be used as a plugin in a standard ab initio atomic relaxation code.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1973

Optimal control of a variable spin speed CMG system for space vehicles

Many future NASA programs require very high accurate pointing stability. These pointing requirements are well beyond anything attempted to date. This paper suggests a control system which has the capability of meeting these requirements. An optimal control law for the suggested system is specified. However, since no direct method of solution is known for this complicated system, a computation technique using successive approximations is used to develop the required solution. The method of calculus of variations is applied for estimating the changes of index of performance as well as those constraints of inequality of state variables and terminal conditions. Thus, an algorithm is obtained by the steepest descent method and/or conjugate gradient method. Numerical examples are given to show the optimal controls.

Liu, T. C.↗
NASA NTRS2004

Putting Priors in Mixture Density Mercer Kernels

This paper presents a new methodology for automatic knowledge driven data mining based on the theory of Mercer Kernels, which are highly nonlinear symmetric positive definite mappings from the original image space to a very high, possibly infinite dimensional feature space. We describe a new method called Mixture Density Mercer Kernels to learn kernel function directly from data, rather than using predefined kernels. These data adaptive kernels can en- code prior knowledge in the kernel using a Bayesian formulation, thus allowing for physical information to be encoded in the model. We compare the results with existing algorithms on data from the Sloan Digital Sky Survey (SDSS). The code for these experiments has been generated with the AUTOBAYES tool, which automatically generates efficient and documented C/C++ code from abstract statistical model specifications. The core of the system is a schema library which contains template for learning and knowledge discovery algorithms like different versions of EM, or numeric optimization methods like conjugate gradient methods. The template instantiation is supported by symbolic- algebraic computations, which allows AUTOBAYES to find closed-form solutions and, where possible, to integrate them into the code. The results show that the Mixture Density Mercer-Kernel described here outperforms tree-based classification in distinguishing high-redshift galaxies from low- redshift galaxies by approximately 16% on test data, bagged trees by approximately 7%, and bagged trees built on a much larger sample of data by approximately 2%.

Srivastava, Ashok N.↗
NASA NTRS2005

An Ensemble Approach to Building Mercer Kernels with Prior Information

This paper presents a new methodology for automatic knowledge driven data mining based on the theory of Mercer Kernels, which are highly nonlinear symmetric positive definite mappings from the original image space to a very high, possibly dimensional feature space. we describe a new method called Mixture Density Mercer Kernels to learn kernel function directly from data, rather than using pre-defined kernels. These data adaptive kernels can encode prior knowledge in the kernel using a Bayesian formulation, thus allowing for physical information to be encoded in the model. Specifically, we demonstrate the use of the algorithm in situations with extremely small samples of data. We compare the results with existing algorithms on data from the Sloan Digital Sky Survey (SDSS) and demonstrate the method's superior performance against standard methods. The code for these experiments has been generated with the AUTOBAYES tool, which automatically generates efficient and documented C/C++ code from abstract statistical model specifications. The core of the system is a schema library which contains templates for learning and knowledge discovery algorithms like different versions of EM, or numeric optimization methods like conjugate gradient methods. The template instantiation is supported by symbolic-algebraic computations, which allows AUTOBAYES to find closed-form solutions and, where possible, to integrate them into the code.

Srivastava, Ashok N.↗
DOE OSTI2013

Approximate error conjugation gradient minimization methods

In one embodiment, a method includes selecting a subset of rays from a set of all rays to use in an error calculation for a constrained conjugate gradient minimization problem, calculating an approximate error using the subset of rays, and calculating a minimum in a conjugate gradient direction based on the approximate error. In another embodiment, a system includes a processor for executing logic, logic for selecting a subset of rays from a set of all rays to use in an error calculation for a constrained conjugate gradient minimization problem, logic for calculating an approximate error using the subset of rays, and logic for calculating a minimum in a conjugate gradient direction based on the approximate error. In other embodiments, computer program products, methods, and systems are described capable of using approximate error in constrained conjugate gradient minimization problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

A Novel Method Of Gradient Forming and Fluid Manipulation in Reduced Gravity Environments

The use of magnetic fields to control the motion and position of non-conducting liquids has received growing interest in recent times. The possibility of using the forces exerted by a nonuniform magnetic field on a ferrofluid to not only achieve fluid manipulation but also to actively control fluid motion makes it an attractive candidate for applications such as heat transfer in space systems. Terrestrial heat transfer equipment often relies on the normal gravitational force to hold liquid in a desired position or to provide a buoyant force to enhance the heat transfer rate. The residual gravitational force present in a space environment may no longer serve these useful functions and other forces, such as surface tension, can play a significant role in determining heat transfer rates. Although typically overwhelmed by gravitational forces in terrestrial applications, the body force induced in a ferrofluid by a nonuniform magnetic field can help to achieve these objectives in a microgravity environment. This paper will address the fluid manipulation aspect and will comprise of results from model fluid experiments and numerical modeling of the problem. Results from a novel method of forming concentration gradients that are applicable to low gravity applications will be presented. The ground based experiments are specifically tailored to demonstrate the magnetic manipulation capability of a ferrofluid and show that gravitational effects can be countered in carefully designed systems. The development of governing equations for the system will be presented along with a sampling of numerical results.

Ramachandran N.↗
DOE OSTI2019

Distributed Energy Resources Coordination over Time-varying Directed Communication Networks

In this paper, we consider an optimal coordination problem for distributed energy resources (DERs) including distributed generators and energy storage devices. We first propose an algorithm based on the push-sum and gradient method to optimally coordinate storage devices and distributed generators in a distributed manner. In the proposed algorithm, each DER only maintains a set of variables and updates them through information exchange with a few neighboring DERs over a time-varying directed communication network. We show that the proposed distributed algorithm with appropriately chosen step-sizes solve the optimal DER coordination problem if the time-varying directed communication network is uniformly jointly strongly connected, which is a mild condition on the connectivity of communication topologies. The proposed distributed algorithm is illustrated and validated by numerical simulations.

Yang, Tao↗
DOE OSTI2021

A Scalable Gradient Free Method for Bayesian Experimental Design with Implicit Models

Bayesian experimental design (BED) is to answer the question that how to choose designs that maximize the information gathering. For implicit models, where the likelihood is intractable but sampling is possible, conventional BED methods have difficulties in efficiently estimating the posterior distribution and maximizing the mutual information (MI) between data and parameters. Recent work proposed the use of gradient ascent to maximize a lower bound on MI to deal with these issues. However, the approach requires a sampling path to compute the pathwise gradient of the MI lower bound with respect to the design variables, and such a pathwise gradient is usually inaccessible for implicit models. In this paper, we propose a novel approach that leverages recent advances in stochastic approximate gradient ascent incorporated with a smoothed variational MI estimator for efficient and robust BED. Without the necessity of pathwise gradients, our approach allows the design process to be achieved through a unified procedure with an approximate gradient for implicit models. Several experiments show that our approach outperforms baseline methods, and significantly improves the scalability of BED in high-dimensional problems.

Zhang, Jiaxin↗
NASA NTRS1979

Induced electric currents in the Alaska oil pipeline measured by gradient, fluxgate, and SQUID magnetometers

The field gradient method for observing the electric currents in the Alaska pipeline provided consistent values for both the fluxgate and SQUID method of observation. These currents were linearly related to the regularly measured electric and magnetic field changes. Determinations of pipeline current were consistent with values obtained by a direct connection, current shunt technique at a pipeline site about 9.6 km away. The gradient method has the distinct advantage of portability and buried- pipe capability. Field gradients due to the pipe magnetization, geological features, or ionospheric source currents do not seem to contribute a measurable error to such pipe current determination. The SQUID gradiometer is inherently sensitive enough to detect very small currents in a linear conductor at 10 meters, or conversely, to detect small currents of one amphere or more at relatively great distances. It is fairly straightforward to achieve imbalance less than one part in ten thousand, and with extreme care, one part in one million or better.

Campbell, W. H.↗
DOE OSTI2018

A new accelerated algorithm for ill-conditioned ridge regression problems

In this work, we study the primal and dual formulations of the regularized least squares problem, in special norm L{sub 2}, named ridge regression. We state the Gradient method for the corresponding primal-dual problem through a fixed point approach. For this formulation, we present an original convergence analysis involving the spectral radius of a suitable iteration matrix. The main contribution of this work is to accelerate the studied method by means of a memory-like strategy. Some connections with the spectral properties of the associated iteration matrices were established to prove the convergence of the accelerated method. Preliminary experiments showed the good performance of the proposed algorithm when applied to solve ill-conditioned problems, providing better results than the conjugate gradient method for this class of problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Nucleon-pair approximation for nuclei from spherical to deformed regions

Here we model low-lying states of atomic nuclei in the nucleon-pair approximation of the shell model, using three approaches to select collective nucleon pairs: the generalized seniority scheme, the conjugate gradient method, and the Hartree-Fock approach. We find the collective pairs obtained from the generalized seniority scheme provide a good description for nearly spherical nuclei, and those from the conjugate gradient method or the Hartree-Fock approach work well for transitional and deformed nuclei. Our NPA calculations using collective pairs with angular momenta 0, 2, and 4 (denoted by S D G pairs) reproduce the nuclear shape evolution in the N = 26 isotones, Ca 46 , Ti 48 , Cr 50 , and Fe 52 , and yield good agreement with full configuration-interaction calculations of low-lying states in medium-heavy transitional and deformed nuclei: Ti 44 – 48 , Cr 48 , Cr 50 , Fe 52 , Zn 60 – 64 , Ge 64 , 66 , Mo 84 , and Xe 108 – 112 . Finally, using the S D G I -pair approximation we describe low-lying states of Ba 112 , 114 , cases difficult to reach by conventional configuration-interaction methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

Comparison of solution approaches for minimum-fuel, low-thrust, power-limited orbital transfers

An initial assessment of the feasibility of a function space gradient method for computing solutions to minimum-fuel power-limited transfers encompassing a wide range of thrust to weight ratios is conducted. Three transfers between coplanar ellipses are used as test cases. The gradient method performs best at the high end of the thrust to weight ratio range. At the lower end, there is reduced sensitivity of the fuel consumption to the control profiles. The minimum fuel consumption and the trajectory are computed quite accurately but the control profiles are in error. An approximate analytical solution, obtained by Edelbaum using the method of averaging, is discussed.

Mease, Kenneth D.↗
DOE OSTI2013

Porous media properties of reticulated shape memory polymer foams and mock embolic coils for aneurysm treatment

Background: Shape memory polymer (SMP) foams are being investigated as an alternative aneurysm treatment method to embolic coils. The goal of both techniques is the reduction of blood flow into the aneurysm and the subsequent formation of a stable thrombus, which prevents future aneurysm rupture. The purpose of this study is to experimentally determine the parameters, permeability and form factor, which are related to the flow resistance imposed by both media when subjected to a pressure gradient. Methods: The porous media properties—permeability and form factor—of SMP foams and mock embolic coils (MECs) were measured with a pressure gradient method by means of an in vitro closed flow loop. We implemented the Forchheimer-Hazen-Dupuit-Darcy equation to calculate these properties. Mechanically-reticulated SMP foams were fabricated with average cell sizes of 0.7E-3 and 1.1E-3 m, while the MECs were arranged with volumetric packing densities of 11-28%. Results: The permeability of the SMP foams was an order of magnitude lower than that of the MECs. The form factor differed by up to two orders of magnitude and was higher for the SMP foams in all cases. The maximum flow rate of all samples tested was within the inertial laminar flow regime, with Reynolds numbers ranging between 1 and 35. Conclusions: The SMP foams impose a greater resistance to fluid flow compared to MECs, which is a result of increased viscous and inertial losses. These results suggest that aneurysms treated with SMP foam will have flow conditions more favorable for blood stasis than those treated with embolic coils having packing densities ≤ 28%.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A gradient based method for modeling baryons and matter in halos of fast simulations

Fast N-body PM simulations with a small number of time steps such as FastPM or COLA have been remarkably successful in modeling the galaxy statistics, but their lack of small scale force resolution and long time steps cannot give accurate halo matter profiles or matter power spectrum. High resolution N-body simulations can improve on this, but lack baryonic effects, which can only be properly included in hydro simulations. Here we present a scheme to calibrate the fast simulations to mimic the precision of the hydrodynamic simulations or high resolution N-body simulations. The scheme is based on a gradient descent of either effective gravitational potential, which mimics the short range force, or of effective enthalpy, which mimics gas hydrodynamics and feedback. The scheme is fast and differentiable, and can be incorporated as a post-processing step into any simulation. It gives very good results for the matter power spectrum for several of the baryonic feedback and dark matter simulations, and also gives improved dark matter halo profiles. The scheme is even able to find the large subhalos, and increase the correlation coefficient between the fast simulations and the high resolution N-body or hydro simulations. It can also be used to add baryonic effects to the high resolution N-body simulations. While the method has free parameters that can be calibrated on various simulations, they can also be viewed as astrophysical nuisance parameters describing baryonic effects that can be marginalized over during the data analysis. In this view these parameters can be viewed as an efficient parametrization of baryonic effects.

79 ASTRONOMY AND ASTROPHYSICS↗