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,495 records · Page 83

Orbital parameters for the two young binaries VSB 111 and VSB 126

We report orbital parameters for two low-mass, pre-main-sequence, double-lined spectroscopic binaries: VSB 111 and VSB 126. These systems were originally identified as single-lined on the basis of visible-light observations. We obtained high-resolution infrared spectra with the 10 m Keck II telescope, detected absorption lines of the secondary stars, and measured radial velocities of both components in the systems. The visible-light spectra were obtained with the 1.5 m Wyeth reflector at the Oak Ridge Observatory, the 1.5 m Tillinghast reflector at the F. L. Whipple Observatory, and the 4.5 m equivalent Multiple Mirror Telescope. The combination of our visible and infrared observations of VSB 111 leads to a period of 902.1 ± 0.9 days, an eccentricity of 0.788 ± 0.008, and a mass ratio of 0.52 ± 0.05. VSB 126 has a period of 12.9244 ± 0.0002 days, an eccentricity of 0.18 ± 0.02, and a mass ratio of 0.29 ± 0.02. Visible-light photometry, using the 0.8 m telescope at Lowell Observatory, provided rotation periods for the primary stars in both systems: 3.74 ± 0.02 days for VSB 111 and 5.71 ± 0.07 days for VSB 126. Both binaries are located in the young, active star-forming cluster NGC 2264 at a distance of ∼800 pc. The difference in the center-of-mass velocities of the two systems is consistent with the radial velocity gradient seen across NGC 2264. To test the evolutionary models for accuracy and consistency, we compare the stellar properties derived from several sets of theoretical calculations for pre-main-sequence evolution with our dynamical results.

79 ASTRONOMY AND ASTROPHYSICS↗

Full-sky formulae for weak lensing power spectra from total angular momentum method

We systematically derive full-sky formulae for the weak lensing power spectra generated by scalar, vector and tensor perturbations from the total angular momentum (TAM) method. Based on both the geodesic and geodesic deviation equations, we first give the gauge-invariant expressions for the deflection angle and Jacobi map as observables of the CMB lensing and cosmic shear experiments. We then apply the TAM method, originally developed in the theoretical studies of CMB, to a systematic derivation of the angular power spectra. The TAM representation, which characterizes the total angular dependence of the spatial modes projected along a line-of-sight, can carry all the information of the lensing modes generated by scalar, vector, and tensor metric perturbations. This greatly simplifies the calculation, and we present a complete set of the full-sky formulae for angular power spectra in both the E-/B-mode cosmic shear and gradient-/curl-mode lensing potential of deflection angle. Based on the formulae, we give illustrative examples of non-vanishing B-mode cosmic shear and curl-mode of deflection angle in the presence of the vector and tensor perturbations, and explicitly compute the power spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Verification and validation of integrated simulation of energetic particles in fusion plasmas

Here, we report on verification and validation of linear simulations of Alf´ven eigenmodes in the current ramp phase of DIII-D L-mode discharge #159243 using gyrokinetic, gyrokinetic-MHD hybrid, and eigenvalue codes. Using a classical fast ion profile, all simulation codes find that reversed shear Alf´ven eigenmodes (RSAE) are the dominant instability. The real frequencies from all codes have a coefficient of variation of less than 5% for the most unstable modes with toroidal mode number n = 4 and 5. Due to uncertainty in the experimentally measured minimum safety factor (q min ) and the sensitivity of the RSAE’s frequency to q min , the average difference in frequency between simulated and measured is df/f = 14.7%, 22.6%, for n = 4 and 5, respectively. If qmin is adjusted slightly, within experimental errors, an agreement of within 4.5% and 3.1%, for n = 4 and 5 respectively, is found for all codes. The simulated growth rates exhibit greater variation, and simulations find that pressure gradients of thermal plasmas make a significant contribution to the growth rates. Mode structures of the dominant modes agree well among all codes. Using a calculated fast ion profile that takes into account the action of the spectrum of multiple unstable modes, a toroidal Alf´ven eigenmode (TAE) with n = 6 is found to be unstable in the outer edge, consistent with the experimental data. Variations of the real frequencies and growth rates of the TAE are slightly larger than those of the RSAE. Lastly, electron temperature fluctuations and radial phase shifts from simulations show no significant differences with the experimental data for the strong n = 4 RSAE, but significant differences for the weak n = 6 TAE. The verification and validation for the linear simulations of Alf´ven eigenmodes is the first step to determine an integrated simulation of energetic particles confinement in burning plasmas incorporating multiple physical processes with kinetic effects of both energetic and thermal particles on an equal footing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A log-barrier Newton-CG method for bound constrained optimization with complexity guarantees

