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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Dose-volume and radiobiological dependence on the calculation grid size in prostate VMAT planning

This study evaluated the effects of dose-volume and radiobiological dependence on the calculation grid size in prostate volumetric-modulated arc therapy (VMAT) planning. Ten patients with prostate cancer were selected for this retrospective treatment planning study. Prostate VMAT plans were created for the patients using the 6 MV photon beam produced by a Varian TrueBEAM linac with the calculation grid size equal to 1, 2, 2.5, 3, 4, and 5 mm. Dose-volume histograms (DVHs) of targets and organs at risk were generated for different grid sizes. We calculated the radiobiological parameters of the tumor control probability (TCP) of clinical target volume (CTV) and planning target volume (PTV), and the normal tissue complication probability (NTCP) of organs at risk (rectal wall, rectum, bladder wall, bladder, left femur, and right femur). The homogeneity, conformity, and gradient indexes of CTV and PTV were calculated for different grid sizes. The TCP of PTV was found decreasing with a rate of 0.06%/mm when the calculation grid size increased from 1 to 5 mm. On the other hand, both NTCPs of rectal wall and rectum were found decreasing with rates of 0.03%/mm and 0.05%/mm, respectively, with an increase of grid size. The homogeneity index of PTV increased with a rate of 0.57/mm of the calculation grid size, whereas the conformity index of PTV decreased with a rate of 0.0075/mm. The gradient index of PTV was found increasing with a rate equal to 0.05/mm. In prostate VMAT planning, variations of dose-volume and radiobiological parameters with calculation grid size on PTV, rectal wall, and rectum were more significant than those of CTV and other organs at risk such as bladder wall, bladder, and femurs. Results in this study are important in the treatment planning quality assurance when the calculation grid size is varied to compromise a shorter dose computing time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Reliable density functional calculations for the electronic structure of thermoelectric material ZnSb

In this paper, we present the results of systematic test calculations for the electronic structure of thermoelectric material ZnSb using a first-principles full-potential all electron computational method. We used a linear combination of atomic orbitals (LACO) formalism, based on density functional theory (DFT). The exchange-correlation interaction potential of the many electron system was described by using a generalized gradient approximation (GGA). We compared the calculated indirect and direct band gaps as well as the effective masses of holes and electrons in ZnSb with experimental measurement results. The calculated indirect band gap of ZnSb is 0.56 eV, which agrees very well with the experimentally measured values of 0.50 eV ~ 0.61 eV. The calculated direct band gap at X point is 0.89 eV. The calculated effective masses of electrons and holes in ZnSb also agree with experimental data. The systematical test calculations as well as the comparisons of the calculated results with experimental measurements show that the obtained electronic structure of ZnSb would be reliable. We did not observe a major deficiency of the first-principles DFT calculation for the electronic structure of ZnSb, using full-potential all electron LACO method. The reported electronic structure of single crystal ZnSb from this work may provide a fundamental knowledge base for further research and applications for this important thermoelectric material.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Determining Equilibrium Position For Acoustical Levitation

Equilibrium position and orientation of acoustically-levitated weightless object determined by calibration technique on Earth. From calibration data, possible to calculate equilibrium position and orientation in presence of Earth gravitation. Sample not levitated acoustically during calibration. Technique relies on Boltzmann-Ehrenfest adiabatic-invariance principle. One converts resonant-frequency-shift data into data on normalized acoustical potential energy. Minimum of energy occurs at equilibrium point. From gradients of acoustical potential energy, one calculates acoustical restoring force or torque on objects as function of deviation from equilibrium position or orientation.

Barmatz, M. B.↗
DOE OSTI2023

Benchmarks and results of the two-band Hubbard model from the Gutzwiller conjugate gradient minimization theory

Ground-state properties, such as energies and double occupancies, of a one-dimensional two-band Hubbard model are calculated using a first-principles Gutzwiller conjugate gradient minimization theory. The favorable agreement with the results from the density matrix renormalization group theory demonstrates the accuracy of our method. A rotationally invariant approach is further incorporated into the method to greatly reduce the computational complexity with a speedup of approximately 50 times. Moreover, we investigate the Mott transition between a metal and a Mott insulator by evaluating the charge gap. In conclusion, with greatly reduced computational effort, our method reproduces the phase diagram in reasonable agreement with the density matrix renormalization group theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Screening tool to evaluate the vulnerability of down-gradient receptors to groundwater contaminants from uncapped landfills

