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 487 records · Page 27

Volatile transport on Venus and implications for surface geochemistry and geology

The high vapor pressure of volatile metal halides and chalcogenides (e.g., of Cu, Zn, Sn, Pb, As, Sb, Bi) at typical Venus surface temperatures, coupled with the altitude-dependent temperature gradient of approximately 8.5 K/km, is calculated to transport volatile metal vapors to the highlands of Venus, where condensation and accumulation will occur. The predicted geochemistry of volatile metals on Venus is supported by observations of CuCl in volcanic gases at Kilauea and Nyiragongo, and large enrichments of these and other volatile elements in terrestrial volcanic aerosols. A one-dimensional finite difference vapor transport model shows the diffusive migration of a thickness of 0.01 to greater than 10 microns/yr of moderately to highly volatile phases (e.g., metal halides and chalcogenides) from the hot lowlands (740 K) to the cold highlands (660 K) on Venus. The diffusive transport of volatile phases on Venus may explain the observed low emissivity of the Venusian highlands, hazes at 6-km altitude observed by two Pioneer Venus entry probes, and the Pioneer Venus entry probe anomalies at 12.5 km.

Brackett, Robert A.↗

Multi-Dimensional Quantum Effect Simulation Using a Density-Gradient Model and Script-Level Programming Techniques

A density-gradient (DG) model is used to calculate quantum-mechanical corrections to classical carrier transport in MOS (Metal Oxide Semiconductor) inversion/accumulation layers. The model is compared to measured data and to a fully self-consistent coupled Schrodinger and Poisson equation (SCSP) solver. Good agreement is demonstrated for MOS capacitors with gate oxide as thin as 21 A. It is then applied to study carrier distribution in ultra short MOSFETs (Metal Oxide Semiconductor Field Effect Transistor) with surface roughness. This work represents the first implementation of the DG formulation on multidimensional unstructured meshes. It was enabled by a powerful scripting approach which provides an easy-to-use and flexible framework for solving the fourth-order PDEs (Partial Differential Equation) of the DG model.

Rafferty, Connor S.↗

Tectonic-Climate Interactions in Action Orogenic Belts: Quantification of Dynamic Topography with SRTM data

This project was undertaken to examine the approach to steady state in collisional mountain belts. Although the primary thrust of this grant was to look at larger collisional mountain belts, such as the Himalaya, the Tien Shan, and Southern Alps, we began by looking at smaller structures represented by growing and propagating folds. Like ranges that are evolving toward a topographic steady state, these folds undergo a series of morphologic changes as they are progressively uplifted and eroded. We wanted to document the nature of these changes and to try to discern some of the underlying controls on them. We initially focused on the Wheeler Ridge anticline in southern California. Subsequently, we progressed to looking at the topographic development and the effects of differential uplift and glaciation on the Kyrgyz Range in the northern Tien Shan. This range is unusual inasmuch as it is transformed along its length from a simple uplift with a largely preserved Mesozoic erosion surface arching across it to a highly dissected and heavily glaciated uplift in the region where uplift has been sustained at higher rates over longer intervals. In efforts to understand the distribution of erosion rates at 10(exp 3) - 10(exp 5) year time scales, cosmogenic radionuclide (CRN) concentrations have been gaining increasingly widespread usage (Brown et al., 1995; Riebe et al., 2004; Riebe et al., 2001; Vance et al., 2003). Most studies to date, however, have been conducted in slowly eroding ranges. In rapidly eroding mountains where landslides deliver most of the sediments to the rivers, we hypothesized that CRN concentrations could be highly perturbed by the stochastic processes of landsliding. Therefore, we undertook the development of a numerical model that simulated the effects of both landsliding and grain-by-grain attrition within fluvial catchments. This modeling effort has shown the effects of catchment size and erosion rate on CRN concentrations and allows a prediction of where to sample to obtain the optimal erosion rate estimates using CRN techniques. Finally, we developed computational techniques to operate on DEMs to extract useful information that would enable quantification of climate-erosion interactions. In particular, we worked on rapid techniques to define catchments of any given range of sizes, to extract channel gradients, to combine precipitation information to calculate discharge, and to utilize various stream-power models to determine the erosional energy within any given catchment within a transect. We briefly describe results from Wheeler Ridge, the Kyrgyz Range, the Nepal Himalaya, and our numerical modeling.

