Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Expanded Equations for Torque and Force on a Cylindrical Permanent Magnet Core in a Large-Gap Magnetic Suspension System

The expanded equations for torque and force on a cylindrical permanent magnet core in a large-gap magnetic suspension system are presented. The core is assumed to be uniformly magnetized, and equations are developed for two orientations of the magnetization vector. One orientation is parallel to the axis of symmetry, and the other is perpendicular to this axis. Fields and gradients produced by suspension system electromagnets are assumed to be calculated at a point in inertial space which coincides with the origin of the core axis system in its initial alignment. Fields at a given point in the core are defined by expanding the fields produced at the origin as a Taylor series. The assumption is made that the fields can be adequately defined by expansion up to second-order terms. Examination of the expanded equations for the case where the magnetization vector is perpendicular to the axis of symmetry reveals that some of the second-order gradient terms provide a method of generating torque about the axis of magnetization and therefore provide the ability to produce six-degree-of-freedom control.

Groom, Nelson J.↗
DOE OSTI2019

Introduction of a bounded penalty function in contact‐assisted simulations of protein structures to omit false restraints

Contact‐assisted simulations, the contacts being predicted or determined experimentally, have become very important in the determination of the structures of proteins and other biological macromolecules. In this work, the effect of contact‐distance restraints on the simulated structures was investigated with the use of multiplexed replica exchange simulations with the coarse‐grained UNRES force field. A modified bounded flat‐bottom restraint function that does not generate a gradient when a restraint cannot be satisfied was implemented. Calculations were run with (i) a set of four small proteins, with contact restraints derived from experimental structures, and (ii) selected CASP11 and CASP12 targets, with restraints as used at prediction time. The bounded penalty function largely omitted false contacts, which were usually inconsistent. It was found that at least 20% of correct contacts must be present in the restraint set to improve model quality with respect to unrestrained simulations. © 2019 Wiley Periodicals, Inc.

Lubecka, Emilia A.↗
DOE OSTI2019

Defect Characterization Through Automated Laser Track Trace Identification in SLM Processes Using Laser Profilometer Data

This paper presents recent developments in statistical image processing to demonstrate feasibility of layer-by-layer nondestructive inspection for selective laser melting (SLM) of metallic parts. A matrix of 10 mm × 10 mm stainless 316 squares was deposited with the first row consisting of just the sintered surface, Row 2 a single layer, Row 3 with 2 layers and so on. These layers were scanned to emulate the layer-by-layer collection of top surface geometry using a laser sensor. The resultant data were utilized to generate ISO 25718 roughness parameters and subsequently to perform track identification. This calculation consisted of two salient steps: (1) computing the gradients of surface roughness, which vanish at the peaks and pits, and (2) extracting the scanning direction by taking Hough transform of the geometry. This algorithm was repeated for all layers, and the alignment of the roughness with the rotating scan was observed to be persistent for all layers. Correlation between surface properties obtained along various scan directions and defect probability is explored.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Plant organic matter inputs exert a strong control on soil organic matter decomposition in a thawing permafrost peatland

Peatlands are climate critical carbon (C) reservoirs that could become a C source under continued warming. A strong relationship between plant tissue chemistry and the soil organic matter (SOM) that fuels C gas emissions is inferred, but rarely examined at the molecular level. Here we compared Fourier transform infrared (FT-IR) spectroscopy measurements of solid phase functionalities in plants and SOM to ultra-high-resolution mass spectrometric analyses of plant and SOM water extracts across a palsa-bog-fen thaw and moisture gradient in an Arctic peatland. From these analyses we calculated the C oxidation state (NOSC), a measure which can be used to assess organic matter quality. Palsa plant extracts had the highest NOSC, indicating high quality, whereas extracts of Sphagnum, which dominated the bog, had the lowest NOSC. The percentage of plant compounds that are less bioavailable and accumulate in the peat, increases from palsa (25%) to fen (41%) to bog (47%), reflecting the pattern of percent Sphagnum cover. The pattern of NOSC in the plant extracts was consistent with the high number of consumed compounds in the palsa and low number of consumed compounds in the bog. However, in the FT-IR analysis of the solid phase bog peat, carbohydrate content was high implying high quality SOM. Here, we explain this discrepancy as the result of low solubilization of bog SOM facilitated by the low pH in the bog which makes the solid phase carbohydrates less available to microbial decomposition. Plant-associated condensed aromatics, tannins, and lignin-like compounds declined in the unsaturated palsa peat indicating decomposition, but lignin-like compounds accumulated in the bog and fen peat where decomposition was presumably inhibited by the anaerobic conditions. A molecular-level comparison of the aboveground C sources and peat SOM demonstrates that climate-associated vegetation shifts in peatlands are important controls on the mechanisms underlying changing C gas emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Unravelling Photocarrier Dynamics beyond the Space Charge Region for Photoelectrochemical Water Splitting

