Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Simple Two-Dimensional Corrections for One-Dimensional Pulse Tube Models

One-dimensional oscillating flow models are very useful for designing pulse tubes. They are simple to use, not computationally intensive, and the physical relationship between temperature, pressure and mass flow are easy to understand when used in conjunction with phasor diagrams. They do not possess, however, the ability to directly calculate thermal and momentum diffusion in the direction transverse to the oscillating flow. To account for transverse effects, lumped parameter corrections, which are obtained though experiment, must be used. Or two-dimensional solutions of the differential fluid equations must be obtained. A linear two-dimensional solution to the fluid equations has been obtained. The solution provides lumped parameter corrections for one-dimensional models. The model accounts for heat transfer and shear flow between the gas and the tube. The complex Nusselt number and complex shear wall are useful in describing these corrections, with phase relations and amplitudes scaled with the Prandtl and Valensi numbers. The calculated ratio, a, between a two-dimensional solution of the oscillating temperature and velocity and a one-dimensional solution for the same shows a scales linearly with Va for Va less than 30. In this region alpha less than 0.5, that is, the enthalpy flow calculated with a two-dimensional model is 50% of a calculation using a one-dimensional model. For Va greater than 250, alpha = 0.8, showing that diffusion is still important even when it is confined to a thing layer near the tube wall.

Lee, J. M.↗
DOE OSTI2025

Spontaneous generation of athermal phonon bursts within bulk silicon causing excess noise, low energy background events, and quasiparticle poisoning in superconducting sensors

Solid state phonon detectors used in the search for dark matter and coherent neutrino nucleus interactions (CE $v$ NS) require excellent energy resolution (eV-scale or below) and low backgrounds. An unknown source of phonon bursts, the low energy excess (LEE), dominates other above-threshold backgrounds and generates excess shot noise from subthreshold bursts. Here, in this paper, we measure these phonon bursts for 12 days after cooldown in two nearly identical 1 cm 2 silicon detectors that differ only in the thickness of their substrate (1 vs 4 mm thick). We find that both the channel-correlated shot noise and near-threshold shared LEE relax with time since cooldown. Additionally, both the correlated shot noise and LEE rates scale linearly with substrate thickness. When combined with previous measurements of other silicon phonon detectors with different substrate geometries and mechanical support strategies, these measurements strongly suggest that the dominant source of both above and below threshold LEE is the bulk substrate. By monitoring the relation between bias power and excess phonon shot noise, we estimate that the energy scale for subthreshold noise events is 0.68 ± 0.38 meV. In our final dataset, we report a world-leading energy resolution of 258.5 ± 0.4 meV in the 1 mm thick detector. Simple calculations suggest that these silicon substrate phonon bursts are likely a significant source of quasiparticle poisoning in superconducting qubits operated in well shielded and vibration free environments.

