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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-D-BRC-01: An Automatic Beam Model Commissioning Method for Monte Carlo Simulations in Pencil-Beam Scanning Proton Therapy

Purpose: Monte Carlo (MC) simulation is typically regarded as the most accurate dose calculation method for proton therapy. Yet for real clinical cases, the overall accuracy also depends on that of the MC beam model. Commissioning a beam model to faithfully represent a real beam requires finely tuning a set of model parameters, which could be tedious given the large number of pencil beams to commmission. This abstract reports an automatic beam-model commissioning method for pencil-beam scanning proton therapy via an optimization approach. Methods: We modeled a real pencil beam with energy and spatial spread following Gaussian distributions. Mean energy, and energy and spatial spread are model parameters. To commission against a real beam, we first performed MC simulations to calculate dose distributions of a set of ideal (monoenergetic, zero-size) pencil beams. Dose distribution for a real pencil beam is hence linear superposition of doses for those ideal pencil beams with weights in the Gaussian form. We formulated the commissioning task as an optimization problem, such that the calculated central axis depth dose and lateral profiles at several depths match corresponding measurements. An iterative algorithm combining conjugate gradient method and parameter fitting was employed to solve the optimization problem. We validated our method in simulation studies. Results: We calculated dose distributions for three real pencil beams with nominal energies 83, 147 and 199 MeV using realistic beam parameters. These data were regarded as measurements and used for commission. After commissioning, average difference in energy and beam spread between determined values and ground truth were 4.6% and 0.2%. With the commissioned model, we recomputed dose. Mean dose differences from measurements were 0.64%, 0.20% and 0.25%. Conclusion: The developed automatic MC beam-model commissioning method for pencil-beam scanning proton therapy can determine beam model parameters with satisfactory accuracy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Greenland Surface Mass Balance as Simulated by the Community Earth System Model. Part II: Twenty-First-Century Changes

This study presents the first twenty-first-century projections of surface mass balance (SMB) changes for the Greenland Ice Sheet (GIS) with the Community Earth System Model (CESM), which includes a new ice sheet component. For glaciated surfaces, CESM includes a sophisticated calculation of energy fluxes, surface albedo, and snowpack hydrology (melt, percolation, refreezing, etc.). To efficiently resolve the high SMB gradients at the ice sheet margins and provide surface forcing at the scale needed by ice sheet models, the SMB is calculated at multiple elevations and interpolated to a finer 5-km ice sheet grid. During a twenty-first-century simulation driven by representative concentration pathway 8.5 (RCP8.5) forcing, the SMB decreases from 372 ± 100 Gt yr -1 in 1980–99 to -78 ± 143 Gt yr -1 in 2080–99. The 2080–99 near-surface temperatures over the GIS increase by 4.7 K (annual mean) with respect to 1980–99, only 1.3 times the global increase (+3.7 K). Snowfall increases by 18%, while surface melt doubles. The ablation area increases from 9% of the GIS in 1980–99 to 28% in 2080–99. Over the ablation areas, summer downward longwave radiation and turbulent fluxes increase, while incoming shortwave radiation decreases owing to increased cloud cover. The reduction in GIS-averaged July albedo from 0.78 in 1980–99 to 0.75 in 2080–99 increases the absorbed solar radiation in this month by 12%. Summer warming is strongest in the north and east of Greenland owing to reduced sea ice cover. In the ablation area, summer temperature increases are smaller due to frequent periods of surface melt.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Electronic Structures and Ferromagnetic Properties of 3d (Cr)-Doped BaSe Barium Selenide

In this study, we have employed the first-principle methods based on spin-polarized density functional theory to investigate the structural parameters, the electronic structures, and the half-metallic ferromagnetic behavior of chromium (Cr)-doped barium selenide (BaSe) such as Ba{sub 1 − x}Cr{sub x}Se at concentrations x = 0.25, 0.5, and 0.75. The exchange and correlation potential is described by the generalized gradient approximation of Wu and Cohen (GGA-WC). The calculated structural parameters of BaSe are in good agreement with theoretical data. Our findings reveal that the p-d exchange coupling is ferromagnetic for Ba{sub 0.75}Cr{sub 0.25}Se and Ba{sub 0.5}Cr{sub 0.5}S, but it becomes anti-ferromagnetic for Ba{sub 0.25}Cr{sub 0.75}S. The electronic structures exhibit that the Ba{sub 1 − x}Cr{sub x}Se materials for all concentrations are half-metallic ferromagnets with spin polarization of 100% and total magnetic moment per Cr atom of 4 μ{sub B}. Therefore, the Ba{sub 1 − x}Cr{sub x}Se compounds are suitable candidates for possible spintronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Effect of Slotted Anode on Gas Bubble Behaviors in Aluminum Reduction Cell

