Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time parallelization”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Solving Unit Commitment Problems with Demand Responsive Loads: Preprint

This paper focuses on using variations of the Frank-Wolfe algorithm for solving unit commitment problems with high volumes of demand responsive loads on the power grid. We present a formulation of the unit commitment problem with demand responsive loads. We then show through reformulation and relaxations of the problem that variations of the Frank-Wolfe algorithm can be used to determine the time series decisions for the demand responsive loads. We show through computational experiments on the RTS-GMLC test system that the time series of demand responsive load decisions obtained through our approach are near optimal and provide details regarding how large scale parallel implementations of our approach can be highly computationally efficient.

demand response↗
DOE OSTI2019

Characterization of NbTiN Thin Film Structures

Approaching the bulk Nb material RF performance limits has urged development of alternative superconducting materials for superconducting radio frequency (SRF) accelerator cavities for further performance enhancement. A promising theory has predicted that thin film structures of superconductor-insulator-superconductor (SIS) [1] can delay magnetic flux penetration in accelerator cavities to higher fields. A candidate superconductor for the SIS structures is NbTiN. A few key aspects of SIS structures development are high quality individual layers, sharp interfaces and optimum thickness for first flux penetration (Hfp) delay. High quality monocrystalline NbTiN films are deposited by reactive DC magnetron sputtering. In a parallel development, interface quality was assessed by depositing bilayers of 3 nm NbTiN with ~1 nm AlN repeated up to 16 times with no increase in roughness of the structure. The stacked layers form a metamaterial, which could exhibit Tc greater than bulk NbTiN [2]. This contribution presents the characterization of the surface, material and superconductivity of NbTiN with concentration on the Hfp enhancement for 200 to 5 nm films and multilayer nanostructures.

Beverstock, David↗
DOE OSTI2015

Multipurpose epithermal neutron beam on new research station at MARIA research reactor in Swierk-Poland

MARIA reactor is an open-pool research reactor what gives the chance to install uranium fission converter on the periphery of the core. It could be installed far enough not to induce reactivity of the core but close enough to produce high flux of fast neutrons. Special design of the converter is now under construction. It is planned to set the research stand based on such uranium converter in the near future: in 2015 MARIA reactor infrastructure should be ready (preparation started in 2013), in 2016 the neutron beam starts and in 2017 opening the stand for material and biological research or for medical training concerning BNCT. Unused for many years, horizontal channel number H2 at MARIA research rector in Poland, is going to be prepared as a part of unique stand. The characteristics of the neutron beam will be significant advantage of the facility. High flux of neutrons at the level of 2x10{sup 9} cm{sup -2}s{sup -1} will be obtainable by uranium neutron converter located 90 cm far from the reactor core fuel elements (still inside reactor core basket between so called core reflectors). Due to reaction of core neutrons with converter U{sub 3}Si{sub 2} material it will produce high flux of fast neutrons. After conversion neutrons will be collimated and moderated in the channel by special set of filters and moderators. At the end of H2 channel i.e. at the entrance to the research room neutron energy will be in the epithermal energy range with neutron intensity at least at the level required for BNCT (2x10{sup 9} cm{sup -2}s{sup -1}). For other purposes density of the neutron flux could be smaller. The possibility to change type and amount of installed filters/moderators which enables getting different properties of the beam (neutron energy spectrum, neutron-gamma ratio and beam profile and shape) is taken into account. H2 channel is located in separate room which is adjacent to two other empty rooms under the preparation for research laboratories (200 m2). It is planned to create fully equipped complex facility possible to perform various experiments on the intensive neutron beam. Epithermal neutron beam enables development across the full spectrum of materials research for example shielding concrete tests or electronic devices construction improvement. Due to recent reports on the construction of the accelerator for the Boron Neutron Capture Therapy (BNCT) it has the opportunity to become useful and successful method in the fight against brain and other types of cancers not treated with well known medical methods. In Europe there is no such epithermal neutron source which could be used throughout the year for training and research for scientist working on BNCT what makes the stand unique in Europe. Also our research group which specializes in mixed radiation dosimetry around nuclear and medical facilities would be able to carry out research on new detectors and methods of measurements for radiological protection and in-beam (therapeutic) dosimetry. Another group of scientists from National Centre for Nuclear Research, where MARIA research reactor is located, is involved in research of gamma detector systems. There is an idea to develop Prompt-gamma Single Photon Emission Computed Tomography (Pg- SPECT). This method could be used as imaging system for compounds emitting gamma rays after nuclear reaction with thermal neutrons e.g. for boron concentration in BNCT. Inside the room, where H2 channel is located, there is another horizontal channel - H1 which is also unused. Simultaneously with the construction of the H2 stand it will be possible to create special pneumatic horizontal mail inside the H1 channel for irradiation material samples in the vicinity of the core i.e. in the distal part of the H1 channel. It might expand the scope of research at the planned neutron station. Secondly it is planned to equip both stands with moveable positioning system, video system and facilities to perform animal experiments (anaesthesia, vital signs control, imaging devices, positioning). These all above make constructed station unique in the world (uranium fission converter-based beam) and the only one of such intense neutron beam in the Europe. Moreover implementation of the station would allow the development of research on a number of issues for researchers from all over the Europe. One of very important advantages of the station is undisturbed exploitation of the reactor and other vertical and horizontal channels. MARIA reactor operates 6000 hours per year and that amount of time will be achievable for research on the neutron station. It have to be underlined that new neutron station will work parallel to all another ventures. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Development of ARCHER towards clinical use: modeling and simulation of Varian Linac for radiation therapy dose calculations - 25449

