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,567 records · Page 87

The effects of orbital and climatic variations on Martian surface heat flow

Large changes in the orbital elements of Mars on timescales of 10(exp 4) to 10(exp 6) years will cause widely varying climate, specifically surface temperatures, as a result of varying insolation. These surface temperature oscillations will produce subsurface thermal gradients which contribute to the total surface heat flux. We investigate the thermal behavior of the Martian regolith on orbital timescales and show that this climatological surface heat flux is spatially variable and contributes significantly to the total surface heat flux at many locations. We model the thermal behavior of the Martian regolith by calculating the mean annual surface temperatures for each epoch (spaced 1000 years apart to resolve orbital variations) for the past 200,000 years at a chosen location on the surface. These temperatures are used as a boundary condition for the deeper regolith and subsurface temperature oscillation are then computed. The surface climatological heat flux due to past climate changes can then be found from the temperature gradient between the surface and about 150 m depth (a fraction of the thermal skin depth on these timescales). This method provides a fairly accurate determination of the climatological heat flow component at a point; however, this method is computationally time consuming and cannot be applied to all points on the globe. To map the spatial variations in the surface heat flow we recognize that the subsurface temperature structure will be largely dominated by the most recent surface temperature oscillations. In fact, the climate component of the surface heat flow will be approximately proportional to the magnitude of the most recent surface temperature change. By calculating surface temperatures at all points globally for the present epoch and an appropriate past epoch, and combining these results with a series of more precise calculations described above, we estimate the global distribution of climatological surface heat flow.

Mellon, Michael T.↗

Multiscale analysis of the CMB temperature derivatives

We study the Planck CMB temperature at different scales through its derivatives up to second order, which allows one to characterize the local shape and isotropy of the field. The problem of having an incomplete sky in the calculation and statistical characterization of the derivatives is addressed in the paper. The analysis confirms the existence of a low variance in the CMB at large scales, which is also noticeable in the derivatives. Moreover, deviations from the standard model in the gradient, curvature and the eccentricity tensor are studied in terms of extreme values on the data. As it is expected, the Cold Spot is detected as one of the most prominent peaks in terms of curvature, but additionally, when the information of the temperature and its Laplacian are combined, another feature with similar probability at the scale of 10{sup o} is also observed. However, the p -value of these two deviations increase above the 6% when they are referred to the variance calculated from the theoretical fiducial model, indicating that these deviations can be associated to the low variance anomaly. Finally, an estimator of the directional anisotropy for spinorial quantities is introduced, which is applied to the spinors derived from the field derivatives. An anisotropic direction whose probability is <1% is detected in the eccentricity tensor.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-447: Electronic Brachytherapy (EBT) Treatment of Cervical Cancer - First Clinical Experience

Purpose: To study the first trial patient in which an electronic brachytherapy (EBT) x-ray source is utilized for treatment of cervical cancer. Methods: During patient treatment, a miniaturized x-ray source was used in combination with a customized titanium tandem and ovoid applicator set. The semi-specialized source was modeled with formalisms outlined by AAMP Task Group 43. Multiple models were used to compensate for variable attenuation conditions as a function of source positions. Varian Brachyvision treatment planning software was utilized on CT data sets for dose calculations prior to treatment delivery. The dose was prescribed to “point A” as defined by American Brachytherapy society. Additional treatments plans were created from those clinically utilized in patient care and were recalculated for an existing Ir-192 source model. Dose volume histograms (DVH) and point dose calculations were compared between the modalities for the clinical condition present in patients treated with EBT. Results: Clinical treatment times, though longer than those typically experienced by Ir-192 users, were manageable. Instantaneous dose rates at personal positions within the treatment vault were lower than those measured during intra operative radiation therapy and breast EBT treatments. Due to lower average photon energy in EBT, dose gradients within the treatment plans were as expected steeper than those observed in Ir-192 based brachytherapy. DVH comparisons between Ir-192 and EBT treatments showed an expected decrease in the integral dose to normal tissues of interest for EBT. In comparing plans created for EBT delivery with those calculated for Ir-192, average dose values for EBT were more than 4%, 11%, and 9% lower at predefined bladder, rectum and “point B” positions, respectively. Conclusion: For the first time, we have demonstrated that the utilizing electronic brachytherapy system for tandem and ovoid based treatment of cancer of the cervix is feasible, and potentially beneficial, in clinical environment.