Semiconductor photoelectrodes for photoelectrochemical (PEC) water splitting require efficient carrier generation, separation, and transport at and beyond the space charge region (SCR) formed at the aqueous interface. The trade-off between photon collection and minority carrier delivery governs the photoelectrode design and implies maximum water splitting efficiency at an electrode thickness equivalent to the light absorption depth. Here, using planar ZnO thin films as a model system, we identify the photocarriers beyond the SCR as another significant source to substantially enhance the PEC performance. The high-quality ZnO films synthesized by pulsed laser deposition feature very few deep trap states and support a long photocarrier lifetime. Combined with photoelectrochemical characterization, ultrafast spectroscopy, and numerical calculations, it is revealed that engineering the exciton concentration gradient by film thickness facilitates the inward diffusion of photocarriers from the neighboring illuminated region to the SCR and, therefore, achieves a record high quantum efficiency over 80% at a thickness far beyond its light absorption depth and the SCR width. Furthermore, these results elucidate the important role of the photocarriers beyond SCR for the PEC process and provide new insight into exploring the full potential for efficient photoelectrode materials with large exciton diffusivity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction

In this study, we have extended our graphical processing unit (GPU)-accelerated direct configuration interaction program to multiple devices, reducing iteration times for configuration spaces of 165 million determinants to only 3 s using NVIDIA P100 GPUs. Similar improvements in the one- and two-particle reduced density matrix formation allow for fast analytical energy gradients and electronic properties. Our parallel algorithm enables the calculation of arbitrarily large configuration spaces (limited only by available system memory), with iteration times of 13 min for an active space of 18 electrons in 18 orbitals (2.4 billion determinants) using six consumer grade NVIDIA 1080Ti GPUs. These advances enable routine molecular dynamics simulations, geometry optimizations, and absorption spectrum calculations for molecules with large configuration spaces, a task that has heretofore required massive computational effort. In this work, we demonstrate the utility of our program by generating the absorption spectrum for diphenyl acetylene at the floating occupation molecular orbital complete active space configuration interaction level of theory. Lastly, several active spaces were investigated to assess the dependence of spectral features on orbital space dimension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Using Ensemble Data Assimilation to Estimate Transient Hydrologic Exchange Flow Under Highly Dynamic Flow Conditions

Abstract Quantifying dynamic hydrologic exchange flows (HEFs) within river corridors that experience high‐frequency flow variations caused by dam regulations is important for understanding the biogeochemical processes at the river water and groundwater interfaces. Heat has been widely used as a tracer to infer steady‐state flow velocities through analytical solutions of heat transport defined by the diurnal temperature signals. Under sub‐daily dynamic flow conditions, however, such analytical solutions are not applicable due to the violation of their fundamental assumptions. In this study, we developed a data assimilation‐based approach to estimate the sub‐daily flux under highly dynamic flow conditions using multi‐depth temperature observations at a 5‐min resolution. If the hydraulic gradient is measured, Darcy's law was used to calculate the flux with permeability estimated from temperature responses below the riverbed. Otherwise, flux was estimated directly by assimilating multi‐depth temperature data at 1‐ or 2‐hr time intervals assuming one‐dimensional flow and heat transport governing equation. By comparing estimated fluxes with model‐generated synthetic truth, we demonstrated that both schemes have robust performance in estimating fluxes under highly dynamic flow conditions. This data assimilation‐based flux estimation method was able to capture the vertical sub‐daily fluxes using multi‐depth high‐resolution temperature data alone, even in the presence of multi‐dimensional flow. This approach has been successfully applied to real field temperature data collected at the Hanford site, which experiences highly dynamic HEFs. Our study shows the promise of adopting distributed 1‐D temperature monitoring to capture spatial and temporal exchange dynamics in river corridors at a watershed scale or beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Computationally predicted energies and properties of defects in GaN