Burbank, Douglas W.↗

Anomalous elastic softening in ferroelectric hafnia under pressure

his study employs first-principles density-functional theory (DFT) calculations to explore the elastic and mechanical properties of ferroelectric hafnia (HfO 2 ) in its polar orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase under varying hydrostatic pressure conditions up to 30 GPa. Utilizing a plane-wave basis set and Perdew-Burke-Ernzerhof generalized-gradient approximation for solids in our DFT calculations, we investigate both pure and yttrium-substituted HfO 2 . Our findings reveal an anomalous reduction in the 𝐶 33 component of the elastic tensor with increasing pressure, which becomes significant above 15 GPa and signals a potential pressure-driven structural phase transition at higher pressure. The analysis of atomic displacements under pressure sheds light on the unusual mechanical behavior and phase stability of this material. Additionally, we observe a transition from an indirect band gap to a direct band gap with increasing pressure, which could have significant implications for optical applications. Here, the effects of yttrium substitution on the mechanical and electronic properties are further examined, revealing that yttrium substitution softens the elastic response of this material and reduces the electronic band gap. These results enhance our understanding of elastic and mechanical responses of ferroelectric hafnia and its potential for applications in microelectronics, piezoelectric devices, and nonvolatile ferroelectric random-access memories. Further experimental validation is recommended to confirm our predictions and explore the practical implications of the observed phase transitions and electronic behavior of the ferroelectric hafnia under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparative energetics of carbon storage molecules in green algae

Several members of the green algae possess the ability to produce lipids and/or high value compounds in significant quantities. While for several of these green algal species induction of increased lipid production has been shown, and cultivation of species for high value molecules occurs at production scale, the molecular mechanisms governing over-accumulation of molecules synthesized from isoprenoid precursors, carotenoids, for example, have received far less attention. Here, we present a calculation of the required ATP equivalencies per carbon atom and reducing power equivalencies as NADH/NADPH (NAD(P)H) per carbon atom for the isoprenoid molecules ..beta..-carotene (C40), astaxanthin (C40), and squalene (C30). We compared energetic requirements of carbohydrates, triacylglycerol, and isoprenoid molecules under a gradient of conditions of cellular stress. Our calculations revealed slightly less ATP and NAD(P)H equivalency per carbon atom between triacylglycerol and the three isoprenoid molecules. Based on our results, we propose that the driving force for differences in accumulation patterns of carotenoids vs. triacylglycerols in algal cells under stress is largely dependent on the presence and regulation of bypass mechanisms at metabolic junction bottlenecks, like pyruvate dehydrogenase (PDH), within particular species. We provide a discussion of several molecular mechanisms that may influence carbon partitioning within different groups of green algae, including metabolic inhibition through accumulation of specific substrates related to ATP and reducing equivalent production (NAD(P)H) as well as cellular compartmentalization. This work contributes to the ongoing discussion of cellular homeostatic regulation during stress, as well as the potential mechanisms driving long-term carbon storage as it relates to energy and redox states within the algal cell.

09 BIOMASS FUELS↗

Molecular dynamics simulation of methane transport in confined organic nanopores with high relative roughness