60 APPLIED LIFE SCIENCES↗

Investigation of Electronic, Magnetic, Thermodynamic, and Thermoelectric Properties of Half-Metallic XLiSn (X = Ce, Nd) Alloys

We have investigated the electronic, magnetic, thermodynamic, and thermoelectric properties of half-Heusler CeLiSn and NdLiSn using density functional theory. The materials crystallize in F 4-bar 3mstructure with ferromagnetic state at optimized lattice parameter of 6.73 Å for CeLiSn and 6.74 Å for NdLiSn. The modified Becke-Johnson potential (mBJ) and generalized gradient approximation (GGA) have been used for computation of band structure, where mBJ shows a perfect half-metallic nature of the bands in both the alloys. The calculated magnetic moment of 1.0 μ{sub B} for CeLiSn and 3.0 μ{sub B} for NdLiSn follow Slater-Pauling rule, with ferromagnetic interactions. The temperature variation of the Seebeck coefficient shows the p-type behavior with room temperature values of 27 μV/K (43 μV/K) for CeLiSn (NdLiSn). Pressure-dependent variation of specific heat, Debye temperature, and Gruneisen parameter were calculated to infer the nature of lattice vibration and stability at high pressures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nonlinear, real-time optimization for actuator management in tokamaks

Experiments in DIII-D have been carried out to test a novel actuator management approach in tokamaks. Here, the actuator management scheme is posed as a nonlinear-optimization problem in which the actuator commands are calculated in real time according to the changing control priorities, plasma state, and actuator availability. Such optimization problem is solved using the augmented Lagrangian method, combined with a gradient projection method and a conjugate-gradient iteration algorithm. The algorithmic approach followed in this work does not depend on the particular control objectives or actuators considered, which facilitates its integration with other independently-designed control components within a plasma-control system. In addition, the actuator-management algorithm is able to handle the optimization problem in a computationally efficient manner, making it suitable for real-time implementations. Initial DIII-D results in the steady-state high-q min scenario have demonstrated the capabilities of the actuator manager to perform both simultaneous multiple mission and repurposing sharing, which will be required in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ion temperature gradient (ITG), trapped electron mode (TEM) and electron temperature gradient (ETG) modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. Finally, it is also confirmed that thermodiffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulating One-Photon Absorption and Resonance Raman Scattering Spectra Using Analytical Excited State Energy Gradients within Time-Dependent Density Functional Theory

A parallel implementation of analytical time-dependent density functional theory gradients is presented for the quantum chemistry program NWChem. The implementation is based on the Lagrangian approach developed by Furche and Ahlrichs. To validate our implementation, we first calculate the Stokes shifts for a range of organic dye molecules using a diverse set of exchange-correlation functionals (traditional density functionals, global hybrids, and range-separated hybrids) followed by simulations of the one-photon absorption and resonance Raman scattering spectrum of the phenoxyl radical, the well-studied dye molecule rhodamine 6G, and a molecular host–guest complex (TTFcCBPQT 4+ ). The study of organic dye molecules illustrates that B3LYP and CAM-B3LYP generally give the best agreement with experimentally determined Stokes shifts unless the excited state is a charge transfer state. Absorption, resonance Raman, and fluorescence simulations for the phenoxyl radical indicate that explicit solvation may be required for accurate characterization. For the host–guest complex and rhodamine 6G, it is demonstrated that absorption spectra can be simulated in good agreement with experimental data for most exchange-correlation functionals. Finally, however, because one-photon absorption spectra generally lack well-resolved vibrational features, resonance Raman simulations are necessary to evaluate the accuracy of the exchange-correlation functional for describing a potential energy surface.

