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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Methodology to Generate the Decay Gamma Source for a Second Target Station Target Wedge

This report details the specifics of applying the Position-Averaged Method methodology for the calculation of the rotating lasagna-style target wedge activation and the resulting decay gamma source terms. The methodology uses MCNP® Code Version 6.2.0 with the RNUCS patch coupled with CINDER2008 from the AARE V1.0 package. The process begins with model preparations that include ensuring impurities are included in the materials in the regions of interest along with segmentation of the geometry in order to capture the decay gamma source gradients. This report also details the physics models and other MCNP®, CINDER2008, and ADVANTG 3.2.0 related options used to calculate the decay gamma source terms. ADVANTG 3.2.0 is used to generate a weight window set to aid in convergence of the neutron fluxes in the initial MCNP® protons-on-target particle transport calculation to calculate the neutron fluxes and spallation products in the regions of interest. The neutron fluxes and spallation products are then used in the CINDER2008 transmutation calculation, and the decay gamma sources are generated by the AARE_GAMMA_SOURCE_SCRIPT as a part of the AARE V1.0 package. These decay gamma sources include the energy distributions as well as the necessary distributions to sample the location of the source using the cell rejection method. The gamma sources that are produced using the methodology may be used for a wide variety of analysis involving the target wedge. This report includes the material specifications, physics options, proton source term description, geometrical configuration, and post-processing details.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Three-dimensional gyrokinetic simulation of the relaxation of a magnetized temperature filament

An electromagnetic, 3D gyrokinetic particle code is used to study the relaxation of a magnetized electron temperature filament embedded in a large, uniform plasma of lower temperature. The study provides insight into the role played by unstable drift-Alfvén waves observed in a basic electron heat transport experiment [D. C. Pace et al., Phys. Plasmas 15, 122304 (2008)] in which anomalous cross-field transport has been documented. The simulation exhibits the early growth of temperature-gradient-driven, drift-Alfvén fluctuations that closely match the eigenmodes predicted by linear theory. At the onset of saturation, the unstable fluctuations display a spiral spatial pattern, similar to that observed in the laboratory, which causes the rearrangement of the temperature profile. After saturation of the linear instability, the system exhibits a markedly different behavior depending on the inclusion in the computation of modes without variation along the magnetic field, i.e., k{sub z} = 0. In their absence, the initial filament evolves into a broadened temperature profile, self-consistent with undamped, finite amplitude drift-Alfvén waves. But the inclusion of k{sub z} = 0 modes causes the destruction of the filament and damping of the drift-Alfvén modes leading to a final state consisting of undamped convective cells and multiple, smaller-scale filaments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Self-consistent implementation of meta-GGA functionals for the ONETEP linear-scaling electronic structure package

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Lastly, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms.

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

Insight into Structural, Electronic, Magnetic, Mechanical, and Thermodynamic Properties of Actinide Perovskite BaPuO{sub 3}

In this paper, we have performed a systematic investigation on structural, magnetic, electronic, mechanical, and thermodynamic properties of BaPuO{sub 3} perovskite within density functional theory (DFT) using generalized gradient approximations (GGA), onsite coulomb repulsion (GGA + U), and modified Becke-Johnson (mBJ). The calculated structural parameters were found in good agreement with the experimental results. A large value of magnetic moment equal to integer value of 4 μ{sub B} was obtained for the compound. The spin-polarized electronic band structure and density of states present 100% of spin polarization at Fermi level, resulting in half-metallic nature for the compound with spin-up states as metallic and spin- down states as a semiconducting. The elastic and mechanical properties have also been predicted. Moreover, we have calculated thermodynamic properties like Debye temperature (𝜃{sub D}), specific heat (C{sub V}), entropy (S), etc. using quasi-harmonic Debye model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Periodic boundary conditions for arbitrary deformations in molecular dynamics simulations

A generalization of the Lees-Edwards periodic boundary conditions (gLE-PBC) for molecular dynamics (MD) simulations is developed to allow for arbitrary deformations to be applied to the domain. The gLE-PBC domain remains a rectangular cuboid regardless of the applied deformation in contrast with the Lagrangian-rhomboid periodic boundary conditions (LR-PBC) where the domain deforms according to the applied deformation. Furthermore, the kinematics of gLE-PBC are validated against pure shear. The gLE-PBC method for interacting systems is then validated against the LR-PBC method and analytical solutions for a solid under isotropic compression, one-dimensional shearing, three-dimensional extension and shearing and for a liquid under Couette flow. Bulk physical properties extracted from the gLE-PBC simulations agree well with values calculated from equilibrium MD simulations. Three dimensional shearing and a deformation with a full velocity gradient matrix are also simulated, showing the range of problems gLE-PBC can explore.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Relationship between Conductivity, Ion Diffusion, and Transference Number in Perfluoropolyether Electrolytes

