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 937 records · Page 52

First principles study of the structural, elastic, electronic and phonon properties of CdX{sub 2}O{sub 4} (X=Al, Ga, In) spinel-type oxides

We have performed ab-initio calculations of the structural, electronic, elastic and dynamical properties for the spinel compounds CdX{sub 2}O{sub 4} (X=Al, Ga, In) using the plane wave pseudo-potential method within the generalized gradient approximation (GGA). The calculated lattice parameters, elastic constants for these compounds are in good agreement with the previous calculated values. The computed direct band gaps of CdAl{sub 2}O{sub 4}, CdGa{sub 2}O{sub 4} and CdIn{sub 2}O{sub 4} are 2.90 eV, 1.92 eV and 1.16 eV, respectively. The lattice vibrations were calculated by direct method. The calculated phonon dispersion curves show that all compounds are dynamically stable in the spinel structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Volatility Studies of Lanthanide Hexafluoroacetylacetone Chelates

Nuclear forensics has become very important to the field of nuclear security, where there is considerable interest in developing rapid separations of radionuclides and isotopes. Although prevention and detection form the foundation of nuclear security science, in the unfortunate event of detonation of a nuclear device, the fallout debris can be analyzed for signature isotopes that may provide evidence on the origin and type of nuclear materials. This is of considerable political importance to many governing bodies and this information is desired in a timely manner. One method of analysis is the rapid gas phase separation of rare earth metals from other components of the debris via isothermal thermochromatography. Lanthanides cover a portion of the fission product curve with varying yields dependent on a number of variables (e.g. neutron energy, fissile parent, etc), motivating the chromatographic separation of these elements. In previous work, 1,1,1,5,5,5 - hexafluoroacetylacetone (hfac) has been shown to produce volatile complexes with lanthanides. Thermochromatography is a gas phase separation method that utilizes temperature gradients on a column to vary adsorption location or constant column temperatures to vary retention time. The adsorption time of a compound on a column, i.e. retention time, is dependent on the thermodynamics, most notably the adsorption enthalpy. Studies of heavy metals have generated empirical correlations between adsorption enthalpy on a quartz column and sublimation enthalpy of the solid material. Thermogravimetric analysis (TGA) is useful in the characterization of the thermodynamics of materials, especially the sublimation enthalpy. Many methods analyze weight loss curves as a function of a constant heating rate. Another method introduced by Ashcroft analyzes weight loss curves over isothermal heating regions. Work presented here focuses on the determination of the sublimation enthalpies of thirteen lanthanide hfac complexes via the latter method. Additionally, there is investigation into the degradation of the compound via constant heating TGA up to higher temperatures. This provides information on the onset and end-set temperatures of degradation stages, as well as the peak temperature at which total sublimation occurs. (authors)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of the displacement field along a surface crack in a flat plate using optical techniques

Structural components may contain cracks, and for these components a crack tip represents a highly singular stress. The evaluation of the associated strain gradient is difficult to achieve with experimental discrete methods. An efficient alternative are the optical techniques which are non-contact, and so delivering continuous information of the displacement fields and associated derivatives for strain evaluation. This paper describes some experimental methods to fully characterize the displacement field near a crack tip existing in flat plates. Three optical field techniques based on image analysis were used in the present work; respectively, Electronic Speckle Pattern Interferometry, Moiré Interferometry, and Digital Image Correlation. These methods allow different resolutions which can be adjusted according to the expected strain gradient. While the first method depends on the laser wavelength and position of the illumination sources, the second depends on the grating pitch and the last on the surface texture or painted speckle. Algorithms to derivate and filtering the displacement field are developed to compute the strain field and will be used for further works.

36 MATERIALS SCIENCE↗

Deep probabilistic direction prediction in 3D with applications to directional dark matter detectors