Highlights: • A spreadsheet-based risk screening tool for groundwater affected by landfills is presented. • Domenico solute transport equations are used to estimate downgradient contaminant concentrations. • Landfills are categorized as presenting high, moderate or low risks. • Analysis of parameter sensitivity and examples of the method’s application are given. • The method has value to regulators and those considering redeveloping closed landfills. - Abstract: A screening tool for quantifying levels of concern for contaminants detected in monitoring wells on or near landfills to down-gradient receptors (streams, wetlands and residential lots) was developed and evaluated. The tool uses Quick Domenico Multi-scenario (QDM), a spreadsheet implementation of Domenico-based solute transport, to estimate concentrations of contaminants reaching receptors under steady-state conditions from a constant-strength source. Unlike most other available Domenico-based model applications, QDM calculates the time for down-gradient contaminant concentrations to approach steady state and appropriate dispersivity values, and allows for up to fifty simulations on a single spreadsheet. Sensitivity of QDM solutions to critical model parameters was quantified. The screening tool uses QDM results to categorize landfills as having high, moderate and low levels of concern, based on contaminant concentrations reaching receptors relative to regulatory concentrations. The application of this tool was demonstrated by assessing levels of concern (as defined by the New Jersey Pinelands Commission) for thirty closed, uncapped landfills in the New Jersey Pinelands National Reserve, using historic water-quality data from monitoring wells on and near landfills and hydraulic parameters from regional flow models. Twelve of these landfills are categorized as having high levels of concern, indicating a need for further assessment. This tool is not a replacement for conventional numerically-based transport model or other available Domenico-based applications, but is suitable for quickly assessing the level of concern posed by a landfill or other contaminant point source before expensive and lengthy monitoring or remediation measures are taken. In addition to quantifying the level of concern using historic groundwater-monitoring data, the tool allows for archiving model scenarios and adding refinements as new data become available.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1988

Viscous/inviscid interaction procedure for high-amplitude oscillating airfoils

A coupling procedure between an inviscid potential theoretical panel method and an unsteady interactive boundary layer calculation procedure has been developed for calculating unsteady airloads on oscillating airfoils in the high-incidence regime. A coordinate transformation into a time-dependent surface coordinate has been incorporated into the boundary layer code to avoid a quasi-steady treatment of the viscous flow equations in the leading edge region. The upper and lower sides of the airfoil are calculated with the complete set of unsteady equations. The interactive region of the boundary layer calculation can be straightforwardly extended over the complete upper surface of the airfoil while accounting for interaction areas in the region of strong adverse pressure gradients with incipient separation. Comparisons of sample calculations with experimental data show reasonably good agreement even for unsteady drag coefficients.

Geissler, W.↗
DOE OSTI2022

Machine learning for impurity charge-state transition levels in semiconductors from elemental properties using multi-fidelity datasets

Quantifying charge-state transition energy levels of impurities in semiconductors is critical to understanding and engineering their optoelectronic properties for applications ranging from solar photovoltaics to infrared lasers. While these transition levels can be measured and calculated accurately, such efforts are time-consuming and more rapid prediction methods would be beneficial. Here, we significantly reduce the time typically required to predict impurity transition levels using multi-fidelity datasets and a machine learning approach employing features based on elemental properties and impurity positions. We use transition levels obtained from low-fidelity (i.e., local-density approximation or generalized gradient approximation) density functional theory (DFT) calculations, corrected using a recently proposed modified band alignment scheme, which well-approximates transition levels from high-fidelity DFT (i.e., hybrid HSE06). Further, the model fit to the large multi-fidelity database shows improved accuracy compared to the models trained on the more limited high-fidelity values. Crucially, in our approach, when using the multi-fidelity data, high-fidelity values are not required for model training, significantly reducing the computational cost required for training the model. Our machine learning model of transition levels has a root mean squared (mean absolute) error of 0.36 (0.27) eV vs high-fidelity hybrid functional values when averaged over 14 semiconductor systems from the II–VI and III–V families. As a guide for use on other systems, we assessed the model on simulated data to show the expected accuracy level as a function of bandgap for new materials of interest. Finally, we use the model to predict a complete space of impurity charge-state transition levels in all zinc blende III–V and II–VI systems.

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

Osmotic Pressure Gradient Effects on Water Diffusion in Porous Rock: Can They Pervert Permeability Tests?

Abstract Generation of a large network of hydraulic cracks is of key importance not only for the success of fracking of shale but also for the recent scheme of sequestration of CO2 in deep formations of basalt and peridotite, which are mafic and ultramafic rocks that combine chemically with CO2. In numerical simulation of the creation of a fracture network in porous rock, an important goal is to enhance the rock permeability. The objective of this article is to calculate the effect of osmotic pressure gradients caused by gradients of concentration of the ions of Ca, Mg, Na, etc. on the effective permeability of the rock. The basic differential equations are formulated, and their explicit solutions for appropriate initial and boundary conditions are obtained under certain plausible simplifications. The main result is explicit approximate formulas for the critical time before which no water permeation through a test specimen can be observed. Depending on various parameters, this time can be unacceptably long, which is manifested as a zero water outflow. The solution may also explain the unreasonably small permeability values reported for some shales.