Monte Carlo particle transport methods are well recognized as the most accurate approach to dose calculations for radiation therapy treatment planning. However, the lengthy calculation time has prevented Monte Carlo methods from being applied to dose calculations in the planning of radiation therapy. Recently, the emerging parallel computing hardware accelerators, such as general-purpose Graphics Processing Unit (GPU), have significantly boosted the performance of Monte Carlo codes. The accuracy of a clinical dose engine depends greatly on the accuracy of source term modeling and the methodology of dose calculations. Hence, in this study, an efficient and accurate source model of the Varian TrueBeam{sup R} Linac is developed and integrated with a fast photon-electron transport module of the recently developed high-performance Monte Carlo code, ARCHER. The validity of the source model and the particle transport physics are verified with established dose calculation codes such as EGSnrc and validated with experimental measurements. To demonstrate the capability of the platform in clinical applications, three clinical treatment plans (the lung, breast, and prostate) are simulated and analyzed. From the perspective of performance, the entire simulation of each clinical plan, starting from the source modeling to the end of particle transport within the patient phantom, can be finished in less than one minute using a single Titan X GPU card with 1% statistical error. In comparison, the same cases take the EGSnrc code hundreds of CPU hours to finish. The results suggest that this new version of ARCHER code has been successfully developed and tested as a clinically deployable dose engine for external photon beam IMRT plans in the clinical setting. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

Short-Term Load Forecasting-Based Automatic Distribution Network Reconfiguration

In a traditional dynamic network reconfiguration study, the optimal topology is determined at every scheduled time point by using the real load data measured at that time. The development of the load forecasting technique can provide an accurate prediction of the load power that will happen in a future time and provide more information about load changes. With the inclusion of load forecasting, the optimal topology can be determined based on the predicted load conditions during a longer time period instead of using a snapshot of the load at the time when the reconfiguration happens; thus, the distribution system operator can use this information to better operate the system reconfiguration and achieve optimal solutions. This paper proposes a short-term load forecasting approach to automatically reconfigure distribution systems in a dynamic and pre-event manner. Specifically, a short-term and high-resolution distribution system load forecasting approach is proposed with a forecaster based on support vector regression and parallel parameters optimization. The network reconfiguration problem is solved by using the forecasted load continuously to determine the optimal network topology with the minimum amount of loss at the future time. The simulation results validate and evaluate the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Equilibrium Properties of the Mixed State in Superconducting Niobium in a Transverse Magnetic Field: Experiment and Theoretical Model

Equilibrium magnetic properties of the mixed state in type II superconductors were studied on high-purity film and single-crystal niobium samples with different Ginzburg-Landau parameters in perpendicular and parallel magnetic fields using dc magnetometry and scanning Hall-probe microscopy. The magnetization curve for samples with unity demagnetizing factor (slabs in perpendicular field) was obtained for the first time. It was found that none of the existing theories is consistent with these new data. To address this problem, a theoretical model is developed and comprehensively validated. The new model describes the mixed state in an averaged limit, i.e., without detailing the samples’ magnetic structure and therefore ignoring the surface current and interactions between the structural units (vortices). At low values of the Ginzburg-Landau parameter, it converts to the model of Peierls and London for the intermediate state in type I superconductors. The model quantitatively describes the magnetization curve for the perpendicular field and provides new insights into the properties of the mixed state, including properties of individual vortices. In particular, it suggests that description of the vortex matter in superconductors of the transverse geometry as a “gas-like” system of non-interacting vortices is more appropriate than the frequently used solid-like scenarios.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

