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 991 records · Page 55

Ab-initio Cu alloy design for high-gradient accelerating structures

Operation of normal conducting accelerator structures at high accelerating gradients is beneficial for many accelerator applications in basic science, industry, medicine, and National Security. RF breakdown is the major factor that limits the achievable accelerating gradients. Previous experiments on copper (Cu) have demonstrated that RF breakdown probability can be significantly decreased by hardening the material and alloying Cu with solutes such as silver (Ag). In this paper, we propose a figure-of-merit (FOM) that characterizes the ability of Cu alloys to withstand high-gradients. The FOM represents a trade-off between hardening through solid solution strengthening and the additional thermal stress induced by incremental RF pulse heating resulting from changes in electronic properties induced by alloying. We performed high-throughput ab initio calculations and computed the FOM for a large number of binary Cu alloys. Several promising candidate alloys for high-gradient accelerating structures were identified, such as CuAg, CuCd, CuHg, CuAu, CuIn, and CuMg. CuAg alloys have previously exhibited low RF breakdown rates in experiments. The results provide guidance for selecting alloys for the future high-gradient normal conducting accelerating structures operating at very high gradients.

36 MATERIALS SCIENCE↗

WE-D-BRA-05: Pseudo In Vivo Patient Dosimetry Using a 3D-Printed Patient-Specific Phantom

Purpose: To test the feasibility of using 3D-printed patient-specific phantoms for intensity-modulated radiation therapy (IMRT) quality assurance (QA). Methods: We created a patient-specific whole-head phantom using a 3D printer. The printer data file was created from high-resolution DICOM computed tomography (CT) images of 3-year old child treated at our institution for medulloblastoma. A custom-modified extruder system was used to create tissue-equivalent materials. For the printing process, the Hounsfield Units from the CT images were converted to proportional volumetric densities. A 5-field IMRT plan was created from the patient CT and delivered to the 3D- phantom. Dose was measured by an ion chamber placed through the eye. The ion chamber was placed at the posterior edge of the planning target volume in a high dose gradient region. CT scans of the patient and 3D-phantom were fused by using commercial treatment planning software (TPS). The patient’s plan was calculated on the phantom CT images. The ion chamber’s active volume was delineated in the TPS; dose per field and total dose were obtained. Measured and calculated doses were compared. Results: The 3D-phantom dimensions and tissue densities were in good agreement with the patient. However, because of a printing error, there was a large discrepancy in the density in the frontal cortex. The calculated and measured treatment plan doses were 1.74 Gy and 1.72 Gy, respectively. For individual fields, the absolute dose difference between measured and calculated values was on average 3.50%. Conclusion: This study demonstrated the feasibility of using 3D-printed patient-specific phantoms for IMRT QA. Such phantoms would be particularly advantageous for complex IMRT treatment plans featuring high dose gradients and/or for anatomical sites with high variation in tissue densities. Our preliminary findings are promising. We anticipate that, once the printing process is further refined, the agreement between measured and calculated doses will improve.

60 APPLIED LIFE SCIENCES↗

Study of Quadrupole Fringe Fields in the Interaction Region of the Hadron Storage Ring of the Electron Ion Collider

