SEARCH · Search NASA

Results for “gradient calculation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,621 records · Page 90

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Coastal inundation regime moderates the short-term effects of sediment and soil additions on seawater oxygen and greenhouse gas dynamics: a microcosm experiment

The frequency and persistence of tidal inundation varies along the coastal terrestrial-aquatic interface, from frequently inundated wetlands to rarely inundated upland forests. This inundation gradient controls soil and sediment biogeochemistry and influence the exchange of soils and sediments from terrestrial to aquatic domains. Although a rich literature exist on studies of the influence of tidal waters on the biogeochemistry of coastal ecosystem soils, few studies have experimentally addressed the reverse question: How do soils (or sediments) from different coastal ecosystems influence the biogeochemistry of the tidal waters that inundate them? To better understand initial responses of coastal waters that flood coastal wetlands and uplands, we conducted short-term laboratory experiments where seawater was amended with sediments and soils collected across regional gradients of inundation exposure (i.e., frequently to rarely inundated) for 14 sites across the Mid-Atlantic, USA. Measured changes in dissolved oxygen and greenhouse gas concentrations were used to calculate gas consumption or production rates occurring during seawater exposure to terrestrial materials. We also measured soil and water physical and chemical properties to explore potential drivers. We observed higher oxygen consumption rates for seawater incubated with soils/sediments from frequently inundated locations and higher carbon dioxide production for seawater incubated with soils from rarely inundated transect locations. Incubations with soil from rarely inundated sites produced the highest global warming potential, primarily driven by carbon dioxide and secondarily by nitrous oxide. We also found environmental drivers of gas rates varied notably between transect locations. Our findings indicate that seawater responses to soil and sediment inputs across coastal terrestrial-aquatic interfaces exhibit some consistent patterns and high intra- and inter-site variability, suggesting potential biogeochemical feedback loops as inundation regimes shift inland.

coastal inundation↗

NASA NTRS1977

Stellar rotation and the thermomagnetic torque

It is noted that any reasonable assumption for the rotation rates of the interstellar clouds which must collapse to produce stars will result in stars which rotate much faster than suggested by observations if angular momentum is conserved throughout the collapse. The contribution of thermomagnetic torque, which can exist when a rarefied gas of polyatomic molecules is subjected to a magnetic field and a nonuniform temperature gradient, to the angular momentum of a collapsing protostellar object is examined in order to determine whether the magnitude of this torque is sufficient to slow down contracting protostellar clouds. Calculations are performed for two times during protostellar collapse when conditions may be favorable for thermomagnetic torque to play a role in angular-momentum considerations. The results indicate that this torque is apparently of no importance in determining or modifying stellar rotation during the early stages of stellar evolution.

Spear, G. G.↗

NASA NTRS1984

Parametric instabilities during electron cyclotron heating of tandem mirrors

Electron cyclotron resonance heating is one of the most commonly used methods of heating electrons in the plugs and in the thermal barriers of tandem mirrors. The intense coherent electromagnetic waves used for such heating are susceptible to parametric decay into other modes. Significant growth rates are found for the decay of either ordinary or extraordinary waves into two magnetized electron plasma waves. This and related effects may result in electron heating mechanisms rather different than those assumed in linear ray-tracing calculations. These results may help explain the unusual effects observed during heating of the Phaedrus tandem mirror device. In the general case, these instabilities may be strongly inhibited by density gradients.

Nicholson, D. R.↗

NASA NTRS1986

SSME thrust chamber modeling with Navier Stokes equations

The capability of predicting two-dimensional, compressible and reacting flow in the combustion chamber and nozzle of the Space Shuttle Main Engine (SSME) is demonstrated. A nonorthogonal body fitted coordinate system has been used to represent the combustor and nozzle geometry. The Navier-Stokes equations are solved for the entire thrust chamber with the k-epsilon turbulence model accounting for compressibility and large pressure gradients effects. Results of the computational test cases reveal all expected features of the transonic nozzle flows including location of sonic line, internal shock and boundary layer build-up. Calculated performance parameters such as thrust, flow rate, and specific impulse are also in reasonble agreement with available data. The results show promising potential of solving full Navier-Stokes equations with heat transfer and two-phase combustion in truly comprehensive modeling of rocket engines.

Przekwas, A. J.↗

NASA NTRS1988

Analysis of notched metal matrix composites under tension loading