High-Performance Molecular Dynamics Simulation for Biological and Materials Sciences: Challenges of Performance Portability, Dallas, TX, November 2018.

Highly-optimized parallel molecular dynamics programs have allowed researchers to achieve ground-breaking results in biological and materials sciences. This type of performance has come at the expense of portability: a significant effort is required for performance optimization on each new architecture. Using a metric that emphasizes speedup, we assess key accelerating programming components of four different best-performing molecular dynamics programs - GROMACS, NAMD, LAMMPS and CP2K-- each having a particular scope of application, for contribution to performance and for portability. We use builds with and without these components, tested on HPC systems. We also analyze the code-bases to determine compliance with portability recommendations. We find that for all four programs, the contributions of the non-portable components to speed are essential to the programs' performances; without them we see a reduction in time-to-solution of a magnitude that is insufferable to domain scientists. This characterizes the performance efficiency that must be approached for good performance portability on a programmatic level, suggesting solutions to this difficult problem, which should come from developers, industry and funding institutions, and possibly new research in programming languages.

Computer Science↗
DOE OSTI2018

Scalable Multilevel Uncertainty Quantification Concepts for Extreme-Scale Multiscale Problems

Our work in this project is aimed at making fundamental advances in multilevel and multiscale methods for problems with uncertainties related to flows in porous media. The main thrust of this research is to develop a systematic multiscale, multilevel parallel approaches for problems with uncertainties and use them in the uncertainty quantification. A key emphasis is on problems without an apparent scale separation in space, time, and uncertainty spaces. Multiscale solution methods are currently under active investigation. However, the develop- ment of multiscale methods for problems with uncertainties where the scales and uncertainties are tightly coupled is a research area that is less explored. During this project, we have made some significant contributions to the developments of multiscale methods, the uncertainty quantification for multiscale methods, multiscale model learning, uncertainty quantification in inverse problems, and applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Ion Manipulation from Liquid Xe to Vacuum: Ba-Tagging for a nEXO Upgrade and Future 0 νββ Experiments

Neutrinoless double beta decay (0𝜈𝛽𝛽) provides a way to probe physics beyond the Standard Model of particle physics. The upcoming nEXO experiment will search for 0𝜈𝛽𝛽 decay in 136 Xe with a projected half-life sensitivity exceeding 10 28 years at the 90% confidence level using a liquid xenon (LXe) Time Projection Chamber (TPC) filled with 5 tonnes of Xe enriched to ∼90% in the 𝛽𝛽-decaying isotope 136 Xe. In parallel, a potential future upgrade to nEXO is being investigated with the aim to further suppress radioactive backgrounds and to confirm 𝛽𝛽-decay events. This technique, known as Ba-tagging, comprises extracting and identifying the 𝛽𝛽-decay daughter 136 Ba ion. One tagging approach being pursued involves extracting a small volume of LXe in the vicinity of a potential 𝛽𝛽-decay using a capillary tube and facilitating a liquid-to-gas phase transition by heating the capillary exit. The Ba ion is then separated from the accompanying Xe gas using a radio-frequency (RF) carpet and RF funnel, conclusively identifying the ion as 136 Ba via laser-fluorescence spectroscopy and mass spectrometry. Simultaneously, an accelerator-driven Ba ion source is being developed to validate and optimize this technique. The motivation for the project, the development of the different aspects, along with the current status and results, are discussed here.

a-tagging↗
DOE OSTI2018

Tomosaic : efficient acquisition and reconstruction of teravoxel tomography data using limited-size synchrotron X-ray beams

X-rays offer high penetration with the potential for tomography of centimetre-sized specimens, but synchrotron beamlines often provide illumination that is only millimetres wide. Here an approach is demonstrated termed Tomosaic for tomographic imaging of large samples that extend beyond the illumination field of view of an X-ray imaging system. This includes software modules for image stitching and calibration, while making use of existing modules available in other packages for alignment and reconstruction. The approach is compatible with conventional beamline hardware, while providing a dose-efficient method of data acquisition. By using parallelization on a distributed computing system, it provides a solution for handling teravoxel-sized or larger datasets that cannot be processed on a single workstation in a reasonable time. Using experimental data, the package is shown to provide good quality three-dimensional reconstruction for centimetre-sized samples with sub-micrometre pixel size.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Characterization of a Pixelated Cadmium Telluride Detector System Using a Polychromatic X-Ray Source and Gold Nanoparticle-Loaded Phantoms for Benchtop X-Ray Fluorescence Imaging

