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,549 records · Page 86

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Optimization of turbine design in wind farms with multiple hub heights, using exact analytic gradients and structural constraints

Wind farms are generally designed with turbines of all the same hub height. If wind farms were designed with turbines of different hub heights, wake interference between turbines could be reduced, lowering the cost of energy (COE). This paper demonstrates a method to optimize onshore wind farms with two different hub heights using exact, analytic gradients. Gradient-based optimization with exact gradients scales well with large problems and is preferable in this application over gradient-free methods. Our model consisted of the following: a version of the FLOw Redirection and Induction in Steady-State wake model that accommodated three-dimensional wakes and calculated annual energy production, a wind farm cost model, and a tower structural model, which provided constraints during optimization. Structural constraints were important to keep tower heights realistic and account for additional mass required from taller towers and higher wind speeds. We optimized several wind farms with tower height, diameter, and shell thickness as coupled design variables. Our results indicate that wind farms with small rotors, low wind shear, and closely spaced turbines can benefit from having two different hub heights. A nine-by-nine grid wind farm with 70-meter rotor diameters and a wind shear exponent of 0.08 realized a 4.9% reduction in COE by using two different tower sizes. If the turbine spacing was reduced to 3 diameters, the reduction in COE decreased further to 11.2%. Allowing for more than two different turbine heights is only slightly more beneficial than two heights and is likely not worth the added complexity.

17 WIND ENERGY↗

DOE OSTI2022

Spatially resolved properties of early-type group-dominant galaxies with MUSE: gas content, ionization mechanisms, and metallicity gradients

ABSTRACT With the goal of a thorough investigation of the ionized gas and its origin in early-type group-dominant galaxies, we present archival MUSE data for 18 galaxies from the Complete Local-Volume Groups Sample (CLoGS). These data allowed us to study the spatially resolved warm gas properties, including the morphology of the ionized gas, EW(H α), and kinematics as well as the gas-phase metallicity (12 + log(O/H)) of these systems. In order to distinguish between different ionization mechanisms, we used the emission-line ratios [O iii]/H β and [N ii]/H α in the BPT diagrams and EW(H α). We find that the ionization sources in our sample have variable impacts at different radii; central regions are more influenced by low-luminosity active galactic nuclei, while extended regions of low-ionization nuclear emission-line region-like emission are ionized by other mechanisms with post-asymptotic giant branch stars photoionization likely contributing significantly. We classified our sample into three H α +[N ii] emission morphology types. We calculate the gas-phase metallicity assuming several methods and ionization sources. In general, 12 + log(O/H) decreases with radius from the centre for all galaxies, independently of nebular morphology type, indicating a metallicity gradient in the abundance profiles. Interestingly, the more extended filamentary structures and all extranuclear star-forming regions present shallow metallicity gradients. Within the uncertainties these extended structures can be considered chemically homogeneous. We suggest that group-dominant galaxies in our sample likely acquired their cold gas in the past as a consequence of one or more mechanisms, e.g. gas-clouds or satellite mergers/accretion and/or cooling flows that contribute to the growth of the ionized gas structures.

Lagos, P. (ORCID:0000000223218657)↗

DOE OSTI2021

Time-dependent experimental identification of inter-ELM microtearing modes in the tokamak edge on DIII-D

In a series of discharges on the DIII-D tokamak, fast vertical plasma jogs are used to induce current perturbations in the steep gradient region of the H-mode edge. These current perturbations directly influence the edge q profile, decoupling the resonant location and instability drive of pedestal-localized microtearing modes (MTMs.) By exploiting this effect, we develop and apply a new experimental technique to track the dynamical frequency evolution of MTMs in the pedestal region, providing a compelling validation of the MTM model. The frequency of potential MTMs is calculated as the Doppler-shifted electron diamagnetic frequency at rational q=m/n surfaces, showing remarkable agreement with chirped frequency behavior of n=3, 4 and 5 modes detected with fast magnetics. Data is collected throughout multiple ELM cycles in order to build robust statistics describing the time-dependent frequency evolution of MTMs, which can be explained by examining the recovery of pedestal gradients after an ELM event. MTMs have a dominant transport contribution in the electron thermal channel, so the presented results indicate that reduced models of pedestal transport must be electromagnetic in nature and constructed with accurate calculations of MTM stability; inclusion of this physics is essential for accurate predictions of the electron temperature pedestal profile. Supporting measurements of mode saturation, propagation direction and transport fingerprints are made to support the dynamic frequency determination, unambiguously and experimentally identifying MTMs in the pedestal region of DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗

