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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Scale-dependent spatial variabilities of hydrological exchange flows and transit time in a large regulated river

Hydrological exchange flows (HEF) across the river-aquifer interface and the associated residence time of river water in the aquifer have important implications for contaminant plume migration and biogeochemical processes in the river corridor. HEFs and residence time are influenced by both subsurface physical features and hydrologic forcing related to the transport process, which can exhibit complex spatial and temporal variations. In this study, we used a massively parallel subsurface flow model and a particle-tracking model to study the influences of different control factors on spatial variability of HEFs and residence time distributions (RTD) in the Hanford Reach of the Columbia River in Washington State. A total number of 100M particles were randomly injected in time and space and then tracked in a model domain that covers a 51-km 2 area (15.1M model cells). We used hourly river stages and groundwater levels to drive the model to provide dynamic velocity fields for the particle tracking in the simulation period that was longer than 2 years. The groundwater flow simulation and particle-tracking results provide the first comprehensive assessment of the spatial distribution of HEFs and residence time in large complex river corridors. Overall, our results show that the aquifer hydrogeological structure has the strongest correlation with the extent and magnitude of exchange flux. The residence time exhibits complex patterns that are impacted by all the river geomorphologic, hydrodynamic, and hydrogeologic factors and are strongly correlated with the downwelling ratio of exchange flux. The new insights gained through this study can be used to support the development of reduced-order models of HEFs and RTDs for large complex river systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

Passive microwave remote sensing of precipitation from space - Impact of cloud microphysics

Numerical experiments to investigate the impact of time-dependent cloud microphysical structure on the transfer to space of passive microwave radiation are carried out, using a vertically and angularly detailed, plane-parallel radiative transfer model for unpolarized microwave radiation. The object is to explore the feasibility of multichannel passive-microwave retrieval techniques to estimate precipitation from space. Results on the time-dependent influence of large graupel-size ice particles on top-of-atmosphere brightness temperature - rainfall rate relationships at 6.6, 10.7, 18, 19.35, 22.235, 37, 89.5, 130, and 231 GHz are given. Frequency-dependent vertical profiles of a generalized brightness-temperature weighting function, introduced to vertically resolve the contributions by individual cloud and precipitation layers to the brightness temperatures, are also shown.

Mugnai, Alberto↗
NASA NTRS1979

Autonomous attitude estimation via star sensing and pattern recognition

Results are reported on the development of an autonomous, onboard, near real time spacecraft attitude estimation technique. The approach uses CCD based star sensors to digitize relative star positions. Three microcomputers are envisioned, configured in parallel, to: (1) determine star image centroids and delete spurious images; (2) identify measured stars with stars in an onboard catalog and determine discrete attitude estimates; (3) integrate gyro rate measurements and determine optimal real time attitude estimates for use in the control system and for feedback to the star identification algorithm. Algorithms for the star identification are presented. The discrete attitude estimation algorithm recovers thermally varying interlock angles between two star sensors. The optimal state estimation process recovers rate gyro biases in addition to real time attitude estimates.

Junkins, J. L.↗
DOE OSTI2020

Real-time optimization of multi-cell industrial evaporative cooling towers using machine learning and particle swarm optimization

Existing electrical generating stations must operate with greater flexibility due to increasing renewable energy penetration on the electrical grid, and many coal-fired power stations have transitioned away from baseload operation to load-following operation to aid in grid stability. In cases where multiple independently controlled cooling tower cells are used in parallel for the cooling purposes of such stations, there is an opportunity to increase plant efficiency through data-driven optimization across their full load ranges. This work presents a novel application of real-time optimization using machine learning and particle swarm optimization on a multi-cell induced-draft cooling tower servicing a coal-fired power station under variable load. This is the first work to demonstrate simultaneous optimization of a multi-cell cooling tower, in addition to using machine learning for closed-loop control on a cooling tower. A novel control configuration is presented that ensures original control logic is not adversely affected and that the overall plant process is not disrupted using only existing hardware and operational data. To verify this methodology, the 12 independent cooling tower cells are simulated in parallel using historic operating data to demonstrate the effectiveness of real-time optimization compared to current practice. An artificial neural network is trained to predict overall cooling tower power consumption using only operational data and ambient conditions with an R2 value of greater than 0.96. The real-time optimization using particle swarm yields 6.7% annual energy usage savings compared to current practices, although the extent of the real-time savings varies greatly with both plant load and environmental conditions. This is particularly significant for a variable load situation because frequent ramping typically results in reduced overall efficiency. Furthermore, this proposed AI-based solution presents an opportunity to improve the overall heat rate of a load-following coal-fired power plant without the need to perform extensive first-principles modeling or add additional hardware to the cooling tower, resulting in more resources conserved and less overall emissions per unit of electricity generated.