In this paper, the imaging dose and scan time have been considered as the two major constraints for routine benchtop x-ray fluorescence computed tomography (XFCT) imaging. One way to address this issue is to acquire x-ray fluorescence (XRF) signals in parallel through a 2D array of single-crystal detectors or a pixelated detector along with the cone-beam x-ray source. To identify a detector system suitable for this purpose, a commercially available, fully spectroscopic cadmium telluride (CdTe) pixelated detector, HEXITEC (High-Energy X-ray Imaging Technology), was tested under the experimental conditions optimized for benchtop XFCT imaging of gold nanoparticles (GNPs). Specifically, two different parallel-hole stainless steel collimators were fabricated and coupled with the detector for seamless integration into our existing benchtop cone-beam XFCT system. After the detector deployment, this benchtop XFCT system was used to detect XRF photons from GNP-loaded phantoms. A pixel-merging algorithm was introduced to enhance the sensitivity of XRF photon detection thereby minimizing the scan time. The effect of pixel-level charge sharing correction algorithms was investigated within the context of benchtop XFCT imaging. The detector energy resolution, in terms of the full width at half maximum (FWHM) values at different gold K-shell XRF energies, was also determined. Of the two charge sharing correction algorithms examined, the charge sharing addition gave better sensitivity than the charge sharing discrimination (csd). On the other hand, under the current experimental conditions, the energy resolution of the HEXITEC detector was the best with the csd and estimated to be 1.56 keV FWHM at 66-69 keV photon energy. Overall, despite some degradation of the detector energy resolution (compared with typical single crystal CdTe detectors), the HEXITEC detector enabled parallel data acquisition under the experimental conditions typical of benchtop XFCT imaging and operated well within our benchtop XFCT setup.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Extreme-scale EV charging infrastructure planning for last-mile delivery using high-performance parallel computing

Here, this paper addresses stochastic charger location and allocation problems under queue congestion for last-mile delivery using electric vehicles (EVs). The objective is to decide where to open charging stations and how many chargers of each type to install, subject to budgetary and waiting-time constraints. We formulate the problem as a mixed-integer non-linear program, where each station-charger pair is modeled as a multiserver queue with stochastic arrivals and service times to capture the notion of waiting in fleet operations. The model is extremely large, with billions of variables and constraints for a typical metropolitan area; even loading the model in solver memory is difficult, let alone solving it. To address this challenge, we develop a Lagrangian-based dual decomposition framework that decomposes the problem by station and leverages parallelization on high-performance computing systems, where the subproblems are solved by using a cutting plane method and their solutions are collected at the master level. We also develop a three-step rounding heuristic to transform the fractional subproblem solutions into feasible integral solutions. Computational experiments on data from the Chicago metropolitan area with hundreds of thousands of households and thousands of candidate stations show that our approach produces high-quality solutions in cases where existing exact methods cannot even load the model in memory. We also analyze various policy scenarios, demonstrating that combining existing depots with newly built stations under multiagency collaboration substantially reduces costs and congestion. These findings offer a scalable and efficient framework for developing sustainable large-scale EV charging networks.

Capacity allocation↗
DOE OSTI2024

Experimental investigation of PDI bifurcation of lower hybrid waves during electron density ramp-up in EAST

Abstract The effect of parametric decay instability (PDI) on the current drive efficiency of 4.6 GHz lower hybrid (LH) waves in EAST is investigated experimentally, showing the PDI channel bifurcation of LH waves for the first time in EAST. First, experiments with three platforms of LH power were performed, achieving the LH power required for the PDI occurrence. Second, PDI bifurcation experiments were further carried out by ramping up the plasma electron density. The loop voltage increases with an increase in density, implying a decrease in the current driven by the LH wave. PDI bifurcation during electron density ramp-up was studied by analyzing the parallel refractive index ( n ∥ ) and the frequency spectrum broadening, which is measured with a radio frequency magnetic probe array recently installed close to the LH antenna. It is observed for the first time that they both first increase with density, then there is not much variation and a clear sideband in the frequency spectrum is also observed when the density is up to 4 × 10 19 m − 3 , suggesting a change in the PDI channel. Calculation of the mode growth rate driven by PDI shows that when the edge electron density is up to 1.9 × 10 18 m − 3 , the growth rate of the ion cyclotron quasi-mode (ICQM) will exceed that of the ion sound quasi-mode (ISQM), quantitively explaining that with an increase in density, the PDI channel partly transits from the ISQM to the ICQM channel. Studies provide a possible way to reduce the power deposition in the edge region and improve drive capability by means of mitigating PDI behavior.

