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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-T-109: An Investigation of Including Variable Relative Biological Effectiveness in Intensity Modulated Proton Therapy Planning Optimization for Head and Neck Cancer Patients

Purpose: The current practice of considering the relative biological effectiveness (RBE) of protons in intensity modulated proton therapy (IMPT) planning is to use a generic RBE value of 1.1. However, RBE is indeed a variable depending on the dose per fraction, the linear energy transfer, tissue parameters, etc. In this study, we investigate the impact of using variable RBE based optimization (vRBE-OPT) on IMPT dose distributions compared by conventional fixed RBE based optimization (fRBE-OPT). Methods: Proton plans of three head and neck cancer patients were included for our study. In order to calculate variable RBE, tissue specific parameters were obtained from the literature and dose averaged LET values were calculated by Monte Carlo simulations. Biological effects were calculated using the linear quadratic model and they were utilized in the variable RBE based optimization. We used a Polak-Ribiere conjugate gradient algorithm to solve the model. In fixed RBE based optimization, we used conventional physical dose optimization to optimize doses weighted by 1.1. IMPT plans for each patient were optimized by both methods (vRBE-OPT and fRBE-OPT). Both variable and fixed RBE weighted dose distributions were calculated for both methods and compared by dosimetric measures. Results: The variable RBE weighted dose distributions were more homogenous within the targets, compared with the fixed RBE weighted dose distributions for the plans created by vRBE-OPT. We observed that there were noticeable deviations between variable and fixed RBE weighted dose distributions if the plan were optimized by fRBE-OPT. For organs at risk sparing, dose distributions from both methods were comparable. Conclusion: Biological dose based optimization rather than conventional physical dose based optimization in IMPT planning may bring benefit in improved tumor control when evaluating biologically equivalent dose, without sacrificing OAR sparing, for head and neck cancer patients. The research is supported in part by National Institutes of Health Grant No. 2U19CA021239-35.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1979

Mathematical model for steady state, simple ampholyte isoelectric focusing: Development, computer simulation and implementation

The elimination of Ampholine from the system by establishing the pH gradient with simple ampholytes is proposed. A mathematical model was exercised at the level of the two-component system by using values for mobilities, diffusion coefficients, and dissociation constants representative of glutamic acid and histidine. The constants assumed in the calculations are reported. The predictions of the model and computer simulation of isoelectric focusing experiments are in direct importance to obtain Ampholine-free, stable pH gradients.

Palusinski, O. A.↗
DOE OSTI2019

Magnetic and Thermal Design of HTS Quadrupole Magnet for Newly Developed Superconducting Proton Cyclotron Beam Line

A magnetic and thermal design of high-temperature superconducting (HTS) quadrupole magnet for newly developed superconducting proton cyclotron was presented in this research. With superconducting technology, the design can reduce the magnet size and remove the heat loads more efficiently. Calculations are conducted with finite element method (FEM) to study the quadrupole magnet’s magnetic and thermal properties. For the magnet cold-mass and cryostat system, the heat load is mainly generated by conduction and radiation. Multilayer thermal insulation and G10 supports are used to restrict them. The magnetic field distribution in the HTS coils is calculated to consider the critical current degradation of the YBCO tapes. The effects on the field gradient quality in relation to the pole tip profile and end chamfer are analyzed. The sixth multipole component in the integral field can be improved with end chamfer. Finally, a 20 T/m HTS quadrupole magnet with integral field uniformity less than 0.05% out to 75% of the inscribed radius is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Crystallographic texture gradient along the wall thickness of an SF-copper tube

Metallic tubes usually show heterogeneities along circumference and through wall thickness. Local texture analysis by synchrotron diffraction and global texture analysis by neutron diffraction was used describing texture variations related to the wall thickness. Test sample was an SF-Cu tube with outer diameter of 140 mm and an average wall thickness of 10 mm. The texture is dominated by the cube component {001} <100>. Other minor components are rotated cube {013} <100>, Goss {011} <100>, Cu {112} <111>, brass {110} <112>, S {123} <634>. The texture gradient over the wall thickness shows a remarkable variation of the sharpness of the cube component. Other components show only little variation. Comparing samples obtained at maximum wall thickness and at minimum wall thickness one can see that the wall thickness has an influence on the global texture as well as texture gradient. Another interesting result is, that global textures measured directly by neutrons and averaged global textures, calculated from the set of local measurements by synchrotron diffraction, agree surprisingly well despite their different resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Toward accurate measurement of electromagnetic field by retrieving and refining the center position of non-uniform diffraction disks in Lorentz 4D-STEM

