Search NASA⌕ Search

SEARCH · Search NASA

Results for “Widely Used”

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 541 records · Page 30

Development of the Architectural Simulation Model for Future Launch Systems and its Application to an Existing Launch Fleet

A significant portion of lifecycle costs for launch vehicles are generated during the operations phase. Research indicates that operations costs can account for a large percentage of the total life-cycle costs of reusable space transportation systems. These costs are largely determined by decisions made early during conceptual design. Therefore, operational considerations are an important part of vehicle design and concept analysis process that needs to be modeled and studied early in the design phase. However, this is a difficult and challenging task due to uncertainties of operations definitions, the dynamic and combinatorial nature of the processes, and lack of analytical models and the scarcity of historical data during the conceptual design phase. Ultimately, NASA would like to know the best mix of launch vehicle concepts that would meet the missions launch dates at the minimum cost. To answer this question, we first need to develop a model to estimate the total cost, including the operational cost, to accomplish this set of missions. In this project, we have developed and implemented a discrete-event simulation model using ARENA (a simulation modeling environment) to determine this cost assessment. Discrete-event simulation is widely used in modeling complex systems, including transportation systems, due to its flexibility, and ability to capture the dynamics of the system. The simulation model accepts manifest inputs including the set of missions that need to be accomplished over a period of time, the clients (e.g., NASA or DoD) who wish to transport the payload to space, the payload weights, and their destinations (e.g., International Space Station, LEO, or GEO). A user of the simulation model can define an architecture of reusable or expendable launch vehicles to achieve these missions. Launch vehicles may belong to different families where each family may have it own set of resources, processing times, and cost factors. The goal is to capture the required resource levels of the major launch elements and their required facilities. The model s output can show whether or not a certain architecture of vehicles can meet the launch dates, and if not, how much the delay cost would be. It will also produce aggregate figures of missions cost based on element procurement cost, processing cost, cargo integration cost, delay cost, and mission support cost. One of the most useful features of this model is that it is stochastic where it accepts statistical distributions to represent the processing times mimicking the stochastic nature of real systems.

Rabadi, Ghaith↗

Pharmacokinetics of Acetaminophen in Hind Limbs Unloaded Mice: A Model System Simulating the Effects of Low Gravity on Astronauts in Space

The pharmacokinetics (PK) of medications administered to astronauts could be altered by the conditions in Space. Low gravity and free floating (and associated hemodynamic changes) could affect the absorption, distribution, metabolism and excretion of the drugs. Knowledge of these alterations is essential for adjusting the dosage and the regimen of drug administration in astronauts. Acquiring of such knowledge has inherent difficulties due to limited opportunities for experimenting in Space. One of the approaches is to use model systems that simulate some of the Space conditions on Earth. In this study we used hind limbs unloaded mice (HLU) to investigate the possible changes in PK of acetaminophen, a widely used analgesic with high probability of use by astronauts. The HLU is recognized as an appropriate model for simulating the effects of low gravity on hemodynamic parameters. Mice were tail suspended (n = 24) for 24-96 hours prior to introduction of acetaminophen (150 - 300 mg/kg). The drug (in aqueous solution containing 10% ethyl alcohol by volume) was given orally by a gavage procedure and after the administration of acetaminophen mice were additionally suspended for 30 min, 1 and 2 hours. Control mice (n = 24) received the same dose of acetaminophen and were kept freely all the time. Blood specimens were obtained either from retroorbital venous sinuses or from heart. Acetaminophen concentration was measured in plasma by the fluorescent polarization immunoassay and the AxSYM analyzer (Abbott Laboratories). In control mice peak acetaminophen concentration was achieved at 30 min. By 1 hour the concentration decreased to less than 50% of the peak level and at 2 hours the drug was almost undetectable in the serum. HLU for 24 hours significantly altered the acetaminophen pharmacokinetic: at 30 min the acetaminophen concentrations were significantly (both statistically and medically significant) lower than in control mice. The concentrations also reduced less significantly after 1 and 2 hours. At 2 hours approximately 20% of the drug still remained in the circulation. After 96 hrs of HLU the changes in acetaminophen PK were less prominent. These data indicate that short term HLU causes significant changes in acetaminophen PK most likely associated with HUL-related hemodynamic changes. However, after 96 hour these changes diminished. This suggests hemodynamic adaptation to the HUL conditions that possibly occurs also in real space conditions.