In the aluminum reduction cells, gas bubbles are generated at the bottom of the anode which eventually reduces the effective current contact area and the system efficiency. To encourage the removal of gas bubbles, slotted anode has been proposed and increasingly adopted by some industrial aluminum reduction cells. Nonetheless, the exact gas bubble removal mechanisms are yet to be fully understood. A three-dimensional (3D) transient, multiphase flow mathematical model coupled with magnetohydrodynamics has been developed to investigate the effect of slotted anode on the gas bubble movement. The Eulerian volume of fluid approach is applied to track the electrolyte (bath)–molten aluminum (metal) interface. Meanwhile, the Lagrangian discrete particle model is employed to handle the dynamics of gas bubbles with considerations of the buoyancy force, drag force, virtual mass force, and pressure gradient force. The gas bubble coalescence process is also taken into account based on the O’Rourke’s algorithm. The two-way coupling between discrete bubbles and fluids is achieved by the inter-phase momentum exchange. Numerical predictions are validated against the anode current variation in an industrial test. Comparing the results using slotted anode with the traditional one, the time-averaged gas bubble removal rate increases from 36 to 63 pct; confirming that the slotted anode provides more escaping ways and shortens the trajectories for gas bubbles. Furthermore, the slotted anode also reduces gas bubble’s residence time and the probability of coalescence. Moreover, the bubble layer thickness in aluminum cell with slotted anode is reduced about 3.5 mm (17.4 pct), so the resistance can be cut down for the sake of energy saving and the metal surface fluctuation amplitude is significantly reduced for the stable operation due to the slighter perturbation with smaller bubbles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Magnetism and finite-temperature effects in UZr 2 : A density functional theory analysis

The structure and the thermophysical properties of d-UZr 2 are investigated using 0 K density functional theory and ab initio molecular dynamics (AIMD). Modeling the true paramagnetic state of this intermetallic compound has been challenging using first-principles calculations. For the first time, we find that the generalized gradient approximation method without applying an on-site Coulomb interaction term (Hubbard U) can result in a ground state that is antiferromagnetic (AFM). We believe that this weak AFM ground state is the closest to the real paramagnetic state. We found that structure optimization at finite temperatures using AIMD is necessary to achieve this ground state instead of the non-magnetic and ferromagnetic states previously reported in the literature. Our findings indicate that applying the Hubbard U on uranium f-orbitals in this metallic system is unnecessary and not recommended, as it leads to a large overestimation of the volume and introduces an unphysical strong spin polarization. Finally, our approach results in atomic volume, thermal expansion, and heat capacities that have strong agreement with experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Assessment of CASL VERA for BWR analysis and application to accident tolerant SiC/SiC channel box

Application of the Virtual Environment for Reactor Applications (VERA) to BWR analysis is assessed in this paper by comparing results to those calculated using other widely-used modeling tools, namely the U.S. Nuclear Regulatory Commission’s PARCS/PATHS and the Serpent Monte Carlo particle transport code. Additionally, VERA is used to calculate 3-D temperature and fast neutron flux distributions in silicon carbide (SiC) fiber-reinforced, SiC matrix composite (SiC/SiC) BWR channel boxes, which are being studied as an Accident Tolerant Fuel core structural material concept. We found that the code-to-code comparisons were favorable, and the SiC/SiC channel box evaluation demonstrates the many advanced modeling features VERA provides while also highlighting the non-uniformity in fast neutron flux distributions that can play a role in potential SiC/SiC channel box deformation. Traditional BWR analysis tools do not have the calculation fidelity necessary for coupled assessment of flux and temperature gradients in a SiC/SiC channel box. VERA is a state-of-the-art modeling environment that was developed to increase the safety and economic competitiveness of nuclear power through improved modeling accuracy. While VERA has already been deployed in the nuclear industry for PWR applications, the current study is a vital initial step in the extensive development, validation, and verification that VERA must go through to be useful for BWR applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Hillock formation in co-deposited thin films of immiscible metal alloy systems