42 ENGINEERING↗
DOE OSTI2023

Online Monitoring of Medium Voltage Cable Systems with Spread Spectrum Time Domain and Frequency Domain Reflectometry

In-service failures of wave energy convertor (WEC) cable systems can have a significant cost and power availability impact. Close parallel research 2019 data showed > 1B£ and 9 Terra-Watt-Hours associated with global off-shore wind (OSW) cable failures (Strang-Moran 2020). OSW is a closely related technology but currently is significantly cheaper than WEC technology. For wave energy to compete, the problem of reliable cable transmission must be mitigated. This project develops isolation technology to allow online high frequency reflectometry testing of medium voltage cables (1 to 10 kV and higher) without arcing or damage to the test instrument. Online spread spectrum time domain reflectometry (SSTDR) testing has been established for low voltage cable systems in the aircraft and rail industry and the ability to detect and locate cable flaws of interest is well understood. Extending reflectometry testing to medium voltage systems could enable detection of cable damage before failures occur thereby allowing repair and replacement of damaged cable segments to be scheduled and managed. The seedling project succeeded to pass and receive high frequency SSTDR signals onto a cable up to 1 kV using a parallel trace isolation circuit board that can be connected onto the test cable. The approach used a novel circuit design for which an invention disclosure has been filed. A proposed sapling project would extend the technology toward the higher operating voltages used by WEC systems, thereby enabling online SSTDR cable monitoring. The goal of the seedling project was to extend the capability of the ARENA cable/motor test bed to address medium voltages and to develop a high pass filter isolation architecture to protect the reflectometry instrument from the low frequency (DC – 60 Hz) line voltage while allowing the high frequency diagnostic signal to pass to and from the test instrument to the live line. Initial efforts focused on passive LCR filter circuits to reduce 60 Hz levels below 10 volts from a 10 kV line while allowing the MHz high frequency chirps to pass onto the cables and for mV signals to be detected. We discovered that the parasitic loss behavior of real high voltage components precluded this approach from working. An alternate approach was adapted for the electric field to couple between two parallel traces on a printed circuit board much like a radio-frequency coupler. The challenge here was and is to have the parallel traces close enough to each other to effectively pass the high frequency chirp onto the live line and receive any reflected signal from any encountered impedance change along the cable. This reflected signal will be in the mV range. The traces however must be far enough apart to not allow arcing on the board. A design with 3 mm spacing was determined to allow the high frequency signal to pass onto the live line and receive the mV signal back into the instrument while reducing the 60 Hz voltage amplitude by >80 dB (more than a factor of 10,000) without allowing arcing from across the parallel traces. This was confirmed by simulation and test.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗
DOE OSTI2021

Enabling rapid COVID-19 small molecule drug design through scalable deep learning of generative models

We improved the quality and reduced the time to produce machine learned models for use in small molecule antiviral design. Our globally asynchronous multi-level parallel training approach strong scales to all of Sierra with up to 97.7% efficiency. We trained a novel, character-based Wasserstein autoencoder that produces a higher quality model trained on 1.613 billion compounds in 23 minutes while the previous state of the art takes a day on 1 million compounds. Reducing training time from a day to minutes shifts the model creation bottleneck from computer job turnaround time to human innovation time. Our implementation achieves 318 PFLOPs for 17.1% of half-precision peak. We will incorporate this model into our molecular design loop enabling the generation of more diverse compounds; searching for novel, candidate antiviral drugs improves and reduces the time to synthesize compounds to be tested in the lab.

