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,675 records · Page 93

MPS solidification model. Analysis and calculation of macrosegregation in a casting ingot

Work performed on several existing solidification models for which computer codes and documentation were developed is presented. The models describe the solidification of alloys in which there is a time varying zone of coexisting solid and liquid phases; i.e., the S/L zone. The primary purpose of the models is to calculate macrosegregation in a casting or ingot which results from flow of interdendritic liquid in this S/L zone during solidification. The flow, driven by solidification contractions and by gravity acting on density gradients in the interdendritic liquid, is modeled as flow through a porous medium. In Model 1, the steady state model, the heat flow characteristics are those of steady state solidification; i.e., the S/L zone is of constant width and it moves at a constant velocity relative to the mold. In Model 2, the unsteady state model, the width and rate of movement of the S/L zone are allowed to vary with time as it moves through the ingot. Each of these models exists in two versions. Models 1 and 2 are applicable to binary alloys; models 1M and 2M are applicable to multicomponent alloys.

Poirier, D. R.↗

Horizontal mixing coefficients for two-dimensional chemical models calculated from National Meteorological Center Data

Calculations of the two-dimensional, species-independent mixing coefficients for two-dimensional chemical models for the troposphere and stratosphere are performed using quasi-geostrophic potential vorticity fluxes and gradients from 4 years of National Meteorological Center data for the four seasons in both hemispheres. Results show that the horizontal mixing coefficient values for the winter lower stratosphere are broadly consistent with those currently employed in two-dimensional models, but the horizontal mixing coefficient values in the northern winter upper stratosphere are much larger than those usually used.

Newman, P. A.↗

A Modified Mixing Length Turbulence Model for Zero and Adverse Pressure Gradients

The modified mixing length (MML) turbulence model was installed in the Proteus Navier-Stokes code, then modified to make it applicable to a wider range of flows typical of aerospace propulsion applications. The modifications are based on experimental data for three flat-plate flows having zero, mild adverse, and strong adverse pressure gradients. Three transonic diffuser test cases were run with the new version of the model in order to evaluate its performance. All results are compared with experimental data and show improvements over calculations made using the Baldwin-Lomax turbulence model, the standard algebraic model in Proteus.

Conley, Julianne M.↗

The Favre-Reynolds Average Distinction and a Consistent Gradient Transport Expression for the Dissipation

Two equation and higher order closures for compressible turbulence fail to capture the compressible wall layers' log scaling. Accounting for the distinction between Favre and Reynolds averaged variables in the compressible moment equations indicate that turbulent transport expressions obtained using the 'variable density approximation' are in error. The error is related to the enstrophy, a Reynolds averaged variable appearing in the equation for the Favre averaged k; recognizing this fact an expression for the transport of dissipation consistent with simple mixing length arguments is obtained. Within the (limited) context of a gradient transport hypothesis a rational form for the turbulent transport of the dissipation is found. Modestly better agreement with the well established compressible Van Driest log scaling is found in k - epsilon calculation.

Ristorcelli, J. R.↗

A Grid-resolved Analysis of Base Flowfield for a Four-Engine Clustered Nozzle Configuration

The objective of this study is to propose a computational methodology that can effectively anchor the base flowfield of a four-engine clustered nozzle configuration. This computational methodology is based on a three-dimensional, viscous flow, pressure-based computational fluid dynamics (CFD) formulation. For efficient CFD calculation, a Prandtl-Meyer solution treatment is applied to the algebraic grid lines for initial plume expansion resolution. As the solution evolves, the computational grid is adapted to the pertinent flow gradients. The CFD model employs an upwind scheme in which second- and fourth-order central differencing schemes with artificial dissipation are used. The computed quantitative base flow properties such as the radial base pressure distributions, model centerline static pressure, Mach number and impact pressure variations, and base pressure characteristic curve agreed reasonably well with those of the measurement.

Wang, Ten-See↗

GPS World, Innovation: Autonomous Navigation at High Earth Orbits

