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

Ion temperature gradient turbulence in helical and axisymmetric RFP plasmas

Turbulence induced by the ion temperature gradient (ITG) is investigated in the helical and axisymmetric plasma states of a reversed field pinch device by means of gyrokinetic calculations. The two magnetic configurations are systematically compared, both linearly and nonlinearly, in order to evaluate the impact of the geometry on the instability and its ensuing transport, as well as on the production of zonal flows. Despite its enhanced confinement, the high-current helical state demonstrates a lower ITG stability threshold compared to the axisymmetric state, and ITG turbulence is expected to become an important contributor to the total heat transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Comparing measurements of limiting current of electrolytes with theoretical predictions up to the solubility limit

Imposing a steady ionic current through an electrolyte results in the formation of salt concentration gradients that compromise battery performance. The limiting current is usually defined as the current at which the salt concentration at the cathode approaches zero. Higher currents cannot be imposed on the cell as larger concentration gradients are unsustainable. We study the limiting current in electrolytes comprising a perfluorinated oligomer, C8-DMC, and lithium bis-(fluorosulfonyl)imide salt in symmetric lithium cells. The time-dependence of the potential, which increases as salt concentration gradients develop, was also measured. Both steady-state and transient behaviors are modeled using Newman's concentrated solution theory; transport and thermodynamic parameters needed to perform the calculations were measured independently. The limiting current is a nonmonotonic function of salt concentration in both theory and experiment. The model shows that at low salt concentrations (below 0.88 mol/kg solvent), the concentration at the cathode approaches zero at limiting current. In contrast, at high salt concentrations (above 0.88 mol/kg solvent), the concentration at the anode approaches the solubility limit (2.03 mol/kg solvent). The experimentally determined salt concentration at which the limiting current is maximized is in excellent agreement with theoretical predictions made without resorting to any adjustable parameters.

Shah, Deep B.↗
DOE OSTI2024

Validation of Hole-Drilling Residual Stress Measurements in Workpieces of Various Thickness

A recent revision to the ASTM E837 standard for near-surface residual stress measurement by the hole-drilling method describes a new thickness-dependent stress calculation procedure applicable to “thin” and “intermediate” workpieces for which strain versus depth response depends on workpiece thickness. This new calculation procedure differs from that of the prior standard, which applies only to thick workpieces with strain versus depth response independent of thickness. Herein we assess the new calculation procedures by performing hole-drilling residual stress measurements in samples with a range of thickness. Near-surface residual stress is measured in a thick aluminum plate containing near-surface residual stress from a uniform shot peening treatment, and in samples of different thickness removed from the plate at the peened surface. A finite element (FE) model is used to assess consistency between measured residual stress across the range of sample thickness. Measured residual stress varies with sample thickness, with thinner samples exhibiting smaller near-surface compressive stress and a larger gradient of subsurface stress. These trends are consistent with both observed bending (curvature) of the removed samples and the trend in FE-calculated expected residual stress. The measured and expected residual stresses are in good agreement for samples of intermediate thickness, but the agreement decreases with sample thickness. Measured residual stress is invariant with gage circle diameter. The new thickness-dependent stress calculation procedure for hole-drilling provides meaningful improvement compared to thick-workpiece calculations.

42 ENGINEERING↗
DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

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

Theory and observations of upward field‐aligned currents at the magnetopause boundary layer

Abstract The dependence of the upward field‐aligned current density ( J || ) at the dayside magnetopause boundary layer is well described by a simple analytic model based on a velocity shear generator. A previous observational survey confirmed that the scaling properties predicted by the analytical model are applicable between 11 and 17 MLT. We utilize the analytic model to predict field‐aligned currents using solar wind and ionospheric parameters and compare with direct observations. The calculated and observed parallel currents are in excellent agreement, suggesting that the model may be useful to infer boundary layer structures. However, near noon, where velocity shear is small, the kinetic pressure gradients and thermal currents, which are not included in the model, could make a small but significant contribution to J || . Excluding data from noon, our least squares fit returns log( J ||,max_cal ) = (0.96 ± 0.04) log( J ||_obs ) + (0.03 ± 0.01) where J ||,max_cal = calculated J ||,max and J ||_obs = observed J || .

Wing, Simon↗
DOE OSTI2016

PROTEUS-SN User Manual