Recent advancement in scanning transmission electron microscopy (STEM) allows the use of 4D-STEM, a technique that captures an electron diffraction pattern at each scan point in STEM, to measure electrostatic and magnetic potential and field in materials. However, accurate measurement, separation of the magnetic and electric signals, and removal of artifacts remain challenging, especially in the presence of complex non-uniform diffraction contrast within the disks. In this work, based on dynamic simulations of 4D-STEM patterns built upon superstructures consisting of millions of atoms to account for different sample thickness and edge geometries, we show how the shape and intensity distribution of the central disk are affected by multiple scattering. We propose a robust refinement procedure through iteration of the spin-sensitive peak position of the disk-center in the circular Hough transform filtered images from experimental Lorentz 4D-STEM dataset after minimizing the possible artifacts, such as those due to the change of thickness, dynamic scattering, and scanning process. We verify that caution must be taken as in practice the rigid-disk-shift model used to reconstruct induction maps can easily break down due to disk-protrusion when there exists a nonconstant phase gradient or thickness within the width of the probe. Through quantitative analysis and comparing experiment with calculation the effect of the non-spin-related intensity distribution inside the disk as well as that causes the disk shift due to the intensity-protrusion can be removed, and high-quality magnetic field mapping is possible.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Precursor Ion–Ion Aggregation in the Brust–Schiffrin Synthesis of Alkanethiol Nanoparticles

Tetraoctylammonium bromide is used in the Brust-Schiffrin nanoparticle synthesis to phase-transfer chloroaurate ions from the aqueous phase to the organic phase. While it is established that the quaternary ammonium complex self-associates in the organic phase, the actual self-assembled structure is debated. We have confirmed the presence of ion-ion aggregates through quantitative 1H Nuclear Magnetic Resonance spectroscopy (NMR), pulsed field gradient, diffusion-ordered NMR (DOSY-NMR) and density functional theory (DFT) based NMR shift calculations. Tetraoctylammonium complexes (TOA-X, where X = Br, Cl, AuCl4-xBrx, AuBr4/Br and AuCl4-xBrx/Br) were investigated to measure the extraction of water into the organic phase. 1H NMR and DFT based NMR shielding calculations indicated that deshielding of water is due to hydration of the anion and not the formation of the aqueous core of a reverse micelle. DOSYNMR results were consistent with the formation of small aggregates at typical Brust-Schiffrin synthesis concentrations. The extent of aggregation correlated with the size and electronegativity of the anion and was analyzed with a modified, isodesmic, indefinite aggregation model. The substitution of bromoauric acid for chlororoauric acid at conditions emulating the Brust-Schiffrin synthesis increased the aggregation of the quaternary ammonium complex. The increase in aggregation corresponded with an increase in the size of the produced nanoparticles from 4.3 to 4.6 nm. Understanding the selfassembly and supramolecular structure of precursors in the Brust-Schiffrin synthesis will enable further refinement of models that predict the growth of noble metal nanoparticles.

Density Functional Theory↗
DOE OSTI2018

Reevaluating growing season length controls on net ecosystem production in evergreen conifer forests

Abstract Growing season length (GSL) is a key unifying concept in ecology that can be estimated from eddy covariance-derived estimates of net ecosystem production (NEP). Previous studies disagree on how increasing GSLs may affect NEP in evergreen coniferous forests, potentially due to the variety of methods used to quantify GSL from NEP. We calculated GSL and GSL-NEP regressions at eleven evergreen conifer sites across a broad climatic gradient in western North America using three common approaches: (1) variable length (3–7 days) regressions of day of year versus NEP, (2) a smoothed threshold approach, and (3) the carbon uptake period, followed by a new approach of a method-averaged ensemble. The GSL and the GSL-NEP relationship differed among methods, resulting in linear relationships with variable sign, slope, and statistical significance. For all combinations of sites and methods, the GSL explained between 6% and 82% of NEP with p -values ranging from 0.45 to < 0.01. These results demonstrate the variability among GSL methods and the importance of selecting an appropriate method to accurately project the ecosystem carbon cycling response to longer growing seasons in the future. To encourage this approach in future studies, we outline a series of best practices for GSL method selection depending on research goals and the annual NEP dynamics of the study site(s). These results contribute to understanding growing season dynamics at ecosystem and continental scales and underscore the potential for methodological variability to influence forecasts of the evergreen conifer forest response to climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Electronic structural and magnetic properties of Mn{sub 5}Ge{sub 3} clusters

