Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gradient”

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 19 records

The molecular gradient using the divide-expand-consolidate resolution of the identity second-order Møller-Plesset perturbation theory: The DEC-RI-MP2 gradient

We report an implementation of the molecular gradient using the divide-expand-consolidate resolution of the identity second-order Møller-Plesset perturbation theory (DEC-RI-MP2). The new DEC-RI-MP2 gradient method combines the precision control as well as the linear-scaling and massively parallel features of the DEC scheme with efficient evaluations of the gradient contributions using the RI approximation. We further demonstrate that the DEC-RI-MP2 gradient method is capable of calculating molecular gradients for very large molecular systems. A test set of supramolecular complexes containing up to 158 atoms and 1960 contracted basis functions has been employed to demonstrate the general applicability of the DEC-RI-MP2 method and to analyze the errors of the DEC approximation. Moreover, the test set contains molecules of complicated electronic structures and is thus deliberately chosen to stress test the DEC-RI-MP2 gradient implementation. Additionally, as a showcase example the full molecular gradient for insulin (787 atoms and 7604 contracted basis functions) has been evaluated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dependence of nuclear quadrupole resonance transitions on the electric field gradient asymmetry parameter for nuclides with half-integer spins

Allowed transition energies and eigenstate expansions have been calculated and tabulated in numerical form as functions of the electric field gradient asymmetry parameter for the zero field Hamiltonian of quadrupolar nuclides with I=3/2,5/2,7/2, and 9/2. These results are essential to interpret nuclear quadrupole resonance (NQR) spectra and extract accurate values of the electric field gradient tensors. Applications of NQR methods to studies of electronic structure in heavy element systems are proposed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Secondary Ion Mass Spectrometry Methodology for Isotopic Ratio Measurement of Micro‐Grains in Thin Sections: True Grain Size Estimation and Deconvolution of Inter‐Grain Size Gradients and Intra‐Grain Radial Gradients

In a thin section, grains that were approximately spherical in situ appear circular in cross section, and the distribution of apparent diameters frequently assumed to be their size distribution. Scanning ion imaging by secondary ion mass spectrometry (SIMS) is capable of providing precise (< 1‰) stable isotope ratio measurements of such grains, but, importantly, also registers their rate of evolution in apparent size as they are ablated by the primary beam. By assessing rates of radius change with depth, the described methodology enables the ‘true’ size of grains to be estimated, as well as the distance of the sectioned surface from the original grain centre. Transects in three dimensions are made possible, and this capability enables better identification (and thus separation) of both inter‐grain chemical signatures as a function of grain size, and intra‐grain radial trends. In this example, we highlight the specific application to pyrite (FeS 2 ) minerals, which are frequently analysed by SIMS to determine their inter‐grain and intra‐grain geochemical variations, particularly in their sulfur stable isotopic ratios (δ 34 S). Benefits of the new methodology over the Faraday cup ‘spot mode’ are described. Data correction algorithms and precision considerations are discussed.

Jones, Clive↗

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix

We present an iterative algorithm for computing an invariant subspace associated with the algebraically smallest eigenvalues of a large sparse or structured Hermitian matrix A. We are interested in the case in which the dimension of the invariant subspace is large (e.g., over several hundreds or thousands) even though it may still be small relative to the dimension of A. These problems arise from, for example, density functional theory (DFT) based electronic structure calculations for complex materials. The key feature of our algorithm is that it performs fewer Rayleigh–Ritz calculations compared to existing algorithms such as the locally optimal block preconditioned conjugate gradient or the Davidson algorithm. It is a block algorithm, and hence can take advantage of efficient BLAS3 operations and be implemented with multiple levels of concurrency. We discuss a number of practical issues that must be addressed in order to implement the algorithm efficiently on a high performance computer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Communication: slow relaxation, spatial mobility gradients, and vitrification in confined films

Two decades of experimental research indicate that spatial confinement of glass-forming molecular and polymeric liquids results in major changes of their slow dynamics beginning at large confinement distances. A fundamental understanding remains elusive given the generic complexity of activated relaxation in supercooled liquids and the major complications of geometric confinement, interfacial effects, and spatial inhomogeneity. For this research, we construct a predictive, quantitative, force-level theory of relaxation in free-standing films for the central question of the nature of the spatial mobility gradient. The key new idea is that vapor interfaces speed up barrier hopping in two distinct, but coupled, ways by reducing near surface local caging constraints and spatially long range collective elastic distortion. Effective vitrification temperatures, dynamic length scales, and mobile layer thicknesses naturally follow. In conclusion, our results provide a unified basis for central observations of dynamic and pseudo-thermodynamic measurements.