Recent developments in theoretical techniques have significantly improved the predictive power of density-functional-based calculations. In this article, we discuss how such advancements have enabled improved understanding of native point defects in GaN. We review the methodologies for the calculation of point defects, and discuss how techniques for overcoming the band-gap problem of density functional theory affect native defect calculations. In particular, we examine to what extent calculations performed with semilocal functionals (such as the generalized gradient approximation), combined with correction schemes, can produce accurate results. The properties of vacancy, interstitial, and antisite defects in GaN are described, as well as their interaction with common impurities. We also connect the first-principles results to experimental observations, and discuss how native defects and their complexes impact the performance of nitride devices. Overall, we find that lower-cost functionals, such as the generalized gradient approximation, combined with band-edge correction schemes can produce results that are qualitatively correct. However, important physics may be missed in some important cases, particularly for optical transitions and when carrier localization occurs.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Oscillator strengths, first-order properties, and nuclear gradients for local ADC(2)

We describe theory and implementation of oscillator strengths, orbital-relaxed first-order properties, and nuclear gradients for the local algebraic diagrammatic construction scheme through second order. The formalism is derived via time-dependent linear response theory based on a second-order unitary coupled cluster model. The implementation presented here is a modification of our previously developed algorithms for Laplace transform based local time-dependent coupled cluster linear response (CC2LR); the local approximations thus are state specific and adaptive. The symmetry of the Jacobian leads to considerable simplifications relative to the local CC2LR method; as a result, a gradient evaluation is about four times less expensive. Test calculations show that in geometry optimizations, usually very similar geometries are obtained as with the local CC2LR method (provided that a second-order method is applicable). As an exemplary application, we performed geometry optimizations on the low-lying singlet states of chlorophyllide a.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Accurate lattice geometrical parameters and bulk moduli from a semilocal density functional

Accurate prediction of lattice constants is very important in applications of density functional theory. In this work, we assess the efficacy of a non-empirical meta-generalized gradient approximation proposed by Tao and Mo (TM) by calculating the lattice constants as well as bulk moduli of 33 crystalline semiconductors within the TM scheme. We find that the TM functional is able to produce very accurate lattice constants, with a mean absolute error of 0.038 Å, and bulk moduli with a mean absolute error of 3.2 GPa, improving upon commonly-used semilocal density functionals, such as the LSDA, PBE, SOGGA, PBEsol, TPSS, M06L, and SCAN. The high computational efficiency and remarkable agreements with the corresponding experimental values suggest that the TM functional can be a very competitive candidate in electronic structure theory. We attribute the accuracy of the TM functional to be the result of its satisfaction of many exact or nearly-exact conditions related to the exchange-correlation energy and the associated hole, leading to an improved description of the short- as well as intermediate-range van der Waals interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A unified theory of zonal flow shears and density corrugations in drift wave turbulence

A unified theory of zonal flow shears and density corrugations in drift wave turbulence is presented. Polarization and density advection beat excitation are studied in combination with modulational response. Noise is driven by two-time flux correlation. While the effective zonal flow eddy viscosity can go negative, the zonal diffusivity is positive definite. There is no inverse cascade of density corrugation. The connection between avalanches and corrugations is discussed. Here, the zonal cross-correlation is identified and calculated. Conditions for alignment of zonal shears and corrugation gradients are determined, and the implications for staircase structure are discussed. We show that the synergy of beat noise and modulational effects is stronger than either alone. Strong zonal flows can be excited well below the modulational instability threshold. In the context of L–H transition, zonal noise quenches turbulence overshoot by eliminating the threshold for zonal flow excitation. The power threshold for L–H transition is lowered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Real-Space Multiple-Scattering Theory and Its Applications at Exascale