Understanding and characterizing the transport of shale gas (methane) through the nanopores of kerogens are critical for the accurate prediction of shale gas recovery. However, the key factors that regulate shale gas transport through highly roughened nanopores of shale kerogens are not fully understood. Here in this work, methane transport in organic nanopores with a high relative roughness is characterized using equilibrium and non-equilibrium molecular dynamics methods. According to our results, the CH 4 mass flux has a linear relationship with the pressure gradient, consistent with previous studies, while the calculated slip lengths and gas fluxes varied with different roughness geometries in the order of sigmoidal ≥ triangular > rectangular. Surface slip flow can be a major contributor to the overall gas flux, but surprisingly, the relative contribution of surface slip flow is independent of the pressure gradient. In contrast, the contributions of both slip flow and the average gas fluxes vary strongly with pore diameters. Typical contributions of the adsorbed layer to the overall gas flux are in the range 20–40% but vary from as high as 74% in a 4-nm pore to as low as 6% in a 16-nm pore. Compared to smooth nanopores, we find that, in nanopores with realistically high degrees of relative roughness, methane confinement in cavities decouples slip flow from the flow in the pore interior, significantly reducing the overall flux.

03 NATURAL GAS↗

Remote Sensing techniques used to characterize soil erosion in southwestern Sao Paulo state

Within randomly sampled squares of a 1 km x 1 km grid, rill/gullies frequency, land cover/land use type and shape of the slopes were extracted from aerial photographs of the Ribeirao Anhumas drainage basin. Mean slope gradient, stream frequency and slope length were calculated on topographic maps. Ground truth data on fine sand/coarse sand ratio and vegetation cover densities were obtained. The MSS-LANDSAT-2 data (CCTs) were analyzed using single-cell, cluster synthesis and slicer algorithms. Graphical and statistical analyses of the data indicate that different slope gradients and land cover/land use types are the most significant factors related to the soil erosion process. The digital analysis of MSS data allowed the association among gray level classes and vegetation cover classes, which defined seven classes. These gray level classes and slope gradient classes were used to rank erosion risk.

Parada, N. D. J.↗

Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities

Here, we present a framework for calibration of parameters in elastoplastic constitutive models that is based on the use of automatic differentiation (AD). The model calibration problem is posed as a partial differential equation-constrained optimization problem where a finite element (FE) model of the coupled equilibrium equation and constitutive model evolution equations serves as the constraint. The objective function quantifies the mismatch between the displacement predicted by the FE model and full-field digital image correlation data, and the optimization problem is solved using gradient-based optimization algorithms. Forward and adjoint sensitivities are used to compute the gradient at considerably less cost than its calculation from finite difference approximations. Through the use of AD, we need only to write the constraints in terms of AD objects, where all of the derivatives required for the forward and inverse problems are obtained by appropriately seeding and evaluating these quantities. We present three numerical examples that verify the correctness of the gradient, demonstrate the AD approach's parallel computation capabilities via application to a large-scale FE model, and highlight the formulation's ease of extensibility to other classes of constitutive models.

42 ENGINEERING↗

Predicting the Rate of Homogeneous Intermetallic Nucleation within Steep Composition Gradients

Simulations of isothermal homogeneous nucleation from deeply undercooled amorphous melts exhibit systematic variations in nucleation behavior depending upon the strength of an imposed composition gradient. Data from molecular dynamics (MD) simulations in a model Ni/Al system permit quantification of the nucleation rate of the NiAl-B2 intermetallic phase and indicate that nucleation proceeds in a polymorphous mode. The nucleation rate decreases with increasing gradient, and nucleation is completely suppressed above a critical gradient. Based on an argument of crystal nucleus stability, a simple estimate provides good prediction of the critical gradient. Here, a modified classical nucleation model parameterized with thermodynamic and kinetic quantities calculated independently predicts the nucleation rates and their variation with the imposed gradient, matching the MD results very well.

36 MATERIALS SCIENCE↗

Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method

We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Buoyant Shear Flows Using Anisotropic Dissipation Rate Models