Abstract We present the first method to probabilistically predict 3D direction in a deep neural network model. The probabilistic predictions are modeled as a heteroscedastic von Mises-Fisher distribution on the sphere S 2 , giving a simple way to quantify aleatoric uncertainty. This approach generalizes the cosine distance loss which is a special case of our loss function when the uncertainty is assumed to be uniform across samples. We develop approximations required to make the likelihood function and gradient calculations stable. The method is applied to the task of predicting the 3D directions of electrons, the most complex signal in a class of experimental particle physics detectors designed to demonstrate the particle nature of dark matter and study solar neutrinos. Using simulated Monte Carlo data, the initial direction of recoiling electrons is inferred from their tortuous trajectories, as captured by the 3D detectors. For 40 keV electrons in a 70% He 30% CO 2 gas mixture at STP, the new approach achieves a mean cosine distance of 0.104 (26 ∘ ) compared to 0.556 (64 ∘ ) achieved by a non-machine learning algorithm. We show that the model is well-calibrated and accuracy can be increased further by removing samples with high predicted uncertainty. This advancement in probabilistic 3D directional learning could increase the sensitivity of directional dark matter detectors.

Computer Science↗

An improved quasi-diabatic representation of the 1, 2, 3 1 A coupled adiabatic potential energy surfaces of phenol in the full 33 internal coordinates

In a previous project we constructed a quasi-diabatic representation, H d , of the 1, 2, 3 1 A adiabatic states of phenol from high level multireference single and double excitation configuration interaction electronic structure data, energies, energy gradients, and derivative couplings. That H d accurately describes surface minima, saddle points, and also regions of strong nonadiabatic interactions, reproducing the locus of conical intersection seams and the coordinate dependence of the derivative couplings. The present work determines the accuracy of H d for describing phenol photodissociation. Additionally, we demonstrate that a modest energetic shift of two diabats yields a quantifiably more accurate H d compared with experimental energetics. The analysis shows that in favorable circumstances it is possible to use single point energies obtained from the most reliable electronic structure methods available, including methods for which the energy gradients and derivative couplings are not available, to improve the quality of a global representation of several coupled potential energy surfaces. Our data indicate an alternative interpretation of kinetic energy release measurements near λ phot ~ 248 nm.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Hybridizable Discontinuous Galerkin solver for the Grad–Shafranov equation

In axisymmetric fusion reactors, the equilibrium magnetic configuration can be expressed in terms of the solution to a semi-linear elliptic equation known as the Grad–Shafranov equation, the solution of which determines the poloidal component of the magnetic field. When the geometry of the confinement region is known, the problem becomes an interior Dirichlet boundary value problem. We propose a high order solver based on the Hybridizable Discontinuous Galerkin method. The resulting algorithm (1) provides high order of convergence for the flux function and its gradient, (2) incorporates a novel method for handling piecewise smooth geometries by extension from polygonal meshes, (3) can handle geometries with non-smooth boundaries and x-points, and (4) deals with the semi-linearity through an accelerated two-grid fixed-point iteration. The effectiveness of the algorithm is verified with computations for cases where analytic solutions are known on configurations similar to those of actual devices (ITER with single null and double null divertor, NSTX, ASDEX upgrade, and Field Reversed Configurations).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A new hybrid algorithm combining a new chaos optimization approach with gradient descent for high dimensional optimization problems

Chaos optimization algorithm (COA) is a recently developed method for global optimization based on chaos theory. It has the features of easy implementation, short execution time and robust mechanisms of escaping from local optima compared with the existing stochastic searching algorithm. By integrating the gradient descent method with a new COA approach using new strategies including symmetrization and levelling, a novel hybrid method, called SCGD, is proposed in this paper. The main idea is to combine the good ability of exploitation of the gradient descent with the high ability of exploration of COA. The purpose is to avoid being trapped into local optima and to improve convergence in large space and high-dimension optimization problems. The simulation results based on a set of 50 benchmark functions with diverse properties and different dimensions show the efficiency and the abilities of the proposed method in finding the global optima compared with the existing methods.

97 MATHEMATICS AND COMPUTING↗

Calculations of structural, elastic, electronic, magnetic and phonon properties of FeNiMnAl by the first principles

