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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Application of a new K-tau model to near wall turbulent flows

A recently developed K-tau model for near wall turbulent flows is applied to two severe test cases. The turbulent flows considered include the incompressible flat plate boundary layer with the adverse pressure gradients and incompressible flow past a backward facing step. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.

Thangam, S.↗
NASA NTRS2001

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗
NASA NTRS2007

Analyzing Radio-Frequency Coverage for the ISS

The Interactive Coverage Analysis Tool (iCAT) is an interactive desktop computer program serving to (1) support planning of coverage, and management of usage of frequencies, of current and proposed radio communication systems on and near the International Space Station (ISS) and (2) enable definition of requirements for development of future such systems. The iCAT can also be used in design trade studies for other (both outer-space and terrestrial) communication systems. A user can enter the parameters of a communication-system link budget in a table in a worksheet. The nominal (onaxis) link values for the bit-to-noise-energy ratio, received isotropic power (RIP), carrier-to-noise ratio (C/N), power flux density (PFD), and link margin of the system are calculated and displayed in the table. Plots of field gradients for the RIP, C/N, PFD, and link margin are constructed in an ISS coordinate system, at a specified link range, for both the forward and return link parameters, and are displayed in worksheets. The forward and reverse link antenna gain patterns are also constructed and displayed. Line-of-sight (LOS) obstructions can be both incorporated into the gradient plots and displayed on separate plots.

Bolen, Steven M.↗
NASA NTRS2011

Relaxing USOS Solar Array Constraints for Russian Vehicle Undocking

With the retirement of Space Shuttle cargo delivery capability and the ten year life extension of the International Space Station (ISS) more emphasis is being put on preservation of the service life of ISS critical components. Current restrictions on the United States Orbital Segment (USOS) Solar Array (SA) positioning during Russian Vehicle (RV) departure from ISS nadir and zenith ports cause SA to be positioned in the plume field of Service Module thrusters and lead to degradation of SAs as well as potential damage to Sun tracking Beta Gimbal Assemblies (BGA). These restrictions are imposed because of the single fault tolerant RV Motion Control System (MCS), which does not meet ISS Safety requirements for catastrophic hazards and dictates 16 degree Solar Array Rotary Joint position, which ensures that ISS and RV relative motion post separation, does lead to collision. The purpose of this paper is to describe a methodology and the analysis that was performed to determine relative motion trajectories of the ISS and separating RV for nominal and contingency cases. Analysis was performed in three phases that included ISS free drift prior to Visiting Vehicle separation, ISS and Visiting Vehicle relative motion analysis and clearance analysis. First, the ISS free drift analysis determined the worst case attitude and attitude rate excursions prior to RV separation based on a series of different configurations and mass properties. Next, the relative motion analysis calculated the separation trajectories while varying the initial conditions, such as docking mechanism performance, Visiting Vehicle MCS failure, departure port location, ISS attitude and attitude rates at the time of separation, etc. The analysis employed both orbital mechanics and rigid body rotation calculations while accounting for various atmospheric conditions and gravity gradient effects. The resulting relative motion trajectories were then used to determine the worst case separation envelopes during the clearance analysis. Analytical models were developed individually for each stage and the results were used to build initial conditions for the following stages. In addition to the analysis approach, this paper also discusses the analysis results, showing worst case relative motion envelopes, the recommendations for ISS appendage positioning and the suggested approach for future analyses.

Menkin, Evgeny↗
NASA NTRS2014

Aerosol Polarimetry Sensor (APS): Design Summary, Performance and Potential Modifications

APS is a mature design that has already been built and has a TRL of 7. Algorithmic and retrieval capabilities continue to improve and make better and more sophisticated used of the data. Adjoint solutions, both in one dimensional and three dimensional are computationally efficient and should be the preferred implementation for the calculation of Jacobians (one dimensional), or cost-function gradients (three dimensional). Adjoint solutions necessarily provide resolution of internal fields and simplify incorporation of active measurements in retrievals, which will be necessary for a future ACE mission. Its best to test these capabilities when you know the answer: OSSEs that are well constrained observationally provide the best place to test future multi-instrument platform capabilities and ensure capabilities will meet scientific needs.

