Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient methods”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

CFD-based design optimization of ducted hydrokinetic turbines

Abstract Hydrokinetic turbines extract kinetic energy from moving water to generate renewable electricity, thus contributing to sustainable energy production and reducing reliance on fossil fuels. It has been hypothesized that a duct can accelerate and condition the fluid flow passing the turbine blades, improving the overall energy extraction efficiency. However, no substantial evidence has been provided so far for hydrokinetic turbines. To investigate this problem, we perform a CFD-based optimization study with a blade-resolved Reynolds-averaged Navier–Stokes (RANS) solver to explore the design of a ducted hydrokinetic turbine that maximizes the efficiency of energy extraction. A gradient-based optimization approach is utilized to effectively deal with the high-dimensional design space of the blade and duct geometry, with gradients being calculated through the adjoint method. The final design is re-evaluated through higher-fidelity unsteady RANS (URANS) simulations. Our optimized ducted turbine achieves an efficiency of about 54% over a range of operating conditions, higher than the typical 46% efficiency of unducted turbines.

13 HYDRO ENERGY↗
DOE OSTI2012

Equation of state of dense plasmas by ab initio simulations: Bridging the gap between quantum molecular dynamics and orbital-free molecular dynamics at high temperature

The applicability of quantum molecular dynamics to the calculation of the equation of state of a dense plasma is limited at high temperature by computational cost. Orbital-free molecular dynamics, based on the Thomas-Fermi semiclassical approximation and possibly on a gradient correction, is the only simulation method currently available at high temperature. We show in the case of a dense boron plasma that the two approaches give pressures differing by a few percent even at temperatures as high as a few tens of electron-volts. We indicate how the pressures obtained by orbital-free molecular dynamics can be corrected in order to appear as a limit of the quantum molecular dynamics results as temperature increases. We thus obtain a method to calculate the equation of state of a dense plasma up to high temperatures where quantum molecular dynamics cannot be directly implemented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Benchmarking gyrokinetic simulations in a toroidal flux-tube

A flux-tube model is implemented in the global turbulence code GEM [Y. Chen and S. E. Parker, J. Comput. Phys. 220, 839 (2007)] in order to facilitate benchmarking with Eulerian codes. The global GEM assumes the magnetic equilibrium to be completely given. The initial flux-tube implementation simply selects a radial location as the center of the flux-tube and a radial size of the flux-tube, sets all equilibrium quantities (B, ∇B, etc.) to be equal to the values at the center of the flux-tube, and retains only a linear radial profile of the safety factor needed for boundary conditions. This implementation shows disagreement with Eulerian codes in linear simulations. An alternative flux-tube model based on a complete local equilibrium solution of the Grad-Shafranov equation [J. Candy, Plasma Phys. Controlled Fusion 51, 105009 (2009)] is then implemented. This results in better agreement between Eulerian codes and the particle-in-cell (PIC) method. The PIC algorithm based on the v{sub ||}-formalism [J. Reynders, Ph.D. dissertation, Princeton University, 1992] and the gyrokinetic ion/fluid electron hybrid model with kinetic electron closure [Y. Chan and S. E. Parker, Phys. Plasmas 18, 055703 (2011)] are also implemented in the flux-tube geometry and compared with the direct method for both the ion temperature gradient driven modes and the kinetic ballooning modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Dirac cones in transition metal doped boron nitride

The transition metal (TM) doped zinc blende boron nitride (c-BN) is studied by using the first principle calculation. TM atoms fill in the interstitials in c-BN and form two-dimensional honeycomb lattice. The generalized gradient approximation and projector augmented wave method are used. The calculated density of states and band structures show that d electrons of TM atoms form impurity bands in the gap of c-BN. When the TM-BN system is in ferromagnetic or non-magnetic state, Dirac cones emerge at the K point in Brillouin zone. When TM is Ti and Co, the Dirac cones are spin polarized and very close to the Fermi level, which makes them promising candidates of Dirac half-metal [H. Ishizuka and Y. Motome, Phys. Rev. Lett. 109, 237207 (2012)]. While TM is Ni and Cu, the system is non-magnetic and Dirac cones located above the Fermi level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Accurate adiabatic singlet-triplet gaps in atoms and molecules employing the third-order spin-flip algebraic diagrammatic construction scheme for the polarization propagator