The electronic, elastic and dynamical properties of the quaternary alloy FeNiMnAl have been investigated using a pseudopotential plane wave method within the generalized gradient approximation (GGA). We determined the lattice parameters and the bulk modulus B. In addition, the elastic properties such as elastic constans (C{sub 11}, C{sub 12} and C{sub 44}), the shear modulus G, the young modulus E, the poisson's ratio σ and the B/G ratio are also given. The FeNiMnAl Heusler alloy exhibit a ferromagnetic half-metallic behavior with the total magnetic moment of 4.02 μ{sub B}. The phonon dispersion of FeNiMnAl has been performed using the density functional theory and the direct method with 2×2×2 supercell.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Equation of state of dense plasmas: Orbital-free molecular dynamics as the limit of quantum molecular dynamics for high-Z elements

The applicability of quantum molecular dynamics to the calculation of the equation of state of a dense plasma is limited at high temperature by computational cost. Orbital-free molecular dynamics, based on a semiclassical approximation and possibly on a gradient correction, is a simulation method available at high temperature. For a high-Z element such as lutetium, we examine how orbital-free molecular dynamics applied to the equation of state of a dense plasma can be regarded as the limit of quantum molecular dynamics at high temperature. For the normal mass density and twice the normal mass density, we show that the pressures calculated with the quantum approach converge monotonically towards those calculated with the orbital-free approach; we observe a faster convergence when the orbital-free approach includes the gradient correction. We propose a method to obtain an equation of state reproducing quantum molecular dynamics results up to high temperatures where this approach cannot be directly implemented. With the results already obtained for low-Z plasmas, the present study opens the way for reproducing the quantum molecular dynamics pressure for all elements up to high temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High gradient lens for charged particle beam

Methods and devices enable shaping of a charged particle beam. A dynamically adjustable electric lens includes a series of alternating a series of alternating layers of insulators and conductors with a hollow center. The series of alternating layers when stacked together form a high gradient insulator (HGI) tube to allow propagation of the charged particle beam through the hollow center of the HGI tube. A plurality of transmission lines are connected to a plurality of sections of the HGI tube, and one or more voltage sources are provided to supply an adjustable voltage value to each transmission line of the plurality of transmission lines. By changing the voltage values supplied to each section of the HGI tube, any desired electric field can be established across the HGI tube. This way various functionalities including focusing, defocusing, acceleration, deceleration, intensity modulation and others can be effectuated on a time varying basis.

47 OTHER INSTRUMENTATION↗

Light Scattering by Small Multilayer Particles: A Generalized Separation of Variables Method

A Rayleigh approximation is constructed for light scattering by small multilayer axisymmetric particles, in which their polarizability is determined by the generalized separation of variables method (SVM). In this method, scalar potentials, the gradients of which yield the electric-field strengths, are represented as expansions in spherical harmonics of the Laplace equation. Unknown coefficients of expansions are determined from the boundary conditions, which are reduced to infinite systems of linear algebraic equations (ISLAEs), since the separation of variables is incomplete. The T matrix of the electrostatic problem, principal element T{sub 11} of which is proportional to the particle polarizability, is determined. The necessary condition for the ISLAEs solvability for the SVM coincides with the condition of the correct application of the extended boundary conditions method (EBCM). However, numerical calculations in which finite-dimensional (i.e., reduced) systems are solved, yield different results in ranges of variation of parameters that are close to the boundary of the range of applicability. An analysis of the numerical calculations of the scattering and absorption cross sections for two-layer confocal spheroids, an exact solution for which can be obtained using spheroidal harmonics, shows that the SVM is preferable to the EBCM. It turned out that the proposed method yields workable results in a wider range of variation of parameters. Even outside the range of applicability, in which it should be regarded as a certain approximate solution, its use in a number of cases is quite acceptable. Additional calculations for three-layer nonconfocal spheroids, as well as for three-layer similar pseudospheroids and Pascal’s snails, which can be obtained from spheroids as a result of the inversion with respect to the coordinate origin and one of the foci, respectively, confirm these inferences. We note that, for certain values of the parameters, the shapes of the latter particles are nonconvex.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Angular momentum desaturation for Skylab using gravity gradient torques