algorithms↗
NASA NTRS2005

Computing Gravitational Fields of Finite-Sized Bodies

A computer program utilizes the classical theory of gravitation, implemented by means of the finite-element method, to calculate the near gravitational fields of bodies of arbitrary size, shape, and mass distribution. The program was developed for application to a spacecraft and to floating proof masses and associated equipment carried by the spacecraft for detecting gravitational waves. The program can calculate steady or time-dependent gravitational forces, moments, and gradients thereof. Bodies external to a proof mass can be moving around the proof mass and/or deformed under thermoelastic loads. An arbitrarily shaped proof mass is represented by a collection of parallelepiped elements. The gravitational force and moment acting on each parallelepiped element of a proof mass, including those attributable to the self-gravitational field of the proof mass, are computed exactly from the closed-form equation for the gravitational potential of a parallelepiped. The gravitational field of an arbitrary distribution of mass external to a proof mass can be calculated either by summing the fields of suitably many point masses or by higher-order Gauss-Legendre integration over all elements surrounding the proof mass that are part of a finite-element mesh. This computer program is compatible with more general finite-element codes, such as NASTRAN, because it is configured to read a generic input data file, containing the detailed description of the finiteelement mesh.

Quadrelli, Marco↗
DOE OSTI2014

Beam-driven acceleration in ultra-dense plasma media

Accelerating parameters of beam-driven wakefield acceleration in an extremely dense plasma column has been analyzed with the dynamic framed particle-in-cell plasma simulator, and compared with analytic calculations. In the model, a witness beam undergoes a TeV/m scale alternating potential gradient excited by a micro-bunched drive beam in a 10{sup 25 }m{sup −3} and 1.6 × 10{sup 28 }m{sup −3} plasma column. The acceleration gradient, energy gain, and transformer ratio have been extensively studied in quasi-linear, linear-, and blowout-regimes. The simulation analysis indicated that in the beam-driven acceleration system a hollow plasma channel offers ∼20% higher acceleration gradient by enlarging the channel radius (r) from 0.2 λ{sub p} to 0.6 λ{sub p} in a blowout regime. This paper suggests a feasibility of TeV/m scale acceleration with a hollow crystalline structure (e.g., nanotubes) of high electron plasma density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Projection-operator methods for classical transport in magnetized plasmas. Part 2. Nonlinear response and the Burnett equations [Projection-operator methods for classical transport in magnetized plasmas. Part 1. Nonlinear response and the Burnett equations]

The time-independent projection-operator formalism for the derivation of Burnett equations is extended and considered in the context of multispecies and magnetized plasmas. The procedure provides specific formulas for transport coefficients in terms of two-time correlation functions involving both two and three phase-space points. It is shown how to calculate those correlation functions in the limit of weak coupling. The results are used to demonstrate, with the aid of a particular non-trivial example, that the Chapman–Enskog methodology employed by Catto & Simakov (CS) (Phys. Plasmas, vol. 11, 2004, pp. 90–102) to calculate the contributions to the parallel viscosity driven by temperature gradients is consistent with formulas previously derived from the two-time formalism by Brey (J. Chem. Phys., vol. 79, 1983, pp. 4585–4598). Here, the work serves to unify previous work on plasma kinetic theory with formalism usually applied to turbulence. Additional contributions include discussions of (i) Braginskii-order interspecies momentum exchange from the point of view of two-time correlations; and (ii) a simple stochastic model, unrelated to many-body theory, that exhibits Burnett effects. Insights from that model emphasize the role of non-Gaussian statistics in the evaluation of Burnett transport coefficients, including the effects calculated by CS that stem from the nonlinear collision operator. Together, Parts 1 and 2 of this series provide an introduction to projection-operator methods that should be broadly useful in theoretical plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

SU-E-T-196: Comparative Analysis of Surface Dose Measurements Using MOSFET Detector and Dose Predicted by Eclipse - AAA with Varying Dose Calculation Grid Size