This paper examines the modeling of two-dimensional homogeneous stratified turbulent shear flows using the Reynolds-stress and Reynolds-heat-flux equations. Several closure models have been investigated-, the emphasis is placed on assessing the effect of modeling the dissipation rate tensor in the Reynolds-stress equation. Three different approaches are considered: one is an isotropic approach while the other two are anisotropic approaches. The isotropic approach is based on Kolmogorov's hypothesis and a dissipation rate equation modified to account for vortex stretching. One of the anisotropic approaches is based on an algebraic representation of the dissipation rate tensor, while another relies on solving a modeled transport equation for this tensor. In addition, within the former anisotropic approach, two different algebraic representations are examined one is a function of the Reynolds-stress anisotropy tensor, and the other is a function of' the mean velocity gradients. The performance of these closure models is evaluated against experimental and direct numerical simulation data of pure shear flows. pure buoyant flows and buoyant shear flows. Calculations have been carried out over a range of Richardson numbers (Ri) and two different Prandtl numbers (Pr); thus the effect of Pr on the development of counter-gradient heat flux in a stratified shear flow can be assessed. At low Ri, the isotropic model performs well in the predictions of stratified shear flows; however, its performance deteriorates as Ri increases. At high Ri, the transport equation model for the dissipation rate tensor gives the best result. Furthermore, the results also lend credence to the algebraic dissipation rate model based on the Reynolds stress anisotropy tensor. Finally, it is found that Pr has an effect on the development of counter-gradient heat flux. The calculations show that, under the action of shear, counter-gradient heat flux does not occur even at Ri = 1 in an air flow.

So, R. M. C.↗

Calculation of Sensitivity Derivatives in an MDAO Framework

During gradient-based optimization of a system, it is necessary to generate the derivatives of each objective and constraint with respect to each design parameter. If the system is multidisciplinary, it may consist of a set of smaller "components" with some arbitrary data interconnection and process work ow. Analytical derivatives in these components can be used to improve the speed and accuracy of the derivative calculation over a purely numerical calculation; however, a multidisciplinary system may include both components for which derivatives are available and components for which they are not. Three methods to calculate the sensitivity of a mixed multidisciplinary system are presented: the finite difference method, where the derivatives are calculated numerically; the chain rule method, where the derivatives are successively cascaded along the system's network graph; and the analytic method, where the derivatives come from the solution of a linear system of equations. Some improvements to these methods, to accommodate mixed multidisciplinary systems, are also presented; in particular, a new method is introduced to allow existing derivatives to be used inside of finite difference. All three methods are implemented and demonstrated in the open-source MDAO framework OpenMDAO. It was found that there are advantages to each of them depending on the system being solved.

Moore, Kenneth T.↗

SU-C-BRA-02: Gradient Based Method of Target Delineation On PET/MR Image of Head and Neck Cancer Patients

Purpose: Validate the consistency of a gradient-based segmentation tool to facilitate accurate delineation of PET/CT-based GTVs in head and neck cancers by comparing against hybrid PET/MR-derived GTV contours. Materials and Methods: A total of 18 head and neck target volumes (10 primary and 8 nodal) were retrospectively contoured using a gradient-based segmentation tool by two observers. Each observer independently contoured each target five times. Inter-observer variability was evaluated via absolute percent differences. Intra-observer variability was examined by percentage uncertainty. All target volumes were also contoured using the SUV percent threshold method. The thresholds were explored case by case so its derived volume matched with the gradient-based volume. Dice similarity coefficients (DSC) were calculated to determine overlap of PET/CT GTVs and PET/MR GTVs. Results: The Levene’s test showed there was no statistically significant difference of the variances between the observer’s gradient-derived contours. However, the absolute difference between the observer’s volumes was 10.83%, with a range from 0.39% up to 42.89%. PET-avid regions with qualitatively non-uniform shapes and intensity levels had a higher absolute percent difference near 25%, while regions with uniform shapes and intensity levels had an absolute percent difference of 2% between observers. The average percentage uncertainty between observers was 4.83% and 7%. As the volume of the gradient-derived contours increased, the SUV threshold percent needed to match the volume decreased. Dice coefficients showed good agreement of the PET/CT and PET/MR GTVs with an average DSC value across all volumes at 0.69. Conclusion: Gradient-based segmentation of PET volume showed good consistency in general but can vary considerably for non-uniform target shapes and intensity levels. PET/CT-derived GTV contours stemming from the gradient-based tool show good agreement with the anatomically and metabolically more accurate PET/MR-derived GTV contours, but tumor delineation accuracy can be further improved with the use PET/MR.