Connecting continuum-scale ion transport properties such as conductivity and cation transference number to microscopic transport properties such as ion dissociation and ion self-diffusivities is an unresolved challenge in characterizing polymer electrolytes. Better understanding of the relationship between microscopic and continuum scale transport properties would enable the rational design of improved electrolytes for applications such as lithium batteries. Here, we present measurements of continuum and microscopic ion transport properties of nonflammable liquid electrolytes consisting of binary mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and perfluoropolyethers (PFPE) with different end groups: diol, dimethyl carbonate, ethoxy-diol, and ethoxy-dimethyl carbonate. The continuum properties, conductivity and cation transference number, were measured by ac impedance spectroscopy and potentiostatic polarization, respectively. The ion self-diffusivities were measured by pulsed field gradient nuclear magnetic resonance spectroscopy (PFG-NMR), and a microscopic cation transference number was calculated from these measurements. The measured ion self-diffusivities did not reflect the measured conductivities; in some cases, samples with high diffusivities exhibited low conductivity. We introduce a nondimensional parameter, β, that combines microscopic diffusivities and conductivity. We show that β is a sensitive function of end-group chemistry. In the ethoxylated electrolytes, β is close to unity, the value expected for electrolytes that obey the Nernst-Einstein equation. In these cases, the microscopic and continuum transference numbers are in reasonable agreement. PFPE electrolytes devoid of ethoxy groups exhibit values of β that are significantly lower than unity. In these cases, there is significant deviation between microscopic and continuum transference numbers. We propose that this may be due to electrostatic coupling of the cation and anion or contributions to the NMR signal from neutral ion pairs.

25 ENERGY STORAGE↗
DOE OSTI2020

Implications of the fractional charge of hydroxide at the electrochemical interface

Rational design of materials that efficiently convert electrical energy into chemical bonds will ultimately depend on a thorough understanding of the electrochemical interface at the atomic level. Towards this goal, the use of density functional theory (DFT) at the generalized gradient approximation (GGA) level has been applied widely in the past 15 years. In the calculation of electrochemical reaction energetics using GGA-DFT, it is frequently implicitly assumed that ions in the Helmholtz plane have unit charge. However, the ion charge is observed to be fractional near the interface through both a capacitor model and through Bader charge partitioning. In this work, we show that this spurious charge transfer can be effectively mitigated by continuum charging of the electrolyte. Finally, we show that, similar to hydronium, the observed fractional charge of hydroxide is not due to a GGA level self-interaction error, as the partial charge is observed even when using hybrid level exchange–correlation functionals.

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

Fermi level pinning induced electrostatic fields and band bending at organic heterojunctions

The energy level alignment at interfaces between organic semiconductors is of direct relevance to understand charge carrier generation and recombination in organic electronic devices. Commonly, work function changes observed upon interface formation are interpreted as interface dipoles. In this study, using ultraviolet and X-ray photoelectron spectroscopy, complemented by electrostatic calculations, we find a huge work function decrease of up to 1.4 eV at the C{sub 60} (bottom layer)/zinc phthalocyanine (ZnPc, top layer) interface prepared on a molybdenum trioxide (MoO{sub 3}) substrate. However, detailed measurements of the energy level shifts and electrostatic calculations reveal that no interface dipole occurs. Instead, upon ZnPc deposition, a linear electrostatic potential gradient is generated across the C{sub 60} layer due to Fermi level pinning of ZnPc on the high work function C{sub 60}/MoO{sub 3} substrate, and associated band-bending within the ZnPc layer. This finding is generally of importance for understanding organic heterojunctions when Fermi level pinning is involved, as induced electrostatic fields alter the energy level alignment significantly.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

On the electronic, structural, and thermodynamic properties of Au supported on α-Fe{sub 2}O{sub 3} surfaces and their interaction with CO