PROTEUS-SN is a three-dimensional, highly scalable, high-fidelity neutron transport code developed at Argonne National Laboratory. The code is applicable to all spectrum reactor transport calculations, particularly those in which a high degree of fidelity is needed either to represent spatial detail or to resolve solution gradients. PROTEUS-SN solves the second order formulation of the transport equation using the continuous Galerkin finite element method in space, the discrete ordinates approximation in angle, and the multigroup approximation in energy. PROTEUS-SN’s parallel methodology permits the efficient decomposition of the problem by both space and angle, permitting large problems to run efficiently on hundreds of thousands of cores. PROTEUS-SN can also be used in serial or on smaller compute clusters (10’s to 100’s of cores) for smaller homogenized problems, although it is generally more computationally expensive than traditional homogenized methodology codes. PROTEUS-SN has been used to model partially homogenized systems, where regions of interest are represented explicitly and other regions are homogenized to reduce the problem size and required computational resources. PROTEUS-SN solves forward and adjoint eigenvalue problems and permits both neutron upscattering and downscattering. An adiabatic kinetics option has recently been included for performing simple time-dependent calculations in addition to standard steady state calculations. PROTEUS-SN handles void and reflective boundary conditions. Multigroup cross sections can be generated externally using the MC2-3 fast reactor multigroup cross section generation code or internally using the cross section application programming interface (API) which can treat the subgroup or resonance table libraries. PROTEUS-SN is written in Fortran 90 and also includes C preprocessor definitions. The code links against the PETSc, METIS, HDF5, and MPICH libraries. It optionally links against the MOAB library and is a part of the SHARP multi-physics suite for coupled multi-physics analysis of nuclear reactors. This user manual describes how to set up a neutron transport simulation with the PROTEUS-SN code. A companion methodology manual describes the theory and algorithms within PROTEUS-SN.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Three-dimensional analysis of 0/90s and 90/0s laminates with a central circular hole

Stress distributions were calculated near a circular hole in laminates, using a three dimensional finite element analysis. These stress distributions were presented three ways: through the thickness at the hole boundary, along radial lines at the 0/90 and 90/0 interfaces, and around the hole at these interfaces. The interlaminar normal stress, and the shear stress, distributions had very steep gradients near the hole boundary, suggesting interlaminar stress singularities. The largest compressive stress occurred at about 60 deg from the load axis. A simple procedure was introduced to calculate interlaminar stresses near the hole boundary. It used stresses calculated by an exact two dimensional analysis of a laminate with a hole as input to a quasi three dimensional model. It produced stresses that agreed closely with those from the three dimensional finite element model.

Raju, I. S.↗
DOE OSTI2024

Evaporation barrier for dark matter in celestial bodies

The minimum testable dark matter (DM) mass for almost all DM signatures in celestial bodies is determined by the rate at which DM evaporates. DM evaporation has previously been calculated assuming a competition between the gravitational potential of the object, and thermal kicks from the celestial-body matter. We point out a new effect, where mediators with a range larger than the interparticle spacing induce a force proportional to the density gradient of celestial objects, forming an evaporation barrier for the DM. This effect can be so significant that evaporation does not occur even for sub-MeV DM, in stark contrast to previous calculations. This opens up a wide range of new light DM searches, many orders of magnitude in DM mass below the sensitivity of direct detection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Blind Modeling Validation Exercises Using the Horizontal Dry Cask Simulator

The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). A dry cask simulator (DCS) was built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the accuracy of the modeling used to predict cladding temperatures. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. A model validation exercise was previously carried out for the DCS in a vertical configuration. Lessons learned during the previous validation exercise have been applied to a new, blind study using a horizontal dry cask simulator (HDCS). Three modeling institutions – the Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) – were granted access to the input parameters from the DCS Handbook, SAND2017-13058R, and results from a limited data set from the horizontal BWR dry cask simulator tests reported in the HDCS update report, SAND2019-11688R. With this information, each institution was tasked to calculate peak cladding temperatures and air mass flow rates for ten HDCS test cases. Axial as well as vertical and horizontal transverse temperature profiles were also calculated. These calculations were done using modeling codes (ANSYS/Fluent, STAR-CCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, the ten test cases of the horizontal dry cask simulator were defined by three independent variables – fuel assembly decay heat (0.5 kW, 1 kW, 2.5 W, and 5 kW), internal backfill pressure (100 kPa and 800 kPa), and backfill gas (helium and air). The plots provided in Chapter 3 of this report show the axial, vertical, and horizontal temperature profiles obtained from the dry cask simulator experiments in the horizontal configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), PCT axial locations, axial temperature profiles, vertical and horizontal temperature profiles at two different axial locations, and air mass flow rates for the ten test cases, normalized by the experimental results. The model results are assigned arbitrary model numbers to retain anonymity. Due to the relatively flat axial temperature profiles, small temperature gradients resulted in large deviations of all models’ PCT axial location from the experimental PCT axial location. When the PCT axial location error is excluded in the calculation of the combined RMS of the normalized errors that considers PCT, the temperature profiles, and the air mass flow rates, the model data fits the experimental data to within 5%. When the vault information is excluded, the model data fits the experimental data to within 2.5%. An error analysis was developed further for one model, using the model and experimental uncertainties in each validation parameter to calculate validation uncertainties. The uncertainties for each parameter were used to define quantifiable validation criteria. For this analysis, the model was considered validated for a given comparison metric if the normalized error in that metric divided by the validation uncertainty was less than or equal to 1. When considering the combined RMS of the normalized errors of all metrics divided by their validation uncertainties, the model was found to have satisfied the criterion for model validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