Fringe fields in quadrupole magnets are usually neglected in studies of beam dynamics at accelerators. However, the extreme optical parameters present in the final focus of a collider such as the Electron–Ion Collider (EIC) may give rise to effects that should not be overlooked. The calculation of quadrupole fringe fields presented in this study follows the procedure outlined in Ref. [1], specialized to the case of a straight reference orbit (i.e., with no dipole field component). A right-handed Cartesian coordinate system is employed, with the $z$-axis aligned with the quadrupole axis and $x$ and $y$ denoting the horizontal and vertical transverse coordinates, respectively. The magnetic quadrupole field gradient, $G = \partial B_y / \partial x$, transitions from its peak value inside the quadrupole—where it is nearly independent of the longitudinal coordinate $z$—to zero at some distance beyond the magnet edge. Consequently, $G$ is treated as a function of $z$. The region over which this variation occurs is defined as the quadrupole fringe field region. The study begins with the development of a description of the magnetic field in the fringe region using a power-series expansion in the transverse coordinates $x$and $y$, consistent with the longitudinally varying gradient. A model for the $z$-dependence of the gradient is then proposed and adjusted to reproduce magnetic field data obtained from three-dimensional field calculations. To evaluate the impact of the fringe fields on beam dynamics, the corresponding vector potential is derived and incorporated into a Hamiltonian formulation of particle motion. The significance of the fringe fields is quantified by calculating the amplitude-dependent tune shift from the Hamiltonian. Using linear beam optics parameters of the Hadron Storage Ring (HSR) of the EIC, the tune shift due to the fringe fields of all quadrupole magnets in the IR-6 interaction region is evaluated. Finally, the resulting tune shifts are compared with those arising from other nonlinear field components present in the HSR.

43 PARTICLE ACCELERATORS↗

Infinite neural network quantum states: entanglement and training dynamics

We study infinite limits of neural network quantum states (∞-NNQS), which exhibit representation power through ensemble statistics, and also tractable gradient descent dynamics. Ensemble averages of entanglement entropies are expressed in terms of neural network correlators, and architectures that exhibit volume-law entanglement are presented. The analytic calculations of entanglement entropy bound are tractable because the ensemble statistics are simplified in the Gaussian process limit. A general framework is developed for studying the gradient descent dynamics of neural network quantum states (NNQS), using a quantum state neural tangent kernel (QS-NTK). For ∞-NNQS the training dynamics is simplified, since the QS-NTK becomes deterministic and constant. An analytic solution is derived for quantum state supervised learning, which allows an ∞-NNQS to recover any target wavefunction. Numerical experiments on finite and infinite NNQS in the transverse field Ising model and Fermi Hubbard model demonstrate excellent agreement with theory. ∞-NNQS opens up new opportunities for studying entanglement and training dynamics in other physics applications, such as in finding ground states.

quantum state supervised learning↗

The AGS Ggamma Meter and Calibrating the Gauss Clock