An angular momentum desaturation method for momentum exchange devices of orbiting spacecraft is described. The specific application of the method is to the Skylab which contains three double-gimbaled control moment gyros for precise attitude control and maneuvering. It is assumed that the attitude reference is inertially fixed and that two of the vehicle principal moments of inertia are much larger than the third. Gravity gradient torques and resultant angular momentum accumulation are developed for small deviations from the reference. The assumed moment-of-inertia distribution allows desaturation about all axes with only two attitude angles each for the two axes with large moments of inertia. The necessary desaturation maneuvers can be decoupled for a special set of orbital coordinates. All maneuvers are made during the night portion of the orbit, and the percentage utilized for desaturation is selectable. Expressions for the attitude angle commands are developed assuming infinite vehicle rates. The effect of finite rates introduces an efficiency into the desaturation. Expressions for this efficiency are developed and means for compensation are treated. Arbitrary misalignments between geometric vehicle axes and principal moment-of-inertia axes are permissible. An angle bias about the sun line minimizes the angular momentum accumulation about the sun line projection into the orbital plane. Adaptive desaturation maneuver limiting consistent with the available maneuver momentum is included.

Kennel, H. F.↗

WE-AB-BRA-07: Operating Room Quality Assurance (ORQA) for Spine Surgery Using Known-Component 3D-2D Image Registration

Purpose: Intraoperative x-ray radiography/fluoroscopy is commonly used to qualitatively assess delivery of surgical devices (e.g., spine pedicle screws) but can fail to reliably detect suboptimal placement (e.g., breach of adjacent critical structures). We present a method wherein prior knowledge of the patient and surgical components is leveraged to match preoperative CT and intraoperative radiographs for quantitative assessment of 3D pose. The method presents a new means of operating room quantitative quality assurance (ORQA) that could improve quality and safety, and reduce the frequency of revision surgeries. Methods: The algorithm (known-component registration, KC-Reg) uses patient-specific preoperative CT and parametrically defined surgical component models within a robust 3D-2D registration method to iteratively optimize gradient similarity using the covariance matrix adaptation evolution strategy. Advances from previous work address key challenges to clinical translation: i) absolving the need for offline geometric calibration of the C-arm; and ii) solving multiple component bodies simultaneously, thereby allowing QA in a single step (e.g., spinal construct with 4–20 screws), rather than sequential QA of each component. Performance was tested in a spine phantom with 10 pedicle screws, and first results from clinical studies are reported. Results: Phantom experiments demonstrated median target registration error (TRE) of (1.0±0.3) mm at the screw tip and (0.7°±0.4°) in angulation. The simultaneous multi-body registration approach improved TRE from the previous (sequential) method by 42%, reduced outliers, and fits into the natural workflow. Initial application of KC-Reg in clinical data shows TRE of (2.5±4.5) mm and (4.7°±0.5°). Conclusion: The KC-Reg algorithm offers a potentially valuable method for quantitative QA of the surgical product, using radiographic systems that are already within the surgical arsenal. For spine surgery, the method offers a near-real-time independent check on the quality of surgical product, facilitating immediate revision if necessary and potentially avoiding postoperative morbidity and/or revision surgery. Gerhard Kleinszig and Sebastian Vogt are employees of Siemens Healthcare.

60 APPLIED LIFE SCIENCES↗

Research, Development and Fabrication of Lithium Solar Cells

This report continues studies on cells made with low temperature (<400degC) single-cycle, lithium diffusion schedules. Measurements of concentration of lithium at the back surface showed that C(sub S) changes with time consisting of a build-up to a maximum value (above 10(exp 17) cm(exp -3)) with gradual decrease at longer diffusion times. Using the estimated C(sub S) values, theoretical lithium profiles were calculated. The calculated concentration near the PN junction was higher than those measured by the capacitance voltage (C-V) technique. The calculated concentration gradients were one-seventh of those measured by C-V methods. Both the concentration and its gradient were proportional to C(sub S). Estimates of C(sub S) obtained by extrapolating from C-V measured values of concentration and gradient near the PN junction, showed a decrease in C(sub S) for longer diffusion times. This fact supported the trend observed in Shipment C-13, where,, at longer times, the spread in values of V(sub oc), lithium concentration, and gradient became greater. The general conclusion was that the I-V characteristics were well controlled, with good cell output, but for longer diffusion times the lithium concentration gradient could have wider spread. The work described has led to greater insight into the lithium distribution dependence on the surface conditions. The surface conditions could be changed by varying the method used to apply lithium. Vacuum evaporation of lithium shows promise, and another exhaustive test of this method is in progress. Shipment C-13 is summarized, and the succeeding lots (C-14, C-15) are described. It is felt that improved monitoring of cell properties is at hand, and that the control of cell fabrication processes has advanced to where the effects of subtle differences in lithium distribution can now be seen. Evaluation is continuing on the best method to scale up the fabrication steps. A different sequence of fabrication is planned for Shipments C-14 and C-15.