TH-A-19A-12: A GPU-Accelerated and Monte Carlo-Based Intensity Modulated Proton Therapy Optimization System

Purpose: To develop a clinically applicable intensity modulated proton therapy (IMPT) optimization system that utilizes more accurate Monte Carlo (MC) dose calculation, rather than analytical dose calculation. Methods: A very fast in-house graphics processing unit (GPU) based MC dose calculation engine was employed to generate the dose influence map for each proton spot. With the MC generated influence map, a modified gradient based optimization method was used to achieve the desired dose volume histograms (DVH). The intrinsic CT image resolution was adopted for voxelization in simulation and optimization to preserve the spatial resolution. The optimizations were computed on a multi-GPU framework to mitigate the memory limitation issues for the large dose influence maps that Result from maintaining the intrinsic CT resolution and large number of proton spots. The dose effects were studied particularly in cases with heterogeneous materials in comparison with the commercial treatment planning system (TPS). Results: For a relatively large and complex three-field bi-lateral head and neck case (i.e. >100K spots with a target volume of ∼1000 cc and multiple surrounding critical structures), the optimization together with the initial MC dose influence map calculation can be done in a clinically viable time frame (i.e. less than 15 minutes) on a GPU cluster consisting of 24 Nvidia GeForce GTX Titan cards. The DVHs of the MC TPS plan compare favorably with those of a commercial treatment planning system. Conclusion: A GPU accelerated and MC-based IMPT optimization system was developed. The dose calculation and plan optimization can be performed in less than 15 minutes on a hardware system costing less than 45,000 dollars. The fast calculation and optimization makes the system easily expandable to robust and multi-criteria optimization. This work was funded in part by a grant from Varian Medical Systems, Inc.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1951

Determination and Use of the Local Recovery Factor for Calculating the Effectiveness Gas Temperature for Turbine Blades / Jack B. Esgar and Alfred L. Lea

In an in experimental investigation of local recovery factors for a blade having a pressure distribution similar to that of a typical reaction-type turbine blade, it a was found that the recovery factors were essentially independent of Mach number, Reynolds number, pressure gradient, and position on the blade surface except for regions where the boundary layer was probably in the transition range from laminar to turbulent. The recommended value of local subsonic recovery factor for use in calculating the effective gas temperature for gas turbine blades was 0.89.

HEAT TRANSFER, CASCADES↗
NASA NTRS1987

NASCRAC - A computer code for fracture mechanics analysis of crack growth

NASCRAC - a computer code for fracture mechanics analysis of crack growth - is described in this paper. The need for such a code is increasing as requirements grow for high reliability and low weight in aerospace components. The code is comprehensive and versatile, as well as user friendly. The major purpose of the code is calculation of fatigue, corrosion fatigue, or stress corrosion crack growth, and a variety of crack growth relations can be selected by the user. Additionally, crack retardation models are included. A very wide variety of stress intensity factor solutions are contained in the code, and extensive use is made of influence functions. This allows complex stress gradients in three-dimensional crack problems to be treated easily and economically. In cases where previous stress intensity factor solutions are not adequate, new influence functions can be calculated by the code. Additional features include incorporation of J-integral solutions from the literature and a capability for estimating elastic-plastic stress redistribution from the results of a corresponding elastic analysis. An example problem is presented which shows typical outputs from the code.

Harris, D. O.↗
DOE OSTI2015

The effect of shear flow on the rotational diffusion of a single axisymmetric particle

Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$ , where $p$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $p$, as $1/D_{0}^{r}p^{4}$ and as $1/D_{0}^{r}p^{2}$ for $p\gg$ 1. For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of ${\approx}2$ and the rotational diffusion by a factor of ${\approx}40$.

42 ENGINEERING↗
DOE OSTI2025

Measurement of phonon angular momentum

In condensed matter, angular momentum is intimately related to the emergence of topological quantum states, including chiral superconductivity, quantum spin liquids and various chiral quasiparticles. Recently, it has been predicted that microscopic lattice excitations, known as phonons, can carry finite angular momentum, leading to specific physical properties of materials. However, phonon angular momentum has not yet been observed directly. Here, in this study, we demonstrate that angular momentum conservation results in a macroscopic mechanical torque when applying a time-reversal symmetry-breaking thermal gradient along the chiral axis of single-crystal tellurium. We probe this torque using a cantilever-based device and establish that it changes sign by flipping the chirality or thermal gradient. This behavior disappears in polycrystalline samples that lack a preferred chirality. Our experimental results align well with theoretical calculations. We provide compelling evidence for phonon angular momentum, which might enable quantum states with potential applications in microelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Local nonlinear rf forces in inhomogeneous magnetized plasmas