Chang, C. L. [Argonne National Laboratory 1 , 9700↗
NASA NTRS1982

Pole placement and order reduction in two-time-scale control systems through Riccati iteration

A transformation of variables taken from singular perturbations may be applied to two-time-scale linear systems in state space form to reduce the system to block-diagonal form with slow and fast modes decoupled. The transformation is easily computed by applying the new Riccati iteration. The iteration yields a solution to the nonsymmetric algebraic Riccati equation obtained by partitioning the original system matrix A. The numerical procedure is initiated with the trivial iterate L(0) = 0, and is globally convergent to the desired unique time scale decoupling solution. After transformation, the decoupled system may be used in controller design to achieve exact closed loop pole placement in the slow subsystem without altering the poles of the fast subsystem. The decoupled form may also be used to reduce system order by wetting a small parameter to zero. Provided the fast subsystem is stable, the order reduction can be expected to yield a good approximation to the original system. These methods are demonstrated using the 16th order linear model of a turbofan engine.

Anderson, L. R.↗
DOE OSTI2022

Hydrogen-Bonding Bottlebrush Networks: Self-Healing Materials from Super-Soft to Stiff

The impact of polymer architecture on network dynamics and self-healing is presented using bottlebrushes containing side chains that are end-functionalized with 2-ureido-4[1H]-pyrimidinone (UPy). The synthesis of these materials is straightforward through a three-step process: (1) synthesizing rubbery poly(4-methylcaprolactone) macromonomers (p4MCL–OH) with a norbornene-based initiator, (2) functionalizing the terminal hydroxyl group with UPy–isocyanate (p4MCL–UPy), and (3) statistically copolymerizing p4MCL–OH and p4MCL–UPy via ring-opening metathesis polymerization (ROMP) to form hydrogen-bonding bottlebrushes having a fraction (p) of side chains functionalized with UPy. Attaching UPy to the free end of bottlebrush side chains dilutes the impact of friction from complementary UPy interactions on segmental dynamics, leading to a much weaker dependence of the glass-transition temperature (T g ) on p than observed in linear analogues, while the activation energy to dissociate UPy–UPy bonds (41–47 kJ/mol) remains mostly unchanged. Longer side chains result in a competition between reducing T g and inducing entanglements that influence hydrogen-bonded network dynamics. Increasing the backbone length extends the sticky Rouse region without affecting the network modulus (G x ) or UPy–UPy dissociation time (τ s ). Gx scales linearly with p and ranges from 27 kPa to 1.6 MPa, while τs remains nearly constant in contrast to linear telechelic ionomers, implying a similar self-healability across bottlebrush networks containing different amounts of UPy. Furthermore, these stretchable networks with p ≤ 0.25 undergo self-healing upon repeated rupture and melt pressing at ≤100 °C while retaining similar tensile properties. In summary, decorating bottlebrush polymers with hydrogen bonds creates opportunities to independently manipulate associative network dynamics and mechanical moduli.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effective cosmic density field reconstruction with convolutional neural network

ABSTRACT We present a cosmic density field reconstruction method that augments the traditional reconstruction algorithms with a convolutional neural network (CNN). Following previous work, the key component of our method is to use the reconstructed density field as the input to the neural network. We extend this previous work by exploring how the performance of these reconstruction ideas depends on the input reconstruction algorithm, the reconstruction parameters, and the shot noise of the density field, as well as the robustness of the method. We build an eight-layer CNN and train the network with reconstructed density fields computed from the Quijote suite of simulations. The reconstructed density fields are generated by both the standard algorithm and a new iterative algorithm. In real space at z = 0, we find that the reconstructed field is 90 per cent correlated with the true initial density out to $k\sim 0.5 \, \mathrm{ h}\, \rm {Mpc}^{-1}$, a significant improvement over $k\sim 0.2 \, \mathrm{ h}\, \rm {Mpc}^{-1}$ achieved by the input reconstruction algorithms. We find similar improvements in redshift space, including an improved removal of redshift space distortions at small scales. We also find that the method is robust across changes in cosmology. Additionally, the CNN removes much of the variance from the choice of different reconstruction algorithms and reconstruction parameters. However, the effectiveness decreases with increasing shot noise, suggesting that such an approach is best suited to high density samples. This work highlights the additional information in the density field beyond linear scales as well as the power of complementing traditional analysis approaches with machine learning techniques.

Astronomy & Astrophysics↗
DOE OSTI2025

Realization and Calibration of Continuously Parameterized Two-Qubit Gates on a Trapped-Ion Quantum Processor

Continuously parameterized two-qubit gates are a key feature of state-of-the-art trapped-ion quantum processors, as they have favorable error scalings and show distinct improvements in circuit performance over more restricted maximally entangling gatesets. In this work, we provide a comprehensive and pedagogical discussion on how to practically implement these continuously parameterized Mølmer–Sørensen gates on the Quantum Scientific Computing Open User Testbed, a low-level trapped-ion processor. To generate the arbitrary entangling angles, θ, we simply scale the amplitude of light used to generate the entanglement. However, doing so requires careful consideration of amplifier saturation as well as the variable light shifts that result. As such, we describe a method to calibrate and cancel the dominant fourth-order effects, followed by a dynamic virtual phase advance during the gate to cancel any residual light shifts, and find a linear scaling between θ and the residual light shift. Once we have considered and calibrated these effects, we demonstrate performance improvement with decreasing θ. Finally, we describe nuances of hardware control to transform the XX-type interaction of the arbitrary-angle Mølmer–Sørensen gate into a phase-agnostic and crosstalk-mitigating ZZ interaction.

Yale, Christopher Gordon [Sandia National Laborato↗
DOE OSTI2020

Breaking Simple Scaling Relations through Metal–Oxide Interactions: Understanding Room-Temperature Activation of Methane on M/CeO 2 (M = Pt, Ni, or Co) Interfaces

The clean activation of methane at low temperatures remains an eminent challenge and a field of competitive research. In particular, on late transition metal surfaces such as Pt(111) or Ni(111), higher temperatures are necessary to activate the hydrocarbon molecule, but a massive deposition of carbon makes the metal surface useless for catalytic activity. However, on very low-loaded M/CeO 2 (M = Pt, Ni, or Co) surfaces, the dissociation of methane occurs at room temperature, which is unexpected considering simple linear scaling relationships. This intriguing phenomenon has been studied using a combination of experimental techniques (ambient-pressure X-ray photoelectron spectroscopy, time-resolved X-ray diffraction, and X-ray absorption spectroscopy) and density functional theory-based calculations. The experimental and theoretical studies show that the size and morphology of the supported nanoparticles together with strong metal–support interactions are behind the deviations from the scaling relations. These findings point toward a possible strategy for circumventing scaling relations, producing active and stable catalysts that can be employed for methane activation and conversion.

03 NATURAL GAS↗
DOE OSTI2023

A high-resolution large-eddy simulation framework for wildland fire predictions using TensorFlow

Background: Wildfires are becoming more severe, so we need improved tools to predict them over a wide range of conditions and scales. One approach towards this goal entails the use of coupled fire/atmosphere modelling tools. Although significant progress has been made in advancing their physical fidelity, existing tools have not taken full advantage of emerging programming paradigms and computing architectures to enable high-resolution wildfire simulations. Aims: The aim of this study was to present a new framework that enables landscape-scale wildfire simulations with physical representation of combustion at an affordable cost. Methods: We developed a coupled fire/atmosphere simulation framework using TensorFlow, which enables efficient and scalable computations on Tensor Processing Units. Key Results: Simulation results for a prescribed fire were compared with experimental data. Predicted fire behavior and statistical analysis for fire spread rate, scar area, and intermittency showed overall reasonable agreement. Scalability analysis was performed, showing close to linear scaling. Conclusions: While mesh refinement was shown to have less impact on global quantities, such as fire scar area and spread rate, it benefits predictions of intermittent fire behavior, buoyancy-driven dynamics, and small-scale turbulent motion. Implications: This new simulation framework is efficient in capturing both global quantities and unsteady dynamics of wildfires at high spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

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

The impact of the fiducial cosmology assumption on BAO distance scale measurements

ABSTRACT Standard analysis pipelines for measurements of Baryon Acoustic Oscillations (BAO) in galaxy surveys make use of a fiducial cosmological model to guide the data compression required to transform from observed redshifts and angles to the measured angular and radial BAO peak positions. In order to remove any dependence on the fiducial cosmology from the results, all models compared to the data should mimic the compression and its dependence on the fiducial model. In practice, approximations are made when testing models: (1) There is assumed to be no residual dependence on the fiducial cosmology after reconstruction, (2) differences in the distance–redshift relationship are assumed to match a linear scaling, and (3) differences in clustering between true and fiducial models are assumed to be removed by the free parameters used to null the non-BAO signal. We test these approximations using the current standard measurement procedure with a set of halo catalogues from the aemulus suite of N-body simulations, which span a range of wCDM cosmological models. We focus on reconstruction of the primordial BAO and locating the BAO. For the range of wCDM cosmologies covered by the aemulus suite, we find no evidence for systematic errors in the measured BAO shift parameters α∥ and α⊥ to $\lt 0.1\%$. However, the measured errors $\sigma _{\alpha _{\parallel }}$ and $\sigma _{\alpha _{\bot }}$ show a notable absolute increase by up to +0.001 and +0.002, respectively, in the case that the fiducial cosmology does not match the truth. These effects on the inferred BAO scale will be important, given the precision of measurements expected from future surveys including DESI, Euclid, and WFIRST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Connecting Hindered Transport in Porous Media across Length Scales: From Single-Pore to Macroscopic

Hindered mass transport is widely observed in various porous media, including filtration membranes and porous rocks. However, there is no universal model capable of predicting transport in porous media due to the heterogeneity of porous structures and the complexity of the underlying microscopic mechanisms, including electrostatic and hydrodynamic effects. Here, we used a highly ordered porous medium (i.e. a silica inverse opal film) as a model system to directly explore the effects of geometric parameters (i.e. pore size, pore throat size and tracer particle size) and microscopic interaction parameters (e.g. controlled by ionic strength) on nanoparticle transport in porous environments using three-dimensional single-particle tracking. We found a linear scaling relation between the macroscopic diffusion coefficient and microscopic diffusion behaviors and a combination of parameters associated with microscopic features and phenomena, including both geometric effects and particle-level interactions. The proportionality coefficient relating micro and macro behaviors was complex and related to the connectivity of the matrix and pore-size variation, which could lead to tortuous diffusion pathways, hindering macroscopic transport. Furthermore, these findings provide important insights into nanoparticles and macromolecules transport in various porous media.

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

Scalable training of graph convolutional neural networks for fast and accurate predictions of HOMO-LUMO gap in molecules

Abstract Graph Convolutional Neural Network (GCNN) is a popular class of deep learning (DL) models in material science to predict material properties from the graph representation of molecular structures. Training an accurate and comprehensive GCNN surrogate for molecular design requires large-scale graph datasets and is usually a time-consuming process. Recent advances in GPUs and distributed computing open a path to reduce the computational cost for GCNN training effectively. However, efficient utilization of high performance computing (HPC) resources for training requires simultaneously optimizing large-scale data management and scalable stochastic batched optimization techniques. In this work, we focus on building GCNN models on HPC systems to predict material properties of millions of molecules. We use HydraGNN, our in-house library for large-scale GCNN training, leveraging distributed data parallelism in PyTorch. We use ADIOS, a high-performance data management framework for efficient storage and reading of large molecular graph data. We perform parallel training on two open-source large-scale graph datasets to build a GCNN predictor for an important quantum property known as the HOMO-LUMO gap. We measure the scalability, accuracy, and convergence of our approach on two DOE supercomputers: the Summit supercomputer at the Oak Ridge Leadership Computing Facility (OLCF) and the Perlmutter system at the National Energy Research Scientific Computing Center (NERSC). We present our experimental results with HydraGNN showing (i) reduction of data loading time up to 4.2 times compared with a conventional method and (ii) linear scaling performance for training up to 1024 GPUs on both Summit and Perlmutter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

First-Generation Jet Propulsion Laboratory "Hockey-Puck" Free-Flying Magnetometers for Distributed In-Situ Multiprobe Measurement of Current Density Filamentation in the Northern Auroral Zone: Enstrophy Mission

The sub-orbital rocket mission was a collaborative project between the University of New Hampshire, Cornell University, and the Jet Propulsion Laboratory (JPL) to study filamentation phenomena in the northern Auroral zone. The Enstrophy mission test flies the JPL Free-Flying Magnetometer (FFM) concept. The FFM technology development task has been funded by NASA develop miniaturized, low-power, integrated "sensorcrafts". JPL's role was to design, integrate, test, and deliver four FFMs for deployment from the sounding rocket, allowing a unique determination of curl-B. This provides a direct measurement of magnetic-field-aligned current density along the rocket trajectory. A miniaturized three-axis fluxgate magnetometer was integrated with a 4-channel 22-bit sigma-delta Analog to Digital Converter (ADC), four temperature sensors, digital control electronics, seven (Li-SOCl2) batteries, two (4 deg x 170 deg field of view) sun-sensors, a fan-shaped-beam laser diode beacon, a (16 MHz) stable Temperature Compensated Crystal Oscillator (TCXO) clock, Radio Frequency (RF) communication subsystem, and an antenna for approximately 15 minutes of operation where data was collected continuously and transmitted in three (3) bursts (approximately 26 seconds each) to ground station antennas at Poker Flat, Alaska. FFMs were stowed within two trays onboard the rocket during the rocket launch and were released simultaneously using the spinning action of the rocket at approximately 300 km altitude (approximately 100 sec. into the flight). FFMs were deployed with spin rate of approximately 17 Hz and approximately 3 m/sec linear velocity with respect to the rocket. For testing purposes while the rocket was in the launch pad and during flight prior to release of FFMs from the rocket, commands (such as "power on", "test", "flight", "power off', and clock "Reset" signal) were transmitted via a infrared Light Emitting Diode to an infrared detector in the FFM. Special attention was paid to low magnetic signature electronic design and choice of materials in packaging. The miniaturized fluxgate magnetometers had a range of 1-60000 nT with 0.1% full-scale linearity. The frequency range of interest for magnetic measurement was 10 mHz - 50 Hz. Digital data from the magnetometer's three axes were placed in a 4MB Static Random Access Memory (SRAM) in data packages (frames) formatted together with time tags and frame ID. After a specified time was elapsed, the data were Viterbi encoded and transmitted at a rate of 100 kbps (BPSK). Each of the four FFMs transmitted at different frequency. These carrier frequencies were in the range of 2200-2300 MHz. The antenna was a single patch on a high dielectric constant substrate covering one end-plate of the hockey-puck-sized unit. The local clocks aboard the FFMs were reset at the start of the mission and stayed synchronized within 3 msec during the mission. Position of each FFM with respect to the rocket is calculated by the knowledge of its release velocity (measured at exit point of the FFM launcher tract) providing an accuracy of 1 m over the maximum range of 3 km. Spatial and temporal nature of observants can be separated to within 3 m in space or 3 msec time interval.

Javadi, H.↗
NASA NTRS2003

Ocular motor responses to abrupt interaural head translation in normal humans

We characterized the interaural translational vestibulo-ocular reflex (tVOR) in 6 normal humans to brief (approximately 200 ms), high-acceleration (0.4-1.4g) stimuli, while they fixed targets at 15 or 30 cm. The latency was 19 +/- 5 ms at 15-cm and 20 +/- 12 ms at 30-cm viewing. The gain was quantified using the ratio of actual to ideal behavior. The median position gain (at time of peak head velocity) was 0.38 and 0.37, and the median velocity gain, 0.52 and 0.62, at 15- and 30-cm viewing, respectively. These results suggest the tVOR scales proportionally at these viewing distances. Likewise, at both viewing distances, peak eye velocity scaled linearly with peak head velocity and gain was independent of peak head acceleration. A saccade commonly occurred in the compensatory direction, with a greater latency (165 vs. 145 ms) and lesser amplitude (1.8 vs. 3.2 deg) at 30- than 15-cm viewing. Even with saccades, the overall gain at the end of head movement was still considerably undercompensatory (medians 0.68 and 0.77 at 15- and 30-cm viewing). Monocular viewing was also assessed at 15-cm viewing. In 4 of 6 subjects, gains were the same as during binocular viewing and scaled closely with vergence angle. In sum the low tVOR gain and scaling of the response with viewing distance and head velocity extend previous results to higher acceleration stimuli. tVOR latency (approximately 20 ms) was lower than previously reported. Saccades are an integral part of the tVOR, and also scale with viewing distance.

Non-NASA Center↗
DOE OSTI2023

DESI mock challenge: constructing DESI galaxy catalogues based on FastPM simulations

Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large-scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a halo occupation distribution (HOD) model upon the FastPM simulations, such that the resulting galaxy clustering reproduces high-resolution N-body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements – monopole and quadrupole – to a precision required by the DESI Year 1 emission line galaxy sample down to non-linear scales, i.e. $k\lt 0.5\, h\, \mathrm{Mpc}^{-1}$ or $s\gt 10\, \mathrm{Mpc}\, h^{-1}$. Furthermore, we compute covariance matrices based on the resulting FastPM galaxy clustering – monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate (e.g. power spectrum) on the covariance matrix of the Fourier counterpart (e.g. correlation function). We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between 2 and 10 per cent. Consequently, using the current pipeline we can precisely reproduce the clustering of N-body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

CGSim: A Simulation Framework for Large Scale Distributed Computing Environment

Large-scale distributed computing infrastructures such as the Worldwide LHC Computing Grid (WLCG) require comprehensive simulation tools for evaluating performance, testing new algorithms, and optimizing resource allocation strategies. However, existing simulators suffer from limited scalability, hardwired algorithms, lack of real-time monitoring, and inability to generate datasets suitable for modern machine learning approaches. We present CGSim, a simulation framework for large-scale distributed computing environments that addresses these limitations. Built upon the validated SimGrid simulation framework, CGSim provides high-level abstractions for modeling heterogeneous grid environments while maintaining accuracy and scalability. Key features include a modular plugin mechanism for testing custom workflow scheduling and data movement policies, interactive real-time visualization dashboards, and automatic generation of event-level datasets suitable for AI-assisted performance modeling. We demonstrate CGSim’s capabilities through a comprehensive evaluation using production ATLAS PanDA workloads, showing significant calibration accuracy improvements across WLCG computing sites. Scalability experiments show near-linear scaling for multi-site simulations, with distributed workloads achieving 6 × better performance compared to single-site execution. The framework enables researchers to simulate WLCG-scale infrastructures with hundreds of sites and thousands of concurrent jobs within practical time budget constraints on commodity hardware.

Vatsavai, Sairam Sri [Brookhaven National Laborato↗
DOE OSTI2023

Scalable DPG multigrid solver for Helmholtz problems: A study on convergence

This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov–Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid preconditioning techniques of Petrides and Demkowicz (Comput. Math. Appl. 87 (2021) pp. 12–26) and extends the convergence results from $\mathscr{O}$(10 7 ) degrees of freedom (DOFs) to $\mathscr{O}$(10 9 ) DOFs using a new scalable parallel MPI/OpenMP implementation. Novel contributions of this paper include an alternative definition of coarse-grid systems based on restriction of fine-grid operators, yielding superior convergence results. In the uniform refinement setting, a detailed convergence study is provided, demonstrating h and p robust convergence and linear scaling with respect to the wave frequency. Finally, the paper concludes with numerical results on hp -adaptive simulations including a large-scale seismic modeling benchmark problem with high material contrast.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Linearized Pair-Density Functional Theory for Vertical Excitation Energies

Multiconfiguration pair-density functional theory (MC-PDFT) is a computationally efficient method that computes the energies of electronic states in a state specific or state average framework via an on-top functional. However, MC-PDFT does not include state interaction among these states since the final energies do not come from the diagonalization of an effective model-space Hamiltonian. Recently, multistate extensions such as linearized PDFT (L-PDFT) have been developed to accurately model the potentials near conical intersections and avoided crossings. However, there has not been any systematic study evaluating their performance for predicting vertical excitations at the equilibrium geometry of a molecule, when the excited states are generally well separated. In this paper, we report the performance of L-PDFT on the extensive QUESTDB data set of vertical excitations using a database of automatically selected active spaces. We show that L-PDFT performs well on all these excitations and successfully reproduces the performance of MC-PDFT. These results further demonstrate the potential of L-PDFT, as its scaling is constant with the number of states included in the state-average manifold, whereas MC-PDFT scales linearly in this regard.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