Jacobs, Sam Ade↗
DOE OSTI2023

Microstructurally driven self-sharpening mechanism in beaver incisor enamel facilitates their capacity to fell trees

Beavers ( Castor ) stand out among mammals for their unique capacity to fell trees using their large, ever-growing incisors. This routine consumption of resistant fodder induces prodigious wear in the lower incisors, despite this blunting effect the incisors maintain a remarkably sharp cutting edge. Notably, the enamel edges of their incisors show a highly complex two-part microstructure of which the biomechanical import is unknown. Here, using fracture analysis, nanoindentation, and wear testing on North American beaver ( C. canadensis ) incisors we test the microstructure's possible contribution to maintaining incisal sharpness. Although comparable in hardness, the inner enamel preferentially fails and readily wears at 2.5 times the rate of the outer enamel. The outer microstructure redirects all fractures in parallel, decreasing fracture coalescence. Conversely, the inner microstructure facilitates crack coalescence increasing the wear rate by isolating layers of enamel prisms that readily fragment. Together these two architectures form a microstructurally driven self-sharpening mechanism contained entirely within the thin enamel shell. Here, our results demonstrate that enamel microstructures exposed at the occlusal surface can markedly influence both enamel crest shape and surface texture in wearing dentitions. The methods introduced here open the door to exploring the biomechanical functionality and evolution of enamel microstructures throughout Mammalia.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

An optical-input Maximum Likelihood Estimation feedback system demonstrated on tokamak horizontal equilibrium control

A readily parallelized Maximum Likelihood Estimation (MLE) algorithm with linear computational complexity is demonstrated in real time using only measurements from an extreme ultraviolet (EUV) diagnostic to control the horizontal position of a tokamak plasma. A set of trial emissivity profiles are parameterized by the control quantity of interest (R m ), and the MLE is identified from the profile which minimizes the signal reconstruction residual. The algorithm depends on an empirically determined likelihood function with exponential form. EUV emission (λ ≈ 15eV-1keV) is captured in a poloidal plane by four 16-channel AXUV diodes mounted at different poloidal angles with radial and angular resolution sufficient to discern plasma equilibrium evolution in HBT-EP. Calculations of the plasma major radius by the system are consistent within diagnostic uncertainty for the majority of the discharge with those of: a weighted average of vertical soft X-ray or EUV chords, magnetic sensors, and an equilibrium reconstruction. The feedback system corrects for a horizontal displacement of the major radius equal to 20% of the plasma minor radius by adjusting the vertical field produced from 40 in-vessel control coils in real time. The MLE calculation is performed on a GPU in a 15 μs cycle, with similar performance in this application to a simple weighted average of vertical chords. Finally, results demonstrate horizontal position control using magnetic actuators and an optical observer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Carburization and decarburization behavior of Grade 91 ferritic-martensitic steel in liquid sodium environments

Here, this paper presents a study of carbon transfer and its effect on microstructure and tensile properties of Grade 91 (G91) ferritic-martensitic steel exposed to sodium at 550–650 °C. Sodium exposure tests were conducted in Argonne's forced convection sodium loops up to exposure times of ~40,000 h. Thermal aging study of G91 steel was conducted in parallel to isolate the thermal aging effect from the sodium effect. It was found that sodium exposures at 650 °C dissolved M 23 C 6 carbides, eliminated the martensite subgrain structure resulting in excessive grain growth and reduced the tensile strength by >50%, while sodium exposures at 550 and 600 °C had an insignificant effect on its microstructure and tensile properties. These effects were attributed to the carburization/decarburization process of G91 steel in sodium environments. Carbon concentrations in sodium were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8–1.2 ppm in the SMT-1 loop and 0.3–0.7 ppm in the SMT-2 loop. Thermodynamic analysis of the carburization – decarburization process was conducted for G91 steels exposed in sodium environments. The carbon activity-concentration relationship for G91 was evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC and VC carbides dominated the process at low carbon activities. The calculated carburization-decarburization boundary showed that G91 would undergo decarburization at 650 °C and carburization at 550 °C in the sodium loop environments, which was consistent with our experimental observations. This experimental and theoretical analysis provided a basis for predicting the effect of carbon transfer on the integrity of reactor components in sodium environments and for the design of new alloys used in sodium-cooled fast reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

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