The local nonlinear forces induced by radio frequency (rf) waves are derived in inhomogeneous magnetized plasmas, where the inhomogeneity exists in the rf fields, in the static magnetic field as well as in the equilibrium density and temperature. The local parallel force is completely resonant, but a novel component dependent on those inhomogeneities is obtained as the result of the inhomogeneous transport of parallel resonant-absorbed momentum by the nonlinear perpendicular drift flux. In the local poloidal force, the component induced by the inhomogeneity of rf power absorption is also confirmed and it can be recognized as the residual effect from the incomplete cancellation between the rate of the diamagnetic poloidal momentum gain and the Lorentz force due to the radial diffusion-like flux. The compact expression for radial force is also obtained for the first time, whose nonresonant component is expressed as the sum of the ponderomotive force on particles and the gradients of the nonresonant perpendicular pressure and of the nonresonant momentum flux due to the finite temperature effect. Numerical calculations in a 1-D slab model show that the resonant component dependent on the inhomogeneities may be significant when the ion absorption dominates the resonant wave-particle interaction. A quantitative estimation shows that the novel component in the parallel force is important to understand the experiments of the ion-cyclotron-frequency mode-conversion flow drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Microtearing modes in tokamak discharges

Microtearing modes (MTMs) have been identified as a source of significant electron thermal transport in tokamak discharges. In order to describe the evolution of these discharges, it is necessary to improve the prediction of electron thermal transport. This can be accomplished by utilizing a model for transport driven by MTMs in whole device predictive modeling codes. The objective of this paper is to develop the dispersion relation that governs the MTM driven transport. A unified fluid/kinetic approach is used in the development of a nonlinear dispersion relation for MTMs. The derivation includes the effects of electrostatic and magnetic fluctuations, arbitrary electron-ion collisionality, electron temperature and density gradients, magnetic curvature, and the effects associated with the parallel propagation vector. An iterative nonlinear approach is used to calculate the distribution function employed in obtaining the nonlinear parallel current and the nonlinear dispersion relation. The third order nonlinear effects in magnetic fluctuations are included, and the influence of third order effects on a multi-wave system is considered. An envelope equation for the nonlinear microtearing modes in the collision dominant limit is introduced in order to obtain the saturation level. In the limit that the mode amplitude does not vary along the field line, slab geometry, and strong collisionality, the fluid dispersion relation for nonlinear microtearing modes is found to agree with the kinetic dispersion relation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Discovering Active Subspaces for High-Dimensional Computer Models

Dimension reduction techniques have long been an important topic in statistics, and active subspaces (AS) have received much attention this past decade in the computer experiments literature. The most common approach towards estimating the AS is to use Monte Carlo with numerical gradient evaluation. While sensible in some settings, this approach has obvious drawbacks. Recent research has demonstrated that active subspace calculations can be obtained in closed form, conditional on a Gaussian process (GP) surrogate, which can be limiting in high-dimensional settings for computational reasons. In this paper, we produce the relevant calculations for a more general case when the model of interest is a linear combination of tensor products. These general equations can be applied to the GP, recovering previous results as a special case, or applied to the models constructed by other regression techniques including multivariate adaptive regression splines (MARS). Furthermore, using a MARS surrogate has many advantages including improved scaling, better estimation of active subspaces in high dimensions and the ability to handle a large number of prior distributions in closed form. In one real-world example, we obtain the active subspace of a radiation-transport code with 240 inputs and 9,372 model runs in under half an hour.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

ELECTRON ACCELERATION BY CASCADING RECONNECTION IN THE SOLAR CORONA. I. MAGNETIC GRADIENT AND CURVATURE DRIFT EFFECTS

We investigate the electron acceleration by magnetic gradient and curvature drift effects in cascading magnetic reconnection of a coronal current sheet via a test particle method in the framework of the guiding center approximation. After several Alfvén transit times, most of the electrons injected at the current sheet are still trapped in the magnetic islands. A small fraction of the injected electrons precipitate into the chromosphere. The acceleration of trapped electrons is dominated by the magnetic curvature drifts, which change the parallel momentum of the electron, and appears to be more efficient than the acceleration of precipitating electrons, which is dominated by the perpendicular momentum change caused by the magnetic gradient drifts. With the resulting trapped energetic electron distribution, the corresponding hard X-ray (HXR) radiation spectra are calculated using an optically thin Bremsstrahlung model. Trapped electrons may explain flare loop top HXR emission as well as the observed bright spots along current sheets trailing coronal mass ejections. The asymmetry of precipitating electrons with respect to the polarity inversion line may contribute to the observed asymmetry of footpoint emission.

79 ASTRONOMY AND ASTROPHYSICS↗