60 APPLIED LIFE SCIENCES↗

Zodiacal emission. II - Dust near ecliptic

The location and brightness of the peak zodiacal emission are derived from IRAS observations of the diffuse infrared background radiation. A uniform ellipsoid geometrical model is fitted to the latitudes of peak brightness, from which the orientation of the surface of maximum dust density near the earth's orbit is derived. Deviations from the uniform ellipsoid indicate that this surface bends toward the invariable plane of the solar system beyond 2.5 AU from the sun. The variations of the peak zodiacal emission brightness with solar elongation and observation date allow two independent determinations of the gradient of the volume emissivity of interplanetary dust. Combined with theoretical calculations of the infrared emission from silicate and graphite grains heated by the sun, the gradient of the dust density is found to be n proportional to r exp (-1.1 + or - 0.2). Systematic brightness variations with ecliptic longitude reveal two potential asymmetries in the dust complex: excess emission following the earth in its orbit, and excess emission along lines of sight with ecliptic longitude lambda = 190 + or - 20 deg and solar lambda about 100 deg.

Reach, William T.↗

Secondary subharmonic instability of boundary layers with pressure gradient and suction

Three-dimensional linear secondary instability is investigated for boundary layers with pressure gradient and suction in the presence of a finite amplitude TS wave. The focus is on principal parametric resonance responsible for a strong growth of subharmonics in a low disturbance environment. Calculations are presented for the effect of pressure gradients and suction on controlling the onset and amplification of the secondary instability.

El-Hady, Nabil M.↗

Hybrid functional pseudopotentials

The consistency between the exchange-correlation functional used in pseudopotential construction and in the actual density functional theory calculation is essential for the accurate prediction of fundamental properties of materials. However, routine hybrid density functional calculations at present still rely on generalized gradient approximation pseudopotentials due to the lack of hybrid functional pseudopotentials. Here, we present a scheme for generating hybrid functional pseudopotentials, and we analyze the importance of pseudopotential density functional consistency for hybrid functionals. For the PBE0 hybrid functional, we benchmark our pseudopotentials for structural parameters and fundamental electronic gaps of the Gaussian-2 (G2) molecular dataset and some simple solids. Our results demonstrate that using our PBE0 pseudopotentials in PBE0 calculations improves agreement with respect to all-electron calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic structure of ferromagnetic semiconductor material on the monoclinic and rhombohedral ordered double perovskites La{sub 2}FeCoO{sub 6}

Double perovskite La{sub 2}FeCoO{sub 6} with monoclinic structure and rhombohedra structure show as ferromagnetic semiconductor based on density functional theory calculation. The ferromagnetic semiconductor state can be well explained by the superexchange interaction. Moreover, the ferromagnetic semiconductor state remains under the generalized gradient approximation (GGA) and GGA plus onsite Coulomb interaction calculation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calculated diffusion coefficients and the growth rate of olivine in a basalt magma

Concentration gradients in glass adjacent to skeletal olivines in a basalt have been examined by electron probe. The glass is depleted in Mg, Fe, and Cr and enriched in Si, Al, Na, and Ca relative to that far from olivine. Ionic diffusion coefficients for the glass compositions are calculated from temperature, ionic radius and melt viscosity, using the Stokes-Einstein relation. At 1170 C, the diffusion coefficient of Mg(2+) ions in the basalt is 4.5 billionths sq cm per sec. Comparison with measured diffusion coefficients in a mugearite suggests this value may be 16 times too small. The concentration gradient data and the diffusion coefficients are used to calculate instantaneous olivine growth rates. Growth necessarily preceded emplacement such that the composition of the crystals plus the enclosing glass need not be that of a melt. The computed olivine growth rates are compatible with the rate of crystallization deduced for the Skaegaard intrusion.

Donaldson, C. H.↗