Purpose: Skin dose can be the limiting factor and fairly common reason to interrupt the treatment, especially for treating head-and-neck with Intensity-modulated-radiation-therapy(IMRT) or Volumetrically-modulated - arc-therapy (VMAT) and breast with tangentially-directed-beams. Aim of this study was to investigate accuracy of near-surface dose predicted by Eclipse treatment-planning-system (TPS) using Anisotropic-Analytic Algorithm (AAA)with varying calculation grid-size and comparing with metal-oxide-semiconductor-field-effect-transistors(MOSFETs)measurements for a range of clinical-conditions (open-field,dynamic-wedge, physical-wedge, IMRT,VMAT). Methods: QUASAR™-Body-Phantom was used in this study with oval curved-surfaces to mimic breast, chest wall and head-and-neck sites.A CT-scan was obtained with five radio-opaque markers(ROM) placed on the surface of phantom to mimic the range of incident angles for measurements and dose prediction using 2mm slice thickness.At each ROM, small structure(1mmx2mm) were contoured to obtain mean-doses from TPS.Calculations were performed for open-field,dynamic-wedge,physical-wedge,IMRT and VMAT using Varian-21EX,6&15MV photons using twogrid-sizes:2.5mm and 1mm.Calibration checks were performed to ensure that MOSFETs response were within ±5%.Surface-doses were measured at five locations and compared with TPS calculations. Results: For 6MV: 2.5mm grid-size,mean calculated doses(MCD)were higher by 10%(±7.6),10%(±7.6),20%(±8.5),40%(±7.5),30%(±6.9) and for 1mm grid-size MCD were higher by 0%(±5.7),0%(±4.2),0%(±5.5),1.2%(±5.0),1.1% (±7.8) for open-field,dynamic-wedge,physical-wedge,IMRT,VMAT respectively.For 15MV: 2.5mm grid-size,MCD were higher by 30%(±14.6),30%(±14.6),30%(±14.0),40%(±11.0),30%(±3.5)and for 1mm grid-size MCD were higher by 10% (±10.6), 10%(±9.8),10%(±8.0),30%(±7.8),10%(±3.8) for open-field, dynamic-wedge, physical-wedge, IMRT, VMAT respectively.For 6MV, 86% and 56% of all measured values agreed better than ±20% for 1mm and 2.5mm grid-sizes respectively. For 18MV, 56% and 18% of all measured-values agreed better than ±20% for 1mm and 2.5mm grid-sizes respectively. Conclusion: Reliable Skin-dose calculations by TPS can be very difficult due to steep dose-gradient and inaccurate beam-modelling in buildup region.Our results showed that Eclipse over-estimates surface-dose.Impact of grid-size is also significant,surface-dose increased up to 40% from 1mm to 2.5mm,however, 1mm calculated-values closely agrees with measurements. Due to large uncertnities in skin-dose predictions from TPS, outmost caution must be exercised when skin dose is evaluated,a sufficiently smaller grid-size(1mm)can improve the accuracy and MOSFETs can be used for verification.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2021

Satellite Remote Sensing Observations of Trans-Atlantic Dust Transport and Deposition: A Multi-Sensor Analysis

We analyze the decade-long (2007-2016) record of aerosol measurements from four distinctive sensors, namely CALIOP, MODIS, MISR, and IASI, to quantify the trans-Atlantic dust transport and deposition. These satellite sensors use different techniques to characterize particle size and shape properties; and we have developed sensor-specific methods (broadly categorized into size-based and shape-based method) to derive dust optical depth (DOD). The size-based DOD from MODIS and IASI generally agrees better with AERONET-derived DOD than the shape-based DOD from CALIOP and MISR does. Overall, the shape-based DOD is smaller than the size-based DOD by about 25%, which is consistent with distinctive ways of accounting for coarse-spherical particles of dust-pollution internal mixture. While such dust-pollution mixtures are counted as dust in the size-based DOD, they are excluded in the shape-based DOD. DOD is not a good proxy for dust deposition. Instead, the dust deposition depends strongly on the gradient of DOD on a monthly basis and can be derived by calculating the meridional and zonal dust mass flux based on the three-dimensional distributions of dust. Among the remote sensing measurements, difference in dust deposition is smaller than that of DOD, suggesting that different satellites characterize the DOD gradient more consistently than DOD itself. Satellite measurements of dust deposition and DOD also provide an accurate estimate of the dust loss frequency (LF) that measures how efficient the dust is removed from the atmosphere. We found that these remote sensing measurements yield similar LF values of 0.078 –0.102 d-1, which however is factors of 2-5 smaller than model simulations. This analysis provides valuable insights into potential deficiencies in models’ emission and transport/removal processes and hence helps guide model improvement.