Presented are techniques based on 3-D finite-element analysis for the examination of continuous fiber reinforced metal matrix composites. Examples are shown for specific metal matrix composites such as boron/aluminum and silicon-carbide/aluminum. Specimen stress-strain behavior and stress at first fiber failure were predicted for boron/aluminum laminates containing circular holes and crack-like slits. The predictions compared very well for (+ or - 45) sub 2s laminates. Mesh configuration was shown to have an effect on the calculation of stresses local to the notch. The presence of thin interface layers of matrix material had a significant influence on the slit-tip stress state, causing sharper stress gradients near the notch. Interface layers reduced the slit-tip fibers stresses in a (+ or - 45) sub s silicon-carbide/aluminum laminate but increased them in a (0/90) sub s laminate.

Bigelow, C. A.↗

NASA NTRS1990

Effect of Suction on Controlling the Secondary Instability of Boundary Layers

The effect of suction on controlling the 3-D secondary instability is investigated for a boundary layer with pressure gradient in the presence of small but finite amplitude Tollmien-Schlichting wave. The focus is on principal parametric resonance responsible for strong growth of subharmonics in low disturbance environment. Calculations are presented for the effect of suction on the onset and amplification of the secondary instability in Blasius and Falkner-Skan flows, as well as its effect on controlling the production of the vortical structure.

El-Hady, Nabil M.↗

NASA NTRS1991

Stochastic inflation lattice simulations - Ultra-large scale structure of the universe

Non-Gaussian fluctuations for structure formation may arise in inflation from the nonlinear interaction of long wavelength gravitational and scalar fields. Long wavelength fields have spatial gradients, a (exp -1), small compared to the Hubble radius, and they are described in terms of classical random fields that are fed by short wavelength quantum noise. Lattice Langevin calculations are given for a toy model with a scalar field interacting with an exponential potential where one can obtain exact analytic solutions of the Fokker-Planck equation. For single scalar field models that are consistent with current microwave background fluctuations, the fluctuations are Gaussian. However, for scales much larger than our observable Universe, one expects large metric fluctuations that are non-Gaussian. This example illuminates non-Gaussian models involving multiple scalar fields which are consistent with current microwave background limits.

Salopek, D. S.↗

NASA NTRS1993

An investigation of the atmospheric sources and sinks of Methyl Bromide

Methyl Bromide (CH3Br) is a ubiquitous component of the atmosphere and has been implicated as an important player in the depletion of stratospheric ozone. Atmospheric CH3Br abundances, interhemispheric gradients, oceanic concentrations, man-made emissions, and removal processes have been analyzed and interpreted with the help of a simple box model and a 2D global photochemical model. Its calculated atmospheric lifetime (T) of 1.7-1.9 years, based on reaction with OH radicals, is consistent with a global source of 90-110 Gg (10 exp 9 g)/year. Consequences of a much shorter lifetime of 1.2 years, due to possible deposition/hydrolysis losses, are also considered. Available data are used to estimate a CH3Br source that is 35 percent (20-50 percent) man-made and 65 percent (80-50 percent) natural. Oceans are substantially supersaturated and provide the most important natural source of about 60 (40-80) Gg/year.

Singh, Hanwant B.↗

NASA NTRS1994

Modeling of the hydrogen maser disk in MWC 349

Maser amplification in a Keplerian circumstellar disk seen edge on-the idea put forward by Gordon (1992), Martin-Pintado, & Serabyn (1992), and Thum, Martin-Pintado, & Bachiller (1992) to explain the millimeter hydrogen recombination lines in MWC 349-is further justified and developed here. The double-peaked (vs. possible triple-peaked) form of the observed spectra is explained by the reduced emission from the inner portion of the disk, the portion responsible for the central ('zero velocity') component of a triple-peaked spectrum. Radial gradient of electron density and/or free-free absorption within the disk are identified as the probable causes of this central 'hole' in the disk and of its opacity. We calculate a set of synthetic maser spectra radiated by a homogeneous Keplerian ring seen edge-on and compare them to the H30-alpha observations of Thum et al., averaged over about 1000 days. We used a simple graphical procedure to solve an inverse problem and deduced the probable values of some basic disk and maser parameters. We find that the maser is essentially unsaturated, and that the most probable values of electron temperature. Doppler width of the microturbulence, and electron density, all averaged along the amplification path are, correspondingly, T(sub e) less than or equal to 11,000 K, V(sub micro) less than or equal to 14 km/s, n(sub e) approx. = (3 +/- 2) x 10(exp 7)/cu cm. The model shows that radiation at every frequency within the spectrum arises in a monochromatic 'hot spot.' The maximum optical depth within the 'hot spot' producing radiation at the spectral peak maximum is tau(sub max) approx. = 6 +/- 1; the effective width of the masing ring is approx. = 0.4-0.7 times its outer diameter; the size of the 'hot spot' responsible for the radiation at the spectral peak frequency is approx. = 0.2-0.3 times the distance between the two 'hot spots' corresponding to two peaks. An important derivation of our model is the dynamical mass of the central star, M(sub *) approx. = 26 solar masses (D/1.2 kpc), D being the distance to the star. Prospects for improving the model are discussed.