Abstract We describe an algorithm based on a logarithmic barrier function, Newton’s method and linear conjugate gradients that seeks an approximate minimizer of a smooth function over the non-negative orthant. We develop a bound on the complexity of the approach, stated in terms of the required accuracy and the cost of a single gradient evaluation of the objective function and/or a matrix-vector multiplication involving the Hessian of the objective. The approach can be implemented without explicit calculation or storage of the Hessian.

O’Neill, Michael↗

SU-E-T-226: Junction Free Craniospinal Irradiation in Linear Accelerator Using Volumetric Modulated Arc Therapy : A Novel Technique Using Dose Tapering

Purpose: Spatially separated fields are required for craniospinal irradiation due to field size limitation in linear accelerator. Field junction shits are conventionally done to avoid hot or cold spots. Our study was aimed to demonstrate the feasibility of junction free irradiation plan of craniospinal irradiation (CSI) for Meduloblastoma cases treated in linear accelerator using Volumetric modulated arc therapy (VMAT) technique. Methods: VMAT was planned using multiple isocenters in Monaco V 3.3.0 and delivered in Elekta Synergy linear accelerator. A full arc brain and 40° posterior arc spine fields were planned using two isocentre for short (<1.3 meter height ) and 3 isocentres for taller patients. Unrestricted jaw movement was used in superior-inferior direction. Prescribed dose to PTV was achieved by partial contribution from adjacent beams. A very low dose gradient was generated to taper the isodoses over a long length (>10 cm) at the conventional field junction. Results: In this primary study five patients were planned and three patients were delivered using this novel technique. As the dose contribution from the adjacent beams were varied (gradient) to create a complete dose distribution, therefore there is no specific junction exists in the plan. The junction were extended from 10–14 cm depending on treatment plan. Dose gradient were 9.6±2.3% per cm for brain and 7.9±1.7 % per cm for spine field respectively. Dose delivery error due to positional inaccuracy was calculated for brain and spine field for ±1mm, ±2mm, ±3mm and ±5 mm were 1%–0.8%, 2%–1.6%, 2.8%–2.4% and 4.3%–4% respectively. Conclusion: Dose tapering in junction free CSI do not require a junction shift. Therefore daily imaging for all the field is also not essential. Due to inverse planning dose to organ at risk like thyroid kidney, heart and testis can be reduced significantly. VMAT gives a quicker delivery than Step and shoot or dynamic IMRT.

60 APPLIED LIFE SCIENCES↗

Simple methods for determining the virtual origin of turbulent boundary layers in hypersonic flow on sharp-edged flat plates

Two methods for determining the virtual origin of turbulent boundary layers in hypersonic flow are evaluated. The results of the analyses are restricted to wind-tunnel models having sharp-edged surfaces with zero or small pressure gradients. Virtual origin and skin friction estimates from these two methods are compared with values from a base method for which the virtual origin is calculated from the measured momentum thickness at a station downstream of boundary layer transition.

Hopkins, E. J.↗

Entopistic scintillators

Trajectories of cosmic rays through a balloon-borne experiment will be determined by pulse height analyzing the signals from each of the two or more phototubes coupled to each scintillator. Thin or tapered scintillators produce sufficiently strong gradients of phototube signal with respect to scintillation position that position uncertainties of about 1 cm can be achieved over 1-sq-m areas. Calculations and some test results for entopistic, or position determining, scintillators are presented.

Arens, J. F.↗

Isothermal blast wave model of supernova remnants