Mechanics↗
NASA NTRS1974

The numerical calculation of laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

Klineberg, J. M.↗
NASA NTRS1974

On laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for a number of separated flows and the accuracy of the method is verified.

Klineberg, J. M.↗
DOE OSTI2016

Theoretical prediction of the band offsets at the ZnO/anatase TiO{sub 2} and GaN/ZnO heterojunctions using the self-consistent ab initio DFT/GGA-1/2 method

The band offsets of the ZnO/anatase TiO{sub 2} and GaN/ZnO heterojunctions are calculated using the density functional theory/generalized gradient approximation (DFT/GGA)-1/2 method, which takes into account the self-energy corrections and can give an approximate description to the quasiparticle characteristics of the electronic structure of semiconductors. We present the results of the ionization potential (IP)-based and interfacial offset-based band alignments. In the interfacial offset-based band alignment, to get the natural band offset, we use the surface calculations to estimate the change of reference level due to the interfacial strain. Based on the interface models and GGA-1/2 calculations, we find that the valence band maximum and conduction band minimum of ZnO, respectively, lie 0.64 eV and 0.57 eV above those of anatase TiO{sub 2}, while lie 0.84 eV and 1.09 eV below those of GaN, which agree well with the experimental data. However, a large discrepancy exists between the IP-based band offset and the calculated natural band offset, the mechanism of which is discussed. Our results clarify band alignment of the ZnO/anatase TiO{sub 2} heterojunction and show good agreement with the GW calculations for the GaN/ZnO heterojunction.

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

Magneto-electronic and Thermodynamic Properties of Quaternary NiFeMnZ (Z = In, Sn) New Spin Gapless Semiconductors

This paper deals with investigation of magneto-electronic and thermodynamic properties of new spin gapless semiconductor quaternary NiFeMnZ (Z = In, Sn) Heusler alloys. In addition, electronic properties such as band structure and density of states were studied in detail. The obtained results are based on first-principles calculations, using the full potential linearized augmented plane wave (FP-LAPW) method within well-known generalized gradient approximation (GGA) coupling with spin-orbit calculation (SOC) in framework of the density functional theory (DFT) implemented in the WIEN2k package. As results, the alloys were found to be new spin gapless semiconductors (SGSs), and both compounds in LiMgPdSn-type crystal structure are stable. The magnetic moment of NiFeMnIn and NiFeMnSn was 4.39 and 4.99 μ{sub B} per formula unit, respectively. This result is in agreement with the Slater-Pauling rule (M{sub tot} = 24 −Z{sub tot}). Furthermore, the thermodynamic properties were also predicted. All these multifunctionalities lead to crucial interest for developing future spintronics device generations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1975

Simple computer method provides contours for radiological images

Computer is provided with information concerning boundaries in total image. Gradient of each point in digitized image is calculated with aid of threshold technique; then there is invoked set of algorithms designed to reduce number of gradient elements and to retain only major ones for definition of contour.

Newell, J. D.↗
DOE OSTI2015

Profile stiffness measurements in the Helically Symmetric experiment and comparison to nonlinear gyrokinetic calculations

Stiffness measurements are presented in the quasi-helically symmetric experiment (HSX), in which the neoclassical transport is comparable to that in a tokamak and turbulent transport dominates throughout the plasma. Electron cyclotron emission is used to measure the local electron temperature response to modulated electron cyclotron resonant heating. The amplitude and phase of the heat wave through the steep electron temperature gradient (ETG) region of the plasma are used to determine a transient electron thermal diffusivity that is close to the steady-state diffusivity. The low stiffness in the region between 0.2 ≤ r/a ≤ 0.4 agrees with the scaling of the steady-state heat flux with temperature gradient in this region. These experimental results are compared to gyrokinetic calculations in a flux-tube geometry using the gyrokinetic electromagnetic numerical experiment code with two kinetic species. Linear simulations show that the ETG mode may be experimentally relevant within r/a ≤ 0.2, while the Trapped Electron Mode (TEM) is the dominant long-wavelength microturbulence instability across most of the plasma. The TEM is primarily driven by the density gradient. Non-linear calculations of the saturated heat flux driven by the TEM and ETG bracket the experimental heat flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Behavior of lithium ions in the turbulent near-wall tokamak plasma under heating of ions and electrons of the main plasma