Calculating a spacecraft's precise location at high orbital altitudes-22,000 miles (35,800 km) and beyond-is an important and challenging problem. New and exciting opportunities become possible if satellites are able to autonomously determine their own orbits. First, the repetitive task of periodically collecting range measurements from terrestrial antennas to high altitude spacecraft becomes less important-this lessens competition for control facilities and saves money by reducing operational costs. Also, autonomous navigation at high orbital altitudes introduces the possibility of autonomous station keeping. For example, if a geostationary satellite begins to drift outside of its designated slot it can make orbit adjustments without requiring commands from the ground. Finally, precise onboard orbit determination opens the door to satellites flying in formation-an emerging concept for many scientific space applications. The realization of these benefits is not a trivial task. While the navigation signals broadcast by GPS satellites are well suited for orbit and attitude determination at lower altitudes, acquiring and using these signals at geostationary (GEO) and highly elliptical orbits is much more difficult. The light blue trace describes the GPS orbit at approximately 12,550 miles (20,200 km) altitude. GPS satellites were designed to provide navigation signals to terrestrial users-consequently the antenna array points directly toward the earth. GEO and HE0 orbits, however, are well above the operational GPS constellation, making signal reception at these altitudes more challenging. The nominal beamwidth of a Block II/IIA GPS satellite antenna array is approximately 42.6 degrees. At GEO and HE0 altitudes, most of these primary beam transmissions are blocked by the Earth, leaving only a narrow region of nominal signal visibility near opposing limbs of the earth. This region is highlighted in gray. If GPS receivers at GEO and HE0 orbits were designed to use these higher power signals only, precise orbit determination would not be practical. Fortunately, the GPS satellite antenna array also produces side lobe signals at much lower power levels. NASA has designed and tested the Navigator, a new GPS receiver that can acquire and track these weaker signals, thereby dramatically increasing the signal visibility at these altitudes. While using much weaker signals is a fundamental requirement for a high orbital altitude GPS receiver, it is certainly not the only challenge. There are other unique characteristics of this application that must also be considered. For example, Position Dilution of Precision (PDOP) figures are much higher at GEO and HE0 altitudes because visible GPS satellites are concentrated in a much smaller area with respect to the spacecraft antenna. These poor PDOP values contribute considerable error to the point solutions calculated by the spacecraft GPS receiver. Finally, spacecraft GPS receivers must be designed to withstand a variety of extreme environmental conditions. Variations in acceleration between launch and booster separation are extreme. Temperature gradients in the space environment are also severe. Furthermore, radiation effects are a major concern-spacecraft-borne GPS receivers must be designed with radiation-hardened electronics to guard against this phenomenon, otherwise they simply will not work. Perhaps most importantly, there are no opportunities to repair or modify any space-borne GPS receiver after it has been launched. Great care must be taken to ensure all performance characteristics have been analyzed prior to liftoff.

Bamford, William↗

Simulation Studies for a Space-Based CO2 Lidar Mission

We report results of initial space mission simulation studies for a laser-based, atmospheric CO2 sounder, which are based on real-time carbon cycle process modelling and data analysis. The mission concept corresponds to the Active Sensing of CO2 Emissions over Nights, Days and Seasons (ASCENDS) recommended by the US National Academy of Sciences' Decadal Survey. As a pre-requisite for meaningful quantitative evaluation, we employ a CO2 model that has representative spatial and temporal gradients across a wide range of scales. In addition, a relatively complete description of the atmospheric and surface state is obtained from meteorological data assimilation and satellite measurements. We use radiative transfer calculations, an instrument model with representative errors and a simple retrieval approach to quantify errors in 'measured' CO2 distributions, which are a function of mission and instrument design specifications along with the atmospheric/surface state. Uncertainty estimates based on the current instrument design point indicate that a CO2 laser sounder can provide data consistent with ASCENDS requirements and will significantly enhance our ability to address carbon cycle science questions. Test of a dawn/dusk orbit deployment, however, shows that diurnal differences in CO2 column abundance, indicative of plant photosynthesis and respiration fluxes, will be difficult to detect

Kawa, S. R.↗

WE-EF-210-06: Ultrasound 2D Strain Measurement of Radiation-Induced Toxicity: Phantom and Ex Vivo Experiments

