Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “APPL”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Direct-drive inertial confinement fusion: A review

In this study, the direct-drive, laser-based approach to inertial confinement fusion (ICF) is reviewed from its inception following the demonstration of the first laser to its implementation on the present generation of high-power lasers. The review focuses on the evolution of scientific understanding gained from target-physics experiments in many areas, identifying problems that were demonstrated and the solutions implemented. The review starts with the basic understanding of laser–plasma interactions that was obtained before the declassification of laser-induced compression in the early 1970s and continues with the compression experiments using infrared lasers in the late 1970s that produced thermonuclear neutrons. The problem of suprathermal electrons and the target preheat that they caused, associated with the infrared laser wavelength, led to lasers being built after 1980 to operate at shorter wavelengths, especially 0.35 um—the third harmonic of the Nd:glass laser—and 0.248 um (the KrF gas laser). The main physics areas relevant to direct drive are reviewed. The primary absorption mechanism at short wavelengths is classical inverse bremsstrahlung. Nonuniformities imprinted on the target by laser irradiation have been addressed by the development of a number of beam-smoothing techniques and imprint-mitigation strategies. The effects of hydrodynamic instabilities are mitigated by a combination of imprint reduction and target designs that minimize the instability growth rates. Several coronal plasma physics processes are reviewed. The two-plasmon–decay instability, stimulated Brillouin scattering (together with cross-beam energy transfer), and (possibly) stimulated Raman scattering are identified as potential concerns, placing constraints on the laser intensities used in target designs, while other processes (self-focusing and filamentation, the parametric decay instability, and magnetic fields), once considered important, are now of lesser concern for mainline direct-drive target concepts. Filamentation is largely suppressed by beam smoothing. Thermal transport modeling, important to the interpretation of experiments and to target design, has been found to be non-local in nature. Advances in shock timing and equation-of-state measurements relevant to direct-drive ICF are reported. Room-temperature implosions have provided an increased understanding of the importance of stability and uniformity. The evolution of cryogenic implosion capabilities, leading to an extensive series carried out on the 60-beam OMEGA laser [T. R. Boehly et al., Opt. Commun. 133, 495 (1997)], is reviewed together with major advances in cryogenic target formation. A polar-drive concept has been developed that will enable direct-drive–ignition experiments to be performed on the National Ignition Facility [C. A. Haynam et al., Appl. Opt. 46 (16), 3276 (2007)]. The advantages offered by the alternative approaches of fast ignition and shock ignition and the issues associated with these concepts are described. The lessons learned from target-physics and implosion experiments are taken into account in ignition and high-gain target designs for laser wavelengths of 1/3 μm and 1/4 μm. Substantial advances in direct-drive inertial fusion reactor concepts are reviewed. Overall, the progress in scientific understanding over the past five decades has been enormous, to the point that inertial fusion energy using direct drive shows significant promise as a future environmentally attractive energy source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Characterization of electrically-active defects in ultraviolet light-emitting diodes with laser-based failure analysis techniques