Whistler wave generation by anisotropic tail electrons during asymmetric magnetic reconnection in space and laboratory

Whistler wave generation near the magnetospheric separatrix during reconnection at the dayside magnetopause is studied with data from the Magnetospheric Multiscale (MMS) mission. The dispersion relation of the whistler mode is measured for the first time near the reconnection region in space, which shows that whistler waves propagate nearly parallel to the magnetic field line. A linear analysis indicates that the whistler waves are generated by temperature anisotropy in the electron tail population. This is caused by loss of electrons with a high velocity parallel to the magnetic field to the exhaust region. Here, there is a positive correlation between activities of whistler waves and the lower–hybrid drift instability (LHDI) both in laboratory and space, indicating the enhanced transport by LHDI may be responsible for the loss of electrons with a high parallel velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Oxygen transport properties estimation by DSMC-CT simulations

Coupling DSMC simulations with classical trajectories calculations is emerging as a powerful tool to improve predictive capabilities of computational rarefied gas dynamics. The considerable increase of computational effort outlined in the early application of the method (Koura,1997) can be compensated by running simulations on massively parallel computers. In particular, GPU acceleration has been found quite effective in reducing computing time (Ferrigni,2012; Norman et al.,2013) of DSMC-CT simulations. The aim of the present work is to study rarefied Oxygen flows by modeling binary collisions through an accurate potential energy surface, obtained by molecular beams scattering (Aquilanti, et al.,1999). The accuracy of the method is assessed by calculating molecular Oxygen shear viscosity and heat conductivity following three different DSMC-CT simulation methods. In the first one, transport properties are obtained from DSMC-CT simulations of spontaneous fluctuation of an equilibrium state (Bruno et al, Phys. Fluids, 23, 093104, 2011). In the second method, the collision trajectory calculation is incorporated in a Monte Carlo integration procedure to evaluate the Taxman’s expressions for the transport properties of polyatomic gases (Taxman,1959). In the third, non-equilibrium zero and one-dimensional rarefied gas dynamic simulations are adopted and the transport properties are computed from the non-equilibrium fluxes of momentum and energy. The three methods provide close values of the transport properties, their estimated statistical error not exceeding 3%. The experimental values are slightly underestimated, the percentage deviation being, again, few percent.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Implementation and benchmark of a long-range corrected functional in the density functional based tight-binding method

Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data.

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

The flush-mounted rail Langmuir probe array designed for the Alcator C-Mod vertical target plate divertor

An array of flush-mounted and toroidally elongated Langmuir probes (henceforth called rail probes) have been specifically designed for the Alcator C-Mod’s vertical target plate divertor and operated over multiple campaigns. The “flush” geometry enables the tungsten electrodes to survive high heat flux conditions in which traditional “proud” tungsten electrodes suffer damage from melting. The toroidally elongated rail-like geometry reduces the influence of sheath expansion, which is an important effect to consider in the design and interpretation of flush-mounted Langmuir probes. The new rail probes successfully operated during C-Mod’s FY2015 and FY2016 experimental campaigns with no evidence of damage, despite being regularly subjected to heat flux densities parallel to the magnetic field exceeding ~1 GW m -2 for short periods of time. A comparison between rail and proud probe data indicates that sheath expansion effects were successfully mitigated by the rail design, extending the use of these Langmuir probes to incident magnetic field line angles as low as 0.5°.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Triplet-pair spin signatures from macroscopically aligned heteroacenes in an oriented single crystal