Purpose: Radiation-induced fibrosis is a common long-term complication affecting many patients following cancer radiotherapy. Standard clinical assessment of subcutaneous fibrosis is subjective and often limited to visual inspection and palpation. Ultrasound strain imaging describes the compressibility (elasticity) of biological tissues. This study’s purpose is to develop a quantitative ultrasound strain imaging that can consistently and accurately characterize radiation-induce fibrosis. Methods: In this study, we propose a 2D strain imaging method based on deformable image registration. A combined affine and B-spline transformation model is used to calculate the displacement of tissue between pre-stress and post-stress B-mode image sequences. The 2D displacement is estimated through a hybrid image similarity measure metric, which is a combination of the normalized mutual information (NMI) and normalized sum-of-squared-differences (NSSD). And 2D strain is obtained from the gradient of the local displacement. We conducted phantom experiments under various compressions and compared the performance of our proposed method with the standard cross-correlation (CC)- based method using the signal-to-noise (SNR) and contrast-to-noise (CNS) ratios. In addition, we conducted ex-vivo beef muscle experiment to further validate the proposed method. Results: For phantom study, the SNR and CNS values of the proposed method were significantly higher than those calculated from the CC-based method under different strains. The SNR and CNR increased by a factor of 1.9 and 2.7 comparing to the CC-based method. For the ex-vivo experiment, the CC-based method failed to work due to large deformation (6.7%), while our proposed method could accurately detect the stiffness change. Conclusion: We have developed a 2D strain imaging technique based on the deformable image registration, validated its accuracy and feasibility with phantom and ex-vivo data. This 2D ultrasound strain imaging technology may be valuable as physicians try to eliminate radiation-induce fibrosis and improve the therapeutic ratio of cancer radiotherapy. This research is supported in part by DOD PCRP Award W81XWH-13-1-0269, and National Cancer Institute (NCI) Grant CA114313.

60 APPLIED LIFE SCIENCES↗

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal hybrid exchange-correlation density functional for improving the description of warm dense matter

Finite-temperature density functional theory (FT-DFT) has become a standard tool for first principles calculations of the properties of warm dense matter (WDM) relevant to high-energy density physics (HEDP) applications. In this work we introduce theoretical grounds of thermal hybrid exchange-correlation (XC) functionals within the generalized Mermin-Kohn-Sham scheme for an improved description of WDM. Building on the previously developed KSDT (Karasiev-Sjostrom-Dufty-Trickey) [Karasiev et al., Phys. Rev. Lett. 112, 076403 (2014)] local density approximation (LDA) and the KDT16 (Karasiev-Dufty-Trickey 2016) [Karasiev et al., Phys. Rev. Lett. 120, 076402 (2018)] generalized-gradient approximation (GGA) XC free-energy density functionals, we construct a new thermal hybrid XC functional, referred to here as KDT0. The KDT0 model at low temperature reduces to the popular ground-state PBE0 hybrid due to properties of the used KDT16 density functional approximation. Application to static calculations of electronic band gap and band structure at a wide range of temperatures for various systems of interest to HEDP show that KDT0 provides a significant improvement to the lower LDA and GGA rung XC functionals and to the ground-state PBE0 hybrid.

36 MATERIALS SCIENCE↗

Hyperfine interactions for small systems including transition-metal elements using self-interaction corrected density-functional theory

The interactions between the electronic magnetic moment and the nuclear spin moment, i.e., magnetic hyperfine (HF) interactions, play an important role in understanding electronic properties of magnetic systems and in realizing platforms for quantum information science applications. We investigate the HF interactions for atomic systems and small molecules, including Ti or Mn, by using Fermi–Löwdin orbital (FLO) based self-interaction corrected (SIC) density-functional theory. We calculate the Fermi contact (FC) and spin-dipole terms for the systems within the local density approximation (LDA) in the FLO-SIC method and compare them with the corresponding values without SIC within the LDA and generalized-gradient approximation (GGA), as well as experimental data. For the moderately heavy atomic systems (atomic number Z ≤ 25), we find that the mean absolute error of the FLO-SIC FC term is about 27 MHz (percentage error is 6.4%), while that of the LDA and GGA results is almost double that. Therefore, in this case, the FLO-SIC results are in better agreement with the experimental data. For the non-transition-metal molecules, the FLO-SIC FC term has the mean absolute error of 68 MHz, which is comparable to both the LDA and GGA results without SIC. For the seven transition-metal-based molecules, the FLO-SIC mean absolute error is 59 MHz, whereas the corresponding LDA and GGA errors are 101 and 82 MHz, respectively. Therefore, for the transition-metal-based molecules, the FLO-SIC FC term agrees better with experiment than the LDA and GGA results. We observe that the FC term from the FLO-SIC calculation is not necessarily larger than that from the LDA or GGA for all the considered systems due to the core spin polarization, in contrast to the expectation that SIC would increase the spin density near atomic nuclei, leading to larger FC terms.