Extensive first principles calculations are carried out to investigate Au monomers and dimers supported on α-Fe{sub 2}O{sub 3}(0001) surfaces in terms of structure optimizations, electronic structure analyses, and ab initio thermodynamics calculations of surface phase diagrams. All computations rely on density functional theory in the generalized gradient approximation (Perdew-Burke-Ernzerhof (PBE)) and account for on-site Coulomb interactions via inclusion of a Hubbard correction (PBE+U). The relative stability of Au monomers/dimers on the stoichiometric termination of α-Fe{sub 2}O{sub 3}(0001) decorated with various vacancies (multiple oxygen vacancies, iron vacancy, and mixed iron-oxygen vacancies) has been computed as a function of the oxygen chemical potential. The charge rearrangement induced by Au at the oxide contact is analyzed in detail and discussed. On one hand, ab initio thermodynamics predicts that under O-rich conditions, structures obtained by replacing a surface Fe atom with a Au atom are thermodynamically stable over a wide range of temperatures. On the other hand, the complex of a CO molecule on a Au atom substituting surface Fe atoms is thermodynamically stable only in a much more narrow range of values of the O chemical potential under O-rich conditions. In the case of a Au dimer, under O-rich conditions, supported Au atoms at an O-Fe di-vacancy are more stable. However, upon CO adsorption, the complex of a CO molecule and 2 Au atoms located at a single Fe vacancy is more favorable.

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

Relationship between thermoelectric figure of merit and energy conversion efficiency

Significance Thermoelectric materials generate electricity from temperature gradients. The dimensionless figure of merit, ZT = S 2 ρ −1 κ −1 T , is calculated from the Seebeck coefficient ( S ), electrical resistivity ( ρ ), and thermal conductivity ( κ ). The calculated efficiency based on ZT using the conventional formula is not reliable in some cases due to the assumption of temperature-independent S , ρ , and κ . We established a new efficiency formula by introducing an engineering figure of merit ( ZT ) eng and an engineering power factor ( PF ) eng to predict reliably and accurately the efficiency of materials at a large temperature difference between the hot and cold sides, unlike the conventional ZT and PF providing performance only at specific temperatures. These new formulas will profoundly impact the search for new thermoelectric materials.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗
DOE OSTI2025

Range-separated thermal hybrid exchange-correlation density functional for accurate band-gap calculations of warm dense matter

Here, in this study, we present a range-separated thermal hybrid exchange-correlation (XC) functional, termed RS-KDT0, which provides highly accurate band-gap predictions at elevated temperature while reducing the computational cost relative to thermal global hybrid functionals. At zero temperature, RS-KDT0 recovers the widely used ground-state range-separated hybrid functional HSE06. To assess its performance across temperature regimes, we conducted static calculations of electronic pressure, band gaps, and band structures. RS-KDT0 exhibits the same level of accuracy as HSE06 in the low-temperature regime and demonstrates improved agreement with finite-temperature GW calculations at elevated temperatures when compared to other thermal XC functionals such as the generalized gradient approximation (GGA)-level KDT16 and the global hybrid KDT0.

density functional theory↗
DOE OSTI2014

Effects of reflector and crystal surface on the performance of a depth-encoding PET detector with dual-ended readout

Purpose: Depth encoding detectors are required to improve the spatial resolution and spatial resolution uniformity of small animal positron emission tomography (PET) scanners, as well as dedicated breast and brain scanners. Depth of interaction (DOI) can be measured by using dual-ended readout of lutetium oxyorthosilicate (LSO) scintillator arrays with position-sensitive avalanche photodiodes. Inter-crystal reflectors and crystal surface treatments play important roles in determining the performance of dual-ended detectors. In this paper, the authors evaluated five LSO arrays made with three different intercrystal reflectors and with either polished or unpolished crystal surfaces. Methods: The crystal size in all arrays was 1.5 mm, which is typical of the detector size used in small animal and dedicated breast scanners. The LSO arrays were measured with dual-ended readout and were compared in terms of flood histogram, energy resolution, and DOI resolution performance. Results: The four arrays using enhanced specular reflector (ESR) and Toray reflector provided similar quality flood histograms and the array using Crystal Wrap reflector gave the worst flood histogram. The two arrays using ESR reflector provided the best energy resolution and the array using Crystal Wrap reflector yielded the worst energy resolution. All arrays except the polished ESR array provided good DOI resolution ranging from 1.9 mm to 2.9 mm. DOI resolution improved as the gradient in light collection efficiency with depth (GLCED) increased. The geometric mean energies were also calculated for these dual-ended readout detectors as an alternative to the conventional summed total energy. It was shown that the geometric mean energy is advantageous in that it provides more uniform photopeak amplitude at different depths for arrays with high GLCED, and is beneficial in event selection by allowing a fixed energy window independent of depth. A new method of DOI calculation that improved the linearity of DOI ratio vs depth and simplifies the DOI calibration procedure also was developed and tested. Conclusions: The results of these studies provide useful guidance in selecting the proper reflectors and crystal surface treatments when LSO arrays are used for high-resolution PET applications in small animal scanners or dedicated breast and brain scanners.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-F-BRD-05: Robustness of Dose Painting by Numbers in Proton Therapy