During AGS Polarized Proton acceleration periods, one output from the AGS Ggamma Meter, namely the energy (or Ggamma) calculated from the magnetic field in the AGS main magnets and the beam radius- both measured in particular instant, is used to figure out the times in the AGS magnet acceleration cycle when the beam passes through a particular set of depolarizing resonances. The resonance set occur whenever a particle’s Ggamma (energy*(G/m) becomes nearly equal to n*Qx (i.e. any integer multiplied by the horizontal betatron tune). This deliverable is why the machinery is referred to as the ''Ggamma Meter'' rather than the AGS energy meter. The Ggamma Meter takes as inputs a set of measurements of frequency (F(t)), radius (r(t)), and gauss clock counts (GCC(t)). The other energy (GgammaBr) assumes the field when the gauss clock starts counting is known. The change in field to time t is given by the measured accumulated gauss clock counts multiplied by the gauss clock calibration (gauss/GCC). In order to deal with experimental data, this calibration factor gets an added ad hoc complication, namely a correction dependent on the rate of change the counting rate. The Ggamma meter takes GCC(t) and together with the past history for this cycle calculates B(t).

43 PARTICLE ACCELERATORS↗

Phase Transitions and Melting in the Venusian Basaltic Crust: Implications for Crustal Recycling

Understanding phase transitions and melting in the Venusian crust and the associated changes in density are critical to constrain crustal thickness, recycling, and remelting processes. For example, average surface conditions on Venus of 92 bars and 460 °C correspond to low-grade metamorphic conditions on Earth equivalent to hornfels- and greenschist-facies. Further, the primary igneous mineralogy in the Venusian crust should be replaced by metamorphic minerals over time and potentially melt, as temperatures and pressures increase with depth along geothermal gradients. Here, we use Perple_X (Connolly, 2005), a Gibbs free energy minimization program to calculates table phase equilibria over a range of pressure and temperature conditions for whole-rock compositions of a dry basalt, alkali basalt, and peridotite. Further, we extract mineral abundancies along five geotherms of 5, 10, 15, 20 and 25 °C/km and calculate the rock density after the extraction of melt. On the coldest geotherm of 5 °C/km, melting would start at ~ 197 km and ~121 km in the basalt and alkali basalt, respectively, and therefore at depths that significantly exceed the estimated thickness of the Venusian crust (8-25 km). A density cross-over, where the crustal density would exceed that of the mantle could potentially induce delamination and is estimated to occur at ~ 40 km for both compositions. On the 10 °C/km geotherm, melting of the basalt starts at 65 km and melt extraction only causes a gradual densification due the formation of initially lower melt proportions. In contrast, the alkali basalt starts to melt at a depth of ~ 57 km and the subsequent densification of the residual composition could trigger delamination at ~ 60km depth. For hotter geotherms (i.e., 25 °C/km gradient), melting would occur at shallower depths of 21 – 23 km. However, the basalt will not get significantly denser than the mantle while the residuum of the alkali basalt would reach a density cross-over at 30 km. Hence, differences in Venusian crustal compositions can significantly influence the thickness of the crust on individual geotherms. Our results demonstrate that phase transitions and melting are strongly dependent on the assumed geotherm and could, combined with geodynamic models, further constrain crustal parameters.

Julia Semprich↗

Effects of Spatial Gradients on Electron Runaway Acceleration

The runaway process is known to accelerate electrons in many laboratory plasmas and has been suggested as an acceleration mechanism in some astrophysical plasmas, including solar flares. Current calculations of the electron velocity distributions resulting from the runaway process are greatly restricted because they impose spatial homogeneity on the distribution. We have computed runaway distributions which include consistent development of spatial gradients in the energetic tail. Our solution for the electron velocity distribution is presented as a function of distance along a finite length acceleration region, and is compared with the equivalent distribution for the infinitely long homogenous system (i.e., no spatial gradients), as considered in the existing literature. All these results are for the weak field regime. We also discuss the severe restrictiveness of this weak field assumption.

MacNeice, Peter↗

Optimization of Depletion Modeling and Simulation for the High Flux Isotope Reactor

Monte Carlo based depletion tools used for the high-fidelity modeling and simulation of the High Flux Isotope Reactor (HFIR) come at a great computational cost; finding sufficient approximations is necessary to make the use of these tools feasible. The optimization of the neutronics and depletion model for the HFIR is based on two factors: (i) the explicit representation of the involute fuel plates with sets of polyhedra and (ii) the treatment of depletion mixtures and control element position during depletion calculations. A very fine representation (i.e., more polyhedra in the involute plate approximation) does not significantly improve simulation accuracy. The recommended representation closely represents the physical plates and ensures sufficient fidelity in regions with high flux gradients. Including the fissile targets in the central flux trap of the reactor as depletion mixtures has the greatest effect on the calculated cycle length, while localized effects (e.g., the burnup of specific isotopes or the power distribution evolution over the cycle) are more noticeable consequences of including a critical control element search or depleting burnable absorbers outside the fuel region.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Atlantic Tropical Cyclone Monitoring with AMSU-A: Estimation of Maximum Sustained Wind Speeds

The first Advanced Microwave Sounding Unit temperature sounder (AMSU-A) was launched on the NOAA-15 satellite on 13 May 1998. The AMSU-A's higher spatial and radiometric resolutions provide more useful information on the strength of the middle- and upper-tropospheric warm cores associated with tropical cyclones than have previous microwave temperature sounders. The gradient wind relationship suggests that the temperature gradient near the core of tropical cyclones increases nonlinearly with wind speed. The gradient wind equation is recast to include AMSU-A-derived variables, Stepwise regression is used to determine which of these variables is most closely related to maximum sustained winds (V(sub max)). The satellite variables investigated include the radially averaged gradients at two spatial resolutions of AMSU-A channels 1-10 T(sub b) data (delta(sub r)T(sub B)), the squares of these gradients, a channel-15-based scattering index (SI(sub 89)), and area-averaged T(sub B). Calculations of T(sub B) and delta(sub r)T(sub B) from mesoscale model simulations of Andrew reveal the effects of the AMSU spatial sampling on the cyclone warm core presentation. Stepwise regression of 66 AMSU-A terms against National Hurricane Center V(sub max) estimates from the 1998 and 1999 Atlantic hurricane season confirms the existence of a nonlinear relationship between wind speed and radially averaged temperature gradients near the cyclone warm core. Of six regression terms, four are dominated by temperature information, and two are interpreted as correcting for hydrometeor contamination. Jackknifed regressions were performed to estimate the algorithm performance on independent data. For the 82 cases that had in situ measurements of V(sub max), the average error standard deviation was 4.7 m/s. For 108 cases without in situ wind data, the average error standard deviation was 7.5 m/s Operational considerations, including the detection of weak cyclones and false alarm reduction, are also discussed.

Spencer, Roy W.↗

Efficient Discovery of Optimal N -Layered TMDC Hetero-Structures

ABSTRACT Vertical hetero-structures made from stacked monolayers of transition metal dichalcogenides (TMDC) are promising candidates for next-generation optoelectronic and thermoelectric devices. Identification of optimal layered materials for these applications requires the calculation of several physical properties, including electronic band structure and thermal transport coefficients. However, exhaustive screening of the material structure space using ab initio calculations is currently outside the bounds of existing computational resources. Furthermore, the functional form of how the physical properties relate to the structure is unknown, making gradient-based optimization unsuitable. Here, we present a model based on the Bayesian optimization technique to optimize layered TMDC hetero-structures, performing a minimal number of structure calculations. We use the electronic band gap and thermoelectric figure of merit as representative physical properties for optimization. The electronic band structure calculations were performed within the Materials Project framework, while thermoelectric properties were computed with BoltzTraP. With high probability, the Bayesian optimization process is able to discover the optimal hetero-structure after evaluation of only ~20% of all possible 3-layered structures. In addition, we have used a Gaussian regression model to predict not only the band gap but also the valence band maximum and conduction band minimum energies as a function of the momentum.

Bassman Oftelie, Lindsay↗

Momentum transport in the vicinity of q{sub min} in reverse shear tokamaks due to ion temperature gradient turbulence

We present an analytic study of momentum transport of tokamak plasmas in the vicinity of minimum safety factor (q) position in reversed magnetic shear configuration. Slab ion temperature gradient modes with an equilibrium flow profile are considered in this study. Quasi-linear calculations of momentum flux clearly show the novel effects of q-curvature on the generation of intrinsic rotation and mean poloidal flow without invoking reflectional symmetry breaking of parallel wavenumber (k{sub ∥}). This q-curvature effect originates from the inherent asymmetry in k{sub ∥} populations with respect to a rational surface due to the quadratic proportionality of k{sub ∥} when q-curvature is taken into account. Discussions are made of possible implications of q-curvature induced plasma flows on internal transport barrier formation in reversed shear tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonempirical Semilocal Free-Energy Density Functional for Matter under Extreme Conditions

The potential for density functional calculations to predict the properties of matter under extreme conditions depends crucially upon having a non-empirical approximate free energy functional valid over a wide range of state conditions. Unlike the ground-state case, no such free-energy exchange- correlation (XC) functional exists. We remedy that with systematic construction of a generalized gradient approximation XC free-energy functional based on rigorous constraints, including the free energy gradient expansion. The new functional provides the correct temperature dependence in the slowly varying regime and the correct zero-T, high-T, and homogeneous electron gas limits. Application in Kohn-Sham calculations for hot electrons in a static fcc Aluminum lattice demon- strates the combined magnitude of thermal and gradient effects handled by this functional. Its accuracy in the increasingly important warm dense matter regime is attested by excellent agreement of the calculated deuterium equation of state with reference path integral Monte Carlo results at intermediate and elevated temperatures and by low density Al calculations over a wide T range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An integral analysis of transonic normal shock wave/turbulent boundary layer interactions in internal flow

An approximate integral viscous-inviscid interaction method is presented for calculating the development of a turbulent boundary layer subjected to a normal shock wave induced adverse pressure gradient in an internal axisymmetric flow. The inflow conditions and the downstream pressure are provided for the computation. In the supersonic region of shock pressure rise, the Prandtl-Meyer function is used to couple the viscous and inviscid flows. An analytical model for the coupling process is postulated and appropriate equations are defined. Downstream of the sonic point, one-dimensional inviscid flow is assumed for coupling with the viscous flow. The turbulent boundary layer is calculated using Green's integral lag-entrainment method. Comparisons of the solutions with the experimental data are made for interactions which are unseparated, near separation and separated. For comparison purposes, solutions to the time-dependent, mass-averaged, Navier-Stokes equations incorporating a two-equation, Wilcox-Rubesin turbulence model are also shown. The computed results from the integral method show good agreement with experimental data for unseparated interactions and reasonable agreement with the trend of the viscous effects when the interaction becomes increasingly separated.

Om, D.↗

Theoretical Study of Mechanical Stability and Physical Properties of Co{sub 2}V{sub 1−x}Zr{sub x}Ga

The mechanical stability and physical properties of Co{sub 2}V{sub 1−x}Zr{sub x}Ga (x = 0, 0.25, 0.50, 0.75, 1) have been predicted by using ab initio calculations based on density functional theory. The exchange-correlation potentials are treated within the generalized gradient approximation-Perdew-Burke and Ernzerhhof (GGA-PBE) and Hubbard coefficient (GGA+U). The investigated minimum lattice parameter of Co{sub 2}VGa is 5.8030 Å, which is in excellent agreement with available experimental and theoretical data and for Co{sub 2}V{sub 1−x}Zr{sub x}Ga (x = 0.25, 0.50, 0.75, 1) are 5.9026, 5.9625, 6.0060 and 6.0021 Å, respectively, which are predicted for the very first time. The total magnetic moments decrease with increasing concentration of x in Co{sub 2}V{sub 1−x}Zr{sub x}Ga, in agreement with the Slater-Pauling rule. Band structure and density of states calculations show that the minority spin channel exhibits band gaps of 0.43, 0.46, 0.52, 0.54 and 1.01 eV for GGA+U scheme around the Fermi level confirming that all the studied composites are half-metallic in nature. In this study, we also studied the elastic constant and it is established that all the materials are mechanically stable and ductile in nature. The Co{sub 2}VGa is stiffer than other materials and all the studied composites have anisotropic behaviour. Moreover, by using a quasi-harmonic Debye model and calculated elastic constant, the Debye temperature and temperature-dependent constant volume heat capacity have been investigated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First-Principles Calculations to Investigate the Structural, Electronic, and Half-Metallic Properties of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} (x = 0, 0.25, 0.5, 0.75, and 1) Quaternary Heusler Alloys