Peterson, Amanda↗

Numerical Predictions of Sonic Boom Signatures for a Straight Line Segmented Leading Edge Model

A sonic boom wind tunnel test was conducted on a straight-line segmented leading edge (SLSLE) model in the NASA Langley 4- by 4- Foot Unitary Plan Wind Tunnel (UPWT). The purpose of the test was to determine whether accurate sonic boom measurements could be obtained while continuously moving the SLSLE model past a conical pressure probe. Sonic boom signatures were also obtained using the conventional move-pause data acquisition method for comparison. The continuous data acquisition approach allows for accurate signatures approximately 15 times faster than a move-pause technique. These successful results provide an incentive for future testing with greatly increased efficiency using the continuous model translation technique with the single probe to measure sonic boom signatures. Two widely used NASA codes, USM3D (Navier-Stokes) and CART3D-AERO (Euler, adjoint-based adaptive mesh), were used to compute off-body sonic boom pressure signatures of the SLSLE model at several different altitudes below the model at Mach 2.0. The computed pressure signatures compared well with wind tunnel data. The effect of the different altitude for signature extraction was evaluated by extrapolating the near field signatures to the ground and comparing pressure signatures and sonic boom loudness levels.

Elmiligui, Alaa A.↗

Empirical and Face Validity of Software Maintenance Defect Models Used at the Jet Propulsion Laboratory

At the Mission Design and Navigation Software Group at the Jet Propulsion Laboratory we make use of finite exponential based defect models to aid in maintenance planning and management for our widely used critical systems. However a number of pragmatic issues arise when applying defect models for a post-release system in continuous use. These include: how to utilize information from problem reports rather than testing to drive defect discovery and removal effort, practical model calibration, and alignment of model assumptions with our environment.

software reliability↗

Harmonizing solar induced fluorescence across spatial scales, instruments, and extraction methods using proximal and airborne remote sensing: A multi-scale study in a soybean field

Solar-induced chlorophyll fluorescence (SIF) has been widely used to track vegetation photosynthesis at different scales ranging from in-situ measurements to satellite products. Airborne platforms sample SIF data at a spatial scale intermediate between in-situ and satellite, matching that of ground measurement (e.g. flux tower footprints and other field sampling), enabling us to explore causes of SIF variation and validate satellite-based SIF products. However, harmonizing SIF across sensors and platforms (correcting for systematic errors to yield a consistent, comparable SIF product) is challenging because SIF can be retrieved in different absorption windows, with different instruments and methods complicating the comparison between different observational levels (i.e., ground, airborne, satellites) and between sites equipped with different instruments with varying optical properties (spectral resolution and sampling intervals, spatial resolution). Additionally, the spatial and temporal variability of atmospheric properties can influence the retrieval of the weak SIF signal. Because of these complications, direct comparisons of airborne and ground SIF across scales are rarely attempted. Here, in this study, we combined airborne SIF data with simultaneous ‘ground truth’ data collected by stationary and mobile platforms in a soybean field in Nebraska, USA. In this effort, we tested several SIF extraction methods, including Fraunhofer Line Discrimination (FLD), improved Fraunhofer Line Discrimination (iFLD), Spectral Fitting Method (SFM), SpecFit, and a Singular Vector Decomposition (SVD) method. The SpecFit method was sensitive to the 715–740 nm water bands and removing the water bands in the fitting process yielded better agreement between the airborne and ground SIF spectra. Accurate estimation of the ground level downwelling irradiance obtained by ground measurements over a calibration target improved agreement between airborne and ground SIF retrievals at the O 2 A band, and allowed us to derive a SIF dataset with improved agreement across platforms and sampling scales. This experimental approach provided a method for generating comparable SIF signals across instruments, methods and platforms, which is critical to understanding the SIF-GPP relationship at different scales and to cross-validate the diversity of platforms used for satellite products calibration and validation.