Chemistry↗

Additively manufactured Ni-20Cr to V functionally graded material: Computational predictions and experimental verification of phase formations

Throughout this work, a database for the Cr-Ni-V system was constructed by modeling the binary Cr-V and ternary Cr-Ni-V systems using the CALPHAD approach aided by density functional theory (DFT)-based first-principles calculations and ab initio molecular dynamics (AIMD) simulations. To validate this new database, a functionally graded material (FGM) using Ni-20Cr and elemental V was fabricated using directed energy deposition additive manufacturing (DED AM) and experimentally characterized. The deposited Ni-20Cr was pure fcc phase, while increasing the amount of V across the gradient resulted in the formation of sigma phase, followed by the bcc phase. The experimentally measured phase data was compared with computational predictions made using a Cr-Ni-V thermodynamic database from the literature as well as the database developed in the present work. The newly developed database was shown to better predict the experimentally observed phases due to its accurate modeling of binary systems within the database and the ternary liquid phase, which is critical for accurate Scheil calculations.

36 MATERIALS SCIENCE↗

Radiation-induced bowing of SiC/SiC composites under neutron flux gradients—integral experimental data for model validation

Here, the radiation-induced swelling of SiC and its composites, including strong dependencies on temperature and dose, can drive significant lateral bowing in the presence of temperature and/or dose gradients. In recent years, simulations have been performed to assess the extent of bowing in SiC composite light-water reactor (LWR) fuel cladding and boiling water reactor (BWR) channel boxes. However, to date, no integral experimental data exist to validate these models. This work provides the first experimental bowing evaluation of three ∼380 mm long SiC composite specimens irradiated under varying neutron dose gradients (∼50°C–60°C, 0.03–0.06 dpa): two tubes (∼9.8 mm diameter) and a miniature BWR channel box (∼30 mm square). The measured radiation-induced length swelling (∼0.3%–0.7% linear) was consistently 10%–21% higher than values obtained from 3D finite element structural analyses with inputs from 3D radiation transport calculations. This discrepancy could be at least partially explained by differences in dose rate (∼10 -8 dpa/s) compared to the literature data (∼10-6 dpa/s) used to establish the dose-to-swelling correlations in the model. Nevertheless, the modeled bowing magnitudes (<2 mm) obtained from finite element analyses and simple analytical equations were within the bounds of the experimental measurements for all specimens. With improved confidence in the ability to predict the structural response and measure the macroscopic deformations, future experiments will target transient bowing under neutron flux gradients at representative LWR temperatures and assess whether grid spacers can mitigate the tens of millimeters of bowing that would otherwise be expected in ∼4 m long LWR components.

bowing↗

Shock formation induced by poloidal flow and its effects on the edge stability in tokamaks

In the high confinement mode of tokamaks, magnitude of the radial electric field increases at the edge. Thus, the poloidal flow inside the transport barrier can be sonic when the edge pressure gradient is not steep enough to make the poloidal flow subsonic. When the poloidal Mach number is close to unity, a shock appears in the low field side and causes a large density perturbation. In this study, we describe a shock induced by the sonic poloidal plasma flow. Then, an entropy production across the shock is calculated. Finally, we introduce a simple model for Type III edge localized modes using the poloidal density variation driven by the sonic poloidal flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fractional occupation numbers and self-interaction correction-scaling methods with the Fermi-Löwdin orbital self-interaction correction approach

In this work, we present a new assessment of the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) approach with an emphasis on its performance for predicting energies as a function of fractional occupation numbers (FONs) for various multielectron systems. Our approach is implemented in the massively parallelized NWChem quantum chemistry software package and has been benchmarked on the prediction of total energies, atomization energies, and ionization potentials of small molecules and relatively large aromatic systems. Within our study, we also derive an alternate expression for the FLO-SIC energy gradient expressed in terms of gradients of the Fermi-orbital eigenvalues and revisit how the FLO-SIC methodology can be seen as a constrained unitary transformation of the canonical Kohn–Sham orbitals. Finally, we conclude with calculations of energies as a function of FONs using various SIC-scaling methods to test the limits of the FLO-SIC formalism on a variety of multielectron systems. We find that these relatively simple scaling methods do improve the prediction of total energies of atomic systems as well as enhance the accuracy of energies as a function of FONs for other multielectron chemical species.