Environmental Molecular Sciences Laboratory↗

Magnetophoretic potential at the movement of cluster products of electrochemical reactions in an inhomogeneous magnetic field

An electric field arises from the influence of a nonuniform static magnetic field on charged colloid particles with magnetic susceptibility different from that of the surrounding liquid. It arises, for example, under the influence of a nonuniform static magnetic field in clusters of electrochemical reaction products created during metal etching, deposition, and corrosion processes without an external electric current passing through an electrolyte near a magnetized electrode surface. The corresponding potential consists of a Nernst potential of inhomogeneous distribution of concentration of colloid particles and a magnetophoretic potential (MPP). This potential has been calculated using a thermodynamic approach based on the equations of thermodynamics of nonequilibrium systems and the Onsager relations for a mass flow of correlated magnetic clusters under a gradient magnetic force in the electrolyte. The conditions under which the MPP contribution to the total electric potential may be significant are discussed with a reference to the example of a corroding spherical ferromagnetic steel electrode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Concordance

Dimension reduction techniques have long been an important topic in statistics, and active subspaces (AS) have received much attention this past decade in the computer experiments literature. The most common approach towards estimating the AS is to use Monte Carlo with numerical gradient evaluation. This software takes a new approach towards estimating the active subspace, in which Multivariate Adaptive Regression Splines (MARS) is used as a surrogate model, and relevant calculation are obtained analytically. Using a MARS surrogate has many advantages including improved scaling, better estimation of active subspaces in high dimensions and the ability to handle a large number of prior distributions in closed form. This software is also capable of conducting a so-called "concordance analysis", using a generalization of ASs which we refer to as "Co-Active" subspaces.

Rumsey, Kelin↗

NRC Multiphysics Analysis Capability Deployment (FY2021 - Part 2)

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” Task 4b was completed for this report. INL developed a sample problem showing how to use Serpent 2 to calculate macroscopic cross sections for use in Griffin for a typical sodium fast reactor Unprotected Loss of Flow (ULOF) transient. The complicating factor is the relatively large axial gradient of the coolant (and hence fuel) temperature and its evolution during the transient. A 3D Griffin model of the Advanced Burner Test Reactor (ABTR) is coupled to the System Analysis Module (SAM) to perform the ULOF transient. The model includes various explicit feedback mechanisms for fast reactors, including Doppler, radial expansion from the displacement of the support plate, and axial expansion from the displacement of the fuel pins. The results for the ULOF transient are consistent with published values. We observed small effects in the transient results that arise from the superhomogenization equivalence correction of uniform and nonuniform temperature data sets, but they diminish as more dominant reactivity mechanisms are added to the model. Potential improvements to the neutronics model include: adding the differential control rod and sodium feedback mechanisms, implementing cross section re-homogenization based on the partial volume of the various materials, and enhancing the kinetics parameters. Furthermore, the fluids model only includes four subchannels to represent the entire active core due to performance limitations in SAM when attempting to model one subchannel per assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Longitudinal waves in a perpendicular collisionless plasma shock. IV - Gradient B.

The consideration of elastic waves in a Vlasov plasma of unmagnetized ions and magnetized electrons undergoing E x B electron drift and gradient B drift, pursued in the earlier three parts, is brought to conclusion in this last part of the longitudinal wave study in a collisionless plasma shock. Detailed calculations of the effects of the beta sub e dimensionless parameter on the E x B electron drift instability are presented. It is shown that the range of propagation of the elastic waves about the perpendicular remains quite narrow, and that, for oblique propagation, the already narrow angular range of unstable waves is decreased by increases in the value of the beta sub e dimensionless parameter. Also, increases in wave number generally reduce the growth rate and the angular range of propagation.

Gary, S. P.↗

The five-minute oscillations as nonradial pulsations of the entire sun.