Hillocks are protruding sections of a surface that can influence the electronic and mechanical properties of a thin film. In monolithic films, these hillocks are formed due to thermally induced stress gradients during deposition. In this study, the surface morphology of co-sputtered immiscible Cu-X (X is a body centered cubic (BCC) group V or VI metal) thin films is characterized to elucidate the conditions that lead to hillock formation. Cu-Ta and Cu-Mo-Ag films were co-deposited with physical vapor deposition (PVD) magnetron sputtering at 25, 400, 600, and 800°C. A significant number of homogeneously distributed surface hillocks were observed only in the 600 and 800°C films. High-angle annular dark field (HAADF) cross-sectional imaging revealed significant Cu agglomerations underneath the protruding features. For the Cu-Ta films, the Cu was enveloped in a Cu-Ta nanocrystalline matrix. For Cu-Mo-Ag, the Cu was surrounded by Mo-Ag concentration modulations. While thermal stress gradients arise during deposition of immiscible metal films, biaxial stress calculations and literature cases reveal that they are not solely responsible for hillock formation. The observed morphologies align with a surface diffusion kinetic model that evaluates diffusion length of adatoms during deposition as a function of deposition temperature. Finally, this indicates a phase separation driving force paired with the constituent elements’ dissimilar mobilities at elevated deposition temperatures contribute to the presence of hillocks.

36 MATERIALS SCIENCE↗
DOE OSTI2013

The transfer between electron bulk kinetic energy and thermal energy in collisionless magnetic reconnection

By performing two-dimensional particle-in-cell simulations, we investigate the transfer between electron bulk kinetic and electron thermal energy in collisionless magnetic reconnection. In the vicinity of the X line, the electron bulk kinetic energy density is much larger than the electron thermal energy density. The evolution of the electron bulk kinetic energy is mainly determined by the work done by the electric field force and electron pressure gradient force. The work done by the electron gradient pressure force in the vicinity of the X line is changed to the electron enthalpy flux. In the magnetic island, the electron enthalpy flux is transferred to the electron thermal energy due to the compressibility of the plasma in the magnetic island. The compression of the plasma in the magnetic island is the consequence of the electromagnetic force acting on the plasma as the magnetic field lines release their tension after being reconnected. Therefore, we can observe that in the magnetic island the electron thermal energy density is much larger than the electron bulk kinetic energy density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Numerically optimized structures for dielectric asymmetric dual-grating laser accelerators

Optical scale dielectric structures are promising candidates to realize future compact, low cost particle accelerators, since they can sustain high acceleration gradients in the range of GeV/m. Here, we present numerical simulation results for a dielectric asymmetric dual-grating accelerator. It was found that the asymmetric dual-grating structures can efficiently modify the laser field to synchronize it with relativistic electrons, therefore increasing the average acceleration gradient by ∼10% in comparison to symmetric structures. The optimum pillar height which was determined by simulation agrees well with that estimated analytically. The effect of the initial kinetic energy of injected electrons on the acceleration gradient is also discussed. Finally, the required laser parameters were calculated analytically and a suitable laser is proposed as energy source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Magnetic islands and singular currents at rational surfaces in three-dimensional magnetohydrodynamic equilibria

Using the recently developed multiregion, relaxed MHD (MRxMHD) theory, which bridges the gap between Taylor's relaxation theory and ideal MHD, we provide a thorough analytical and numerical proof of the formation of singular currents at rational surfaces in non-axisymmetric ideal MHD equilibria. These include the force-free singular current density represented by a Dirac δ-function, which presumably prevents the formation of islands, and the Pfirsch-Schlüter 1/x singular current, which arises as a result of finite pressure gradient. An analytical model based on linearized MRxMHD is derived that can accurately (1) describe the formation of magnetic islands at resonant rational surfaces, (2) retrieve the ideal MHD limit where magnetic islands are shielded, and (3) compute the subsequent formation of singular currents. The analytical results are benchmarked against numerical simulations carried out with a fully nonlinear implementation of MRxMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The density gradient effect on quantum Weibel instability