Purpose: Proton range uncertainties may cause important dose perturbations within the target volume, especially when steep dose gradients are present as in dose painting. The aim of this study is to assess the robustness against setup and range errors for high heterogeneous dose prescriptions (i.e., dose painting by numbers), delivered by proton pencil beam scanning. Methods: An automatic workflow, based on MATLAB functions, was implemented through scripting in RayStation (RaySearch Laboratories). It performs a gradient-based segmentation of the dose painting volume from 18FDG-PET images (GTVPET), and calculates the dose prescription as a linear function of the FDG-uptake value on each voxel. The workflow was applied to two patients with head and neck cancer. Robustness against setup and range errors of the conventional PTV margin strategy (prescription dilated by 2.5 mm) versus CTV-based (minimax) robust optimization (2.5 mm setup, 3% range error) was assessed by comparing the prescription with the planned dose for a set of error scenarios. Results: In order to ensure dose coverage above 95% of the prescribed dose in more than 95% of the GTVPET voxels while compensating for the uncertainties, the plans with a PTV generated a high overdose. For the nominal case, up to 35% of the GTVPET received doses 5% beyond prescription. For the worst of the evaluated error scenarios, the volume with 5% overdose increased to 50%. In contrast, for CTV-based plans this 5% overdose was present only in a small fraction of the GTVPET, which ranged from 7% in the nominal case to 15% in the worst of the evaluated scenarios. Conclusion: The use of a PTV leads to non-robust dose distributions with excessive overdose in the painted volume. In contrast, robust optimization yields robust dose distributions with limited overdose. RaySearch Laboratories is sincerely acknowledged for providing us with RayStation treatment planning system and for the support provided.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1972

Diffuse galactic X-rays from discrete sources.

The galactic ridge is interpreted as a population of intrinsically weak X-ray sources, and calculations are presented of the predicted intensity as a function of scale height and radial gradient. Observations of diffuse galactic X-radiation above 1 keV impose severe constraints on discrete-source interpretations of the galactic-plane excess observed at 0.25 keV.

Silk, J.↗
NASA NTRS1974

Control system design using frequency domain models and parameter optimization, with application to supersonic inlet controls

A technique is described for designing feedback control systems using frequency domain models, a quadratic cost function, and a parameter optimization computer program. FORTRAN listings for the computer program are included. The approach is applied to the design of shock position controllers for a supersonic inlet. Deterministic or random system disturbances, and the presence of random measurement noise are considered. The cost function minimization is formulated in the time domain, but the problem solution is obtained using a frequency domain system description. A scaled and constrained conjugate gradient algorithm is used for the minimization. The approach to a supersonic inlet included the calculations of the optimal proportional-plus integral (PI) and proportional-plus-integral-plus-derivative controllers. A single-loop PI controller was the most desirable of the designs considered.

Seidel, R. C.↗
NASA NTRS1986

Analyses of convective flows and spatial gradients in energetic ion observations

Based on the measured anisotropy of three-dimensional energetic ion distributions obtained from the Medium Energy Particle Instrument experiment of the ISEE 1 spacecraft, bulk flow, epsilon x B convection, and pressure gradients in magnetospheric plasmas are measured with two techniques. In the first technique, using a Compton-Getting transform into the appropriate rest frame in an iterative fashion and converging to the unique velocity that removes the first-order anisotropy from the flux angular distribution, the actual measured angular distribution, and spectrum at each angle, are employed. The second technique uses a first moment of the distribution function in which the shift in velocity space of the distribution function is obtained at two phase space density values. The convection contribution to the velocity shift is assumed to be perpendicular to the pressure gradient contribution if an MHD generator condition is assumed, rendering the two shifts separable, and a gradient scale and a lower limit to the convection electric field are derived in one calculation.

Mitchell, D. G.↗