DOE OSTI2015

SU-D-303-03: Impact of Uncertainty in T1 Measurements On Quantification of Dynamic Contrast Enhanced MRI

Purpose: Quantification of dynamic contrast enhanced (DCE) MRI requires native longitudinal relaxation time (T1) measurement. This study aimed to assess uncertainty in T1 measurements using two different methods. Methods and Materials: Brain MRI scans were performed on a 3T scanner in 9 patients who had low grade/benign tumors and partial brain radiotherapy without chemotherapy at pre-RT, week-3 during RT (wk-3), end-RT, and 1, 6 and 18 months after RT. T1-weighted images were acquired using gradient echo sequences with 1) 2 different flip angles (50 and 150), and 2) 5 variable TRs (100–2000ms). After creating quantitative T1 maps, average T1 was calculated in regions of interest (ROI), which were distant from tumors and received a total of accumulated radiation doses < 5 Gy at wk-3. ROIs included left and right normal Putamen and Thalamus (gray matter: GM), and frontal and parietal white matter (WM). Since there were no significant or even a trend of T1 changes from pre-RT to wk-3 in these ROIs, a relative repeatability coefficient (RC) of T1 as a measure of uncertainty was estimated in each ROI using the data pre-RT and at wk-3. The individual T1 changes at later time points were evaluated compared to the estimated RCs. Results: The 2-flip angle method produced small RCs in GM (9.7–11.7%) but large RCs in WM (12.2–13.6%) compared to the saturation-recovery (SR) method (11.0–17.7% for GM and 7.5–11.2% for WM). More than 81% of individual T1 changes were within T1 uncertainty ranges defined by RCs. Conclusion: Our study suggests that the impact of T1 uncertainty on physiological parameters derived from DCE MRI is not negligible. A short scan with 2 flip angles is able to achieve repeatability of T1 estimates similar to a long scan with 5 different TRs, and is desirable to be integrated in the DCE protocol. Present study was supported by National Institute of Health (NIH) under grant numbers; UO1 CA183848 and RO1 NS064973.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2013

High Power Co-Axial Coupler

A superconducting RF (SRF) power coupler capable of handling 500 kW CW RF power at 750 MHz is required for present and future storage rings and linacs. There are over 35 coupler designs for SRF cavities ranging in frequency from 325 to 1500 MHz. Coupler windows vary from cylinders to cones to disks and RF power couplers will always be limited by the ability of ceramic windows and their matching systems to withstand the stresses due to non-uniform heating from dielectric and wall losses, multipactor, and mechanical flexure. In the Phase II project, we built a double window coaxial system with materials that would not otherwise be useable due to individual VSWRs. Double window systems can be operated such that one is cold (LN2) and one is warm. They can have different materials and still have a good match without using matching elements that create problematic multipactor bands. The match of the two windows will always result from the cancellation of the two window’s reflections when they are located approximately a quarter wavelength apart or multiples of a quarter wavelength. The window assemblies were carefully constructed to put the window material and its braze joint in compression at all times. This was done using explosion bonding techniques which allow for inexpensive fabrication of the vacuum / compression ring out of stainless steel with copper plating applied to the inner surface. The EIA 3-1/8” double window assembly was then successfully baked out and tested to 12 kW in a 3-1/8” co-axial system. The thermal gradient across the window was measured to be 90 C which represents about 15 ksi tensile stress in an uncompressed window. In our design the compression was calculated to be about 25 ksi, so the net compressive force was 5 ksi at full power.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2018