remote sensing↗
DOE OSTI2020

Constitutional under-potential plating (CUP) – New insights for predicting the morphological stability of deposited lithium anodes in lithium metal batteries

Lithium (Li) cycling in lithium metal batteries (LMBs) inevitably witnesses a non-planar Li/electrolyte interface (cellular/dendritic microstructure) affecting the performance, safety and reliability of LMBs. Herein, the “constitutional under-potential plating” (CUP) concept is presented for the first time to understand the evolution, formation and predict the morphological stability/instability of the Li/electrolyte interface. Here, the CUP condition, a measure of the driving force for interface perturbation, occurs when the true potential gradient at the interface is smaller than the equilibrium electrodeposition potential gradient at the Li/electrolyte interface resulting in the electrolyte ahead of the Li-electrolyte interface to be in an under-potential state. Detailed CUP criterion calculations conducted herein establishes the relationship of battery operating conditions (e.g. current density, potential gradient) with the electrolyte properties (e.g. Li-ion transference number and diffusion coefficient) predicting the morphological stability condition at the interface. In addition to the driving force for interface perturbation, the capillarity effect or the Gibbs Thomson effect, a measure of the resistance to perturbation of the perturbed interface is discussed for understanding the influence of interfacial energy/adhesion energy between the current collector and the deposited Li on the overall morphological stability of the Li/electrolyte interface, a critical criterion for the unwanted dendritic morphology.

25 ENERGY STORAGE↗
NASA NTRS1995

Understanding the large-scale structure of magnetic fields in the corona and heliosphere

To understand the non-radial nature of the southwestern streamer in the 3 November 1994 eclipse and the near absence of a latitudinal gradient of the interplanetary magnetic field strength seen by Ulysses over the southern pole, we calculate the coronal and inter planetary magnetic field from WSO photospheric field data using both the horizontal current-current sheet model (Zhao and Hoeksema, 1994) and the current sheet-source surface model (Zhao and Hoeksema, 1995). These models account for the effects of large-currents and current sheets in the corona. Comparison of the calculated results with white light, X-ray, and IMF observations show that non-radial coronal streamers and a near-zero latitudinal gradient of the radial component of the IMF may be explained by a large-scale current sheet. To best match the observations, the height of the cusp point (above which current sheets form) differs from one streamer to another. While heliospheric volume currents flowing above the source surface certainly affect the interplanetary field, they have little or no influence on the structure of the corona.

Zhao, X. P.↗
NASA NTRS1990

Spectra of circularly polarized radiation from astrophysical OH masers

A striking feature of astrophysical masers is the tendency for either one or the other of the circular polarizations to dominate in the radiation from the strong, widely observed masing transitions of OH at 18 cm. Spectral line profiles are calculated for polarized maser radiation due to the combined effects of a velocity gradient and, as is indicated for these transitions, a Zeeman splitting that is at least comparable with the thermal contributions to the breadths of the spectral lines. The resulting spectral features are similar in appearance, including the presence of large net circular polarization and narrow line breadths, to the commonly observed spectra of OH masers in molecular clouds. The calculations presented here are performed as a function of frequency without making the approximations of a large velocity gradient. Rapid cross relaxation, which has been advocated by others for the OH masers, is assumed.

Nedoluha, Gerald E.↗
NASA NTRS1973

Unsteady turbulent boundary layers in two-dimensional, incompressible flow.

The paper presents the results of some calculations aimed at studying the characteristics of turbulent boundary layers in unsteady flow. The calculations were performed using the method of Singleton and Nash, which is based on the turbulent kinetic energy equation. Emphasis is placed on the approach to zero wall stress, and the results show the extent to which onset of this condition can be delayed by the effects of time dependence. The delay is not associated merely with alleviation of the pressure gradients, and quasi-steady methods cannot readily be used to calculate this type of flow.