Turbulent dynamics of the near-wall tokamak plasma is simulated by numerically solving the nonlinear reduced Braginskii magnetohydrodynamic equations with allowance for a lithium ion admixture. The effects of turbulence and radiation of the admixture are analyzed in the framework of a self-consistent approach. The radial distributions of the radiative loss power and the density of Li{sup 0} atoms and Li{sup +1} ions are obtained as functions of the electron and ion temperatures of the main plasma in the near-wall layer. The results of numerical simulations show that supply of lithium ions into the low-temperature near-wall plasma substantially depends on whether the additional power is deposited into the electron or ion component of the main plasma. If the electron temperature in the layer increases (ECR heating), then the ion density drops. At the same time, an increase in the temperature of the main ions (ICR heating) leads to an increase in the density of Li{sup +1} ions. The results of numerical simulations are explained by the different influence of the electron and ion temperatures on the atomic processes governing the accumulation and loss of particles in the balance equations for neutral Li{sup 0} atoms and Li{sup +1} ions in the admixture. The radial profile of the electron temperature and the corresponding distribution of the radiative loss power for different densities of neutral Li{sup 0} atoms on the wall are obtained. The calculations show that the presence of Li{sup +1} ions affects turbulent transport of the main ions. In this case, the electron heat flux increases by 20–30% with increasing Li{sup +1} density, whereas the flux of the main ions drops by nearly the same amount. The radial profile of the turbulent flux of lithium ions is obtained. It is demonstrated that the appearance of the pinch effect is related to the positive density gradient of lithium ions across the calculation layer. For the parameters of the T-10 tokamak, the effect of radiative cooling of the near-wall plasma layer becomes appreciable when the near-wall density of neutral lithium atoms exceeds 7 × 10{sup 11} cm{sup −3}. In this case, the density of radiative loss power in the center of the layer is estimated to be about 500–600 kW/m{sup 3}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

SU-G-BRB-05: Automation of the Photon Dosimetric Quality Assurance Program of a Linear Accelerator

Purpose: To develop an automated method to calculate a linear accelerator (LINAC) photon radiation field size, flatness, symmetry, output and beam quality in a single delivery for flattened (FF) and flattening-filter-free (FFF) beams using an ionization chamber array. Methods: The proposed method consists of three control points that deliver 30×30, 10×10 and 5×5cm{sup 2} fields (FF or FFF) in a step-and-shoot sequence where the number of monitor units is weighted for each field size. The IC Profiler (Sun Nuclear Inc.) with 5mm detector spacing was used for this study. The corrected counts (CCs) were calculated and the locations of the maxima and minima values of the first-order gradient determined data of each sub field. Then, all CCs for each field size are summed in order to obtain the final profiles. For each profile, the radiation field size, symmetry, flatness, output factor and beam quality were calculated. For field size calculation, a parameterized gradient method was used. For method validation, profiles were collected in the detector array both, individually and as part of the step-and-shoot plan, with 9.9cm buildup for FF and FFF beams at 90cm source-to-surface distance. The same data were collected with the device (plus buildup) placed on a movable platform to achieve a 1mm resolution. Results: The differences between the dosimetric quantities calculated from both deliveries, individually and step-and-shoot, were within 0.31±0.20% and 0.04±0.02mm. The differences between the calculated field sizes with 5mm and 1mm resolution were ±0.1mm. Conclusion: The proposed single delivery method proved to be simple and efficient in automating the photon dosimetric monthly and annual quality assurance.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

A shortcut to gradient-corrected magnon dispersion: exchange-only case

Ab initio calculations of the magnon dispersion in ferromagnetic materials typically rely on the adiabatic local density approximation (ALDA) in which the effective exchange-correlation field is everywhere parallel to the magnetization. These calculations, however, tend to overestimate the “magnon stiffness”, defined as the curvature of the magnon frequency vs. wave vector relation evaluated at zero wave vector. Here we suggest a simple procedure to improve the magnon dispersion by taking into account gradient corrections to the ALDA at the exchange-only level. We find that this gradient correction always reduces the magnon stiffness. The surprisingly large size of these corrections (~30%) greatly improves the agreement between the calculated and the observed magnon stiffness for cobalt and nickel, which are known to be overestimated within the ALDA.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Engineering calculations for communications satellite systems planning

An extended gradient search code for broadcasting satellite service (BSS) spectrum/orbit assignment synthesis is discussed. Progress is also reported on both single-entry and full synthesis computational aids for fixed satellite service (FSS) spectrum/orbit assignment purposes.

Levis, C. A.↗