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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

DOE OSTI2016

West Flank Coso FORGE: 3D Temperature Model

x,y,z data of the 3D temperature model for the West Flank Coso FORGE site. Model grid spacing is 250m. The temperature model for the Coso geothermal field used over 100 geothermal production sized wells and intermediate-depth temperature holes. At the near surface of this model, two boundary temperatures were assumed: (1) areas with surface manifestations, including fumaroles along the northeast striking normal faults and northwest striking dextral faults with the hydrothermal field, a temperature of ~104 deg C was applied to datum at +1066 meters above sea level elevation, and (2) a near-surface temperature at about 10 meters depth, of 20 deg C was applied below the diurnal and annual conductive temperature perturbations. These assumptions were based on heat flow studies conducted at the CVF and for the Mojave Desert. On the edges of the hydrothermal system, a 73 deg C/km temperature gradient contour was established using conductive gradient data from shallow and intermediate-depth temperature holes. This contour was continued to all elevation datums between the 20 deg C surface and -1520 meters below mean sea level. Because the West Flank is outside of the geothermal field footprint, during Phase 1, the three wells inside the FORGE site were incorporated into the preexisting temperature model. To ensure a complete model was built based on all the available data sets, measured bottom-hole temperature gradients in certain wells were downward extrapolated to the next deepest elevation datum (or a maximum of about 25% of the well depth where conductive gradients are evident in the lower portions of the wells). After assuring that the margins of the geothermal field were going to be adequately modelled, the data was contoured using the Kriging method algorithm. Although the extrapolated temperatures and boundary conditions are not rigorous, the calculated temperatures are anticipated to be within ~6 deg C (20 deg F), or one contour interval, of the observed data within the Coso geothermal field. Based on a lack of temperature data west of 74-2TCH, the edges of this model still seem to have an effect on West Flank modeled temperatures.

15 GEOTHERMAL ENERGY↗

DOE OSTI2016

Modification and validation of an analytical source model for external beam radiotherapy Monte Carlo dose calculations

Purpose: A dose calculation tool, which combines the accuracy of the dose planning method (DPM) Monte Carlo code and the versatility of a practical analytical multisource model, which was previously reported has been improved and validated for the Varian 6 and 10 MV linear accelerators (linacs). The calculation tool can be used to calculate doses in advanced clinical application studies. One shortcoming of current clinical trials that report dose from patient plans is the lack of a standardized dose calculation methodology. Because commercial treatment planning systems (TPSs) have their own dose calculation algorithms and the clinical trial participant who uses these systems is responsible for commissioning the beam model, variation exists in the reported calculated dose distributions. Today’s modern linac is manufactured to tight specifications so that variability within a linac model is quite low. The expectation is that a single dose calculation tool for a specific linac model can be used to accurately recalculate dose from patient plans that have been submitted to the clinical trial community from any institution. The calculation tool would provide for a more meaningful outcome analysis. Methods: The analytical source model was described by a primary point source, a secondary extra-focal source, and a contaminant electron source. Off-axis energy softening and fluence effects were also included. The additions of hyperbolic functions have been incorporated into the model to correct for the changes in output and in electron contamination with field size. A multileaf collimator (MLC) model is included to facilitate phantom and patient dose calculations. An offset to the MLC leaf positions was used to correct for the rudimentary assumed primary point source. Results: Dose calculations of the depth dose and profiles for field sizes 4 × 4 to 40 × 40 cm agree with measurement within 2% of the maximum dose or 2 mm distance to agreement (DTA) for 95% of the data points tested. The model was capable of predicting the depth of the maximum dose within 1 mm. Anthropomorphic phantom benchmark testing of modulated and patterned MLCs treatment plans showed agreement to measurement within 3% in target regions using thermoluminescent dosimeters (TLD). Using radiochromic film normalized to TLD, a gamma criteria of 3% of maximum dose and 2 mm DTA was applied with a pass rate of least 85% in the high dose, high gradient, and low dose regions. Finally, recalculations of patient plans using DPM showed good agreement relative to a commercial TPS when comparing dose volume histograms and 2D dose distributions. Conclusions: A unique analytical source model coupled to the dose planning method Monte Carlo dose calculation code has been modified and validated using basic beam data and anthropomorphic phantom measurement. While this tool can be applied in general use for a particular linac model, specifically it was developed to provide a singular methodology to independently assess treatment plan dose distributions from those clinical institutions participating in National Cancer Institute trials.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗

DOE OSTI2019

Interaction of Strain Refined Precipitates and Recrystallized Grains in Nb-Ti Microalloyed Steel during Continuous Casting, Hot-Core Heavy Reduction Rolling, and Reheating Process

Hot-core heavy reduction rolling (HHR{sup 2}) is an energy-saving process that directly rolls continuously cast slabs at the end of solidification that have a large temperature gradient. Here, we elucidated the interaction of strain refined precipitates and recrystallized grains in Nb-Ti microalloyed steel during continuous casting, hot-core heavy reduction rolling, and reheating process. Furthermore, theoretical models were established to calculate recrystallization driving force and pinning force at temperatures in the range of 850-1300 °C by a two-stage interrupted compression test using a Gleeble-3800 thermomechanical simulator. The study indicated that a random dispersion of strain refined precipitates with a size range of ~ 6-10 nm was obtained in the hot-core heavy reduction rolled slab, which was different from the coarse and uneven distribution of precipitates in a conventional slab. Upon reheating, the austenite grains in the hot-core heavy reduction rolled slab were significantly refined and evenly distributed along the thickness direction. This is attributed to the random dispersion of strain refined precipitates formed in the HHR{sup 2} process.

36 MATERIALS SCIENCE↗

DOE OSTI2015

SU-E-T-536: LINAC-Based Single Isocenter Frameless SRT for Brain Metastases

Purpose: Single-isocenter Stereotactic Radiotherapy of multiple brain metastases with Varian 21 IX LINAC, using Aktina Pinpoint system for patient setup. Methods: In 2014, five single-isocenter RapidArc SRT plans were delivered to five patients with 2 to 8 brain metastases using Varian 21 IX. Aktina Pinpoint system was used for setup and 2mm PTV margin were used. CBCT was acquired before and after the beam delivery. The prescription is 2100 cGy in 3 fractions. Eclipse planning system was used for treatment planning. Depending on the number of metastases and their locations, 1 to 5 coplanar or non coplanar arcs were used. Typically, 2 or 3 arcs are used. IMRT QAs were performed by comparing an A1SL ion chamber point dose measurement in solid water phantom to point dose of the plan; also, based on EPID measurement, 3D spatial dose was calculated using DosimetryCheck software package from MathResolutions Inc. The EPID system has an active area of 40cm by 30cm with 1024 by 768 photodiodes, which corresponds to a resolution of 0.4mm by 0.4mm pixel dimension. Results: for all the plans, at least 95% PTV coverage was achieved for full prescription dose, with plan normalization > 75%. RTOG conformity indices are less than 1.1 and Paddick gradient indices are less than 4.5. The distance from prescription IDL to 50% IDL increases as the number of metastases increases, and it ranges from 0.6mm to 0.8mm. Treatment time varies from 10mins to 30mins, depending on the number of arcs and if the arcs are coplanar. IMRT QA shows that the ion chamber measurement agree with the eclipse calculation within 3%, and 95% of the points passed Gamma, using 3% dose difference and 3mm DTA Conclusion: High quality single isocenter RapidArc SRT plan can be optimized and accurately delivered using Eclipse and Varian 21IX.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2023

A perfused multi-well bioreactor platform to assess tumor organoid response to a chemotherapeutic gradient