Laser-based failure analysis techniques demonstrate the ability to quickly and non-intrusively screen deep ultraviolet light-emitting diodes (LEDs) for electrically-active defects. In particular, two laser-based techniques, light-induced voltage alteration and thermally-induced voltage alteration, generate applied voltage maps (AVMs) that provide information on electrically-active defect behavior including turn-on bias, density, and spatial location. Here, multiple commercial LEDs were examined and found to have dark defect signals in the AVM indicating a site of reduced resistance or leakage through the diode. The existence of the dark defect signals in the AVM correlates strongly with an increased forward-bias leakage current. This increased leakage is not present in devices without AVM signals. Transmission electron microscopy analysis of a dark defect signal site revealed a dislocation cluster through the pn junction. The cluster included an open core dislocation. Even though LEDs with few dark AVM defect signals did not correlate strongly with power loss, direct association between increased open core dislocation densities and reduced LED device performance has been presented elsewhere [M. W. Moseley et al., J. Appl. Phys. 117, 095301 (2015)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Precise parameterization of the recombination velocity at passivated phosphorus doped surfaces

We investigate the surface recombination velocity S{sub p} at the silicon-dielectric interface of phosphorus-doped surfaces for two industrially relevant passivation schemes for crystalline silicon solar cells. A broad range of surface dopant concentrations together with a high accuracy of evaluating the latter is achieved by incremental back-etching of the surface. The analysis of lifetime measurements and the simulation of the surface recombination consistently apply a set of well accepted models, namely, the Auger recombination by Richter et al. [Phys. Rev. B 86, 1–14 (2012)], the carrier mobility by Klaassen [Solid-State Electron. 35, 953–959 (1992); 35, 961–967 (1992)], the intrinsic carrier concentration for undoped silicon by Altermatt et al. [J. Appl. Phys. 93, 1598–1604 (2003)], and the band-gap narrowing by Schenk [J. Appl. Phys. 84, 3684–3695 (1998)]. The results show an increased S{sub p} at textured in respect to planar surfaces. The obtained parameterizations are applicable in modern simulation tools such as EDNA [K. R. McIntosh and P. P. Altermatt, in Proceedings of the 35th IEEE Photovoltaic Specialists Conference, Honolulu, Hawaii, USA (2010), pp. 1–6], PC1Dmod [Haug et al., Sol. Energy Mater. Sol. Cells 131, 30–36 (2014)], and Sentaurus Device [Synopsys, Sentaurus TCAD, Zürich, Switzerland] as well as in the analytical solution under the assumption of local charge neutrality by Cuevas et al. [IEEE Trans. Electron Devices 40, 1181–1183 (1993)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Modulated electromagnetic fields in inhomogeneous media, hyperbolic pseudoanalytic functions, and transmutations

The time-dependent Maxwell system describing electromagnetic wave propagation in inhomogeneous isotropic media in the one-dimensional case reduces to a Vekua-type equation for bicomplex-valued functions of a hyperbolic variable, see Kravchenko and Ramirez [Adv. Appl. Cliord Algebr. 21(3), 547–559 (2011)]. Using this relation, we solve the problem of the transmission through an inhomogeneous layer of a normally incident electromagnetic time-dependent plane wave. The solution is written in terms of a pair of Darboux-associated transmutation operators [Kravchenko, V. V. and Torba, S. M., J. Phys. A: Math. Theor. 45, 075201 (2012)], and combined with the recent results on their construction [Kravchenko, V. V. and Torba, S. M., Complex Anal. Oper. Theory 9, 379-429 (2015); Kravchenko, V. V. and Torba, S. M., J. Comput. Appl. Math. 275, 1–26 (2015)] can be used for efficient computation of the transmitted modulated signals. We develop the corresponding numerical method and illustrate its performance with examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Laser-driven hydrothermal process studied with excimer laser pulses

Previously, we discovered [Mariella et al., J. Appl. Phys. 114, 014904 (2013)] that modest-fluence/modest-intensity 351-nm laser pulses, with insufficient fluence/intensity to ablate rock, mineral, or concrete samples via surface vaporization, still removed the surface material from water-submerged target samples with confinement of the removed material, and then dispersed at least some of the removed material into the water as a long-lived suspension of nanoparticles. We called this new process, which appears to include the generation of larger colorless particles, “laser-driven hydrothermal processing” (LDHP) [Mariella et al., J. Appl. Phys. 114, 014904 (2013)]. Now, we report that we have studied this process using 248-nm and 193-nm laser light on submerged concrete, quartzite, and obsidian, and, even though light at these wavelengths is more strongly absorbed than at 351 nm, we found that the overall efficiency of LDHP, in terms of the mass of the target removed per Joule of laser-pulse energy, is lower with 248-nm and 193-nm laser pulses than with 351-nm laser pulses. Given that stronger absorption creates higher peak surface temperatures for comparable laser fluence and intensity, it was surprising to observe reduced efficiencies for material removal. We also measured the nascent particle-size distributions that LDHP creates in the submerging water and found that they do not display the long tail towards larger particle sizes that we had observed when there had been a multi-week delay between experiments and the date of measuring the size distributions. This is consistent with transient dissolution of the solid surface, followed by diffusion-limited kinetics of nucleation and growth of particles from the resulting thin layer of supersaturated solution at the sample surface.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A photoemission moments model using density functional and transfer matrix methods applied to coating layers on surfaces: Theory

Recent experimental measurements of a bulk material covered with a small number of graphene layers reported by Yamaguchi et al. [NPJ 2D Mater. Appl. 1, 12 (2017)] (on bialkali) and Liu et al.[Appl. Phys. Lett. 110, 041607 (2017)] (on copper) and the needs of emission models in beam optics codes have lead to substantial changes in a Moments model of photoemission. The changes account for (i) a barrier profile and density of states factor based on density functional theory (DFT) evaluations, (ii) a Drude-Lorentz model of the optical constants and laser penetration depth, and (iii) a transmission probability evaluated by an Airy Transfer Matrix Approach. Importantly, the DFT results lead to a surface barrier profile of a shape similar to both resonant barriers and reflectionless wells: the associated quantum mechanical transmission probabilities are shown to be comparable to those recently required to enable the Moments (and Three Step) model to match experimental data but for reasons very different than the assumption by conventional wisdom that a barrier is responsible. The substantial modifications of the Moments model components, motivated by computational materials methods, are developed. The results prepare the Moments model for use in treating heterostructures and discrete energy level systems (e.g., quantum dots) proposed for decoupling the opposing metrics of performance that undermine the performance of advanced light sources like the x-ray Free Electron Laser. The consequences of the modified components on quan-tum yield, emittance, and emission models needed by beam optics codes are discussed. Published by AIP Publishing. https://doi.org/10.1063/1.5008600

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sound speed measurements in lithium fluoride single crystals shock compressed to 168 GPa along [100]

The shock wave response of [100] lithium fluoride (LiF) single crystals at high stresses is of long-standing interest due to their extensive use as optical windows in dynamic compression experiments. The report of melting in shock compressed LiF single crystals between 134 GPa and 152 GPa – based on a single sound speed datum [J. Appl. Phys. 117, 045901 (2015)] – was surprising because good optical transmission was previously demonstrated in LiF shock compressed to ~200 GPa [J. Appl. Phys. 116, 033515 (2014)]. To address these apparent differences, we report on plate impact experiments on [100] LiF single crystals shock compressed to 168 GPa. Wave profiles were measured using laser interferometry to determine Hugoniot states and longitudinal sound speeds in shock compressed LiF. The measured Hugoniot states are in good agreement with those measured in previous studies. However, the measured sound speeds presented here show no evidence of melting up to 168 GPa. In particular, the abrupt drop reported previously in sound speed at 152 GPa was not observed in the present work. Furthermore, our results establish a lower bound of 168 GPa for the onset of melting in shock compressed LiF single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effect of surface defects on damage and failure of dynamically driven cylindrical shells

Here, in this work, expanding cylinder experiments are performed using a gas gun on the aluminum 6061 alloy to understand the dynamic deformation and failure behavior of this material in the presence of surface defects in the form of scratches. The experiments are then modeled with the modified Tepla (for tensile plasticity) model, which is equipped with the damage yield function of Gurson [“Continuum theory of ductile rupture by void nucleation and growth: Part I—Yield criteria and flow rules for porous ductile media,” (1977)] and Tvergaard and Needleman [Acta Metall. 32(1), 157–169 (1984)], over-stress viscosity [F. Addessio and J. Johnson, J. Appl. Phys. 74(3), 1640–1648 (1993)], and the underlying viscoplastic model of Preston et al. [J. Appl. Phys. 93(1), 211–220 (2003)] (Preston–Tonks–Wallace model) for the dense material. The Tepla model parameters were calibrated using the Bayesian approach against experimentally measured free-surface velocity data obtained from plate impact experiments on the aluminum 6061 alloy. The evolution of radial expansion and shapes of the deformed cylinder from the simulations on a pristine cylinder wall and a cylinder wall with a longitudinal surface defect shows excellent agreement with the corresponding experimental results. Finally, the role of surface defect size and cylinder wall geometry on the dynamic strength of the cylinder wall was also investigated through a series of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2023

The Random Ray Method Versus Multigroup Monte Carlo: The Method of Characteristics in OpenMC and SCONE

The Random Ray Method (TRRM) is a recently developed approach to solving neutral particle transport problems based on the Method of Characteristics. While the method previously has been implemented only in closed-source or limited-functionality codes, this work describes its implementation in two open-source Monte Carlo codes: OpenMC and SCONE. The random ray implementations required small modifications to the existing Multigroup Monte Carlo (MGMC) solvers, offering a rare venue for redundant, fine-grained, "apples-to-apples" speed and accuracy comparisons between transport methods. To this end, TRRM and MGMC solvers are evaluated against each other using each code's native capabilities on reactor eigenvalue problems with different degrees of energy discretization. On the C5G7 benchmark (featuring only seven energy groups), TRRM achieves a maximum pin power error comparable to or lower than that of MGMC for a given run time. On a problem with 69 energy groups, MGMC is found to scale more efficiently, obtaining a lower pin power error for a given run time. However, the defining difference between the two transport methods is found to be their vastly different uncertainty distributions. Specifically, TRRM is found to maintain similar levels of accuracy and uncertainty throughout the simulation domain whereas MGMC can exhibit orders-of-magnitude greater errors in areas of the problem that feature low neutron flux. For instance, TRRM provided an up to 373 times speed advantage compared with MGMC for computing the flux in low-flux regions in the moderator surrounding the C5G7 core.

42 ENGINEERING↗
DOE OSTI2025

A comparison of effective field theory models of redshift space galaxy power spectra for DESI 2024 and future surveys

In preparation for the next generation of galaxy redshift surveys, and in particular the year-one data release from the Dark Energy Spectroscopic Instrument (DESI), we investigate the consistency of a variety of effective field theory models that describe the galaxy-galaxy power spectra in redshift space into the quasi-linear regime using 1-loop perturbation theory. These models are employed in the pipelines velocileptors, PyBird, and Folpsν. While these models have been validated independently, a detailed comparison with consistent choices has not been attempted. After briefly discussing the theoretical differences between the models we describe how to provide a more apples-to-apples comparison between them. We present the results of fitting mock spectra from the AbacusSummit suite of N-body simulations provided in three redshift bins to mimic the types of dark time tracers targeted by the DESI survey. We show that the theories behave similarly and give consistent constraints in both the forward-modeling and ShapeFit compressed fitting approaches. We additionally generate (noiseless) synthetic data from each pipeline to be fit by the others, varying the scale cuts in order to show that the models agree within the range of scales for which we expect 1-loop perturbation theory to be applicable. Finally, this work lays the foundation of Full-Shape analysis with DESI Y1 galaxy samples where in the tests we performed, we found no systematic error associated with the modeling of the galaxy redshift space power spectrum for this volume.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

X-ray back-diffraction: can we further increase the energy resolution by tuning the energy slightly below that of exact backscattering?

X-ray beams at energies tuned slightly below that of exact backscattering (extreme conditions, where X-ray back-diffraction is almost extinguished - called residual XBD) are better focused if the experiment is carried out at lower energies in order to avoid multiple-beam diffraction effects. Following previous work by the authors [Honnicke, Conley, Cusatis, Kakuno, Zhou, Bouet, Marques & Vicentin (2014). J. Appl. Cryst.47, 1658-1665], herein efforts are directed towards characterizing the residual XBD beam of an ultra-thin Si 220 crystal (UTSiXTAL) at similar to 3.2 keV. To achieve the residual XBD condition the UTSiXTAL was cooled from 310 to 273 K. The results indicate that under this extreme condition the energy resolution can be further improved. Issues with the energy resolution measurements due to incoming beam divergence and the ultra-thin crystal flatness are discussed.

low-energy X-ray back-diffraction↗
DOE OSTI2019

Partitioned Coupling for Structural Acoustics

Here, we expand the second-order fluid–structure coupling scheme of Farhat et al. (1998, “Load and Motion Transfer Algorithms for 19 Fluid/Structure Interaction Problems With Non-Matching Discrete Interfaces: Momentum and Energy Conservation, Optimal Discretization and Application to Aeroelasticity,” Comput. Methods Appl. Mech. Eng., 157(1–2), pp. 95–114; 2006, “Provably Second-Order Time-Accurate Loosely-Coupled Solution Algorithms for Transient Nonlinear Computational Aeroelasticity,” Comput. Methods Appl. Mech. Eng., 195(17), pp. 1973–2001) to structural acoustics. The staggered structural acoustics solution method is demonstrated to be second-order accurate in time, and numerical results are compared to a monolithically coupled system. The partitioned coupling method is implemented in the Sierra Mechanics software suite, allowing for the loose coupling of time domain acoustics in SIERRA/SD to structural dynamics ( SIERRA/SD ) or solid mechanics ( SIERRA/SD ). The coupling is demonstrated to work for nonconforming meshes. Results are verified for a one-dimensional piston, and the staggered and monolithic results are compared to an exact solution. Huang, H. (1969, “Transient Interaction of Plane Acoustic Waves With a Spherical Elastic Shell,” J. Acoust. Soc. Am., 45(3), pp. 661–670) sphere scattering problem with a spherically spreading acoustic load demonstrates parallel capability on a complex problem. Our numerical results compare well for a bronze plate submerged in water and sinusoidally excited (Fahnline and Shepherd, 2017, “Transient Finite Element/Equivalent Sources Using Direct Coupling and Treating the Acoustic Coupling Matrix as Sparse,” J. Acoust. Soc. Am., 142(2), pp. 1011–1024).

42 ENGINEERING↗
DOE OSTI2014

A modern Monte Carlo investigation of the TG-43 dosimetry parameters for an {sup 125}I seed already having AAPM consensus data

Purpose: To investigate potential causes for differences in TG-43 brachytherapy dosimetry parameters in the existent literature for the model IAI-125A{sup 125}I seed and to propose new standard dosimetry parameters. Methods: The MCNP5 code was used for Monte Carlo (MC) simulations. Sensitivity of dose distributions, and subsequently TG-43 dosimetry parameters, was explored to reproduce historical methods upon which American Association of Physicists in Medicine (AAPM) consensus data are based. Twelve simulation conditions varying{sup 125}I coating thickness, coating mass density, photon interaction cross-section library, and photon emission spectrum were examined. Results: Varying{sup 125}I coating thickness, coating mass density, photon cross-section library, and photon emission spectrum for the model IAI-125A seed changed the dose-rate constant by up to 0.9%, about 1%, about 3%, and 3%, respectively, in comparison to the proposed standard value of 0.922 cGy h{sup −1} U{sup −1}. The dose-rate constant values by Solberg et al. [“Dosimetric parameters of three new solid core {sup 125}I brachytherapy sources,” J. Appl. Clin. Med. Phys. 3, 119–134 (2002)], Meigooni et al. [“Experimental and theoretical determination of dosimetric characteristics of IsoAid ADVANTAGE™ {sup 125}I brachytherapy source,” Med. Phys. 29, 2152–2158 (2002)], and Taylor and Rogers [“An EGSnrc Monte Carlo-calculated database of TG-43 parameters,” Med. Phys. 35, 4228–4241 (2008)] for the model IAI-125A seed and Kennedy et al. [“Experimental and Monte Carlo determination of the TG-43 dosimetric parameters for the model 9011 THINSeed™ brachytherapy source,” Med. Phys. 37, 1681–1688 (2010)] for the model 6711 seed were +4.3% (0.962 cGy h{sup −1} U{sup −1}), +6.2% (0.98 cGy h{sup −1} U{sup −1}), +0.3% (0.925 cGy h{sup −1} U{sup −1}), and −0.2% (0.921 cGy h{sup −1} U{sup −1}), respectively, in comparison to the proposed standard value. Differences in the radial dose functions between the current study and both Solberg et al. and Meigooni et al. were <10% for r ≤ 5 cm, and increased for r > 5 cm with a maximum difference of 29% at r = 9 cm. In comparison to Taylor and Rogers, these differences were lower (maximum of 2% at r = 9 cm). For the similarly designed model 6711 {sup 125}I seed, differences did not exceed 0.5% for 0.5 ≤ r ≤ 10 cm. Radial dose function values varied by 1% as coating thickness and coating density were changed. Varying the cross-section library and source spectrum altered the radial dose function by 25% and 12%, respectively, but these differences occurred at r = 10 cm where the dose rates were very low. The 2D anisotropy function results were most similar to those of Solberg et al. and most different to those of Meigooni et al. The observed order of simulation condition variables from most to least important for influencing the 2D anisotropy function was spectrum, coating thickness, coating density, and cross-section library. Conclusions: Several MC radiation transport codes are available for calculation of the TG-43 dosimetry parameters for brachytherapy seeds. The physics models in these codes and their related cross-section libraries have been updated and improved since publication of the 2007 AAPM TG-43U1S1 report. Results using modern data indicated statistically significant differences in these dosimetry parameters in comparison to data recommended in the TG-43U1S1 report. Therefore, it seems that professional societies such as the AAPM should consider reevaluating the consensus data for this and others seeds and establishing a process of regular evaluations in which consensus data are based upon methods that remain state-of-the-art.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Iterated Gauss-Seidel GMRES

The GMRES algorithm of Saad and Schultz [SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869] is an iterative method for approximately solving linear systems Ax = b, with initial guess x0 and residual r0 = b Ax0. The algorithm employs the Arnoldi process to generate the Krylov basis vectors (the columns of Vk ). It is well known that this process can be viewed as a QR factorization of the matrix Bk = [r0, AVk] at each iteration. Despite an O (..epsilon..)..kappa.. (Bk ) loss of orthogonality, for unit roundoff ..epsilon..and condition number ..kappa.. , the modified Gram-Schmidt formulation was shown to be backward stable in the seminal paper by Paige et al. [SIAM J. Matrix Anal.Appl., 28 (2006), pp. 264-284]. We present an iterated Gauss-Seidel formulation of the GMRES algorithm (IGS-GMRES) based on the ideas of Ruhe [Linear Algebra Appl., 52 (1983), pp. 591-601] and Swirydowicz et al. [Numer. Linear Algebra Appl., 28 (2020), pp. 1-20]. IGS-GMRES maintains orthogonality to the level O (..epsilon..)..kappa.. (Bk ) or O (..epsilon..), depending on the choice of one or two iterations; for two Gauss-Seidel iterations, the computed Krylov basis vectors remain orthogonal to working accuracy and the smallest singular value of Vk remains close to one. The resulting GMRES method is thus backward stable. We show that IGS-GMRES can be implemented with only a single synchronization point per iteration, making it relevant to large-scale parallel computing environments. We also demonstrate that, unlike MGS-GMRES, in IGS-GMRES the relative Arnoldi residual corresponding to the computed approximate solution no longer stagnates above machine precision even for highly nonnormal systems.

Arnoldi-QR↗
DOE OSTI2016

Evaluation of Graph Pattern Matching Workloads in Graph Analysis Systems

Graph analysis has emerged as a powerful method for data scientists to represent, integrate, query, and explore heterogeneous data sources. As a result, graph data management and mining became a popular area of research, and led to the development of plethora of systems in recent years. Unfortunately, the number of emerging graph analysis systems and the wide range of applications, coupled with a lack of apples-to-apples comparisons, make it difficult to understand the trade-offs between different systems and the graph operations for which they are designed. A fair comparison of these systems is a challenging task for the following reasons: multiple data models, non-standardized serialization formats, various query interfaces to users, and diverse environments they operate in. To address these key challenges, in this paper we present a new benchmark suite by extending the Lehigh University Benchmark (LUBM) to cover the most common capabilities of various graph analysis systems. We provide the design process of the benchmark, which generalizes the workflow for data scientists to conduct the desired graph analysis on different graph analysis systems. Equipped with this extended benchmark suite, we present performance comparison for nine subgraph pattern retrieval operations over six graph analysis systems, namely NetworkX, Neo4j, Jena, Titan, GraphX, and uRiKA. Through the proposed benchmark suite, this study reveals both quantitative and qualitative findings in (1) implications in loading data into each system; (2) challenges in describing graph patterns for each query interface; and (3) different sensitivity of each system to query selectivity. We envision that this study will pave the road for: (i) data scientists to select the suitable graph analysis systems, and (ii) data management system designers to advance graph analysis systems.

Hong, Seokyong↗
DOE OSTI2026

Emerging Flexible Designs for Geospatial Multimodal Foundation Models

Foundation models are rapidly transforming Earth observation by enabling scalable pretraining across diverse unlabeled geospatial modalities. However, their architectural diversity—ranging from encoder-only to encoder-decoder and masked autoencoding paradigms—makes it challenging to assess performance trade-offs in a consistent manner. In this work, we present an apples-to-apples comparison of leading FM architectures designed for geospatial multimodal reasoning, with a particular focus on flexibility across varied spectral band configurations. We standardize pretraining using identical self-supervised learning objectives and training datasets, and evaluate all models under consistent parameterization on the GEOBench benchmark across classification and segmentation tasks. Our results offer new insights into the design trade-offs between model flexibility, modality alignment, and downstream task performance. By highlighting architectural strengths and limitations under controlled conditions, this study provides practical guidance for building next-generation geospatial foundation models capable of robust multimodal reasoning.

Ambrozio Dias, Philipe [ORNL] (ORCID:0000000194277↗
DOE OSTI2016

GraphBench

GraphBench is a benchmark suite for graph pattern mining and graph analysis systems. The benchmark suite is a significant addition to conducting apples-apples comparison of graph analysis software (databases, in-memory tools, triple stores, etc.)

Sukumar, SreenivasR.↗