Theoretical understanding of the stability, ferromagnetism, and spin polarization of Mn{sub 5}Ge{sub 3} clusters has been performed by using the density functional theory with generalized gradient approximation for exchange and correlation. The magnetic moments and magnetic anisotropy energy (MAE) have been calculated for both bulk and clusters, and the enhanced magnetic moment as well as the enlarged MAE have been identified in clusters. The most attractive achievement is that Mn{sub 5}Ge{sub 3} clusters show a fine half-metallic character with large energy scales. The present results may have important implications for potential applications of small Mn{sub 5}Ge{sub 3} clusters as both emerging spintronics and next-generation data-storage technologies.

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

Adsorption of silver dimer on graphene - A DFT study

We performed a systematic density functional theory (DFT) study of the adsorption of silver dimer (Ag{sub 2}) on graphene using SIESTA (Spanish Initiative for Electronic Simulations with Thousands of Atoms) package, in the generalized gradient approximation (GGA). The adsorption energy, geometry, and charge transfer of Ag2-graphene system are calculated. The minimum energy configuration for a silver dimer is parallel to the graphene sheet with its two atoms directly above the centre of carbon-carbon bond. The negligible charge transfer between the dimer and the surface is also indicative of a weak bond. The methodology demonstrated in this paper may be applied to larger silver clusters on graphene sheet.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Structural, bonding, and electronic properties of the hexagonal ferroelectric and paraelectric phases of LuMnO{sub 3} compound: A density functional theory study

We have investigated the structural, bonding, and electronic properties of both ferroelectric (FE) and paraelectric (PE) phases of the hexagonal LuMnO{sub 3} compound using calculations based on density functional theory. The structural properties have been determined by employing the generalized gradient approximation with Perdew-Burke-Ernzerhof and Wu-Cohen parameterization. The bonding and electronic properties have been treated by recently developed modified Becke-Johnson exchange potential, which succeeded to open a band gap for both PE and FE phases, in agreement with experimental predictions. The Bader’s topological analysis of electronic density showed that the character of the Lu–O axial bonds changes when the crystal exhibits the PE → FE structural transition. This fact is in agreement with experimental findings. The covalent character of the Lu–O bond significantly increases due to orbital hybridization between the Lu 5d{sub z}{sup 2} and O 2p{sub z}-states. This bonding mechanism causes the ferroelectricity in the hexagonal LuMnO{sub 3} compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD{sub 3}CN treated with a parallel multi-state EVB model

We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD{sub 3}CN → DF + CD{sub 2}CN reaction in CD{sub 3}CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD{sub 3}CN solvent, equilibrium power spectra of DF in CD{sub 3}CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol{sup −1} localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD{sub 3}CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.

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

Methods of improving spatial resolution for IR spectroscopy in atmospheric-pressure plasma systems

High-temperature, atmospheric pressure plasma systems operated in molecular gases present complex diagnostic challenges. Infrared spectroscopy has been used to make measurements of the absorbance spectrum of complex molecular gas mixtures and thereby calculate the concentrations and species temperatures in these systems. For atmospheric-pressure systems, high spatial gradients arise and high spatial-resolution measurements are thus desirable. Some systems have achieved increased spatial resolution by reducing the beam diameter. However, this increase in spatial resolution comes at the expense of the optical throughput. Here we propose modifying a commercial Fourier Transform, Infrared spectrometer system with a set of simple optical elements. This design allows for scannable, high spatial resolution absorbance spectrum measurements of a complex molecular gas mix. We analyze the signal-to-noise ratio for this design and compare it with the diameter reduction method. For sufficiently low plasma radiation levels, this design improves the signal-to-noise ratio relative to reducing the beam diameter size while transmitting nearly all of the beam power. Furthermore, whether this design can achieve the desired sub-mm spatial resolution depends on the radiating area of the infrared light source in the spectrometer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Enhanced pedestal H-mode at low edge ion collisionality on NSTX