The structural, electronic, and magnetic properties of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} (x = 0, 0.25, 0.5, 0.75, and 1) quaternary Heusler alloys have been investigated using first-principles density functional theory within generalized gradient approximation (GGA) and local spin-density approximation with Hubbard-U corrections (LSDA+U). The calculations have been performed using the full-potential linearized augmented Plane wave plus local orbitals method. The crystal structures of Ti{sub 2}RhSn{sub 1-x}Si{sub x}, Ti{sub 2}RhSn{sub 1-x}Ge{sub x}, and Ti{sub 2}RhGe{sub 1-x}Si{sub x} were found to be cubic in (F-43m) symmetry for x = 0, 1 and cubic in (P-43m) symmetry for x = 0.25, 0.75, while tetragonal in (P-4m2) symmetry for x = 0.5. The results show that all Ti{sub 2}RhSn1-xSix, Ti{sub 2}RhGe{sub 0.75}Si{sub 0.25}, and Ti{sub 2}RhGe alloys are half-metallic ferrimagnetic. These new materials are good candidates for potential applications in spintronic.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Density functional theory investigation of the LiIn 1- x Ga x Se 2 solid solution

Here, the electronic structure and optical properties of the LiIn 1-x Ga x Se 2 (x=0, 0.25, 0.5, 0.75, 1) solid solution were studied by density functional theory (DFT) with pure functionals. The exchange-correlation is treated within the local density approximation (LDA) and generalized-gradient approximation (GGA). The electronic structures for each respective compound are discussed in detail. Calculations reveal that gallium incorporation can be used to tune the optical-electrical properties of the solid solution and correlates with the lattice parameter. The band gap trend of the LiIn 1-x Ga x Se 2 system follows a nonlinear behavior between the LiInSe 2 and LiGaSe 2 ternary boundaries. The bowing parameter is estimated to be on the order of 0.1- 0.3 eV at the point. Low-temperature optical absorption revealed a 30% change in the temperature dependence of the band gap for the intermediate compound LiIn 0.6 Ga 0.4 Se 2 compared to ternary boundaries and suggests the heat capacity to be another control element through strain.