In recent decades, the ab initio methods based on density functional theory (DFT) (Hohenberg and Kohn 1964, Kohn and Sham 1965) have become a widely used tool in computational materials science, which allows theoretical prediction of physical properties of materials from the first principles and theoretical interpretation of new physical phenomena found in experiments. In the framework of DFT, the original problem that requires solving a quantum mechanical equation for a many-electron system is reduced to a one-electron problem that involves an electron moving in an effective field, while the effective field potential is made up of an electrostatic potential, also known as Hartree potential, arising from the electronic and ion charge distribution in space and an exchange–correlation potential, which is a function of the electron density and encapsulates the exchange and correlation effects of the many-electron system. Even though the exact functional form of the exchange-correlation potential is formally unknown, a local density approximation (LDA) or a generalized gradient approximation (GGA) is usually applied so that the calculation of the exchange–correlation potential, as well as the exchange–correlation energy, becomes tractable while a required accuracy is retained. Based on DFT, ab initio electronic structure calculations for a material generally involve a self-consistent process that iterates between two computational tasks: (1) solving an one-electron Schrödinger equation, also known as Kohn–Sham equation, to obtain the electron density and, if needed, the magnetic moment density, and (2) solving the Poisson equation to obtain the electrostatic potential corresponding to the electron density and constructing the effective potential by adding the exchange–correlation potential to the electrostatic potential. This self-consistent process proceeds until a convergence criteria is reached.

Eisenbach, Markus↗
DOE OSTI2020

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Reviews and syntheses: Australian vegetation phenology: new insights from satellite remote sensing and digital repeat photography

Phenology is the study of periodic biological occurrences and can provide important insights into the influence of climatic variability and change on ecosystems. Understanding Australia's vegetation phenology is a challenge due to its diverse range of ecosystems, from savannas and tropical rainforests to temperate eucalypt woodlands, semi-arid scrublands, and alpine grasslands. These ecosystems exhibit marked differences in seasonal patterns of canopy development and plant life-cycle events, much of which deviates from the predictable seasonal phenological pulse of temperate deciduous and boreal biomes. Many Australian ecosystems are subject to irregular events (i.e. drought, flooding, cyclones, and fire) that can alter ecosystem composition, structure, and functioning just as much as seasonal change. We show how satellite remote sensing and ground-based digital repeat photography (i.e. phenocams) can be used to improve understanding of phenology in Australian ecosystems. First, we examine temporal variation in phenology on the continental scale using the enhanced vegetation index (EVI), calculated from MODerate resolution Imaging Spectroradiometer (MODIS) data. Spatial gradients are revealed, ranging from regions with pronounced seasonality in canopy development (i.e. tropical savannas) to regions where seasonal variation is minimal (i.e. tropical rainforests) or high but irregular (i.e. arid ecosystems). Next, we use time series colour information extracted from phenocam imagery to illustrate a range of phenological signals in four contrasting Australian ecosystems. These include greening and senescing events in tropical savannas and temperate eucalypt understorey, as well as strong seasonal dynamics of individual trees in a seemingly static evergreen rainforest. We also demonstrate how phenology links with ecosystem gross primary productivity (from eddy covariance) and discuss why these processes are linked in some ecosystems but not others. We conclude that phenocams have the potential to greatly improve the current understanding of Australian ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1973

A method for calculating unsteady turbulent boundary layers in two- and three-dimensional flows.

The governing equations for an unsteady turbulent boundary layer on a swept infinite cylinder, composed of a continuity equation, a pair of momentum equations and a pair of rate equations for the shear stress, based on the turbulent kinetic-energy equation, are solved numerically. These rate equations for the shear-stress express a balance between the convection, production, dissipation and diffusion of the turbulent shear stress components. Due to the physical model of turbulence employed, this system of equations is hyperbolic. An explicit, second-order accurate, conditionally stable difference scheme for this system, together with appropriate boundary and initial conditions, is formulated and developed into a computer program. Calculations for oscillating free stream flows with no pressure gradient show significant unsteady effects on the turbulent boundary layer.

Singleton, R. E.↗
NASA NTRS1974

Diffusion with chemical reaction: An attempt to explain number density anomalies in experiments involving alkali vapor

The mutual diffusion of two reacting gases is examined which takes place in a bath of inert gas atoms. Solutions are obtained between concentric spheres, each sphere acting as a source for one of the reactants. The calculational model is used to illustrate severe number density gradients observed in absorption experiments with alkali vapor. Severe gradients result when sq root k/D R is approximately 5 where k, D, and R are respectively the second order rate constant, the multicomponent diffusion constant, and the geometrical dimension of the experiment.

Snow, W. L.↗