36 MATERIALS SCIENCE↗

Theory of activated glassy relaxation, mobility gradients, surface diffusion, and vitrification in free standing thin films

We have constructed a quantitative, force level, statistical mechanical theory for how confinement in free standing thin films introduces a spatial mobility gradient of the alpha relaxation time as a function of temperature, film thickness, and location in the film. The crucial idea is that relaxation speeds up due to the reduction of both near-surface barriers associated with the loss of neighbors in the local cage and the spatial cutoff and dynamical softening near the vapor interface of the spatially longer range collective elasticity cost for large amplitude hopping. These two effects are fundamentally coupled. Quantitative predictions are made for how an apparent glass temperature depends on the film thickness and experimental probe technique, the emergence of a two-step decay and mobile layers in time domain measurements, signatures of confinement in frequency-domain dielectric loss experiments, the dependence of film-averaged relaxation times and dynamic fragility on temperature and film thickness, surface diffusion, and the relationship between kinetic experiments and pseudo-thermodynamic measurements such as ellipsometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Testing thermal gradient driving force for grain boundary migration using molecular dynamics simulations

Strong thermal gradients in low-thermal-conductivity ceramics may drive extended defects, such as grain boundaries and voids, to migrate in preferential directions. In this work, molecular dynamics simulations are conducted to study thermal gradient driven grain boundary migration and to verify a previously proposed thermal gradient driving force equation, using uranium dioxide as a model system. It is found that a thermal gradient drives grain boundaries to migrate up the gradient and the migration velocity increases under a constant gradient owing to the increase in mobility with temperature. Different grain boundaries migrate at very different rates due to their different intrinsic mobilities. The extracted mobilities from the thermal gradient driven simulations are compared with those calculated from two other well-established methods and good agreement between the three different methods is found, demonstrating that the theoretical equation of the thermal gradient driving force is valid, although a correction of one input parameter should be made. The discrepancy in the grain boundary mobilities between modeling and experiments is also discussed.

36 MATERIALS SCIENCE↗

Field-gradient partitioning for fracture and frictional contact in the material point method: Field-gradient partitioning for fracture and frictional contact in the material point method [Fracture and frictional contact in material point method using damage-field gradients for velocity-field partitioning]

Contact and fracture in the material point method require grid-scale enrichment or partitioning of material into distinct velocity fields to allow for displacement or velocity discontinuities at a material interface. We present a new method which a kernel-based damage field is constructed from the particle data. The gradient of this field is used to dynamically repartition the material into contact pairs at each node. Our approach avoids the need to construct and evolve explicit cracks or contact surfaces and is therefore well suited to problems involving complex 3-D fracture with crack branching and coalescence. A straightforward extension of this approach permits frictional ‘self-contact’ between surfaces that are initially part of a single velocity field, enabling more accurate simulation of granular flow, porous compaction, fragmentation, and comminution of brittle materials. Finally, numerical simulations of self contact and dynamic crack propagation are presented to demonstrate the accuracy of the approach.

97 MATHEMATICS AND COMPUTING↗

Stabilization of electron-scale turbulence by electron density gradient in national spherical torus experiment

Theory and experiments have shown that electron temperature gradient (ETG) turbulence on the electron gyro-scale, k{sub ⊥}ρ{sub e} ≲ 1, can be responsible for anomalous electron thermal transport in NSTX. Electron scale (high-k) turbulence is diagnosed in NSTX with a high-k microwave scattering system [D. R. Smith et al., Rev. Sci. Instrum. 79, 123501 (2008)]. Here we report on stabilization effects of the electron density gradient on electron-scale density fluctuations in a set of neutral beam injection heated H-mode plasmas. We found that the absence of high-k density fluctuations from measurements is correlated with large equilibrium density gradient, which is shown to be consistent with linear stabilization of ETG modes due to the density gradient using the analytical ETG linear threshold in F. Jenko et al. [Phys. Plasmas 8, 4096 (2001)] and linear gyrokinetic simulations with GS2 [M. Kotschenreuther et al., Comput. Phys. Commun. 88, 128 (1995)]. We also found that the observed power of electron-scale turbulence (when it exists) is anti-correlated with the equilibrium density gradient, suggesting density gradient as a nonlinear stabilizing mechanism. Higher density gradients give rise to lower values of the plasma frame frequency, calculated based on the Doppler shift of the measured density fluctuations. Linear gyrokinetic simulations show that higher values of the electron density gradient reduce the value of the real frequency, in agreement with experimental observation. Nonlinear electron-scale gyrokinetic simulations show that high electron density gradient reduces electron heat flux and stiffness, and increases the ETG nonlinear threshold, consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrated Design of Downwind Land-Based Wind Turbines using Analytic Gradients