54 ENVIRONMENTAL SCIENCES↗

Investigation of the effect of structural defects from hydride precipitation on superconducting properties of high purity SRF cavity Nb using magneto-optical and electron imaging methods

Abstract Nb is an elemental superconductor with a critical temperature of 9.3 K and is widely used to fabricate superconducting radiofrequency (SRF) cavities for particle accelerators. However, microstructural defects in Nb, such as grain boundaries (GBs) and dislocations, can act as pinning centers for magnetic flux that can degrade SRF cavity performance. Hydrogen contamination is also detrimental to SRF cavity performance due to the formation of normal conducting hydrides during cool down. In this study, disc shaped Nb bi-crystals extracted from high-purity large-grain Nb slices were investigated to study the effects of GBs, hydrogen, and dislocations on superconducting properties. Grain orientation and GB misorientation were measured using Laue x-ray diffraction and electron backscattered diffraction (EBSD) analyses. Cryogenic magneto-optical imaging was used to directly observe magnetic flux penetration below T c = 9.3 K. Damage caused by low temperature precipitation of hydrides and their dissolution upon reheating after cryogenic cycles was examined using electron channeling contrast imaging, and EBSD. The relationships between hydride formation, dislocation content, GBs, cryo-cooling, heat treatment (HT), and flux penetration indicate that both GB character and hydrogen content affect magnetic flux penetration. Such flux penetration could be facilitated by dislocation structures and low angle GBs resulting from hydride precipitation and HT.

36 MATERIALS SCIENCE↗

High speed corner and gap-seal computations using an LU-SGS scheme

The hybrid Lower-Upper Symmetric Gauss-Seidel (LU-SGS) algorithm was added to a widely used series of 2D/3D Euler/Navier-Stokes solvers and was demonstrated for a particular class of high-speed flows. A limited study was conducted to compare the hybrid LU-SGS for approximate Newton iteration and diagonalized Beam-Warming (DBW) schemes on a work and convergence history basis. The hybrid LU-SGS algorithm is more efficient and easier to implement than the DBW scheme originally present in the code for the cases considered. The code was validated for the hypersonic flow through two mutually perpendicular flat plates and then used to investigate the flow field in and around a simplified scramjet module gap seal configuration. Due to the similarities, the gap seal flow was compared to hypersonic corner flow at the same freestream conditions and Reynolds number.

Coirier, William J.↗

High speed corner and gap-seal computations using an LU-SGS scheme

The hybrid Lower-Upper Symmetric Gauss-Seidel (LU-SGS) algorithm was added to a widely used series of 2D/3D Euler/Navier-Stokes solvers and was demonstrated for a particular class of high-speed flows. A limited study was conducted to compare the hybrid LU-SGS for approximate Newton iteration and diagonalized Beam-Warming (DBW) schemes on a work and convergence history basis. The hybrid LU-SGS algorithm is more efficient and easier to implement than the DBW scheme originally present in the code for the cases considered. The code was validated for the hypersonic flow through two mutually perpendicular flat plates and then used to investigate the flow field in and around a simplified scramjet module gap seal configuration. Due to the similarities, the gap seal flow was compared to hypersonic corner flow at the same freestream conditions and Reynolds number.

Coirier, William J.↗

Automating CPM-GOMS