Ponomarev, Victor O.↗

NASA NTRSDate not supplied

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗

NASA NTRS1978

Adhesion of a bimetallic interface

The Hohenberg-Kohn and Kohn-Sham formalisms are used to examine binding (binding energy as a function of separation) for combinations of the simple metals Al(111), Zn(0001), Mg(0001), and Na(110) in contact. Similar metal contacts between Al, Zn, Mg, and Na are examined self-consistently in an ab initio calculation using the Kohn-Sham formalism. Crystallinity is included using the Aschroft pseudopotential via first order perturbation theory for the electron-ion interaction; and the ion-ion interaction is included exactly via a lattice sum. Binding energy was determined both in the local-density approximation and including gradient corrections to the exchange and correlation energy. Binding was found in all cases. In dissimilar metal contacts, interfacial bonding was greater than that in the weaker material predicting the possibility of metallic transfer. The nonzero position of the energy minimum in like metal contacts is explained in terms of consistency between the Ashcroft pseudopotential and the bulk charge density. Good agreement with experimental surface energies is obtained in the self-consistent calculation when nonlocal terms are included.

Ferrante, J.↗

DOE OSTI2016

SU-F-T-316: A Model to Deal with Dosimetric and Delivery Uncertainties in Radiotherapy Treatment Planning

Purpose The conventional way of dealing with uncertainties resulting from dose calculation or beam delivery in IMRT, is to do verification measurements for the plan in question. Here we present an alternative based on recommendations given in the AAPM 142 report and treatment specific parameters that model the uncertainties for the plan delivery. Methods Basis of the model is the assignment of uncertainty parameters to all segment fields or control point sequences of a plan. The given field shape is analyzed for complexity, dose rate, number of MU, field size related output as well as factors for in/out field position and penumbra regions. Together with depth related uncertainties, a 3D matrix is generated by a projection algorithm. Patient anatomy is included as uncertainty CT data set as well. Therefore, object density is classified in 4 categories close to water, lung, bone and gradient regions with additional uncertainties. The result is then exported as a DICOM dose file by the software tool (written in IDL, Exelis), having the given resolution and target point. Results Uncertainty matrixes for several patient cases have been calculated and compared side by side in the planning system. The result is not quite always intuitive but it clearly indicates high and low uncertainties related to OARs and target volumes as well as to measured gamma distributions.ConclusionThe imported uncertainty datasets may help the treatment planner to understand the complexity of the treatment plan. He then might decide to change the plan to produce a more suited uncertainty distribution, e.g. by changing the beam angles the high uncertainty spots can be influenced or try to use another treatment setup, resulting in a plan with lower uncertainties. A next step could be to include such a model into the optimization algorithm to add a new dose uncertainty constraint.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

First-principles study of L -shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Ab initio phonon dispersion in crystalline naphthalene using van der Waals density functionals

Acene molecular crystals are of current interest in organic optoelectronics, both as active materials and for exploring and understanding new phenomena. Phonon scattering can be an important facilitator and dissipation mechanism in charge separation and carrier transport processes. Here, we carry out density functional theory (DFT) calculations of the structure and the full phonon dispersion of crystalline naphthalene, a well-characterized acene crystal for which detailed neutron-diffraction measurements, as well as infrared and Raman spectroscopy, are available. We evaluate the performance, relative to experiments, of DFT within the local density approximation (LDA); the generalized gradient approximation of Perdew, Burke, and Ernzerhof (PBE); and a recent van der Waals-corrected nonlocal correlation (vdW-DF-cx) functional. We find that the vdW-DF-cx functional accurately predicts lattice parameters of naphthalene within 1%. Intermolecular and intramolecular phonon frequencies across the Brillouin zone are reproduced within 7.8% and 1%, respectively. As expected, LDA (PBE) underestimates (overestimates) the lattice parameters and overestimates (underestimates) phonon frequencies, demonstrating their shortcomings for predictive calculations of weakly bound materials. If the unit cell is fixed to the experimental lattice parameters, PBE is shown to lead to improved phonon frequencies. Our study provides a detailed understanding of the phonon spectrum of naphthalene, and highlights the importance of including van der Waals dispersion interactions in predictive calculations of lattice parameters and phonon frequencies of molecular crystals and related organic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Theoretical prediction of nuclear magnetic shieldings and indirect spin-spin coupling constants in 1,1-, cis-, and trans-1,2-difluoroethylenes