There is an urgent need to develop new therapies for colorectal cancer that has metastasized to the liver and, more fundamentally, to develop improved preclinical platforms of colorectal cancer liver metastases (CRCLM) to screen therapies for efficacy. To this end, we developed a multi-well perfusable bioreactor capable of monitoring CRCLM patient-derived organoid response to a chemotherapeutic gradient. CRCLM patient-derived organoids were cultured in the multi-well bioreactor for 7 days and the subsequently established gradient in 5-fluorouracil (5-FU) concentration resulted in a lower IC 50 in the region near the perfusion channel versus the region far from the channel. We compared behaviour of organoids in this platform to two commonly used PDO culture models: organoids in media and organoids in a static (no perfusion) hydrogel. The bioreactor IC 50 values were significantly higher than IC 50 values for organoids cultured in media whereas only the IC 50 for organoids far from the channel were significantly different than organoids cultured in the static hydrogel condition. Using finite element simulations, we showed that the total dose delivered, calculated using area under the curve (AUC) was similar between platforms, however normalized viability was lower for the organoid in media condition than in the static gel and bioreactor. Our results highlight the utility of our multi-well bioreactor for studying organoid response to chemical gradients and demonstrate that comparing drug response across these different platforms is nontrivial.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Theoretical Investigation of Half-Metallic Ferromagnetism in Sodium-Based Fluoro-perovskite NaXF{sub 3} (X = V, Co)

We report on the structural, mechanical, magneto-electronic and thermodynamic properties of cubic sodium based fluoro-perovskite NaXF3 (X = V, Co) through density functional theory (DFT). The generalized gradient approximation of Perdew, Burke, and Ernzerh (GGA-PBE) and the Tran-Blaha modified Becke-Johnson potential (TB-mBJ) are used for modeling exchange-correlation effects. The calculated lattice parameters agree well with the experimental measurements. Elastic properties show that these compounds are mechanically stable. The results of spin-polarized band structure show a half metallic behavior of NaVF{sub 3} and NaCoF{sub 3}. The magnetic study shows that the nature of these compounds is ferromagnetic and the integer value of the total magnetic moment confirms the half metallicity for both materials. The thermodynamic characters are also investigated through quasi-harmonic approximation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Assessing dual continuum method for multicomponent reactive transport

The subsurface flow and reactive transport code PFLOTRAN can accommodate fracture–matrix interaction for multicomponent reactive transport using a Dual Continuum formulation. In this work we present two examples to assess this approach. Through these examples, we also illustrate how this Dual Continuum model can be parameterised and what kind of grid spacing is needed to accurately capture reactive transport in the secondary continuum. The results of the calculations show that when geochemical reactions in the matrix (or the secondary continuum) are included, a very fine grid spacing close to the fracture–matrix interface is required to properly capture geochemical gradients at the interface. We also show that one can reduce the computational cost by using a variable grid spacing, that limits the total number in the secondary continuum, and yet obtain accurate results. Furthermore, this feature makes the Dual Continuum approach suitable for large-scale reactive transport modelling in fractured crystalline rock.

58 GEOSCIENCES↗

DOE OSTI2015

A first-principles study of cementite (Fe{sub 3}C) and its alloyed counterparts: Elastic constants, elastic anisotropies, and isotropic elastic moduli