Calculations of the stability coefficient show that the Sun should be pulsating as a unit in nonradial modes of high order. The pulsations are driven by the superadiabatic gradient of the low photosphere and by the same sensitive changes in opacity that are known to be important in variable stars. The existence of solar pulsations can explain many of the large scale features observed in the well known 5-minute oscillations of the solar atmosphere.

Wolff, C. L.↗

Kinetics of boron deposition by the hydrogen reduction of boron trichloride

The chemical vapor deposition of boron filament is generally considered to be limited by the diffusion of the reactive species to the hot substrate. From concentration gradient considerations and from experimental measurements it is concluded that the product species (HCl) is more likely to be involved in limiting the deposition rate. Comparison of theoretically calculated and experimentally measured deposition rates show that the rate is generally partially limited by diffusion and partially by the surface reaction rate (desorption) of HCl. The relative amounts of surface reaction or diffusion rate impedance were determined as a function of temperature and gas velocity. The temperature coefficient for the desorption of HCl was found to be 52 Kcal/mole.

Mehalso, R. M.↗

Persistent sunward flow of approximately 1.6 MeV protons at 1 AU

The paper is concerned with the anisotropy of 1.3 to 2.3 MeV protons measured with an electron/isotope spectrometer aboard the IMP 7 satellite during 317 six-hour periods, excluding periods when prompt solar particle events are observed. It is found that the diffusive streaming is predominantly toward the sun, indicating an essentially positive radial gradient in the particle density. The long-term average radial component of this anisotropy and its dependence on the solar wind speed are consistent with the anisotropy calculated with a steady-state propagation model including adiabatic energy loss.

Marshall, F. E.↗

Flow velocities and streamlines

Computer program calculates subsonic or transonic flow on hubshroud, midchannel, stream surface of single-blade row of turbomachine. Program uses finite-different and quasi-orthogonal (velocity-gradient) methods. Program is reported in two volumes: Part I is User's Manual, Part II is Programmer's Manual.

Katsanis, T.↗

Decay of the zincate concentration gradient at an alkaline zinc cathode after charging

The transport of the zincate ion to the alkaline zinc cathode was studied by observing the decay of the zincate concentration gradient at a horizontal zinc cathode after charging. This decay was found to approximate first order kinetics as expected from a proposed boundary layer model. The concentrations were calculated from polarization voltages. The decay half life was shown to be a linear function of the thickness of porous zinc deposit on the cathode indicating a very rapid transport of zincate through porous zinc metal. The rapid transport is attributed to an electrochemical mechanism. From the linear dependence of the half life on the thickness the boundary layer thickness was found to be about 0.010 cm when the cathode was at the bottom of the cell. No significant dependence of the boundary layer thickness on the viscosity of electrolyte was observed. The data also indicated a relatively sharp transition between the diffusion and convection transport regions. When the cathode was at the top of the cell, the boundary layer thickness was found to be roughly 0.080 cm. The diffusion of zincate ion through asbestos submerged in alkaline electrolyte was shown to be comparable with that predicted from the bulk diffusion coefficient of the zincate ion in alkali.

Kautz, H. E.↗

The relative importance of aerosol scattering and absorption in remote sensing

The relative importance of aerosol optical thickness and absorption is illustrated through computing radiances for radiative transfer models. The radiance of sunlight reflected from models of the earth-atmosphere system is computed as a function of the aerosol optical thickness and its albedo of single scattering; it is noted that the albedo varies from 0.6 in urban environment to nearly 1 in areas with low graphitic carbon content. The calculations are applied to the example of satellite measurements of biomass. It is found that when surface classifications are made by means of clustering techniques the presence of gradients in the aerosol optical properties results in the dispersion of points in the plot correlating radiances viewed in two different directions. Finally, though such a remote sensing parameter as contrast is weakly affected by aerosol absorption, it is highly dependent on its optical thickness.

Fraser, R. S.↗