A theoretical prediction of nuclear magnetic shieldings and indirect spin-spin coupling constants in 1,1-, cis- and trans-1,2-difluoroethylenes is reported. The results obtained using density functional theory (DFT) combined with large basis sets and gauge-independent atomic orbital calculations were critically compared with experiment and conventional, higher level correlated electronic structure methods. Accurate structural, vibrational, and NMR parameters of difluoroethylenes were obtained using several density functionals combined with dedicated basis sets. B3LYP/6-311++G(3df,2pd) optimized structures of difluoroethylenes closely reproduced experimental geometries and earlier reported benchmark coupled cluster results, while BLYP/6-311++G(3df,2pd) produced accurate harmonic vibrational frequencies. The most accurate vibrations were obtained using B3LYP/6-311++G(3df,2pd) with correction for anharmonicity. Becke half and half (BHandH) density functional predicted more accurate {sup 19}F isotropic shieldings and van Voorhis and Scuseria's τ-dependent gradient-corrected correlation functional yielded better carbon shieldings than B3LYP. A surprisingly good performance of Hartree-Fock (HF) method in predicting nuclear shieldings in these molecules was observed. Inclusion of zero-point vibrational correction markedly improved agreement with experiment for nuclear shieldings calculated by HF, MP2, CCSD, and CCSD(T) methods but worsened the DFT results. The threefold improvement in accuracy when predicting {sup 2}J(FF) in 1,1-difluoroethylene for BHandH density functional compared to B3LYP was observed (the deviations from experiment were −46 vs. −115 Hz)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1982

Numerical calculations relevant to the initial expansion of the polar wind

The results of a numerical simulation of the expansion of an H(+)-0(+) electron plasma into a vacuum are reported. Variations were introduced in the ionic density ratios, the scale lengths for the density gradient at the plasma-vacuum interface, and the initial electron-ion temperature ratio. Conditions of a gravitationally bound O(+) plasma were also considered. It was found that if H(+) is a major ion, the H(+) density profile in the expansion region is always concave, with a phase-space divided into three regions, i.e., an undisturbed plasma, a rarefaction region, and an expansion region in which ions are accelerated. When the H(+) is a minor ion, the O(+) ions accelerate the H(+) ions, the H(+) phase-space is divided into five regions, and the increased electron temperatures produce enhanced electrostatic potentials and H(+) drifts in the expansion region.

Singh, N.↗

DOE OSTI2018

A DFT Study on Structural, Electronic Mechanical and Thermodynamic Properties of 5f-Electron System BaAmO{sub 3}

The structural, electronic, mechanical and thermodynamic properties of the perovskite oxide BaAmO{sub 3} have been predicted using the full-potential linearized augmented plane wave (FP-LAPW) method. The equilibrium lattice constant, bulk modulus and pressure derivative were computed using different exchange correlations. The optimization of structure was carried out in ferromagnetic, anti-ferromagnetic and non-magnetic states, and the compound was found to be stable in the ferromagnetic state. A systematic study on the band structure and density of states was accomplished using generalized gradient approximation (GGA), Hubbard approximation (GGA + U) and modified Becke–Johnson exchange potential (mBJ),and the compound was found to have a half-metallic nature in all the approximations. The calculated total spin magnetic moment was found to be 5 μ{sub B} in all the approximations used. The second-order elastic constants, Young modulus, shear modulus, Poisson ratio and anisotropic factor have also been calculated. In order to have a complete understanding of BaAmO{sub 3}, the thermodynamic properties were studied in the pressure range of 0 to 40 GPa and the temperature range extending from 0 to 600 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Design and Application of CSC-Roll for Heavy Reduction of the Bloom Continuous Casting Process

Center quality control is a key issue affecting large-section bloom production and product quality. The central convex roll could avoid the bloom side part and yield a high reduction amount, which benefit the implement of heavy reduction (HR) for bloom continuous casting process without upgrading the mechanical and hydraulic system of withdrawal and straightening machine. In this work, a curving surface convex roll (CSC-Roll), which could significantly improve bloom center quality, was designed by finite element simulation for a 360 mm × 450 mm #45 steel continuous casting bloom. The optimal roll surface shape of the CSC-Roll was determined based on an analysis of the distribution of the reduction amount in the bloom center and the strain of the bloom centerline and inner surface. The calculation results showed that the optimal size of the CSC-Roll was 200 mm in convex platform region length, 30 mm in convex platform region height, and 80 mm in gradient curvature region length. The bloom heat transfer and deformation behavior in the CSC-Roll reduction process were then analyzed. The results showed that in the CSC-Roll HR process, 61.71 pct of the reduction amount could be transmitted to the inner side of the porosity zone. Macrographs of the 360 mm × 450 mm #45 steel continuous casting bloom proved that the bloom center quality was significantly improved after CSC-Roll HR was applied.

36 MATERIALS SCIENCE↗