For the calculation of adiabatic singlet-triplet gaps (STG) in diradicaloid systems the spin-flip (SF) variant of the algebraic diagrammatic construction (ADC) scheme for the polarization propagator in third order perturbation theory (SF-ADC(3)) has been applied. Due to the methodology of the SF approach the singlet and triplet states are treated on an equal footing since they are part of the same determinant subspace. This leads to a systematically more accurate description of, e.g., diradicaloid systems than with the corresponding non-SF single-reference methods. Furthermore, using analytical excited state gradients at ADC(3) level, geometry optimizations of the singlet and triplet states were performed leading to a fully consistent description of the systems, leading to only small errors in the calculated STGs ranging between 0.6 and 2.4 kcal/mol with respect to experimental references.

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

Detection of topological materials with machine learning

Databases compiled using ab initio and symmetry-based calculations now contain tens of thousands of topological insulators and topological semimetals. This makes the application of modern machine learning methods to topological materials possible. Using gradient boosted trees, we show how to construct a machine learning model which can predict the topology of a given existent material with an accuracy of 90%. Such predictions are orders of magnitude faster than actual ab initio calculations. In this work, we use machine learning models to probe how different material properties affect topological features. Notably, we observe that topology is mostly determined by the “coarse-grained” chemical composition and crystal symmetry and depends little on the particular positions of atoms in the crystal lattice. We identify the sources of our model's errors and we discuss approaches to overcome them.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Prospects for the expansion of standing wave ambient pressure photoemission spectroscopy to reactions at elevated temperatures

Standing wave ambient pressure photoemission spectroscopy (SWAPPS) is a promising method to investigate chemical and potential gradients across solid-vapor and solid-liquid interfaces under close-to-realistic environmental conditions, far away from high vacuum. Until now, these investigations have been performed only near room temperature, but for a wide range of interfacial processes, chief among them being heterogeneous catalysis, measurements at elevated temperatures are required. One concern in these investigations is the temperature stability of the multilayer mirrors, which generate the standing wave field. At elevated temperatures, degradation of the multilayer mirror due to, for example, interdiffusion between the adjacent layers, decreases the modulation of the standing wave field, thus rendering SWAPPS experiments much harder to perform. Here, we show that multilayer mirrors consisting of alternate B4C and W layers are stable at temperatures exceeding 600 °C and are, thus, promising candidates for future studies of surface and subsurface species in heterogeneous catalytic reactions using SWAPPS.

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

Correction of MRI-induced geometric distortions in whole-body small animal PET-MRI

Purpose: The fusion of positron emission tomography (PET) and magnetic resonance imaging (MRI) data can be a challenging task in whole-body PET-MRI. The quality of the registration between these two modalities in large field-of-views (FOV) is often degraded by geometric distortions of the MRI data. The distortions at the edges of large FOVs mainly originate from MRI gradient nonlinearities. This work describes a method to measure and correct for these kind of geometric distortions in small animal MRI scanners to improve the registration accuracy of PET and MRI data. Methods: The authors have developed a geometric phantom which allows the measurement of geometric distortions in all spatial axes via control points. These control points are detected semiautomatically in both PET and MRI data with a subpixel accuracy. The spatial transformation between PET and MRI data is determined with these control points via 3D thin-plate splines (3D TPS). The transformation derived from the 3D TPS is finally applied to real MRI mouse data, which were acquired with the same scan parameters used in the phantom data acquisitions. Additionally, the influence of the phantom material on the homogeneity of the magnetic field is determined via field mapping. Results: The spatial shift according to the magnetic field homogeneity caused by the phantom material was determined to a mean of 0.1 mm. The results of the correction show that distortion with a maximum error of 4 mm could be reduced to less than 1 mm with the proposed correction method. Furthermore, the control point-based registration of PET and MRI data showed improved congruence after correction. Conclusions: The developed phantom has been shown to have no considerable negative effect on the homogeneity of the magnetic field. The proposed method yields an appropriate correction of the measured MRI distortion and is able to improve the PET and MRI registration. Furthermore, the method is applicable to whole-body small animal imaging routines including different standard MRI sequences.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-163: Characterization of a Novel High Resolution 1D Silicon Monolithic Array for Small Field Commissioning and Quality Assurance