CPM-GOMS is a modeling method that combines the task decomposition of a GOMS analysis with a model of human resource usage at the level of cognitive, perceptual, and motor operations. CPM-GOMS models have made accurate predictions about skilled user behavior in routine tasks, but developing such models is tedious and error-prone. We describe a process for automatically generating CPM-GOMS models from a hierarchical task decomposition expressed in a cognitive modeling tool called Apex. Resource scheduling in Apex automates the difficult task of interleaving the cognitive, perceptual, and motor resources underlying common task operators (e.g. mouse move-and-click). Apex's UI automatically generates PERT charts, which allow modelers to visualize a model's complex parallel behavior. Because interleaving and visualization is now automated, it is feasible to construct arbitrarily long sequences of behavior. To demonstrate the process, we present a model of automated teller interactions in Apex and discuss implications for user modeling. available to model human users, the Goals, Operators, Methods, and Selection (GOMS) method [6, 21] has been the most widely used, providing accurate, often zero-parameter, predictions of the routine performance of skilled users in a wide range of procedural tasks [6, 13, 15, 27, 28]. GOMS is meant to model routine behavior. The user is assumed to have methods that apply sequences of operators and to achieve a goal. Selection rules are applied when there is more than one method to achieve a goal. Many routine tasks lend themselves well to such decomposition. Decomposition produces a representation of the task as a set of nested goal states that include an initial state and a final state. The iterative decomposition into goals and nested subgoals can terminate in primitives of any desired granularity, the choice of level of detail dependent on the predictions required. Although GOMS has proven useful in HCI, tools to support the construction of GOMS models have not yet come into general use.

GOMS↗

Performance evaluation of a Terry GS-2 steam impulse turbine with air-water mixtures

Terry steam turbines are widely used in various industries because of their robust design. Within the nuclear power generation industry, they are used in the Reactor Core Isolation Cooling System to remove decay heat during reactor isolation events. During the Fukushima Daiichi nuclear power station disaster in Japan in 2011, the Reactor Core Isolation Cooling System and associated Terry turbine operated for over 70 hours in Unit 2; this runtime is well beyond the expected operating duration. Theories suggest the turbine was subjected to a two-phase inlet flow, which could degrade the turbine performance. In this work, an experimental test rig was constructed to test a full-scale Terry model GS-2 steam turbine under two-phase air/water flows. Steady-state efficiency and torque performance maps of the turbine were developed over a range of turbine inlet pressures (1.38–4.83 bar or 20–70 psia), air mass fractions (0.05–1.0) and rotational speeds up to 4000 RPM. Furthermore, turbine performance followed expected trends with torque varying linearly and efficiency varying quadratically with rotational speed. In addition, high-speed images of the two-phase flow entering the turbine were also analyzed to understand how changes in inlet pressure and air mass fraction affect the flow regime and homogenization. The present tests with air–water two-phase mixtures are an important step towards providing an understanding of the full-scale Terry turbine’s behavior and performance curves under two-phase conditions. The results of this work will be combined with air/water and steam/water data gathered using a small-scale Terry ZS-1 steam turbine in order to understand the scaling relationship between large and small size Terry turbines and fluid pairs. The combined data set will enable further development of analytical models over a wide range of conditions and may be used to provide technical justification for expanded use of the Terry turbines in nuclear power plant safety systems and other systems.

42 ENGINEERING↗

Correction for nonlinear photon counting effects in lidar systems