A comprehensive computational study of elastic properties of cementite (Fe{sub 3}C) and its alloyed counterparts (M{sub 3}C (M = Al, Co, Cr, Cu, Fe, Hf, Mn, Mo, Nb, Ni, Si, Ta, Ti, V, W, Zr, Cr{sub 2}FeC and CrFe{sub 2}C) having the crystal structure of Fe{sub 3}C is carried out employing electronic density-functional theory (DFT), all-electron PAW pseudopotentials and the generalized gradient approximation for the exchange-correlation energy (GGA). Specifically, as a part of our systematic study of cohesive properties of solids and in the spirit of materials genome, following properties are calculated: (i) single-crystal elastic constants, C{sub ij}, of above M{sub 3}Cs; (ii) anisotropies of bulk, Young’s and shear moduli, and Poisson’s ratio based on calculated C{sub ij}s, demonstrating their extreme anisotropies; (iii) isotropic (polycrystalline) elastic moduli (bulk, shear, Young’s moduli and Poisson’s ratio) of M{sub 3}Cs by homogenization of calculated C{sub ij}s; and (iv) acoustic Debye temperature, θ{sub D}, of M{sub 3}Cs based on calculated C{sub ij}s. We provide a critical appraisal of available data of polycrystalline elastic properties of alloyed cementite. Calculated single crystal properties may be incorporated in anisotropic constitutive models to develop and test microstructure-processing-property-performance links in multi-phase materials where cementite is a constituent phase.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Ab-initio study of electronic structure and elastic properties of ZrC

The electronic and elastic properties of ZrC have been investigated using the linear combination of atomic orbitals method within the framework of density functional theory. Different exchange-correlation functionals are taken into account within generalized gradient approximation. We have computed energy bands, density of states, elastic constants, bulk modulus, shear modulus, Young’s modulus, Poisson’s ratio, lattice parameters and pressure derivative of the bulk modulus by calculating ground state energy of the rock salt structure type ZrC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Direct measurements of the 3D plasma velocity in single-helical-axis RFP plasmas

The first local velocity measurements of helical equilibrium plasmas in the Reversed Field Pinch (RFP) Single Helical Axis (SHAx) state using a Charge Exchange Recombination Spectroscopy (CHERS) diagnostic are presented. Measurements show strong axisymmetric and non-axisymmetric flow, with n ¼ 5 components of flow related to the ðm; nÞ¼ð1; 5Þ dominant magnetic mode on the order of the axisymmetric flow in certain regions of the plasma, as well as significant n > 5 flow. Flow measurements are compared with NIMROD simulations of visco-resistive, single-fluid MHD in toroidal and cylindrical geometries with limited axial periodicity. Both measurements and the simulation with toroidal geometry show stronger inboard flows relative to the outboard flows, which is attributed to the toroidal geometry of the device. In the experiment, the n ¼ 5 component of flow is phase shifted from the reconnection-like flow pattern observed in the single-fluid simulations, possibly due to decoupling of the ion and electron fluids over much of the plasma. Finally, the strength of the helical angular flow shear relative to the critical shear necessary to disrupt nonlinear coupling between tearing modes is calculated around the helical magnetic axis. The shear in the measured flow is on the order of the theoretical critical threshold needed to nonlinearly decouple modes, but the measurement uncertainty in the gradient of the flow is large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Phase-field-based calculations of the disregistry fields of static extended dislocations in FCC metals

In the continuum context, the displacements of atoms induced by a dislocation can be approximated by a continuum disregistry field. In this paper, two phase-field (PF)-based approaches and their variants are employed to calculate the disregistry fields of static, extended dislocations of pure edge and pure screw character in two face-centred cubic metals: Au and Al, which have distinct stable stacking fault energy and elastic anisotropy. A new truncated Fourier series form is developed to approximate the generalised stacking fault energy (GSFE) surface, which shows significant improvement over the previously employed Fourier series form. By measuring the intrinsic stacking fault (ISF) width and partial dislocation core size in different ways, the PF-based disregistry fields are quantitatively compared against those predicted by molecular statics. In particular, two new measures for the ISF widths are proposed and shown to overcome drawbacks of the more commonly used standards in the literature. Our calculations also show that continuum formulation of the elastic energy and the GSFE for a homogeneous surface can successfully characterise the core structure. Finally, our comparisons highlight the significance of including the gradient energy in the free energy formulation when an accurate description of the dislocation core structure is desired.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Ground state properties of 3d metals from self-consistent GW approach

The self consistent GW approach (scGW) has been applied to calculate the ground state properties (equilibrium Wigner–Seitz radius S WZ and bulk modulus B) of 3d transition metals Sc, Ti, V, Fe, Co, Ni, and Cu. The approach systematically underestimates S WZ with average relative deviation from the experimental data of about 1% and it overestimates the calculated bulk modulus with relative error of about 25%. We show that scGW is superior in accuracy as compared to the local density approximation but it is less accurate than the generalized gradient approach for the materials studied. If compared to the random phase approximation, scGW is slightly less accurate, but its error for 3d metals looks more systematic. Lastly, the systematic nature of the deviation from the experimental data suggests that the next order of the perturbation theory should allow one to reduce the error.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

First-principles study of an S = 1 quasi one-dimensional quantum molecular magnetic material

We use density functional theory to study the structural, magnetic, and electronic structures of the organometallic quantum magnet NiCl 2-4 SC (NH 2 ) 2 (DTN). Recent work has demonstrated the quasi one-dimensional nature of the molecular crystal and studied its quantum phase transitions at low temperatures. The system includes a magnetoelectric (ME) coupling and, when doped with Br, the presence of an exotic Bose-glass state. Using the generalized gradient approximation with inclusion of a van der Waals term to account for weak intermolecular forces and by introducing a Hubbard U term to the total energy, we systematically show that our calculations reproduce the magnetic anisotropy, the intermolecular exchange coupling strength, and the magnetoelectric effect in DTN, which have been observed in previous experiments. Further analysis of the electronic structure gives insight into the underlying magnetic interactions, including what mechanisms may be causing the ME effect. Using this computationally efficient model, we predict what effect applying an electric field might have on the magnetic properties of this quantum magnet.

1-dimensional spin chains↗

NASA NTRS1973

Minimum impulse three-body trajectories.

A rapid and accurate method of calculating optimal impulsive transfers in the restricted problem of three bodies has been developed. The technique combines a multi-conic method of trajectory integration with primer vector theory and an accelerated gradient method of trajectory optimization. A unique feature is that the state transition matrix and the primer vector are found analytical without additional integrations or differentiations. The method has been applied to the determination of optimal two and three impulse transfers between the L2 libration point and circular orbits about both the earth and the moon.

D'Amario, L.↗

NASA NTRS1980

Plasma collection by high voltage spacecraft at low earth orbit

A computer model of the three-dimensional sheath formation and plasma current collection by high voltage spacecraft has been developed. By using new space charge density and plasma collection algorithms, it is practical to perform calculations for large, complex spacecraft. The model uses NASCAP compatible objects and geometries. Results indicate that ion focusing observed in the laboratory during high voltage collection experiments is probably due to voltage gradients on the collecting surfaces.

Katz, I.↗

NASA NTRS1980

Numerical models of solar modulation of galactic cosmic rays, including drifts

Numerical solutions are presented of the steady-state modulation of galactic cosmic rays for a realistic two-dimensional model which includes drifts. Calculations of energy spectra and the spatial density variation of protons are shown for a model that incorporates the solar-minimum magnetic field and nominal interplanetary propagation parameters. The calculated values are relatively insensitive to the assumed size of the modulating region and to the magnitude of the particle diffusion coefficient. The solutions show a broad interior 'plateau' in which radial gradients are small. A simple model that neglects diffusion but includes drifts and adiabatic cooling accounts for these results. It is concluded that in this model and the range of parameters, diffusion plays a role secondary to that of drifts and adiabatic cooling.

Jokipii, J. R.↗

NASA NTRS1982

Features of a reattaching turbulent shear layer subject to an adverse pressure gradient

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step with superimposed adverse pressure gradient. Mean velocities, Reynolds stresses and triple-products measured by a laser Doppler velocimeter are presented for two cases of adverse pressure gradient. Mixing lengths, eddy viscosities, production, convection, turbulent diffusion, and dissipation terms are extracted from the data. These data are compared with various mixing length and eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and the algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model

Driver, D. M.↗