Purpose: To study the suitability of a novel 1D silicon monolithic array for dosimetry of small radiation fields and for QA of high dose gradient treatment modalities (IMRT and SBRT). Methods: A 1D array composed of 4 monolithic silicon modules of 64 mm length and 1 mm pixel pitch was developed by IBA Dosimetry. Measurements were carried out for 6MV and 15MV photons on two commercial different linacs (TrueBeam and Clinac iX, Varian Medical Systems, Palo Alto, CA) and for a CyberKnife G4 (Accuray Inc., Sunnyvale, CA). The 1D array was used to measure output factors (OF), profiles and off axis correction factors (OACF) for the Iris CyberKnife variable collimator (5–60 mm). In addition, dose profiles (at the isocenter plane) were measured for multiple IMRT and SBRT treatment plans and compared with those obtained using EDR2radiographic film (Carestream Health, Rochester NY), a commercial 2D diode array and with the dose distribution calculated using a commercial TPS (Eclipse, Varian Medical Systems, Palo Alto, CA). Results: Due to the small pixel pitch of the detector, IMRT and SBRT plan profiles deviate from film measurements by less than 2%. Similarly, the 1D array exhibits better performance than the 2D diode array due to the larger (7 mm) pitch of that device. Iris collimator OFs measured using the 1D silicon array are in good agreement with the commissioning values obtained using a commercial stereotactic diode as well as with published data. Maximum deviations are < 3% for the smallest field (5 and 7.5mm) and below 1% for all other dimensions. Conclusion: We have demonstrated good performances of the array for commissioning of small photon fields and in patient QA, compared with diodes and film typically used in these clinical applications. The technology compares favorably with existing commercial solutions The presenting author is founded by a Marie Curie Early Initial Training Network Fellowship of the European Communitys Seventh Framework Programme under contract number (PITN-GA-2011-289198-ARDENT). The research activity is hosted by IBA Dosimetry, Gmbh.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-IeP1-08: MR Geometric Distortion Dependency On Imaging Sequence, Acquisition Orientation and Receiver Bandwidth of a Dedicated 1.5T MR-Simulator

Purpose: To investigate the 3D geometric distortion of four potential MR sequences for radiotheraptic applications, and its dependency on sequence-type, acquisition-orientation and receiver-bandwidth from a dedicated 1.5T 700mm-wide bore MR-simulator (Magnetom-Aera, Sienmens Healthcare, Erlangen, Germany), using a large customized geometric accuracy phantom. Methods: This work studied 3D gradient-echo (VIBE) and spin-echo (SPACE) sequences for anatomical imaging; a specific ultra-short-TE sequence (PETRA) potentially for bone imaging and MR-based dosimetry; and a motion-insensitive sequence (BLADE) for dynamic applications like 4D-MRI. Integrated geometric-correction was employed, three orthogonal acquisition-orientations and up to three receiver-bandwidths were used, yielding 27 acquisitions for testing (Table 1a).A customized geometric accuracy phantom (polyurethane, MR/CT invisible, W×L×H:55×55×32.5cm3) was constructed and filled with 3892 spherical markers (6mm diameter, MR/CT visible) arranged on a 25mm-interval 3D isotropic-grid (Fig.1). The marker positions in MR images were quantitatively calculated and compared against those in the CT-reference using customized MatLab scripts. Results: The average distortion within various diameter-of-spherical-volumes (DSVs) and the usable DSVs under various distortion limits were measured (Tables 1b-c). It was observed that distortions fluctuated when sequence-type, acquisition-orientation or receiver-bandwidth changed (e.g. within 300mm-DSV, the lowest/highest average distortions of VIBE were 0.40mm/0.59mm, a 47.5% difference). According to AAPM-TG66 (<1mm distortion, left-most column of Table 1c), PETRA (Largest-DSV:253.9mm) has the potential on brain treatment, while BLADE (Largest-DSV:207.2mm) may need improvement for thoracic/abdominal applications. The results of VIBE (Largest-DSVs:294.3mm, the best among tested acquisitions) and SPACE (Largest-DSVs:267.7mm) suggests their potentials on head and neck applications. These Largest-DSVs were attained on different acquisition-orientations and receiver-bandwidths. Conclusion: Geometric distortion was shown to be dependent on sequence-type, acquisition-orientation and receiver-bandwidth. In the experiment, no configuration in any one of these factors could consistently reduce distortion while the others were varying. The distortion analysis result is a valuable guideline for sequence selection and optimization for MR-aided radiotherapy applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