Iles, Peter A.↗

Gradient Flow for Parton Distribution Functions: First Application to the Pion

Parton distribution functions (PDFs) are central to precision QCD phenomenology. Their Mellin moments can be computed on the lattice, but direct determinations using local operators, besides ⟨𝑥⟩, face severe challenges from reduced hypercubic symmetry, limiting results to the lowest moments. A recently proposed method resolves these issues using gradient flow. We demonstrate the efficacy of this method by computing ratios of flavor nonsinglet pion PDF moments up to ⟨𝑥 5 ⟩ on four lattice spacings at 𝑚 𝜋 ≃ 411 MeV. The moments and reconstructed PDF agree quantitatively with recent phenomenological extractions.

Francis, Anthony [National Yang Ming Chiao Tung Un↗

An adaptive discontinuous Petrov-Galerkin method for the Grad-Shafranov equation

In this work, we propose and develop an arbitrary-order adaptive discontinuous Petrov--Galerkin (DPG) method for the nonlinear Grad--Shafranov equation. An ultraweak formulation of the DPG scheme for the equation is given based on a minimal residual method. The DPG scheme has the advantage of providing more accurate gradients compared to conventional finite element methods, which is desired for numerical solutions to the Grad--Shafranov equation. The numerical scheme is augmented with an adaptive mesh refinement approach, and a criterion based on the residual norm in the minimal residual method is developed to achieve dynamic refinement. Nonlinear solvers for the resulting system are explored and a Picard iteration with Anderson acceleration is found to be efficient to solve the system. Finally, the proposed algorithm is implemented in parallel on MFEM using a domain-decomposition approach, and our implementation is general, supporting arbitrary order of accuracy and general meshes. Furthermore, numerical results are presented to demonstrate the efficiency and accuracy of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗

Enhancing Multi-Step Reservoir Inflow Forecasting: A Time-Variant Encoder–Decoder Approach

Accurate reservoir inflow forecasting is vital for effective water resource management. Reliable forecasts enable operators to optimize storage and release strategies to meet competing sectoral demands—such as water supply, irrigation, and hydropower scheduling—while also mitigating flood and drought risks. To address this need, in this study, we propose a novel time-variant encoder–decoder (ED) model designed specifically to improve multi-step reservoir inflow forecasting, enabling accurate predictions of reservoir inflows up to seven days ahead. Unlike conventional ED-LSTM and recursive ED-LSTM models, which use fixed encoder parameters or recursively propagate predictions, our model incorporates an adaptive encoder structure that dynamically adjusts to evolving conditions at each forecast horizon. Additionally, we introduce the Expected Baseline Integrated Gradients (EB-IGs) method for variable importance analysis, enhancing interpretability of inflow by incorporating multiple baselines to capture a broader range of hydrometeorological conditions. The proposed methods are demonstrated at several diverse reservoirs across the United States. Our results show that they outperform traditional methods, particularly at longer lead times, while also offering insights into the key drivers of inflow forecasting. These advancements contribute to enhanced reservoir management through improved forecasting accuracy and practical decision-making insights under complex hydroclimatic conditions.

58 GEOSCIENCES↗

Spectral Gradients and u-v Sampling

I describe a systematic method for characterizing the effects of differences in u-v coverage in terms of inferred spectral gradients in interferometric images. This method is directly applicable to optical and radio interferometry; it is useful in situations when observations using scale arrays are not feasible.

u-v sampling spectral gradients interferometric im↗