Simulation of Conduction Cooled SRF Cavity

SRF cavities are used heavily in high energy physics research, and simplifying operation of SRF cavities could enable their use in many fields of industry. Liquid helium coolant systems are a notable complexity of operating SRF cavities, and replacing liquid helium with cryocooler conduction cooling is an important research and development interest. As part of this effort, a system consisting of a 650 MHz single cell cavity, a commercial grade pulse tube refrigerator, and an aluminum thermal link was simulated through finite element analysis. The electromagnetic fields inside the cavity, heat flow through the system, and temperature distribution of the thermal link and cavity were calculated. Based on the simulations, the pulse tube refrigerator is capable of supplying adequate cooling power to a cavity with \ce{Nb3Sn} coating operating with an accelerating gradient between 6 and 12 MV/m. The maximum simulated temperature of the cavity was below 7 K, and the oscilla tions in temperature introduced by the pulse tube refrigerator were shown to dampen in amplitude before reaching the cavity surface.

43 PARTICLE ACCELERATORS↗

NASA NTRS1986

A search for the prewetting line

This paper describes efforts to locate the prewetting line in a binary liquid system (isopropanol-perfluoromethylcyclohexane) at the vapor-liquid interface. Tight upper bounds were placed on the temperature separation (0.2 K) between the prewetting line and the line of bulk liquid phase separation. The prewetting line in systems at equilibrium was not detected. Experimental signatures indicative of the prewetting line occurred only in nonequilibrium situations. Several theories predict that the adsorption of one of the components (the fluorocarbon, in this case) at the liquid-vapor interface should increase abruptly, at a temperature sightly above the temperature at which the mixture separates into two liquid phases. A regular solution calculation indicates that this prewetting line should have been easily detectable with the instruments used in this experiment. Significant features of the experiment are: (1) low-gradient thermostatting, (2) in situ stirring, (3) precision ellipsometry from the vapor-liquid interface, (4) high resolution differential index of refraction measurements using a novel cell design, and (5) computer control.

Schmidt, J. W.↗

NASA NTRS1986

Explicit finite-volume time-marching calculations of total temperature distributions in turbulent flow

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. This paper investigates the introduction of a new formulation of the energy equation which gives improved transient behavior as the calculation converges. The effect of variable Prandtl number on the temperature distribution through the boundary layer is also investigated. A turbulent boundary layer in an adverse pressure gradient (M = 0.55) is used to demonstrate the improved transient temperature distribution obtained when the new formulation of the energy equation is used. A flat plate turbulent boundary layer with a supersonic free-stream Mach number of 2.8 is used to investigate the effect of Prandtl number on the distribution of properties through the boundary layer. The computed total temperature distribution and recovery factor agree well with the measurements when a variable Prandtl number is used through the boundary layer.

Nicholson, Stephen↗

NASA NTRS1987

Explicit finite-volume time-marching calculations of total temperature distributions in turbulent flow

Nicholson, Stephen↗

NASA NTRS1992

Determination of refractive properties of fluids for dual-wavelength interferometry

Methods to calculate the refractive properties of solutions at different wavelengths are described by using experimental data at just two wavelengths. The properties are the refractive index and its gradients with temperature and concentration. Cauchy's equation is used to determine the refractive indices. The gradients versus temperature and concentration are then determined by using the Murphy-Alpert and the Lorentz-Lorenz equation, respectively. Finally, the particular case of a triglycine sulfate aqueous solution is considered as an example. The approach should provide the desired information for fringe analysis when dual-wavelength holographic or other interferometry is used for solving heat and mass transfer problems in fluids during crystal-growth experiments.

Vikram, Chandra S.↗

NASA NTRS2000

An Initial Study of the Sensitivity of Aircraft Vortex Spacing System (AVOSS) Spacing Sensitivity to Weather and Configuration Input Parameters

A study has been performed on a computer code modeling an aircraft wake vortex spacing system during final approach. This code represents an initial engineering model of a system to calculate reduced approach separation criteria needed to increase airport productivity. This report evaluates model sensitivity toward various weather conditions (crosswind, crosswind variance, turbulent kinetic energy, and thermal gradient), code configurations (approach corridor option, and wake demise definition), and post-processing techniques (rounding of provided spacing values, and controller time variance).