STEMDL

The code "STEMDL" provides highly-optimized multi-gpu implementations of DL models of applicability in electron microscopy via the use of mixed-precision arithmetic, native graph construction and evaluation in Tensorflow, and adaptive scaling of learning parameters to avoid numerical instabilities due to reduced precision. STEMDL also uses low-level methods from the Horovod library for gradient communication and reduction and adds an additional optimization by allowing the user to scheduling the gradient reduction operations based on their position within the layers of a neural net. STEMDL provides optimized I/O routines for data staging and data reading/writing in both Tensorflow and Pytorch using multi-processing and asynchronous transfer between host and device. .

Laanait, Nouamane↗
DOE OSTI2024

Stochastic minibatch approach to the ptychographic iterative engine

The ptychographic iterative engine (PIE) is a widely used algorithm that enables phase retrieval at nanometer-scale resolution over a wide range of imaging experiment configurations. By analyzing diffraction intensities from multiple scanning locations where a probing wavefield interacts with a sample, the algorithm solves a difficult optimization problem with constraints derived from the experimental geometry as well as sample properties. The effectiveness at which this optimization problem is solved is highly dependent on the ordering in which we use the measured diffraction intensities in the algorithm, and random ordering is widely used due to the limited ability to escape from stagnation in poor-quality local solutions. In this study, we introduce an extension to the PIE algorithm that uses ideas popularized in recent machine learning training methods, in this case minibatch stochastic gradient descent. Our results demonstrate that these new techniques significantly improve the convergence properties of the PIE numerical optimization problem.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Scalable Second Order Optimization for Machine Learning

Many machine learning (ML) training tasks are essentially optimization processes that would at first glance appear eminently parallelizable and scalable. However, effective acceleration of these tasks with scalable parallel hardware has proven to be elusive. While standard methods for machine learning, e.g., stochastic gradient descent (SGD) for DNNs, tend to be resource efficient, they appear to be fundamentally sequential in nature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2010

ISOTHERMAL AIR INGRESS VALIDATION EXPERIMENTS AT IDAHO NATIONAL LABORATORY: DESCRIPTION AND SUMMARY OF DATA

Idaho National Laboratory performed air ingress experiments as part of validating computational fluid dynamics code (CFD). An isothermal stratified flow experiment was designed and set to understand stratified flow phenomena in the very high temperature gas cooled reactor (VHTR) and to provide experimental data for validating computer codes. The isothermal experiment focused on three flow characteristics unique in the VHTR air-ingress accident: stratified flow in the horizontal pipe, stratified flow expansion at the pipe and vessel junction, and stratified flow around supporting structures. Brine and sucrose were used as heavy fluids and water was used as light fluids. The density ratios were changed between 0.87 and 0.98. This experiment clearly showed that a stratified flow between heavy and light fluids is generated even for very small density differences. The code was validated by conducting blind CFD simulations and comparing the results to the experimental data. A grid sensitivity study was also performed based on the Richardson extrapolation and the grid convergence index method for modeling confidence. As a result, the calculated current speed showed very good agreement with the experimental data, indicating that the current CFD methods are suitable for predicting density gradient stratified flow phenomena in the air-ingress accident.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Prairie Strips and Lower Land Use Intensity Increase Biodiversity and Ecosystem Services