The Weibel instability plays an important role in stopping the hot electrons and energy deposition mechanism in the fast ignition of inertial fusion process. In this paper, the effects of the density gradient and degeneracy on Weibel instability growth rate are investigated. Calculations show that decreasing the density degenerate in the plasma corona, near the relativistic electron beam emitting region by 8.5% leads to a 92% reduction in the degeneracy parameter and about 90% reduction in Weibel instability growth rate. Also, decreasing the degenerate density near the fuel core by 8.5% leads to 1% reduction in the degeneracy parameter and about 8.5% reduction in Weibel instability growth rate. The Weibel instability growth rate shrinks to zero and the deposition condition of relativistic electron beam energy can be shifted to the fuel core for a suitable ignition by increasing the degeneracy parameter in the first layer of plasma corona.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

FP-LAPW study of structural, electronic, elastic, mechanical and thermal properties of AlFe intermetallic

The structural, electronic, elastic, mechanical and thermal properties of AlFe intermetallic compound in B{sub 2}-type (CsCl) structure have been investigated using first-principles calculations. The exchange-correlation term was treated within generalized gradient approximation. Ground state properties i.e. lattice constants (a{sub 0}), bulk modulus (B) and first-order pressure derivative of bulk modulus (B’) are presented. The density of states are derived which show the metallic character of present compound. Our results for C{sub 11}, C{sub 12} and C{sub 44} agree well with previous theoretical data. Using Pugh’s criteria (B/G{sub H} < 1.75), brittle character of AlFe is satisfied. In addition shear modulus (G{sub H}), Young’s modulus (E), sound wave velocities and Debye temperature (θ{sub D}) have also been estimated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

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

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ≤ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased PNBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i /T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Unsupervised Learning Based Interaction Force Model for Nonspherical Particles in Incompressible Flows

This project provides a neural network-based interaction force model for gas-solid flows from low to intermediate Reynolds numbers and concentration, which can be linked to MFiX-DEM. We have constructed a database of the interaction force between the irregular-shaped particles using a spherical harmonic method and the fluid phase based on the particle-resolved direct numerical simulation (PR-DNS) with immersed boundary-based gas kinetic scheme. Unsupervised learning method, i.e., variational auto-encoder (VAE) has been applied to extract the primitive shape factors determining the drag force, lifting forces, and torque. The interaction force model has been trained and validated with a simple but effective multi-layer feed-forward neural network: multi-layer perceptron (MLP), which will be concatenated after the encoder of the previously trained VAE for geometry feature extraction for single, irregular particles. We have trained transpose convolutional neural networks with the PR-DNS data to predict the velocity and pressure gradient of the single particle systems and utilized them to calculate drag force of multi-particle systems. This model can provide high computational efficiency because it does not require collecting multiparticle system data from PR-DNS.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Towards 2+1 Flavor Lattice QCD Results for the Heavy Quark Diffusion Coefficient

We apply and extend a novel approach to non-perturbatively estimate the heavy-quark momentum diffusion coefficient κ, which is a key input for the theoretical description of heavy quarkonium production in heavy ion collisions, and is important for the understanding of the elliptic flow and nuclear suppression factor of heavy flavor hadrons. In the heavy-quark limit, this coefficient is encoded in the spectral functions of color-electric and color-magnetic correlators that we calculate on the lattice to high pre cision by applying gradient flow. In a recent study we have considered quenched QCD at 1.5 T c , where we performed a detailed study of the lat tice spacing and flow time dependence of the color-electric correlator, and, using theoretically well-established model fits for the spectral reconstruc tion, we estimated the heavy-quark diffusion coefficient. Equipped with the experience obtained in quenched QCD, we estimate $κ$ from 2+1 flavor QCD ensembles at small but finite lattice spacing and flow time without increasing systematic errors significantly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1974

Evolution of a moderately short turbulent boundary layer in a severe pressure gradient

The early and intermediate development of a highly accelerated (or decelerated) turbulent boundary layer is analyzed. For sufficiently large accelerations (or pressure gradients) and for total normal strains which are not excessive, the equation for the Reynolds shear stress simplifies to give a stress that remains approximately constant as it is convected along streamlines. The theoretical results for the evolution of the mean velocity in favourable and adverse pressure gradients agree well with experiment for the cases considered. A calculation which includes mass injection at the wall is also given.

Deissler, R. G.↗