The Enhanced Pedestal (EP) H-mode regime is an attractive wide-pedestal high-βp scenario for the National Spherical Torus Experiment Upgrade (NSTX-U) and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (βN > 5), and significant bootstrap fraction (fBS > 0.6) at Ip/BT = 2 MA/T. This regime is realized when the edge ion collisionality becomes sufficiently small that a positive feedback interaction occurs between a reduction in the ion neoclassical energy transport and an increase in the particle transport from pressure-driven edge instabilities. The EP H-mode was most often observed as a transition following a large edge-localized mode in conditions with low edge neutral recycling. It is hypothesized that the onset of pressure-driven instabilities prior to the full recovery of the neutral density leads to a temporary period with an elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy and magnetic spectroscopy in order to describe the evolution of the edge particle transport mechanisms during the ELM recovery and the saturated EP H-mode state. The observations are consistent with the hypothesis that the onset of pressure-driven edge instabilities, such as the kinetic ballooning mode and kink-peeling mode, can be responsible for the increased particle transport in EP H-mode.

Battaglia, D. J. (ORCID:0000000188979740)↗
DOE OSTI2025

Molybdenum density profiles in Alcator C-Mod plasmas

Molybdenum (Z = 42) density spatial profiles have been obtained for Alcator C-Mod L-, I-, and H-mode plasmas from x-ray observations of the neon-like Mo 32+ 4D transition. Profile shapes range from flat to modestly peaked in the core. The observed peaking factors for molybdenum are in good agreement with those calculated by TGLF across all confinement regimes, and there is little dependence seen with core density gradient, sawtooth amplitude, or rotation velocity. Predictions for tungsten peaking in SPARC are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Relativistic Bondi–Hoyle–LYttleton accretion onto a rotating black hole: density gradients

In this work, for the first time, we present a numerical study of Bondi–Hoyle accretion with density gradients in the fully relativistic regime. In this context, we consider the accretion onto a Kerr black hole (BH) of a supersonic ideal gas with density gradients perpendicular to the relative motion. The parameters of interest in this study are the Mach number, M, the spin of the BH, a, and the density-gradient parameter of the gas, ϵ{sub ρ}. We show that, unlike in the Newtonian case, all of the studied cases, especially those with a density gradient, approach a stationary flow pattern. To illustrate that the system reaches a steady state, we calculate the mass and angular momentum accretion rates on a spherical surface almost located at the event horizon. In the particular case of M=1, ϵ{sub ρ} = 0.5, and BH spin a = 0.5, we observe a disk-like configuration surrounding the BH. Finally, we present the gas morphology and some of its properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Viable cosmology in bimetric theory

We study cosmological perturbations in bimetric theory with two fluids each of which is coupled to one of the two metrics. Focusing on a healthy branch of background solutions, we clarify the stability of the cosmological perturbations. For this purpose, we extend the condition for the absence of the so-called Higuchi ghost, and show that the condition is guaranteed to be satisfied on the healthy branch. We also calculate the squared propagation speeds of perturbations and derive the conditions for the absence of the gradient instability. To avoid the gradient instability, we find that the model parameters are weakly constrained.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

PFPO plasma scenarios for exploration of long pulse operation in ITER

Long Pulse Scenarios (LPS) in ITER foreseen during the Pre-Fusion Power Operation (PFPO) phase of the ITER Research Plan (IRP) are assessed using 1.5D transport simulations within the ASTRA framework. Such assessment is required to predict the operational space for LPS operation in PFPO, as well as to evaluate which physics processes for LPS operation during Fusion Power Operation (FPO) could be studied during PFPO. An important aspect in the development of LPSs in PFPO is to minimize lifetime consumption of the Central Solenoid (CS) for these scenarios. The maximum pulse length achievable for LPSs in PFPO with no consumption of CS lifetime (currents in CS coils $&LessSlantEqual;$30 kA per turn) has been assessed for a range of heating schemes and heating mixes, confinement regimes (L-mode and H-mode) and for helium and hydrogen plasmas. The operational space of LPS and pulse length has been explored through density scans with the Heating and Current Drive mix required for the FPO Q $&GreaterSlantEqual;$ 5 steady-state plasma scenario (namely Neutral Beam Injection and Electron Cyclotron Heating) including acceptable shine through losses on the first wall for both helium and hydrogen plasmas. Fast particle physics aspects that are common between FPO plasmas and LPS PFPO H-mode plasmas at low densities are studied including MHD stability analysis with the KINX code and non-perturbative critical gradient model based on high-n Toroidal Alfven Eigenmodes (TAE) stability kinetic ballooning code HINST calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies

We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke- Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.

36 MATERIALS SCIENCE↗