Riddick, Stephen E.↗

NASA NTRS2003

Dissociative Ionization of Aromatic and Heterocyclic Molecules

Space radiation poses a major health hazard to humans in space flight. The high-energy charged particles in space radiation ranging from protons to high atomic number, high-energy (HZE) particles, and the secondary species they produce, attack DNA, cells, and tissues. Of the potential hazards, long-term health effects such as carcinogenesis are likely linked to the DNA lesions caused by secondary electrons in the 1 - 30 eV range. Dissociative ionization (DI) is one of the electron collision processes that can damage the DNA, either directly by causing a DNA lesion, or indirectly by producing radicals and cations that attack the DNA. To understand this process, we have developed a theoretical model for DI. Our model makes use of the fact that electron motion is much faster than nuclear motion and assumes DI proceeds through a two-step process. The first step is electron-impact ionization resulting in a particular state of the molecular ion in the geometry of the neutral molecule. In the second step the ion undergoes unimolecular dissociation. Thus the DI cross section sigma(sup DI)(sub a) for channel a is given by sigma(sup DI)(sub a) = sigma(sup I)(sub a) P(sub D) with sigma(sup I)(sub a) the ionization cross section of channel a and P(sub D) the dissociation probability. This model has been applied to study the DI of H2O, NH3, and CH4, with results in good agreement with experiment. The ionization cross section sigma(sup I)(sub a) was calculated using the improved binary encounter-dipole model and the unimolecular dissociation probability P(sub D) obtained by following the minimum energy path determined by the gradients and Hessians of the electronic energy with respect to the nuclear coordinates of the ion. This model is used to study the DI from the low-lying channels of benzene and pyridine to understand the different product formation in aromatic and heterocyclic molecules. DI study of the DNA base thymine is underway. Solvent effects will also be discussed.

Huo, Winifred M.↗

NASA NTRS2012

Mathematical Simulation for Integrated Linear Fresnel Spectrometer Chip

A miniaturized solid-state optical spectrometer chip was designed with a linear gradient-gap Fresnel grating which was mounted perpendicularly to a sensor array surface and simulated for its performance and functionality. Unlike common spectrometers which are based on Fraunhoffer diffraction with a regular periodic line grating, the new linear gradient grating Fresnel spectrometer chip can be miniaturized to a much smaller form-factor into the Fresnel regime exceeding the limit of conventional spectrometers. This mathematical calculation shows that building a tiny motionless multi-pixel microspectrometer chip which is smaller than 1 cubic millimter of optical path volume is possible. The new Fresnel spectrometer chip is proportional to the energy scale (hc/lambda), while the conventional spectrometers are proportional to the wavelength scale (lambda). We report the theoretical optical working principle and new data collection algorithm of the new Fresnel spectrometer to build a compact integrated optical chip.

Park, Yeonjoon↗

NASA NTRS1983

Engineering calculations for communications satellite systems planning

Progress is reported on a computer code to improve the efficiency of spectrum and orbit utilization for the Broadcasting Satellite Service in the 12 GHz band for Region 2. It implements a constrained gradient search procedure using an exponential objective function based on aggregate signal to noise ratio and an extended line search in the gradient direction. The procedure is tested against a manually generated initial scenario and appears to work satisfactorily. In this test it was assumed that alternate channels use orthogonal polarizations at any one satellite location.

Martin, C. H.↗

DOE OSTI2013

Experimentally studied dynamic dose interplay does not meaningfully affect target dose in VMAT SBRT lung treatments