Physics↗
DOE OSTI2023

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

On the mobilization of metals by CO 2 leakage into shallow aquifers: exploring release mechanisms by modeling field and laboratory experiments

The dissolution of CO 2 in water leads to a pH decrease and a carbonate content increase in affected groundwater, which in turn can drive the mobilization of metals from sediments. The mechanisms of metal release postulated in various field and laboratory studies often differ. Drawing primarily on previously published results, we examine contrasting metal mobilization behaviors at two field tests and in one laboratory study, to investigate whether the same mechanisms could explain metal releases in these different experiments. Numerical modeling of the two field tests reveals that fast Ca-driven cation exchange (from calcite dissolution) can explain the release of most major and trace metal cations at both sites, and their parallel concentration trends. The dissolution of other minerals reacting more slowly (superimposed on cation exchange) also contributes to metal release over longer time frames, but can be masked by fast ambient groundwater velocities. Therefore, the magnitude and extent of mobilization depends not only on metal-mineral associations and sediment pH buffering characteristics, but also on groundwater flow rates, thus on the residence time of CO 2 -impacted groundwater relative to the rates of metal-release reactions. Sequential leaching laboratory tests modeled using the same metal-release concept as postulated from field experiments show that both field and laboratory data can be explained by the same processes. Finally, the reversibility of metal release upon CO 2 degassing by de-pressurization is also explored using simple geochemical models, and shows that the sequestration of metals by resorption and re-precipitation upon CO 2 exsolution is quite plausible and may warrant further attention.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kinetic Monte Carlo simulations of structural evolution during anneal of additively manufactured materials

Our experiments indicated that upon a post-processing anneal, an additively manufactured 316L stainless steel exhibits cubic grains rather than the conventional equiaxed grains. In this work, we have used kinetic Monte Carlo simulations to explore the origin of these cubic grains. First, we implemented a new kinetic Monte Carlo model in parallel code SPPARKS to simulate grain growth and recrystallization under a residual energy distribution. Our model incorporates physical properties and real-time, as opposed to generic properties and relative time. We further validated that our SPPARKS simulations reproduced the expected kinetic behavior of single-grain evolution. We then used the validated approach to simulate the anneal of an additively manufactured material under the same conditions used in our experiments. We found that the cubic grains can origin from a periodically varying residual energy that may be present in additively manufactured materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The data acquisition system of the LZ dark matter detector: FADR

The Data Acquisition System (DAQ) for the LUX-ZEPLIN (LZ) dark matter detector is described. The signals from 745 PMTs, distributed across three subsystems, are sampled with 100-MHz 32-channel digitizers (DDC-32s). A basic waveform analysis is carried out on the on-board Field Programmable Gate Arrays (FPGAs) to extract information about the observed scintillation and electroluminescence signals. This information is used to determine if the digitized waveforms should be preserved for offline analysis. The system is designed around the Kintex-7 FPGA. In addition to digitizing the PMT signals and providing basic event selection in real time, the flexibility provided by the use of FPGAs allows us to monitor the performance of the detector and the DAQ in parallel to normal data acquisition. Furthermore, the hardware and software/firmware of this FPGA-based Architecture for Data acquisition and Realtime monitoring (FADR) are discussed and performance measurements are described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

A parallel strategy for density functional theory computations on accelerated nodes

Using the Löwdin orthonormalization of tall-skinny matrices as a proxy-app for wavefunction-based Density Functional Theory solvers, we investigate a distributed memory parallel strategy focusing on Graphics Processing Unit (GPU)-accelerated nodes as available on some of the top ranked supercomputers at the present time. Here we present numerical results in the strong limit regime, as it is particularly relevant for First-Principles Molecular Dynamics. We also examine how matrix product-based iterative solvers provide a competitive alternative to dense eigensolvers on GPUs, allowing to push the strong scaling limit of these computations to a larger number of distributed tasks. Our strategy, which relies on replicated Gram matrices and efficient collective communications using the NCCL library, leads to a time-to-solution under 0.5 s for the Löwdin orthonormalization of a tall-skinny matrix of 3000 columns on Summit at Oak Ridge Leadership Facility (OLCF). Given the similarity in computational operations between one iteration of a DFT solver and this proxy-app, this shows the possibility of solving accurately the DFT equations well under a minute for 3000 electronic wave functions, and thus perform First-Principles molecular dynamics of physical systems much larger than traditionally solved on CPU systems.

97 MATHEMATICS AND COMPUTING↗