36 MATERIALS SCIENCE↗

Ab Initio Investigation of Structural Stability and Electronic and Magnetic Properties of the Half-Heusler Alloys: MTiSb (M = Fe, Co, and Ni)

The structural, electronic, and magnetic properties of the half-Heusler MTiSb (M = Fe, Co, Ni) compounds have been examined using the first-principles calculations with the full-potential linear augmented plane wave (FP-LAPW) method within the density functional theory. The electronic structures and magnetic properties of MTiSb (M = Fe, Co, Ni) are studied. It is found that the calculated lattice parameter is in good agreement with the theoretical values. Using the general gradient approximation (GGA), we observe that for the three compounds, the α phase is more stable than β and γ phases. CoTiSb is a semiconductor, NiTiSb is a metal, while the FeTiSb compound is a half-metallic, and the total magnetic moments are 0.00085, − 0.26074, and 0.99353 μ B, respectively, which is in agreement with the Slater–Pauling rule, Furthermore, the origin for the appearance of the half-metallic band gap in the FeTiSb compound was also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Half-Metallic and Half-Semiconductor Gaps in Cr-Based Chalcogenides: DFT + U Calculations

Using full-potential linear-augmented plane waves plus local orbitals (FP-LAPW) method, within the generalized gradient approximation (GGA) and GGA + U (Hubbard Coulomb onsite correction) as well as the modified Becke-Johnson approach (mBJ-GGA and mBJ-GGA + U), we have carried out a systematic investigation of the structural, electronic, and magnetic properties of CrS, CrSe, and CrTe compounds in two competing phases: the zinc-blende (ZB) and wurtzite (WZ) structures. We obtained half-metallic (HM) gaps of 1.35 (2.32), 1.65 (2.49), and 1.52 (1.80) for the CrS, CrSe, and CrTe compounds, respectively, in the ZB phase within the GGA + U (mBJ-GGA + U) approach. The HM gaps are 1.78, 2.11, and 1.53 for CrS, CrSe, and CrTe compounds, respectively, in the WZ phase using the GGA + U. Notably, the CrSe and CrTe using the mBJ-GGA + U approach are a half-semiconductor (HSC) in the WZ phase. The HSC gaps are 0.94 and 0.25 eV for CrSe and CrTe compounds, respectively, while CrS has a HM gap of 2.30 eV. We observed a HSC-to-HM transition for the CrSe and CrTe crystal structures on application of strain. The CrS, CrSe, and CrTe compounds have a total magnetic moment of 4.0 µ{sub B} per formula unit with the magnetization predominantly from the Cr ion. The main features observed from the density of states and strain evolution should motivate further experimental exploration on the possible application of the CrS, CrSe, and CrTe compounds for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