Nash, J. F.↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

Thermal energy and charge currents in multi-terminal nanorings

We study in experiment and theory thermal energy and charge transfer close to the quantum limit in a ballistic nanodevice, consisting of multiply connected one-dimensional electron waveguides. The fabricated device is based on an AlGaAs/GaAs heterostructure and is covered by a global top-gate to steer the thermal energy and charge transfer in the presence of a temperature gradient, which is established by a heating current. The estimate of the heat transfer by means of thermal noise measurements shows the device acting as a switch for charge and thermal energy transfer. The wave-packet simulations are based on the multi-terminal Landauer-Büttiker approach and confirm the experimental finding of a mode-dependent redistribution of the thermal energy current, if a scatterer breaks the device symmetry.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

MO-F-16A-05: Spot-Scanning Proton Therapy Patient-Specific Quality Assurance: Our Methodology and Results From 295 Treatment Plans

Purpose: To report on the patient-specific quality assurance (PSQA) results for 295 spot-scanning proton therapy treatment plans from the MD Anderson PTC-Houston. We show how the results differed by treatment site and how they were affected by the treatment plan optimization method and by a range shifter in the treatment field. We also discuss some causes of PSQA problems. Methods: The PSQA procedure, which is designed to verify both the accuracy of the treatment planning system's (Eclipse™ v8.9) dose calculations and the dose delivery of the Hitachi PROBEAT synchrotron, consists of (1) an end-to-end test in which the beam is delivered and measured at the prescribed gantry angle, and (2) additional dose plane measurements made from gantry angle 270°. HPlusQA™ software automatically performs the gamma analysis with criteria 3% (dose tolerance), 3 mm (distance-to-agreement, DTA) and 2%, 2 mm. Passing is defined as at least 90% of the pixels having a gamma score less than 1. Results: The PSQA gamma passing rate was 96.2% for 3%, 3 mm, and 85.3% for 2%, 2 mm. The rate depended on the treatment site. For example, the 3%, 3 mm passing rate was 95% for head and neck plans, vs 100% for prostate plans. The passing rates of multi- vs. single-field optimization plans did not significantly differ. However, the rate for fields with range shifters was 94.8±0.6%, vs 99.0±0.6% for those without (p = 0.002). Longitudinal dose gradients caused most of the low scores. Overestimation of the calculated dose proximal to the spread-out Bragg peak (SOBP) caused many of the others. Conclusion: The planned and delivered doses consistently agreed within tolerance levels. Minor dose modeling deficiencies remain proximal to the SOBP. The 3% dose tolerance, 3 mm DTA, with 90% pixel passing rate is a reasonable action level for 2D gamma comparisons.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1998

Operational Data Reduction Procedure for Determining Density and Vertical Structure of the Martian Upper Atmosphere from Mars Global Surveyor Accelerometer Measurements

The success of aerobraking by the Mars Global Surveyor (MGS) spacecraft was partly due to the analysis of MGS accelerometer data. Accelerometer data was used to determine the effect of the atmosphere on each orbit, to characterize the nature of the atmosphere, and to predict the atmosphere for future orbits. To interpret the accelerometer data, a data reduction procedure was developed to produce density estimations utilizing inputs from the spacecraft, the Navigation Team, and pre-mission aerothermodynamic studies. This data reduction procedure was based on the calculation of aerodynamic forces from the accelerometer data by considering acceleration due to gravity gradient, solar pressure, angular motion of the MGS, instrument bias, thruster activity, and a vibration component due to the motion of the damaged solar array. Methods were developed to calculate all of the acceleration components including a 4 degree of freedom dynamics model used to gain a greater understanding of the damaged solar array. The total error inherent to the data reduction procedure was calculated as a function of altitude and density considering contributions from ephemeris errors, errors in force coefficient, and instrument errors due to bias and digitization. Comparing the results from this procedure to the data of other MGS Teams has demonstrated that this procedure can quickly and accurately describe the density and vertical structure of the Martian upper atmosphere.

Cancro, George J.↗