fractional occupation number↗

Statistical Shape Modeling of Cam Femoroacetabular Impingement

In this study, statistical shape modeling (SSM) was used to quantify three-dimensional (3D) variation and morphologic differences between femurs with and without cam femoroacetabular impingement (FAI). 3D surfaces were generated from CT scans of femurs from 41 controls and 30 cam FAI patients. SSM correspondence particles were optimally positioned on each surface using a gradient descent energy function. Mean shapes for control and patient groups were defined from the resulting particle configurations. Morphological differences between group mean shapes and between the control mean and individual patients were calculated. Principal component analysis was used to describe anatomical variation present in both groups. The first 6 modes (or principal components) captured statistically significant shape variations, which comprised 84% of cumulative variation among the femurs. Shape variation was greatest in femoral offset, greater trochanter height, and the head-neck junction. The mean cam femur shape protruded above the control mean by a maximum of 3.3 mm with sustained protrusions of 2.5-3.0 mm along the anterolateral head-neck junction and distally along the anterior neck, corresponding well with reported cam lesion locations and soft-tissue damage. This study provides initial evidence that SSM can describe variations in femoral morphology in both controls and cam FAI patients and may be useful for developing new measurements of pathological anatomy. SSM may also be applied to characterize cam FAI severity and provide templates to guide patient-specific surgical resection of bone.

statistical shape modeling↗

Meso-scale anisotropic hydrogen segregation near grain-boundaries in polycrystalline nickel characterized by EBSD/SIMS

To study anisotropic hydrogen segregation and diffusion in nickel polycrystalline, Secondary Ion Mass Spectrometry (SIMS) and Electron Back Scattered Diffraction (EBSD) are integrated to investigate hydrogen distribution around grain boundaries. Hydrogen distribution in pre-charged samples were correlated with grain boundary character by integrating high-resolution grain microstructure from EBSD inverse pole figure map and low-resolution hydrogen concentration profile map from SIMS. This multimodal imaging instrumentation shows that grain boundaries in nickel can be categorized into two families based on behavior of hydrogen distribution crossing grain boundary: the first one includes random grain boundaries with fast hydrogen diffusivity, showing a sharp gap for hydrogen concentration profile cross the grain boundaries. The second family are special Σ3n grain boundaries with low hydrogen diffusivity, showing a smooth gradient of hydrogen concentration cross the grain boundary. Heterogeneous hydrogen distributions due to grain boundary family revealed by SIMS/EBSD on mesoscale further validate the recent hydrogen permeation data and anisotropic ab-initio calculations in nanoscale. The results highlight the fact that grain boundaries character impacts hydrogen distribution significantly.

Electron Back Scattered Diffraction↗

Topological Dirac Semimetal Phase in Bismuth Based Anode Materials for Sodium-Ion Batteries

Bismuth has recently attracted interest in connection with Na-ion battery anodes due to its high volumetric capacity. It reacts with Na to form Na 3 Bi which is a prototypical Dirac semimetal with a nontrivial electronic structure. Density-functional-theory based first-principles calculations are playing a key role in understanding the fascinating electronic structure of Na 3 Bi and other topological materials. In particular, the strongly-constrained-and-appropriately-normed (SCAN) meta-generalized-gradient-approximation (meta-GGA) has shown significant improvement over the widely used generalized-gradient-approximation (GGA) scheme in capturing energetic, structural, and electronic properties of many classes of materials. Here, we discuss the electronic structure of Na 3 Bi within the SCAN framework and show that the resulting Fermi velocities and s-band shift around the Γ point are in better agreement with experiments than the corresponding GGA predictions. SCAN yields a purely spin-orbit-coupling (SOC) driven Dirac semimetal state in Na 3 Bi in contrast with the earlier GGA results. Our analysis reveals the presence of a topological phase transition from the Dirac semimetal to a trivial band insulator phase in Na 3 Bi x Sb 1-x alloys as the strength of the SOC varies with Sb content, and gives insight into the role of the SOC in modulating conduction properties of Na 3 Bi.

36 MATERIALS SCIENCE↗