Wind turbines are complex systems where component-level changes can have significant system-level effects. Effective wind turbine optimization generally requires an integrated analysis approach with a large number of design variables. Optimizing across large variable sets is orders of magnitude more efficient with gradient-based methods as compared with gradient-free method, particularly when using exact gradients. We have developed a wind turbine analysis set of over 100 components where 90% of the models provide numerically exact gradients through symbolic differentiation, automatic differentiation, and adjoint methods. This framework is applied to a specific design study focused on downwind land-based wind turbines. Downwind machines are of potential interest for large wind turbines where the blades are often constrained by the stiffness required to prevent a tower strike. The mass of these rotor blades may be reduced by utilizing a downwind configuration where the constraints on tower strike are less restrictive. The large turbines of this study range in power rating from 5-7MW and in diameter from 105m to 175m. The changes in blade mass and power production have important effects on the rest of the system, and thus the nacelle and tower systems are also optimized. For high-speed wind sites, downwind configurations do not appear advantageous. The decrease in blade mass (10%) is offset by increases in tower mass caused by the bending moment from the rotor-nacelle-assembly. For low-wind speed sites, the decrease in blade mass is more significant (25-30%) and shows potential for modest decreases in overall cost of energy (around 1-2%).

17 WIND ENERGY↗

Up-gradient particle flux in a drift wave-zonal flow system

We report a net inward, up-gradient turbulent particle flux in a cylindrical plasma when collisional drift waves generate a sufficiently strong sheared azimuthal flow that drives positive (negative) density fluctuations up (down) the background density gradient, resulting in a steepening of the mean density gradient. The results show the existence of a saturation mechanism for drift-turbulence driven sheared flows that can cause up-gradient particle transport and density profile steepening.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of Wind Farm Layout Optimization Results Using a Simple Wake Model and Gradient-Based Optimization to Large-Eddy Simulations: Preprint

Comparison of wind farm layout optimization results with simulations or wind tunnel tests is important because the models used during optimization use simplifying assumptions that can alter the design space. Some of the characteristics of the simple models may negatively influence the resulting layouts. In this paper we perform wind farm layout optimization and compare the resulting improvements to large eddy simulation (LES) results. We begin by describing the models used, including changes specific for use with gradient-based optimization. We then compare our version of the models with the original versions and large eddy simulation LES) data from previous publications. Using the models described, we performed gradient-based wind farm layout optimization, using exact gradients, and compare the results to LES simulations. The models predicted an improvement in annual energy production (AEP) of 7.4%, while the LES reported an AEP improvement of 9:9%. We conclude that the improvements found by optimizing with the simple models are not just an artifact of the models, but are real improvements.

17 WIND ENERGY↗

Accelerating gradient descent and Adam via fractional gradients

Here we propose a class of novel fractional-order optimization algorithms. We define a fractional-order gradient via the Caputo fractional derivatives that generalizes integer-order gradient. We refer it to as the Caputo fractional-based gradient, and develop an efficient implementation to compute it. A general class of fractional-order optimization methods is then obtained by replacing integer-order gradients with the Caputo fractional-based gradients. To give concrete algorithms, we consider gradient descent (GD) and Adam, and extend them to the Caputo fractional GD (CfGD) and the Caputo fractional Adam (CfAdam). We demonstrate the superiority of CfGD and CfAdam on several large scale optimization problems that arise from scientific machine learning applications, such as ill-conditioned least squares problem on real-world data and the training of neural networks involving non-convex objective functions. Numerical examples show that both CfGD and CfAdam result in acceleration over GD and Adam, respectively. We also derive error bounds of CfGD for quadratic functions, which further indicate that CfGD could mitigate the dependence on the condition number in the rate of convergence and results in significant acceleration over GD.

97 MATHEMATICS AND COMPUTING↗