The validity of the 'adiabatic' assumption in supernova-remnant calculations is examined, and the alternative extreme of an isothermal blast wave is explored. It is concluded that, because of thermal conductivity, the large temperature gradients predicted by the adiabatic model probably are not maintained in nature. Self-similar solutions to the hydrodynamic equations for an isothermal blast wave have been found and studied. These solutions are then used to determine the relationship between X-ray observations and inferred parameters of supernova remnants. A comparison of the present results with those for the adiabatic model indicates differences which are less than present observational uncertainties. It is concluded that most parameters of supernova remnants inferred from X-ray measurements are relatively insensitive to the specifics of the blast-wave model.

Solinger, A.↗

The dynamics and control of large flexible space structures-IV

The effects of solar radiation pressure as the main environmental disturbance torque were incorporated into the model of the rigid orbiting shallow shell and computer simulation results indicate that within the linear range the rigid modal amplitudes are excited in proportion to the area to mass ratio. The effect of higher order terms in the gravity-gradient torque expressions previously neglected was evaluated and found to be negligible for the size structures under consideration. A graph theory approach was employed for calculating the eigenvalues of a large flexible system by reducing the system (stiffness) matrix to lower ordered submatrices. The related reachability matrix and term rank concepts are used to verify controllability and can be more effective than the alternate numerical rank tests. Control laws were developed for the shape and orientation control of the orbiting flexible shallow shell and numerical results presented.

Bainum, P. M.↗

Melting and solidification of acoustically levitated drops

A single axis acoustic levitator has been used to suspend a drop of commercial grade succinonitrile in air. The drop is subjected to thermal cycling in order to observe melting and solidification of the drop in the acoustically levitated state. Melting takes place over a temperature rise of a few degrees, which is typical for an impure material. Solidification initially progresses along the drop surface because of the temperature gradient developed in the drop. The resultant solid drop shows a dendritic surface feature with a cavity in some circumstances. An attempt is made to calculate the growth velocity of the dendrites.

Ohsaka, K.↗

General contraction of Gaussian basis sets. Part 2: Atomic natural orbitals and the calculation of atomic and molecular properties

A recently proposed scheme for using natural orbitals from atomic configuration interaction (CI) wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outmost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital (ANO) sets.

Almloef, Jan↗

Downward continuation of the free-air gravity anomalies to the ellipsoid using the gradient solution and terrain correction: An attempt of global numerical computations

The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere-at least to the ellipsoid. The goal is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (the g sub 1 term). The terrain correction was also computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. The fast Fourier transformation was applied to the computations.

Wang, Y. M.↗

General contraction of Gaussian basis sets. II - Atomic natural orbitals and the calculation of atomic and molecular properties

A recently proposed scheme for using natural orbitals from atomic configuration interaction wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outermost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital sets.

Almlof, Jan↗

Microwave Imaging on Metal Objects

This final report for the project discusses the attempts to model, using different methods, microwave image reconstruction. Maximum Entropy Method was not successful. Attempts to use Singular Value Decomposition (SVD) got some good results after initial failure. SVD is based upon a theory of linear algebra, to the effect that any M X N Matrix A whose number of rows M is greater than or equal to its number of columns, N can be written as the product of an M X N column-orthogonal matrix U, an N X N diagonal Matrix, W, with m positive or zero elements (the singular values) and the transposition of an N X N orthogonal matrix V. In microwave imaging, the scattered fields can be expressed by the induced current distribution. The SVD method required more contiguous computer memory than was available. Work was also done on the Conjugate Gradient Method (CGM), which didn't work well when tried earlier. It was found that separation of the imaginary part and the real part during calculation may work. This work was considered incomplete as of the end of the grant period.

Tolliver, C. L.↗

A Numerical Study of Low-Reynolds-Number Separation Bubbles

The present study uses two dimensional numerical simulations to study unsteady low-Reynolds-number separation bubbles. The numerical study is in two parts: (1) a two dimensional time-accurate Navier-Stokes solver is used to simulate flows over the APEX airfoil, and (2) a numerical procedure is developed for localized simulations of transitional separation bubbles. The 2-D computations of flow over the APEX airfoil show that the flow is unsteady with periodic vortex shedding. A linear stability analysis of the separated flow shows that the vortex shedding is caused due to the instability of the separated flow. For transonic flows over the APEX airfoil the vortex shedding is additionally influenced by the presence of shocks. The flowfield has two characteristic time scales, one corresponding to the vortex shedding and another corresponding to the movement of the shocks. The two dimensional (2-D) airfoil simulations also showed the presence of nonlinear effects in the separated region. To better understand the characteristics of separation bubbles a numerical procedure has been developed for localized separation bubble calculations. This procedure is used to perform computations for a flat plate separation bubble test case. The separation bubble is induced by specifying a velocity gradient in the freestream. The growth of disturbances in the separation bubble is analyzed by introducing disturbances upstream of the separation bubble.