The photo-driven process of singlet fission generates coupled triplet pairs (TT) with fundamentally intriguing and potentially useful properties. The quintet 5 TT 0 sublevel is particularly interesting for quantum information because it is highly entangled, is addressable with microwave pulses, and could be detected using optical techniques. Previous theoretical work on a model Hamiltonian and nonadiabatic transition theory, called the JDE model, has determined that this sublevel can be selectively populated if certain conditions are met. Among the most challenging, the molecules within the dimer undergoing singlet fission must have their principal magnetic axes parallel to one another and to an applied Zeeman field. Here, we present time-resolved electron paramagnetic resonance (TR-EPR) spectroscopy of a single crystal sample of a tetracenethiophene compound featuring arrays of dimers aligned in this manner, which were mounted so that the orientation of the field relative to the molecular axes could be controlled. The observed spin sublevel populations in the paired TT and unpaired (T+T) triplets are consistent with predictions from the JDE model, including preferential 5 TT 0 formation at z ‖ B 0 , with one caveat—two 5 TT spin sublevels have little to no population. This may be due to crossings between the 5 TT and 3 TT manifolds in the field range investigated by TR-EPR, consistent with the intertriplet exchange energy determined by monitoring photoluminescence at varying magnetic fields.

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

SAIGE-GPU: accelerating genome- and phenome-wide association studies using GPUs

Genome-wide association studies (GWAS) at biobank scale are computationally intensive, especially for admixed populations requiring robust statistical models. SAIGE is a widely used method for generalized linear mixed-model GWAS but is limited by its CPU-based implementation, making phenome-wide association studies impractical for many research groups. We developed SAIGE-GPU, a GPU-accelerated version of SAIGE that replaces CPU-intensive matrix operations with GPU-optimized kernels. The core innovation is distributing genetic relationship matrix calculations across GPUs and communication layers. Applied to 2068 phenotypes from 635 969 participants in the Million Veteran Program, including diverse and admixed populations, SAIGE-GPU achieved a 5-fold speedup in mixed model fitting on supercomputing infrastructure and cloud platforms. We further optimized the variant association testing step through multi-core and multi-trait parallelization. Deployed on Google Cloud Platform and Azure, the method provided substantial cost and time savings. Source code and binaries are available for download at https://github.com/saigegit/SAIGE/tree/SAIGE-GPU-1.3.3. A code snapshot is archived at Zenodo for reproducibility (DOI: [10.5281/zenodo.17642591]). SAIGE-GPU is available in a containerized format for use across HPC and cloud environments and is implemented in R/C++ and runs on Linux systems.

Rodriguez, Alex [Argonne National Laboratory (ANL)↗
DOE OSTI2025

Entanglement-enhanced ac magnetometry in the presence of Markovian noise

Entanglement is a resource to improve the sensitivity of quantum sensors. In an ideal case, using an entangled state as a probe to detect target fields, we can beat the standard quantum limit by which all classical sensors are bounded. However, since entanglement is fragile against decoherence, it is unclear whether entanglement-enhanced metrology is useful in a noisy environment. Its benefit is indeed limited when estimating the amplitude of dc magnetic fields under the effect of parallel Markovian decoherence, where the noise operator is parallel to the target field. In this paper, on the contrary, we show an advantage to using an entanglement over the classical strategy under the effect of parallel Markovian decoherence when we try to detect ac magnetic fields. We consider a scenario to induce a Rabi oscillation of the qubits with the target ac magnetic fields. Although we can, in principle, estimate the amplitude of the ac magnetic fields from the Rabi oscillation, the signal becomes weak if the qubit frequency is significantly detuned from the frequency of the ac magnetic field. We show that, by using the Greenberger-Horne-Zeilinger (GHZ) states, we can significantly enhance the signal of the detuned Rabi oscillation even under the effect of parallel Markovian decoherence. Further, our method is based on the fact that the interaction time between the GHZ states and ac magnetic fields scales as 1/L to mitigate the decoherence effect, where L is the number of qubits, which contributes to improving the bandwidth of the detectable frequencies of the ac magnetic fields. Our results pave the way for new applications of entanglement-enhanced ac magnetometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