Agricultural landscapes can be managed to protect biodiversity and maintain ecosystem services. One approach to achieve this is to restore native perennial vegetation within croplands. Where rowcrops have displaced prairie, as in the US Midwest, restoration of native perennial vegetation can align with crops in so called “prairie strips.” We tested the effect of prairie strips in addition to other management practices on a variety of taxa and on a suite of ecosystem services. To do so, we worked within a 33-year-old experiment that included treatments that varied methods of agricultural management across a gradient of land use intensity. In the two lowest intensity crop management treatments, we introduced prairie strips that occupied 5% of crop area. We addressed three questions: (1) What are the effects of newly established prairie strips on the spillover of biodiversity and ecosystem services into cropland? (2) How does time since prairie strip establishment affect biodiversity and ecosystem services? (3) What are the tradeoffs and synergies among biodiversity conservation, non-provisioning ecosystem services, and provisioning ecosystem services (crop yield) across a land use intensity gradient (which includes prairie strips)? Within prairie strip treatments, where sampling effort occurred within and at increasing distance from strips, dung beetle abundance, spider abundance and richness, active carbon, decomposition, and pollination decreased with distance from prairie strips, and this effect increased between the first and second year. Across the entire land use intensity gradient, treatments with prairie strips and reduced chemical inputs had higher butterfly abundance, spider abundance, and pollination services. In addition, soil organic carbon, butterfly richness, and spider richness increased with a decrease in land use intensity. Crop yield in one treatment with prairie strips was equal to that of the highest intensity management, even while including the area taken out of production. We found no effects of strips on ant biodiversity and greenhouse gas emissions (N 2 O and CH 4 ). Our results show that, even in early establishment, prairie strips and lower land use intensity can contribute to the conservation of biodiversity and ecosystem services without a disproportionate loss of crop yield.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

kemmerl4/LTER_prairie_strip_2019-2020: Data for: Prairie strips and lower land use intensity increase biodiversity and ecosystem services

Agricultural landscapes can be managed to protect biodiversity and maintain ecosystem services. One approach to achieve this is to restore native perennial vegetation within croplands. Where rowcrops have displaced prairie, as in the US Midwest, restoration of native perennial vegetation can align with crops in so called “prairie strips.” We tested the effect of prairie strips in addition to other management practices on a variety of taxa and on a suite of ecosystem services. To do so, we worked within a 33-year-old experiment that included treatments that varied methods of agricultural management across a gradient of land use intensity. In the two lowest intensity crop management treatments, we introduced prairie strips that occupied 5% of crop area. We addressed three questions: (1) What are the effects of newly established prairie strips on the spillover of biodiversity and ecosystem services into cropland? (2) How does time since prairie strip establishment affect biodiversity and ecosystem services? (3) What are the tradeoffs and synergies among biodiversity conservation, non-provisioning ecosystem services, and provisioning ecosystem services (crop yield) across a land use intensity gradient (which includes prairie strips)? Within prairie strip treatments, where sampling effort occurred within and at increasing distance from strips, dung beetle abundance, spider abundance and richness, active carbon, decomposition, and pollination decreased with distance from prairie strips, and this effect increased between the first and second year. Across the entire land use intensity gradient, treatments with prairie strips and reduced chemical inputs had higher butterfly abundance, spider abundance, and pollination services. In addition, soil organic carbon, butterfly richness, and spider richness increased with a decrease in land use intensity. Crop yield in one treatment with prairie strips was equal to that of the highest intensity management, even while including the area taken out of production. We found no effects of strips on ant biodiversity and greenhouse gas emissions (N 2 O and CH 4 ). Our results show that, even in early establishment, prairie strips and lower land use intensity can contribute to the conservation of biodiversity and ecosystem services without a disproportionate loss of crop yield.

54 ENVIRONMENTAL SCIENCES↗