Photomultiplier tubes (PMT's) employed in the photon counting (PC) mode of operation are widely used as detectors in lidar systems. In our laboratory, we have developed a versatile Nd:YAG lidar which is used for measurement of both the middle atmosphere and the troposphere. With this system, we encounter a very wide range of signal levels ranging from the extremely weak signals from the top of the mesosphere to the very strong returns from low level clouds. Although the system is capable of operating the PMT's in either the analog detection or photon counting mode, we find that often when we use photon counting we have portions of our lidar return which contain very useful information but are not within the linear operating regime of the PC system. We report the results of our efforts to explore the extent to which such high intensity PC signals can be quantitatively analyzed. In particular, a useful model relating the mean 'true' count rate and the observed count rate is presented and it's application to our system demonstrated. This model takes into account the variation in height of the PMT output pulses and the effect of the pulse height discrimination threshold.

Donovan, D. P.↗

New Dimensions in the Theory of Excited States and X-ray Spectra (Final Report)

This Final Technical Report briefly summarizes the achievements during the lifetime of our DOE BES grant DE-FG02-97ER45623. The long-term goal of this project has been the development of quantitative theories of the interaction between radiation and matter, with a focus on x-ray spectroscopies. X-ray spectra have long been among the most important probes of atomic-scale structure and properties of matter, ranging from atoms and molecular systems to condensed matter and exotic states. These spectroscopies are widely used in investigations at the major DOE synchrotron x-ray facilities and related centers world-wide. In addition to fundamental theory, a major goal of our project has been the development of computational software that implements the theory for calculations of x-ray spectra of various materials throughout the periodic table. Due to the complex nature of x-ray spectra, quantitative theory is essential for its interpretation. The theory is challenging since it involves excited state electronic structure and many-body correlation effects that go beyond independent particle approximations like DFT or Hartree-Fock. Moreover, the experimental investigations typically involve a broad range of energy, time, and temperature scales, from the UV-Vis to hard x-ray energies of order 10 4 eV, and temperatures T from ambient up to the warm-dense-matter regime where the Fermi energy kBTF is of order a few eV, i.e., temperatures of order 105 Kelvin. This broad range of experimental conditions has fostered many novel theoretical approaches and computational techniques, many of which we have developed systematically over the duration of the grant. In contrast to the traditional wave-function approach of quantum theory and electronic structure methods, our theoretical approach is based on modern Green's function techniques. This approach is better suited for aperiodic structures, excited states, and broad spectral ranges, since it avoids the computational bottlenecks of sum-over-states approaches, as in the Fermi golden rule. This theoretical framework has been incorporated into efficient, user-friendly x-ray spectroscopy software which is now used routinely worldwide to simulate and analyze spectra. These theoretical tools provide an essential complement to synchrotron and next-generation light sources, which are used to investigate complex materials with ever increasing precision. Moreover, the synergism between theory, computation and experiment contributed by our research enhances scientific understanding and creates opportunities for innovations in materials and energy science and in many fields. As documented in this Report, this research grant has been remarkably successful in achieving these goals. In particular, this grant has supported the development of the x-ray spectroscopy software suite known as FEFF (named for an effective scattering amplitude f eff in the theory). The FEFF codes have become one of the premier tools for quantitative simulations of x-ray spectra as documented by many thousands of citations in the Web of Science and Google-Scholar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Non-OH Chemistry in Oxidation Flow Reactors for the Study of Atmospheric Chemistry Systematically Examined by Modeling

Oxidation flow reactors (OFRs) using low-pressure Hg lamp emission at 185 and 254 nm produce OH radicals efficiently and are widely used in atmospheric chemistry and other fields. However, knowledge of detailed OFR chemistry is limited, allowing speculation in the literature about whether some non-OH reactants, including several not relevant for tropospheric chemistry, may play an important role in these OFRs. These non-OH reactants are UV radiation, O(1D), O(3P), and O3. In this study, we investigate the relative importance of other reactants to OH for the fate of reactant species in OFR under a wide range of conditions via box modeling. The relative importance of non-OH species is less sensitive to UV light intensity than to relative humidity (RH) and external OH reactivity (OHRext), as both non-OH reactants and OH scale roughly proportional to UV intensity. We show that for field studies in forested regions and also the urban area of Los Angeles, reactants of atmospheric interest are predominantly consumed by OH. We find that O(1D), O(3P), and O3 have relative contributions to VOC consumption that are similar or lower than in the troposphere. The impact of O atoms can be neglected under most conditions in both OFR and troposphere. Under pathological OFR conditions of low RH and/or high OHRext, the importance of non-OH reactants is enhanced because OH is suppressed. Some biogenics can have substantial destructions by O3, and photolysis at non-tropospheric wavelengths (185 and 254 nm) may also play a significant role in the degradation of some aromatics under pathological conditions. Working under low O2 with the OFR185 mode allows OH to completely dominate over O3 reactions even for the biogenic species most reactive with O3. Non-tropospheric VOC photolysis may have been a problem in some laboratory and source studies, but can be avoided or lessened in future studies by diluting source emissions and working at lower precursor concentrations in lab studies, and by humidification. SOA photolysis is shown to be insignificant for most functional groups, except for nitrates and especially aromatics, which may be photolyzed at high UV flux settings. Our work further establishes the OFR's usefulness as a tool to study atmospheric chemistry and enables better experiment design and interpretation, as well as improved future reactor design.

Peng, Zhe↗

Non-OH Chemistry in Oxidation Flow Reactors for the Study of Atmospheric Chemistry Systematically Examined by Modeling

Oxidation flow reactors (OFRs) using low-pressure Hg lamp emission at 185 and 254 nm produce OH radicals efficiently and are widely used in atmospheric chemistry and other fields. However, knowledge of detailed OFR chemistry is limited, allowing speculation in the literature about whether some non-OH reactants, including several not relevant for tropospheric chemistry, may play an important role in these OFRs. These non-OH reactants are UV radiation, O(1D), O(3P), and O3. In this study, we investigate the relative importance of other reactants to OH for the fate of reactant species in OFR under a wide range of conditions via box modeling. The relative importance of non-OH species is less sensitive to UV light intensity than to relative humidity (RH) and external OH reactivity (OHRext), as both non-OH reactants and OH scale roughly proportional to UV intensity. We show that for field studies in forested regions and also the urban area of Los Angeles, reactants of atmospheric interest are predominantly consumed by OH. We find that O(1D), O(3P), and O3 have relative contributions to VOC consumption that are similar or lower than in the troposphere. The impact of O atoms can be neglected under most conditions in both OFR and troposphere. Under pathological OFR conditions of low RH and/or high OHRext, the importance of non-OH reactants is enhanced because OH is suppressed. Some biogenics can have substantial destructions by O3, and photolysis at non-tropospheric wavelengths (185 and 254 nm) may also play a significant role in the degradation of some aromatics under pathological conditions. Working under low O2 with the OFR185 mode allows OH to completely dominate over O3 reactions even for the biogenic species most reactive with O3. Non-tropospheric VOC photolysis may have been a problem in some laboratory and source studies, but can be avoided or lessened in future studies by diluting source emissions and working at lower precursor concentrations in lab studies, and by humidification. SOA photolysis is shown to be insignificant for most functional groups, except for nitrates and especially aromatics, which may be photolyzed at high UV flux settings. Our work further establishes the OFR's usefulness as a tool to study atmospheric chemistry and enables better experiment design and interpretation, as well as improved future reactor design.

Peng, Zhe↗

Establishment of a genome editing tool using CRISPR-Cas9 ribonucleoprotein complexes in the non-model plant pathogen Sphaerulina musiva

CRISPR-Cas9 is a versatile genome editing system widely used since 2013 to introduce site-specific modifications into the genomes of model and non-model species. This technology is used in various applications, from gene knock-outs, knock-ins, and over-expressions to more precise changes, such as the introduction of nucleotides at a targeted locus. CRISPR-Cas9 has been demonstrated to be easy to establish in new species and highly efficient and specific compared to previous gene editing strategies such as Zinc finger nucleases and transcription activator-like effector nucleases. Grand challenges for emerging CRISPR-Cas9 tools in filamentous fungi are developing efficient transformation methods for non-model organisms. In this paper, we have leveraged the establishment of CRISPR-Cas9 genome editing tool that relies on Cas9/sgRNA ribonucleoprotein complexes (RNPs) in the model species Trichoderma reesei and developed the first protocol to efficiently transform the non-model species, Sphaerulina musiva. This fungal pathogen constitutes a real threat to the genus Populus, a foundational bioenergy crop used for biofuel production. Herein, we highlight the general considerations to design sgRNAs and their computational validation. We also describe the use of isolated protoplasts to deliver the CRISPR-Cas9 RNP components in both species and the screening for targeted genome editing events. The development of engineering tools in S. musiva can be used for studying genes involved in diverse processes such as secondary metabolism, establishment, and pathogenicity, among many others, but also for developing genetic mitigation approaches. The approach described here provides guidance for potential development of transformation systems in other non-model spore-bearing ascomycetes.

59 BASIC BIOLOGICAL SCIENCES↗

Interfacing HDF5 with a scalable object‐centric storage system on hierarchical storage

Summary Object storage technologies that take advantage of multitier storage on HPC systems are emerging. However, to use these technologies at present, applications have to be modified significantly from current I/O libraries. HDF5, a widely used I/O middleware on HPC systems, provides a virtual object layer (VOL) that allows applications to connect to different storage mechanisms transparently without requiring significant code modifications. We recently designed the proactive data containers (PDC) object‐centric storage system that provides the capabilities of transparent, asynchronous, and autonomous data movement taking advantage of multiple storage tiers—a decision that has so far been left upon the user on most current systems. To enable PDC's features through HDF5 without modifying application codes, we have developed an HDF5 VOL connector that interfaces with PDC. We present in this article the connector interface and evaluate its performance on Cori, a Cray XC40 supercomputer located at the National Energy Research Scientific Computing Center (NERSC). Our evaluation demonstrates up to an 8× improvement compared with HDF5 that has the most recent optimizations.

Mu, Jingqing↗

Compactly‐Supported Nonstationary Kernels for Computing Exact Gaussian Processes on Big Data

The Gaussian process (GP) is a widely used method for analyzing large-scale data sets, including spatio-temporal measurements of nonlinear processes that are now commonplace in the environmental sciences. Traditional implementations of GPs involve stationary kernels (also termed covariance functions) that limit their flexibility, and exact methods for inference that prevent application to data sets with more than about 10,000 points. Modern approaches to address stationarity assumptions generally fail to accommodate large data sets, while all attempts to address scalability focus on approximating the Gaussian likelihood, which can involve subjectivity and lead to inaccuracies. In this work, we explicitly derive an alternative kernel that can discover and encode both sparsity and nonstationarity. We embed the kernel within a fully Bayesian GP model and leverage high-performance computing resources to enable the analysis of massive data sets. We demonstrate the favorable performance of our novel kernel relative to existing exact and approximate GP methods across a variety of synthetic data examples. Furthermore, we conduct space–time prediction based on more than 1 million measurements of daily maximum temperature and verify that our results outperform state-of-the-art methods in the Earth sciences. More broadly, having access to exact GPs that use ultra-scalable, sparsity-discovering, nonstationary kernels allows GP methods to truly compete with a wide variety of machine learning methods.

Gaussian processes↗

Parallel performance of algebraic multigrid domain decomposition

Algebraic multigrid (AMG) is a widely used scalable solver and preconditioner for large-scale linear systems resulting from the discretization of a wide class of elliptic PDEs. While AMG has optimal computational complexity, the cost of communication has become a significant bottleneck that limits its scalability as processor counts continue to grow on modern machines. This article examines the design, implementation, and parallel performance of a novel algorithm, algebraic multigrid domain decomposition (AMG-DD), designed specifically to limit communication. The goal of AMG-DD is to provide a low-communication alternative to standard AMG V-cycles by trading some additional computational overhead for a significant reduction in communication cost. Numerical results show that AMG-DD achieves superior accuracy per communication cost compared with AMG, and speedup over AMG is demonstrated on a large GPU cluster.

97 MATHEMATICS AND COMPUTING↗