Tatineni, Mahidhar↗

Remote Sensing of the Urban Heat Island Effect Across Biomes in the Continental USA

Impervious surface area (ISA) from the Landsat TM-based NLCD 2001 dataset and land surface temperature (LST) from MODIS averaged over three annual cycles (2003-2005) are used in a spatial analysis to assess the urban heat island (UHI) skin temperature amplitude and its relationship to development intensity, size, and ecological setting for 38 of the most populous cities in the continental United States. Development intensity zones based on %ISA are defined for each urban area emanating outward from the urban core to the nonurban rural areas nearby and used to stratify sampling for land surface temperatures and NDVI. Sampling is further constrained by biome and elevation to insure objective intercomparisons between zones and between cities in different biomes permitting the definition of hierarchically ordered zones that are consistent across urban areas in different ecological setting and across scales. We find that ecological context significantly influences the amplitude of summer daytime UHI (urban-rural temperature difference) the largest (8 C average) observed for cities built in biomes dominated by temperate broadleaf and mixed forest. For all cities combined, ISA is the primary driver for increase in temperature explaining 70% of the total variance in LST. On a yearly average, urban areas are substantially warmer than the non-urban fringe by 2.9 C, except for urban areas in biomes with arid and semiarid climates. The average amplitude of the UHI is remarkably asymmetric with a 4.3 C temperature difference in summer and only 1.3 C in winter. In desert environments, the LST's response to ISA presents an uncharacteristic "U-shaped" horizontal gradient decreasing from the urban core to the outskirts of the city and then increasing again in the suburban to the rural zones. UHI's calculated for these cities point to a possible heat sink effect. These observational results show that the urban heat island amplitude both increases with city size and is seasonally asymmetric for a large number of cities across most biomes. The implications are that for urban areas developed within forested ecosystems the summertime UHI can be quite high relative to the wintertime UHI suggesting that the residential energy consumption required for summer cooling is likely to increase with urban growth within those biomes.

Imhoff, Marc L.↗

Temperature-aware task scheduling and proactive power management

Systems, apparatuses, and methods for performing temperature-aware task scheduling and proactive power management. A SoC includes a plurality of processing units and a task queue storing pending tasks. The SoC calculates a thermal metric for each pending task to predict an amount of heat the pending task will generate. The SoC also determines a thermal gradient for each processing unit to predict a rate at which the processing unit's temperature will change when executing a task. The SoC also monitors a thermal margin of how far each processing unit is from reaching its thermal limit. The SoC minimizes non-uniform heat generation on the SoC by scheduling pending tasks from the task queue to the processing units based on the thermal metrics for the pending tasks, the thermal gradients of each processing unit, and the thermal margin available on each processing unit.

97 MATHEMATICS AND COMPUTING↗

Determination of the Thicknesses and Visualization of Ion-Exchange Waveguides in Glasses by Scanning Electron Microscopy

The formation of a surface layer with silver and sodium concentration gradients in K8-grade glasses as a result of ion exchange is demonstrated by scanning electron microscopy (SEM) and electron probe X-ray microanalysis. The silver ions enriched layers with different thicknesses, depending on the duration of the ion exchange from the silver melt, are visualized in the secondary electrons mode. SEM data on the thickness of the silver enriched layer are in good agreement with the results of calculation from the diffusion equation for silver in glass. The dependence of the silver concentration on the diffusion depth is obtained using elemental composition mapping over the glass cross section. The presence of a layer with a silver concentration gradient leads to formation of a gradient waveguide. Discrete peaks in silver concentration profiles, caused by specific features of mutual diffusion of sodium and silver ions in glasses, are found.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