Purpose: The effects of respiratory motion on the tumor dose can be divided into the gradient and interplay effects. While the interplay effect is likely to average out over a large number of fractions, it may play a role in hypofractionated [stereotactic body radiation therapy (SBRT)] treatments. This subject has been extensively studied for intensity modulated radiation therapy but less so for volumetric modulated arc therapy (VMAT), particularly in application to hypofractionated regimens. Also, no experimental study has provided full four-dimensional (4D) dose reconstruction in this scenario. The authors demonstrate how a recently described motion perturbation method, with full 4D dose reconstruction, is applied to describe the gradient and interplay effects during VMAT lung SBRT treatments.Methods: VMAT dose delivered to a moving target in a patient can be reconstructed by applying perturbations to the treatment planning system-calculated static 3D dose. Ten SBRT patients treated with 6 MV VMAT beams in five fractions were selected. The target motion (motion kernel) was approximated by 3D rigid body translation, with the tumor centroids defined on the ten phases of the 4DCT. The motion was assumed to be periodic, with the period T being an average from the empirical 4DCT respiratory trace. The real observed tumor motion (total displacement ≤8 mm) was evaluated first. Then, the motion range was artificially increased to 2 or 3 cm. Finally, T was increased to 60 s. While not realistic, making T comparable to the delivery time elucidates if the interplay effect can be observed. For a single fraction, the authors quantified the interplay effect as the maximum difference in the target dosimetric indices, most importantly the near-minimum dose (D{sub 99%}), between all possible starting phases. For the three- and five-fractions, statistical simulations were performed when substantial interplay was found.Results: For the motion amplitudes and periods obtained from the 4DCT, the interplay effect is negligible (<0.2%). It is also small (0.9% average, 2.2% maximum) when the target excursion increased to 2–3 cm. Only with large motion and increased period (60 s) was a significant interplay effect observed, with D{sub 99%} ranging from 16% low to 17% high. The interplay effect was statistically significantly lower for the three- and five-fraction statistical simulations. Overall, the gradient effect dominates the clinical situation.Conclusions: A novel method was used to reconstruct the volumetric dose to a moving tumor during lung SBRT VMAT deliveries. With the studied planning and treatment technique for realistic motion periods, regardless of the amplitude, the interplay has nearly no impact on the near-minimum dose. The interplay effect was observed, for study purposes only, with the period comparable to the VMAT delivery time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2020

Thermochemically-Informed Mass Transport Model for Zr in U-Zr Fuel

Recent improvements to the coupled Thermochimica-MOOSE/BISON code system have enabled efficient calculations of species transport based on direct evaluation of composition and temperature dependent chemical potentials of the species. This presents an alternative to the traditional approach to species transport in nuclear fuels, which has been to employ a diffusion formulation that combines concentration-gradient driven Fickian diffusion with a Soret term based on a heat of transport fit to experimental data. Here we describe the application of the coupled code system to the diffusion of Zr in U-Zr metallic fuel. New classes implemented in BISON to solve this problem are documented. The Zr concentration profile after 50 years of diffusion is found to be strongly dependent on the assumptions made pertaining how to mobility is calculated in multi-phase regions of the fuel element. Two assumptions are compared (simple averaging of mobilities and using the majority phase mobility), and good qualitative agreement with experimental measurements is obtained using the majority phase assumption.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Selected results from full-core hydrogen redistribution asymptotic analysis in YH-moderated heat-pipe cooled microreactor

Yttrium hydride is the main candidate for moderation of high temperature nuclear microreactors . This is due to its high thermal stability and relatively high hydrogen retention at temperatures exceeding 870 C. One of the main issues associated with the use of yttrium hydride is that, when exposed to temperature, stress, or concentration gradients, the hydrogen contained in the metallic matrix tends to redistribute and leak from the moderating elements, potentially leading to reactivity losses and power swings. This paper is focused on presenting selected results from the asymptotic hydrogen redistribution analysis performed on the Simplified Microreactor Benchmark Assessment (SiMBA) problem, a microreactor core conceptual design developed at INL. The analysis is based on full core analysis, relies uniquely on Bison to perform hydrogen redistribution calculations, and discusses the physical causes of the reactivity feedback. It was found that the feedback is negative, but it is one order of magnitude lower than the one found for the empire reactor unit-cell. This is due to the lower axial temperature gradient together with the effect of the reflector.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